cbc_plots_sim4.C File Reference

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

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 ()
cout<< endl;if(write_ascii){fclose(fev);for(int l=0;l< nfactor;l++){fclose(fev_single[l]);}}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;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 ()
hcandle	Draw ("CANDLE")
mdc	Draw (sel,"")
sim	Draw (sel, newcut,"")
sim	Draw (sel, newcut2,"same")
sim	Draw (sel, newcut2)
	exit (0)
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[0], mass[1])
D_dchirp_rec	Fill (chirp[0], chirp[0]-chirp[1], range[1])
dchirp_rec	Fill (chirp[0], chirp[0]-chirp[1])
	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 (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("--d-distr")!="")
	if (MDC.GetInspiralOption("--dchirp-distr")!="")
	if (write_ascii)
	if (range[1]==0.0)
	if ((l >=NBINS_mass2)||(l< 0))
	if ((m >=NBINS_mass)||(m< 0))
	if (DDistrVolume)
else	if (DDistrUniform)
else	if (DDistrChirpMass)
	if (FixedFiducialVolume)
	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")
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 (0)
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/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 (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
CWB::config *	cfg = new CWB::config
float	chi [3]
float	chi2 = penalty>0 ? log10(penalty) : 0
float	CYS = 31560000.0
int	cz
TH3F *	D_dchirp_rec = new TH3F("Distance vs dchirp rec.","",5,10.0,50.,100,-50,50.,.5000,MINDISTANCE/1000.,MAXDISTANCE/1000)
TH2F *	D_Mchirp_inj = new TH2F("Distance vs MChirp inj.","",1000,MINCHIRP,MAXCHIRP,5000,MINDISTANCE/1000.,MAXDISTANCE/1000.)
TH2F *	D_Mchirp_rec = new TH2F("Distance vs MChirp rec.","",1000,MINCHIRP,MAXCHIRP,5000,MINDISTANCE/1000.,MAXDISTANCE/1000)
TH2F *	D_Mtot_inj = new TH2F("Distance vs Mtot inj.","",1000,MINMtot,MAXMtot,5000,MINDISTANCE/1000.,MAXDISTANCE/1000.)
TH2F *	D_Mtot_rec = new TH2F("Distance vs Mtot rec.","",1000,MINMtot,MAXMtot,5000,MINDISTANCE/1000.,MAXDISTANCE/1000.)
TH2F *	D_q_inj = new TH2F("Distance vs q inj.","",1000,MINRATIO,MAXRATIO,5000,MINDISTANCE/1000.,MAXDISTANCE/1000.)
TH2F *	D_q_rec = new TH2F("Distance vs q rec.","",1000,MINRATIO,MAXRATIO,5000,MINDISTANCE/1000.,MAXDISTANCE/1000)
TH2F *	dchirp_rec = new TH2F("dchirp rec.","Chirp Mass estimate",5,xq,100,yq)
bool	DDistrChirpMass
bool	DDistrUniform
bool	DDistrVolume =1
float	distance
TH1F *	Dt = new TH1F("Dt","",1000,-0.5,0.5)
	dV = shell_volume[ifactor]
	Temporay patch for redshifted mass distributions, i.e. point-like in source frame and spread over multiple bins in the detector frame. More...
	dV1 = dV/nevts
	else {cout << "No defined distance distribution?????! WARNING: Assuming uniform in volume"<<endl
	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_cut, 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)
float	FACTORS [nfactor]
cout<< "Injected signals: " << mdc.GetEntries()<< endl;char cut[512];sprintf(cut,"rho[%d]> %f && %s", pp_irho, T_cut, ch2);float 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("ecor",&ecor);sim.SetBranchAddress("penalty",&penalty);sim.SetBranchAddress("time", time);sim.SetBranchAddress("iSNR", iSNR);sim.SetBranchAddress("sSNR", sSNR);sim.SetBranchAddress("spin", spin);sim.SetBranchAddress("frequency", frequency);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;char fname[1024];if(write_ascii){sprintf(fname,"%s/recovered_signals.txt", netdir);FILE *fev=fopen(fname,"w");fprintf(fev,"#GPS@L1 factor rho frequency iSNR sSNR \n");FILE *fev_single[nfactor];for(int l=1;l< nfactor+1;l++){sprintf(fname,"%s/recovered_signals_%d.txt", netdir, l);fev_single[l-1]=fopen(fname,"w");fprintf(fev_single[l-1],"#GPS@L1 factor rho frequency iSNR sSNR \n");}}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, nevts, internal_volume;int nT;int countv=0;int cntfreq=0;bool bcut=false;for(int g=0;g<(int) sim.GetEntries();g++){sim.GetEntry(g);if(rho[pp_irho] <=T_cut){countv++;continue;}if(netcc[0] <=T_cor){countv++;continue;}if((time[0]-time[nIFO]) <-T_win||(time[0]-time[nIFO])> 2 *T_win){countv++;continue;}if(T_vED > 0){if(neted/ecor >=T_vED){countv++;continue;}}if(T_pen > 0){if(penalty<=T_vED){countv++;continue;}}if(T_pen > 0){if(penalty<=T_vED){countv++;continue;}}bcut=false;for(int j=0;j< nFCUT;j++){if((frequency[0]> lowFCUT[j])&&(frequency[0] <=highFCUT[j])) bcut=true;}if(bcut){countv++;cntfreq++;continue;}if(++cnt%1000==0){cout << cnt<< " - ";}Dt->	Fill (time[0]-time[nIFO])
bool	FixedFiducialVolume
cout<< "nmdc : "<< nmdc<< endl;TH2F *htemp=(TH2F *) gPad->	GetPrimitive ("hist")
int	gIFACTOR =0
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 *	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 = 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[0]
	m2 =mass[1]
float	mass [2]
int	MAX_AXIS_Z = inj_events->GetBinContent(inj_events->GetMaximumBin()) + 1
float	max_mass1 = MAX_MASS
float	max_mass2 = MAX_MASS
float	maxDistance [nfactor]
float	maxDistanceXML [nfactor]
float	maxMChirp [nfactor]
float	maxMtot [nfactor]
float	maxRatio [nfactor]
float	mchirp
float	min_mass1 = MIN_MASS
float	min_mass2 = MIN_MASS
float	minDistance [nfactor]
float	minDistanceXML [nfactor]
float	minMChirp [nfactor]
float	minMtot [nfactor]
float	minRatio [nfactor]
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]
double	NEVTS =0.0
	nevts =0
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]
cout<<"MDC set: "<< gIFACTOR << endl;cout<<"xml conf: waveform="<< waveform[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;}if(bminMtot){float MINMtot=0.99 *minMtot[gIFACTOR-1];}else{float MINMtot=0.0;}if(bmaxMtot){float MAXMtot=MAX_MASS;int 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);}if(bminDistance){float MINDISTANCE=0.9 *ShellminDistance;}else{cout <<"Undefined minimal distance!! Defined float MINDISTANCE"<< endl;float MINDISTANCE=0.0;}if(bmaxDistance){float MAXDISTANCE=1.1 *ShellmaxDistance;}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=0;cout <<"Beware: different fiducial volumes for different factors!!"<< endl;exit(1);}}cout <<"Plotting bounds: MINMtot=" << MINMtot<<" MAXMtot=" << MAXMtot<<" MINRATIO=" << MINRATIO<<" MAXRATIO=" << MAXRATIO<<" MINDISTANCE=" << MINDISTANCE<<" MAXDISTANCE=" << MAXDISTANCE<<" MINCHIRP=" << MINCHIRP<<" MAXCHIRP=" << MAXCHIRP<< endl;if(strlen(veto_not_vetoed)==0){sprintf(veto_not_vetoed,"%s", ch2);}TCanvas *c1=new TCanvas("c1","c1", 3, 47, 1000, 802);c1->	Range (-1.216392,-477.6306, 508.8988, 2814.609)
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 }
bool	Redshift
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	shell_volume [nfactor]
float	ShellmaxDistance =0.0
float	ShellminDistance =9999999.0
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]
	volume [m][l] = dV
