Functions
gROOT	Macro ("$VMCWORKDIR/gconfig/rootlogon.C")

timer	Start ()

t	SetBranchAddress ("LMDHitsStrip",&hit_array)

T	SetBranchAddress ("LMDPoint",&mc_array)

h2	SetFillColor (kRed)

hisxy	SetTitle ("LMD Reconstructed points, xy view;x / cm;y / cm")

hisrz	SetTitle ("LMD Reconstructed points, rz view;z / cm;r/ cm")

mchisxy	SetTitle ("LMD MC points, xy view;x / cm;y / cm")

mchisrz	SetTitle ("LMD MC points, rz view;z / cm;r/ cm")

hisDiffXY	SetTitle ("MC - RECO Hit coordinates xy view;#Deltax / cm;#Deltay / cm")

hisDiffRZ	SetTitle ("MC - RECO Hit coordinates rz view;#Deltaz / cm;#Deltar / cm")

hisDiffX	SetTitle ("MC - RECO Hit coordinate x;x / cm;")

hisDiffY	SetTitle ("MC - RECO Hit coordinate y;y / cm;")

hisDiffZ	SetTitle ("MC - RECO Hit coordinate z;z / cm;")

T	GetEntry (j)

Variables
TStopwatch	timer

TFile *	f = new TFile("/private/huagen/simdata/Lmd_Reco_DPM_elastic_6.2_1.9mrad_5M_2.root")

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

TClonesArray *	hit_array =new TClonesArray("PndSdsHit")

TFile *	F = new TFile("/private/huagen/simdata/Lmd_DPM_elastic_6.2_1.9mrad_5M_2.root")

TTree *	T =(TTree *) F->Get("pndsim")

TClonesArray *	mc_array =new TClonesArray("PndSdsMCPoint")

TH3F *	h2 = new TH3F("c2","Reconstructed Points", 200, -10, 10, 100, -10, 10, 100, 1070,1150)

TH3F *	h3 = new TH3F("c3","MC Points", 200, -10, 10, 100, -10, 10, 100, 1070,1150)

TH2D *	hisxy = new TH2D("hisxy","",200,-10.,10.,200,-10.,10.)

TH2D *	hisrz = new TH2D("hisrz","",800,1070.,1150.,100,-10.,10.)

TH2D *	mchisxy = new TH2D("mchisxy","",200,-10.,10.,200,-10.,10.)

TH2D *	mchisrz = new TH2D("mchisrz","",800,1070.,1150.,100,-10.,10.)

TH2D *	hisDiffXY = new TH2D("hisdiffxy","",100,-0.05,0.05,100,-0.05,0.05)

TH2D *	hisDiffRZ = new TH2D("hisdiffrz","",100,-0.1.,0.1.,100,-0.05,0.05)

TH1F *	hisDiffX = new TH1F("hisdiffx","",400,-0.02.,0.02.)

TH1F *	hisDiffY = new TH1F("hisdiffy","",400,-0.02.,0.02)

TH1F *	hisDiffZ = new TH1F("hisdiffz","",400,-0.02,0.02)

TH1D *	RecoCharge = new TH1D("RecoCharge","",200,0.,1000)

TH1D *	DigiCharge = new TH1D("DigiCharge","",200,0.,1000)

TH1D *	MCCharge = new TH1D("MCCharge","",200,0.,1000)

TVector3	vecs

TVector3	vecmc

TVector3	vecdiif

double	mcX =0

double	mcY =0

double	mcZ =0

double	reX =0

double	reY =0

double	reZ =0

double	diffX =0

double	diffY =0

double	diffZ =0

TVector3	RecoPosition

TVector3	MCposition

TVector3	vecdiff

