cbc_plots_sim4_ifar.C File Reference

#include "cwb.hh"
#include "config.hh"
#include "network.hh"
#include "wavearray.hh"
#include "TString.h"
#include "TObjArray.h"
#include "TObjString.h"
#include "TRandom.h"
#include "TComplex.h"
#include "TMath.h"
#include "mdc.hh"
#include <vector>
#include "GToolbox.hh"

Include dependency graph for cbc_plots_sim4_ifar.C:

Macros
#define	CONTOURS   7
#define	LIV_FILE_NAME   "live.txt"
#define	NCont   99
#define	NRGBs   6
#define	RHO_BIN   0.1
#define	RHO_MIN   5.0
#define	RHO_NBINS   5000
#define	XIFO   4

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")
leg_D	AddEntry (D_Mtot_rec, lab,"p")
leg_D	AddEntry (D_Mtot_inj, lab,"p")
leg_D2	AddEntry (t0,"IFAR<1 yr","p")
leg_D2	AddEntry (t1,"1 yr < IFAR < 10 yr","p")
leg_D2	AddEntry (t10,"10 yr < IFAR < 100 yr","p")
leg_D2	AddEntry (t100,"IFAR > 100 yr","p")
cout<< "Length of ifar/volume vector: "<< vV.size()<< endl;for(int i=0;i< vV.size();i++){veV[i]=TMath::Sqrt(veV[i]);vfar[i]*=TRIALS *CYS;vefar[i]*=TRIALS *CYS;}c1->	Clear ()
efficiency	Divide (inj_events)
gr	Draw ("ap")
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 ()
hcandle	Draw ("CANDLE")
mdc	Draw (sel,"")
sim	Draw (sel, newcut,"")
sim	Draw (sel, newcut2,"same")
sim	Draw (sel, newcut2)
D_Mtot_inj	Draw ("p")
D_Mtot_rec	Draw ("p same")
D_inj	Draw ("lp")
D_rec	Draw ("lp same")
efficiency	Draw ("colz text")
	exit (1)
gSystem	Exit (0)
	ROC Curve. More...
D_Mtot_inj	Fill (mass[1]+mass[0], distance)
D_Mchirp_inj	Fill (mchirp, distance)
D_q_inj	Fill (mass[0]/mass[1], distance)
factor_events_inj[ifactor]	Fill (mass[0], mass[1])
	for (int i=0;i< nfactor;i++)
	for (int i=1;i<=NBINS_mass;i++)
gROOT	ForceStyle ()
sim	GetHistogram () -> Fit("gaus")
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 (ShellminDistance > minDistanceXML[gIFACTOR-1])
	if (MDC.GetInspiralOption("--max-distance")!="")
	if (ShellmaxDistance< maxDistanceXML[gIFACTOR-1])
	if (MDC.GetInspiralOption("--min-mratio")!="")
	if (MDC.GetInspiralOption("--max-mratio")!="")
	if ((MDC.GetInspiralOption("--min-mass1")!="")&&(MDC.GetInspiralOption("--max-mass2")!=""))
	if ((MDC.GetInspiralOption("--max-mass1")!="")&&(MDC.GetInspiralOption("--min-mass2")!=""))
	if (MDC.GetInspiralOption("--d-distr")!="")
	if (MDC.GetInspiralOption("--dchirp-distr")!="")
	if (sim.GetListOfBranches() ->FindObject("ifar"))
	if (minchi)
	if ((ifactor > nfactor-1)||(ifactor< 0))
	if (write_ascii)
	if (Redshift)
	if (if() c1->SetLogy(kTRUE) pp_rho_log)
cout<< "Number of Factors:" << nfactor<< endl;for(int l=0;l < nfactor;l++){gIFACTOR=l+1;FACTORS[l]=gIFACTOR;gROOT->	Macro (configPlugin.GetTitle())
TChain	mdc ("mdc")
finj	open (fname3, std::ofstream::out)
vV	push_back (vdv[mindex[0]])
veV	push_back (pow(vdv[mindex[0]], 2))
vsifar	push_back (vifar[mindex[0]])
vfar	push_back (1./vifar[mindex[0]])
vefar	push_back (TMath::Sqrt(TMath::Nint(liveTot/vifar[mindex[0]]))/liveTot)
c1	SaveAs (fname)
c1	SetBorderMode (0)
c1	SetBorderSize (2)
p_inj	SetBorderSize (0)
c1	SetBottomMargin (0.1450777)
mdc	SetBranchAddress ("time", mytime)
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)
D_inj	SetFillColor (2)
D_rec	SetFillColor (4)
htemp2	SetFillStyle (3017)
gStyle	SetFrameBorderMode (0)
c1	SetGridx ()
c1	SetGridy ()
gr	SetLineColor (kBlue)
htemp	SetLineColor (kBlack)
htemp2	SetLineColor (kRed)
D_inj	SetLineColor (2)
D_rec	SetLineColor (4)
htemp	SetLineWidth (3)
htemp2	SetLineWidth (2)
D_inj	SetLineWidth (4)
c1	SetLogx (1)
c1	SetLogx (0)
c1	SetLogx (false)
c1	SetLogy (1)
c1	SetLogy (0)
D_Mtot_inj	SetMarkerColor (2)
D_Mtot_rec	SetMarkerColor (4)
rhocc	SetMarkerColor (1)
sim	SetMarkerColor (kRed)
sim	SetMarkerColor (kBlue)
t0	SetMarkerColor (kMagenta)
t1	SetMarkerColor (kCyan)
t100	SetMarkerColor (kYellow)
D_Mtot_inj	SetMarkerSize (0.5)
gr	SetMarkerSize (1.0)
D_Mtot_inj	SetMarkerStyle (20)
D_Mtot_inj	SetName ("D_Mtotinj")
D_Mchirp_inj	SetName ("D_Chirp_inj")
D_Chi_inj	SetName ("D_Chi_inj")
D_Chi_injx	SetName ("D_Chi_injx")
D_Chi_injy	SetName ("D_Chi_injy")
D_q_inj	SetName ("D_q_inj")
efficiency	SetName ("efficiency")
gStyle	SetNumberContours (256)
gStyle	SetOptFit (1)
gStyle	SetOptFit (0)
gStyle	SetOptStat (kFALSE)
cout<< "Selection: "<< sel << endl;gStyle->	SetOptStat (1)
gStyle	SetOptStat (0)
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", myptitle.Data());htemp->	SetTitle (title)
D_inj	SetTitle (0)
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 (fname3,"%s/injected_signals.txt", netdir)
	sprintf (line,"#GPS@L1 mass1 mass2 distance spinx1 spiny1 spinz1 ""spinx2 spiny2 spinz2 \n")
	sprintf (fname,"%s/ROCV.png", netdir)
	sprintf (fname,"%s/ROC.png", netdir)
	sprintf (inj_title,"Injected events")
	sprintf (fname,"%s/Injected_mass1_mass2.eps", netdir)
	sprintf (fname,"%s/Dchirp_candle.png", netdir)
	sprintf (fname,"%s/Dchirp_candle2.png", 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_cut)
	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 (fname,"%s/Relative_snr_Loss.png", netdir)
	sprintf (lab,"Injections: %i",(int) mdc.GetEntries())
	sprintf (lab,"found: %i", cnt)
	sprintf (lab,"missed: %i",(int) mdc.GetEntries()-cnt)
	sprintf (fname,"%s/Distance_vs_total_mass.eps", netdir)
	sprintf (fname,"%s/Distance_vs_total_mass_ifar.eps", netdir)
	sprintf (fname,"%s/Distance_vs_chirp_mass.eps", netdir)
	sprintf (fname,"%s/Distance_vs_mass_ratio.eps", netdir)
	sprintf (fname,"%s/Injected_distances_distribution.png", netdir)
	sprintf (fname,"%s/Total_mass_distribution.png", netdir)
	sprintf (fname,"%s/Chirp_mass_distribution.png", netdir)
	sprintf (fname,"%s/Mass_ratio_distribution.png", netdir)
	sprintf (eff_title,"Efficiency")
	sprintf (radius_title,"%s : Effective radius (Mpc)", networkname)
	sprintf (fname,"%s/Effective_radius.png", netdir)
	strcpy (cfg->tmp_dir,"tmp")