TString	waveform [nfactor]
bool	write_ascii = false
const Float_t	xq [6] = {8.0, 15.0, 25.0, 35.0, 45.0, 55.0}
const Float_t *	yq = new Float_t[101]

Macro Definition Documentation

#define CONTOURS   7

Definition at line 5 of file cbc_plots_sim4.C.

#define LIV_FILE_NAME   "live.txt"

Definition at line 8 of file cbc_plots_sim4.C.

#define NCont   99

#define NRGBs   6

#define RHO_BIN   0.1

Definition at line 3 of file cbc_plots_sim4.C.

#define RHO_MIN   5.0

Definition at line 2 of file cbc_plots_sim4.C.

#define RHO_NBINS   5000

Definition at line 4 of file cbc_plots_sim4.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.

cout<< endl; if(write_ascii){ fclose(fev); for(int l=0;l<nfactor;l++){fclose(fev_single[l]);} } 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; 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

(

)

dchirp_rec Draw

(

"CANDLE"

)

mdc Draw	(	sel	,
		""
	)

sim Draw	(	sel	,
		newcut	,
		""
	)

sim Draw	(	sel	,
		newcut2	,
		"same"
	)

sim Draw	(	sel	,
		newcut2
	)

exit

(

0

)

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	[0],
		mass	[1]
	)

