cwb_report_prod_2.C File Reference

#include "GToolbox.hh"
#include <map>

Include dependency graph for cwb_report_prod_2.C:

Functions
wave	Add (net_file_name)
live	Add (liv_file_name)
c1	Clear ()
out_far	close ()
c1	Divide (1, nIFO)
wave	Draw (TString::Format("rho[%d]", pp_irho).Data(),"","goff")
pf_vED	Draw ("colz")
cutF	Draw ("")
gr2	Draw ("AP")
gSystem	Exec ("date")
cout<< "End CWB::Toolbox::getLiveTime : ";cout.flush();gSystem->	Exec ("/bin/date")
	fclose (flive)
	fclose (ftrig)
	fclose (filelags)
	fclose (frate)
	fclose (fratev)
	fclose (fratem)
	for (int n=0;n< nIFO;n++)
gROOT	ForceStyle ()
	fprintf (flive,"nonzero lags live=%10.1f \n", liveTot)
	fprintf (ftrig,"# correlation threshold = %f \n", T_cor)
	fprintf (ftrig,"# -/+ - not passed/passed final selection cuts\n")
	fprintf (ftrig,"# 1 - effective correlated amplitude rho \n")
	fprintf (ftrig,"# 2 - correlation coefficient 0/1 (1G/2G)\n")
	fprintf (ftrig,"# 3 - correlation coefficient 2\n")
	fprintf (ftrig,"# 4 - correlation coefficient 3\n")
	fprintf (ftrig,"# 5 - correlated amplitude \n")
	fprintf (ftrig,"# 6 - time shift \n")
	fprintf (ftrig,"# 7 - time super shift \n")
	fprintf (ftrig,"# 8 - likelihood \n")
	fprintf (ftrig,"# 9 - penalty factor \n")
	fprintf (ftrig,"# 10 - energy disbalance \n")
	fprintf (ftrig,"# 11 - central frequency \n")
	fprintf (ftrig,"# 12 - bandwidth \n")
	fprintf (ftrig,"# 13 - duration \n")
	fprintf (ftrig,"# 14 - number of pixels \n")
	fprintf (ftrig,"# 15 - frequency resolution \n")
	fprintf (ftrig,"# 16 - cwb run number \n")
gr2	GetHistogram () -> SetXTitle("lag")
xCor	GetXaxis () -> SetTitle("network correlation")
xCor	GetYaxis () -> SetTitle("events")
	if (nIFO==1)
	if (!slagFound)
	if (flive==NULL)
	if ((cwb_lag_number!=0)&&(cwb_slag_number!=0))
	if (TString(analysis)=="1G")
	if (pp_rho_vs_rate_no_multiplicity)
	if (!out_far.good())
	if (!out_efar.good())
	if (delete c1)
TChain	live ("liveTime")
TIter	next (wave.GetListOfBranches())
out_far	open (fname, ios::out)
ofstream	out_lf (fname, ios::out)
	printf ("total live time: non-zero lags = %10.1f \n", liveTot)