c1	Update ()

Variables
int	actual_nfactor = 0
cout<< "Injected signals: " << mdc.GetEntries()<< endl;cout << "Injected signals in histogram factor_events_inj: " << NEVTS<< endl;float myifar, ecor, m1, m2, netcc[3], neted, penalty;float rho[2];float chirp[6];float range[2];float frequency[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("ifar",&myifar);sim.SetBranchAddress("ecor",&ecor);sim.SetBranchAddress("penalty",&penalty);sim.SetBranchAddress("time", mytime);sim.SetBranchAddress("iSNR", iSNR);sim.SetBranchAddress("sSNR", sSNR);sim.SetBranchAddress("spin", spin);sim.SetBranchAddress("frequency", frequency);float **volume=new float *[NBINS_mass1];float **volume_first_shell=new float *[NBINS_mass1];float **radius=new float *[NBINS_mass1];float **error_volume=new float *[NBINS_mass1];float **error_volume_first_shell=new float *[NBINS_mass1];float **error_radius=new float *[NBINS_mass1];for(int i=0;i < NBINS_mass1;i++){volume[i]=new float[NBINS_mass2];volume_first_shell[i]=new float[NBINS_mass2];radius[i]=new float[NBINS_mass2];error_volume[i]=new float[NBINS_mass2];error_volume_first_shell[i]=new float[NBINS_mass2];error_radius[i]=new float[NBINS_mass2];for(int j=0;j < NBINS_mass2;j++){volume[i][j]=0.;volume_first_shell[i][j]=0.;radius[i][j]=0.;error_volume[i][j]=0.;error_volume_first_shell[i][j]=0.;error_radius[i][j]=0.;}}float **spin_mtot_volume=new float *[NBINS_MTOT+1];float **spin_mtot_radius=new float *[NBINS_MTOT+1];float **error_spin_mtot_volume=new float *[NBINS_MTOT+1];float **error_spin_mtot_radius=new float *[NBINS_MTOT+1];for(int i=0;i< NBINS_MTOT+1;i++){spin_mtot_volume[i]=new float[NBINS_SPIN+1];spin_mtot_radius[i]=new float[NBINS_SPIN+1];error_spin_mtot_volume[i]=new float[NBINS_SPIN+1];error_spin_mtot_radius[i]=new float[NBINS_SPIN+1];for(int j=0;j < NBINS_SPIN+1;j++){spin_mtot_volume[i][j]=0.;error_spin_mtot_volume[i][j]=0.;spin_mtot_radius[i][j]=0.;error_spin_mtot_radius[i][j]=0.;}}char fname[1024];sprintf(fname,"%s/recovered_signals.txt", netdir);ofstream fev;fev.open(fname, std::ofstream::out);sprintf(line,"#GPS@L1 FAR[Hz] eFAR[Hz] Pval ""ePval factor rho frequency iSNR sSNR \n");fev<< line<< endl;ofstream *fev_single=new ofstream[nfactor];for(int l=1;l< nfactor+1;l++){sprintf(fname,"%s/recovered_signals_%d.txt", netdir, l);fev_single[l-1].open(fname, std::ofstream::out);fev_single[l-1] << line<< endl;}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, nevts, internal_volume;int nT;int countv=0;int cnt=0;int cnt2=0;int cntfreq=0;bool bcut=false;double liveTot=sim.GetMaximum("ifar");double BKG_LIVETIME_yr=liveTot/CYS;double BKG_LIVETIME_Myr=BKG_LIVETIME_yr/(1.e6);cout.precision(14);cout << "Total live time --->	background
double	blue [NRGBs] = {0.97, 0.97, 0.00, 0.00, 0.00, 0.00}
bool	bmaxDistance
bool	bmaxMass1
bool	bmaxMass2
bool	bmaxMtot
bool	bmaxRatio
bool	bminDistance
bool	bminMass1
bool	bminMass2
bool	bminMtot
bool	bminRatio
CWB::config *	cfg = new CWB::config
float	chi [3]
float	CYS
int	cz
TH2F *	D_Chi_injx = (TH2F *)D_Chi_inj->Clone("D_Chi_injx")
TH2F *	D_Chi_injy = (TH2F *)D_Chi_inj->Clone("D_Chi_injy")
TH2F *	D_Chi_rec
TH2F *	D_Chi_recx = (TH2F *)D_Chi_rec->Clone("D_Chi_recx")
TH2F *	D_Chi_recy = (TH2F *)D_Chi_rec->Clone("D_Chi_recy")
TH3F *	D_dchirp_rec
TH1D *	D_inj = D_Mtot_inj->ProjectionY()
TH2F *	D_Mchirp_inj
TH2F *	D_Mchirp_rec
TH2F *	D_Mtot_inj
TH2F *	D_Mtot_rec
TH3F *	D_Mtot_rec3
TH2F *	D_q_inj
TH2F *	D_q_rec
TH1D *	D_rec = D_Mtot_rec->ProjectionY()
TH2F *	dchirp_rec
bool	DDistrChirpMass
bool	DDistrUniform
bool	DDistrVolume
float	distance
TH1F *	Dt = new TH1F("Dt", "", 1000, -0.5, 0.5)
char	eff_title [256]
TH2F *	efficiency = (TH2F *)rec_events->Clone()
TH2F *	efficiency_first_shell = (TH2F *)factor_events_rec->Clone()
	else
delete[]	error_radius
delete[]	error_spin_mtot_radius
delete[]	error_spin_mtot_volume
delete[]	error_volume
delete[]	error_volume_first_shell
TExec *	ex1 = new TExec("ex1", "gStyle->SetPaintTextFormat(\".2f\");")
TExec *	ex2 = new TExec("ex2", "gStyle->SetPaintTextFormat(\".0f\");")
TF1 *	f2
float	factor
TH2F *	factor_events_inj = new TH2F[nfactor]
TH2F *	factor_events_rec
delete[]	factor_events_spin_mtot_inj
float *	FACTORS = new float[nfactor]
finj<< line<< endl;ofstream *finj_single=new ofstream[nfactor];TString xml;for(int l=0;l< nfactor;l++){sprintf(fname3,"%s/injected_signals_%d.txt", netdir, l+1);finj_single[l].open(fname3, std::ofstream::out);finj_single[l] << line<< endl;}for(int g=0;g < (int) mdc.GetEntries();g++){mdc.GetEntry(g);if(Redshift){mass[0]=SFmasses[factor-1][0];mass[1]=SFmasses[factor-1][1];mchirp=pow(mass[0]*mass[1], 3./5.)/pow(mass[0]+mass[1], 1./5.);}inj_events->	Fill (mass[0], mass[1])
ofstream	finj
bool	FixedFiducialVolume
char	fname3 [2048]
cout<< "nmdc : "<< nmdc<< endl;TH2F *htemp=(TH2F *) gPad->	GetPrimitive ("hist")
int	gIFACTOR =-1
TGraphErrors *	gr
TGraphErrors *	gr2
double	green [NRGBs] = {0.97, 0.30, 0.40, 0.97, 0.00, 0.00}
TH2F *	h_radius
TH1 *	hcandle = D_dchirp_rec->Project3D("yx")
TH2F *	htemp4 = (TH2F *)gPad->GetPrimitive("hist4")
TH2F *	htemp5 = (TH2F *)gPad->GetPrimitive("hist5")
int	ifactor = (int)factor - 1
TH2F *	inj_events
char	inj_title [256]
TString	Insp = MDC.GetInspiral()
int	l
char	lab [256]
	leg_D
	leg_D2
	leg_snr = new TLegend(0.6, 0.755, 0.885, 0.923, "", "brNDC")
char	line [1024]
int	m
float	M1
float	M2
float	mass [2]
int	massbins = 0
int	MAX_AXIS_Z
float	max_mass1 = MAX_MASS
float	max_mass2 = MAX_MASS
float *	maxDistance = new float[nfactor]
float *	maxDistanceXML = new float[nfactor]
double	MAXIFAR = 0.0
float *	maxMChirp = new float[nfactor]
float *	maxMtot = new float[nfactor]
float *	maxRatio = new float[nfactor]
float	mchirp
TH1D *	Mchirp_inj = D_Mchirp_inj->ProjectionX()
TH1D *	Mchirp_rec = D_Mchirp_rec->ProjectionX()
int	mcount = 0
CWB::mdc	MDC (net)
float	min_mass1 = MIN_MASS
float	min_mass2 = MIN_MASS
bool	minchi
float *	minDistance = new float[nfactor]
float *	minDistanceXML = new float[nfactor]
float *	minMChirp = new float[nfactor]
float *	minMtot = new float[nfactor]
float *	minRatio = new float[nfactor]
TH1D *	Mtot_inj = D_Mtot_inj->ProjectionX()
TH1D *	Mtot_rec = D_Mtot_rec->ProjectionX()
int	mtt
double	mytime [6]
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)
int	NBINS_SPIN = (int)((MAXCHI - MINCHI) / CHI_BIN)
network *	net = NULL
char	networkname [256]
double	NEVTS = 0.0
Int_t	nGraphs = 0
int	nmdc = mdc.GetSelectedRows()
int	nwave = sim.GetSelectedRows()
int	nwave_final = sim.GetSelectedRows()
TPaveText *	p_eff
TPaveText *	p_inj
TPaveText *	p_radius
int	point
TH1D *	q_inj = D_q_inj->ProjectionX()
TH1D *	q_rec = D_q_rec->ProjectionX()
delete[]	radius
char	radius_title [256]
cout<< "MDC set: "<< gIFACTOR << endl;cout<< "xml conf: waveform="<< waveforms[gIFACTOR-1] << " minMtot="<< minMtot[gIFACTOR-1] << " maxMtot="<< maxMtot[gIFACTOR-1] << " minDistance=" << minDistance[gIFACTOR-1] << " maxDistance=" << maxDistance[gIFACTOR-1] << " minRatio="<< minRatio[gIFACTOR-1] << " maxRatio="<< maxRatio[gIFACTOR-1] << endl;}float MINMtot=0.0;if(bminMtot){float MINMtot=0.99 *minMtot[gIFACTOR-1];}float MAXMtot=0.0;int NBINS_MTOT=0;if(bmaxMtot){MAXMtot=MAX_MASS;NBINS_MTOT=TMath::FloorNint((MAX_MASS-MIN_MASS)/MASS_BIN/2.);cout << "NBINS_MTOT: "<< NBINS_MTOT << endl;}else{cout << "Undefined maximal total mass!! Define float MAXMtot" << endl;exit(1);}float MINDISTANCE=0.0;if(bminDistance){MINDISTANCE=0.9 *ShellminDistance;cout << "MINDISTANCE = " << MINDISTANCE<< endl;}else{cout << "Undefined minimal distance!! Defined float MINDISTANCE"<< endl;}float MAXDISTANCE=0.0;if(bmaxDistance){MAXDISTANCE=TMath::Max(ShellmaxDistance, maxDistance[gIFACTOR-1]);cout << "MAXDISTANCE = " << MAXDISTANCE<< endl;}else{cout << "Undefined maximal distance!! Define float MAXDISTANCE"<< endl;exit(1);}float MINCHIRP=100.0;float MAXCHIRP=0.0;float MINRATIO=1.0;float MAXRATIO=1.0;if((bminRatio)&&(bmaxRatio)){for(int l=0;l< nfactor;l++){if(MINRATIO > minRatio[l]){MINRATIO=minRatio[l];}if(MAXRATIO < maxRatio[l]){MAXRATIO=maxRatio[l];}if(MINCHIRP > minMtot[l]*pow(maxRatio[l], 3./5.)/pow(1+minRatio[l], 6./5.)){MINCHIRP=minMtot[l]*pow(maxRatio[l], 3./5.)/pow(1+minRatio[l], 6./5.);};if(MAXCHIRP< maxMtot[l]*pow(minRatio[l], 3./5.)/pow(1+minRatio[l], 6./5.)){MAXCHIRP=maxMtot[l]*pow(minRatio[l], 3./5.)/pow(1+minRatio[l], 6./5.);};}}else{cout << "Undefined minRatio.. " << endl;exit(1);}for(int l=0;l < nfactor-1;l++){if((minDistanceXML[l]==minDistanceXML[l+1])&&(maxDistanceXML[l]==maxDistanceXML[l+1])){FixedFiducialVolume=1;}else{FixedFiducialVolume=1;cout << "Beware: different fiducial volumes for ""different factors!!"<< endl;}}if(strlen(veto_not_vetoed)==0){sprintf(veto_not_vetoed,"%s", ch2);}cout<< "ERROR! Missing zero lag livetime...Please add:"<< endl;cout<< "\"#define LIVE_ZERO XXXXXXXXXX\" where XXXXXXXXXX is the ""livetime in seconds to your config/user_pparameters.C file...."<< endl;exit(1);TCanvas *c1=new TCanvas("c1","c1", 3, 47, 1000, 802);c1->	Range (-1.216392,-477.6306, 508.8988, 2814.609)
TH2F *	rec_events
double	red [NRGBs] = {0.00, 0.00, 0.00, 0.97, 0.97, 0.10}
bool	Redshift
TH2F *	rho_pf
TH2F *	rhocc
float *	shell_volume = new float[nfactor]
float	ShellmaxDistance = 0.0
float	ShellminDistance = 9999999.0
float	spin [6]
delete[]	spin_mtot_radius
delete[]	spin_mtot_volume
double	stops [NRGBs] = {0.10, 0.25, 0.45, 0.60, 0.75, 1.00}
TH1 *	t0 = D_Mtot_rec3->Project3D("yx_0")
TH1 *	t1 = D_Mtot_rec3->Project3D("yx_1")
TH1 *	t10 = D_Mtot_rec3->Project3D("yx_10")
TH1 *	t100 = D_Mtot_rec3->Project3D("yx_100")
TText *	text = p_inj->AddText(inj_title)
TLatex	Tl
TLatex	Tl2
Int_t	TotalConts = 0
float	Tscale = 1.0
double	V0 = 0.0
char	val [20]
std::vector< double >	vdv
std::vector< double >	vefar
std::vector< double >	veR
std::vector< double >	veV
std::vector< double >	vfar
std::vector< double >	vifar
delete[]	volume
delete[]	volume_first_shell
std::vector< double >	vR
std::vector< double >	vsifar
cout<< endl;if(write_ascii){fev.close();for(int l=0;l< nfactor;l++){fev_single[l].close();}}cout << "Recovered entries: "<< cnt << endl;cout<< "Recovered entries: "<< cnt2<< endl;cout << "Recovered entries cut by frequency: "<< cntfreq<< endl;cout << "Recovered entries vetoed: " << countv<< endl;cout<< "dV : " << dV<< " dV1 : "<< dV1<< endl;internal_volume=4./3.*TMath::Pi()*pow(minDistance[0]/1000., 3.);if(INCLUDE_INTERNAL_VOLUME){for(int i=0;i< vdv.size();i++){if(vdv[i] > 0.0){vdv[i]+=internal_volume;}}for(int i=0;i< RHO_NBINS;i++){if(Vrho[i] > 0.0){Vrho[i]+=internal_volume;}}for(int i=0;i< NBINS_MTOT+1;i++){for(int j=0;j< NBINS_SPIN+1;j++){if(spin_mtot_volume[i][j] > 0.0){spin_mtot_volume[i][j]+=internal_volume;}}}for(int i=0;i< NBINS_mass;i++){for(int j=0;j< NBINS_mass2;j++){if(volume[i][j] > 0.0){volume[i][j]+=internal_volume;}}}}Int_t *mindex=new Int_t[vdv.size()];TMath::Sort((int) vdv.size(),&vifar[0], mindex, true);std::vector< double >	vV
TString *	waveforms = new TString[nfactor]
bool	write_ascii = false
Double_t	x0
Float_t	xq [6] = {8.0, 15.0, 25.0, 35.0, 45.0, 55.0}
Double_t	Y0
Float_t *	yq = new Float_t[101]
Double_t	z0

