Logo Coherent WaveBurst  
Reference Guide
Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Macros | Functions | Variables
cwb_mkeff.C File Reference

Macros

#define NINJ_MAX   50
 
#define NMDC_MAX   64
 
#define NTYPE_MAX   20
 

Functions

 exit (0)
 
 for (int i=0;i< ninj;i++)
 
gROOT ForceStyle ()
 
 if (ninj==0)
 
 if (!in2.good()) break
 
 if ()
 
gStyle SetCanvasColor (kWhite)
 
gStyle SetFrameBorderMode (0)
 
gStyle SetLabelFont (42,"X")
 
gStyle SetLabelFont (42,"Y")
 
gStyle SetLineColor (kWhite)
 
cout<< ewaveform[k].Data()
<< endl;double par0=TMath::Log10(hrss50[iset][k]);fFit[k]=new
TF1("logNfit", logNfit, pow(10.0,
inf), pow(10.0, sup), 5);fFit[k]-> 
SetNpx (100000)
 
gStyle SetNumberContours (256)
 
gStyle SetOptStat (kFALSE)
 
fFit[kSetParameters (par0, par1[k], par2[k], par3[k], pp_factor2distance)
 
gStyle SetStatBorderSize (1)
 
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.0,"X")
 
gStyle SetTitleOffset (1.2,"Y")
 
gStyle SetTitleW (0.95)
 
gStyle SetTitleY (0.98)
 

Variables

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]
 
 hrss10 [iset][k] =fFit[k]->GetX(.1,pow(10.0,inf),pow(10.0,sup))
 
 hrss50_bis [k] =fFit[k]->GetX(.5,pow(10.0,inf),pow(10.0,sup))
 
 hrss90 [iset][k] =fFit[k]->GetX(.9,pow(10.0,inf),pow(10.0,sup))
 
double imdc_fbandwidth [NMDC_MAX]
 
double imdc_fcentral [NMDC_MAX]
 
size_t imdc_index [NMDC_MAX]
 
size_t imdc_iset [NMDC_MAX]
 
char imdc_name [NMDC_MAX][128]
 
char imdc_set [NMDC_MAX][128]
 
TStringimdc_set_name = new TString[ninj]
 
size_t imdc_type [NMDC_MAX]
 
int ninj
 
int nset =0
 

Macro Definition Documentation

#define NINJ_MAX   50
#define NMDC_MAX   64
#define NTYPE_MAX   20

Function Documentation

exit ( )
for ( )

Definition at line 58 of file cwb_mkeff.C.

gROOT ForceStyle ( )
if ( ninj  = =0)

Definition at line 54 of file cwb_mkeff.C.

if ( !in2.  good())
if ( )

Definition at line 124 of file cwb_mkeff.C.

gStyle SetCanvasColor ( kWhite  )
gStyle SetFrameBorderMode ( )
gStyle SetLabelFont ( 42  ,
"X"   
)
gStyle SetLabelFont ( 42  ,
"Y"   
)
gStyle SetLineColor ( kWhite  )
cout<< ewaveform[k].Data() << endl; double par0=TMath::Log10(hrss50[iset][k]); fFit[k] = new TF1("logNfit",logNfit,pow(10.0,inf),pow(10.0,sup),5); fFit[k]-> SetNpx ( 100000  )
gStyle SetNumberContours ( 256  )
gStyle SetOptStat ( kFALSE  )
fFit [k] SetParameters ( par0  ,
par1  [k],
par2  [k],
par3  [k],
pp_factor2distance   
)
gStyle SetStatBorderSize ( )
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.  0,
"X"   
)
rho_pf SetTitleOffset ( 1.  2,
"Y"   
)
gStyle SetTitleW ( 0.  95)
gStyle SetTitleY ( 0.  98)

Variable Documentation

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]

Definition at line 114 of file cwb_mkeff.C.

hrss10[iset][k] =fFit[k]->GetX(.1,pow(10.0,inf),pow(10.0,sup))

Definition at line 121 of file cwb_mkeff.C.

hrss50_bis[k] =fFit[k]->GetX(.5,pow(10.0,inf),pow(10.0,sup))

Definition at line 122 of file cwb_mkeff.C.

hrss90[iset][k] =fFit[k]->GetX(.9,pow(10.0,inf),pow(10.0,sup))

Definition at line 123 of file cwb_mkeff.C.

double imdc_fbandwidth[NMDC_MAX]

Definition at line 48 of file cwb_mkeff.C.

double imdc_fcentral[NMDC_MAX]

Definition at line 47 of file cwb_mkeff.C.

size_t imdc_index[NMDC_MAX]

Definition at line 49 of file cwb_mkeff.C.

size_t imdc_iset[NMDC_MAX]

Definition at line 50 of file cwb_mkeff.C.

char imdc_name[NMDC_MAX][128]

Definition at line 46 of file cwb_mkeff.C.

char imdc_set[NMDC_MAX][128]

Definition at line 44 of file cwb_mkeff.C.

TString* imdc_set_name = new TString[ninj]

Definition at line 56 of file cwb_mkeff.C.

size_t imdc_type[NMDC_MAX]

Definition at line 45 of file cwb_mkeff.C.

int ninj
Initial value:
char mdc_inj_file[1024]
#define NMDC_MAX
char imdc_name[NMDC_MAX][128]
Definition: cwb_mkeff.C:46
int ReadInjType(TString ifName, int ntype_max, char set[][128], size_t type[], char name[][128], double fcentral[], double fbandwidth[])
Definition: Toolfun.hh:740
size_t imdc_type[NMDC_MAX]
Definition: cwb_mkeff.C:45
double imdc_fcentral[NMDC_MAX]
Definition: cwb_mkeff.C:47
char imdc_set[NMDC_MAX][128]
Definition: cwb_mkeff.C:44
double imdc_fbandwidth[NMDC_MAX]
Definition: cwb_mkeff.C:48

Definition at line 52 of file cwb_mkeff.C.

int nset =0

Definition at line 57 of file cwb_mkeff.C.