cbc_plots.C File Reference

Macros
#define	CONTOURS   20
#define	LIV_FILE_NAME   "live.txt"
#define	NCont   99
#define	NRGBs   6
#define	RHO_BIN   0.1
#define	RHO_MIN   6.0
#define	RHO_NBINS   1000

Functions
sim	Add (sim_file_name)
mdc	Add (mdc_file_name)
leg_snr	AddEntry ("", lab,"a")
leg_snr	AddEntry (htemp, lab,"l")
leg_snr	AddEntry (htemp2, lab,"l")
leg_snr	AddEntry (htemp3, lab,"l")
c1	Clear ()
efficiency_first_shell	Divide (factor_events_inj[nfactor-1])
cout<< endl;fclose(fev);cout << "Recovered entries: "<< cnt << endl;cout<< "Recovered entries: "<< cnt2<< endl;cout << "Recovered entries vetoed: " << countv<< endl;cout<<"dV : " << dV<< " dV1 : "<< dV1<< endl;for(int i=0;i< RHO_NBINS;i++){eVrho[i]=TMath::Sqrt(eVrho[i]);}cout << "Vrho[0] = "<< Vrho[0] <<" +/- "<< eVrho[0]<< endl;cout << "Vrho[RHO_NBINS-1] = " << Vrho[RHO_NBINS-1] <<" +/- "<< eVrho[RHO_NBINS-1] << endl;inj_events->	Draw ("colz")
p_inj	Draw ()
mdc	Draw (sel,"")
sim	Draw (sel, newcut,"")
sim	Draw (sel, newcut2,"same")
sim	Draw (sel, newcut2)
	exit (1)
	exit (0)
	fclose (finj)
D_Mtot_rec	Fill (mass[1]+mass[0], range[1])
D_Mchirp_rec	Fill (chirp[0], range[1])
D_q_rec	Fill (mass[0]/mass[1], range[1])
rhocc	Fill (netcc[0], rho[pp_irho])
rho_pf	Fill (chi2, rho[pp_irho])
rec_events	Fill (mass[1], mass[0])
	for (int i=0;i< nfactor;i++)
	for (int i=1;i<=NBINS_mass;i++)
gROOT	ForceStyle ()
sim	GetHistogram ()
inj_events	GetXaxis () -> SetRangeUser(MIN_plot_mass1, MAX_plot_mass1)
inj_events	GetYaxis () -> SetRangeUser(MIN_plot_mass2, MAX_plot_mass2)
inj_events	GetZaxis () -> SetTickLength(0.01)
	if (strlen(ifo[0]), 0)
	if (MDC.GetInspiralOption("--waveform")!="")
	if (MDC.GetInspiralOption("--min-mtotal")!="")
	if (MDC.GetInspiralOption("--max-mtotal")!="")
	if (MDC.GetInspiralOption("--min-distance")!="")
	if (MDC.GetInspiralOption("--max-distance")!="")
	if (MDC.GetInspiralOption("--min-mratio")!="")
	if (MDC.GetInspiralOption("--max-mratio")!="")
	if (MDC.GetInspiralOption("--d-distr")!="")
	if (bminMtot)
	if (bmaxMtot)
	if (bminDistance)
	if (bmaxDistance)
	if (strlen(veto_not_vetoed)==0)
	if (range[1]==0.0)
	if (l >NBINS_mass2+1)
	if (m >NBINS_mass+1)
	if (rho[pp_irho]<=T_out)
	if (if() c1->SetLogy(kTRUE) pp_rho_log)