Macro Definition Documentation

#define CONTOURS   7

Definition at line 5 of file cbc_plots_sim4_ifar.C.

#define LIV_FILE_NAME   "live.txt"

Definition at line 7 of file cbc_plots_sim4_ifar.C.

#define NCont   99

#define NRGBs   6

#define RHO_BIN   0.1

Definition at line 3 of file cbc_plots_sim4_ifar.C.

#define RHO_MIN   5.0

Definition at line 2 of file cbc_plots_sim4_ifar.C.

#define RHO_NBINS   5000

Definition at line 4 of file cbc_plots_sim4_ifar.C.

#define XIFO   4

Definition at line 9 of file cbc_plots_sim4_ifar.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"
	)

leg_D AddEntry	(	D_Mtot_rec	,
		lab	,
		"p"
	)

leg_D AddEntry	(	D_Mtot_inj	,
		lab	,
		"p"
	)

leg_D2 AddEntry

(

t0

)

leg_D2 AddEntry

(

t1

)

leg_D2 AddEntry

(

t10

)

leg_D2 AddEntry	(	t100	,
		"	IFAR,
		100 yr"	,
		"p"
	)

cout<< "Length of ifar/volume vector: " << vV.size() << endl; for (int i = 0; i < vV.size(); i++) { veV[i] = TMath::Sqrt(veV[i]); vfar[i] *= TRIALS * CYS; vefar[i] *= TRIALS * CYS; } c1-> Clear