next	Reset ()
c1	SaveAs (fname)
c1	SetBorderMode (0)
c1	SetBorderSize (2)
c1	SetBottomMargin (0.103939)
gStyle	SetCanvasColor (kWhite)
wave	SetEstimate (nsize)
c1	SetFillColor (0)
cutF	SetFillColor (2)
gStyle	SetFrameBorderMode (0)
c1	SetGridx ()
c1	SetGridy ()
gStyle	SetLabelFont (42,"X")
gStyle	SetLabelFont (42,"Y")
gStyle	SetLabelOffset (0.010,"X")
gStyle	SetLabelOffset (0.010,"Y")
gStyle	SetLabelSize (0.06,"X")
gStyle	SetLabelSize (0.06,"Y")
gStyle	SetLineColor (kWhite)
gr2	SetLineColor (1)
c1	SetLogx (kFALSE)
c1	SetLogy (kFALSE)
c1	SetLogy (kTRUE)
rhocc	SetMarkerColor (1)
rhocc	SetMarkerSize (0.5)
gr2	SetMarkerSize (1)
rhocc	SetMarkerStyle (20)
gStyle	SetNumberContours (256)
gStyle	SetOptStat (kFALSE)
gStyle	SetOptStat (kTRUE)
gStyle	SetPalette (1, 0)
c1	SetRightMargin (0.1517039)
gStyle	SetStatBorderSize (1)
cutF	SetStats (kFALSE)
gStyle	SetTextFont (12)
xCor	SetTitle ("")
Mchirp	SetTitle ("After pp cuts")
rhocc	SetTitle ("0 < cc < 1")
rho_pf	SetTitle ("chi2")
rho_mchirp	SetTitle ("rho vs Mchirp(after pp cuts)")
hrho[0]	SetTitle ("full events (black), after pp cuts (green)")
gr2	SetTitle (_title)
gStyle	SetTitleColor (kBlue,"D")
gStyle	SetTitleFillColor (kWhite)
gStyle	SetTitleFont (42,"X")
gStyle	SetTitleFont (42,"Y")
gStyle	SetTitleFont (12,"D")
gStyle	SetTitleH (0.060)
gStyle	SetTitleOffset (1.4,"X")
gStyle	SetTitleOffset (2.0,"Y")
xCor	SetTitleOffset (1.3,"Y")
gStyle	SetTitleSize (0.05,"X")
gStyle	SetTitleSize (0.04,"Y")
gStyle	SetTitleW (0.95)
gStyle	SetTitleY (0.99)
c1	SetTopMargin (0.0772727)
	sprintf (fname,"%s/live.txt", netdir)
	sprintf (fname,"%s/events.txt", netdir)
	sprintf (fname,"%s/lags.txt", netdir)
	sprintf (fname,"%s/far_rho.txt", netdir)
	sprintf (fname,"%s/efar_rho.txt", netdir)
	sprintf (fname,"%s/rate_threshold.txt", netdir)
	sprintf (fname,"%s/rate_threshold_veto.txt", netdir)
	sprintf (fname,"%s/rate_threshold_multi.txt", netdir)
	sprintf (fname,"%s/sigma.txt", netdir)
	sprintf (fname,"%s/penalty_vED.eps", netdir)
	sprintf (fname,"%s/rho_L.eps", netdir)
	sprintf (fname,"%s/penalty_ed.eps", netdir)
	sprintf (fname,"%s/penalty_cc.eps", netdir)
	sprintf (fname,"%s/fraction_time.eps", netdir)
	sprintf (fname,"%s/fraction_frequency.eps", netdir)
	sprintf (fname,"%s/frequency.eps", netdir)
	sprintf (fname,"%s/rate_time.eps", netdir)
	sprintf (_title,"after pp-cuts & rho > %3.2f", T_out)
	sprintf (fname,"%s/rate_lag.eps", netdir)
	sprintf (fname,"%s/density.eps", netdir)
	sprintf (fname,"%s/penalty.eps", netdir)
	sprintf (fname,"%s/vED.eps", netdir)
	sprintf (fname,"%s/netcor.eps", netdir)
c1	Update ()
TChain	wave ("waveburst")
	while ((branch=(TBranch *) next()))

Variables
char	_title [1024]
double	acor
double	asnr
bool	bcat2 = false
bool	bcat3 = false
bool	bhveto = false
TBranch *	branch
TCanvas *	c1 = new TCanvas("c","C",0,0,600,600)
char	ch1 [256]
char	ch3 [256]
cout<< fname<< endl;out_lf.close();sprintf(fname,"%s/sigma_veto.txt", netdir);ofstream out_veto(fname, ios::out);for(j=0;j < pp_rho_bin;j++){int elag=0;double mean=0;double sigma=0;double sum=0;for(int k=1;k< nlag-1;k++){float rate=0.0;if(Tlag[k]>0) rate=y_lag_veto[k][j]/Tlag[k];mean+=rate;sum+=rate *rate;elag++;}if(elag >0){mean=mean/elag;sigma=sqrt((sum-elag *mean *mean)/elag);}sprintf(output,"%3.2f %4.3e %4.3e\n", Xcut[j], Ycut_veto[j], sigma);out_veto<< output;}cout<< fname<< endl;out_veto.close();CWB::Toolbox::doPoissonPlot(nIFO, Wlag, Tlag, Rlag, netdir);c1->	Clear ()
int	countVETO [100]
TH1F *	cutF = new TH1F("cutF","",24,60,2048.)
TH1F *	dens = new TH1F("dens","",100,0.,3)
TH1F *	ecor = new TH1F("ecor","",200,0,10.)
TH1F *	ECOR = new TH1F("ECOR","",200,0,10.)
TH2F *	ECOR_cc = new TH2F("ECOR_cc","",100,0.,1.,500,0,15.)
TH2F *	ecor_cc = new TH2F("ecor_cc","",100,0.,1.,500,0,15.)
wavearray< double >	Elag = Tlag
	else