else	if (DDistrUniform)
gROOT	Macro (configPlugin.GetTitle())
TChain	mdc ("mdc")
c1	Range (-1.216392,-477.6306, 508.8988, 2814.609)
c1	SaveAs (fname)
c1	SetBorderMode (0)
c1	SetBorderSize (2)
p_inj	SetBorderSize (0)
c1	SetBottomMargin (0.1450777)
mdc	SetBranchAddress ("time", time)
mdc	SetBranchAddress ("mass", mass)
mdc	SetBranchAddress ("factor",&factor)
mdc	SetBranchAddress ("distance",&distance)
mdc	SetBranchAddress ("mchirp",&mchirp)
mdc	SetBranchAddress ("spin", spin)
gStyle	SetCanvasColor (kWhite)
inj_events	SetContour (NCont)
h_radius	SetEntries (1)
c1	SetFillColor (0)
htemp2	SetFillColor (kRed)
htemp3	SetFillColor (kBlue)
htemp2	SetFillStyle (3017)
gStyle	SetFrameBorderMode (0)
c1	SetGridx ()
c1	SetGridy ()
htemp	SetLineColor (kBlack)
htemp2	SetLineColor (kRed)
htemp3	SetLineColor (kBlue)
htemp	SetLineWidth (3)
htemp2	SetLineWidth (2)
c1	SetLogx (false)
c1	SetLogy (1)
D_Mtot_inj	SetMarkerColor (2)
D_Mtot_rec	SetMarkerColor (4)
rhocc	SetMarkerColor (1)
sim	SetMarkerColor (kRed)
sim	SetMarkerColor (kBlue)
D_Mtot_inj	SetMarkerSize (0.5)
D_Mtot_inj	SetMarkerStyle (20)
D_Mtot_inj	SetName ("D_Mtotinj")
D_Mchirp_inj	SetName ("D_Chirp_inj")
D_q_inj	SetName ("D_q_inj")
gStyle	SetNumberContours (256)
gStyle	SetOptFit (1)
gStyle	SetOptStat (kFALSE)
cout<< "Selection: "<< sel << endl;gStyle->	SetOptStat (1)
c1	SetRightMargin (0.1154618)
gStyle	SetStatBorderSize (1)
rhocc	SetStats (kFALSE)
p_inj	SetTextColor (1)
p_inj	SetTextFont (32)
p_inj	SetTextSize (0.045)
inj_events	SetTitle ("")
rhocc	SetTitle ("0 < cc < 1")
rho_pf	SetTitle ("chi2")
cout<< "nwave_final : " << nwave_final<< endl;sprintf(title,"%s", newcut);sprintf(title,"%s", ptitle.Data());htemp->	SetTitle (title)
cout<< "cbc_plots.C starts..." << endl;cout<< "Mass1 : [" << min_mass1<<","<< max_mass1 <<"] with "<< NBINS_mass1 <<" bins"<< endl;cout<< "Mass2 : [" << min_mass2<<","<< max_mass2 <<"] with "<< NBINS_mass2 <<" bins"<< endl;CWB::Toolbox TB;CWB::CBCTool cbcTool;TB.checkFile(gSystem-> Getenv("CWB_ROOTLOGON_FILE"));TB.checkFile(gSystem-> Getenv("CWB_PARAMETERS_FILE"));TB.checkFile(gSystem-> Getenv("CWB_UPARAMETERS_FILE"));TB.checkFile(gSystem-> Getenv("CWB_PPARAMETERS_FILE"));TB.checkFile(gSystem-> Getenv("CWB_UPPARAMETERS_FILE"));TB.checkFile(gSystem-> Getenv("CWB_EPPARAMETERS_FILE"));TB.mkDir(netdir, true);gStyle->	SetTitleFillColor (kWhite)
rhocc	SetTitleOffset (1.3,"Y")
c1	SetTopMargin (0.07642487)
TChain	sim ("waveburst")
	sprintf (inj_title,"Injected events")
	sprintf (fname,"%s/Injected_mass1_mass2.eps", netdir)
	sprintf (fname,"%s/Delta_t.png", netdir)
	sprintf (sel,"%s)>>hist(500)", sel)
	sprintf (newcut,"(((time[0]-time[%d])>-%g) || (time[0]-time[%d])<%g) && rho[%d]> %g", nIFO, T_win, nIFO, 2 *T_win, pp_irho, T_out)
	sprintf (sel,"sqrt(iSNR[%d]", 0)
	sprintf (lab,"Injections Average SNR: %g", htemp->GetMean())
	sprintf (lab,"Injected: %i", nmdc)
	sprintf (lab,"Found(minimal cuts): %i", nwave)
	sprintf (lab,"Found(final cuts): %i", nwave_final)
	sprintf (fname,"%s/Injected_snr_distributions.png", netdir)
	sprintf (fname,"%s/Estimated_snr_vs_Injected_snr.eps", netdir)
	sprintf (radius_title,"%s : Effective radius (Mpc)", networkname)
	sprintf (fname,"%s/Effective_radius.png", netdir)
	strcpy (cfg->tmp_dir,"tmp")
c1	Update ()

