30 if(gSystem->Getenv(
"CWB_DOC_URL")!=NULL) {
31 cwb_doc_url=
TString(gSystem->Getenv(
"CWB_DOC_URL"));
43 TString NetCC,Rho,RhOut,RhAcor,svED,sPEN,sHRSS,sfLow,sfHigh,sfResample,spp_rho_min,spp_rho_max;
60 TTree* wave_tree = (TTree*)psp->Get(
"waveburst");
67 cout << fileout << endl;
69 if (!out.good()) {cout <<
"Error Opening File : " << fileout << endl;
exit(1);}
81 for(
int i=0;
i<epsList.size();
i++) {
83 ofile.ReplaceAll(
".eps",
".gif");
85 sprintf(cmd,
"convert %s %s",epsList[
i].Data(),ofile.Data());
88 sprintf(cmd,
"rm %s",epsList[
i].Data());
95 sprintf(fileliv,
"%s/live.txt",PP_DATA_PATH.Data());
96 cout << fileliv << endl;
98 double OLIVETIME =
GetLiveTime(fileliv,0,0,countlag);
99 double LIVETIME =
GetLiveTime(fileliv,-1,-1,countlag);
103 cout << OLIVETIME <<
" " << LIVETIME << endl;
113 bkgin.getline(istring,256);
114 if (!bkgin.good())
break;
117 ostringa.ReplaceAll(
"<!--VETO",
"");
118 ostringa.ReplaceAll(
"VETO-->",
"");
120 ostringa.ReplaceAll(
"<!--NOVETO",
"");
121 ostringa.ReplaceAll(
"VETONO-->",
"");
124 ostringa.ReplaceAll(
"<!--RHO_PF",
"");
125 ostringa.ReplaceAll(
"RHO_PF-->",
"");
127 ostringa.ReplaceAll(
"<!--RHO_SUBNET",
"");
128 ostringa.ReplaceAll(
"RHO_SUBNET-->",
"");
130 if(pp_rho_vs_rate_no_multiplicity) {
131 ostringa.ReplaceAll(
"<!--MULTI",
"");
132 ostringa.ReplaceAll(
"MULTI-->",
"");
134 ostringa.ReplaceAll(
"FLOW",sfLow);
135 ostringa.ReplaceAll(
"FHIGH",sfHigh);
136 ostringa.ReplaceAll(
"PP_RHO_MIN",spp_rho_min);
137 ostringa.ReplaceAll(
"PP_RHO_MAX",spp_rho_max);
138 ostringa.ReplaceAll(
"PP_DATA_DIR",PP_DATA_DIR);
139 ostringa.ReplaceAll(
"NETCC",NetCC);
140 ostringa.ReplaceAll(
"RHO",Rho);
141 ostringa.ReplaceAll(
"RH_OUT",RhOut);
142 ostringa.ReplaceAll(
"RH_ACOR",RhAcor);
143 ostringa.ReplaceAll(
"svED",svED);
145 ostringa.ReplaceAll(
"RATEBACKGROUND",Rate);
146 ostringa.ReplaceAll(
"SUBTITLE",
subtitle);
147 ostringa.ReplaceAll(
"TITLE",
title);
149 sprintf(ltime,
"zero lag : <font color=\"red\"> %.2f sec = %.2f days</font>",
150 OLIVETIME,OLIVETIME/86400.);
151 ostringa.ReplaceAll(
"LIVETIME",ltime);
152 sprintf(ltime,
"non-zero lags : <font color=\"red\"> %i lags - %.2f sec = %.2f days = %.1f years </fonts>",
153 countlag,LIVETIME,LIVETIME/(
float)86400.,LIVETIME/86400./365.);
154 ostringa.ReplaceAll(
"LIVETIM2",ltime);
156 if(cwb_doc_url!=
"") {
157 ostringa.ReplaceAll(
"<!--CWB_DOC_URL",
"");
158 ostringa.ReplaceAll(
"CWB_DOC_URL-->",
"");
159 ostringa.ReplaceAll(
"XCWB_DOC_URL",cwb_doc_url.Data());
162 ostringa.ReplaceAll(
"</html>",
"");
163 out << ostringa.Data() << endl;
174 sprintf(file_ev,
"%s/events.txt",PP_DATA_PATH.Data());
175 ifstream for_ev(file_ev);
176 char file_header[1024];
177 sprintf(file_header,
"%s/events_header.txt",PP_DATA_PATH.Data());
178 ofstream for_header(file_header);
179 char file_notsorted[1024];
180 sprintf(file_notsorted,
"%s/events_notsorted.txt",PP_DATA_PATH.Data());
181 ofstream for_notsorted(file_notsorted);
184 for_ev.getline(in_ev,100000);
185 if (!for_ev.good())
break;
187 if (sin_ev.Contains(
"#")) for_header << sin_ev.Data() << endl;
188 else for_notsorted << sin_ev.Data() << endl;
191 for_notsorted.close();
192 char file_sorted[1024];
193 sprintf(file_sorted,
"%s/events_sorted.txt",PP_DATA_PATH.Data());
194 sprintf(exec,
"sort -g -r --key=3 %s > %s",file_notsorted, file_sorted);
196 sprintf(exec,
"cat %s %s > %s/EVENTS.txt",file_header, file_sorted, PP_DATA_PATH.Data());
199 sprintf(file_notsorted,
"%s/events_notsorted.txt",PP_DATA_PATH.Data());
200 ifstream
f_ev(file_sorted);
203 float icc,
irho,
iacor,
ilag,
islag,
ilik,
ipen,
icHH,
ivHH,
ivED,
icc2,
icc3;
214 out <<
"<table border=0 cellpadding=2 class=\"datagrid\">" << endl;
215 out <<
"<tr align=\"center\">"<< endl;
216 out <<
"<td>ID</td>"<< endl;
217 if(bhveto) out <<
"<td>KW</td>"<< endl;
218 if(bcat3) out <<
"<td>cat3 </td>"<< endl;
219 out <<
"<td>rho["<<
pp_irho<<
"]</td>"<< endl;
220 out <<
"<td>cc["<<
pp_inetcc<<
"]</td>"<< endl;
221 if(
TString(analysis)==
"2G") out <<
"<td>subnet</td>"<< endl;
223 out <<
"<td>lag</td>"<< endl;
224 out <<
"<td>slag</td>"<< endl;
225 out <<
"<td>SNRnet</td>"<< endl;
226 if(
TString(analysis)==
"1G") out <<
"<td>pf</td>"<< endl;
227 if(
TString(analysis)==
"1G") out <<
"<td>vED</td>"<< endl;
228 out <<
"<td>freq</td>"<< endl;
229 out <<
"<td>bw</td>"<< endl;
230 out <<
"<td>dur</td>"<< endl;
231 out <<
"<td>size</td>"<< endl;
232 if(
optim==
true) out <<
"<td>res</td>"<< endl;
233 out <<
"<td>run</td>"<< endl;
234 for (
int nn=0;nn<
nIFO;nn++) out <<
"<td>GPS " << _ifo[nn] <<
"</td>"<< endl;
235 for (
int nn=0;nn<
nIFO;nn++) out <<
"<td>SNR " << _ifo[nn] <<
"</td>"<< endl;
242 out <<
"</tr>"<< endl;
247 f_ev>>pm>>c3>>irho>>icc>>icc2>>icc3>>iacor>>ilag>>islag>>ilik>>ipen>>icHH>>ifreq>>iband>>idur>>isize>>irate>>irun>>time0>>time1>>time2>>time3>>SNR0>>SNR1>>SNR2>>SNR3>>hrss0>>hrss1>>hrss2>>hrss3>>phi>>theta>>
psi;
249 f_ev>>pm>>c3>>irho>>icc>>icc2>>icc3>>iacor>>ilag>>islag>>ilik>>ipen>>icHH>>ifreq>>iband>>idur>>isize>>irate>>irun>>time0>>time1>>time2>>SNR0>>SNR1>>SNR2>>hrss0>>hrss1>>hrss2>>phi>>theta>>
psi;
251 f_ev>>pm>>c3>>irho>>icc>>icc2>>icc3>>iacor>>ilag>>islag>>ilik>>ipen>>icHH>>ifreq>>iband>>idur>>isize>>irate>>irun>>time0>>time1>>SNR0>>SNR1>>hrss0>>hrss1>>phi>>theta>>
psi;
253 if (!
f_ev.good())
break;
258 cout << ecount <<
" ";cout.flush();
259 out <<
"<tr align=\"center\">"<< endl;
265 start[0]=
GetStart(wave_tree, 0, irun, irho, time0, time1, analysis,
pp_irho);
268 for (
int nn=0; nn<
nIFO; nn++)
sprintf(namedir,
"%s%s",namedir,_ifo[nn]);
269 for (
int nn=0; nn<
nIFO; nn++) {
270 start[nn]=
GetStart(wave_tree, nn, irun, irho, time0, time1, analysis,
pp_irho);
271 sprintf(namedir,
"%s_%.3f",namedir,start[nn]);
274 sprintf(os,
"<td><a href=\"ced/ced_%s/%s\" target=\"_blank\">%i</a></td>",s_ced.Data(),namedir,
ecount);
277 sprintf(os,
"<td>%s</td>",
ppm.Remove(0,1).Data());
281 sprintf(os,
"<td>%s</td>",
cc3.Remove(0,1).Data());
284 sprintf(os,
"<td>%.2f</td>",irho);
286 sprintf(os,
"<td>%.2f</td>",icc);
289 sprintf(os,
"<td>%.2f</td>",icc2);
296 sprintf(os,
"<td>%i</td>",(
int)ilag);
298 sprintf(os,
"<td>%i</td>",(
int)islag);
300 sprintf(os,
"<td>%3.1f</td>",sqrt(ilik));
302 sprintf(os,
"<td>%.3f</td>",ipen);
303 if(
TString(analysis)==
"1G") out << os << endl;
304 sprintf(os,
"<td>%.3f</td>",icHH);
305 if(
TString(analysis)==
"1G") out << os << endl;
306 sprintf(os,
"<td><font color=\"black\">%i</font></td>",ifreq);
308 sprintf(os,
"<td>%i</td>",iband);
310 sprintf(os,
"<td>%.3f</td>",idur);
312 sprintf(os,
"<td>%i</td>",isize);
314 sprintf(os,
"<td>%i</td>",irate);
315 if(
optim==
true) out << os << endl;
316 sprintf(os,
"<td>%i</td>",irun);
318 sprintf(os,
"<td>%.2f</td>",time0);
321 sprintf(os,
"<td>%.2f</td>",time1);
325 sprintf(os,
"<td>%.2f</td>",time2);
329 sprintf(os,
"<td>%.2f</td>",time3);
332 sprintf(os,
"<td>%3.1f</td>",sqrt(SNR0));
335 sprintf(os,
"<td>%3.1f</td>",sqrt(SNR1));
339 sprintf(os,
"<td>%3.1f</td>",sqrt(SNR2));
343 sprintf(os,
"<td>%3.1f</td>",sqrt(SNR3));
366 out <<
"</tr>" << endl;
372 out <<
"</table>" << endl;
373 out <<
"<p>" << endl;
376 cout << exec << endl;
379 sprintf(file_notsorted,
"%s/events_notsorted.txt",PP_DATA_PATH.Data());
380 sprintf(exec,
"rm %s",file_notsorted);
381 cout << exec << endl;
383 sprintf(exec,
"rm %s",file_header);
384 cout << exec << endl;
387 out <<
"</table>" << endl;
388 out <<
"<p>" << endl;
void Export(TString fname="")
sprintf(netdir,"%s/%s", pp_dir, pp_data_dir)
char html_body_prod_template[1024]
char channelNamesRaw[NIFO_MAX][50]
vector< TString > epsList
cout<< "baudline_FFL : "<< baudline_FFL<< endl;ofstream out;out.open(baudline_FFL, ios::out);if(!out.good()){cout<< "Error Opening File : "<< baudline_FFL<< endl;exit(1);}ifstream in;in.open(frFiles[ifoID], ios::in);if(!in.good()){cout<< "Error Opening File : "<< frFiles[ifoID]<< endl;exit(1);}TString pfile_path="";char istring[1024];while(1){in > istring
void Import(TString umacro="")
ifstream f_ev(file_sorted)
strcpy(RunLabel, RUN_LABEL)
detectorParams detParms[4]
void SetSingleDetectorMode()