16 cout <<
"You are going to produce the report for lag zero" << endl;
17 cout <<
"Do you want to procede ? (y/n) ";
20 }
while ((strcmp(answer,
"y")!=0)&&(strcmp(answer,
"n")!=0));
21 if (strcmp(answer,
"n")==0)
exit(1);
30 if(
pp_rho_min<0.) {cout <<
"Error : pp_rho_min must be >=0" << endl;
exit(1);}
36 if(
pp_rho_bin<=0) {cout <<
"Error : pp_rho_bin must be >0, check pp_rho_min, pp_rho_max, pp_drho" << endl;
exit(1);}
38 if(
TString(gSystem->Getenv(
"CWB_REPORT_PE"))!=
"") {
55 ifLow=200;ifHigh=2048;
57 #ifdef PP_1600_5000_HZ
59 ifLow=1600;ifHigh=5000;
90 else {cout <<
"Error : T_win must be >0" << endl;
exit(1);}}
92 sT_ifar.ReplaceAll(
".",
"d");
111 cout <<
"cwb_epparameters.C : Error - wrong pp_irho=" <<
pp_irho <<
" value "
112 <<
"must be 0/1" << endl;
116 cout <<
"cwb_epparameters.C : Error - wrong pp_irho=" <<
pp_inetcc <<
" value "
117 <<
"must be 0/1" << endl;
128 #ifdef PP_200_2048_HZ
131 #ifdef PP_1600_5000_HZ
139 cout<<
"Error in Frequency cuts : lowFCUT imust be < highFCUT"
154 sprintf(epp_cmd,
"freqLow = %f;",freqLow);
EXPORT(
double,freqLow,epp_cmd)
155 sprintf(epp_cmd,"freqHigh = %
f;",freqHigh);
EXPORT(
double,freqHigh,epp_cmd)
160 vector<waveSegment>
iseg; iseg.clear();
169 for(
int n=0;
n<oseg.size()-1;
n++) {
170 sprintf(freqs_label,
"%s_freq%d_%d",freqs_label,
int(oseg[
n].
stop),
int(oseg[
n+1].
start));
173 sfreqs_label.ReplaceAll(
".",
"d");
191 sprintf(
title,
"SIMULATION - Network %s [%4.1f-%4.1f]",xtitle,ifLow,ifHigh);
194 sprintf(
title,
"BACKGROUND NON ZERO LAG - Network %s [%4.1f-%4.1f]",xtitle,ifLow,ifHigh);
251 sprintf(ch2,
"netcc[%d]>%f&&lag[%d]==%d&&slag[%d]==%d",
269 for(
int i=0;i<
nFCUT;i++) {
289 for(
int i=0;i<nvdqf;i++) {
292 char veto_cat2[1024] =
"";
293 for(
int i=0;i<nvdqf;i++) {
295 if(strlen(veto_cat2)==0) {
302 cout <<
"veto_cat2 : " << veto_cat2 << endl;
303 sprintf(ch2,
"%s&&(%s)",ch2,veto_cat2);
306 #if defined (HVETO_VETO)
309 for(
int i=0;i<nvdqf;i++) {
326 #if defined (CAT3_VETO)
329 for(
int i=0;i<nvdqf;i++) {
345 #if defined (PEM_VETO)
349 #if defined (USER_VETO)
353 cout <<
"ch2 : " << ch2 << endl;
355 #if defined (HVETO_VETO) || defined (CAT3_VETO)
vector< waveSegment > oseg
#define EXPORT(TYPE, VAR, CMD)
vector< waveSegment > iseg
strcpy(RunLabel, RUN_LABEL)
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
char veto_not_vetoed[1024]
detectorParams detParms[4]