Logo coherent WaveBurst  
Library Reference Guide
Logo
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++)
 pointers to detectors More...
 
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 [kewaveform [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

◆ NINJ_MAX

#define NINJ_MAX   50

◆ NMDC_MAX

#define NMDC_MAX   64

◆ NTYPE_MAX

#define NTYPE_MAX   20

Function Documentation

◆ exit()

exit ( )

◆ for()

for ( )

pointers to detectors

Definition at line 76 of file cwb_mkeff.C.

◆ ForceStyle()

gROOT ForceStyle ( )

◆ if() [1/3]

if ( ninj  = =0)

Definition at line 72 of file cwb_mkeff.C.

◆ if() [2/3]

if ( !in2.  good())

◆ if() [3/3]

if ( )

Definition at line 142 of file cwb_mkeff.C.

◆ SetCanvasColor()

gStyle SetCanvasColor ( kWhite  )

◆ SetFrameBorderMode()

gStyle SetFrameBorderMode ( )

◆ SetLabelFont() [1/2]

gStyle SetLabelFont ( 42  ,
"X"   
)

◆ SetLabelFont() [2/2]

gStyle SetLabelFont ( 42  ,
"Y"   
)

◆ SetLineColor()

gStyle SetLineColor ( kWhite  )

◆ SetNpx()

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  )

◆ SetNumberContours()

gStyle SetNumberContours ( 256  )

◆ SetOptStat()

gStyle SetOptStat ( kFALSE  )

◆ SetParameters()

fFit [k] SetParameters ( par0  ,
par1  [k],
par2  [k],
par3  [k],
pp_factor2distance   
)

◆ SetStatBorderSize()

gStyle SetStatBorderSize ( )

◆ SetTextFont()

gStyle SetTextFont ( 12  )

◆ SetTitleColor()

gStyle SetTitleColor ( kBlue  ,
"D"   
)

◆ SetTitleFillColor()

gStyle SetTitleFillColor ( kWhite  )

◆ SetTitleFont() [1/3]

gStyle SetTitleFont ( 42  ,
"X"   
)

◆ SetTitleFont() [2/3]

gStyle SetTitleFont ( 42  ,
"Y"   
)

◆ SetTitleFont() [3/3]

gStyle SetTitleFont ( 12  ,
"D"   
)

◆ SetTitleH()

gStyle SetTitleH ( 0.  050)

◆ SetTitleOffset() [1/2]

gStyle SetTitleOffset ( 1.  0,
"X"   
)

◆ SetTitleOffset() [2/2]

rho_mchirp SetTitleOffset ( 1.  2,
"Y"   
)

◆ SetTitleW()

gStyle SetTitleW ( 0.  95)

◆ SetTitleY()

gStyle SetTitleY ( 0.  98)

Variable Documentation

◆ ewaveform

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 132 of file cwb_mkeff.C.

◆ hrss10

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

Definition at line 139 of file cwb_mkeff.C.

◆ hrss50_bis

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

Definition at line 140 of file cwb_mkeff.C.

◆ hrss90

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

Definition at line 141 of file cwb_mkeff.C.

◆ imdc_fbandwidth

double imdc_fbandwidth[NMDC_MAX]

Definition at line 66 of file cwb_mkeff.C.

◆ imdc_fcentral

double imdc_fcentral[NMDC_MAX]

Definition at line 65 of file cwb_mkeff.C.

◆ imdc_index

size_t imdc_index[NMDC_MAX]

Definition at line 67 of file cwb_mkeff.C.

◆ imdc_iset

size_t imdc_iset[NMDC_MAX]

Definition at line 68 of file cwb_mkeff.C.

◆ imdc_name

char imdc_name[NMDC_MAX][128]

Definition at line 64 of file cwb_mkeff.C.

◆ imdc_set

char imdc_set[NMDC_MAX][128]

Definition at line 62 of file cwb_mkeff.C.

◆ imdc_set_name

TString* imdc_set_name = new TString[ninj]

Definition at line 74 of file cwb_mkeff.C.

◆ imdc_type

size_t imdc_type[NMDC_MAX]

Definition at line 63 of file cwb_mkeff.C.

◆ ninj

int ninj
Initial value:
#define NMDC_MAX
char mdc_inj_file[1024]
Definition: cwb_dump_inj.C:98
char imdc_name[NMDC_MAX][128]
Definition: cwb_mkeff.C:64
int ReadInjType(TString ifName, int ntype_max, char set[][128], size_t type[], char name[][128], double fcentral[], double fbandwidth[])
Definition: Toolfun.hh:816
size_t imdc_type[NMDC_MAX]
Definition: cwb_mkeff.C:63
double imdc_fcentral[NMDC_MAX]
Definition: cwb_mkeff.C:65
char imdc_set[NMDC_MAX][128]
Definition: cwb_mkeff.C:62
double imdc_fbandwidth[NMDC_MAX]
Definition: cwb_mkeff.C:66

Definition at line 70 of file cwb_mkeff.C.

◆ nset

int nset =0

Definition at line 75 of file cwb_mkeff.C.