58 pp_pe_mdc.ReplaceAll(
"/",
"");
59 if(!pp_pe_mdc.IsDigit() && pp_pe_mdc!=
"-1") {
60 cout <<
"cwb_report_pe.C : Error - bad option '--mdc type' in pp_pe parameter" << endl;
61 cout <<
" mdc type must be a integer number [-1:N] or i/j/.../k=[1:N]" << endl;
69 for(
int i=0;
i<token->GetEntries();
i++) {
70 int itype = ((TObjString*)token->At(
i))->
GetString().Atoi();
71 mtype.push_back(itype);
76 for(
int i=0;
i<mtype.size();
i++) {
77 if(mtype[
i]==-1 || mtype[
i]==0) {
78 cout <<
"cwb_report_pe.C : Error - bad option '--mdc type' in pp_pe parameter" << endl;
79 cout <<
" multiple mdc types must be in the range [1:N]" << endl;
86 cout <<
"cwb_report_pe.C : Error - bad option '--mdc type' in pp_pe parameter" << endl;
87 cout <<
" mdc type must be in the range [-1:N] " << endl;
105 imdc_name,imdc_fcentral,imdc_fbandwidth);
107 cout <<
"cwb_report_pe.C : Error - no injection - terminated" << endl;
111 for(
int i=0;
i<mtype.size();
i++) {
113 if(mtype[
i]==(imdc_type[
j]+1)) {
114 vtype.push_back(imdc_type[
j]+1);
115 vname.push_back(imdc_name[
j]);
121 vtype.push_back(imdc_type[
j]+1);
122 vname.push_back(imdc_name[
j]);
127 vname.push_back(
"ALL");
130 Color_t
colors[16] = {2, 3, 6, 4, 8, 43, 7, 8, 4, 5, 2, 43, 1, 3, 2, 1};
136 cout << listFile << endl;
138 if(!out.good()) {cout <<
"cwb_report_pe.C : Error Opening File : " << listFile << endl;
exit(1);}
139 for(
int i=0;
i<vtype.size();
i++) {
140 char wave_file_name[1024];
142 out << wave_file_name <<
"\t" << vname[
i] <<
"\t" << vtype[
i]
143 <<
"\t" << (
int)colors[
i%16] <<
"\t" << 0 << endl;
148 gROOT->LoadMacro(gSystem->ExpandPathName(
"$HOME_CWB/macros/DrawWRC.C"));
149 gROOT->LoadMacro(gSystem->ExpandPathName(
"$HOME_CWB/macros/DrawRECvsINJ.C"));
150 gROOT->LoadMacro(gSystem->ExpandPathName(
"$HOME_CWB/macros/DrawSkyDistributionPRC.C"));
151 gROOT->LoadMacro(gSystem->ExpandPathName(
"$HOME_CWB/macros/DrawSearchAreaPRC.C"));
152 gROOT->LoadMacro(gSystem->ExpandPathName(
"$HOME_CWB/macros/DrawCosOmegaPRC.C"));
153 gROOT->LoadMacro(gSystem->ExpandPathName(
"$HOME_CWB/macros/DrawCoverageVsPercentagePRC.C"));
154 gROOT->LoadMacro(gSystem->ExpandPathName(
"$HOME_CWB/macros/DrawMedianPRCvsSNR.C"));
160 cout <<
"pe polarization options : " <<
parPlugin << endl;
162 for(
int j=0;
j<token->GetEntries();
j++) {
164 TObjString* tok = (TObjString*)token->At(
j);
165 TString stok = tok->GetString();
167 if(stok.Contains(
"--pe_polarization=")) {
169 polarization.Remove(0,polarization.Last(
'=')+1);
170 polarization.ToUpper();;
205 for(
int n=1;n<13;n++) {
206 char gtype[16];
sprintf(gtype,
"mch%d",n);
sprintf(listFile,"%s/DrawMedianPRCvsSNR.lst", netdir)
void Export(TString fname="")
DrawCosOmegaPRC(data_label, netdir, mlabel, nIFO, T_win, pp_inetcc, T_cor, pp_irho, T_cut, T_vED, T_pen, T_ifar)
DrawWRC("nre", data_label, netdir, mlabel, nIFO, T_win, pp_inetcc, T_cor, pp_irho, T_cut, T_vED, T_pen, T_ifar)
void Import(TString umacro="")
size_t imdc_iset[NMDC_MAX]
DrawCoverageVsPercentagePRC("prc", data_label, netdir, mlabel, nIFO, T_win, pp_inetcc, T_cor, pp_irho, T_cut, T_vED, T_pen, T_ifar)
size_t imdc_index[NMDC_MAX]
DrawSkyDistributionPRC(data_label, netdir, mlabel, IFO, detParms, T_win, pp_inetcc, T_cor, pp_irho, T_cut, T_vED, T_pen, T_ifar, true, polarization, true)
char imdc_name[NMDC_MAX][128]
DrawRECvsINJ("snr", data_label, netdir, mlabel, nIFO, T_win, pp_inetcc, T_cor, pp_irho, T_cut, T_vED, T_pen, T_ifar)
DrawSearchAreaPRC(data_label, netdir, mlabel, nIFO, T_win, pp_inetcc, T_cor, pp_irho, T_cut, T_vED, T_pen, T_ifar)
DrawMedianPRCvsSNR(listFile,"","MEDIAN50", TString(netdir)+"/median50_vs_snr.png", T_win, pp_inetcc, T_cor, pp_irho, T_cut, T_vED, T_pen, T_ifar)
size_t imdc_type[NMDC_MAX]
detectorParams detParms[4]
double imdc_fcentral[NMDC_MAX]
char imdc_set[NMDC_MAX][128]
void SetSingleDetectorMode()
double imdc_fbandwidth[NMDC_MAX]