compare_bkg.C File Reference

Macros
#define	BINMAX   305.
#define	BINMIN   5.
#define	NBINS   3000
#define	REBIN_MIN_EVTS_PER_BIN   10.0

Functions
wave	Add (net_file_name)
live	Add (liv_file_name)
leg	AddEntry (curh, entry,"f")
leg	AddEntry (comphr, entry,"f")
curh	AndersonDarlingTest (comph,"D")
cout<<"After Rebinning"<< endl;cout <<"Anderson-Darling test (no rescaling)"<< endl;hnew2->	AndersonDarlingTest (hnew3,"D")
cout<<"Anderson-Darling 2-sample test significance distribution vs Normal (no rescaling)"<< endl;Sh4->	AndersonDarlingTest (Shn,"D")
cout<<"Plotting pad 1"<< endl;canvas1->	cd (1)
cout<<"Plotting pad 2"<< endl;canvas1->	cd (2)
cout<<"Plotting pad 3"<< endl;canvas1->	cd (3)
cout<<"Plotting pad 4"<< endl;canvas1->	cd (4)
cout<<"Plotting pad 5"<< endl;canvas1->	cd (5)
cout<<"Plotting pad 6"<< endl;canvas1->	cd (6)
canvas2	cd ()
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"))
cout<< "Setting up the canvas and the pads..."<< endl;TCanvas *canvas1=new TCanvas("canvas1","C1", 0, 0, 1200, 800);canvas1->	Divide (2, 3)
curh	Draw ()
comphr	Draw ("same")
hnew3	Draw ("SAME")
gr	Draw ("ALP")
	for (int i=1;i< nq+1;i++)
	for (int i=0;i<(int) Sh4->GetEntries();i++)
