33 gStyle->SetTitleOffset(1.0,
"X");
34 gStyle->SetTitleOffset(1.2,
"Y");
37 gStyle->SetLabelFont(42,
"X");
38 gStyle->SetLabelFont(42,
"Y");
39 gStyle->SetTitleFont(42,
"X");
40 gStyle->SetTitleFont(42,
"Y");
44 gStyle->SetTitleH(0.050);
45 gStyle->SetTitleW(0.95);
46 gStyle->SetTitleY(0.98);
47 gStyle->SetTitleFont(12,
"D");
48 gStyle->SetTitleColor(kBlue,
"D");
49 gStyle->SetTextFont(12);
50 gStyle->SetTitleFillColor(kWhite);
51 gStyle->SetLineColor(kWhite);
52 gStyle->SetNumberContours(256);
53 gStyle->SetCanvasColor(kWhite);
54 gStyle->SetStatBorderSize(1);
55 gStyle->SetOptStat(kFALSE);
58 gStyle->SetFrameBorderMode(0);
71 imdc_name,imdc_fcentral,imdc_fbandwidth);
72 if(ninj==0) {cout <<
"Error - no injection - terminated" << endl;
exit(1);}
78 for(
int j=0;
j<
nset;
j++)
if(imdc_set[
i]==imdc_set_name[
j]) bnew=
false;
79 if(bnew) imdc_set_name[nset++]=imdc_set[
i];
81 cout <<
"nset : " << nset << endl;
84 for(
int j=0;
j<
ninj;
j++)
if(imdc_set[
j]==imdc_set_name[
i]) imdc_iset[
j]=
i;
86 for (
int iset=0;iset<
nset;iset++) cout << iset <<
" " << imdc_set_name[iset].Data() << endl;
106 sup = log10(pp_factor2distance/
factors[0]);
110 for (
int iset=0;iset<
nset;iset++) {
113 sprintf(file,
"%s/fit_parameters_%s.txt",
netdir, imdc_set_name[iset].Data());
114 cout << file << endl;
118 if (!in2.good()) {cout <<
"Error Opening File : " << file << endl;
exit(1);}
131 in2 >> ecount[
k] >> chi2[
k] >> hrss50[iset][
k] >> piumeno[
k]
132 >> err[
k] >> par1[
k] >> par2[
k] >> par3[
k] >> ewaveform[
k];
133 if (!in2.good())
break;
134 cout << ewaveform[
k].Data() << endl;
135 double par0=TMath::Log10(hrss50[iset][k]);
136 fFit[
k] =
new TF1(
"logNfit",
logNfit,pow(10.0,inf),pow(10.0,sup),5);
137 fFit[
k]->SetNpx(100000);
138 fFit[
k]->SetParameters(par0,par1[k],par2[k],par3[k],pp_factor2distance);
139 hrss10[iset][
k]=fFit[
k]->GetX(.1,pow(10.0,inf),pow(10.0,sup));
140 hrss50_bis[
k] =fFit[
k]->GetX(.5,pow(10.0,inf),pow(10.0,sup));
141 hrss90[iset][
k]=fFit[
k]->GetX(.9,pow(10.0,inf),pow(10.0,sup));
142 if(fFit[k]->Eval(hrss90[iset][k])<0.89) hrss90[iset][
k]=1
e-10;
144 par0=TMath::Log10(hrss50[iset][k]+err[k]);
145 fFit[
k]->SetParameters(par0,par1[k],par2[k],par3[k],pp_factor2distance);
146 ehrss10[
k]=
fabs(fFit[k]->GetX(.1,pow(10.0,inf),pow(10.0,sup))-hrss10[iset][k]);
147 ehrss50[
k]=
fabs(fFit[k]->GetX(.5,pow(10.0,inf),pow(10.0,sup))-hrss50[iset][k]);
148 ehrss90[
k]=
fabs(fFit[k]->GetX(.9,pow(10.0,inf),pow(10.0,sup))-hrss90[iset][k]);
150 cout << hrss10[iset][
k] <<
" " << hrss50[iset][
k] <<
" " << hrss90[iset][
k] << endl;
162 for(
int iset=0;iset<
nset;iset++) {
164 sprintf(etitle,
"%s",imdc_set_name[iset].Data());
166 mg[iset] =
new TMultiGraph();
167 legend[iset] =
new TLegend(0.732,0.125,0.956,0.354,NULL,
"brNDC");
170 gr10[iset] =
new TGraphErrors();
171 gr10[iset]->SetLineColor(1);
172 gr10[iset]->SetLineWidth(1);
173 gr10[iset]->SetMarkerColor(1);
174 gr10[iset]->SetMarkerStyle(20);
175 gr10[iset]->SetLineStyle(7);
177 gr50[iset] =
new TGraphErrors();
178 gr50[iset]->SetLineColor(2);
179 gr50[iset]->SetLineWidth(1);
180 gr50[iset]->SetMarkerColor(2);
181 gr50[iset]->SetMarkerStyle(20);
182 gr50[iset]->SetLineStyle(7);
184 gr90[iset] =
new TGraphErrors();
185 gr90[iset]->SetLineColor(kBlue);
186 gr90[iset]->SetLineWidth(1);
187 gr90[iset]->SetMarkerColor(kBlue);
188 gr90[iset]->SetMarkerStyle(20);
189 gr90[iset]->SetLineStyle(7);
193 if(imdc_iset[
i]!=iset)
continue;
194 if(pp_show_eff_fit_curve) {
195 gr10[iset]->SetPoint(npoint,imdc_fcentral[
i],hrss10[iset][i]);
196 gr50[iset]->SetPoint(npoint,imdc_fcentral[i],hrss50[iset][i]);
197 gr90[iset]->SetPoint(npoint,imdc_fcentral[i],hrss90[iset][i]);
199 gr10[iset]->SetPoint(npoint,imdc_fcentral[i],0.);
200 gr50[iset]->SetPoint(npoint,imdc_fcentral[i],0.);
201 gr90[iset]->SetPoint(npoint,imdc_fcentral[i],0.);
206 mg[iset]->Add(gr10[iset]);
207 mg[iset]->Add(gr50[iset]);
208 mg[iset]->Add(gr90[iset]);
211 sprintf(namecanvas,
"c%i",iset);
212 canvas[iset] =
new TCanvas(namecanvas,namecanvas,125+iset*10,82,976,576);
213 canvas[iset]->Clear();
214 canvas[iset]->ToggleEventStatus();
215 canvas[iset]->SetLogy();
216 #ifdef CWB_MKEFF_LINX 217 canvas[iset]->SetLogx(
false);
219 canvas[iset]->SetLogx(
true);
221 canvas[iset]->SetGridx();
222 canvas[iset]->SetGridy();
223 canvas[iset]->SetFillColor(kWhite);
226 mg[iset]->SetTitle(etitle);
227 mg[iset]->Paint(
"alp");
228 mg[iset]->GetHistogram()->GetXaxis()->SetTitle(
"Frequency (Hz)");
229 mg[iset]->GetHistogram()->GetXaxis()->CenterTitle(
true);
230 mg[iset]->GetHistogram()->GetXaxis()->SetLabelFont(42);
231 mg[iset]->GetHistogram()->GetXaxis()->SetTitleFont(42);
232 mg[iset]->GetHistogram()->GetYaxis()->SetLabelFont(42);
233 mg[iset]->GetHistogram()->GetYaxis()->SetTitleFont(42);
235 mg[iset]->GetHistogram()->GetYaxis()->SetTitle(
"distance (Kpc)");
236 mg[iset]->GetHistogram()->GetYaxis()->SetRangeUser(pp_factor2distance/
factors[
nfactor-1],pp_factor2distance/
factors[0]);
238 mg[iset]->GetHistogram()->GetYaxis()->SetTitle(
"snr");
241 mg[iset]->GetHistogram()->GetYaxis()->SetTitle(
"hrss (1/Hz^{-1/2})");
245 mg[iset]->Draw(
"alp");
247 legend[iset]->SetBorderSize(1);
248 legend[iset]->SetTextAlign(22);
249 legend[iset]->SetTextFont(12);
250 legend[iset]->SetLineColor(1);
251 legend[iset]->SetLineStyle(1);
252 legend[iset]->SetLineWidth(1);
253 legend[iset]->SetFillColor(0);
254 legend[iset]->SetFillStyle(1001);
255 legend[iset]->SetTextSize(0.04);
256 legend[iset]->SetLineColor(kBlack);
257 legend[iset]->SetFillColor(kWhite);
260 legend[iset]->AddEntry(gr90[iset],
"dstance (Kpc) 90%",
"lp");
261 legend[iset]->AddEntry(gr50[iset],
"dstance (Kpc) 50%",
"lp");
262 legend[iset]->AddEntry(gr10[iset],
"dstance (Kpc) 10%",
"lp");
264 legend[iset]->AddEntry(gr90[iset],
"snr 90%",
"lp");
265 legend[iset]->AddEntry(gr50[iset],
"snr 50%",
"lp");
266 legend[iset]->AddEntry(gr10[iset],
"snr 10%",
"lp");
268 legend[iset]->AddEntry(gr90[iset],
"Hrss 90%",
"lp");
269 legend[iset]->AddEntry(gr50[iset],
"Hrss 50%",
"lp");
270 legend[iset]->AddEntry(gr10[iset],
"Hrss 10%",
"lp");
272 legend[iset]->Draw();
274 sprintf(ofile,
"%s/eff_freq_%s.gif",
netdir,imdc_set_name[iset].Data());
275 cout << ofile << endl;
276 canvas[iset]->SaveAs(ofile);
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 iset=0;iset< nset;iset++) cout<< iset<< " "<< imdc_set_name[iset].Data()<< endl;char etitle[256];char ofile[256];TCanvas *canvas[NTYPE_MAX];int ecount[NINJ_MAX];TString piumeno[NINJ_MAX];float chi2[NINJ_MAX], err[NINJ_MAX], par1[NINJ_MAX], par2[NINJ_MAX], par3[NINJ_MAX];double ehrss10[NINJ_MAX], ehrss50[NINJ_MAX], ehrss90[NINJ_MAX];double hrss50_bis[NINJ_MAX];TString ewaveform[NINJ_MAX];TF1 *fFit[NINJ_MAX];double hrss50[NTYPE_MAX][NINJ_MAX], hrss90[NTYPE_MAX][NINJ_MAX], hrss10[NTYPE_MAX][NINJ_MAX];double inf=simulation==2 ? log10(factors[0]) :-25;double sup=simulation==2 ? log10(factors[nfactor-1]) :-18.5;if(simulation==1 &&pp_factor2distance) { inf=log10(pp_factor2distance/factors[nfactor-1]);sup=log10(pp_factor2distance/factors[0]);} int k=0;for(int iset=0;iset< nset;iset++) { char file[256];sprintf(file,"%s/fit_parameters_%s.txt", netdir, imdc_set_name[iset].Data());cout<< file<< endl;ifstream in2;in2.open(file, ios::in);if(!in2.good()) {cout<< "Error Opening File : "<< file<< endl;exit(1);} for(int j=0;j< NINJ_MAX;j++) { hrss50_bis[j]=0;hrss10[iset][j]=0;hrss50[iset][j]=0;hrss90[iset][j]=0;ecount[j]=0;ewaveform[j]="";} for(int l=0;l< NINJ_MAX;l++) { in2 > ecount [k] chi2 [k] hrss50 [iset][k] piumeno [k] err [k] par1 [k] par2 [k] par3 [k] ewaveform[k]
size_t imdc_iset[NMDC_MAX]
size_t imdc_index[NMDC_MAX]
char imdc_name[NMDC_MAX][128]
double fabs(const Complex &x)
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
size_t imdc_type[NMDC_MAX]
double imdc_fcentral[NMDC_MAX]
char imdc_set[NMDC_MAX][128]
double imdc_fbandwidth[NMDC_MAX]