Variables
double	blue [NRGBs] = { 0.97, 0.97, 0.00, 0.00, 0.00, 0.00 }
bool	bmaxDistance
bool	bmaxMtot
bool	bmaxRatio
bool	bminDistance
bool	bminMtot
bool	bminRatio
TCanvas *	c1 = new TCanvas("c1","c1",3,47,1000,802)
CWB::config *	cfg = new CWB::config
float	chi [3]
float	chi2 = penalty>0 ? log10(penalty) : 0
cout<< "Shell volume: " << shell_volume<< ", internal sphere volume: " << volume_internal_sphere << endl;TH2F *efficiency_first_shell=(TH2F *) factor_events_rec->	Clone ()
int	countv =0
int	cz
TH2F *	D_Mchirp_inj = new TH2F("Distance vs MChirp inj.","",1000,MINCHIRP,MAXCHIRP,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000./FACTORS[0])
TH2F *	D_Mchirp_rec = new TH2F("Distance vs MChirp rec.","",1000,MINCHIRP,MAXCHIRP,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000/FACTORS[0])
TH2F *	D_Mtot_inj = new TH2F("Distance vs Mtot inj.","",1000,MINMtot,MAXMtot,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000./FACTORS[0])
TH2F *	D_Mtot_rec = new TH2F("Distance vs Mtot rec.","",1000,MINMtot,MAXMtot,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000/FACTORS[0])
TH2F *	D_q_inj = new TH2F("Distance vs q inj.","",1000,MINRATIO,MAXRATIO,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000./FACTORS[0])
TH2F *	D_q_rec = new TH2F("Distance vs q rec.","",1000,MINRATIO,MAXRATIO,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000/FACTORS[0])
bool	DDistrUniform
bool	DDistrVolume
float	distance
TH1F *	Dt = new TH1F("Dt","",1000,-0.5,0.5)
	else {float MINMtot = 0.0
	error_volume [m][l] = pow(dV,2)
TExec *	ex2 = new TExec("ex2","gStyle->SetPaintTextFormat(\".0f\");")
TF1 *	f2 = cbcTool.doRangePlot(RHO_NBINS, Trho, Rrho, eRrho, RHO_MIN, T_out, c1, networkname, netdir, write_ascii)
float	factor
TH2F *	factor_events_inj [nfactor]
TH2F *	factor_events_rec = new TH2F("factor_events_rec","",NBINS_mass,MIN_MASS,MAX_MASS,NBINS_mass2,min_mass2,max_mass2)
cout<< "Injected signals: " << mdc.GetEntries()<< endl;char cut[512];sprintf(cut,"rho[%d]> %f && %s", pp_irho, T_out, ch2);float ecor, m1, m2, netcc[3], neted, penalty;float rho[2];float chirp[6];float range[2];float iSNR[3], sSNR[3];sim.SetBranchAddress("mass", mass);sim.SetBranchAddress("factor",&factor);sim.SetBranchAddress("range", range);sim.SetBranchAddress("chirp", chirp);sim.SetBranchAddress("rho", rho);sim.SetBranchAddress("netcc", netcc);sim.SetBranchAddress("neted",&neted);sim.SetBranchAddress("ecor",&ecor);sim.SetBranchAddress("penalty",&penalty);sim.SetBranchAddress("time", time);sim.SetBranchAddress("iSNR", iSNR);sim.SetBranchAddress("sSNR", sSNR);sim.SetBranchAddress("spin", spin);float volume[NBINS_mass][NBINS_mass2], error_volume[NBINS_mass][NBINS_mass2];float volume_first_shell[NBINS_mass][NBINS_mass2], error_volume_first_shell[NBINS_mass][NBINS_mass2];float radius[NBINS_mass][NBINS_mass2], error_radius[NBINS_mass][NBINS_mass2];int cnt=0;int cnt2=0;for(int i=0;i < NBINS_mass;i++){for(int j=0;j < NBINS_mass2;j++){volume[i][j]=0.;error_volume[i][j]=0.;volume_first_shell[i][j]=0.;error_volume_first_shell[i][j]=0.;radius[i][j]=0.;error_radius[i][j]=0.;}}double Vrho[RHO_NBINS], eVrho[RHO_NBINS], Rrho[RHO_NBINS], eRrho[RHO_NBINS], Trho[RHO_NBINS];for(int i=0;i < RHO_NBINS;i++){Vrho[i]=0.;eVrho[i]=0.;Rrho[i]=0.;eRrho[i]=0.;Trho[i]=RHO_MIN+i *RHO_BIN;}double dV, dV1, dV_spin_mtot;int nT;for(int g=0;g <(int) sim.GetEntries();g++){sim.GetEntry(g);if(rho[pp_irho] <=RHO_MIN){continue;}if(netcc[0] <=T_cor){continue;}if((time[0]-time[nIFO]) <-T_win||(time[0]-time[nIFO])> 2 *T_win){continue;}if(T_vED > 0){if(neted/ecor >=T_vED){continue;}}if(T_pen > 0){if(penalty<=T_vED){continue;}}if(++cnt%1000==0){cout << cnt<< " - ";}Dt->	Fill (time[0]-time[nIFO])
char	fname [1024]
cout<< "nmdc : "<< nmdc<< endl;TH2F *htemp=(TH2F *) gPad->	GetPrimitive ("hist")
double	green [NRGBs] = { 0.97, 0.30, 0.40, 0.97, 0.00, 0.00 }
TH2F *	h_radius = new TH2F("h_radius","",NBINS_mass,MIN_MASS,MAX_MASS,NBINS_mass2,min_mass2,max_mass2)
TH2F *	htemp4 = (TH2F*)gPad->GetPrimitive("hist4")
TH2F *	htemp5 = (TH2F*)gPad->GetPrimitive("hist5")
TH2F *	inj_events = new TH2F("inj_events","D_Minj",NBINS_mass,MIN_MASS,MAX_MASS,NBINS_mass2,min_mass2,max_mass2)
char	inj_title [256]
TString	Insp = MDC.GetInspiral()
int	l = TMath::FloorNint((m2-min_mass2)/MASS_BIN)
char	lab [256]
	leg_snr = new TLegend(0.6,0.755,0.885,0.923,"","brNDC")
int	m = TMath::FloorNint((m1-MIN_MASS)/MASS_BIN)
	m1 =mass[1]
	m2 =mass[0]
float	mass [2]
int	massbins = 0
int	MAX_AXIS_Z = inj_events->GetBinContent(inj_events->GetMaximumBin()) + 1
float	max_mass1 = MAX_MASS
float	max_mass2 = MAX_MASS
float	MAXCHIRP = 1.1*pow(max_mass1*max_mass2,3./5.)/pow(minMtot,1./5.)
float	maxDistance
float	maxMtot
float	maxRatio
float	MAXRATIO = 1.
float	mchirp
float	min_mass1 = MIN_MASS
float	min_mass2 = MIN_MASS
float	MINCHIRP = 0.9*pow(min_mass1*min_mass2,3./5.)/pow(maxMtot,1./5.)
float	minDistance
float	minMtot
float	minRatio
float	MINRATIO = 0.0
int	mt
int	NBINS_mass = (int)((MAX_MASS-MIN_MASS)/MASS_BIN)
int	NBINS_mass1 = (int)((max_mass1-min_mass1)/MASS_BIN)
int	NBINS_mass2 = (int)((max_mass2-min_mass2)/MASS_BIN)
network *	net =NULL
char	networkname [256]
int	nmdc = mdc.GetSelectedRows()
	nT =TMath::Min(TMath::Floor((rho[pp_irho]-RHO_MIN)/RHO_BIN),(double)RHO_NBINS)+1
int	nwave = sim.GetSelectedRows()
int	nwave_final = sim.GetSelectedRows()
TPaveText *	p_inj = new TPaveText(0.325301,0.926166,0.767068,0.997409,"blNDC")
TPaveText *	p_radius = new TPaveText(0.325301,0.926166,0.767068,0.997409,"blNDC")
char	radius_title [256]
TH2F *	rec_events = new TH2F("rec_events","D_Mrec",NBINS_mass,MIN_MASS,MAX_MASS,NBINS_mass2,min_mass2,max_mass2)
double	red [NRGBs] = { 0.00, 0.00, 0.00, 0.97, 0.97, 0.10 }
TH2F *	rho_pf = new TH2F("rho_pf","",100,-1.,2.,100,pp_rho_min,pp_rho_max)
TH2F *	rhocc = new TH2F("rhocc","",100,0.,1.,100,pp_rho_min,pp_rho_max)
float	spin [6]
double	stops [NRGBs] = { 0.10, 0.25, 0.45, 0.60, 0.75, 1.00 }
TText *	text = p_inj->AddText(inj_title)
double	time [6]
double	V0 = 0.0
	volume [m][l] = dV
TString	waveform
bool	write_ascii = false

Macro Definition Documentation

#define CONTOURS   20

Definition at line 5 of file cbc_plots.C.

#define LIV_FILE_NAME   "live.txt"

Definition at line 8 of file cbc_plots.C.

#define NCont   99

#define NRGBs   6

#define RHO_BIN   0.1

Definition at line 3 of file cbc_plots.C.

#define RHO_MIN   6.0

Definition at line 2 of file cbc_plots.C.

#define RHO_NBINS   1000

Definition at line 4 of file cbc_plots.C.

Function Documentation

sim Add

(

sim_file_name

)

mdc Add

(

mdc_file_name

)

leg_snr AddEntry	(	""	,
		lab	,
		"a"
	)

leg_snr AddEntry	(	htemp	,
		lab	,
		"l"
	)

leg_snr AddEntry	(	htemp2	,
		lab	,
		"l"
	)

leg_snr AddEntry	(	htemp3	,
		lab	,
		"l"
	)

c1 Clear

(

)

Definition at line 1261 of file CWB_Plugin_xWRC.C.

efficiency_first_shell Divide

(

factor_events_inj

[nfactor-1]

)

cout<< endl;fclose(fev); cout << "Recovered entries: " << cnt << endl; cout << "Recovered entries: " << cnt2 << endl; cout<< "Recovered entries vetoed: "<<countv<<endl; cout <<"dV : "<<dV<< " dV1 : "<<dV1<<endl;for(int i=0; i<RHO_NBINS; i++){eVrho[i]=TMath::Sqrt(eVrho[i]);}cout << "Vrho[0] = "<<Vrho[0] <<" +/- "<<eVrho[0]<<endl;cout<< "Vrho[RHO_NBINS-1] = "<<Vrho[RHO_NBINS-1] <<" +/- "<<eVrho[RHO_NBINS-1] <<endl; inj_events-> Draw

(

"colz"

)

p_inj Draw

(

)

mdc Draw	(	sel	,
		""
	)

sim Draw	(	sel	,
		newcut	,
		""
	)

sim Draw	(	sel	,
		newcut2	,
		"same"
	)

sim Draw	(	sel	,
		newcut2
	)

exit

(

1

)

exit

(

0

)

fclose

(

finj

)

D_Mtot_rec Fill	(	mass+	mass[1][0],
		range	[1]
	)

D_Mchirp_rec Fill	(	chirp	[0],
		range	[1]
	)

D_q_rec Fill	(	mass/	mass[0][1],
		range	[1]
	)

rhocc Fill	(	netcc	[0],
		rho	[pp_irho]
	)

rho_pf Fill	(	chi2	,
		rho	[pp_irho]
	)

rec_events Fill	(	mass	[1],
		mass	[0]
	)

for

(

)

Definition at line 263 of file cbc_plots.C.

for

(

int

i = 1; i<=NBINS_mass; i++

)

Definition at line 1345 of file cbc_plots.C.

gROOT ForceStyle

(

)

sim GetHistogram

(

)

Definition at line 913 of file cbc_plots.C.

inj_events GetXaxis

(

)

-> SetRangeUser(MIN_plot_mass1, MAX_plot_mass1)

Definition at line 212 of file cbc_plots.C.

inj_events GetYaxis

(

)

-> SetRangeUser(MIN_plot_mass2, MAX_plot_mass2)

Definition at line 213 of file cbc_plots.C.

inj_events GetZaxis

(

)

-> SetTickLength(0.01)

if	(	strlen(ifo[0])	,
		0
	)

Definition at line 96 of file cbc_plots.C.

if

(

MDC.GetInspiralOption("--waveform")!

= ""

)

Definition at line 118 of file cbc_plots.C.

if

(

MDC.GetInspiralOption("--min-mtotal")!

= ""

)

Definition at line 119 of file cbc_plots.C.

if

(

MDC.GetInspiralOption("--max-mtotal")!

= ""

)

Definition at line 120 of file cbc_plots.C.

if

(

MDC.GetInspiralOption("--min-distance")!

= ""

)

Definition at line 121 of file cbc_plots.C.

if

(

MDC.GetInspiralOption("--max-distance")!

= ""

)

Definition at line 122 of file cbc_plots.C.

if

(

MDC.GetInspiralOption("--min-mratio")!

= ""

)

Definition at line 123 of file cbc_plots.C.

if

(

MDC.GetInspiralOption("--max-mratio")!

= ""

)

Definition at line 124 of file cbc_plots.C.

if

(

MDC.GetInspiralOption("--d-distr")!

= ""

)

Definition at line 125 of file cbc_plots.C.

if

(

bminMtot

)

Definition at line 165 of file cbc_plots.C.

if

(

bmaxMtot

)

Definition at line 168 of file cbc_plots.C.

if

(

bminDistance

)

Definition at line 176 of file cbc_plots.C.

if

(

bmaxDistance

)

Definition at line 179 of file cbc_plots.C.

if

(

strlen(veto_not_vetoed)

= = 0

)

Definition at line 190 of file cbc_plots.C.

if

(

range

[1] = =0.0

)

Definition at line 596 of file cbc_plots.C.

if	(	l	,
		NBINS_mass2+	1
	)

Definition at line 609 of file cbc_plots.C.

if	(	m	,
		NBINS_mass+	1
	)

Definition at line 611 of file cbc_plots.C.

if

(

rho<=

T_out[pp_irho]

)

Definition at line 670 of file cbc_plots.C.

if

(

if() c1->SetLogy(kTRUE)

pp_rho_log

)

Definition at line 746 of file cbc_plots.C.

else if

(

DDistrUniform

)

Definition at line 1281 of file cbc_plots.C.

gROOT Macro

(

configPlugin.

GetTitle()

)

TChain mdc

(

"mdc"

)

c1 Range	(	-1.	216392,
		-477.	6306,
		508.	8988,
		2814.	609
	)

c1 SaveAs

(

fname

)

c1 SetBorderMode

(

0

)

c1 SetBorderSize

(

2

)

p_radius SetBorderSize

(

0

)

c1 SetBottomMargin

(

0.

1450777

)

mdc SetBranchAddress	(	"time"	,
		time
	)

mdc SetBranchAddress	(	"mass"	,
		mass
	)

mdc SetBranchAddress	(	"factor"	,
		&	factor
	)

mdc SetBranchAddress	(	"distance"	,
		&	distance
	)

mdc SetBranchAddress	(	"mchirp"	,
		&	mchirp
	)

mdc SetBranchAddress	(	"spin"	,
		spin
	)

gStyle SetCanvasColor

(

kWhite

)

h_radius SetContour

(

NCont

)

h_radius SetEntries

(

1

)

c1 SetFillColor

(

0

)

htemp2 SetFillColor

(

kRed

)

htemp3 SetFillColor

(

kBlue

)

q_rec SetFillStyle

(

3017

)

gStyle SetFrameBorderMode

(

0

)

c1 SetGridx

(

)

c1 SetGridy

(

)

htemp SetLineColor

(

kBlack

)

htemp2 SetLineColor

(

kRed

)

htemp3 SetLineColor

(

kBlue

)

htemp SetLineWidth

(

3

)

htemp2 SetLineWidth

(

2

)

c1 SetLogx

(

false

)

c1 SetLogy

(

1

)

D_Mtot_inj SetMarkerColor

(

2

)

D_Mtot_rec3 SetMarkerColor

(

4

)

rhocc SetMarkerColor

(

1

)

sim SetMarkerColor

(

kRed

)

sim SetMarkerColor

(

kBlue

)

D_Mtot_inj SetMarkerSize

(

0.

5

)

D_Mtot_inj SetMarkerStyle

(

20

)

D_Mtot_inj SetName

(

"D_Mtotinj"

)

D_Mchirp_inj SetName

(

"D_Chirp_inj"

)

D_q_inj SetName

(

"D_q_inj"

)

gStyle SetNumberContours

(

256

)

gStyle SetOptFit

(

1

)

gStyle SetOptStat

(

kFALSE

)

cout<< "Selection: "<< sel <<endl; gStyle-> SetOptStat

(

1

)

c1 SetRightMargin

(

0.

1154618

)

gStyle SetStatBorderSize

(

1

)

rhocc SetStats

(

kFALSE

)

p_radius SetTextColor

(

1

)

p_radius SetTextFont

(

32

)

p_inj SetTextSize

(

0.

045

)

inj_events SetTitle

(

""

)

rhocc SetTitle

(

)

rho_pf SetTitle

(

"chi2"

)

cout<< "nwave_final : " << nwave_final << endl; sprintf(title,"%s",newcut); sprintf(title,"%s",ptitle.Data()); htemp-> SetTitle

(

title

)

cout<< "cbc_plots.C starts..."<<endl; cout << "Mass1 : ["<<min_mass1<<","<<max_mass1<<"] with "<<NBINS_mass1<<" bins"<<endl; cout << "Mass2 : ["<<min_mass2<<","<<max_mass2<<"] with "<<NBINS_mass2<<" bins"<<endl; CWB::Toolbox TB; CWB::CBCTool cbcTool; TB.checkFile(gSystem->Getenv("CWB_ROOTLOGON_FILE")); TB.checkFile(gSystem->Getenv("CWB_PARAMETERS_FILE")); TB.checkFile(gSystem->Getenv("CWB_UPARAMETERS_FILE")); TB.checkFile(gSystem->Getenv("CWB_PPARAMETERS_FILE")); TB.checkFile(gSystem->Getenv("CWB_UPPARAMETERS_FILE")); TB.checkFile(gSystem->Getenv("CWB_EPPARAMETERS_FILE")); TB.mkDir(netdir,true); gStyle-> SetTitleFillColor

(

kWhite

)

rhocc SetTitleOffset	(	1.	3,
		"Y"
	)

c1 SetTopMargin

(

0.

07642487

)

TChain sim

(

"waveburst"

)

sprintf	(	inj_title	,
		"Injected events"
	)

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

sprintf	(	fname	,
		"%s/Delta_t.png"	,
		netdir
	)

sprintf	(	sel	,
		"%	s
	)

sprintf	(	newcut	,
		"(((time[0]-time[%d])>-%g) || (time[0]-time[%d])<%g) &&	rho[%d],
		%g"	,
		nIFO	,
		T_win	,
		nIFO	,
		2 *	T_win,
		pp_irho	,
		T_out
	)

sprintf

(

sel

)

sprintf	(	lab	,
		"Injections Average SNR: %g"	,
		htemp->	GetMean()
	)

sprintf	(	lab	,
		"Injected: %i"	,
		nmdc
	)

sprintf	(	lab	,
		"Found(minimal cuts): %i"	,
		nwave
	)

sprintf	(	lab	,
		"Found(final cuts): %i"	,
		nwave_final
	)

sprintf	(	fname	,
		"%s/Injected_snr_distributions.png"	,
		netdir
	)

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

sprintf	(	radius_title	,
		"%s : Effective radius (Mpc)"	,
		networkname
	)

sprintf	(	fname	,
		"%s/Effective_radius.png"	,
		netdir
	)

strcpy	(	cfg->	tmp_dir,
		"tmp"
	)

c1 Update

(

)

Variable Documentation

double blue[NRGBs] = { 0.97, 0.97, 0.00, 0.00, 0.00, 0.00 }

Definition at line 26 of file cbc_plots.C.

bool bmaxDistance

Definition at line 114 of file cbc_plots.C.

bool bmaxMtot

Definition at line 114 of file cbc_plots.C.

bool bmaxRatio

Definition at line 114 of file cbc_plots.C.

bool bminDistance

Definition at line 114 of file cbc_plots.C.

bool bminMtot

Definition at line 114 of file cbc_plots.C.

bool bminRatio

Definition at line 114 of file cbc_plots.C.

TCanvas* c1 = new TCanvas("c1","c1",3,47,1000,802)

Definition at line 197 of file cbc_plots.C.

CWB::config* cfg = new CWB::config

Definition at line 106 of file cbc_plots.C.

float chi[3]

Definition at line 439 of file cbc_plots.C.

float chi2 = penalty>0 ? log10(penalty) : 0

Definition at line 603 of file cbc_plots.C.

cout<< "Shell volume: " << shell_volume << ", internal sphere volume: " << volume_internal_sphere << endl; TH2F *efficiency_first_shell = (TH2F*) factor_events_rec-> Clone()

Definition at line 1290 of file cbc_plots.C.

int countv =0

Definition at line 587 of file cbc_plots.C.

int cz

Definition at line 434 of file cbc_plots.C.

TH2F* D_Mchirp_inj = new TH2F("Distance vs MChirp inj.","",1000,MINCHIRP,MAXCHIRP,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000./FACTORS[0])

Definition at line 296 of file cbc_plots.C.

TH2F* D_Mchirp_rec = new TH2F("Distance vs MChirp rec.","",1000,MINCHIRP,MAXCHIRP,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000/FACTORS[0])

Definition at line 317 of file cbc_plots.C.

TH2F* D_Mtot_inj = new TH2F("Distance vs Mtot inj.","",1000,MINMtot,MAXMtot,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000./FACTORS[0])

Definition at line 269 of file cbc_plots.C.

TH2F* D_Mtot_rec = new TH2F("Distance vs Mtot rec.","",1000,MINMtot,MAXMtot,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000/FACTORS[0])

Definition at line 290 of file cbc_plots.C.

TH2F* D_q_inj = new TH2F("Distance vs q inj.","",1000,MINRATIO,MAXRATIO,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000./FACTORS[0])

Definition at line 372 of file cbc_plots.C.

TH2F* D_q_rec = new TH2F("Distance vs q rec.","",1000,MINRATIO,MAXRATIO,5000,MINDISTANCE/1000./FACTORS[nfactor-1],MAXDISTANCE/1000/FACTORS[0])

Definition at line 393 of file cbc_plots.C.

bool DDistrUniform

Definition at line 114 of file cbc_plots.C.

bool DDistrVolume

Definition at line 114 of file cbc_plots.C.

float distance

Definition at line 436 of file cbc_plots.C.

TH1F* Dt = new TH1F("Dt","",1000,-0.5,0.5)

Definition at line 398 of file cbc_plots.C.

else {float MINMtot = 0.0

Definition at line 166 of file cbc_plots.C.

error_volume[m][l] = pow(dV,2)

Definition at line 679 of file cbc_plots.C.

TExec* ex2 = new TExec("ex2","gStyle->SetPaintTextFormat(\".0f\");")

Definition at line 1359 of file cbc_plots.C.

TF1* f2 = cbcTool.doRangePlot(RHO_NBINS, Trho, Rrho, eRrho, RHO_MIN, T_out, c1, networkname, netdir, write_ascii)

Definition at line 1710 of file cbc_plots.C.

float factor

Definition at line 436 of file cbc_plots.C.

TH2F* factor_events_inj[nfactor]

Definition at line 262 of file cbc_plots.C.

TH2F* factor_events_rec = new TH2F("factor_events_rec","",NBINS_mass,MIN_MASS,MAX_MASS,NBINS_mass2,min_mass2,max_mass2)

Definition at line 266 of file cbc_plots.C.

factor_events_rec Fill

(

time-

time[0][nIFO]

)

Definition at line 586 of file cbc_plots.C.

char fname[1024]

Definition at line 720 of file cbc_plots.C.

cout<< "nwave : " << nwave << endl; sprintf(title,"%s",newcut); sprintf(newcut2,"%s && %s",newcut,veto_not_vetoed); cout << "final cut " << newcut2 << endl; TString sel_fin = sel; sel_fin.ReplaceAll("hist2","hist3"); sim.Draw(sel_fin,newcut2,"same"); TH2F *htemp3 = (TH2F*)gPad-> GetPrimitive("hist3")

Definition at line 794 of file cbc_plots.C.

double green[NRGBs] = { 0.97, 0.30, 0.40, 0.97, 0.00, 0.00 }

Definition at line 25 of file cbc_plots.C.

TH2F* h_radius = new TH2F("h_radius","",NBINS_mass,MIN_MASS,MAX_MASS,NBINS_mass2,min_mass2,max_mass2)

Definition at line 1322 of file cbc_plots.C.

TH2F* htemp4 = (TH2F*)gPad->GetPrimitive("hist4")

Definition at line 874 of file cbc_plots.C.

TH2F* htemp5 = (TH2F*)gPad->GetPrimitive("hist5")

Definition at line 910 of file cbc_plots.C.

TH2F* inj_events = new TH2F("inj_events","D_Minj",NBINS_mass,MIN_MASS,MAX_MASS,NBINS_mass2,min_mass2,max_mass2)

Definition at line 209 of file cbc_plots.C.

char inj_title[256]

Definition at line 704 of file cbc_plots.C.

TString Insp = MDC.GetInspiral()

Definition at line 111 of file cbc_plots.C.

l = TMath::FloorNint((m2-min_mass2)/MASS_BIN)

Definition at line 434 of file cbc_plots.C.

char lab[256]

Definition at line 844 of file cbc_plots.C.

leg_snr = new TLegend(0.6,0.755,0.885,0.923,"","brNDC")

Definition at line 845 of file cbc_plots.C.

m = TMath::FloorNint((m1-MIN_MASS)/MASS_BIN)

Definition at line 434 of file cbc_plots.C.

m1 =mass[1]

Definition at line 606 of file cbc_plots.C.

m2 =mass[0]

Definition at line 607 of file cbc_plots.C.

float mass[2]

Definition at line 437 of file cbc_plots.C.

int massbins = 0

Definition at line 1299 of file cbc_plots.C.

int MAX_AXIS_Z = inj_events->GetBinContent(inj_events->GetMaximumBin()) + 1

Definition at line 700 of file cbc_plots.C.

float max_mass1 = MAX_MASS

Definition at line 51 of file cbc_plots.C.

float max_mass2 = MAX_MASS

Definition at line 39 of file cbc_plots.C.

float MAXCHIRP = 1.1*pow(max_mass1*max_mass2,3./5.)/pow(minMtot,1./5.)

Definition at line 187 of file cbc_plots.C.

float maxDistance

Definition at line 113 of file cbc_plots.C.

float maxMtot

Definition at line 113 of file cbc_plots.C.

float maxRatio

Definition at line 113 of file cbc_plots.C.

float MAXRATIO = 1.

Definition at line 151 of file cbc_plots.C.

float mchirp

Definition at line 436 of file cbc_plots.C.

float min_mass1 = MIN_MASS

Definition at line 45 of file cbc_plots.C.

float min_mass2 = MIN_MASS

Definition at line 57 of file cbc_plots.C.

float MINCHIRP = 0.9*pow(min_mass1*min_mass2,3./5.)/pow(maxMtot,1./5.)

Definition at line 186 of file cbc_plots.C.

float minDistance

Definition at line 113 of file cbc_plots.C.

float minMtot

Definition at line 113 of file cbc_plots.C.

float minRatio

Definition at line 113 of file cbc_plots.C.

float MINRATIO = 0.0

Definition at line 143 of file cbc_plots.C.

int mt

Definition at line 434 of file cbc_plots.C.

int NBINS_mass = (int)((MAX_MASS-MIN_MASS)/MASS_BIN)

Definition at line 33 of file cbc_plots.C.

int NBINS_mass1 = (int)((max_mass1-min_mass1)/MASS_BIN)

Definition at line 62 of file cbc_plots.C.

int NBINS_mass2 = (int)((max_mass2-min_mass2)/MASS_BIN)

Definition at line 63 of file cbc_plots.C.

network* net =NULL

Definition at line 105 of file cbc_plots.C.

char networkname[256]

Definition at line 95 of file cbc_plots.C.

int nmdc = mdc.GetSelectedRows()

Definition at line 791 of file cbc_plots.C.

nT =TMath::Min(TMath::Floor((rho[pp_irho]-RHO_MIN)/RHO_BIN),(double)RHO_NBINS)+1

Definition at line 659 of file cbc_plots.C.

int nwave = sim.GetSelectedRows()

Definition at line 820 of file cbc_plots.C.

int nwave_final = sim.GetSelectedRows()

Definition at line 838 of file cbc_plots.C.

TPaveText* p_inj = new TPaveText(0.325301,0.926166,0.767068,0.997409,"blNDC")

Definition at line 711 of file cbc_plots.C.

TPaveText* p_radius = new TPaveText(0.325301,0.926166,0.767068,0.997409,"blNDC")

Definition at line 1368 of file cbc_plots.C.

char radius_title[256]

Definition at line 1363 of file cbc_plots.C.

TH2F* rec_events = new TH2F("rec_events","D_Mrec",NBINS_mass,MIN_MASS,MAX_MASS,NBINS_mass2,min_mass2,max_mass2)

Definition at line 238 of file cbc_plots.C.

double red[NRGBs] = { 0.00, 0.00, 0.00, 0.97, 0.97, 0.10 }

Definition at line 24 of file cbc_plots.C.

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

Definition at line 414 of file cbc_plots.C.

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

Definition at line 404 of file cbc_plots.C.

float spin[6]

Definition at line 438 of file cbc_plots.C.

double stops[NRGBs] = { 0.10, 0.25, 0.45, 0.60, 0.75, 1.00 }

Definition at line 23 of file cbc_plots.C.

TText * text = p_inj->AddText(inj_title)

Definition at line 717 of file cbc_plots.C.

double time[6]

Definition at line 435 of file cbc_plots.C.

double V0 = 0.0

Definition at line 1300 of file cbc_plots.C.

volume[m][l] = dV

Definition at line 678 of file cbc_plots.C.

TString waveform

Definition at line 115 of file cbc_plots.C.

bool write_ascii = false

Definition at line 447 of file cbc_plots.C.