gROOT	ForceStyle ()
cout<<"Anderson-Darling 1-sample test significance distribution vs Normal (no rescaling) : "<< pvalueAD_1 << endl;cout<<"Starting Q-Q plot ..."<< endl;TGraphQQ *gr;if(entriesy > =entriesx){gr=new TGraphQQ(entriesy, rhoy, entriesx, rhox);gr->	GetHistogram () -> GetXaxis() ->SetTitle("Quantiles current histogram")
canvas1	GetPad (1) -> SetGridx()
canvas1	GetPad (2) -> SetLogx(kTRUE)
canvas1	GetPad (3) -> SetGridx()
canvas1	GetPad (4) -> SetGridx()
canvas1	GetPad (5) -> SetLogy(kTRUE)
canvas1	GetPad (6) -> SetLogy(kTRUE)
cout<<"nq ="<< nq<< endl;Double_t *xq=new Double_t[nq];Double_t *yq=new Double_t[nq];Double_t *yq2=new Double_t[nq+1];for(Int_t i=0;i< nq;i++) xq[i]=Float_t(i+1)/nq;curh->	GetQuantiles (nq, yq, xq)
comph	GetXaxis () -> SetRangeUser(BINMIN, BINMAX)
curh	GetYaxis () -> SetTitle("#events")
	if (strlen(veto_not_vetoed)==0)
	if (gSystem->Getenv("CWB_COMPARE_TREE")==NULL)
cout<<"Un-binned KS test" << endl;TMath::KolmogorovTest(entriesx, rhox, entriesy, rhoy,"D");cout <<"Binned KS test (no rescaling)"<< endl;curh->	KolmogorovTest (comph,"D")
TChain	live ("liveTime")
Compare_live	ReplaceAll ("wave_","live_")
comphr	Scale (K)
hnew3	Scale (liveTot/cliveTot)
gStyle	SetCanvasColor (kWhite)
curh	SetFillColor (kBlue)
comphr	SetFillColor (kRed)
comphr	SetFillStyle (3001)
gStyle	SetFrameBorderMode (0)
canvas2	SetGridx ()
canvas2	SetGridy ()
leg	SetHeader (entry)
comphr	SetLineColor (kRed)
line	SetLineColor (kGreen)
line	SetLineStyle (2)
line	SetLineWidth (2)
gr	SetMarkerColor (kRed)
gr	SetMarkerSize (1.0)
gr	SetMarkerStyle (20)
gStyle	SetNumberContours (256)
g1	SetParameter (0, count/sqrt(TMath::TwoPi()))
g1	SetParameter (1, 0.)
g1	SetParameter (2, 1.)
P	SetSeed (1234.)
gStyle	SetStatBorderSize (1)
leg	SetTextSize (0.03)
curh	SetTitle ("Distribution of rho")
gr	SetTitle ("Quantile-Quantile Plot")
gStyle	SetTitleFillColor (kWhite)
	sprintf (cnet_file_name,"%s", Compare_tree.Data())
	sprintf (cliv_file_name,"%s", Compare_live.Data())
	sprintf (st,"rho[%d]>>hcur(%d,%f,%f)", pp_irho, NBINS, BINMIN, BINMAX)
	sprintf (entry,"Rescaling factor K = %f", K)
	sprintf (entry,"Current BKG")
	sprintf (entry,"Rescaled Comparison BKG")
TChain	wave ("waveburst")

Variables
double	binmax = TMath::Min(gr2->GetX()[nq-1],gr2->GetY()[nq-1])
double	binmin = TMath::Max(gr2->GetX()[0],gr2->GetY()[0])
TCanvas *	canvas2 = new TCanvas("canvas2","C2",0,0,800,600)
TCanvas *	canvas3 = new TCanvas("canvas3","C3",0,0,800,600)
char	cliv_file_name [1024]
char	cnet_file_name [1024]
TString	Compare_live = Compare_tree
TString	Compare_tree
TH1D *	comphr =(TH1D*)comph->Clone("comphr")
int	count =0
	else
Long64_t	entriesx = cnet.Draw(st,mycut,"goff")
Long64_t	entriesy = wave.Draw(st,mycut,"goff")
char	entry [256]
TF1 *	g1 = new TF1("g1","gaus",-10.,10.)
cout<<"Comparison bkg events above threshold: "<< entriesx << endl;double *rhox=cnet.GetV1();comph=(TH1D *) gDirectory->	Get ("hcomp")
ROOT::Math::GoFTest *	goftest_1 = new ROOT::Math::GoFTest(count, w, ROOT::Math::GoFTest::kGaussian)
TGraphQQ *	gr2 = new TGraphQQ(nq,w,g1)
TH1D *	hnew2 = new TH1D("hnew","rebinned",nq,yq2)
TH1D *	hnew3 = new TH1D("hnew3","rebinned",nq,yq2)
TLegend *	leg = new TLegend(0.65,0.7,0.99,0.96)
TLine *	line = new TLine(BINMIN, BINMIN,BINMAX,BINMAX)
double	maxy = TMath::Max(comphr->GetMaximum(),curh->GetMaximum())
	mycut = T_cut
const Int_t	nq = entriesy/REBIN_MIN_EVTS_PER_BIN
TRandom3	P
Double_t	pvalueAD_1 = goftest_1-> AndersonDarlingTest()
TH1D *	Sh3 = new TH1D("Sh3","Significance vs rho bin",nq,yq2)
TH1D *	Sh4 = new TH1D("Sh4","Distribution of significance",200,-10.,10.)
TH1D *	Shn = new TH1D("Shn","Normal Distribution",200,-10.,10.)
double	w [nq+2]
double	w_tmp
	yq2 [0] = BINMIN

Macro Definition Documentation

#define BINMAX   305.

#define BINMIN   5.

#define NBINS   3000

#define REBIN_MIN_EVTS_PER_BIN   10.0

Function Documentation

wave Add

(

net_file_name

)

live Add

(

liv_file_name

)

leg AddEntry	(	curh	,
		entry	,
		"f"
	)

leg AddEntry	(	comphr	,
		entry	,
		"f"
	)

curh AndersonDarlingTest	(	comph	,
		"D"
	)

cout<<"After Rebinning"<<endl;cout <<"Anderson-Darling test (no rescaling)"<<endl;hnew2-> AndersonDarlingTest	(	hnew3	,
		"D"
	)

cout<<"Anderson-Darling 2-sample test significance distribution vs Normal (no rescaling)"<<endl;Sh4-> AndersonDarlingTest	(	Shn	,
		"D"
	)

cout<<"Plotting pad 1"<<endl;canvas1-> cd

(

1

)

cout<<"Plotting pad 2"<<endl;canvas1-> cd

(

2

)

cout<<"Plotting pad 3"<<endl;canvas1-> cd

(

3

)

cout<<"Plotting pad 4"<<endl;canvas1-> cd

(

4

)

cout<<"Plotting pad 5"<<endl;canvas1-> cd

(

5

)

cout<<"Plotting pad 6"<<endl;canvas1-> cd

(

6

)

canvas2 cd

(

)

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"

)