D_dchirp_rec Fill	(	chirp	[0],
		chirp-	chirp[0][1],
		range	[1]
	)

dchirp_rec Fill	(	chirp	[0],
		chirp-	chirp[0][1]
	)

for

(

)

Definition at line 313 of file cbc_plots_sim4.C.

for

(

int

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

)

Definition at line 1851 of file cbc_plots_sim4.C.

gROOT ForceStyle

(

)

sim GetHistogram

(

)

Definition at line 1394 of file cbc_plots_sim4.C.

inj_events GetXaxis

(

)

-> SetRangeUser(MIN_plot_mass1, MAX_plot_mass1)

Definition at line 261 of file cbc_plots_sim4.C.

inj_events GetYaxis

(

)

-> SetRangeUser(MIN_plot_mass2, MAX_plot_mass2)

Definition at line 262 of file cbc_plots_sim4.C.

inj_events GetZaxis

(

)

-> SetTickLength(0.01)

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

Definition at line 95 of file cbc_plots_sim4.C.

if

(

MDC.GetInspiralOption("--waveform")!

= ""

)

Definition at line 122 of file cbc_plots_sim4.C.

if

(

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

= ""

)

Definition at line 123 of file cbc_plots_sim4.C.

if

(

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

= ""

)

Definition at line 124 of file cbc_plots_sim4.C.

if

(

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

= ""

)

Definition at line 125 of file cbc_plots_sim4.C.

if	(	ShellminDistance	,
		minDistanceXML	[gIFACTOR-1]
	)

Definition at line 126 of file cbc_plots_sim4.C.

if

(

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

= ""

)

Definition at line 127 of file cbc_plots_sim4.C.

if

(

)

Definition at line 128 of file cbc_plots_sim4.C.

if

(

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

= ""

)

Definition at line 129 of file cbc_plots_sim4.C.

if

(

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

= ""

)

Definition at line 130 of file cbc_plots_sim4.C.

if

(

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

= ""

)

Definition at line 131 of file cbc_plots_sim4.C.

if

(

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

= ""

)

Definition at line 144 of file cbc_plots_sim4.C.

if

(

write_ascii

)

Definition at line 741 of file cbc_plots_sim4.C.

if

(

range

[1] = =0.0

)

Definition at line 890 of file cbc_plots_sim4.C.

if

(

(l >=NBINS_mass2)||(l< 0)

)

Definition at line 903 of file cbc_plots_sim4.C.

if

(

(m >=NBINS_mass)||(m< 0)

)

Definition at line 905 of file cbc_plots_sim4.C.

if

(

DDistrVolume

)

Definition at line 927 of file cbc_plots_sim4.C.

else if

(

DDistrUniform

)

Definition at line 932 of file cbc_plots_sim4.C.

else if

(

DDistrChirpMass

)

Definition at line 937 of file cbc_plots_sim4.C.

if

(

FixedFiducialVolume

)