double	far_drho = double(far_rho_max-far_rho_min)/double(far_rho_bin)
int	far_rho_bin = float(far_rho_max-far_rho_min)/far_rho_wbin
float	far_rho_max = TMath::Max(pp_rho_max,(double)TMath::Nint(rho_max+0.5))
float	far_rho_min = TMath::Min(pp_rho_min,(double)TMath::Nint(rho_min-0.5))
float	far_rho_wbin = 0.01
FILE *	filelags = fopen(fname,"w")
FILE *	flive = fopen(fname,"w")
char	fname [1024]
FILE *	frate = fopen(fname,"w")
FILE *	fratem = fopen(fname,"w")
FILE *	fratev = fopen(fname,"w")
FILE *	ftrig = fopen(fname,"w")
double	gps
TGraphErrors *	gr2 =new TGraphErrors(Wlag[nIFO].size(),Tdlag.data,Rlag.data,Olag.data,Elag.data)
TH1F *	hF [NIFO_MAX+1]
TH1F *	hF2 [NIFO_MAX+1]
TH1F *	hpf = new TH1F("hpf","",100,0.,1.)
TH1F *	hR [3]
int	hR_col [3] = {2,1,3}
TH1F *	hvED = new TH1F("hvED","",200,0,2.)
int	i = 16
char	IFO [NIFO_MAX][32]
delete[]	iy_lag
int	j
int	k
wavearray< int >	LAG_PASS (ntrg) =0
cout<< "Start CWB::Toolbox::getLiveTime : be patient, it takes a while ..." << endl;liveTot=CWB::Toolbox::getLiveTime(nIFO, live, Trun, Wlag, Wslag, Tlag, Tdlag, cwb_lag_number, cwb_slag_number);std::map < std::pair< int, int >, int >	lagslag
TH1F *	Like = new TH1F("Like","",100,1,4.)
double	Live
double	liveTot = 0.
double	liveZero = 0.
int	m
int	max_lag =0
TH1F *	McErr = new TH1F("McErr","",100,0.,1.)
TH1F *	Mchirp = new TH1F("Mchirp","",100,0.,40.)
int	min_lag =kMaxInt
int	n
TString	nameVETO [100]
int const	nlag =Wlag[nIFO].size()
int	nsize = wave.GetEntries()
size_t	ntype = 0
int	nXDQF =0
wavearray< double >	Olag = Tlag
ofstream	out_efar
ofstream	out_far
char	output [256]
TH2F *	pf_cc = new TH2F("pf_cc","",100,0.,1.,100,0,1.)
TH2F *	pf_ed = new TH2F("pf_ed","",100,0.,1.,100,0.,1.)
TH2F *	pf_vED = new TH2F("pf_vED","",100,0.,1.,100,0,1.)
double	rho
TH2F *	rho_L = new TH2F("rho_L","",500,0.,5.,500,0,5.)
double	rho_max = TMath::MaxElement(nsize, wave.GetV1())
TH2F *	rho_mchirp = new TH2F("rho_mchirp","",200,0.,30.,100,pp_rho_min,pp_rho_max)
double	rho_min = TMath::MinElement(nsize, wave.GetV1())
TH2F *	rho_pf
TH2F *	rho_subnet = new TH2F("rho_subnet","",100,0.,1.,100,pp_rho_min,pp_rho_max)
TH2F *	rho_vED = new TH2F("rho_vED","",100,0.,1.,100,pp_rho_min,pp_rho_max)
TH2F *	rhocc = new TH2F("rhocc","",100,0.,1.,100,pp_rho_min,pp_rho_max)
wavearray< double >	Rlag = Tlag
Int_t	rUn
bool	save
wavearray< int >	SAVE (ntrg) =0
bool	save_veto
bool	saveVETO [100]
bool	slagFound =false
skymap	sm_theta_phi (3)
double	sTARt
double	sTOp
wavearray< double >	Tdlag
wavearray< double >	Tlag
wavearray< double >	Trun (500000) = 0.
double	vED
UChar_t	VETO [100]
netevent W &	wave
wavearray< double >	Wlag [NIFO_MAX+1]
wavearray< int >	Wsel (ntrg)
wavearray< double >	Wslag [NIFO_MAX+1]
double	xCOR
double	xcor
TH1F *	xCor = new TH1F("xCor","",100,0.4,1.)
double	xcor2
double	xcor3
wavearray< double >	Xcut (pp_rho_bin)
dqfile *	XDQF = new dqfile[nVDQF]
wavearray< double >	Xerr (pp_rho_bin) = 0.
wavearray< double >	Xfar (far_rho_bin)
float	xlag
double	xlive
Int_t	xrun
wavearray< double >	Ycut (pp_rho_bin) = 0.
wavearray< double >	yCUT (pp_rho_bin) = 0.
wavearray< double >	Ycut_multi (pp_rho_bin) = 0.
wavearray< double >	Ycut_veto (pp_rho_bin) = 0.
wavearray< double >	Yerr (pp_rho_bin)
wavearray< double >	yERR (pp_rho_bin)
wavearray< double >	Yerr_multi (pp_rho_bin)
wavearray< double >	Yerr_veto (pp_rho_bin)
wavearray< double >	Yfar (far_rho_bin) = 0.