(

)

Definition at line 1261 of file CWB_Plugin_xWRC.C.

efficiency Divide

(

inj_events

)

gr2 Draw

(

"ap"

)

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

(

)

hcandle Draw

(

"CANDLE"

)

mdc Draw	(	sel	,
		""
	)

sim Draw	(	sel	,
		newcut	,
		""
	)

sim Draw	(	sel	,
		newcut2	,
		"same"
	)

sim Draw	(	sel	,
		newcut2
	)

D_q_inj Draw

(

"p"

)

D_q_rec Draw

(

"p same"

)

q_inj Draw

(

"lp"

)

q_rec Draw

(

"lp same"

)

efficiency Draw

(

"colz text"

)

exit

(

1

)

gSystem Exit

(

0

)

ROC Curve.

D_Mtot_inj Fill	(	mass+	mass[1][0],
		distance
	)

D_Mchirp_inj Fill	(	mchirp	,
		distance
	)

D_q_inj Fill	(	mass/	mass[0][1],
		distance
	)

factor_events_inj [ifactor] Fill	(	mass	[0],
		mass	[1]
	)

for

(

)

Definition at line 539 of file cbc_plots_sim4_ifar.C.

for

(

int

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

)

Definition at line 2311 of file cbc_plots_sim4_ifar.C.

gROOT ForceStyle

(

)

sim GetHistogram

(

)

-> Fit("gaus")

Definition at line 913 of file cbc_plots.C.

inj_events GetXaxis

(

)

-> SetRangeUser(MIN_plot_mass1, MAX_plot_mass1)

Definition at line 488 of file cbc_plots_sim4_ifar.C.

inj_events GetYaxis

(

)

-> SetRangeUser(MIN_plot_mass2, MAX_plot_mass2)

Definition at line 489 of file cbc_plots_sim4_ifar.C.

inj_events GetZaxis

(

)

-> SetTickLength(0.01)

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

Definition at line 114 of file cbc_plots_sim4_ifar.C.

if

(

MDC.GetInspiralOption("--waveform")!

= ""

)

Definition at line 168 of file cbc_plots_sim4_ifar.C.

if

(

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

= ""

)

Definition at line 172 of file cbc_plots_sim4_ifar.C.

if

(

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

= ""

)

Definition at line 177 of file cbc_plots_sim4_ifar.C.

if

(

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

= ""

)

Definition at line 183 of file cbc_plots_sim4_ifar.C.

if	(	ShellminDistance	,
		minDistanceXML	[gIFACTOR-1]
	)

Definition at line 189 of file cbc_plots_sim4_ifar.C.

if

(

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

= ""

)

Definition at line 192 of file cbc_plots_sim4_ifar.C.

if

(

)

Definition at line 198 of file cbc_plots_sim4_ifar.C.

if

(

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

= ""

)

Definition at line 202 of file cbc_plots_sim4_ifar.C.

if

(

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

= ""

)

Definition at line 207 of file cbc_plots_sim4_ifar.C.

if

(

(MDC.GetInspiralOption("--min-mass1")!="")&&(MDC.GetInspiralOption("--max-mass2")!="")

)

Definition at line 212 of file cbc_plots_sim4_ifar.C.

if

(

(MDC.GetInspiralOption("--max-mass1")!="")&&(MDC.GetInspiralOption("--min-mass2")!="")

)

Definition at line 219 of file cbc_plots_sim4_ifar.C.

if

(

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

= ""

)

Definition at line 226 of file cbc_plots_sim4_ifar.C.

if

(

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

= ""

)

Definition at line 260 of file cbc_plots_sim4_ifar.C.

if

(

sim.

GetListOfBranches) ->FindObject("ifar"

)

Definition at line 785 of file cbc_plots_sim4_ifar.C.

if

(

minchi

)

Definition at line 867 of file cbc_plots_sim4_ifar.C.

if

(

(ifactor > nfactor-1)||(ifactor< 0)

)

Definition at line 882 of file cbc_plots_sim4_ifar.C.

if

(

write_ascii

)

Definition at line 895 of file cbc_plots_sim4_ifar.C.

if

(

Redshift

)

Definition at line 1452 of file cbc_plots_sim4_ifar.C.

if

(

if () c1->SetLogy(kTRUE)

pp_rho_log

)

Definition at line 1616 of file cbc_plots_sim4_ifar.C.