Definition at line 947 of file cbc_plots_sim4.C.

if

(

Redshift

)

Definition at line 971 of file cbc_plots_sim4.C.

if

(

if() c1->SetLogy(kTRUE)

pp_rho_log

)

Definition at line 1229 of file cbc_plots_sim4.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"

)

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

(

0

)

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

(

)

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/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	(	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 25 of file cbc_plots_sim4.C.

bool bmaxDistance

Definition at line 110 of file cbc_plots_sim4.C.

bool bmaxMtot

Definition at line 110 of file cbc_plots_sim4.C.

bool bmaxRatio

Definition at line 110 of file cbc_plots_sim4.C.

bool bminDistance

Definition at line 110 of file cbc_plots_sim4.C.

bool bminMtot

Definition at line 110 of file cbc_plots_sim4.C.

bool bminRatio

Definition at line 110 of file cbc_plots_sim4.C.

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

Definition at line 105 of file cbc_plots_sim4.C.

float chi[3]

Definition at line 678 of file cbc_plots_sim4.C.

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

Definition at line 897 of file cbc_plots_sim4.C.

float CYS = 31560000.0

Definition at line 103 of file cbc_plots_sim4.C.

int cz

Definition at line 673 of file cbc_plots_sim4.C.

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

Definition at line 496 of file cbc_plots_sim4.C.

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

Definition at line 347 of file cbc_plots_sim4.C.

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

Definition at line 368 of file cbc_plots_sim4.C.

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

Definition at line 320 of file cbc_plots_sim4.C.

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

Definition at line 341 of file cbc_plots_sim4.C.

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

Definition at line 423 of file cbc_plots_sim4.C.

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

Definition at line 444 of file cbc_plots_sim4.C.

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

Definition at line 478 of file cbc_plots_sim4.C.

bool DDistrChirpMass

Definition at line 110 of file cbc_plots_sim4.C.

bool DDistrUniform

Definition at line 110 of file cbc_plots_sim4.C.

DDistrVolume =1

Definition at line 110 of file cbc_plots_sim4.C.

float distance

Definition at line 675 of file cbc_plots_sim4.C.

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

Definition at line 449 of file cbc_plots_sim4.C.

dV = shell_volume[ifactor]

Temporay patch for redshifted mass distributions, i.e. point-like in source frame and spread over multiple bins in the detector frame.

Definition at line 945 of file cbc_plots_sim4.C.

dV1 = dV/nevts

Definition at line 955 of file cbc_plots_sim4.C.

else {cout << "No defined distance distribution?????! WARNING: Assuming uniform in volume"<<endl

Definition at line 943 of file cbc_plots_sim4.C.

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

Definition at line 1013 of file cbc_plots_sim4.C.

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

Definition at line 1865 of file cbc_plots_sim4.C.

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

Definition at line 2224 of file cbc_plots_sim4.C.

float factor

Definition at line 675 of file cbc_plots_sim4.C.

TH2F* factor_events_inj[nfactor]

Definition at line 311 of file cbc_plots_sim4.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 317 of file cbc_plots_sim4.C.

float FACTORS[nfactor]

Definition at line 112 of file cbc_plots_sim4.C.

cout<< "Injected signals: " << mdc.GetEntries() << endl; char cut[512]; sprintf(cut,"rho[%d]>%f && %s",pp_irho,T_cut,ch2); float 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("ecor",&ecor); sim.SetBranchAddress("penalty",&penalty); sim.SetBranchAddress("time",time); sim.SetBranchAddress("iSNR",iSNR); sim.SetBranchAddress("sSNR",sSNR); sim.SetBranchAddress("spin",spin); sim.SetBranchAddress("frequency",frequency); 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; char fname[1024]; if(write_ascii){ sprintf(fname,"%s/recovered_signals.txt",netdir); FILE* fev = fopen(fname,"w"); fprintf(fev,"#GPS@L1 factor rho frequency iSNR sSNR \n"); FILE *fev_single[nfactor]; for(int l=1;l<nfactor+1;l++){ sprintf(fname,"%s/recovered_signals_%d.txt",netdir,l); fev_single[l-1] = fopen(fname,"w"); fprintf(fev_single[l-1],"#GPS@L1 factor rho frequency iSNR sSNR \n"); }} 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,nevts,internal_volume; int nT; int countv=0; int cntfreq=0; bool bcut=false; for(int g=0;g<(int)sim.GetEntries();g++){ sim.GetEntry(g); if(rho[pp_irho]<=T_cut){countv++;continue;} if(netcc[0]<=T_cor){countv++;continue;} if((time[0]-time[nIFO])<-T_win || (time[0]-time[nIFO])>2*T_win) {countv++;continue;} if (T_vED>0) {if(neted/ecor>=T_vED){countv++;continue;}} if (T_pen>0) {if(penalty<=T_vED){countv++;continue;}} if (T_pen>0) {if(penalty<=T_vED){countv++;continue;}} bcut=false; for(int j=0;j<nFCUT;j++) { if((frequency[0]>lowFCUT[j])&&(frequency[0]<=highFCUT[j])) bcut=true; } if(bcut) {countv++;cntfreq++;continue;} if (++cnt%1000==0) {cout << cnt << " - ";} Dt-> Fill(time[0]-time[nIFO])

Definition at line 881 of file cbc_plots_sim4.C.

bool FixedFiducialVolume

Definition at line 110 of file cbc_plots_sim4.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 1276 of file cbc_plots_sim4.C.

int gIFACTOR =0

Definition at line 115 of file cbc_plots_sim4.C.

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

Definition at line 24 of file cbc_plots_sim4.C.

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

Definition at line 1828 of file cbc_plots_sim4.C.

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

Definition at line 1201 of file cbc_plots_sim4.C.

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

Definition at line 1355 of file cbc_plots_sim4.C.

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

Definition at line 1391 of file cbc_plots_sim4.C.

ifactor = (int)factor-1

Definition at line 679 of file cbc_plots_sim4.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 258 of file cbc_plots_sim4.C.

char inj_title[256]

Definition at line 1173 of file cbc_plots_sim4.C.

TString Insp = MDC.GetInspiral()

Definition at line 121 of file cbc_plots_sim4.C.

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

Definition at line 673 of file cbc_plots_sim4.C.

char lab[256]

Definition at line 1326 of file cbc_plots_sim4.C.

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

Definition at line 1327 of file cbc_plots_sim4.C.

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

Definition at line 673 of file cbc_plots_sim4.C.

m1 =mass[0]

Definition at line 900 of file cbc_plots_sim4.C.

m2 =mass[1]

Definition at line 901 of file cbc_plots_sim4.C.

float mass[2]

Definition at line 676 of file cbc_plots_sim4.C.

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

Definition at line 1169 of file cbc_plots_sim4.C.

float max_mass1 = MAX_MASS

Definition at line 50 of file cbc_plots_sim4.C.

float max_mass2 = MAX_MASS

Definition at line 38 of file cbc_plots_sim4.C.

float maxDistance[nfactor]

Definition at line 109 of file cbc_plots_sim4.C.

float maxDistanceXML[nfactor]

Definition at line 109 of file cbc_plots_sim4.C.

float maxMChirp[nfactor]

Definition at line 109 of file cbc_plots_sim4.C.

float maxMtot[nfactor]

Definition at line 109 of file cbc_plots_sim4.C.

float maxRatio[nfactor]

Definition at line 109 of file cbc_plots_sim4.C.

float mchirp

Definition at line 675 of file cbc_plots_sim4.C.

float min_mass1 = MIN_MASS

Definition at line 44 of file cbc_plots_sim4.C.

float min_mass2 = MIN_MASS

Definition at line 56 of file cbc_plots_sim4.C.

float minDistance[nfactor]

Definition at line 109 of file cbc_plots_sim4.C.

float minDistanceXML[nfactor]

Definition at line 109 of file cbc_plots_sim4.C.

float minMChirp[nfactor]

Definition at line 109 of file cbc_plots_sim4.C.

float minMtot[nfactor]

Definition at line 109 of file cbc_plots_sim4.C.

float minRatio[nfactor]

Definition at line 109 of file cbc_plots_sim4.C.

int mt

Definition at line 673 of file cbc_plots_sim4.C.

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

Definition at line 32 of file cbc_plots_sim4.C.

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

Definition at line 61 of file cbc_plots_sim4.C.

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

Definition at line 62 of file cbc_plots_sim4.C.

network* net =NULL

Definition at line 104 of file cbc_plots_sim4.C.

char networkname[256]

Definition at line 94 of file cbc_plots_sim4.C.

NEVTS =0.0

Definition at line 680 of file cbc_plots_sim4.C.

nevts =0

Definition at line 969 of file cbc_plots_sim4.C.

int nmdc = mdc.GetSelectedRows()

Definition at line 1273 of file cbc_plots_sim4.C.

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

Definition at line 982 of file cbc_plots_sim4.C.

int nwave = sim.GetSelectedRows()

Definition at line 1302 of file cbc_plots_sim4.C.

int nwave_final = sim.GetSelectedRows()

Definition at line 1320 of file cbc_plots_sim4.C.

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

Definition at line 1180 of file cbc_plots_sim4.C.

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

Definition at line 1874 of file cbc_plots_sim4.C.

char radius_title[256]

Definition at line 1869 of file cbc_plots_sim4.C.

cout<<"MDC set: "<<gIFACTOR<<endl; cout<<"xml conf: waveform="<< waveform[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; } if(bminMtot){float MINMtot = 0.99*minMtot[gIFACTOR-1];} else {float MINMtot = 0.0;} if(bmaxMtot){ float MAXMtot = MAX_MASS; int 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);} if(bminDistance) {float MINDISTANCE = 0.9*ShellminDistance;} else {cout <<"Undefined minimal distance!! Defined float MINDISTANCE"<<endl;float MINDISTANCE = 0.0;} if(bmaxDistance) {float MAXDISTANCE = 1.1*ShellmaxDistance;} 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=0; cout<<"Beware: different fiducial volumes for different factors!!"<<endl; exit(1); } }cout<<"Plotting bounds: MINMtot="<<MINMtot<<" MAXMtot="<<MAXMtot<<" MINRATIO="<< MINRATIO <<" MAXRATIO="<<MAXRATIO<<" MINDISTANCE="<<MINDISTANCE<<" MAXDISTANCE="<<MAXDISTANCE<<" MINCHIRP="<<MINCHIRP<<" MAXCHIRP="<<MAXCHIRP<<endl; if (strlen(veto_not_vetoed) == 0){sprintf(veto_not_vetoed,"%s",ch2);} TCanvas *c1 = new TCanvas("c1","c1",3,47,1000,802); c1-> Range(-1.216392,-477.6306, 508.8988, 2814.609)