Function Documentation

wave Add

(

net_file_name

)

live Add

(

liv_file_name

)

c1 Clear

(

)

Definition at line 1261 of file CWB_Plugin_xWRC.C.

out_far close

(

)

c1 Divide	(	1	,
		nIFO
	)

wave Draw	(	TString::	Format"rho[%d]", pp_irho).Data(,
		""	,
		"goff"
	)

pf_cc Draw

(

"colz"

)

xCor Draw

(

""

)

FarPlot2 Draw

(

"AP"

)

gSystem Exec

(

"date"

)

cout<< "End CWB::Toolbox::getLiveTime : ";cout.flush();gSystem-> Exec

(

"/bin/date"

)

fclose

(

flive

)

fclose

(

ftrig

)

fclose

(

filelags

)

fclose

(

frate

)

fclose

(

fratev

)

fclose

(

fratem

)

for

(

)

Definition at line 61 of file cwb_report_prod_2.C.

gROOT ForceStyle

(

)

fprintf

(

flive

)

fprintf

(

ftrig

)

fprintf	(	ftrig	,
		"# -/+ - not passed/passed final selection cuts\n"
	)

fprintf	(	ftrig	,
		"# 1 - effective correlated amplitude rho \n"
	)

fprintf	(	ftrig	,
		"# 2 - correlation coefficient 0/1 (1G/2G)\n"
	)

fprintf	(	ftrig	,
		"# 3 - correlation coefficient 2\n"
	)

fprintf	(	ftrig	,
		"# 4 - correlation coefficient 3\n"
	)

fprintf	(	ftrig	,
		"# 5 - correlated amplitude \n"
	)

fprintf	(	ftrig	,
		"# 6 - time shift \n"
	)

fprintf	(	ftrig	,
		"# 7 - time super shift \n"
	)

fprintf	(	ftrig	,
		"# 8 - likelihood \n"
	)

fprintf	(	ftrig	,
		"# 9 - penalty factor \n"
	)

fprintf	(	ftrig	,
		"# 10 - energy disbalance \n"
	)

fprintf	(	ftrig	,
		"# 11 - central frequency \n"
	)

fprintf	(	ftrig	,
		"# 12 - bandwidth \n"
	)

fprintf	(	ftrig	,
		"# 13 - duration \n"
	)

fprintf	(	ftrig	,
		"# 14 - number of pixels \n"
	)

fprintf	(	ftrig	,
		"# 15 - frequency resolution \n"
	)

fprintf	(	ftrig	,
		"# 16 - cwb run number \n"
	)

FarPlot2 GetHistogram

(

)

-> SetXTitle("lag")

Definition at line 913 of file cbc_plots.C.

xCor GetXaxis

(

)

-> SetTitle("network correlation")

Definition at line 321 of file cwb_report_prod_2.C.

xCor GetYaxis

(

)

-> SetTitle("events")

Definition at line 213 of file cbc_plots.C.

if

(

nIFO

= =1

)

Definition at line 52 of file cwb_report_prod_2.C.

if

(

!

slagFound

)

Definition at line 114 of file cwb_report_prod_2.C.

if

(

flive

= =NULL

)

Definition at line 223 of file cwb_report_prod_2.C.

if

(

(cwb_lag_number!=0)&&(cwb_slag_number!=0)

)