cout<< "Number of Factors:" << nfactor << endl; for (int l = 0; l < nfactor; l++) { gIFACTOR = l + 1; FACTORS[l] = gIFACTOR; gROOT-> Macro

(

configPlugin.

GetTitle()

)

TChain mdc

(

"mdc"

)

finj open	(	fname3	,
		std::ofstream::out
	)

vV push_back

(

vdv]

[mindex[0]

)

veV push_back

(

pow(vdv[mindex[0]], 2)

)

vsifar push_back

(

vifar]

[mindex[0]

)

vfar push_back

(

1./vifar]

[mindex[0]

)

vefar push_back

(

TMath::Sqrt(TMath::Nint(liveTot/vifar[mindex[0]]))/

liveTot

)

c1 SaveAs

(

fname

)

c1 SetBorderMode

(

0

)

c1 SetBorderSize

(

2

)

p_inj SetBorderSize

(

0

)

c1 SetBottomMargin

(

0.

1450777

)

mdc SetBranchAddress	(	"time"	,
		mytime
	)

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

)

D_inj SetFillColor

(

2

)

q_rec SetFillColor

(

4

)

htemp2 SetFillStyle

(

3017

)

gStyle SetFrameBorderMode

(

0

)

c1 SetGridx

(

)

c1 SetGridy

(

)

gr SetLineColor

(

kBlue

)

htemp SetLineColor

(

kBlack

)

htemp2 SetLineColor

(

kRed

)

D_inj SetLineColor

(

2

)

q_rec SetLineColor

(

4

)

htemp SetLineWidth

(

3

)

htemp2 SetLineWidth

(

2

)

q_rec SetLineWidth

(

4

)

c1 SetLogx

(

1

)

c1 SetLogx

(

0

)

c1 SetLogx

(

false

)

c1 SetLogy

(

1

)

c1 SetLogy

(

0

)

D_Mtot_inj SetMarkerColor

(

2

)

D_Mtot_rec SetMarkerColor

(

4

)

rhocc SetMarkerColor

(

1

)

sim SetMarkerColor

(

kRed

)

sim SetMarkerColor

(

kBlue

)

t0 SetMarkerColor

(

kMagenta

)

t1 SetMarkerColor

(

kCyan

)

t100 SetMarkerColor

(

kYellow

)

D_Mtot_inj SetMarkerSize

(

0.

5

)

gr SetMarkerSize

(

1.

0

)

D_Mtot_inj SetMarkerStyle

(

20

)

D_Mtot_inj SetName

(

"D_Mtotinj"

)

D_Mchirp_inj SetName

(

"D_Chirp_inj"

)

D_Chi_inj SetName

(

"D_Chi_inj"

)

D_Chi_injx SetName

(

"D_Chi_injx"

)

D_Chi_injy SetName

(

"D_Chi_injy"

)

D_q_inj SetName

(

"D_q_inj"

)

efficiency SetName

(

"efficiency"

)

gStyle SetNumberContours

(

256

)

gStyle SetOptFit

(

1

)

gStyle SetOptFit

(

0

)

gStyle SetOptStat

(

kFALSE

)

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

(

1

)

gStyle SetOptStat

(

0

)

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

(

)

rho_pf SetTitle

(

"chi2"

)

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

(

title

)

q_rec SetTitle

(

0

)

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	(	fname3	,
		"%s/injected_signals.txt"	,
		netdir
	)

sprintf	(	line	,
		"#GPS@L1 mass1 mass2 distance spinx1 spiny1 spinz1 ""spinx2 spiny2 spinz2 \n"
	)

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

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

sprintf	(	inj_title	,
		"Injected events"
	)

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

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

sprintf	(	fname	,
		"%s/Dchirp_candle2.png"	,
		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_cut
	)

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	(	fname	,
		"%s/Relative_snr_Loss.png"	,
		netdir
	)

sprintf	(	lab	,
		"Injections: %i"	,
		(int) mdc.GetEntries()
	)

sprintf	(	lab	,
		"found: %i"	,
		cnt
	)

sprintf	(	lab	,
		"missed: %i"	,
		(int) mdc.GetEntries()-	cnt
	)

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

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

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

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

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

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

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

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

sprintf	(	eff_title	,
		"Efficiency"
	)

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

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

strcpy	(	cfg->	tmp_dir,
		"tmp"
	)

c1 Update

(

)

Variable Documentation

int actual_nfactor = 0

Definition at line 1482 of file cbc_plots_sim4_ifar.C.

cout<< "Total live time ---> background

Definition at line 1070 of file cbc_plots_sim4_ifar.C.

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

Definition at line 38 of file cbc_plots_sim4_ifar.C.

bool bmaxDistance

Definition at line 152 of file cbc_plots_sim4_ifar.C.

bool bmaxMass1

Definition at line 152 of file cbc_plots_sim4_ifar.C.

bool bmaxMass2

Definition at line 152 of file cbc_plots_sim4_ifar.C.

bool bmaxMtot

Definition at line 152 of file cbc_plots_sim4_ifar.C.

bool bmaxRatio

Definition at line 152 of file cbc_plots_sim4_ifar.C.

bool bminDistance

Definition at line 152 of file cbc_plots_sim4_ifar.C.

bool bminMass1

Definition at line 152 of file cbc_plots_sim4_ifar.C.

bool bminMass2

Definition at line 152 of file cbc_plots_sim4_ifar.C.

bool bminMtot

Definition at line 152 of file cbc_plots_sim4_ifar.C.

bool bminRatio

Definition at line 152 of file cbc_plots_sim4_ifar.C.

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

Definition at line 134 of file cbc_plots_sim4_ifar.C.

float chi[3]

Definition at line 808 of file cbc_plots_sim4_ifar.C.

float CYS

Initial value:

=
            86400. * 365.25

Definition at line 131 of file cbc_plots_sim4_ifar.C.

int cz

Definition at line 803 of file cbc_plots_sim4_ifar.C.

TH2F* D_Chi_injx = (TH2F *)D_Chi_inj->Clone("D_Chi_injx")

Definition at line 652 of file cbc_plots_sim4_ifar.C.

TH2F* D_Chi_injy = (TH2F *)D_Chi_inj->Clone("D_Chi_injy")

Definition at line 655 of file cbc_plots_sim4_ifar.C.

TH2F* D_Chi_rec

Initial value:

=
       new TH2F("Distance vs Chi rec.", "", 1000, MINCHI, MAXCHI, 5000,
           MINDISTANCE / 1000., MAXDISTANCE / 1000)

Definition at line 659 of file cbc_plots_sim4_ifar.C.

TH2F* D_Chi_recx = (TH2F *)D_Chi_rec->Clone("D_Chi_recx")

Definition at line 666 of file cbc_plots_sim4_ifar.C.

TH2F* D_Chi_recy = (TH2F *)D_Chi_rec->Clone("D_Chi_recy")

Definition at line 667 of file cbc_plots_sim4_ifar.C.

TH3F* D_dchirp_rec

Initial value:

=
            new TH3F("Distance vs dchirp rec.", "", 5, 10.0, 50., 100, -50, 50.,
                     5000, MINDISTANCE / 1000., MAXDISTANCE / 1000)

Definition at line 750 of file cbc_plots_sim4_ifar.C.

TH1D* D_inj = D_Mtot_inj->ProjectionY()

Definition at line 1968 of file cbc_plots_sim4_ifar.C.

TH2F* D_Mchirp_inj

Initial value:

=
            new TH2F("Distance vs MChirp inj.", "", 1000, MINCHIRP, MAXCHIRP,
                     5000, MINDISTANCE / 1000., MAXDISTANCE / 1000.)

Definition at line 582 of file cbc_plots_sim4_ifar.C.

TH2F* D_Mchirp_rec

Initial value:

=
            new TH2F("Distance vs MChirp rec.", "", 1000, MINCHIRP, MAXCHIRP,
                     5000, MINDISTANCE / 1000., MAXDISTANCE / 1000)

