cwb_report_sim.C File Reference

Macros
#define	NMDC_MAX   64

Functions
cout<< fname<< endl;FILE *fp=fopen(fname,"w");if(fp==NULL){cout << "cwb_report_sim.C - Error opening file : "<< fname << endl;gSystem->	Exit (1)
	exit (0)
	fclose (fp)
	fclose (in_all)
	fclose (evt_all)
	for (int n=0;n< nIFO;n++)
gROOT	ForceStyle ()
	if (nIFO==1)
	if (pp_eff_vs_thr_quit=="TRUE")
gStyle	SetCanvasColor (kWhite)
gStyle	SetFrameBorderMode (0)
gStyle	SetLabelFont (42,"X")
gStyle	SetLabelFont (42,"Y")
gStyle	SetLabelOffset (0.014,"X")
gStyle	SetLabelOffset (0.010,"Y")
gStyle	SetLabelSize (0.03,"X")
gStyle	SetLabelSize (0.03,"Y")
gStyle	SetNumberContours (256)
gStyle	SetStatBorderSize (1)
gStyle	SetStatColor (0)
gStyle	SetStatH (0.160)
gStyle	SetStatW (0.200)
gStyle	SetStatX (0.878)
gStyle	SetStatY (0.918)
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.4,"X")
gStyle	SetTitleOffset (1.2,"Y")
gStyle	SetTitleW (0.95)
gStyle	SetTitleY (0.98)
	sprintf (fname,"%s/eff_%d_threshold_factors.txt", netdir, pp_eff, imdc_name[i])
	sprintf (fname,"%s/fit_parameters_ALL.txt", netdir)
	sprintf (fname,"%s/evt_parameters_ALL.txt", netdir)
pp_eff_vs_thr_quit	ToUpper ()

