Logo Coherent WaveBurst  
Reference Guide 		Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C
Go to the documentation of this file.
1 
2 #define nFAR_MAX 32
3 
4 int nFAR;
5 TString FAR_PLOT_NAME;
6 TString FAR_TITLE;
7 TString FAR_NAME[nFAR_MAX];
8 int LINE_COLOR[nFAR_MAX];
9 int LINE_STYLE[nFAR_MAX];
10 int LINE_MARKER[nFAR_MAX];
11 double OBS_TIME[nFAR_MAX];
12 TString FAR_PATH[nFAR_MAX];
13 double ZERO_LAG_TIME[nFAR_MAX];
14 
15 int readParameters(TString fname, wavearray<double>& RHO, wavearray<double>& eRHO, wavearray<double>& FAR, wavearray<double>& eFAR);
16 int writeParameters(TString fname, wavearray<double>& RHO, wavearray<double>& FAR);
17 void PlotFAR(int nfar, TString far_plot_name, double rho_min=-1, double rho_max=-1);
18 
19 void
20 Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a() {
21 
22 
23  Style_t markers[32]= {20,21,22,23,29,29,21,22,20,21,22,23,29,20,21,22,
24  21,22,23,29,29,21,22,20,21,22,23,29,20,21,22,20 };
25 
26  Color_t colors[32] = { 8, 0, 2, 8,43, 7, 8, 4, 4, 2,43, 1, 3, 1, 6, 7,
27  6, 3, 2, 8,43, 7, 8, 4, 4, 2,43, 1, 3, 1, 6, 7 };
28 
29  for(int i=0;i<nFAR_MAX;i++) {
30  LINE_COLOR[i] = colors[i];
31  LINE_STYLE[i] = 1;
32  LINE_MARKER[i] = markers[i];
33  ZERO_LAG_TIME[i] = -1;
34  OBS_TIME[i] = 1;
35  }
36 
37  FAR_TITLE="";
38  FAR_PLOT_NAME="";
39 
40  //nFAR = 2;
41  //FAR_PATH[0]="chirp_low_bkg_high_bins.txt";
42  //FAR_PATH[1]="chirp_high_bkg_high_bins.txt";
43 
44  //FAR_PATH[0]="chirp_full_bkg_high_bins.txt";
45  //FAR_PATH[0]="chirp_full_bkg_high_bins_0d01.txt";
46  //FAR_PATH[0]="chirp_full_ext_bkg_high_bins_0d01.txt";
47 /*
48  nFAR = 3;
49  FAR_PATH[0]="unmodeled_full_ext_bkg_high_bins_0d01_M5.txt";
50  FAR_PATH[1]="constrained_full_ext_bkg_high_bins_0d01_M5.txt";
51  FAR_PATH[2]="chirp_full_ext_bkg_high_bins_0d01_M5.txt";
52 */
53 /*
54  nFAR = 3;
55  FAR_PATH[0]="unmodeled_full_ext_bkg_high_bins_0d01_C01_M1.txt";
56  FAR_PATH[1]="constrained_full_ext_bkg_high_bins_0d01_C01_M1.txt";
57  FAR_PATH[2]="chirp_full_ext_bkg_high_bins_0d01_C01_M1.txt";
58 */
59 
60 // ===========================================================================
61 // C01+C00
62 // ===========================================================================
63 
64 /*
65  #define MODE "inclusive"
66  nFAR = 6;
67  FAR_PATH[0]="FARvsRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_unmodeled.txt";
68  FAR_PATH[1]="FARvsRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_constrained.txt";
69  FAR_PATH[2]="FARvsRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_chirp.txt";
70  FAR_PATH[3]="FARvsLowRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_unmodeled.txt";
71  FAR_PATH[4]="FARvsLowRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_constrained.txt";
72  FAR_PATH[5]="FARvsLowRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_chirp.txt";
73 */
74 /*
75  #define MODE "exclusive"
76  nFAR = 6;
77  FAR_PATH[0]="FARvsRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_eunmodeled.txt";
78  FAR_PATH[1]="FARvsRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_econstrained.txt";
79  FAR_PATH[2]="FARvsRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_chirp.txt";
80  FAR_PATH[3]="FARvsLowRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_eunmodeled.txt";
81  FAR_PATH[4]="FARvsLowRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_econstrained.txt";
82  FAR_PATH[5]="FARvsLowRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_chirp.txt";
83 */
84 
85 // ===========================================================================
86 // C01+C01
87 // ===========================================================================
88 
89 /*
90  #define MODE "inclusive"
91  nFAR = 6;
92  FAR_PATH[0]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_unmodeled.txt";
93  FAR_PATH[1]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_constrained.txt";
94  FAR_PATH[2]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_chirp.txt";
95  FAR_PATH[3]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_unmodeled.txt";
96  FAR_PATH[4]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_constrained.txt";
97  FAR_PATH[5]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_chirp.txt";
98 */
99 /*
100  #define MODE "exclusive"
101  nFAR = 6;
102  FAR_PATH[0]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_eunmodeled.txt";
103  FAR_PATH[1]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_econstrained.txt";
104  FAR_PATH[2]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_chirp.txt";
105  FAR_PATH[3]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_eunmodeled.txt";
106  FAR_PATH[4]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_econstrained.txt";
107  FAR_PATH[5]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_chirp.txt";
108 */
109 
110 // ===========================================================================
111 // C01+C01 h1_glitch_removed
112 // ===========================================================================
113 
114 /*
115  #define MODE "inclusive"
116  nFAR = 6;
117  FAR_PATH[0]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_h1_glitch_removed_unmodeled.txt";
118  FAR_PATH[1]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_h1_glitch_removed_constrained.txt";
119  FAR_PATH[2]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_h1_glitch_removed_chirp.txt";
120  FAR_PATH[3]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_unmodeled.txt";
121  FAR_PATH[4]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_constrained.txt";
122  FAR_PATH[5]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_chirp.txt";
123 */
124 
125  #define MODE "exclusive"
126  nFAR = 6;
127  FAR_PATH[0]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_h1_glitch_removed_eunmodeled.txt";
128  FAR_PATH[1]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_h1_glitch_removed_econstrained.txt";
129  FAR_PATH[2]="FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_h1_glitch_removed_chirp.txt";
130  FAR_PATH[3]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_eunmodeled.txt";
131  FAR_PATH[4]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_econstrained.txt";
132  FAR_PATH[5]="FARvsLowRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a_chirp.txt";
133 
134 
135  PlotFAR(nFAR, FAR_PLOT_NAME, -1, -1);
136 
137 }
138 
139 void PlotFAR(int nfar, TString far_plot_name, double rho_min, double rho_max) {
140 
141  // create plots
142  gStyle->SetFrameBorderMode(0); // remove the red box around canvas
143  gROOT->ForceStyle();
144 
145  gStyle->SetTitleFont(72);
146  gStyle->SetMarkerColor(50);
147  gStyle->SetLineColor(kWhite);
148  gStyle->SetTitleW(0.98);
149  gStyle->SetTitleH(0.05);
150  gStyle->SetTitleY(0.98);
151  gStyle->SetFillColor(kWhite);
152  gStyle->SetLineColor(kWhite);
153  gStyle->SetTitleFont(12,"D");
154 
155  TCanvas *canvas = new TCanvas("roc", "roc", 300,40, 800, 500);
156  canvas->Clear();
157  canvas->ToggleEventStatus();
158  canvas->SetLogy();
159  canvas->SetGridx();
160  canvas->SetGridy();
161  canvas->SetFillColor(kWhite);
162 
163  double far_min=1e10;
164  double rmin=1e20;
165  double rmax=0;
166  char fname[1024];
167  int nGR=0;
168  int line_style[3*nFAR_MAX];
169  int line_marker[3*nFAR_MAX];
170  int line_color[3*nFAR_MAX];
171  TString FAR_name[3*nFAR_MAX];
172  bool sigma_lines[3*nFAR_MAX];
173  TGraphErrors* gr[3*nFAR_MAX];
174  for(int n=0;n<nFAR;n++) {
175  cout << n << " OBS_TIME : " << OBS_TIME[n] << endl;
176  if(OBS_TIME[n]>0) if(far_min>1./OBS_TIME[n]) far_min=1./OBS_TIME[n];
177  vector<double> rho_far;
178  vector<double> far;
179  // read FAR
180  sprintf(fname,"%s",FAR_PATH[n].Data());
181  cout << "rate_threshold : " << fname << endl;
182  wavearray<double> RHO;
183  wavearray<double> FAR;
184  wavearray<double> eRHO;
185  wavearray<double> eFAR;
186  int far_size = readParameters(fname, RHO, eRHO, FAR, eFAR);
187 #if ( MODE == "exclusive")
188  TString ofname = fname;
189  ofname.ReplaceAll(".txt","_comp_paper.txt");
190  writeParameters(ofname, RHO, FAR);
191 #endif
192  if(rmin<RHO[0]) rmin=RHO[0];
193  if(rmax>RHO[far_size-1]) rmax=RHO[far_size-1];
194  sigma_lines[nGR] = false;
195  line_style[nGR]=LINE_STYLE[n];
196  line_marker[nGR]=LINE_MARKER[n];
197  line_color[nGR]=LINE_COLOR[n];
198  FAR_name[nGR]=FAR_NAME[n];
199  gr[nGR++] = new TGraphErrors(far_size,RHO.data,FAR.data,eRHO.data,eFAR.data);
200  gr[nGR-1]->SetMarkerStyle(20);
201  gr[nGR-1]->SetMarkerSize(0.35);
202  if(n==0) gr[nGR-1]->SetMarkerColor(8);
203  if(n==1) gr[nGR-1]->SetMarkerColor(kBlack);
204  if(n==2) gr[nGR-1]->SetMarkerColor(kRed);
205  if(n==3) gr[nGR-1]->SetMarkerColor(8);
206  if(n==4) gr[nGR-1]->SetMarkerColor(kBlack);
207  if(n==5) gr[nGR-1]->SetMarkerColor(kRed);
208  gr[nGR-1]->SetName(TString::Format("gr%d",nGR-1));
209 
210  }
211 
212  // GW150914
213  RHO.data[0]=14.1*sqrt(2); FAR.data[0]=1./(1369200.00/(365*24*3600));
214  eRHO.data[0]=0; eFAR.data[0]=0;
215  gr[nGR] = new TGraphErrors(1,RHO.data,FAR.data,eRHO.data,eFAR.data);
216  gr[nGR]->SetMarkerColor(4);
217  gr[nGR]->SetMarkerStyle(22);
218  gr[nGR]->SetMarkerSize(2);
219  FAR_name[nGR]="GW150914";
220  nGR++;
221 
222 
223  TMultiGraph* mg = new TMultiGraph();
224  char gTitle[256];
225  if(FAR_TITLE!="") sprintf(gTitle,FAR_TITLE.Data());
226  else sprintf(gTitle,"FAR Comparison");
227 // PAPER
228 //sprintf(gTitle,"ER8b/O1 (C00) background of unmodeled/constrained/chirp data sets (Sep 12 - Oct 20)");
229 #if ( MODE == "inclusive")
230  //sprintf(gTitle,"ER8b/O1 (C01) background of unmodeled/constrained/chirp data sets (Sep 12 - Oct 20)");
231  sprintf(gTitle,"ER8b/O1 (C0101) background of unmodeled/constrained/chirp data sets (Sep 12 - Oct 20)");
232 #endif
233 #if ( MODE == "exclusive")
234  //sprintf(gTitle,"ER8b/O1 (C01) background of eunmodeled/econstrained/chirp data sets (Sep 12 - Oct 20)");
235  //sprintf(gTitle,"ER8b/O1 (C01) background of edge/bulk/chirp data sets (Sep 12 - Oct 20)");
236  sprintf(gTitle,"ER8b/O1 (C0101) background of edge/bulk/chirp data sets (Sep 12 - Oct 20)");
237 #endif
238  sprintf(gTitle,"");
239  mg->SetName("mg");
240  mg->SetTitle(gTitle);
241  for(int n=0;n<nGR;n++) mg->Add(gr[n]);
242  //mg->Paint("APL");
243  mg->Paint("AL");
244 
245  mg->GetHistogram()->GetXaxis()->SetLabelSize(0.05);
246  mg->GetHistogram()->GetYaxis()->SetLabelSize(0.05);
247  mg->GetHistogram()->GetXaxis()->SetTitleSize(0.05);
248  mg->GetHistogram()->GetYaxis()->SetTitleSize(0.05);
249  mg->GetHistogram()->GetXaxis()->SetLabelFont(42);
250  mg->GetHistogram()->GetYaxis()->SetLabelFont(42);
251  mg->GetHistogram()->GetYaxis()->SetLabelOffset(0.01);
252  mg->GetHistogram()->GetYaxis()->SetTitleOffset(1.5);
253 // mg->GetHistogram()->SetMinimum(far_min/2.);
254 
255  if(rho_min>=0) rmin=rho_min*sqrt(2);
256  if(rho_max>=0) rmax=rho_max*sqrt(2);
257  rmax=16*sqrt(2);
258  mg->GetHistogram()->GetXaxis()->SetRangeUser(rmin,rmax);
259 
260  mg->GetXaxis()->SetTitle(gr[0]->GetXaxis()->GetTitle());
261  mg->GetXaxis()->SetLabelFont(42);
262  mg->GetYaxis()->SetLabelFont(42);
263  mg->GetXaxis()->SetTitleFont(42);
264  mg->GetYaxis()->SetTitleFont(42);
265  mg->GetXaxis()->SetTitleOffset(1.20);
266  mg->GetYaxis()->SetTitleOffset(1.20);
267  mg->GetXaxis()->SetTitleSize(0.04);
268  mg->GetYaxis()->SetTitleSize(0.04);
269  mg->GetXaxis()->CenterTitle(true); //PAPER
270  mg->GetYaxis()->CenterTitle(true); //PAPER
271  mg->GetXaxis()->SetTitle("Coherent Network SNR ( #eta_{c} )");
272  //mg->GetXaxis()->SetTitle("#eta_{c}");
273  mg->GetYaxis()->SetTitle("FAR ( yr^{-1} )");
274 
275  mg->Draw("ALP");
276 
277  // draw label GW150914
278  //TLatex *pS1 = new TLatex(RHO[0], FAR[0], "GW150914"); //PAPER
279  //TLatex *pS1 = new TLatex(18.30,2.14, "GW150914"); //PAPER
280  TLatex *pS1 = new TLatex(18.93,2.61, "GW150914"); //PAPER
281  pS1->SetTextFont(52);
282  pS1->SetTextSize(0.040); //PAPER
283  pS1->SetLineWidth(2);
284  pS1->SetTextColor(1);
285 // pS1->Draw();
286 
287  TPad *pad = canvas->GetPad(0);;
288  pad->cd();
289  TFile *froot = new TFile("chirp_set_low_high_rho_paper.root", "RECREATE");
290  pad->Write("pad");
291  froot->Close();
292 
293  // draw the legend
294 
295  TLegend* leg;
296  //leg = new TLegend(0.65,0.24,0.87,0.40,NULL,"brNDC");
297  //leg = new TLegend(0.55,0.14,0.87,0.50,NULL,"brNDC");
298 // leg = new TLegend(0.55,0.14,0.87,0.40,NULL,"brNDC");
299 // leg = new TLegend(0.404,0.495,0.704,0.855,NULL,"brNDC"); // 3 sets
300 // leg = new TLegend(0.5791457,0.6751055,0.879397,0.8586498,NULL,"brNDC");
301 // leg = new TLegend(0.3278894,0.6603376,0.6281407,0.8438819,NULL,"brNDC");
302  leg = new TLegend(0.4120603,0.6265823,0.6306533,0.8607595,NULL,"brNDC");
303 
304 
305  leg->SetBorderSize(1);
306  leg->SetTextAlign(22);
307  leg->SetTextFont(12);
308  leg->SetLineColor(1);
309  leg->SetLineStyle(1);
310  leg->SetLineWidth(1);
311  leg->SetFillColor(0);
312  leg->SetFillStyle(1001);
313  leg->SetTextSize(0.03);
314  leg->SetLineColor(kBlack);
315  leg->SetFillColor(kWhite);
316 
317 cout << nGR << endl;
318  for(int n=0;n<nGR;n++) {
319  if(n==3) continue;
320  if(n==4) continue;
321  if(n==5) continue;
322  char legLabel[256];
323 #if ( MODE == "inclusive")
324  //if(n==0) strcpy(legLabel,"inclusive UnModeled");
325  //if(n==1) strcpy(legLabel,"inclusive Constrained");
326  if(n==0) strcpy(legLabel,"UnModeled");
327  if(n==1) strcpy(legLabel,"Constrained");
328 #endif
329 #if ( MODE == "exclusive")
330  //if(n==0) strcpy(legLabel,"exclusive UnModeled");
331  //if(n==1) strcpy(legLabel,"exclusive Constrained");
332  //if(n==0) strcpy(legLabel,"Edge");
333  //if(n==1) strcpy(legLabel,"Bulk");
334  if(n==0) strcpy(legLabel,"B1");
335  if(n==1) strcpy(legLabel,"B2");
336 #endif
337  if(n==2) strcpy(legLabel,"Chirp");
338  //if(n==2) strcpy(legLabel,"B3");
339  if(n==6) strcpy(legLabel,"GW150914");
340  //sprintf(legLabel,"%s",FAR_name[n].Data());
341  TLegendEntry* eleg = leg->AddEntry(gr[n],legLabel,"lp");
342  //eleg->SetMarkerSize(5);
343  //eleg->SetMarkerColor(kBlack);
344  }
345  leg->SetBorderSize(1);
346  leg->SetTextAlign(22);
347  leg->SetTextFont(12);
348  leg->SetLineColor(1);
349  leg->SetLineStyle(1);
350  leg->SetLineWidth(1);
351  leg->SetFillColor(0);
352  leg->SetFillStyle(1001);
353  leg->SetTextSize(0.04);
354  leg->SetLineColor(kBlack);
355  leg->Draw();
356 
357  // draw markers in the legend (the defaults are too small)
358 
359  TMarker *mP1 = new TMarker(13.2658,2512.81, 20);
360  mP1->SetMarkerStyle(20);
361  mP1->SetMarkerSize(1);
362  mP1->SetMarkerColor(kGreen);
363  mP1->Draw();
364 
365  TMarker *mP2 = new TMarker(13.2658,479.617, 20);
366  mP2->SetMarkerStyle(20);
367  mP2->SetMarkerSize(1);
368  mP2->SetMarkerColor(kBlack);
369  mP2->Draw();
370 
371  TMarker *mP3 = new TMarker(13.2658,91.7268, 20);
372  mP3->SetMarkerStyle(20);
373  mP3->SetMarkerSize(1);
374  mP3->SetMarkerColor(kRed);
375  mP3->Draw();
376 
377  // save plot
378 #if ( MODE == "inclusive")
379  //far_plot_name = "CompareFAR_C0101_UnModeled_Constrained_Chirp_run1ato10a_inclusive";
380  far_plot_name = "CompareFAR_C0101_UnModeled_Constrained_Chirp_run1ato10a_h1_glitch_removed_inclusive";
381 #endif
382 #if ( MODE == "exclusive")
383  //far_plot_name = "CompareFAR_C0101_eUnModeled_eConstrained_Chirp_run1ato10a_exclusive";
384  far_plot_name = "CompareFAR_C0101_eUnModeled_eConstrained_Chirp_run1ato10a_h1_glitch_removed_exclusive";
385 #endif
386  if(far_plot_name!="") {
387  char gfileName[1024];
388  sprintf(gfileName,"%s.gif",far_plot_name.Data());
389  canvas->Print(gfileName);
390  TString pfileName=gfileName;
391  pfileName.ReplaceAll(".gif",".png");
392  char cmd[1024];
393  sprintf(cmd,"convert %s %s",gfileName,pfileName.Data());
394  cout << cmd << endl;
395  gSystem->Exec(cmd);
396  sprintf(cmd,"rm %s",gfileName);
397  cout << cmd << endl;
398  gSystem->Exec(cmd);
399 // gSystem->Exit(0);
400  }
401 
402  return;
403 }
404 
405 int readParameters(TString fname, wavearray<double>& RHO, wavearray<double>& eRHO, wavearray<double>& FAR, wavearray<double>& eFAR) {
406 
407  double rho;
408  double erho;
409  double far;
410  double efar;
411 
412  double year = (24.*3600.*365.);
413 
414  RHO.resize(1000000);
415  eRHO.resize(1000000);
416  FAR.resize(1000000);
417  eFAR.resize(1000000);
418 
419  ifstream in;
420  in.open(fname.Data(),ios::in);
421  if (!in.good()) {cout << "Error Opening File : " << fname.Data() << endl;exit(1);}
422 
423  int n=0;
424  while (1) {
425  //in >> rho >> far >> erho >> efar;
426  in >> rho >> far;
427  if (!in.good()) break;
428 rho*=sqrt(2.);
429 if(rho>22.) continue;
430 far*=year;
431  RHO[n] = rho;
432  eRHO[n] = erho;
433  FAR[n] = far;
434 // eFAR[n] = efar;
435  eFAR[n] = 0;
436 
437  n++;
438  }
439  in.close();
440 
441  RHO.resize(n);
442  eRHO.resize(n);
443  FAR.resize(n);
444  eFAR.resize(n);
445 
446  return n;
447 }
448 
449 int writeParameters(TString fname, wavearray<double>& RHO, wavearray<double>& FAR) {
450 
451  ofstream out;
452  out.open(fname.Data(),ios::out);
453  if (!out.good()) {cout << "Error Opening File : " << fname.Data() << endl;exit(1);}
454 
455  for(int i=0;i<RHO.size();i++) {
456  out << RHO[i] << "\t" << FAR[i] << endl;
457  }
458 
459  out.close();
460 }
461 
462 
rho
double rho
Definition: cwb_report_prod_2.C:149
wavearray::size
virtual size_t size() const
Definition: wavearray.hh:127
FAR
vector< double > FAR
Definition: cwb_mkfad.C:59
cmd
char cmd[1024]
Definition: cwb_condor_cleanup.C:20
eRHO
vector< double > eRHO
Definition: cwb_mkfad.C:56
n
int n
Definition: cwb_net.C:10
rho_min
double rho_min
Definition: cwb_report_prod_2.C:177
TString
TString("c")
Definition: cwb_report_skymap.C:111
leg
TLegend * leg
Definition: compare_bkg.C:246
RHO
vector< double > RHO
Definition: cwb_mkfad.C:55
canvas
return wmap canvas
Definition: DrawClusterWithSkyplot.C:80
LINE_MARKER
int LINE_MARKER[nFAR_MAX]
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:10
colors
Color_t colors[16]
Definition: cwb_report_pe.C:112
eFAR
vector< double > eFAR
Definition: cwb_mkfad.C:60
nFAR_MAX
#define nFAR_MAX
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:2
OBS_TIME
double OBS_TIME[nFAR_MAX]
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:11
i
i drho i
Definition: cwb_epparameters.C:85
ZERO_LAG_TIME
double ZERO_LAG_TIME[nFAR_MAX]
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:13
LINE_STYLE
int LINE_STYLE[nFAR_MAX]
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:9
out
ofstream out
Definition: cwb_merge.C:196
FAR_NAME
TString FAR_NAME[nFAR_MAX]
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:7
gr
TGraph * gr
Definition: cwb_draw_sensitivity.C:4
Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a
void Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a()
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:20
FAR_PATH
TString FAR_PATH[nFAR_MAX]
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:12
rho_max
double rho_max
Definition: cwb_report_prod_2.C:178
FAR_PLOT_NAME
TString FAR_PLOT_NAME
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:5
fname
char fname[1024]
Definition: cwb_report_prod_1.C:110
PlotFAR
void PlotFAR(int nfar, TString far_plot_name, double rho_min=-1, double rho_max=-1)
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:139
FAR_TITLE
TString FAR_TITLE
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:6
writeParameters
int writeParameters(TString fname, wavearray< double > &RHO, wavearray< double > &FAR)
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:449
froot
TFile * froot
Definition: TestReadWriteDetectorObject.C:11
GetXaxis
cout<<"Number of Entries: "<< num<< endl;double *slag1=new double[slag_entries];slag1=wave.GetV1();double *slag2=new double[slag_entries];slag2=wave.GetV2();char mytitle[256];double SlagMax=wave.GetMaximum("slag")+segLen/2.;double SlagMin=wave.GetMinimum("slag")-segLen/2.;int NSlag=TMath::FloorNint((SlagMax-SlagMin)/segLen);cout<< "SLAG MAX : "<< wave.GetMaximum("slag")<< " s SLAG MIN : "<< wave.GetMinimum("slag")<< " s #SLAGS : "<< NSlag-1<< endl;if(NSlag==1){cout<<"Just one slag....Skipping further execution!"<< endl;exit(0);}sprintf(mytitle,"FAR distribution over slags (post cat3 & rho>%f)", T_cut);TH2F *Slag=new TH2F("SLAG", mytitle, NSlag, SlagMin/86400., SlagMax/86400., NSlag, SlagMin/86400., SlagMax/86400.);Slag-> GetXaxis() -> SetTitle("slag[1] shift [day]")
Definition: slag.C:91
nFAR
int nFAR
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:4
in
ifstream in
Definition: CWB_Plugin_HEN_BKG_Config.C:20
LINE_COLOR
int LINE_COLOR[nFAR_MAX]
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:8
strcpy
strcpy(RunLabel, RUN_LABEL)
sprintf
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
wavearray::data
DataType_t * data
Definition: wavearray.hh:301
wavearray::resize
virtual void resize(unsigned int)
Definition: wavearray.cc:445
mg
TMultiGraph * mg
Definition: cwb_draw_sensitivity.C:5
readParameters
int readParameters(TString fname, wavearray< double > &RHO, wavearray< double > &eRHO, wavearray< double > &FAR, wavearray< double > &eFAR)
Definition: Draw_FARvsRHO_ER8b_12Sep20Oct_C0101_BKG_LF_rMRA_run1ato10a.C:405
exit
exit(0)