Definition at line 605 of file cbc_plots_sim4_ifar.C.

TH2F* D_Mtot_inj

Initial value:

=
            new TH2F("Distance vs Mtot inj.", "", 1000, MINMtot, MAXMtot, 5000,
                     MINDISTANCE / 1000., MAXDISTANCE / 1000.)

Definition at line 551 of file cbc_plots_sim4_ifar.C.

TH2F* D_Mtot_rec

Initial value:

=
            new TH2F("Distance vs Mtot rec.", "", 1000, MINMtot, MAXMtot, 5000,
                     MINDISTANCE / 1000., MAXDISTANCE / 1000.)

Definition at line 574 of file cbc_plots_sim4_ifar.C.

TH3F* D_Mtot_rec3

Initial value:

= new TH3F(
            "Distance vs Mtot rec. vs ifar", "", 100, MINMtot, MAXMtot, 1000,
            MINDISTANCE / 1000., MAXDISTANCE / 1000., 1000, 0., MAXIFAR)

Definition at line 795 of file cbc_plots_sim4_ifar.C.

TH2F* D_q_inj

Initial value:

=
            new TH2F("Distance vs q inj.", "", 1000, MINRATIO, MAXRATIO, 5000,
                     MINDISTANCE / 1000., MAXDISTANCE / 1000.)

Definition at line 669 of file cbc_plots_sim4_ifar.C.

TH2F* D_q_rec

Initial value:

=
            new TH2F("Distance vs q rec.", "", 1000, MINRATIO, MAXRATIO, 5000,
                     MINDISTANCE / 1000., MAXDISTANCE / 1000)

Definition at line 692 of file cbc_plots_sim4_ifar.C.

TH1D* D_rec = D_Mtot_rec->ProjectionY()

Definition at line 1991 of file cbc_plots_sim4_ifar.C.

TH2F* dchirp_rec

Initial value:

=
            new TH2F("dchirp rec.", "Chirp Mass estimate", 5, xq, 100, yq)

Definition at line 731 of file cbc_plots_sim4_ifar.C.

bool DDistrChirpMass

Definition at line 152 of file cbc_plots_sim4_ifar.C.

bool DDistrUniform

Definition at line 152 of file cbc_plots_sim4_ifar.C.

bool DDistrVolume

Definition at line 152 of file cbc_plots_sim4_ifar.C.

float distance

Definition at line 805 of file cbc_plots_sim4_ifar.C.

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

Definition at line 699 of file cbc_plots_sim4_ifar.C.

char eff_title[256]

Definition at line 2212 of file cbc_plots_sim4_ifar.C.

TH2F* efficiency = (TH2F *)rec_events->Clone()

Definition at line 2201 of file cbc_plots_sim4_ifar.C.

TH2F* efficiency_first_shell = (TH2F *)factor_events_rec->Clone()

Definition at line 2238 of file cbc_plots_sim4_ifar.C.

else

Initial value:

{
                cout << "Missing ifar branch: either use cbc_plots or add it "
                        "to wave tree..."
                     << endl

Definition at line 788 of file cbc_plots_sim4_ifar.C.

delete [] error_radius

Definition at line 2761 of file cbc_plots_sim4_ifar.C.

delete [] error_spin_mtot_radius

Definition at line 2772 of file cbc_plots_sim4_ifar.C.

delete [] error_spin_mtot_volume

Definition at line 2771 of file cbc_plots_sim4_ifar.C.

delete [] error_volume

Definition at line 2759 of file cbc_plots_sim4_ifar.C.

delete [] error_volume_first_shell

Definition at line 2760 of file cbc_plots_sim4_ifar.C.

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

Definition at line 2209 of file cbc_plots_sim4_ifar.C.

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

Definition at line 2326 of file cbc_plots_sim4_ifar.C.

TF1* f2

Initial value:

=
            cbcTool.doRangePlot(RHO_NBINS, Trho, Rrho, eRrho, RHO_MIN, T_cut,
                                c1, networkname, netdir, write_ascii)

Definition at line 2728 of file cbc_plots_sim4_ifar.C.

float factor

Definition at line 805 of file cbc_plots_sim4_ifar.C.

delete [] factor_events_inj = new TH2F[nfactor]

Definition at line 537 of file cbc_plots_sim4_ifar.C.

TH2F* factor_events_rec

Initial value:

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

Definition at line 547 of file cbc_plots_sim4_ifar.C.

delete [] factor_events_spin_mtot_inj

Definition at line 2739 of file cbc_plots_sim4_ifar.C.

delete [] FACTORS = new float[nfactor]

Definition at line 156 of file cbc_plots_sim4_ifar.C.

finj<< line << endl; ofstream *finj_single = new ofstream[nfactor]; TString xml; for (int l = 0; l < nfactor; l++) { sprintf(fname3, "%s/injected_signals_%d.txt", netdir, l+1); finj_single[l].open(fname3,std::ofstream::out); finj_single[l] << line << endl; } for (int g = 0; g < (int)mdc.GetEntries(); g++) { mdc.GetEntry(g); if (Redshift) { mass[0] = SFmasses[factor - 1][0]; mass[1] = SFmasses[factor - 1][1]; mchirp = pow(mass[0] * mass[1], 3. / 5.) / pow(mass[0] + mass[1], 1. / 5.); } inj_events-> Fill(mass[0], mass[1])

Definition at line 864 of file cbc_plots_sim4_ifar.C.

ofstream finj

Definition at line 827 of file cbc_plots_sim4_ifar.C.

bool FixedFiducialVolume

Definition at line 152 of file cbc_plots_sim4_ifar.C.

char fname3[2048]

Definition at line 819 of file cbc_plots_sim4_ifar.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 1671 of file cbc_plots_sim4_ifar.C.

int gIFACTOR =-1

Definition at line 141 of file cbc_plots_sim4_ifar.C.

TGraphErrors* gr

Initial value:

=
            new TGraphErrors(vV.size(), &vfar[0], &vV[0], &vefar[0], &veV[0])

Definition at line 1449 of file cbc_plots_sim4_ifar.C.

TGraphErrors* gr2

Initial value:

=
            new TGraphErrors(vR.size(), &vfar[0], &vR[0], &vefar[0], &veR[0])

Definition at line 1505 of file cbc_plots_sim4_ifar.C.

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

Definition at line 37 of file cbc_plots_sim4_ifar.C.

TH2F* h_radius

Initial value:

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

Definition at line 2288 of file cbc_plots_sim4_ifar.C.

TH1* hcandle = D_dchirp_rec->Project3D("yx")

Definition at line 1589 of file cbc_plots_sim4_ifar.C.

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

Definition at line 1759 of file cbc_plots_sim4_ifar.C.

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

Definition at line 1803 of file cbc_plots_sim4_ifar.C.

ifactor = (int)factor - 1

Definition at line 809 of file cbc_plots_sim4_ifar.C.

TH2F* inj_events

Initial value:

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

Definition at line 483 of file cbc_plots_sim4_ifar.C.

char inj_title[256]

Definition at line 1560 of file cbc_plots_sim4_ifar.C.

TString Insp = MDC.GetInspiral()

Definition at line 167 of file cbc_plots_sim4_ifar.C.

int l

Definition at line 803 of file cbc_plots_sim4_ifar.C.

char lab[256]

Definition at line 1727 of file cbc_plots_sim4_ifar.C.

leg_D

Initial value:

=
            new TLegend(0.679719, 0.156425, 0.997992, 0.327409, "", "brNDC")

Definition at line 1873 of file cbc_plots_sim4_ifar.C.

leg_D2

Initial value:

=
            new TLegend(0.679719, 0.156425, 0.997992, 0.327409, "", "brNDC")

Definition at line 1907 of file cbc_plots_sim4_ifar.C.

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

Definition at line 1728 of file cbc_plots_sim4_ifar.C.

char line[1024]

Definition at line 820 of file cbc_plots_sim4_ifar.C.

int m

Definition at line 803 of file cbc_plots_sim4_ifar.C.

float M1

Definition at line 2342 of file cbc_plots_sim4_ifar.C.

float M2

Definition at line 2342 of file cbc_plots_sim4_ifar.C.

float mass[2]

Definition at line 806 of file cbc_plots_sim4_ifar.C.

int massbins = 0

Definition at line 2244 of file cbc_plots_sim4_ifar.C.

int MAX_AXIS_Z

Initial value:

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

Definition at line 1556 of file cbc_plots_sim4_ifar.C.

float max_mass1 = MAX_MASS

Definition at line 62 of file cbc_plots_sim4_ifar.C.

float max_mass2 = MAX_MASS

Definition at line 50 of file cbc_plots_sim4_ifar.C.

delete [] maxDistance = new float[nfactor]

Definition at line 149 of file cbc_plots_sim4_ifar.C.

delete [] maxDistanceXML = new float[nfactor]

Definition at line 148 of file cbc_plots_sim4_ifar.C.

double MAXIFAR = 0.0

Definition at line 784 of file cbc_plots_sim4_ifar.C.

delete [] maxMChirp = new float[nfactor]

Definition at line 147 of file cbc_plots_sim4_ifar.C.

delete [] maxMtot = new float[nfactor]

Definition at line 146 of file cbc_plots_sim4_ifar.C.

delete [] maxRatio = new float[nfactor]

Definition at line 151 of file cbc_plots_sim4_ifar.C.

float mchirp

Definition at line 805 of file cbc_plots_sim4_ifar.C.

TH1D* Mchirp_inj = D_Mchirp_inj->ProjectionX()

Definition at line 2049 of file cbc_plots_sim4_ifar.C.

TH1D* Mchirp_rec = D_Mchirp_rec->ProjectionX()

Definition at line 2072 of file cbc_plots_sim4_ifar.C.

int mcount = 0

Definition at line 1419 of file cbc_plots_sim4_ifar.C.

CWB::mdc MDC(net)

float min_mass1 = MIN_MASS

Definition at line 56 of file cbc_plots_sim4_ifar.C.

float min_mass2 = MIN_MASS

Definition at line 68 of file cbc_plots_sim4_ifar.C.

bool minchi

Definition at line 152 of file cbc_plots_sim4_ifar.C.

delete [] minDistance = new float[nfactor]

Definition at line 149 of file cbc_plots_sim4_ifar.C.

delete [] minDistanceXML = new float[nfactor]

Definition at line 148 of file cbc_plots_sim4_ifar.C.

delete [] minMChirp = new float[nfactor]

Definition at line 147 of file cbc_plots_sim4_ifar.C.

delete [] minMtot = new float[nfactor]

Definition at line 146 of file cbc_plots_sim4_ifar.C.

delete [] minRatio = new float[nfactor]

Definition at line 150 of file cbc_plots_sim4_ifar.C.

TH1D* Mtot_inj = D_Mtot_inj->ProjectionX()

Definition at line 2015 of file cbc_plots_sim4_ifar.C.

TH1D* Mtot_rec = D_Mtot_rec->ProjectionX()

Definition at line 2037 of file cbc_plots_sim4_ifar.C.

int mtt

Definition at line 803 of file cbc_plots_sim4_ifar.C.

double mytime[6]

Definition at line 804 of file cbc_plots_sim4_ifar.C.

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

Definition at line 45 of file cbc_plots_sim4_ifar.C.

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

Definition at line 71 of file cbc_plots_sim4_ifar.C.

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

Definition at line 72 of file cbc_plots_sim4_ifar.C.

int NBINS_SPIN = (int)((MAXCHI - MINCHI) / CHI_BIN)

Definition at line 613 of file cbc_plots_sim4_ifar.C.

network* net = NULL

Definition at line 133 of file cbc_plots_sim4_ifar.C.

char networkname[256]

Definition at line 113 of file cbc_plots_sim4_ifar.C.

NEVTS = 0.0

Definition at line 810 of file cbc_plots_sim4_ifar.C.

Int_t nGraphs = 0

Definition at line 2345 of file cbc_plots_sim4_ifar.C.

int nmdc = mdc.GetSelectedRows()

Definition at line 1668 of file cbc_plots_sim4_ifar.C.

int nwave = sim.GetSelectedRows()

Definition at line 1703 of file cbc_plots_sim4_ifar.C.

int nwave_final = sim.GetSelectedRows()

Definition at line 1721 of file cbc_plots_sim4_ifar.C.

TPaveText* p_eff

Initial value:

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

Definition at line 2219 of file cbc_plots_sim4_ifar.C.

TPaveText* p_inj

Initial value:

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

Definition at line 1567 of file cbc_plots_sim4_ifar.C.

p_radius

Initial value:

=
            new TPaveText(0.125301, 0.926166, 0.567068, 0.997409, "blNDC")

Definition at line 1999 of file cbc_plots_sim4_ifar.C.

int point

Definition at line 2349 of file cbc_plots_sim4_ifar.C.

TH1D* q_inj = D_q_inj->ProjectionX()

Definition at line 2159 of file cbc_plots_sim4_ifar.C.

TH1D* q_rec = D_q_rec->ProjectionX()

Definition at line 2182 of file cbc_plots_sim4_ifar.C.

delete [] radius

Definition at line 2758 of file cbc_plots_sim4_ifar.C.

char radius_title[256]

Definition at line 2329 of file cbc_plots_sim4_ifar.C.

cout<< "MDC set: " << gIFACTOR << endl; cout << "xml conf: waveform=" << waveforms[gIFACTOR - 1] << " minMtot=" << minMtot[gIFACTOR - 1] << " maxMtot=" << maxMtot[gIFACTOR - 1] << " minDistance=" << minDistance[gIFACTOR - 1] << " maxDistance=" << maxDistance[gIFACTOR - 1] << " minRatio=" << minRatio[gIFACTOR - 1] << " maxRatio=" << maxRatio[gIFACTOR - 1] << endl; } float MINMtot = 0.0; if (bminMtot) { float MINMtot = 0.99 * minMtot[gIFACTOR - 1]; } float MAXMtot = 0.0; int NBINS_MTOT = 0; if (bmaxMtot) { MAXMtot = MAX_MASS; NBINS_MTOT = TMath::FloorNint((MAX_MASS - MIN_MASS) / MASS_BIN / 2.); cout << "NBINS_MTOT: " << NBINS_MTOT << endl; } else { cout << "Undefined maximal total mass!! Define float MAXMtot" << endl; exit(1); } float MINDISTANCE = 0.0; if (bminDistance) { MINDISTANCE = 0.9 * ShellminDistance; cout << "MINDISTANCE = " << MINDISTANCE << endl; } else { cout << "Undefined minimal distance!! Defined float MINDISTANCE" << endl; } float MAXDISTANCE = 0.0; if (bmaxDistance) { MAXDISTANCE = TMath::Max(ShellmaxDistance, maxDistance[gIFACTOR - 1]); cout << "MAXDISTANCE = " << MAXDISTANCE << endl; } else { cout << "Undefined maximal distance!! Define float MAXDISTANCE" << endl; exit(1); } float MINCHIRP = 100.0; float MAXCHIRP = 0.0; float MINRATIO = 1.0; float MAXRATIO = 1.0; if ((bminRatio) && (bmaxRatio)) { for (int l = 0; l < nfactor; l++) { if (MINRATIO > minRatio[l]) { MINRATIO = minRatio[l]; } if (MAXRATIO < maxRatio[l]) { MAXRATIO = maxRatio[l]; } if (MINCHIRP > minMtot[l] * pow(maxRatio[l], 3. / 5.) / pow(1 + minRatio[l], 6. / 5.)) { MINCHIRP = minMtot[l] * pow(maxRatio[l], 3. / 5.) / pow(1 + minRatio[l], 6. / 5.); }; if (MAXCHIRP < maxMtot[l] * pow(minRatio[l], 3. / 5.) / pow(1 + minRatio[l], 6. / 5.)) { MAXCHIRP = maxMtot[l] * pow(minRatio[l], 3. / 5.) / pow(1 + minRatio[l], 6. / 5.); }; } } else { cout << "Undefined minRatio.. " << endl; exit(1); } for (int l = 0; l < nfactor - 1; l++) { if ((minDistanceXML[l] == minDistanceXML[l + 1]) && (maxDistanceXML[l] == maxDistanceXML[l + 1])) { FixedFiducialVolume = 1; } else { FixedFiducialVolume = 1; cout << "Beware: different fiducial volumes for " "different factors!!" << endl; } } if (strlen(veto_not_vetoed) == 0) { sprintf(veto_not_vetoed, "%s", ch2); } cout << "ERROR! Missing zero lag livetime...Please add:" << endl; cout << "\"#define LIVE_ZERO XXXXXXXXXX\" where XXXXXXXXXX is the " "livetime in seconds to your config/user_pparameters.C file...." << endl; exit(1); TCanvas *c1 = new TCanvas("c1", "c1", 3, 47, 1000, 802); c1-> Range(-1.216392,-477.6306, 508.8988, 2814.609)

Definition at line 473 of file cbc_plots_sim4_ifar.C.

TH2F* rec_events

Initial value:

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

Definition at line 512 of file cbc_plots_sim4_ifar.C.

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

Definition at line 36 of file cbc_plots_sim4_ifar.C.

bool Redshift

Definition at line 152 of file cbc_plots_sim4_ifar.C.

TH2F* rho_pf

Initial value:

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

Definition at line 715 of file cbc_plots_sim4_ifar.C.

TH2F* rhocc

Initial value:

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

Definition at line 704 of file cbc_plots_sim4_ifar.C.

delete [] shell_volume = new float[nfactor]

Definition at line 151 of file cbc_plots_sim4_ifar.C.

float ShellmaxDistance = 0.0

Definition at line 158 of file cbc_plots_sim4_ifar.C.

float ShellminDistance = 9999999.0

Definition at line 157 of file cbc_plots_sim4_ifar.C.

float spin[6]

Definition at line 807 of file cbc_plots_sim4_ifar.C.

delete [] spin_mtot_radius

Definition at line 2770 of file cbc_plots_sim4_ifar.C.

delete [] spin_mtot_volume

Definition at line 2769 of file cbc_plots_sim4_ifar.C.

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

Definition at line 35 of file cbc_plots_sim4_ifar.C.

TH1* t0 = D_Mtot_rec3->Project3D("yx_0")

Definition at line 1892 of file cbc_plots_sim4_ifar.C.

TH1* t1 = D_Mtot_rec3->Project3D("yx_1")

Definition at line 1896 of file cbc_plots_sim4_ifar.C.

TH1* t10 = D_Mtot_rec3->Project3D("yx_10")

Definition at line 1900 of file cbc_plots_sim4_ifar.C.

TH1* t100 = D_Mtot_rec3->Project3D("yx_100")

Definition at line 1904 of file cbc_plots_sim4_ifar.C.

text = p_inj->AddText(inj_title)

Definition at line 1574 of file cbc_plots_sim4_ifar.C.

TLatex Tl

Definition at line 2347 of file cbc_plots_sim4_ifar.C.

TLatex Tl2

Definition at line 2348 of file cbc_plots_sim4_ifar.C.

Int_t TotalConts = 0

Definition at line 2346 of file cbc_plots_sim4_ifar.C.

float Tscale = 1.0

Definition at line 1481 of file cbc_plots_sim4_ifar.C.

double V0 = 0.0

Definition at line 2245 of file cbc_plots_sim4_ifar.C.

char val[20]

Definition at line 2344 of file cbc_plots_sim4_ifar.C.

std::vector<double> vdv

Definition at line 1076 of file cbc_plots_sim4_ifar.C.

std::vector<double> vefar

Definition at line 1410 of file cbc_plots_sim4_ifar.C.

std::vector<double> veR

Definition at line 1480 of file cbc_plots_sim4_ifar.C.

std::vector<double> veV

Definition at line 1407 of file cbc_plots_sim4_ifar.C.

std::vector<double> vfar

Definition at line 1409 of file cbc_plots_sim4_ifar.C.

std::vector<double> vifar

Definition at line 1077 of file cbc_plots_sim4_ifar.C.

delete [] volume

Definition at line 2756 of file cbc_plots_sim4_ifar.C.

delete [] volume_first_shell

Definition at line 2757 of file cbc_plots_sim4_ifar.C.

std::vector<double> vR

Definition at line 1479 of file cbc_plots_sim4_ifar.C.

std::vector<double> vsifar

Definition at line 1408 of file cbc_plots_sim4_ifar.C.

cout<< endl; if (write_ascii) { fev.close(); for (int l = 0; l < nfactor; l++) { fev_single[l].close(); } } cout << "Recovered entries: " << cnt << endl; cout << "Recovered entries: " << cnt2 << endl; cout << "Recovered entries cut by frequency: " << cntfreq << endl; cout << "Recovered entries vetoed: " << countv << endl; cout << "dV : " << dV << " dV1 : " << dV1 << endl; internal_volume = 4. / 3. * TMath::Pi() * pow(minDistance[0] / 1000., 3.); if (INCLUDE_INTERNAL_VOLUME) { for (int i = 0; i < vdv.size(); i++) { if (vdv[i] > 0.0) { vdv[i] += internal_volume; } } for (int i = 0; i < RHO_NBINS; i++) { if (Vrho[i] > 0.0) { Vrho[i] += internal_volume; } } for (int i = 0; i < NBINS_MTOT + 1; i++) { for (int j = 0; j < NBINS_SPIN + 1; j++) { if (spin_mtot_volume[i][j] > 0.0) { spin_mtot_volume[i][j] += internal_volume; } } } for (int i = 0; i < NBINS_mass; i++) { for (int j = 0; j < NBINS_mass2; j++) { if (volume[i][j] > 0.0) { volume[i][j] += internal_volume; } } } } Int_t *mindex = new Int_t[vdv.size()]; TMath::Sort((int)vdv.size(), &vifar[0], mindex, true); std::vector<double> vV

Definition at line 1406 of file cbc_plots_sim4_ifar.C.

delete [] waveforms = new TString[nfactor]

Definition at line 155 of file cbc_plots_sim4_ifar.C.

bool write_ascii = false

Definition at line 818 of file cbc_plots_sim4_ifar.C.

Double_t x0

Definition at line 2343 of file cbc_plots_sim4_ifar.C.

Float_t xq[6] = {8.0, 15.0, 25.0, 35.0, 45.0, 55.0}

Definition at line 726 of file cbc_plots_sim4_ifar.C.

Double_t Y0

Definition at line 2343 of file cbc_plots_sim4_ifar.C.

Float_t* yq = new Float_t[101]

Definition at line 727 of file cbc_plots_sim4_ifar.C.

Double_t z0

Definition at line 2343 of file cbc_plots_sim4_ifar.C.