Variables
FILE *	evt_all = fopen(fname,"w")
cout<<"cwb_report_sim.C starts..."<< endl;Color_t colors[NMDC_MAX]={6, 3, 2, 8, 43, 7, 8, 4, 4, 2, 43, 1, 3, 1, 6, 7, 6, 3, 2, 8, 43, 7, 8, 4, 4, 2, 43, 1, 3, 1, 6, 7, 6, 3, 2, 8, 43, 7, 8, 4, 4, 2, 43, 1, 3, 1, 6, 7, 6, 3, 2, 8, 43, 7, 8, 4, 4, 2, 43, 1, 3, 1, 6, 7};Style_t markers[NMDC_MAX]={20, 21, 22, 23, 29, 29, 21, 22, 20, 21, 22, 23, 29, 20, 21, 22, 21, 22, 23, 29, 29, 21, 22, 20, 21, 22, 23, 29, 20, 21, 22, 20, 21, 22, 23, 29, 29, 21, 22, 20, 21, 22, 23, 29, 20, 21, 22, 20, 21, 22, 23, 29, 29, 21, 22, 20, 21, 22, 23, 29, 20, 21, 22, 20};TH2F *HH[NIFO_MAX][NMDC_MAX];TH2F *NRE[NIFO_MAX][NMDC_MAX];TH1F *hT[NMDC_MAX];TH1F *hF[NMDC_MAX];TGraphErrors *EFF[NMDC_MAX];for(int i=0;i < NMDC_MAX;i++){for(int j=0;j < nIFO;j++) HH[j][i]=NULL;for(int j=0;j< nIFO;j++) NRE[j][i]=NULL;hT[i]=NULL;hF[i]=NULL;EFF[i]=NULL;}T_ifar *=(365.*24.*3600.);dqfile *XDQF=new dqfile[nVDQF];char imdc_set[NMDC_MAX][128];size_t imdc_type[NMDC_MAX];char imdc_name[NMDC_MAX][128];double imdc_flow[NMDC_MAX];double imdc_fhigh[NMDC_MAX];double imdc_fcentral[NMDC_MAX];double imdc_fbandwidth[NMDC_MAX];size_t imdc_index[NMDC_MAX];size_t imdc_iset[NMDC_MAX];size_t imdc_tset[NMDC_MAX];int ninj=ReadInjType(mdc_inj_file, NMDC_MAX, imdc_set, imdc_type, imdc_name, imdc_fcentral, imdc_fbandwidth);if(ninj==0){cout << endl;cout << "cwb_report_sim.C : Error - no injection types or bad format"<< endl;cout<< "mdc injection file list : " << mdc_inj_file<< endl;cout << "format must be : mdc_set mdc_type mdc_name mdc_fcentral mdc_fbandwidth"<< endl;cout << "process terminated !!!" << endl;exit(1);}for(int i=0;i < NMDC_MAX;i++){imdc_flow[i]=imdc_fcentral[i]-imdc_fbandwidth[i];if(imdc_flow[i] < fLow) imdc_flow[i]=fLow;if(imdc_flow[i]> fHigh) imdc_flow[i]=fLow;imdc_fhigh[i]=imdc_fcentral[i]+imdc_fbandwidth[i];if(imdc_fhigh[i] < fLow) imdc_fhigh[i]=fHigh;if(imdc_fhigh[i]> fHigh) imdc_fhigh[i]=fHigh;}for(int i=0;i< NMDC_MAX;i++) imdc_index[i]=NMDC_MAX;for(int j=0;j< ninj;j++) imdc_index[imdc_type[j]]=j;for(int j=0;j< ninj;j++){cout<< j << " "<< imdc_index[j]<< endl;if(imdc_index[j]> =ninj){cout<< endl;cout << "cwb_report_sim.C : Error - mdc type must be < num max type inj = "<< ninj<< endl;cout << "check mdc_type injection file list : "<< mdc_inj_file << endl;cout<< "format : mdc_set mdc_type mdc_name mdc_fcentral mdc_fbandwidth" << endl;cout<< "process terminated !!!"<< endl;exit(1);}}TString *imdc_set_name=new TString[ninj];int nset=0;for(int i=0;i< ninj;i++){bool bnew=true;for(int j=0;j< nset;j++) if(imdc_set[i]==imdc_set_name[j]) bnew=false;if(bnew) imdc_set_name[nset++]=imdc_set[i];}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 i=0;i < ninj;i++){for(int j=0;j < nset;j++) if(imdc_set[i]==imdc_set_name[j]) imdc_tset[imdc_type[i]]=j;}for(int i=0;i< ninj;i++){cout<< i << " "<< imdc_name[i] << " "<< imdc_set[i] << " "<< imdc_iset[i] << " "<< imdc_tset[i] << " "<< imdc_type[i] << " "<< imdc_index[i]<< endl;}TString pp_eff_vs_thr_mode=CWB::Toolbox::getParameter(pp_eff_vs_thr,"--mode");if((pp_eff_vs_thr_mode!="DISABLED")&&(pp_eff_vs_thr_mode!="disabled")){TString pp_eff_vs_thr_eff=CWB::Toolbox::getParameter(pp_eff_vs_thr,"--eff");int pp_eff=50;if(pp_eff_vs_thr_eff.IsFloat()) pp_eff=pp_eff_vs_thr_eff.Atoi();if(pp_eff < 0) pp_eff=0;if(pp_eff > 100) pp_eff=100;TChain sim("waveburst");TChain mdc("mdc");sim.Add(sim_file_name);mdc.Add(mdc_file_name);char fname[1024];wavearray< double > factorXX(ninj);factorXX=0;if((pp_eff_vs_thr_mode!="AUTO")&&(pp_eff_vs_thr_mode!="auto")){sprintf(fname,"%s", pp_eff_vs_thr_mode.Data());FILE *fp=fopen(fname,"r");if(fp==NULL){cout << "cwb_report_sim.C - Error opening file : "<< fname << endl;gSystem->	Exit (1)
char	fname [4096]
char	IFO [NIFO_MAX][32]
FILE *	in_all = fopen(fname,"w")
TString	pp_eff_vs_thr_quit = CWB::Toolbox::getParameter(pp_eff_vs_thr,"--quit")
float	xfactors
char	xmdc_name [256]
int	xmdc_type

Macro Definition Documentation

#define NMDC_MAX   64

Function Documentation

cout<< fname << endl; FILE* fp = fopen(fname,"w"); if(fp==NULL) { cout << "cwb_report_sim.C - Error opening file : " << fname << endl; gSystem-> Exit

(

1

)

exit

(

0

)

fclose

(

fp

)

fclose

(

in_all

)

fclose

(

evt_all

)

for

(

)

Definition at line 29 of file cwb_report_sim.C.

gROOT ForceStyle

(

)

if

(

nIFO

= =1

)

Definition at line 20 of file cwb_report_sim.C.

if

(

pp_eff_vs_thr_quit

= ="TRUE"

)

Definition at line 313 of file cwb_report_sim.C.

gStyle SetCanvasColor

(

kWhite

)

gStyle SetFrameBorderMode

(

0

)

gStyle SetLabelFont	(	42	,
		"X"
	)

gStyle SetLabelFont	(	42	,
		"Y"
	)

gStyle SetLabelOffset	(	0.	014,
		"X"
	)

gStyle SetLabelOffset	(	0.	010,
		"Y"
	)

gStyle SetLabelSize	(	0.	03,
		"X"
	)

gStyle SetLabelSize	(	0.	03,
		"Y"
	)

gStyle SetNumberContours

(

256

)

gStyle SetStatBorderSize

(

1

)

gStyle SetStatColor

(

0

)

gStyle SetStatH

(

0.

160

)

gStyle SetStatW

(

0.

200

)

gStyle SetStatX

(

0.

878

)

gStyle SetStatY

(

0.

918

)

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.	4,
		"X"
	)

gStyle SetTitleOffset	(	1.	2,
		"Y"
	)

gStyle SetTitleW

(

0.

95

)

gStyle SetTitleY

(

0.

98

)

sprintf	(	fname	,
		"%s/eff_%d_threshold_factors.txt"	,
		netdir	,
		pp_eff	,
		imdc_name	[i]
	)

sprintf	(	fname	,
		"%s/fit_parameters_ALL.txt"	,
		netdir
	)

sprintf	(	fname	,
		"%s/evt_parameters_ALL.txt"	,
		netdir
	)

pp_eff_vs_thr_quit ToUpper

(

)

Variable Documentation

FILE* evt_all = fopen(fname,"w")

Definition at line 323 of file cwb_report_sim.C.

cout<<"cwb_report_sim.C starts..."<<endl; Color_t colors[NMDC_MAX] = { 6, 3, 2, 8,43, 7, 8, 4, 4, 2,43, 1, 3, 1, 6, 7, 6, 3, 2, 8,43, 7, 8, 4, 4, 2,43, 1, 3, 1, 6, 7, 6, 3, 2, 8,43, 7, 8, 4, 4, 2,43, 1, 3, 1, 6, 7, 6, 3, 2, 8,43, 7, 8, 4, 4, 2,43, 1, 3, 1, 6, 7 }; Style_t markers[NMDC_MAX]= {20,21,22,23,29,29,21,22,20,21,22,23,29,20,21,22, 21,22,23,29,29,21,22,20,21,22,23,29,20,21,22,20, 21,22,23,29,29,21,22,20,21,22,23,29,20,21,22,20, 21,22,23,29,29,21,22,20,21,22,23,29,20,21,22,20 }; TH2F* HH[NIFO_MAX][NMDC_MAX]; TH2F* NRE[NIFO_MAX][NMDC_MAX]; TH1F* hT[NMDC_MAX]; TH1F* hF[NMDC_MAX]; TGraphErrors* EFF[NMDC_MAX]; for(int i=0;i<NMDC_MAX;i++) { for(int j=0;j<nIFO;j++) HH[j][i]=NULL; for(int j=0;j<nIFO;j++) NRE[j][i]=NULL; hT[i]=NULL; hF[i]=NULL; EFF[i]=NULL; } T_ifar*=(365.*24.*3600.); dqfile* XDQF = new dqfile[nVDQF]; char imdc_set[NMDC_MAX][128]; size_t imdc_type[NMDC_MAX]; char imdc_name[NMDC_MAX][128]; double imdc_flow[NMDC_MAX]; double imdc_fhigh[NMDC_MAX]; double imdc_fcentral[NMDC_MAX]; double imdc_fbandwidth[NMDC_MAX]; size_t imdc_index[NMDC_MAX]; size_t imdc_iset[NMDC_MAX]; size_t imdc_tset[NMDC_MAX]; int ninj=ReadInjType(mdc_inj_file,NMDC_MAX,imdc_set,imdc_type, imdc_name,imdc_fcentral,imdc_fbandwidth); if(ninj==0) { cout << endl; cout << "cwb_report_sim.C : Error - no injection types or bad format" << endl; cout << "mdc injection file list : " << mdc_inj_file << endl; cout << "format must be : mdc_set mdc_type mdc_name mdc_fcentral mdc_fbandwidth" << endl; cout << "process terminated !!!" << endl; exit(1); } for(int i=0; i<NMDC_MAX; i++) { imdc_flow[i] = imdc_fcentral[i]-imdc_fbandwidth[i]; if(imdc_flow[i]<fLow) imdc_flow[i]=fLow; if(imdc_flow[i]>fHigh) imdc_flow[i]=fLow; imdc_fhigh[i] = imdc_fcentral[i]+imdc_fbandwidth[i]; if(imdc_fhigh[i]<fLow) imdc_fhigh[i]=fHigh; if(imdc_fhigh[i]>fHigh) imdc_fhigh[i]=fHigh; } for(int i=0; i<NMDC_MAX; i++) imdc_index[i] = NMDC_MAX; for(int j=0;j<ninj;j++) imdc_index[imdc_type[j]]=j; for(int j=0;j<ninj;j++) { cout << j << " " << imdc_index[j] << endl; if(imdc_index[j]>=ninj) { cout << endl; cout << "cwb_report_sim.C : Error - mdc type must be < num max type inj = " << ninj << endl; cout << "check mdc_type injection file list : " << mdc_inj_file << endl; cout << "format : mdc_set mdc_type mdc_name mdc_fcentral mdc_fbandwidth" << endl; cout << "process terminated !!!" << endl; exit(1); } } TString* imdc_set_name = new TString[ninj]; int nset=0; for(int i=0;i<ninj;i++) { bool bnew=true; for(int j=0;j<nset;j++) if(imdc_set[i]==imdc_set_name[j]) bnew=false; if(bnew) imdc_set_name[nset++]=imdc_set[i]; } 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 i=0;i<ninj;i++) { for(int j=0;j<nset;j++) if(imdc_set[i]==imdc_set_name[j]) imdc_tset[imdc_type[i]]=j; } for(int i=0;i<ninj;i++) { cout << i << " " << imdc_name[i] << " " << imdc_set[i] << " " << imdc_iset[i] << " " << imdc_tset[i] << " " << imdc_type[i] << " " << imdc_index[i] << endl; } TString pp_eff_vs_thr_mode = CWB::Toolbox::getParameter(pp_eff_vs_thr,"--mode"); if((pp_eff_vs_thr_mode!="DISABLED")&&(pp_eff_vs_thr_mode!="disabled")) { TString pp_eff_vs_thr_eff = CWB::Toolbox::getParameter(pp_eff_vs_thr,"--eff"); int pp_eff = 50; if(pp_eff_vs_thr_eff.IsFloat()) pp_eff = pp_eff_vs_thr_eff.Atoi(); if(pp_eff<0) pp_eff=0; if(pp_eff>100) pp_eff=100; TChain sim("waveburst"); TChain mdc("mdc"); sim.Add(sim_file_name); mdc.Add(mdc_file_name); char fname[1024]; wavearray<double> factorXX(ninj); factorXX=0; if((pp_eff_vs_thr_mode!="AUTO")&&(pp_eff_vs_thr_mode!="auto")) { sprintf(fname,"%s",pp_eff_vs_thr_mode.Data()); FILE* fp = fopen(fname,"r"); if(fp==NULL) { cout << "cwb_report_sim.C - Error opening file : " << fname << endl; gSystem-> Exit(1)

Definition at line 195 of file cwb_report_sim.C.

char fname[4096]

Definition at line 319 of file cwb_report_sim.C.

char IFO[NIFO_MAX][32]

Definition at line 28 of file cwb_report_sim.C.

FILE* in_all = fopen(fname,"w")

Definition at line 321 of file cwb_report_sim.C.

TString pp_eff_vs_thr_quit = CWB::Toolbox::getParameter(pp_eff_vs_thr,"--quit")

Definition at line 311 of file cwb_report_sim.C.

float xfactors

Definition at line 199 of file cwb_report_sim.C.

char xmdc_name[256]

Definition at line 197 of file cwb_report_sim.C.

int xmdc_type

Definition at line 198 of file cwb_report_sim.C.