Definition at line 247 of file cbc_plots_sim4.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 287 of file cbc_plots_sim4.C.

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

Definition at line 23 of file cbc_plots_sim4.C.

bool Redshift

Definition at line 110 of file cbc_plots_sim4.C.

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

Definition at line 465 of file cbc_plots_sim4.C.

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

Definition at line 455 of file cbc_plots_sim4.C.

float shell_volume[nfactor]

Definition at line 109 of file cbc_plots_sim4.C.

float ShellmaxDistance =0.0

Definition at line 114 of file cbc_plots_sim4.C.

float ShellminDistance =9999999.0

Definition at line 113 of file cbc_plots_sim4.C.

float spin[6]

Definition at line 677 of file cbc_plots_sim4.C.

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

Definition at line 22 of file cbc_plots_sim4.C.

TText* text = p_inj->AddText(inj_title)

Definition at line 1186 of file cbc_plots_sim4.C.

double time[6]

Definition at line 674 of file cbc_plots_sim4.C.

volume[m][l] = dV

Definition at line 1012 of file cbc_plots_sim4.C.

TString waveform[nfactor]

Definition at line 111 of file cbc_plots_sim4.C.

bool write_ascii = false

Definition at line 688 of file cbc_plots_sim4.C.

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

Definition at line 475 of file cbc_plots_sim4.C.

const Float_t* yq = new Float_t[101]

Definition at line 476 of file cbc_plots_sim4.C.