Definition at line 229 of file cwb_report_prod_2.C.

if

(

TString(analysis)

= ="1G"

)

Definition at line 389 of file cwb_report_prod_2.C.

if

(

pp_rho_vs_rate_no_multiplicity

)

Definition at line 811 of file cwb_report_prod_2.C.

if

(

!out_far.

good()

)

Definition at line 835 of file cwb_report_prod_2.C.

if

(

!out_efar.

good()

)

Definition at line 839 of file cwb_report_prod_2.C.

if

(

delete

c1

)

Definition at line 947 of file cwb_report_prod_2.C.

TChain live

(

"liveTime"

)

TIter next

(

wave.

GetListOfBranches()

)

out_efar open	(	fname	,
		ios::out
	)

ofstream out_lf	(	fname	,
		ios::out
	)

printf

(

)

next Reset

(

)

c1 SaveAs

(

fname

)

c1 SetBorderMode

(

0

)

c1 SetBorderSize

(

2

)

c1 SetBottomMargin

(

0.

103939

)

gStyle SetCanvasColor

(

kWhite

)

wave SetEstimate

(

nsize

)

c1 SetFillColor

(

0

)

q_inj SetFillColor

(

2

)

gStyle SetFrameBorderMode

(

0

)

c1 SetGridx

(

)

c1 SetGridy

(

)

gStyle SetLabelFont	(	42	,
		"X"
	)

gStyle SetLabelFont	(	42	,
		"Y"
	)

gStyle SetLabelOffset	(	0.	010,
		"X"
	)

gStyle SetLabelOffset	(	0.	010,
		"Y"
	)

gStyle SetLabelSize	(	0.	06,
		"X"
	)

gStyle SetLabelSize	(	0.	06,
		"Y"
	)

gStyle SetLineColor

(

kWhite

)

legl SetLineColor

(

1

)

c1 SetLogx

(

kFALSE

)

c2 SetLogy

(

kFALSE

)

c1 SetLogy

(

kTRUE

)

rhocc SetMarkerColor

(

1

)

rhocc SetMarkerSize

(

0.

5

)

FarPlot2 SetMarkerSize

(

1

)

rhocc SetMarkerStyle

(

20

)

gStyle SetNumberContours

(

256

)

gStyle SetOptStat

(

kFALSE

)

gStyle SetOptStat

(

kTRUE

)

gStyle SetPalette	(	1	,
		0
	)

c1 SetRightMargin

(

0.

1517039

)

gStyle SetStatBorderSize

(

1

)

cutF SetStats

(

kFALSE

)

gStyle SetTextFont

(

12

)

h_radius SetTitle

(

""

)

McErr SetTitle

(

"After pp cuts"

)

cout<< endlSh3-> GetXaxis() -> SetTitle

(

)

rho_pf SetTitle

(

"chi2"

)

rho_mchirp SetTitle

(

"rho vs Mchirp(after pp cuts)"

)

hrho [0] SetTitle	(	"full events	black,
		after pp cuts(green)"
	)

gr2 SetTitle

(

_title

)

gStyle SetTitleColor	(	kBlue	,
		"D"
	)

gStyle SetTitleFillColor

(

kWhite

)

gStyle SetTitleFont	(	42	,
		"X"
	)

gStyle SetTitleFont	(	42	,
		"Y"
	)

gStyle SetTitleFont	(	12	,
		"D"
	)

gStyle SetTitleH

(

0.

060

)

gStyle SetTitleOffset	(	1.	4,
		"X"
	)

gStyle SetTitleOffset	(	2.	0,
		"Y"
	)

xCor SetTitleOffset	(	1.	3,
		"Y"
	)

gStyle SetTitleSize	(	0.	05,
		"X"
	)

gStyle SetTitleSize	(	0.	04,
		"Y"
	)

gStyle SetTitleW

(

0.

95

)

gStyle SetTitleY

(

0.

99

)

c1 SetTopMargin

(

0.

0772727

)

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

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

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

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

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

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

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

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

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

sprintf	(	fname	,
		"%s/penalty_vED.eps"	,
		netdir
	)

sprintf	(	fname	,
		"%s/rho_L.eps"	,
		netdir
	)

sprintf	(	fname	,
		"%s/penalty_ed.eps"	,
		netdir
	)

sprintf	(	fname	,
		"%s/penalty_cc.eps"	,
		netdir
	)

