runPULL2.C File Reference

Go to the source code of this file.

Functions
gStyle	SetOptFit (1)
gROOT	LoadMacro ("$VMCWORKDIR/gconfig/basiclibs.C")
	basiclibs ()
gSystem	Load ("libGeoBase")
gSystem	Load ("libParBase")
gSystem	Load ("libBase")
gSystem	Load ("libField")
gSystem	Load ("libPndData")
gSystem	Load ("libGeane")
gSystem	Load ("libPlane")
gSystem	Load ("libGeaneEx")
gSystem	Load ("libTrkBase")
simtree	SetBranchAddress ("GeaneTrackFinal",&fTrackParFinal)
simtree	SetBranchAddress ("GeaneTrackPar",&fTrackParGeane)
	for (Int_t i=0;i< Nevents;i++)
c	Divide (5, 2)
c	cd (1)
h1	Draw ()
h1	Fit ("gaus")
c	cd (2)
c	cd (3)
c	cd (4)
c	cd (5)
c	cd (6)
c	cd (7)
c	cd (8)
c	cd (9)
c	cd (10)
c	cd ()

Variables
TFile *	f =new TFile("ex2g.root")
TTree *	simtree =f->Get("pndsim")
	fTrackParGeane = new TClonesArray("FairTrackParP")
	fTrackParIni = new TClonesArray("FairTrackParP")
	fTrackParFinal = new TClonesArray("FairTrackParP")
TH1F *	h1 =new TH1F("h1","V ",100,-10.,10.)
TH1F *	h2 =new TH1F("h2","W",100,-10,10)
TH1F *	h3 =new TH1F("h3","TV",100,-10,10)
TH1F *	h4 =new TH1F("h4","TW",100,-10,10)
TH1F *	h5 =new TH1F("h5","Qp",100,-10,10)
TH1F *	h6 =new TH1F("h6","Y",100, -10,10)
TH1F *	h7 =new TH1F("h7","Z",100,-10,10)
TH1F *	h8 =new TH1F("h8","Px",100,-10,10)
TH1F *	h9 =new TH1F("h9","Py",100,-10,10)
TH1F *	h10 =new TH1F("h10","Pz",100,-10,10)
FairTrackParP *	fTrkF
FairTrackParP *	fTrkG
Int_t	Nevents = simtree->GetEntriesFast()
TCanvas *	c =new TCanvas("c","c",900,600)

Function Documentation

basiclibs

(

)

c cd

(

1

)

c cd

(

2

)

c cd

(

3

)

c cd

(

4

)

c cd

(

5

)

c cd

(

6

)

c cd

(

7

)

c cd

(

8

)

c cd

(

9

)

c cd

(

10

)

c cd

(

)

c Divide	(	5	,
		2
	)

h1 Draw

(

)

Definition at line 64 of file NHitsPerEvent.C.

           {

        std::cout << "mvdHitArray: " << mvdHitArray->GetEntries() << " projectedHitArray " << mvdProjectedArray->GetEntries() << std::endl;

//      gHitPoints->Clear();
        ExtractData(mvdHitArray, gHitPoints);
        ExtractData(mvdProjectedArray, gProjections);
        //if (gHitPoints->GetN() != x.size()) gHitPoints->Expand(x.size());

        std::cout << "NHits: " << gHitPoints->GetN() << std::endl;
        gHitPoints->Draw("apl");
        gProjections->Draw("pl");
}

h1 Fit

(

"gaus"

)

for

(

)

Definition at line 44 of file runPULL2.C.