cout<<"************Event No********** " << j<< endl;{if(hit_array-> GetEntriesFast()!=mc_array-> GetEntriesFast()) continue;PndSdsHit hit=(PndSdsHit ) hit_array->	At (i)

PndSdsMCPoint *	point =(PndSdsMCPoint*)mc_array->At(i)

Int_t	ID_hit = hit.GetSensorID()

Int_t	ID_mc = point.GetSensorID()

Function Documentation

T GetEntry ( j )

gROOT Macro ( "$VMCWORKDIR/gconfig/rootlogon.C" )

t SetBranchAddress	(	"LMDHitsStrip"	,
		&	hit_array
	)

T SetBranchAddress	(	"LMDPoint"	,
		&	mc_array
	)

h3 SetFillColor ( kRed )

hisxy SetTitle	(	"LMD Reconstructed	points,
		xy view;x/cm;y/cm"
	)

hisrz SetTitle	(	"LMD Reconstructed	points,
		rz view;z/cm;r/cm"
	)

mchisxy SetTitle	(	"LMD MC	points,
		xy view;x/cm;y/cm"
	)

mchisrz SetTitle	(	"LMD MC	points,
		rz view;z/cm;r/cm"
	)

hisDiffXY SetTitle ( "MC - RECO Hit coordinates xy view;#Deltax / cm;#Deltay / cm" )

hisDiffRZ SetTitle ( "MC - RECO Hit coordinates rz view;#Deltaz / cm;#Deltar / cm" )

hisDiffX SetTitle ( "MC - RECO Hit coordinate x;x / cm;" )

hisDiffY SetTitle ( "MC - RECO Hit coordinate y;y / cm;" )

hisDiffZ SetTitle ( "MC - RECO Hit coordinate z;z / cm;" )

timer Start ( )

Variable Documentation

cout<<"**************Event No************ "<<j<<endl; { if(hit_array->GetEntriesFast()!= mc_array->GetEntriesFast())continue; PndSdsHit *hit=(PndSdsHit*)hit_array-> At(i)

Definition at line 93 of file anaLmdReco.C.

double diffX =0

Definition at line 76 of file anaLmdReco.C.

Referenced by PndLmdSigCleanTask::CheckXY(), and PndLmdTrksFilterTask::Exec().

double diffY =0

Definition at line 76 of file anaLmdReco.C.

Referenced by PndLmdSigCleanTask::CheckXY(), and PndLmdTrksFilterTask::Exec().

double diffZ =0

Definition at line 76 of file anaLmdReco.C.

TH1D* DigiCharge = new TH1D("DigiCharge","",200,0.,1000)

Definition at line 70 of file anaLmdReco.C.

TFile* f = new TFile("/private/huagen/simdata/Lmd_Reco_DPM_elastic_6.2_1.9mrad_5M_2.root")

Definition at line 24 of file anaLmdReco.C.

TFile* F = new TFile("/private/huagen/simdata/Lmd_DPM_elastic_6.2_1.9mrad_5M_2.root")

Definition at line 31 of file anaLmdReco.C.

Referenced by anaGemSmearing().

TH3F* h2 = new TH3F("c2","Reconstructed Points", 200, -10, 10, 100, -10, 10, 100, 1070,1150)

Definition at line 41 of file anaLmdReco.C.

TH3F* h3 = new TH3F("c3","MC Points", 200, -10, 10, 100, -10, 10, 100, 1070,1150)

Definition at line 43 of file anaLmdReco.C.

TH2D* hisDiffRZ = new TH2D("hisdiffrz","",100,-0.1.,0.1.,100,-0.05,0.05)

Definition at line 59 of file anaLmdReco.C.

Referenced by anaclust(), and anaGemSmearing().

TH1F* hisDiffX = new TH1F("hisdiffx","",400,-0.02.,0.02.)

Definition at line 62 of file anaLmdReco.C.

Referenced by anaclust(), and anaGemSmearing().

TH2D* hisDiffXY = new TH2D("hisdiffxy","",100,-0.05,0.05,100,-0.05,0.05)

Definition at line 57 of file anaLmdReco.C.

Referenced by anaclust(), and anaGemSmearing().

TH1F* hisDiffY = new TH1F("hisdiffy","",400,-0.02.,0.02)

Definition at line 64 of file anaLmdReco.C.

Referenced by anaclust(), and anaGemSmearing().

TH1F* hisDiffZ = new TH1F("hisdiffz","",400,-0.02,0.02)

Definition at line 66 of file anaLmdReco.C.

Referenced by anaclust(), and anaGemSmearing().

TH2D* hisrz = new TH2D("hisrz","",800,1070.,1150.,100,-10.,10.)

Definition at line 48 of file anaLmdReco.C.

TH2D* hisxy = new TH2D("hisxy","",200,-10.,10.,200,-10.,10.)

Definition at line 46 of file anaLmdReco.C.

TClonesArray* hit_array =new TClonesArray("PndSdsHit")

Definition at line 27 of file anaLmdReco.C.

Int_t ID_hit = hit.GetSensorID()

Definition at line 95 of file anaLmdReco.C.

Int_t ID_mc = point.GetSensorID()

Definition at line 96 of file anaLmdReco.C.

TClonesArray* mc_array =new TClonesArray("PndSdsMCPoint")

Definition at line 33 of file anaLmdReco.C.

TH1D* MCCharge = new TH1D("MCCharge","",200,0.,1000)

Definition at line 71 of file anaLmdReco.C.

TH2D* mchisrz = new TH2D("mchisrz","",800,1070.,1150.,100,-10.,10.)

Definition at line 54 of file anaLmdReco.C.

TH2D* mchisxy = new TH2D("mchisxy","",200,-10.,10.,200,-10.,10.)

Definition at line 52 of file anaLmdReco.C.

TVector3 MCposition

Definition at line 77 of file anaLmdReco.C.

Referenced by PndTrkTracking2::Exec().

double mcX =0

Definition at line 74 of file anaLmdReco.C.

double mcY =0

Definition at line 74 of file anaLmdReco.C.

double mcZ =0

Definition at line 74 of file anaLmdReco.C.

PndSdsMCPoint* point =(PndSdsMCPoint*)mc_array->At(i)

Definition at line 94 of file anaLmdReco.C.

TH1D* RecoCharge = new TH1D("RecoCharge","",200,0.,1000)

Definition at line 69 of file anaLmdReco.C.

TVector3 RecoPosition

Definition at line 77 of file anaLmdReco.C.

double reX =0

Definition at line 75 of file anaLmdReco.C.

double reY =0

Definition at line 75 of file anaLmdReco.C.

double reZ =0

Definition at line 75 of file anaLmdReco.C.

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

Definition at line 26 of file anaLmdReco.C.

TTree* T =(TTree *) F->Get("pndsim")

Definition at line 32 of file anaLmdReco.C.

Referenced by TTracksCatCounters< int >::AddCounter(), TTracksCatCounters< int >::AddCounters(), PndCAInternal::Allocator< T, alignment >::Alloc(), PndFTSInternal::Allocator< T, alignment >::Alloc(), PndCAInternal::Allocator< T, 0 >::Alloc(), PndFTSInternal::Allocator< T, 0 >::Alloc(), nsL1::SimdAlloc< T >::allocate(), anaGemSmearing(), Vec< T >::at(), BinaryStoreRead(), BinaryStoreWrite(), RhoKalmanVtxFitter::Calculate(), PndFTSCAParameters::CALocalToGlobal(), nsL1::SimdAlloc< T >::construct(), PndFTSInternal::AlignedData< typename PndFTSInternal::TypeForAlignmentHelper< T, alignment >::Type, Size::Size, alignment >::ConstructAlignedData(), PndCAInternal::AlignedData< typename PndCAInternal::TypeForAlignmentHelper< T, alignment >::Type, Size::Size, alignment >::ConstructAlignedData(), nsL1::SimdAlloc< T >::deallocate(), RhoLorentzVectorErr::DetermineChisq(), GFMaterialEffects::energyLossBrems(), PndCAParameters::GlobalToCALocal(), PndFTSCAParameters::GlobalToCALocal(), Photospp::PhotosRandom::initialize(), nsL1::SimdAlloc< T >::max_size(), FTSCAElementsOnStation< T >::operator+(), Vec< T >::operator[](), PndCAFixedArray< T, Size, alignment >::PndCAFixedArray(), PndFTSFixedArray< T, Size, alignment >::PndFTSFixedArray(), PndSttTrackFinderReal::PndSttTrkFindCircles(), RhoLorentzVectorErr::PrintOn(), CAMath::Recip(), PndFTSArray< PndFTSInternal::TypeForAlignmentHelper< PndFTSCAGBHit, 0 >::Type, 1 >::ReinterpretCast(), PndCAArray< PndCAInternal::TypeForAlignmentHelper< T, alignment >::Type, Size::Dim >::ReinterpretCast(), CAMath::RSqrt(), sign(), TTracksCatCounters< int >::TTracksCatCounters(), PndFTSInternal::AlignedData< typename PndFTSInternal::TypeForAlignmentHelper< T, alignment >::Type, Size::Size, alignment >::~AlignedData(), and PndCAInternal::AlignedData< typename PndCAInternal::TypeForAlignmentHelper< T, alignment >::Type, Size::Size, alignment >::~AlignedData().

TStopwatch timer

Definition at line 21 of file anaLmdReco.C.

TVector3 vecdiff

Definition at line 77 of file anaLmdReco.C.

Referenced by anaclust(), anaGemSmearing(), PndAnalysis::Propagator(), and QAmacro_mvd_ana().

TVector3 vecdiif

Definition at line 73 of file anaLmdReco.C.

Referenced by anaGemSmearing().

TVector3 vecmc

Definition at line 73 of file anaLmdReco.C.

TVector3 vecs

Definition at line 73 of file anaLmdReco.C.

Functions

Variables

Function Documentation

Variable Documentation