66 vector<TString> vfpvalue;
67 vector<TString>
vname;
70 for(
int n=0;
n<nHM;
n++) {
73 gNAME[
n].ReplaceAll(
"*",
" ");
75 vector<float> valpha,
vpvalue,vre;
86 if(valpha[
i]>=0.2 && valpha[
i]<3.5) {
87 cout << valpha[
i] <<
"\t" << vpvalue[
i] << endl;
90 wALPHA[
n][
j] = valpha[
i];
91 wPVALUE[
n][
j] = vpvalue[
i];
107 gStyle->SetFrameBorderMode(0);
110 gStyle->SetMarkerColor(50);
111 gStyle->SetLineColor(kWhite);
112 gStyle->SetTitleW(0.98);
113 gStyle->SetTitleH(0.08);
114 gStyle->SetTitleY(0.98);
115 gStyle->SetFillColor(kWhite);
116 gStyle->SetLineColor(kWhite);
117 gStyle->SetTitleFont(12,
"D");
120 TCanvas *
canvas =
new TCanvas(
"Higher Modes",
"Higher Modes", 300, 40, 800, 500);
122 canvas->ToggleEventStatus();
127 canvas->SetFillColor(kWhite);
134 for(
int n=0;
n<nHM;
n++) {
135 #ifdef SHOW_RESIDUALS 136 gr[
n] =
new TGraph(
wSIZE[
n],wALPHA[n].data,wRE[n].data);
138 gr[
n] =
new TGraph(
wSIZE[n],wALPHA[n].data,wPVALUE[n].data);
141 if(wALPHA[n][
i]<xmin) xmin=wALPHA[
n][
i];
142 if(wALPHA[n][
i]<xmax) xmax=wALPHA[
n][
i];
143 if(wPVALUE[n][
i]<ymin) ymin=wPVALUE[
n][
i];
144 if(wPVALUE[n][
i]<ymax) ymax=wPVALUE[
n][
i];
148 for(
int n=0;
n<nHM;
n++) {
149 gr[
n]->SetLineWidth(2);
151 gr[
n]->SetMarkerStyle(20);
152 gr[
n]->SetMarkerSize(0.8);
153 gr[
n]->SetLineColor(
gCOLOR[n]);
154 gr[
n]->SetLineStyle(1);
157 TMultiGraph*
mg =
new TMultiGraph();
159 sprintf(gTitle,
"Higher Modes");
160 mg->SetTitle(gTitle);
161 for(
int n=0;
n<nHM;
n++) mg->Add(gr[
n]);
165 mg->GetHistogram()->GetXaxis()->SetLabelSize(0.04);
166 mg->GetHistogram()->GetYaxis()->SetLabelSize(0.04);
167 mg->GetHistogram()->GetXaxis()->SetTitleSize(0.04);
168 mg->GetHistogram()->GetYaxis()->SetTitleSize(0.04);
169 mg->GetHistogram()->GetXaxis()->SetLabelFont(42);
170 mg->GetHistogram()->GetYaxis()->SetLabelFont(42);
171 mg->GetHistogram()->GetYaxis()->SetLabelOffset(0.01);
172 mg->GetHistogram()->GetYaxis()->SetTitleOffset(1.5);
174 mg->GetHistogram()->GetXaxis()->SetRangeUser(0.001,7.5);
175 #ifndef SHOW_RESIDUALS 176 mg->GetHistogram()->GetYaxis()->SetRangeUser(0.001,1);
179 mg->GetHistogram()->GetXaxis()->SetNdivisions(509);
180 mg->GetHistogram()->GetYaxis()->SetNdivisions(509);
182 mg->GetXaxis()->SetTitle(gr[0]->
GetXaxis()->GetTitle());
183 mg->GetXaxis()->SetLabelFont(42);
184 mg->GetYaxis()->SetLabelFont(42);
185 mg->GetXaxis()->SetTitleFont(42);
186 mg->GetYaxis()->SetTitleFont(42);
187 mg->GetXaxis()->SetTitleOffset(0.80);
188 mg->GetYaxis()->SetTitleOffset(0.90);
189 mg->GetXaxis()->SetTitleSize(0.05);
190 mg->GetYaxis()->SetTitleSize(0.05);
191 mg->GetXaxis()->SetTitle(
"alpha");
192 #ifdef SHOW_RESIDUALS 193 mg->GetYaxis()->SetTitle(
"residuals");
195 mg->GetYaxis()->SetTitle(
"p-value");
197 mg->GetXaxis()->CenterTitle(
true);
198 mg->GetYaxis()->CenterTitle(
true);
205 double hleg = 0.18+nHM*0.06;
206 leg =
new TLegend(0.593985,0.128692,0.8834586,hleg,NULL,
"brNDC");
208 leg->SetBorderSize(1);
209 leg->SetTextAlign(22);
210 leg->SetTextFont(12);
211 leg->SetLineColor(1);
212 leg->SetLineStyle(1);
213 leg->SetLineWidth(1);
214 leg->SetFillColor(0);
215 leg->SetFillStyle(1001);
216 leg->SetTextSize(0.05);
217 leg->SetLineColor(kBlack);
218 leg->SetFillColor(kWhite);
220 for(
int n=0;n<nHM;n++) {
223 leg->AddEntry(gr[n],legLabel,
"lp");
228 if(ofname!=
"") canvas->Print(ofname);
235 float a1,
a2,pvalue,re;
239 if (!in.good()) {cout <<
"Error Opening File : " << fname.Data() << endl;
exit(1);}
246 in >> dy >> dy >> dy >> re >> dy >> pvalue >> dy >> dy >> dy >> a1 >> dy >> a2 >> dy >> dy >> dy >> dy;
247 if (!in.good())
break;
248 valpha.push_back((a1+a2)/2);
249 vpvalue.push_back(pvalue);
254 return vpvalue.size();
264 if (!in.good()) {cout <<
"Error Opening File : " << fname.Data() << endl;
exit(1);}
267 in >> fpvalue >>
name;
268 if (!in.good())
break;
269 if(
TString(fpvalue[0])==
"#")
continue;
270 vfpvalue.push_back(fpvalue);
271 vname.push_back(name);
wavearray< double > wALPHA[nHM_MAX]
wavearray< double > wPVALUE[nHM_MAX]
wavearray< double > wRE[nHM_MAX]
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]")
void DrawPvalueHM(TString ifname, TString ofname="")
double vpvalue[NTRIALS][NPVALUES]
void PlotPvalueHM(int nroc, TString ofname="")
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
int readPvalue(TString fname, vector< float > &valpha, vector< float > &vpvalue, vector< float > &vre)
int readConfig(TString fname, vector< TString > &vfpvalue, vector< TString > &vname)
virtual void resize(unsigned int)