Logo coherent WaveBurst  
Library Reference Guide
Logo
Macros | Functions | Variables
cwb_report_sim.C File Reference

Macros

#define EPSILON   0.001
 
#define NMDC_MAX   64
 

Functions

cout<< fname<< endl;FILE *fp=fopen(fname,"w");if(fp==NULL) { cout<< "cwb_report_sim.C - Error opening file : "<< fname<< endl;gSystem-> Exit (1)
 
 exit (0)
 
 fclose (fp)
 
 fclose (in_all)
 
 fclose (evt_all)
 
 for (int n=0;n< nIFO;n++)
 pointers to detectors More...
 
gROOT ForceStyle ()
 
 if (nIFO==1)
 
 if (pp_eff_vs_thr_quit=="TRUE")
 
gStyle SetCanvasColor (kWhite)
 
gStyle SetFrameBorderMode (0)
 
gStyle SetLabelFont (42,"X")
 
gStyle SetLabelFont (42,"Y")
 
gStyle SetLabelOffset (0.014,"X")
 
gStyle SetLabelOffset (0.010,"Y")
 
gStyle SetLabelSize (0.03,"X")
 
gStyle SetLabelSize (0.03,"Y")
 
gStyle SetNumberContours (256)
 
gStyle SetStatBorderSize (1)
 
gStyle SetStatColor (0)
 
gStyle SetStatH (0.160)
 
gStyle SetStatW (0.200)
 
gStyle SetStatX (0.878)
 
gStyle SetStatY (0.918)
 
gStyle SetTextFont (12)
 
gStyle SetTitleColor (kBlue,"D")
 
gStyle SetTitleFillColor (kWhite)
 
gStyle SetTitleFont (42,"X")
 
gStyle SetTitleFont (42,"Y")
 
gStyle SetTitleFont (12,"D")
 
gStyle SetTitleH (0.050)
 
gStyle SetTitleOffset (1.4,"X")
 
gStyle SetTitleOffset (1.2,"Y")
 
gStyle SetTitleW (0.95)
 
gStyle SetTitleY (0.98)
 
 sprintf (fname,"%s/eff_%d_threshold_factors.txt", netdir, pp_eff)
 
 sprintf (fname,"%s/fit_parameters_ALL.txt", netdir)
 
 sprintf (fname,"%s/evt_parameters_ALL.txt", netdir)
 
pp_eff_vs_thr_quit ToUpper ()
 

Variables

FILE * evt_all = fopen(fname,"w")
 
cout<<"cwb_report_sim.C starts..."<< endl;Color_t colors[NMDC_MAX]={ 6, 3, 2, 8, 43, 7, 8, 4, 4, 2, 43, 1, 3, 1, 6, 7, 6, 3, 2, 8, 43, 7, 8, 4, 4, 2, 43, 1, 3, 1, 6, 7, 6, 3, 2, 8, 43, 7, 8, 4, 4, 2, 43, 1, 3, 1, 6, 7, 6, 3, 2, 8, 43, 7, 8, 4, 4, 2, 43, 1, 3, 1, 6, 7 };Style_t markers[NMDC_MAX]={20, 21, 22, 23, 29, 29, 21, 22, 20, 21, 22, 23, 29, 20, 21, 22, 21, 22, 23, 29, 29, 21, 22, 20, 21, 22, 23, 29, 20, 21, 22, 20, 21, 22, 23, 29, 29, 21, 22, 20, 21, 22, 23, 29, 20, 21, 22, 20, 21, 22, 23, 29, 29, 21, 22, 20, 21, 22, 23, 29, 20, 21, 22, 20 };TH2F *HH[NIFO_MAX][NMDC_MAX];TH2F *NRE[NIFO_MAX][NMDC_MAX];TH1F *hT[NMDC_MAX];TH1F *hF[NMDC_MAX];TGraphErrors *EFF[NMDC_MAX];for(int i=0;i< NMDC_MAX;i++) { for(int j=0;j< nIFO;j++) HH[j][i]=NULL;for(int j=0;j< nIFO;j++) NRE[j][i]=NULL;hT[i]=NULL;hF[i]=NULL;EFF[i]=NULL;} T_ifar *=(365.*24.*3600.);dqfile *XDQF=new dqfile[nVDQF];char imdc_set[NMDC_MAX][128];size_t imdc_type[NMDC_MAX];char imdc_name[NMDC_MAX][128];double imdc_flow[NMDC_MAX];double imdc_fhigh[NMDC_MAX];double imdc_fcentral[NMDC_MAX];double imdc_fbandwidth[NMDC_MAX];size_t imdc_index[NMDC_MAX];size_t imdc_iset[NMDC_MAX];size_t imdc_tset[NMDC_MAX];int ninj=ReadInjType(mdc_inj_file, NMDC_MAX, imdc_set, imdc_type, imdc_name, imdc_fcentral, imdc_fbandwidth);if(ninj==0) { cout<< endl;cout<< "cwb_report_sim.C : Error - no injection types or bad format"<< endl;cout<< "mdc injection file list : "<< mdc_inj_file<< endl;cout<< "format must be : mdc_set mdc_type mdc_name mdc_fcentral mdc_fbandwidth"<< endl;cout<< "process terminated !!!"<< endl;exit(1);} for(int i=0;i< NMDC_MAX;i++) { imdc_flow[i]=imdc_fcentral[i]-imdc_fbandwidth[i];if(imdc_flow[i]< fLow) imdc_flow[i]=fLow;if(imdc_flow[i]>fHigh) imdc_flow[i]=fLow;imdc_fhigh[i]=imdc_fcentral[i]+imdc_fbandwidth[i];if(imdc_fhigh[i]< fLow) imdc_fhigh[i]=fHigh;if(imdc_fhigh[i]>fHigh) imdc_fhigh[i]=fHigh;} for(int i=0;i< NMDC_MAX;i++) imdc_index[i]=NMDC_MAX;for(int j=0;j< ninj;j++) imdc_index[imdc_type[j]]=j;for(int j=0;j< ninj;j++) { cout<< j<< " "<< imdc_index[j]<< endl;if(imdc_index[j]>=ninj) { cout<< endl;cout<< "cwb_report_sim.C : Error - mdc type must be < num max type inj = "<< ninj<< endl;cout<< "check mdc_type injection file list : "<< mdc_inj_file<< endl;cout<< "format : mdc_set mdc_type mdc_name mdc_fcentral mdc_fbandwidth"<< endl;cout<< "process terminated !!!"<< endl;exit(1);} } TString *imdc_set_name=new TString[ninj];int nset=0;for(int i=0;i< ninj;i++) { bool bnew=true;for(int j=0;j< nset;j++) if(imdc_set[i]==imdc_set_name[j]) bnew=false;if(bnew) imdc_set_name[nset++]=imdc_set[i];} cout<< "nset : "<< nset<< endl;for(int i=0;i< nset;i++) { for(int j=0;j< ninj;j++) if(imdc_set[j]==imdc_set_name[i]) imdc_iset[j]=i;} for(int i=0;i< ninj;i++) { for(int j=0;j< nset;j++) if(imdc_set[i]==imdc_set_name[j]) imdc_tset[imdc_type[i]]=j;} for(int i=0;i< ninj;i++) { cout<< i<< " "<< imdc_name[i]<< " "<< imdc_set[i]<< " "<< imdc_iset[i]<< " "<< imdc_tset[i]<< " "<< imdc_type[i]<< " "<< imdc_index[i]<< endl;} TString pp_eff_vs_thr_mode=CWB::Toolbox::getParameter(pp_eff_vs_thr,"--mode");if((pp_eff_vs_thr_mode!="DISABLED")&&(pp_eff_vs_thr_mode!="disabled")) { TString pp_eff_vs_thr_eff=CWB::Toolbox::getParameter(pp_eff_vs_thr,"--eff");int pp_eff=50;if(pp_eff_vs_thr_eff.IsFloat()) pp_eff=pp_eff_vs_thr_eff.Atoi();if(pp_eff< 0) pp_eff=0;if(pp_eff >100) pp_eff=100;TChain sim("waveburst");TChain mdc("mdc");sim.Add(sim_file_name);mdc.Add(mdc_file_name);char fname[1024];wavearray< double > factorXX(ninj);factorXX=0;if((pp_eff_vs_thr_mode!="AUTO")&&(pp_eff_vs_thr_mode!="auto")) { sprintf(fname,"%s", pp_eff_vs_thr_mode.Data());FILE *fp=fopen(fname,"r");if(fp==NULL) { cout<< "cwb_report_sim.C - Error opening file : "<< fname<< endl;gSystem-> Exit (1)
 
char fname [4096]
 
char IFO [NIFO_MAX][32]
 
FILE * in_all = fopen(fname,"w")
 
TString pp_eff_vs_thr_quit = CWB::Toolbox::getParameter(pp_eff_vs_thr,"--quit")
 
float xfactors
 
char xmdc_name [256]
 
int xmdc_type
 

Macro Definition Documentation

◆ EPSILON

#define EPSILON   0.001

◆ NMDC_MAX

#define NMDC_MAX   64

Function Documentation

◆ Exit()

cout<< fname << endl; FILE* fp = fopen(fname,"w"); if(fp==NULL) { cout << "cwb_report_sim.C - Error opening file : " << fname << endl; gSystem-> Exit ( )

◆ exit()

exit ( )

◆ fclose() [1/3]

fclose ( fp  )

◆ fclose() [2/3]

fclose ( in_all  )

◆ fclose() [3/3]

fclose ( evt_all  )

◆ for()

for ( )

pointers to detectors

Definition at line 48 of file cwb_report_sim.C.

◆ ForceStyle()

gROOT ForceStyle ( )

◆ if() [1/2]

if ( nIFO  = =1)

Definition at line 39 of file cwb_report_sim.C.

◆ if() [2/2]

if ( pp_eff_vs_thr_quit  = ="TRUE")

Definition at line 332 of file cwb_report_sim.C.

◆ SetCanvasColor()

gStyle SetCanvasColor ( kWhite  )

◆ SetFrameBorderMode()

gStyle SetFrameBorderMode ( )

◆ SetLabelFont() [1/2]

gStyle SetLabelFont ( 42  ,
"X"   
)

◆ SetLabelFont() [2/2]

gStyle SetLabelFont ( 42  ,
"Y"   
)

◆ SetLabelOffset() [1/2]

gStyle SetLabelOffset ( 0.  014,
"X"   
)

◆ SetLabelOffset() [2/2]

gStyle SetLabelOffset ( 0.  010,
"Y"   
)

◆ SetLabelSize() [1/2]

gStyle SetLabelSize ( 0.  03,
"X"   
)

◆ SetLabelSize() [2/2]

gStyle SetLabelSize ( 0.  03,
"Y"   
)

◆ SetNumberContours()

gStyle SetNumberContours ( 256  )

◆ SetStatBorderSize()

gStyle SetStatBorderSize ( )

◆ SetStatColor()

gStyle SetStatColor ( )

◆ SetStatH()

gStyle SetStatH ( 0.  160)

◆ SetStatW()

gStyle SetStatW ( 0.  200)

◆ SetStatX()

gStyle SetStatX ( 0.  878)

◆ SetStatY()

gStyle SetStatY ( 0.  918)

◆ SetTextFont()

gStyle SetTextFont ( 12  )

◆ SetTitleColor()

gStyle SetTitleColor ( kBlue  ,
"D"   
)

◆ SetTitleFillColor()

gStyle SetTitleFillColor ( kWhite  )

◆ SetTitleFont() [1/3]

gStyle SetTitleFont ( 42  ,
"X"   
)

◆ SetTitleFont() [2/3]

gStyle SetTitleFont ( 42  ,
"Y"   
)

◆ SetTitleFont() [3/3]

gStyle SetTitleFont ( 12  ,
"D"   
)

◆ SetTitleH()

gStyle SetTitleH ( 0.  050)

◆ SetTitleOffset() [1/2]

gStyle SetTitleOffset ( 1.  4,
"X"   
)

◆ SetTitleOffset() [2/2]

gStyle SetTitleOffset ( 1.  2,
"Y"   
)

◆ SetTitleW()

gStyle SetTitleW ( 0.  95)

◆ SetTitleY()

gStyle SetTitleY ( 0.  98)

◆ sprintf() [1/3]

sprintf ( fname  ,
"%s/eff_%d_threshold_factors.txt"  ,
netdir  ,
pp_eff   
)

◆ sprintf() [2/3]

sprintf ( fname  ,
"%s/fit_parameters_ALL.txt"  ,
netdir   
)

◆ sprintf() [3/3]

sprintf ( fname  ,
"%s/evt_parameters_ALL.txt"  ,
netdir   
)

◆ ToUpper()

pp_eff_vs_thr_quit ToUpper ( )

Variable Documentation

◆ evt_all

FILE* evt_all = fopen(fname,"w")

Definition at line 342 of file cwb_report_sim.C.

◆ Exit

cout<<"cwb_report_sim.C starts..."<<endl; Color_t colors[NMDC_MAX] = { 6, 3, 2, 8,43, 7, 8, 4, 4, 2,43, 1, 3, 1, 6, 7, 6, 3, 2, 8,43, 7, 8, 4, 4, 2,43, 1, 3, 1, 6, 7, 6, 3, 2, 8,43, 7, 8, 4, 4, 2,43, 1, 3, 1, 6, 7, 6, 3, 2, 8,43, 7, 8, 4, 4, 2,43, 1, 3, 1, 6, 7 }; Style_t markers[NMDC_MAX]= {20,21,22,23,29,29,21,22,20,21,22,23,29,20,21,22, 21,22,23,29,29,21,22,20,21,22,23,29,20,21,22,20, 21,22,23,29,29,21,22,20,21,22,23,29,20,21,22,20, 21,22,23,29,29,21,22,20,21,22,23,29,20,21,22,20 }; TH2F* HH[NIFO_MAX][NMDC_MAX]; TH2F* NRE[NIFO_MAX][NMDC_MAX]; TH1F* hT[NMDC_MAX]; TH1F* hF[NMDC_MAX]; TGraphErrors* EFF[NMDC_MAX]; for(int i=0;i<NMDC_MAX;i++) { for(int j=0;j<nIFO;j++) HH[j][i]=NULL; for(int j=0;j<nIFO;j++) NRE[j][i]=NULL; hT[i]=NULL; hF[i]=NULL; EFF[i]=NULL; } T_ifar*=(365.*24.*3600.); dqfile* XDQF = new dqfile[nVDQF]; char imdc_set[NMDC_MAX][128]; size_t imdc_type[NMDC_MAX]; char imdc_name[NMDC_MAX][128]; double imdc_flow[NMDC_MAX]; double imdc_fhigh[NMDC_MAX]; double imdc_fcentral[NMDC_MAX]; double imdc_fbandwidth[NMDC_MAX]; size_t imdc_index[NMDC_MAX]; size_t imdc_iset[NMDC_MAX]; size_t imdc_tset[NMDC_MAX]; int ninj=ReadInjType(mdc_inj_file,NMDC_MAX,imdc_set,imdc_type, imdc_name,imdc_fcentral,imdc_fbandwidth); if(ninj==0) { cout << endl; cout << "cwb_report_sim.C : Error - no injection types or bad format" << endl; cout << "mdc injection file list : " << mdc_inj_file << endl; cout << "format must be : mdc_set mdc_type mdc_name mdc_fcentral mdc_fbandwidth" << endl; cout << "process terminated !!!" << endl; exit(1); } for(int i=0; i<NMDC_MAX; i++) { imdc_flow[i] = imdc_fcentral[i]-imdc_fbandwidth[i]; if(imdc_flow[i]<fLow) imdc_flow[i]=fLow; if(imdc_flow[i]>fHigh) imdc_flow[i]=fLow; imdc_fhigh[i] = imdc_fcentral[i]+imdc_fbandwidth[i]; if(imdc_fhigh[i]<fLow) imdc_fhigh[i]=fHigh; if(imdc_fhigh[i]>fHigh) imdc_fhigh[i]=fHigh; } for(int i=0; i<NMDC_MAX; i++) imdc_index[i] = NMDC_MAX; for(int j=0;j<ninj;j++) imdc_index[imdc_type[j]]=j; for(int j=0;j<ninj;j++) { cout << j << " " << imdc_index[j] << endl; if(imdc_index[j]>=ninj) { cout << endl; cout << "cwb_report_sim.C : Error - mdc type must be < num max type inj = " << ninj << endl; cout << "check mdc_type injection file list : " << mdc_inj_file << endl; cout << "format : mdc_set mdc_type mdc_name mdc_fcentral mdc_fbandwidth" << endl; cout << "process terminated !!!" << endl; exit(1); } } TString* imdc_set_name = new TString[ninj]; int nset=0; for(int i=0;i<ninj;i++) { bool bnew=true; for(int j=0;j<nset;j++) if(imdc_set[i]==imdc_set_name[j]) bnew=false; if(bnew) imdc_set_name[nset++]=imdc_set[i]; } cout << "nset : " << nset << endl; for(int i=0;i<nset;i++) { for(int j=0;j<ninj;j++) if(imdc_set[j]==imdc_set_name[i]) imdc_iset[j]=i; } for(int i=0;i<ninj;i++) { for(int j=0;j<nset;j++) if(imdc_set[i]==imdc_set_name[j]) imdc_tset[imdc_type[i]]=j; } for(int i=0;i<ninj;i++) { cout << i << " " << imdc_name[i] << " " << imdc_set[i] << " " << imdc_iset[i] << " " << imdc_tset[i] << " " << imdc_type[i] << " " << imdc_index[i] << endl; } TString pp_eff_vs_thr_mode = CWB::Toolbox::getParameter(pp_eff_vs_thr,"--mode"); if((pp_eff_vs_thr_mode!="DISABLED")&&(pp_eff_vs_thr_mode!="disabled")) { TString pp_eff_vs_thr_eff = CWB::Toolbox::getParameter(pp_eff_vs_thr,"--eff"); int pp_eff = 50; if(pp_eff_vs_thr_eff.IsFloat()) pp_eff = pp_eff_vs_thr_eff.Atoi(); if(pp_eff<0) pp_eff=0; if(pp_eff>100) pp_eff=100; TChain sim("waveburst"); TChain mdc("mdc"); sim.Add(sim_file_name); mdc.Add(mdc_file_name); char fname[1024]; wavearray<double> factorXX(ninj); factorXX=0; if((pp_eff_vs_thr_mode!="AUTO")&&(pp_eff_vs_thr_mode!="auto")) { sprintf(fname,"%s",pp_eff_vs_thr_mode.Data()); FILE* fp = fopen(fname,"r"); if(fp==NULL) { cout << "cwb_report_sim.C - Error opening file : " << fname << endl; gSystem-> Exit(1)

Definition at line 214 of file cwb_report_sim.C.

◆ fname

char fname[4096]

Definition at line 338 of file cwb_report_sim.C.

◆ IFO

char IFO[NIFO_MAX][32]

Definition at line 47 of file cwb_report_sim.C.

◆ in_all

FILE* in_all = fopen(fname,"w")

Definition at line 340 of file cwb_report_sim.C.

◆ pp_eff_vs_thr_quit

TString pp_eff_vs_thr_quit = CWB::Toolbox::getParameter(pp_eff_vs_thr,"--quit")

Definition at line 330 of file cwb_report_sim.C.

◆ xfactors

float xfactors

Definition at line 218 of file cwb_report_sim.C.

◆ xmdc_name

char xmdc_name[256]

Definition at line 216 of file cwb_report_sim.C.

◆ xmdc_type

int xmdc_type

Definition at line 217 of file cwb_report_sim.C.