cWB: /home/cWB_docker/git/cWB/library/tools/cwb/macros/cwb_report_prod

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

Include dependency graph for cwb_report_prod_2.C:

Macros
#define	YEAR (24.3600.365.)

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 ("HIST")

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++)
	pointers to detectors More...

gROOT	ForceStyle ()

	fprintf (flive,"nonzero lags live=%10.1f \, liveTot)

	fprintf (ftrig,"# correlation threshold = %f \, T_cor)

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

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

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

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

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

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

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

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

	fprintf (ftrig,"# 8 - likelihood \)

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

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

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

	fprintf (ftrig,"# 12 - bandwidth \)

	fprintf (ftrig,"# 13 - duration \)

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

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

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

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 \, 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\", 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.

Macro Definition Documentation

◆ YEAR

#define YEAR (24.*3600.*365.)

Function Documentation

◆ Add() [1/2]

wave Add ( net_file_name )

◆ Add() [2/2]

live Add ( liv_file_name )

◆ Clear()

c1 Clear ( )

Definition at line 1279 of file CWB_Plugin_xWRC.C.

◆ close()

out_far close ( )

◆ Divide()

c1 Divide	(	1	,
		nIFO
	)

◆ Draw() [1/4]

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

◆ Draw() [2/4]

pf_vED Draw ( "colz" )

◆ Draw() [3/4]

xCor Draw ( "HIST" )

◆ Draw() [4/4]

gr2 Draw ( "AP" )

◆ Exec() [1/2]

gSystem Exec ( "date" )

◆ Exec() [2/2]

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

◆ fclose() [1/6]

fclose ( flive )

◆ fclose() [2/6]

fclose ( ftrig )

◆ fclose() [3/6]

fclose ( filelags )

◆ fclose() [4/6]

fclose ( frate )

◆ fclose() [5/6]

fclose ( fratev )

◆ fclose() [6/6]

fclose ( fratem )

◆ for()

for ( )

pointers to detectors

Definition at line 79 of file cwb_report_prod_2.C.

◆ ForceStyle()

gROOT ForceStyle ( )

◆ fprintf() [1/19]

fprintf	(	flive	,
		"nonzero lags	live = `%10.1f \`,
		liveTot
	)

◆ fprintf() [2/19]

fprintf	(	ftrig	,
		"# correlation	threshold = `%f \`,
		T_cor
	)

◆ fprintf() [3/19]

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

◆ fprintf() [4/19]

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

◆ fprintf() [5/19]

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

◆ fprintf() [6/19]

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

◆ fprintf() [7/19]

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

◆ fprintf() [8/19]

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

◆ fprintf() [9/19]

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

◆ fprintf() [10/19]

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

◆ fprintf() [11/19]

fprintf	(	ftrig	,
		"# 8 - likelihood \
	)

◆ fprintf() [12/19]

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

◆ fprintf() [13/19]

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

◆ fprintf() [14/19]

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

◆ fprintf() [15/19]

fprintf	(	ftrig	,
		"# 12 - bandwidth \
	)

◆ fprintf() [16/19]

fprintf	(	ftrig	,
		"# 13 - duration \
	)

◆ fprintf() [17/19]

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

◆ fprintf() [18/19]

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

◆ fprintf() [19/19]

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

◆ GetHistogram()

gr2 GetHistogram ( ) -> SetXTitle("lag")

◆ GetXaxis()

xCor GetXaxis ( ) -> SetTitle("network correlation")

Definition at line 339 of file cwb_report_prod_2.C.

◆ GetYaxis()

xCor GetYaxis ( ) -> SetTitle("events")

Definition at line 570 of file cwb_report_cbc.C.

◆ if() [1/9]

if ( nIFO = =1 )

Definition at line 70 of file cwb_report_prod_2.C.

◆ if() [2/9]

if ( ! slagFound )

Definition at line 132 of file cwb_report_prod_2.C.

◆ if() [3/9]

if ( flive = =NULL )

Definition at line 241 of file cwb_report_prod_2.C.

◆ if() [4/9]

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

Definition at line 247 of file cwb_report_prod_2.C.

◆ if() [5/9]

if ( TString(analysis) = ="1G" )

Definition at line 407 of file cwb_report_prod_2.C.

◆ if() [6/9]

if ( pp_rho_vs_rate_no_multiplicity )

Definition at line 829 of file cwb_report_prod_2.C.

◆ if() [7/9]

if ( !out_far. good() )

Definition at line 853 of file cwb_report_prod_2.C.

◆ if() [8/9]

if ( !out_efar. good() )

Definition at line 857 of file cwb_report_prod_2.C.

◆ if() [9/9]

if ( delete c1 )

Definition at line 965 of file cwb_report_prod_2.C.

◆ live()

TChain live ( "liveTime" )

◆ next()

TIter next ( wave. GetListOfBranches() )

◆ open()

out_efar open	(	fname	,
		ios::out
	)

◆ out_lf()

ofstream out_lf	(	fname	,
		ios::out
	)

◆ printf()

printf	(	"total live time: non-zero	lags = `%10.1f \`,
		liveTot
	)

◆ Reset()

next Reset ( )

◆ SaveAs()

c1 SaveAs ( fname )

◆ SetBorderMode()

c1 SetBorderMode ( 0 )

◆ SetBorderSize()

c1 SetBorderSize ( 2 )

◆ SetBottomMargin()

c1 SetBottomMargin ( 0. 103939 )

◆ SetCanvasColor()

gStyle SetCanvasColor ( kWhite )

◆ SetEstimate()

wave SetEstimate ( nsize )

◆ SetFillColor() [1/2]

c1 SetFillColor ( 0 )

◆ SetFillColor() [2/2]

cutF SetFillColor ( 2 )

◆ SetFrameBorderMode()

gStyle SetFrameBorderMode ( 0 )

◆ SetGridx()

c1 SetGridx ( )

◆ SetGridy()

c1 SetGridy ( )

◆ SetLabelFont() [1/2]

gStyle SetLabelFont	(	42	,
		"X"
	)

◆ SetLabelFont() [2/2]

gStyle SetLabelFont	(	42	,
		"Y"
	)

◆ SetLabelOffset() [1/2]

gStyle SetLabelOffset	(	0.	010,
		"X"
	)

◆ SetLabelOffset() [2/2]

gStyle SetLabelOffset	(	0.	010,
		"Y"
	)

◆ SetLabelSize() [1/2]

gStyle SetLabelSize	(	0.	06,
		"X"
	)

◆ SetLabelSize() [2/2]

gStyle SetLabelSize	(	0.	06,
		"Y"
	)

◆ SetLineColor() [1/2]

gStyle SetLineColor ( kWhite )

◆ SetLineColor() [2/2]

legl SetLineColor ( 1 )

◆ SetLogx()

c1 SetLogx ( kFALSE )

◆ SetLogy() [1/2]

c1 SetLogy ( kFALSE )

◆ SetLogy() [2/2]

c1 SetLogy ( kTRUE )

◆ SetMarkerColor()

rhocc SetMarkerColor ( 1 )

◆ SetMarkerSize() [1/2]

rhocc SetMarkerSize ( 0. 5 )

◆ SetMarkerSize() [2/2]

gr2 SetMarkerSize ( 1 )

◆ SetMarkerStyle()

rhocc SetMarkerStyle ( 20 )

◆ SetNumberContours()

gStyle SetNumberContours ( 256 )

◆ SetOptStat() [1/2]

gStyle SetOptStat ( kFALSE )

◆ SetOptStat() [2/2]

gStyle SetOptStat ( kTRUE )

◆ SetPalette()

gStyle SetPalette	(	1	,
		0
	)

◆ SetRightMargin()

c1 SetRightMargin ( 0. 1517039 )

◆ SetStatBorderSize()

gStyle SetStatBorderSize ( 1 )

◆ SetStats()

cutF SetStats ( kFALSE )

◆ SetTextFont()

gStyle SetTextFont ( 12 )

◆ SetTitle() [1/7]

xCor SetTitle ( "" )

◆ SetTitle() [2/7]

McErr SetTitle ( "After pp cuts" )

◆ SetTitle() [3/7]

rhocc SetTitle ( )

◆ SetTitle() [4/7]

rho_pf SetTitle ( "chi2" )

◆ SetTitle() [5/7]

rho_mchirp SetTitle ( "rho vs Mchirp(after pp cuts)" )

◆ SetTitle() [6/7]

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

◆ SetTitle() [7/7]

gr2 SetTitle ( _title )

◆ 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. 060 )

◆ SetTitleOffset() [1/3]

gStyle SetTitleOffset	(	1.	4,
		"X"
	)

◆ SetTitleOffset() [2/3]

gStyle SetTitleOffset	(	2.	0,
		"Y"
	)

◆ SetTitleOffset() [3/3]

xCor SetTitleOffset	(	1.	3,
		"Y"
	)

◆ SetTitleSize() [1/2]

gStyle SetTitleSize	(	0.	05,
		"X"
	)

◆ SetTitleSize() [2/2]

gStyle SetTitleSize	(	0.	04,
		"Y"
	)

◆ SetTitleW()

gStyle SetTitleW ( 0. 95 )

◆ SetTitleY()

gStyle SetTitleY ( 0. 99 )

◆ SetTopMargin()

c1 SetTopMargin ( 0. 0772727 )

◆ sprintf() [1/23]

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

◆ sprintf() [2/23]

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

◆ sprintf() [3/23]

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

◆ sprintf() [4/23]

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

◆ sprintf() [5/23]

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

◆ sprintf() [6/23]

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

◆ sprintf() [7/23]

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

◆ sprintf() [8/23]

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

◆ sprintf() [9/23]

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

◆ sprintf() [10/23]

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

◆ sprintf() [11/23]

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

◆ sprintf() [12/23]

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

◆ sprintf() [13/23]

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

◆ sprintf() [14/23]

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

◆ sprintf() [15/23]

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

◆ sprintf() [16/23]

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

◆ sprintf() [17/23]

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

◆ sprintf() [18/23]

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

◆ sprintf() [19/23]

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

◆ sprintf() [20/23]

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

◆ sprintf() [21/23]

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

◆ sprintf() [22/23]

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

◆ sprintf() [23/23]

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

◆ Update()

c1 Update ( )

◆ wave()

TChain wave ( "waveburst" )

◆ while()

while ( (branch=(TBranch *) next()) )

Definition at line 128 of file cwb_report_prod_2.C.

Variable Documentation

◆ _title

char _title[1024]

Definition at line 174 of file cwb_report_prod_2.C.

◆ acor

double acor

Definition at line 167 of file cwb_report_prod_2.C.

◆ asnr

double asnr

Definition at line 167 of file cwb_report_prod_2.C.

◆ bcat2

bool bcat2 = false

Definition at line 106 of file cwb_report_prod_2.C.

◆ bcat3

bool bcat3 = false

Definition at line 108 of file cwb_report_prod_2.C.

◆ bhveto

bool bhveto = false

Definition at line 107 of file cwb_report_prod_2.C.

◆ branch

TBranch* branch

Definition at line 125 of file cwb_report_prod_2.C.

◆ c1

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

Definition at line 87 of file cwb_report_prod_2.C.

◆ ch1

char ch1[256]

Definition at line 175 of file cwb_report_prod_2.C.

◆ ch3

char ch3[256]

Definition at line 176 of file cwb_report_prod_2.C.

◆ Clear

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()

Definition at line 945 of file cwb_report_prod_2.C.

◆ countVETO

int countVETO[100]

Definition at line 140 of file cwb_report_prod_2.C.

◆ cutF

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

Definition at line 348 of file cwb_report_prod_2.C.

◆ dens

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

Definition at line 300 of file cwb_report_prod_2.C.

◆ ecor

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

Definition at line 318 of file cwb_report_prod_2.C.

◆ ECOR

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

Definition at line 324 of file cwb_report_prod_2.C.

◆ ECOR_cc

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

Definition at line 368 of file cwb_report_prod_2.C.

◆ ecor_cc

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

Definition at line 377 of file cwb_report_prod_2.C.

◆ Elag

Elag = Tlag

Definition at line 223 of file cwb_report_prod_2.C.

◆ else

else

Initial value:

{

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

pp_rho_min

double pp_rho_min

Definition: cwb_pparameters.C:10

rho_pf

TH2F * rho_pf

Definition: cwb_report_prod_2.C:406

pp_rho_max

double pp_rho_max

Definition: cwb_pparameters.C:11

Definition at line 411 of file cwb_report_prod_2.C.

◆ far_drho

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

Definition at line 202 of file cwb_report_prod_2.C.

◆ far_rho_bin

int far_rho_bin = float(far_rho_max-far_rho_min)/far_rho_wbin

Definition at line 201 of file cwb_report_prod_2.C.

◆ far_rho_max

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

Definition at line 199 of file cwb_report_prod_2.C.

◆ far_rho_min

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

Definition at line 198 of file cwb_report_prod_2.C.

◆ far_rho_wbin

float far_rho_wbin = 0.01

Definition at line 200 of file cwb_report_prod_2.C.

◆ filelags

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

Definition at line 818 of file cwb_report_prod_2.C.

◆ flive

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

Definition at line 240 of file cwb_report_prod_2.C.

◆ fname

char fname[1024]

Definition at line 173 of file cwb_report_prod_2.C.

◆ frate

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

Definition at line 873 of file cwb_report_prod_2.C.

◆ fratem

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

Definition at line 888 of file cwb_report_prod_2.C.

◆ fratev

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

Definition at line 880 of file cwb_report_prod_2.C.

◆ ftrig

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

Definition at line 564 of file cwb_report_prod_2.C.

◆ gps

double gps

Definition at line 212 of file cwb_report_prod_2.C.

◆ gr2

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

Definition at line 1041 of file cwb_report_prod_2.C.

◆ hF

TH1F* hF[NIFO_MAX+1]

Definition at line 509 of file cwb_report_prod_2.C.

◆ hF2

TH1F* hF2[NIFO_MAX+1]

Definition at line 510 of file cwb_report_prod_2.C.

◆ hpf

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

Definition at line 306 of file cwb_report_prod_2.C.

◆ hR

TH1F* hR[3]

Definition at line 545 of file cwb_report_prod_2.C.

◆ hR_col

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

Definition at line 546 of file cwb_report_prod_2.C.

◆ hvED

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

Definition at line 312 of file cwb_report_prod_2.C.

◆ i

i = 16

Definition at line 168 of file cwb_report_prod_2.C.

◆ IFO

char IFO[NIFO_MAX][32]

Definition at line 78 of file cwb_report_prod_2.C.

◆ iy_lag

delete [] iy_lag

Definition at line 745 of file cwb_report_prod_2.C.

◆ j

int j

Definition at line 168 of file cwb_report_prod_2.C.

◆ k

int k

Definition at line 168 of file cwb_report_prod_2.C.

◆ LAG_PASS

LAG_PASS =0

Definition at line 604 of file cwb_report_prod_2.C.

◆ lagslag

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 230 of file cwb_report_prod_2.C.

◆ Like

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

Definition at line 330 of file cwb_report_prod_2.C.

◆ Live

double Live

Definition at line 215 of file cwb_report_prod_2.C.

◆ liveTot

double liveTot = 0.

Definition at line 213 of file cwb_report_prod_2.C.

◆ liveZero

double liveZero = 0.

Definition at line 214 of file cwb_report_prod_2.C.

◆ m

int m

Definition at line 168 of file cwb_report_prod_2.C.

◆ max_lag

int max_lag =0

Definition at line 276 of file cwb_report_prod_2.C.

◆ McErr

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

Definition at line 342 of file cwb_report_prod_2.C.

◆ Mchirp

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

Definition at line 336 of file cwb_report_prod_2.C.

◆ min_lag

int min_lag =kMaxInt

Definition at line 275 of file cwb_report_prod_2.C.

◆ n

int n

Definition at line 168 of file cwb_report_prod_2.C.

◆ nameVETO

TString nameVETO[100]

Definition at line 141 of file cwb_report_prod_2.C.

◆ nlag

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

Definition at line 274 of file cwb_report_prod_2.C.

◆ nsize

int nsize = wave.GetEntries()

Definition at line 192 of file cwb_report_prod_2.C.

◆ ntype

size_t ntype = 0

Definition at line 172 of file cwb_report_prod_2.C.

◆ nXDQF

int nXDQF =0

Definition at line 144 of file cwb_report_prod_2.C.

◆ Olag

Olag = Tlag

Definition at line 224 of file cwb_report_prod_2.C.

◆ out_efar

ofstream out_efar

Definition at line 854 of file cwb_report_prod_2.C.

◆ out_far

ofstream out_far

Definition at line 850 of file cwb_report_prod_2.C.

◆ output

char output[256]

Definition at line 894 of file cwb_report_prod_2.C.

◆ pf_cc

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

Definition at line 423 of file cwb_report_prod_2.C.

◆ pf_ed

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

Definition at line 441 of file cwb_report_prod_2.C.

◆ pf_vED

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

Definition at line 432 of file cwb_report_prod_2.C.

◆ rho

double rho

Definition at line 167 of file cwb_report_prod_2.C.

◆ rho_L

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

Definition at line 450 of file cwb_report_prod_2.C.

◆ rho_max

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

Definition at line 196 of file cwb_report_prod_2.C.

◆ rho_mchirp

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

Definition at line 459 of file cwb_report_prod_2.C.

◆ rho_min

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

Definition at line 195 of file cwb_report_prod_2.C.

◆ rho_pf

TH2F* rho_pf

Definition at line 406 of file cwb_report_prod_2.C.

◆ rho_subnet

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

Definition at line 386 of file cwb_report_prod_2.C.

◆ rho_vED

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

Definition at line 396 of file cwb_report_prod_2.C.

◆ rhocc

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

Definition at line 358 of file cwb_report_prod_2.C.

◆ Rlag

Rlag = Tlag

Definition at line 222 of file cwb_report_prod_2.C.

◆ rUn

Int_t rUn

Definition at line 211 of file cwb_report_prod_2.C.

◆ save

bool save

Definition at line 169 of file cwb_report_prod_2.C.

◆ SAVE

SAVE =0

Definition at line 605 of file cwb_report_prod_2.C.

◆ save_veto

bool save_veto

Definition at line 170 of file cwb_report_prod_2.C.

◆ saveVETO

bool saveVETO[100]

Definition at line 138 of file cwb_report_prod_2.C.

◆ slagFound

bool slagFound =false

Definition at line 126 of file cwb_report_prod_2.C.

◆ sm_theta_phi

skymap sm_theta_phi(3)

◆ sTARt

double sTARt

Definition at line 212 of file cwb_report_prod_2.C.

◆ sTOp

double sTOp

Definition at line 212 of file cwb_report_prod_2.C.

◆ Tdlag

wavearray<double> Tdlag

Definition at line 220 of file cwb_report_prod_2.C.

◆ Tlag

wavearray<double> Tlag

Definition at line 221 of file cwb_report_prod_2.C.

◆ Trun

Trun = 0.

Definition at line 217 of file cwb_report_prod_2.C.

◆ vED

double vED

Definition at line 167 of file cwb_report_prod_2.C.

◆ VETO

UChar_t VETO[100]

Definition at line 139 of file cwb_report_prod_2.C.

◆ wave

netevent W& wave

Definition at line 103 of file cwb_report_prod_2.C.

◆ Wlag

wavearray<double> Wlag[NIFO_MAX+1]

Definition at line 218 of file cwb_report_prod_2.C.

◆ Wsel

wavearray<int> Wsel(ntrg)

◆ Wslag

wavearray<double> Wslag[NIFO_MAX+1]

Definition at line 219 of file cwb_report_prod_2.C.

◆ xCOR

double xCOR

Definition at line 167 of file cwb_report_prod_2.C.

◆ xcor

double xcor

Definition at line 167 of file cwb_report_prod_2.C.

◆ xCor

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

Definition at line 294 of file cwb_report_prod_2.C.

◆ xcor2

double xcor2

Definition at line 167 of file cwb_report_prod_2.C.

◆ xcor3

double xcor3

Definition at line 167 of file cwb_report_prod_2.C.

◆ Xcut

wavearray<double> Xcut(pp_rho_bin)

◆ XDQF

dqfile* XDQF = new dqfile[nVDQF]

Definition at line 145 of file cwb_report_prod_2.C.

◆ Xerr

Xerr = 0.

Definition at line 181 of file cwb_report_prod_2.C.

◆ Xfar

wavearray<double> Xfar(far_rho_bin)

◆ xlag

float xlag

Definition at line 210 of file cwb_report_prod_2.C.

◆ xlive

double xlive

Definition at line 212 of file cwb_report_prod_2.C.

◆ xrun

Int_t xrun

Definition at line 211 of file cwb_report_prod_2.C.

◆ Ycut

Ycut = 0.

Definition at line 179 of file cwb_report_prod_2.C.

◆ yCUT

yCUT = 0.

Definition at line 180 of file cwb_report_prod_2.C.

◆ Ycut_multi

Ycut_multi = 0.

Definition at line 186 of file cwb_report_prod_2.C.

◆ Ycut_veto

Ycut_veto = 0.

Definition at line 184 of file cwb_report_prod_2.C.

◆ Yerr

wavearray<double> Yerr(pp_rho_bin)

◆ yERR

wavearray<double> yERR(pp_rho_bin)

◆ Yerr_multi

wavearray<double> Yerr_multi(pp_rho_bin)

◆ Yerr_veto

wavearray<double> Yerr_veto(pp_rho_bin)

◆ Yfar

Yfar = 0.

Definition at line 205 of file cwb_report_prod_2.C.

Macros

Functions

Variables

Macro Definition Documentation

◆ YEAR

Function Documentation

◆ Add() [1/2]

◆ Add() [2/2]

◆ Clear()

◆ close()

◆ Divide()

◆ Draw() [1/4]

◆ Draw() [2/4]

◆ Draw() [3/4]

◆ Draw() [4/4]

◆ Exec() [1/2]

◆ Exec() [2/2]

◆ fclose() [1/6]

◆ fclose() [2/6]

◆ fclose() [3/6]

◆ fclose() [4/6]

◆ fclose() [5/6]

◆ fclose() [6/6]

◆ for()

◆ ForceStyle()

◆ fprintf() [1/19]

◆ fprintf() [2/19]

◆ fprintf() [3/19]

◆ fprintf() [4/19]

◆ fprintf() [5/19]

◆ fprintf() [6/19]

◆ fprintf() [7/19]

◆ fprintf() [8/19]

◆ fprintf() [9/19]

◆ fprintf() [10/19]

◆ fprintf() [11/19]

◆ fprintf() [12/19]

◆ fprintf() [13/19]

◆ fprintf() [14/19]

◆ fprintf() [15/19]

◆ fprintf() [16/19]

◆ fprintf() [17/19]

◆ fprintf() [18/19]

◆ fprintf() [19/19]

◆ GetHistogram()

◆ GetXaxis()

◆ GetYaxis()

◆ if() [1/9]

◆ if() [2/9]

◆ if() [3/9]

◆ if() [4/9]

◆ if() [5/9]

◆ if() [6/9]

◆ if() [7/9]

◆ if() [8/9]

◆ if() [9/9]

◆ live()

◆ next()

◆ open()

◆ out_lf()

◆ printf()

◆ Reset()

◆ SaveAs()

◆ SetBorderMode()

◆ SetBorderSize()

◆ SetBottomMargin()

◆ SetCanvasColor()

◆ SetEstimate()

◆ SetFillColor() [1/2]

◆ SetFillColor() [2/2]

◆ SetFrameBorderMode()

◆ SetGridx()

◆ SetGridy()

◆ SetLabelFont() [1/2]

◆ SetLabelFont() [2/2]

◆ SetLabelOffset() [1/2]

◆ SetLabelOffset() [2/2]

◆ SetLabelSize() [1/2]

◆ SetLabelSize() [2/2]

◆ SetLineColor() [1/2]