anaMvdSim.C File Reference

Go to the source code of this file.

Functions
gROOT	LoadMacro ("../Tools.C")
	LoadPandaStyle ()
timer	Start ()
t	SetBranchAddress ("MVDPoint",&hit_array)
t	SetBranchAddress ("MCTrack",&mc_array)
	if (verbose)
can1	Divide (2, 2)
can1	cd (1)
	DrawNice2DHisto (hisxy)
can1	cd (2)
	DrawNice2DHisto (hisrz)
can1	cd (3)
gPad	SetLogy ()
hisde	DrawCopy ("")
can1	cd (4)
can1	Print ("outAnaMvdSim.ps")
timer	Stop ()

Variables
TStopwatch	timer
PndFileNameCreator	namecreator ("../data/Lars/MvdDtsSim.root")
std::string	inFile = namecreator.GetSimFileName(false)
TFile *	f = new TFile(inFile.c_str())
TTree *	t =(TTree *) f->Get("pndsim")
TClonesArray *	hit_array =new TClonesArray("PndSdsMCPoint")
TClonesArray *	mc_array =new TClonesArray("PndMCTrack")
TH2D *	hisxy = new TH2D("hisxy","MVD MC Points, xy view",100,-5.,5.,100,-5.,5.)
TH2D *	hisrz = new TH2D("hisrz","MVD MC Points, rz view",100,-5.,5.,100,-5.,5.)
TH1D *	hisde = new TH1D("hisde","MVD MC Points, Energyloss",100,0.,0.002)
int	nEvents = 1000
bool	verbose = false
TVector3	vecs
TVector3	veco
TCanvas *	can1 = new TCanvas("can1","MCHit view in MVD",0,0,800,800)
Double_t	rtime = timer.RealTime()
Double_t	ctime = timer.CpuTime()

Function Documentation

can1 cd

(

1

)

can1 cd

(

2

)

can1 cd

(

3

)

can1 cd

(

4

)

can1 Divide	(	2	,
		2
	)

hismom DrawCopy

(

""

)

DrawNice2DHisto

(

hisxy

)

DrawNice2DHisto

(

hisrz

)

if

(

verbose

)

Definition at line 38 of file anaMvdSim.C.

References hit, i, layer, and mcpdg.

    {
      if(verbose) cout<<"Point No "<<i<<endl;
      PndSdsMCPoint *hit=(PndSdsMCPoint*)hit_array->At(i);
      int mcpdg = -1;

      //PndMCTrack *mctruth = (PndMCTrack*)mc_array->At(hit->GetTrackID());
      //mcpdg = mctruth->GetPdgCode();
      //cout<<"mcpdg="<<mcpdg<<endl;

      vecs.SetXYZ(hit->GetX(), hit->GetY(), hit->GetZ());
      Int_t layer = Int_t(10.*vecs->Mag());
      if(verbose) cout<<vecs.x()<<"\t"<<vecs.y()<<"\t"<<vecs.z()<<endl;

      hisxy->Fill(vecs.x(),vecs.y());
      hisrz->Fill(vecs.z(),((vecs.y()>0.)?1.:-1.)*vecs.Perp());
      hisde->Fill(hit->GetEnergyLoss());
    }//end for i (points in event)

gROOT LoadMacro

(

"../Tools.C"

)

LoadPandaStyle

(

)

can1 Print

(

"outAnaMvdSim.ps"

)

t SetBranchAddress	(	"MVDPoint"	,
		&	hit_array
	)

t SetBranchAddress	(	"MCTrack"	,
		&	mc_array
	)

gPad SetLogy

(

)

timer Start

(

)

timer Stop

(

)

Variable Documentation

TCanvas* can1 = new TCanvas("can1","MCHit view in MVD",0,0,800,800)

Definition at line 59 of file anaMvdSim.C.

Double_t ctime = timer.CpuTime()

Definition at line 77 of file anaMvdSim.C.

TFile* f = new TFile(inFile.c_str())

Definition at line 17 of file anaMvdSim.C.

TH1D* hisde = new TH1D("hisde","MVD MC Points, Energyloss",100,0.,0.002)

Definition at line 27 of file anaMvdSim.C.

TH2D* hisrz = new TH2D("hisrz","MVD MC Points, rz view",100,-5.,5.,100,-5.,5.)

Definition at line 26 of file anaMvdSim.C.

TH2D* hisxy = new TH2D("hisxy","MVD MC Points, xy view",100,-5.,5.,100,-5.,5.)

Definition at line 25 of file anaMvdSim.C.

TClonesArray* hit_array =new TClonesArray("PndSdsMCPoint")

Definition at line 19 of file anaMvdSim.C.

std::string inFile = namecreator.GetSimFileName(false)

Definition at line 15 of file anaMvdSim.C.

TClonesArray* mc_array =new TClonesArray("PndMCTrack")

Definition at line 22 of file anaMvdSim.C.

PndFileNameCreator namecreator("../data/Lars/MvdDtsSim.root")

int nEvents = 1000

Definition at line 30 of file anaMvdSim.C.

Double_t rtime = timer.RealTime()

Definition at line 76 of file anaMvdSim.C.

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

Definition at line 18 of file anaMvdSim.C.

TStopwatch timer

Definition at line 10 of file anaMvdSim.C.

TVector3 veco

Definition at line 33 of file anaMvdSim.C.

Referenced by ana_MCOpt(), ana_MCpid(), and anaRadLength().

TVector3 vecs

Definition at line 33 of file anaMvdSim.C.

bool verbose = false

Definition at line 31 of file anaMvdSim.C.