testChi2Ordering.C File Reference

Functions
h_chi2_v_ttheta	Draw ("colz")
h_chi2	Draw ()
	for (int j=0;j< 1;j++)
ali	Print ()
new	TCanvas ()

Variables
EdbPVRec *	ali = (EdbPVRec*)f->Get("EdbPVRec")
TObjArray *	arr_s = new TObjArray()
Float_t	chi2Normalized =0
TH1F *	h_chi2 = new TH1F("h_chi2","h_chi2",100,0,2)
TH2F *	h_chi2_v_ttheta = new TH2F("h_chi2_v_ttheta","h_chi2_v_ttheta",100,0,1,100,0,2)
TH1F *	h_chi2Normalized = new TH1F("h_chi2Normalized","h_chi2Normalized",100,0,5)
TH1F *	h_WTilde = new TH1F("h_WTilde","h_WTilde",40,0,0.1)
TH2F *	h_WTilde_v_ttheta = new TH2F("h_WTilde_v_ttheta","h_WTilde_v_ttheta",100,0,1,40,0,0.1)
TH1F *	h_WTildeNormalized = new TH1F("h_WTildeNormalized","h_WTildeNormalized",100,0,5)
TH1F *	h_X2 = new TH1F("h_X2","h_X2",100,0,5)
TProfile *	po = new TProfile("name","title",20,0,1)
TProfile *	po2 = new TProfile("name2","title2",20,0,1)
	return
Float_t	TTheta =0
Float_t	WTilde =0
Float_t	WTildeNormalized =0
Float_t	X2 =0

Function Documentation

Draw() [1/2]

h_chi2_v_ttheta Draw

(

"colz"

)

Draw() [2/2]

h_chi2 Draw

(

)

for()

for

(

)

                       {
pat = ali->GetPattern(j);
//----------------------------------------------
for (int i=0; i<30000;i++) {
EdbSegP* s = (EdbSegP*)pat->GetSegment(i);
TTheta = TMath::Sqrt(s->TX()*s->TX()+s->TY()*s->TY());
// cout << abs(TTheta-0.05) << endl;
// if (TMath::Abs(TTheta-0.2)>0.05) continue;
 
h_chi2_v_ttheta->Fill(TMath::Sqrt(s->TX()*s->TX()+s->TY()*s->TY()),s->Chi2());
h_chi2->Fill(s->Chi2());
WTilde=1./(s->W());
h_WTilde_v_ttheta->Fill(TMath::Sqrt(s->TX()*s->TX()+s->TY()*s->TY()),WTilde);
h_WTilde->Fill(WTilde);
 
 
} // i
// // // // // // // // // // // // // // // // // // // // // 
 
Float_t m_chi2=h_chi2->GetMean();
Float_t s_chi2=h_chi2->GetRMS();
Float_t m_WTilde=h_WTilde->GetMean();
Float_t s_WTilde=h_WTilde->GetRMS();
 
for (int i=0; i<30000;i++) {
EdbSegP* s = (EdbSegP*)pat->GetSegment(i);
TTheta = TMath::Sqrt(s->TX()*s->TX()+s->TY()*s->TY());
 
chi2Normalized=TMath::Power((s->Chi2()-m_chi2)/s_chi2,2);
h_chi2Normalized->Fill(chi2Normalized);
WTilde=1./(s->W());
WTildeNormalized=TMath::Power((WTilde-m_WTilde)/s_WTilde,2);
h_WTildeNormalized->Fill(WTildeNormalized);
// s->PrintNice();
// cout << TMath::Power((WTilde-mWTilde)/s_WTilde,2) << endl;
X2=chi2Normalized + WTildeNormalized;
h_X2->Fill(X2);
} // i
// // // // // // // // // // // // // // // // // // // // // // // // // // // // 
 
 
//----------------------------------------------
/*
// // Loop over bin entries of the profile histogram:
for (int jj=1; jj<2;jj++) {
Float_t mu= po->GetBinContent(jj);
Float_t sigma= po->GetBinError(jj);
for (int i=0; i<30000;i++) {
EdbSegP* s = (EdbSegP*)pat->GetSegment(i);
if (s->Chi2() > mu+1*sigma) continue;
po2->Fill(TMath::Sqrt(s->TX()*s->TX()+s->TY()*s->TY()),s->Chi2());
} // i
}
*/
 
 
 
} // j

Print()

ali Print

(

)

TCanvas()

new TCanvas

(

)

Variable Documentation

ali

EdbPVRec* ali = (EdbPVRec*)f->Get("EdbPVRec")

arr_s

TObjArray* arr_s = new TObjArray()

chi2Normalized

Float_t chi2Normalized =0

h_chi2

TH1F* h_chi2 = new TH1F("h_chi2","h_chi2",100,0,2)

h_chi2_v_ttheta

TH2F* h_chi2_v_ttheta = new TH2F("h_chi2_v_ttheta","h_chi2_v_ttheta",100,0,1,100,0,2)

h_chi2Normalized

TH1F* h_chi2Normalized = new TH1F("h_chi2Normalized","h_chi2Normalized",100,0,5)

h_WTilde

TH1F* h_WTilde = new TH1F("h_WTilde","h_WTilde",40,0,0.1)

h_WTilde_v_ttheta

TH2F* h_WTilde_v_ttheta = new TH2F("h_WTilde_v_ttheta","h_WTilde_v_ttheta",100,0,1,40,0,0.1)

h_WTildeNormalized

TH1F* h_WTildeNormalized = new TH1F("h_WTildeNormalized","h_WTildeNormalized",100,0,5)

h_X2

TH1F* h_X2 = new TH1F("h_X2","h_X2",100,0,5)

po

TProfile* po = new TProfile("name","title",20,0,1)

po2

TProfile* po2 = new TProfile("name2","title2",20,0,1)

return

return

TTheta

Float_t TTheta =0

WTilde

Float_t WTilde =0

WTildeNormalized

Float_t WTildeNormalized =0

X2

Float_t X2 =0