References cemc, cluster_energy, digi, digi_ene, PndEmcCluster::DigiList(), Double_t, PndEmcCluster::energy(), fTrackParFinal, fTrackParGeane, PndMCTrack::Get4Momentum(), PndEmcHit::GetCopy(), PndEmcDigi::GetCopy(), PndEmcHit::GetCrystal(), PndEmcDigi::GetCrystal(), PndEmcDigi::GetDetectorId(), PndEmcHit::GetEnergy(), PndEmcDigi::GetEnergy(), PndEmcHit::GetModule(), PndEmcDigi::GetModule(), PndMCTrack::GetMomentum(), PndEmcDigi::GetPhi(), PndEmcDigi::GetPhiInt(), PndEmcHit::GetRow(), PndEmcDigi::GetRow(), PndEmcDigi::GetTheta(), PndEmcDigi::GetThetaInt(), PndEmcHit::GetXPad(), PndEmcDigi::GetXPad(), PndEmcHit::GetYPad(), PndEmcDigi::GetYPad(), hit(), hit_ene, i, mctrack, nhits, Pi, PndEmcCluster::where(), PndEmcDigi::where(), and PndEmcCluster::z().

                                      {
          fTrackParGeane->Delete();
          fTrackParFinal->Delete();
          simtree->GetEntry(i);
          for (Int_t k=0; k<fTrackParGeane->GetEntriesFast(); k++)      {
            fTrkF = (FairTrackParP *)fTrackParFinal->At(k);
            fTrkG = (FairTrackParP *)fTrackParGeane->At(k);
            if(!fTrkF || !fTrkG) continue;
            if(fTrkF &&fTrkG ){
    if(fTrkG->GetDV())  h1->Fill((fTrkF->GetV()  - fTrkG->GetV())/fTrkG->GetDV());
              if(fTrkG->GetDW())  h2->Fill((fTrkF->GetW()  - fTrkG->GetW())/fTrkG->GetDW());
              if(fTrkG->GetDTV()) h3->Fill((fTrkF->GetTV() - fTrkG->GetTV())/fTrkG->GetDTV());
              if(fTrkG->GetDTW()) h4->Fill((fTrkF->GetTW() - fTrkG->GetTW())/fTrkG->GetDTW());
              if(fTrkG->GetDQp())h5->Fill((fTrkF->GetQp() - fTrkG->GetQp()) /fTrkG->GetDQp());
              if(fTrkG->GetDY()) h6->Fill((fTrkF->GetY()  - fTrkG->GetY()) /fTrkG->GetDY());
              if(fTrkG->GetDZ()) h7->Fill((fTrkF->GetZ()  - fTrkG->GetZ()) /fTrkG->GetDZ());
              if(fTrkG->GetDPx()) h8->Fill((fTrkF->GetPx() - fTrkG->GetPx())/fTrkG->GetDPx());
              if(fTrkG->GetDPy()) h9->Fill((fTrkF->GetPy() - fTrkG->GetPy())/fTrkG->GetDPy());
              if(fTrkG->GetDPz()) h10->Fill((fTrkF->GetPz()- fTrkG->GetPz())/fTrkG->GetDPz());
            }
          }
        }

gSystem Load

(

"libGeoBase"

)

gSystem Load

(

"libParBase"

)

gSystem Load

(

"libBase"

)

gSystem Load

(

"libField"

)

gSystem Load

(

"libPndData"

)

gSystem Load

(

"libGeane"

)

gSystem Load

(

"libPlane"

)

gSystem Load

(

"libGeaneEx"

)

gSystem Load

(

"libTrkBase"

)

gROOT LoadMacro

(

"$VMCWORKDIR/gconfig/basiclibs.C"

)

simtree SetBranchAddress	(	"GeaneTrackFinal"	,
		&	fTrackParFinal
	)

simtree SetBranchAddress	(	"GeaneTrackPar"	,
		&	fTrackParGeane
	)

gStyle SetOptFit

(

1

)

Variable Documentation

TCanvas* c =new TCanvas("c","c",900,600)

Definition at line 66 of file runPULL2.C.

TFile* f =new TFile("ex2g.root")

Definition at line 20 of file runPULL2.C.

fTrackParFinal = new TClonesArray("FairTrackParP")

Definition at line 25 of file runPULL2.C.

Referenced by for().

fTrackParGeane = new TClonesArray("FairTrackParP")

Definition at line 23 of file runPULL2.C.

Referenced by for().

fTrackParIni = new TClonesArray("FairTrackParP")

Definition at line 24 of file runPULL2.C.

FairTrackParP* fTrkF

Definition at line 40 of file runPULL2.C.

FairTrackParP* fTrkG

Definition at line 41 of file runPULL2.C.

TH1F* h1 =new TH1F("h1","V ",100,-10.,10.)

Definition at line 27 of file runPULL2.C.

TH1F* h10 =new TH1F("h10","Pz",100,-10,10)

Definition at line 36 of file runPULL2.C.

TH1F* h2 =new TH1F("h2","W",100,-10,10)

Definition at line 28 of file runPULL2.C.

TH1F* h3 =new TH1F("h3","TV",100,-10,10)

Definition at line 29 of file runPULL2.C.

TH1F* h4 =new TH1F("h4","TW",100,-10,10)

Definition at line 30 of file runPULL2.C.

TH1F* h5 =new TH1F("h5","Qp",100,-10,10)

Definition at line 31 of file runPULL2.C.

TH1F* h6 =new TH1F("h6","Y",100, -10,10)

Definition at line 32 of file runPULL2.C.

TH1F* h7 =new TH1F("h7","Z",100,-10,10)

Definition at line 33 of file runPULL2.C.

TH1F* h8 =new TH1F("h8","Px",100,-10,10)

Definition at line 34 of file runPULL2.C.

TH1F* h9 =new TH1F("h9","Py",100,-10,10)

Definition at line 35 of file runPULL2.C.

Int_t Nevents = simtree->GetEntriesFast()

Definition at line 43 of file runPULL2.C.

TTree* simtree =f->Get("pndsim")

Definition at line 21 of file runPULL2.C.