13 #define ODIR "brst_o1_lf" // output plot directory 19 gROOT->SetBatch(
true);
30 double F_GA[4] = {0.0001,0.001,0.0025,0.004};
31 for(
int n=0;
n<4;
n++) {
32 par[0].name=
"duration"; par[0].value=F_GA[
n];
41 double F_Q3[3] = {70,235,849};
42 for(
int n=0;
n<3;
n++) {
43 par[0].name=
"frequency"; par[0].value=F_Q3[
n];
44 par[1].name=
"Q"; par[1].value=3.;
53 double F_Q8d9[8] = {70,100,153,235,361,554,849,1053};
54 double Q[8] = {8.9,8.9,8.9,8.9,8.9,8.9,8.9,9};
55 for(
int n=0;
n<8;
n++) {
56 par[0].name=
"frequency"; par[0].value=F_Q8d9[
n];
57 par[1].name=
"Q"; par[1].value=Q[
n];
58 if((F_Q8d9[
n]==70)||(F_Q8d9[
n]==100)||(F_Q8d9[
n]==235)||(F_Q8d9[
n]==361)||(F_Q8d9[
n]==1053)) {
70 double F_Q100[3] = {70,235,849};
71 for(
int n=0;
n<3;
n++) {
72 par[0].name=
"frequency"; par[0].value=F_Q100[
n];
73 par[1].name=
"Q"; par[1].value=100.;
82 double F_WNB[2] = {100,250};
83 double B_WNB[2] = {100,100};
84 double D_WNB[2] = {0.1,0.1};
85 for(
int n=0;
n<2;
n++) {
86 par[0].name=
"frequency"; par[0].value=F_WNB[
n];
87 par[1].name=
"bandwidth"; par[1].value=B_WNB[
n];
88 par[2].name=
"duration"; par[2].value=D_WNB[
n];
89 for(
int m=0;
m<30;
m++) {
90 par[3].name=
"pseed"; par[3].value=100000+
n*100+
m;
91 par[4].name=
"xseed"; par[4].value=100001+
n*100+
m;
92 par[5].name=
"mode"; par[5].value=0;
102 for(
int n=0;
n<wfList.size();
n++) {
108 cout <<
"Error : Waveform " << wf.
name <<
" not exist in the MDC pool !!!" << endl << endl;
112 cout <<
n <<
" " << wf.
name <<
" size : " << wf.
hp.
size() <<
" rate : " << wf.
hp.
rate()
120 gSystem->mkdir(TString::Format(
"%s",
ODIR));
121 gSystem->mkdir(TString::Format(
"%s/plots",
ODIR));
126 if(plot) plot->
graph[0]->SetTitle(TString::Format(
"%s : h+(black), hx(red))",wf.
name.Data()));
127 if(plot) plot->
canvas->SaveAs(TString::Format(
"%s/plots/%s_TIME.png",
ODIR,wf.
name.Data()));
133 plot->
graph[0]->GetHistogram()->GetXaxis()->SetRangeUser(32,1024);
134 if(plot) plot->
graph[0]->SetTitle(TString::Format(
"%s : h+(black), hx(red))",wf.
name.Data()));
135 if(plot) plot->
canvas->SaveAs(TString::Format(
"%s/plots/%s_FFT.png",
ODIR,wf.
name.Data()));
138 TString texiName=TString::Format(
"%s/html_index.texi",
ODIR);
140 if(!overwrite) gSystem->Exit(0);
143 out <<
"@c include texi macros" << endl;
144 out <<
"@include macros.texi" << endl;
145 out <<
"@include mathjax.texi" << endl << endl;
146 for(
int i=0;
i<wfList.size();
i++) {
147 out <<
"@center @txtfont{"<<wfList[
i].name<<
", red, h1}" << endl;
148 out <<
"@multitable @columnfractions .5 .5" << endl;
149 out <<
"@item @center @displayimage{../plots,"<<TString::Format(
"%s_TIME",wfList[
i].
name.Data())<<
",480}"<<endl;
150 out <<
"@tab @center @displayimage{../plots,"<<TString::Format(
"%s_FFT",wfList[
i].
name.Data())<<
",480}"<<endl;
151 out <<
"@end multitable" << endl;
159 TString exec_cmd = TString::Format(
"%s/cwb_mkhtml.csh %s wheader",
160 cwb_scripts.Data(),texiName.Data());
161 int ret=gSystem->Exec(exec_cmd);
163 cout <<
"Make_BRST_O1_LF.C : Error while executing cwb_mkhtml html_index.texi !!!" << endl;
166 cout << endl <<
"New html file created : " <<
ODIR <<
"/html_index/index.html" << endl << endl;
watplot * Draw(TString name, int id=0, TString polarization="hp", MDC_DRAW type=MDC_TIME, TString options="ALP", Color_t color=kBlack)
virtual void rate(double r)
std::vector< TGraph * > graph
vector< waveform > GetWaveformList()
virtual void start(double s)
mdcid AddWaveform(MDC_TYPE mdc_type, vector< mdcpar > par, TString uname="")
virtual size_t size() const
waveform GetWaveform(int ID, int id=0)