cout<< "Setting up the canvas and the pads..."<<endl;TCanvas *canvas1 = new TCanvas("canvas1","C1",0,0,1200,800);canvas1-> Divide	(	2	,
		3
	)

p_radius Draw

(

)

g2 Draw

(

"same"

)

hnew3 Draw

(

"SAME"

)

gr Draw

(

"ALP"

)

for

(

)

Definition at line 288 of file compare_bkg.C.

for

(

int

i = 0;i<(int)Sh4->GetEntries();i++

)

Definition at line 452 of file compare_bkg.C.

gROOT ForceStyle

(

)

cout<<"Anderson-Darling 1-sample test significance distribution vs Normal (no rescaling) : "<<pvalueAD_1<<endl;cout<<"Starting Q-Q plot ..."<<endl;TGraphQQ* gr;if (entriesy>=entriesx){ gr = new TGraphQQ(entriesy,rhoy,entriesx,rhox); gr-> GetHistogram

(

)

-> GetXaxis() ->SetTitle("Quantiles current histogram")

Definition at line 913 of file cbc_plots.C.

canvas1 GetPad

(

1

)

-> SetGridx()

canvas1 GetPad

(

2

)

-> SetLogx(kTRUE)

canvas1 GetPad

(

3

)

-> SetGridx()

canvas1 GetPad

(

4

)

-> SetGridx()

canvas1 GetPad

(

5

)

-> SetLogy(kTRUE)

canvas1 GetPad

(

6

)

-> SetLogy(kTRUE)

cout<<"nq ="<<nq<<endl;Double_t* xq = new Double_t[nq]; Double_t* yq = new Double_t[nq]; Double_t* yq2 = new Double_t[nq+1]; for (Int_t i=0;i<nq;i++) xq[i] = Float_t(i+1)/nq;curh-> GetQuantiles	(	nq	,
		yq	,
		xq
	)

comph GetXaxis

(

)

-> SetRangeUser(BINMIN, BINMAX)

Definition at line 321 of file cwb_report_prod_2.C.

curh GetYaxis

(

)

-> SetTitle("#events")

Definition at line 213 of file cbc_plots.C.

if

(

strlen(veto_not_vetoed)

= = 0

)

Definition at line 73 of file compare_bkg.C.

if

(

gSystem->

Getenv"CWB_COMPARE_TREE" = =NULL

)

Definition at line 100 of file compare_bkg.C.

cout<<"Un-binned KS test"<<endl;TMath::KolmogorovTest(entriesx,rhox,entriesy,rhoy,"D");cout<<"Binned KS test (no rescaling)"<<endl;curh-> KolmogorovTest	(	comph	,
		"D"
	)

TChain live

(

"liveTime"

)

Compare_live ReplaceAll	(	"wave_"	,
		"live_"
	)

comphr Scale

(

K

)

hnew3 Scale

(

liveTot/

cliveTot

)

gStyle SetCanvasColor

(

kWhite

)

Sh4 SetFillColor

(

kBlue

)

hnew3 SetFillColor

(

kRed

)

Sh4 SetFillStyle

(

3001

)

gStyle SetFrameBorderMode

(

0

)

canvas2 SetGridx

(

)

canvas2 SetGridy

(

)

leg SetHeader

(

entry

)

hnew3 SetLineColor

(

kRed

)

line SetLineColor

(

kGreen

)

line SetLineStyle

(

2

)

htemp3 SetLineWidth

(

2

)

gr2 SetMarkerColor

(

kRed

)

gr2 SetMarkerSize

(

1.

0

)

gr SetMarkerStyle

(

20

)

gStyle SetNumberContours

(

256

)

g1 SetParameter	(	0	,
		count/	sqrtTMath::TwoPi()
	)

g1 SetParameter	(	1	,
		0.
	)

g1 SetParameter	(	2	,
		1.
	)

P SetSeed

(

1234.

)

gStyle SetStatBorderSize

(

1

)

p_radius SetTextSize

(

0.

03

)

hnew2 SetTitle

(

"Distribution of rho"

)

gr2 SetTitle

(

"Quantile-Quantile Plot"

)

gStyle SetTitleFillColor

(

kWhite

)

sprintf	(	cnet_file_name	,
		"%s"	,
		Compare_tree.	Data()
	)

sprintf	(	cliv_file_name	,
		"%s"	,
		Compare_live.	Data()
	)

sprintf	(	st	,
		"rho>>hcur(%d,%f,%f)"	[%d],
		pp_irho	,
		NBINS	,
		BINMIN	,
		BINMAX
	)

sprintf

(

entry

)

sprintf	(	entry	,
		"Current BKG"
	)

sprintf	(	entry	,
		"Rescaled Comparison BKG"
	)

TChain wave

(

"waveburst"

)

Variable Documentation

double binmax = TMath::Min(gr2->GetX()[nq-1],gr2->GetY()[nq-1])

Definition at line 519 of file compare_bkg.C.

double binmin = TMath::Max(gr2->GetX()[0],gr2->GetY()[0])

Definition at line 518 of file compare_bkg.C.

TCanvas* canvas2 = new TCanvas("canvas2","C2",0,0,800,600)

Definition at line 57 of file compare_bkg.C.

TCanvas* canvas3 = new TCanvas("canvas3","C3",0,0,800,600)

Definition at line 64 of file compare_bkg.C.

char cliv_file_name[1024]

Definition at line 97 of file compare_bkg.C.

char cnet_file_name[1024]

Definition at line 96 of file compare_bkg.C.

TString Compare_live = Compare_tree

Definition at line 117 of file compare_bkg.C.

TString Compare_tree

Definition at line 99 of file compare_bkg.C.

TH1D* comphr =(TH1D*)comph->Clone("comphr")

Definition at line 234 of file compare_bkg.C.

int count =0

Definition at line 373 of file compare_bkg.C.

else

Initial value:

{
    Compare_tree=TString(gSystem->Getenv("CWB_COMPARE_TREE"))

Definition at line 102 of file compare_bkg.C.

Long64_t entriesx = cnet.Draw(st,mycut,"goff")

Definition at line 197 of file compare_bkg.C.

Long64_t entriesy = wave.Draw(st,mycut,"goff")

Definition at line 205 of file compare_bkg.C.

char entry[256]

Definition at line 250 of file compare_bkg.C.

TF1* g1 = new TF1("g1","gaus",-10.,10.)

Definition at line 505 of file compare_bkg.C.

cout<<"Current bkg events above threshold: "<< entriesy <<endl; double* rhoy = wave.GetV1(); double RHOMAX=TMath::Max(TMath::MaxElement(entriesy,rhoy),TMath::MaxElement(entriesx,rhox)); cout<<"Rho max: "<<RHOMAX<<endl; curh = (TH1D*)gDirectory-> Get("hcur")

Definition at line 201 of file compare_bkg.C.

ROOT::Math::GoFTest* goftest_1 = new ROOT::Math::GoFTest(count, w, ROOT::Math::GoFTest::kGaussian)

Definition at line 456 of file compare_bkg.C.

TGraphQQ* gr2 = new TGraphQQ(nq,w,g1)

Definition at line 510 of file compare_bkg.C.

TH1D* hnew2 = new TH1D("hnew","rebinned",nq,yq2)

Definition at line 291 of file compare_bkg.C.

TH1D* hnew3 = new TH1D("hnew3","rebinned",nq,yq2)

Definition at line 292 of file compare_bkg.C.

TLegend* leg = new TLegend(0.65,0.7,0.99,0.96)

Definition at line 246 of file compare_bkg.C.

TLine * line = new TLine(BINMIN, BINMIN,BINMAX,BINMAX)

Definition at line 482 of file compare_bkg.C.

double maxy = TMath::Max(comphr->GetMaximum(),curh->GetMaximum())

Definition at line 236 of file compare_bkg.C.

mycut = T_cut

Definition at line 187 of file compare_bkg.C.

const Int_t nq = entriesy/REBIN_MIN_EVTS_PER_BIN

Definition at line 281 of file compare_bkg.C.

TRandom3 P

Definition at line 296 of file compare_bkg.C.

Double_t pvalueAD_1 = goftest_1-> AndersonDarlingTest()

Definition at line 457 of file compare_bkg.C.

TH1D* Sh3 = new TH1D("Sh3","Significance vs rho bin",nq,yq2)

Definition at line 293 of file compare_bkg.C.

TH1D* Sh4 = new TH1D("Sh4","Distribution of significance",200,-10.,10.)

Definition at line 294 of file compare_bkg.C.

TH1D* Shn = new TH1D("Shn","Normal Distribution",200,-10.,10.)

Definition at line 295 of file compare_bkg.C.

double w[nq+2]

Definition at line 372 of file compare_bkg.C.

double w_tmp

Definition at line 371 of file compare_bkg.C.

yq2[nq] = BINMIN

Definition at line 289 of file compare_bkg.C.