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 |
#define EPSILON 0.001 |
#define NMDC_MAX 64 |
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 | ( | ) |
pointers to detectors
Definition at line 48 of file cwb_report_sim.C.
gROOT ForceStyle | ( | ) |
if | ( | nIFO | = =1 | ) |
Definition at line 39 of file cwb_report_sim.C.
if | ( | pp_eff_vs_thr_quit | = ="TRUE" | ) |
Definition at line 332 of file cwb_report_sim.C.
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 | ) |
pp_eff_vs_thr_quit ToUpper | ( | ) |
FILE* evt_all = fopen(fname,"w") |
Definition at line 342 of file cwb_report_sim.C.
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.
char fname[4096] |
Definition at line 338 of file cwb_report_sim.C.
char IFO[NIFO_MAX][32] |
Definition at line 47 of file cwb_report_sim.C.
FILE* in_all = fopen(fname,"w") |
Definition at line 340 of file cwb_report_sim.C.
TString pp_eff_vs_thr_quit = CWB::Toolbox::getParameter(pp_eff_vs_thr,"--quit") |
Definition at line 330 of file cwb_report_sim.C.
float xfactors |
Definition at line 218 of file cwb_report_sim.C.
char xmdc_name[256] |
Definition at line 216 of file cwb_report_sim.C.
int xmdc_type |
Definition at line 217 of file cwb_report_sim.C.