piddraw.C File Reference

Go to the source code of this file.

Functions
gROOT	Macro ("$VMCWORKDIR/gconfig/rootlogon.C")
timer	Start ()
gStyle	SetHistFillColor (9)
gStyle	SetCanvasColor (0)
gStyle	SetLabelSize (0.025,"X")
gStyle	SetLabelSize (0.025,"Y")
gStyle	SetTitleSize (0.03,"Y")
gStyle	SetBarOffset (10)
gStyle	SetTitleOffset (1,"X")
gStyle	SetTitleOffset (1.55,"Y")
gStyle	SetTitleSize (0.03,"X")
gStyle	SetPalette (1, 0)
gStyle	SetOptFit (1111)
tree	SetBranchAddress ("MVDPoint",&pointlist)
hist10	SetYTitle ("(dE/dx)/(GeV/cm)")
hist10	SetXTitle ("p/(GeV/c)")
	for (int j=0;j< nEvents &&j< tree->GetEntriesFast();j++)
hist10	DrawCopy ("col")
timer	Stop ()

Variables
TStopwatch	timer
TFile *	inFile = new TFile("../data/MvdMC_Pid_10k.root","READ")
TTree *	tree = (TTree *)inFile->Get("pndsim")
TClonesArray *	pointlist =new TClonesArray("PndSdsMCPoint")
TH2D *	hist10 = new TH2D("specEnergyLossVsP","dE/dx(p)",200,0.0,1,200,0,0.0005)
TVector3	vecFront
TVector3	vecBack
TVector3	vecP
Int_t	nEvents = 5000
Double_t	dx
Double_t	dE
Double_t	p
Double_t	dEdX
TCanvas *	can1 = new TCanvas("test","Energy in MVD",20,20,800,600)
Double_t	rtime = timer.RealTime()
Double_t	ctime = timer.CpuTime()

Function Documentation

hist10 DrawCopy

(

"col"

)

for

(

)

Definition at line 40 of file piddraw.C.

References dx, PndSdsMCPoint::GetPxOut(), PndSdsMCPoint::GetPyOut(), PndSdsMCPoint::GetPzOut(), PndSdsMCPoint::GetXOut(), PndSdsMCPoint::GetYOut(), PndSdsMCPoint::GetZOut(), i, point, and vecFront.

        {
                tree->GetEntry(j);
                for(int i=0;i<pointlist->GetEntriesFast();i++)
                {
                  PndSdsMCPoint* point=(PndSdsMCPoint*)pointlist->At(i);

                  vecFront.SetXYZ(point->GetX(),point->GetY(),point->GetZ());
                  vecBack.SetXYZ(point->GetXOut(),point->GetYOut(),point->GetZOut());
                  vecP.SetXYZ(point->GetPxOut(),point->GetPyOut(),point->GetPzOut());

        dx=(vecBack-vecFront).Mag();
      dE=point->GetEnergyLoss();
      p=vecP.Mag();

      if(dx!=0)
      {
        dEdX=(dE/dx);
        hist10->Fill(p,dE);
      }
    }
  }

gROOT Macro

(

"$VMCWORKDIR/gconfig/rootlogon.C"

)

gStyle SetBarOffset

(

10

)

tree SetBranchAddress	(	"MVDPoint"	,
		&	pointlist
	)

gStyle SetCanvasColor

(

0

)

gStyle SetHistFillColor

(

9

)

gStyle SetLabelSize	(	0.	025,
		"X"
	)

gStyle SetLabelSize	(	0.	025,
		"Y"
	)

gStyle SetOptFit

(

1111

)

gStyle SetPalette	(	1	,
		0
	)

gStyle SetTitleOffset	(	1	,
		"X"
	)

gStyle SetTitleOffset	(	1.	55,
		"Y"
	)

gStyle SetTitleSize	(	0.	03,
		"Y"
	)

gStyle SetTitleSize	(	0.	03,
		"X"
	)

hist10 SetXTitle

(

"p/(GeV/c)"

)

hist10 SetYTitle

(

"(dE/dx)/(GeV/cm)"

)

timer Start

(

)

timer Stop

(

)

Variable Documentation

TCanvas* can1 = new TCanvas("test","Energy in MVD",20,20,800,600)

Definition at line 63 of file piddraw.C.

Double_t ctime = timer.CpuTime()

Definition at line 69 of file piddraw.C.

Double_t dE

Definition at line 37 of file piddraw.C.

Double_t dEdX

Definition at line 37 of file piddraw.C.

Double_t dx

Definition at line 37 of file piddraw.C.

Referenced by for().

TH2D* hist10 = new TH2D("specEnergyLossVsP","dE/dx(p)",200,0.0,1,200,0,0.0005)

Definition at line 30 of file piddraw.C.

TFile* inFile = new TFile("../data/MvdMC_Pid_10k.root","READ")

Definition at line 24 of file piddraw.C.

Int_t nEvents = 5000

Definition at line 36 of file piddraw.C.

Double_t p

Definition at line 37 of file piddraw.C.

TClonesArray* pointlist =new TClonesArray("PndSdsMCPoint")

Definition at line 26 of file piddraw.C.

Double_t rtime = timer.RealTime()

Definition at line 68 of file piddraw.C.

TStopwatch timer

Definition at line 7 of file piddraw.C.

TTree* tree = (TTree *)inFile->Get("pndsim")

Definition at line 25 of file piddraw.C.

TVector3 vecBack

Definition at line 35 of file piddraw.C.

TVector3 vecFront

Definition at line 35 of file piddraw.C.

Referenced by for().

TVector3 vecP

Definition at line 35 of file piddraw.C.