sprintf	(	fname	,
		"%s/fraction_time.eps"	,
		netdir
	)

sprintf	(	fname	,
		"%s/fraction_frequency.eps"	,
		netdir
	)

sprintf	(	fname	,
		"%s/frequency.eps"	,
		netdir
	)

sprintf	(	fname	,
		"%s/rate_time.eps"	,
		netdir
	)

sprintf	(	_title	,
		"after pp-cuts &	rho,
		%3.2f"	,
		T_out
	)

sprintf	(	fname	,
		"%s/rate_lag.eps"	,
		netdir
	)

sprintf	(	fname	,
		"%s/density.eps"	,
		netdir
	)

sprintf	(	fname	,
		"%s/penalty.eps"	,
		netdir
	)

sprintf	(	fname	,
		"%s/vED.eps"	,
		netdir
	)

sprintf	(	fname	,
		"%s/netcor.eps"	,
		netdir
	)

c1 Update

(

)

TChain wave

(

"waveburst"

)

while

(

(branch=(TBranch *) next())

)

Definition at line 110 of file cwb_report_prod_2.C.

Variable Documentation

char _title[1024]

Definition at line 156 of file cwb_report_prod_2.C.

double acor

Definition at line 149 of file cwb_report_prod_2.C.

double asnr

Definition at line 149 of file cwb_report_prod_2.C.

bool bcat2 = false

Definition at line 88 of file cwb_report_prod_2.C.

bool bcat3 = false

Definition at line 90 of file cwb_report_prod_2.C.

bool bhveto = false

Definition at line 89 of file cwb_report_prod_2.C.

TBranch* branch

Definition at line 107 of file cwb_report_prod_2.C.

TCanvas* c1 = new TCanvas("c","C",0,0,600,600)

Definition at line 69 of file cwb_report_prod_2.C.

char ch1[256]

Definition at line 157 of file cwb_report_prod_2.C.

char ch3[256]

Definition at line 158 of file cwb_report_prod_2.C.

c1 Clear

(

)

Definition at line 927 of file cwb_report_prod_2.C.

int countVETO[100]

Definition at line 122 of file cwb_report_prod_2.C.

TH1F* cutF = new TH1F("cutF","",24,60,2048.)

Definition at line 330 of file cwb_report_prod_2.C.

TH1F* dens = new TH1F("dens","",100,0.,3)

Definition at line 282 of file cwb_report_prod_2.C.

TH1F* ecor = new TH1F("ecor","",200,0,10.)

Definition at line 300 of file cwb_report_prod_2.C.

TH1F* ECOR = new TH1F("ECOR","",200,0,10.)

Definition at line 306 of file cwb_report_prod_2.C.

TH2F* ECOR_cc = new TH2F("ECOR_cc","",100,0.,1.,500,0,15.)

Definition at line 350 of file cwb_report_prod_2.C.

TH2F* ecor_cc = new TH2F("ecor_cc","",100,0.,1.,500,0,15.)

Definition at line 359 of file cwb_report_prod_2.C.

Elag = Tlag

Definition at line 205 of file cwb_report_prod_2.C.

else

Initial value:

{                       
    rho_pf = new TH2F("rho_pf","",100,-1.,2.,100,pp_rho_min,pp_rho_max)

Definition at line 393 of file cwb_report_prod_2.C.

double far_drho = double(far_rho_max-far_rho_min)/double(far_rho_bin)

Definition at line 184 of file cwb_report_prod_2.C.

int far_rho_bin = float(far_rho_max-far_rho_min)/far_rho_wbin

Definition at line 183 of file cwb_report_prod_2.C.

float far_rho_max = TMath::Max(pp_rho_max,(double)TMath::Nint(rho_max+0.5))

Definition at line 181 of file cwb_report_prod_2.C.

float far_rho_min = TMath::Min(pp_rho_min,(double)TMath::Nint(rho_min-0.5))

Definition at line 180 of file cwb_report_prod_2.C.

float far_rho_wbin = 0.01

Definition at line 182 of file cwb_report_prod_2.C.

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

Definition at line 800 of file cwb_report_prod_2.C.

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

Definition at line 222 of file cwb_report_prod_2.C.

char fname[1024]

Definition at line 155 of file cwb_report_prod_2.C.

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

Definition at line 855 of file cwb_report_prod_2.C.

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

Definition at line 870 of file cwb_report_prod_2.C.

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

Definition at line 862 of file cwb_report_prod_2.C.

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

Definition at line 546 of file cwb_report_prod_2.C.

double gps

Definition at line 194 of file cwb_report_prod_2.C.

gr2 =new TGraphErrors(Wlag[nIFO].size(),Tdlag.data,Rlag.data,Olag.data,Elag.data)

Definition at line 1023 of file cwb_report_prod_2.C.

TH1F* hF[NIFO_MAX+1]

Definition at line 491 of file cwb_report_prod_2.C.

TH1F* hF2[NIFO_MAX+1]

Definition at line 492 of file cwb_report_prod_2.C.

TH1F* hpf = new TH1F("hpf","",100,0.,1.)

Definition at line 288 of file cwb_report_prod_2.C.

TH1F* hR[3]

Definition at line 527 of file cwb_report_prod_2.C.

int hR_col[3] = {2,1,3}

Definition at line 528 of file cwb_report_prod_2.C.

TH1F* hvED = new TH1F("hvED","",200,0,2.)

Definition at line 294 of file cwb_report_prod_2.C.

i = 16

Definition at line 150 of file cwb_report_prod_2.C.

char IFO[NIFO_MAX][32]

Definition at line 60 of file cwb_report_prod_2.C.

delete [] iy_lag

Definition at line 727 of file cwb_report_prod_2.C.

int j

Definition at line 150 of file cwb_report_prod_2.C.

int k

Definition at line 150 of file cwb_report_prod_2.C.

LAG_PASS =0

Definition at line 586 of file cwb_report_prod_2.C.

cout<< "Start CWB::Toolbox::getLiveTime : be patient, it takes a while ..." << endl; liveTot=CWB::Toolbox::getLiveTime(nIFO,live,Trun,Wlag,Wslag,Tlag,Tdlag,cwb_lag_number,cwb_slag_number); std::map<std::pair<int,int>, int> lagslag

Definition at line 212 of file cwb_report_prod_2.C.

TH1F* Like = new TH1F("Like","",100,1,4.)

Definition at line 312 of file cwb_report_prod_2.C.

double Live

Definition at line 197 of file cwb_report_prod_2.C.

double liveTot = 0.

Definition at line 195 of file cwb_report_prod_2.C.

double liveZero = 0.

Definition at line 196 of file cwb_report_prod_2.C.

int m

Definition at line 150 of file cwb_report_prod_2.C.

int max_lag =0

Definition at line 258 of file cwb_report_prod_2.C.

TH1F* McErr = new TH1F("McErr","",100,0.,1.)

Definition at line 324 of file cwb_report_prod_2.C.

TH1F* Mchirp = new TH1F("Mchirp","",100,0.,40.)

Definition at line 318 of file cwb_report_prod_2.C.

int min_lag =kMaxInt

Definition at line 257 of file cwb_report_prod_2.C.

int n

Definition at line 150 of file cwb_report_prod_2.C.

TString nameVETO[100]

Definition at line 123 of file cwb_report_prod_2.C.

int const nlag =Wlag[nIFO].size()

Definition at line 256 of file cwb_report_prod_2.C.

int nsize = wave.GetEntries()

Definition at line 174 of file cwb_report_prod_2.C.

size_t ntype = 0

Definition at line 154 of file cwb_report_prod_2.C.

int nXDQF =0

Definition at line 126 of file cwb_report_prod_2.C.

Olag = Tlag

Definition at line 206 of file cwb_report_prod_2.C.

ofstream out_efar

Definition at line 836 of file cwb_report_prod_2.C.

ofstream out_far

Definition at line 832 of file cwb_report_prod_2.C.

char output[256]

Definition at line 876 of file cwb_report_prod_2.C.

TH2F* pf_cc = new TH2F("pf_cc","",100,0.,1.,100,0,1.)

Definition at line 405 of file cwb_report_prod_2.C.

TH2F* pf_ed = new TH2F("pf_ed","",100,0.,1.,100,0.,1.)

Definition at line 423 of file cwb_report_prod_2.C.

TH2F* pf_vED = new TH2F("pf_vED","",100,0.,1.,100,0,1.)

Definition at line 414 of file cwb_report_prod_2.C.

double rho

Definition at line 149 of file cwb_report_prod_2.C.

TH2F* rho_L = new TH2F("rho_L","",500,0.,5.,500,0,5.)

Definition at line 432 of file cwb_report_prod_2.C.

double rho_max = TMath::MaxElement(nsize, wave.GetV1())

Definition at line 178 of file cwb_report_prod_2.C.

TH2F* rho_mchirp = new TH2F("rho_mchirp","",200,0.,30.,100,pp_rho_min,pp_rho_max)

Definition at line 441 of file cwb_report_prod_2.C.

double rho_min = TMath::MinElement(nsize, wave.GetV1())

Definition at line 177 of file cwb_report_prod_2.C.

TH2F* rho_pf

Definition at line 388 of file cwb_report_prod_2.C.

TH2F* rho_subnet = new TH2F("rho_subnet","",100,0.,1.,100,pp_rho_min,pp_rho_max)

Definition at line 368 of file cwb_report_prod_2.C.

TH2F* rho_vED = new TH2F("rho_vED","",100,0.,1.,100,pp_rho_min,pp_rho_max)

Definition at line 378 of file cwb_report_prod_2.C.

TH2F* rhocc = new TH2F("rhocc","",100,0.,1.,100,pp_rho_min,pp_rho_max)

Definition at line 340 of file cwb_report_prod_2.C.

Rlag = Tlag

Definition at line 204 of file cwb_report_prod_2.C.

Int_t rUn

Definition at line 193 of file cwb_report_prod_2.C.

bool save

Definition at line 151 of file cwb_report_prod_2.C.

SAVE =0

Definition at line 587 of file cwb_report_prod_2.C.

bool save_veto

Definition at line 152 of file cwb_report_prod_2.C.

bool saveVETO[100]

Definition at line 120 of file cwb_report_prod_2.C.

bool slagFound =false

Definition at line 108 of file cwb_report_prod_2.C.

skymap sm_theta_phi(3)

double sTARt

Definition at line 194 of file cwb_report_prod_2.C.

double sTOp

Definition at line 194 of file cwb_report_prod_2.C.

wavearray<double> Tdlag

Definition at line 202 of file cwb_report_prod_2.C.

wavearray<double> Tlag

Definition at line 203 of file cwb_report_prod_2.C.

Trun = 0.

Definition at line 199 of file cwb_report_prod_2.C.

double vED

Definition at line 149 of file cwb_report_prod_2.C.

UChar_t VETO[100]

Definition at line 121 of file cwb_report_prod_2.C.

netevent W& wave

Definition at line 85 of file cwb_report_prod_2.C.

wavearray<double> Wlag[NIFO_MAX+1]

Definition at line 200 of file cwb_report_prod_2.C.

wavearray<int> Wsel(ntrg)

wavearray<double> Wslag[NIFO_MAX+1]

Definition at line 201 of file cwb_report_prod_2.C.

double xCOR

Definition at line 149 of file cwb_report_prod_2.C.

double xcor

Definition at line 149 of file cwb_report_prod_2.C.

TH1F* xCor = new TH1F("xCor","",100,0.4,1.)

Definition at line 276 of file cwb_report_prod_2.C.

double xcor2

Definition at line 149 of file cwb_report_prod_2.C.

double xcor3

Definition at line 149 of file cwb_report_prod_2.C.

wavearray<double> Xcut(pp_rho_bin)

dqfile* XDQF = new dqfile[nVDQF]

Definition at line 127 of file cwb_report_prod_2.C.

Xerr = 0.

Definition at line 163 of file cwb_report_prod_2.C.

wavearray<double> Xfar(far_rho_bin)

float xlag

Definition at line 192 of file cwb_report_prod_2.C.

double xlive

Definition at line 194 of file cwb_report_prod_2.C.

Int_t xrun

Definition at line 193 of file cwb_report_prod_2.C.

Ycut = 0.

Definition at line 161 of file cwb_report_prod_2.C.

yCUT = 0.

Definition at line 162 of file cwb_report_prod_2.C.

Ycut_multi = 0.

Definition at line 168 of file cwb_report_prod_2.C.

Ycut_veto = 0.

Definition at line 166 of file cwb_report_prod_2.C.

wavearray<double> Yerr(pp_rho_bin)

wavearray<double> yERR(pp_rho_bin)

wavearray<double> Yerr_multi(pp_rho_bin)

wavearray<double> Yerr_veto(pp_rho_bin)

Yfar = 0.

Definition at line 187 of file cwb_report_prod_2.C.