Functions
gROOT	LoadMacro ("$VMCWORKDIR/gconfig/basiclibs.C")

	basiclibs ()

gSystem	Load ("libGeoBase")

gSystem	Load ("libParBase")

gSystem	Load ("libBase")

gSystem	Load ("libPndData")

gSystem	Load ("libField")

gSystem	Load ("libGen")

gSystem	Load ("libPassive")

gSystem	Load ("libgenfit")

gSystem	Load ("libtrackrep")

gSystem	Load ("libtpc")

gSystem	Load ("libtpcreco")

gSystem	Load ("librecotasks")

gSystem	Load ("libMvd")

gSystem	Load ("libMvdReco")

timer	Start ()

t	SetBranchAddress ("MVDCluster",&hit_array)

T	SetBranchAddress ("MVDPoint",&mc_array)

hisxy	SetTitle ("MVD MC Cluster, xy view;x / cm;y / cm")

hisrz	SetTitle ("MVD MC Cluster, 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;#Deltax / cm;")

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

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

hisDiffLocalX	SetTitle ("(MC - RECO) Local Hit coordinate x;#Deltax_{L} / cm;")

hisDiffLocalY	SetTitle ("(MC - RECO) Local Hit coordinate y;#Deltay_{L} / cm;")

hisDiffLocalZ	SetTitle ("(MC - RECO) Local Hit coordinate z;#Deltaz_{L} / cm;")

hisDiffLocalXY	SetTitle ("(MC - RECO) Local Hit XY;#Deltax_{L} / cm;#Deltay_{L} / cm")

hisDiffLocalRZ	SetTitle ("(MC - RECO) Local Hit RZ;#Deltaz_{L} / cm;#Deltar_{L} / cm")

hisLocalXY	SetTitle ("(MC - RECO) Local Hit XY;x_{L} / cm;y_{L} / cm")

hisLocalRZ	SetTitle ("(MC - RECO) Local Hit RZ;z_{L} / cm;r_{L} / cm")

hisLocalYZ	SetTitle ("(MC - RECO) Local Hit YZ;z_{L} / cm;y_{L} / cm")

hisLocalXZ	SetTitle ("(MC - RECO) Local Hit XZ;z_{L} / cm;x_{L} / cm")

	for (Int_t j=0;j< nEvents &&j< t->GetEntriesFast();j++)

gStyle	SetOptFit ()

total	SetLineColor (kRed)

total	SetLineWidth (1)

hisDiffX	Fit (g1,"R")

hisDiffX	Fit (g2,"R+")

total	SetParameters (par)

hisDiffX	Fit (total,"R")

hisDiffLocalX	Fit ("gaus")

hisDiffLocalX	GetFunction ("gaus") -> SetLineColor(kRed)

can1	Divide (a, b)

can1	cd (1)

hisxy	DrawCopy ("col")

can1	cd (2)

can1	cd (3)

hisSensxy	DrawCopy ("BOX")

can1	cd (4)

can1	cd (5)

hisPLUV	DrawCopy ()

can1	cd (6)

hisDiffX	DrawCopy ("")

can1	cd (7)

can1	cd (8)

can1	cd (9)

can1	cd (10)

can1	cd (11)

can1	cd (12)

can1	cd (13)

can1	cd (14)

can1	cd (15)

can1	cd (16)

can1	cd (17)

can1	cd (18)

can1	cd (19)

can1	cd (20)

can1	Print ("testoutput.ps")

timer	Stop ()

Variables
bool	verbose = false

TStopwatch	timer

TFile *	f = new TFile("../data/testMCCluster.root")

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

TClonesArray *	hit_array =new TClonesArray("MvdCluster")

TFile *	F = new TFile("testMC.root")

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

TClonesArray *	mc_array =new TClonesArray("MvdMCPoint")

TGeoManager *	geoMan = (TGeoManager*) gDirectory->Get("FAIRGeom")

TH2D *	hisxy = new TH2D("hisxy","",400,-15.,15.,400,-15.,15.)

TH2D *	hisrz = new TH2D("hisrz","",400,-20.,20.,400,-15.,25.)

TH2D *	hisSensxy = new TH2D("hisSensxy","",400,-15.,15.,400,-15.,15.)

TH2D *	hisSensrz = new TH2D("hisSensrz","",400,-20.,20.,400,-15.,25.)

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

TH2D *	hisDiffRZ = new TH2D("hisdiffrz","",100,-0.05,0.05,100,-0.00,0.10)

TH1D *	hisDiffX = new TH1D("hisdiffx","",100,-0.05,0.05)

TH1D *	hisDiffY = new TH1D("hisdiffy","",100,-0.05,0.05)

TH1D *	hisDiffZ = new TH1D("hisdiffz","",100,-0.05,0.05)

TH1D *	hisDiffLocalX = new TH1D("hisdiffLocalx","",100,-0.05,0.05)

TH1D *	hisDiffLocalY = new TH1D("hisdiffLocaly","",100,-0.05,0.05)

TH1D *	hisDiffLocalZ = new TH1D("hisdiffLocalz","",100,-0.05,0.05)

TH2D *	hisDiffLocalXY = new TH2D("hisdiffLocalxy","",100,-0.05,0.05,100,-0.05,0.05)

TH2D *	hisDiffLocalRZ = new TH2D("hisdiffLocalrz","",100,-0.05,0.05,100,-0.00,0.10)

TH2D *	hisLocalXY = new TH2D("hisLocalxy","",100,-5.,5.,100,-5.,5.)

TH2D *	hisLocalRZ = new TH2D("hisLocalrz","",100,-5.,5.,100,-0.,10.)

TH2D *	hisLocalYZ = new TH2D("hisLocalyz","",100,-5.,5.,100,-5.,5.)

TH2D *	hisLocalXZ = new TH2D("hisLocalxz","",100,-5.,5.,100,-5.,5.)

TH1D *	hisPLUV = new TH1D("hispluv","sensor unit vectors",100,0.9999,1.0001)

TVector3	vecs

TVector3	vecmc

TVector3	vecdiff

Double_t	tmpx

Double_t	tmpy

Double_t	tmpz

TVector2	locals

TVector2	localmc

TVector2	localdiff

Double_t	par [6]

TF1 *	g1 = new TF1("g1","gaus",-0.007,0.007)

TF1 *	g2 = new TF1("g2","gaus",-0.03,0.03)

TF1 *	total = new TF1("total","gaus(0)+gaus(3)",-0.03,0.03)

Int_t	a = 5

Int_t	b = 4

TCanvas *	can1 = new TCanvas("MvdTestPlot","MCHit view in MVD",0,0,a250,b250)

Double_t	rtime = timer.RealTime()

Double_t	ctime = timer.CpuTime()

Function Documentation

basiclibs ( )

can1 cd ( 1 )

can1 cd ( 2 )

can1 cd ( 3 )

can1 cd ( 4 )

can1 cd ( 5 )

can1 cd ( 6 )

can1 cd ( 7 )

can1 cd ( 8 )

can1 cd ( 9 )

can1 cd ( 10 )

can1 cd ( 11 )

can1 cd ( 12 )

can1 cd ( 13 )

can1 cd ( 14 )

can1 cd ( 15 )

can1 cd ( 16 )

can1 cd ( 17 )

can1 cd ( 18 )

can1 cd ( 19 )

can1 cd ( 20 )

can1 Divide	(	a	,
		b
	)

hisLocalRZ DrawCopy ( "col" )

hisSensrz DrawCopy ( "BOX" )

hisPLUV DrawCopy ( )

hisDiffX DrawCopy ( "" )

hisDiffX Fit	(	g1	,
		"R"
	)

hisDiffX Fit	(	g2	,
		"R+"
	)

hisDiffZ Fit	(	total	,
		"R"
	)

hisDiffLocalX Fit ( "gaus" )

for ( )

Definition at line 91 of file anaplaneclust.C.

References hit(), i, mcpdg, point, and vecs.

   {
     t->GetEntry(j);
     T->GetEntry(j);
     if(verbose) cout<<"Event No "<<j<<endl;
     for (Int_t i=0; i<hit_array->GetEntriesFast(); i++)
     {
       if(verbose) cout<<"Point No "<<i<<endl;
       PndSdsCluster *hit=(PndSdsCluster*)hit_array->At(i);
       PndSdsMCPoint *point=(PndSdsMCPoint*)mc_array->At(i);
       int mcpdg = -1;
 
       geoMan->cd(point->GetDetName());
 //       currentTransMat = geoMan->GetCurrentMatrix();
 
       vecs.SetXYZ(hit->GetX(), hit->GetY(), hit->GetZ());
       vecmc.SetXYZ(point->GetX(),point->GetY(),point->GetZ());
       vecdiff = vecmc - vecs;
 
       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());
       if(vecs.y() > 0.) hisrz->Fill(vecs.z(),vecs.Perp());
       else hisrz->Fill(vecs.z(),-1.*vecs.Perp());
 
       hisDiffXY->Fill(vecdiff.x(),vecdiff.y());
       hisDiffRZ->Fill(vecdiff.z(),vecdiff.Perp());
       hisDiffX->Fill(vecdiff.x());
       hisDiffY->Fill(vecdiff.y());
       hisDiffZ->Fill(vecdiff.z());
 
       // --- Now move vectors to the local sensor system
       DetPlane plane = hit->getDetPlane();
       TVector3 vnorm = plane->getNormal(); //vnorm->SetMag(1.0);
       TVector3 plO = plane->getO();
       TVector3 plU = plane->getU(); //plU->SetMag(1.0);
       TVector3 plV = plane->getV(); //plV->SetMag(1.0);
       hisSensxy->Fill(plO.x(),plO.y());
       if(plO.y()>0.)hisSensrz->Fill(plO.z(),plO.Perp());
       else hisSensrz->Fill(plO.z(),-1.*plO.Perp());
 
       vecs -= plO;
       tmpz = vecs.Dot(plU);
       tmpx = vecs.Dot(plV);
       tmpy = vecs.Dot(vnorm);
       vecs.SetXYZ(tmpx,tmpy,tmpz);
       hisLocalXY->Fill(vecs.x(),vecs.y());
       hisLocalXZ->Fill(vecs.z(),vecs.x());
       hisLocalYZ->Fill(vecs.z(),vecs.y());
       hisLocalRZ->Fill(vecs.z(),vecs.Perp());
 
       vecmc -= plO;
       tmpz = vecmc.Dot(plU);
       tmpx = vecmc.Dot(plV);
       tmpy = vecmc.Dot(vnorm);
       vecmc.SetXYZ(tmpx,tmpy,tmpz);
       vecdiff = vecmc - vecs;
 
       tmpx = plU.Dot(plU);tmpy = plU.Dot(plV);tmpz = plU.Dot(vnorm);
       hisPLUV->Fill(tmpx);hisPLUV->Fill(tmpy);hisPLUV->Fill(tmpz);
 
       tmpx = plV.Dot(plU);tmpy = plV.Dot(plV);tmpz = plV.Dot(vnorm);
       hisPLUV->Fill(tmpx);hisPLUV->Fill(tmpy);hisPLUV->Fill(tmpz);
 
       tmpx = vnorm.Dot(plU);tmpy = vnorm.Dot(plV);tmpz = vnorm.Dot(vnorm);
       hisPLUV->Fill(tmpx);hisPLUV->Fill(tmpy);hisPLUV->Fill(tmpz);
 
       hisDiffLocalX->Fill(vecdiff.X());
       hisDiffLocalY->Fill(vecdiff.Y());
       hisDiffLocalZ->Fill(vecdiff.Z());
       hisDiffLocalXY->Fill(vecdiff.X(),vecdiff.Y());
       hisDiffLocalRZ->Fill(vecdiff.Z(),vecdiff.Perp());
 
     }//end for i (points in event)
   }// end for j (events)

hisDiffLocalZ GetFunction ( "gaus" ) -> SetLineColor(kRed)

gSystem Load ( "libGeoBase" )

gSystem Load ( "libParBase" )

gSystem Load ( "libBase" )

gSystem Load ( "libPndData" )

gSystem Load ( "libField" )

gSystem Load ( "libGen" )

gSystem Load ( "libPassive" )

gSystem Load ( "libgenfit" )

gSystem Load ( "libtrackrep" )

gSystem Load ( "libtpc" )

gSystem Load ( "libtpcreco" )

gSystem Load ( "librecotasks" )

gSystem Load ( "libMvd" )

gSystem Load ( "libMvdReco" )

gROOT LoadMacro ( "$VMCWORKDIR/gconfig/basiclibs.C" )

can1 Print ( "testoutput.ps" )

t SetBranchAddress	(	"MVDCluster"	,
		&	hit_array
	)

T SetBranchAddress	(	"MVDPoint"	,
		&	mc_array
	)

total SetLineColor ( kRed )

total SetLineWidth ( 1 )

gStyle SetOptFit ( )

total SetParameters ( par )

hisxy SetTitle	(	"MVD MC	Cluster,
		xy view;x/cm;y/cm"
	)

hisrz SetTitle	(	"MVD MC	Cluster,
		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;#Deltax / cm;" )

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

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

hisDiffLocalX SetTitle ( "(MC - RECO) Local Hit coordinate x;#Deltax_{L} / cm;" )

hisDiffLocalY SetTitle ( "(MC - RECO) Local Hit coordinate y;#Deltay_{L} / cm;" )

hisDiffLocalZ SetTitle ( "(MC - RECO) Local Hit coordinate z;#Deltaz_{L} / cm;" )

hisDiffLocalXY SetTitle ( "(MC - RECO) Local Hit XY;#Deltax_{L} / cm;#Deltay_{L} / cm" )

hisDiffLocalRZ SetTitle ( "(MC - RECO) Local Hit RZ;#Deltaz_{L} / cm;#Deltar_{L} / cm" )

hisLocalXY SetTitle ( "(MC - RECO) Local Hit XY;x_{L} / cm;y_{L} / cm" )

hisLocalRZ SetTitle ( "(MC - RECO) Local Hit RZ;z_{L} / cm;r_{L} / cm" )

hisLocalYZ SetTitle ( "(MC - RECO) Local Hit YZ;z_{L} / cm;y_{L} / cm" )

hisLocalXZ SetTitle ( "(MC - RECO) Local Hit XZ;z_{L} / cm;x_{L} / cm" )

timer Start ( )

timer Stop ( )

Variable Documentation

Int_t a = 5

Definition at line 213 of file anaplaneclust.C.

Int_t b = 4

Definition at line 213 of file anaplaneclust.C.

TCanvas* can1 = new TCanvas("MvdTestPlot","MCHit view in MVD",0,0,a*250,b*250)

Definition at line 215 of file anaplaneclust.C.

Double_t ctime = timer.CpuTime()

Definition at line 249 of file anaplaneclust.C.

TFile* f = new TFile("../data/testMCCluster.root")

Definition at line 32 of file anaplaneclust.C.

TFile* F = new TFile("testMC.root")

Definition at line 37 of file anaplaneclust.C.

TF1* g1 = new TF1("g1","gaus",-0.007,0.007)

Definition at line 171 of file anaplaneclust.C.

TF1* g2 = new TF1("g2","gaus",-0.03,0.03)

Definition at line 172 of file anaplaneclust.C.

TGeoManager* geoMan = (TGeoManager*) gDirectory->Get("FAIRGeom")

Definition at line 42 of file anaplaneclust.C.

TH2D* hisDiffLocalRZ = new TH2D("hisdiffLocalrz","",100,-0.05,0.05,100,-0.00,0.10)

Definition at line 73 of file anaplaneclust.C.

TH1D* hisDiffLocalX = new TH1D("hisdiffLocalx","",100,-0.05,0.05)

Definition at line 65 of file anaplaneclust.C.

TH2D* hisDiffLocalXY = new TH2D("hisdiffLocalxy","",100,-0.05,0.05,100,-0.05,0.05)

Definition at line 71 of file anaplaneclust.C.

TH1D* hisDiffLocalY = new TH1D("hisdiffLocaly","",100,-0.05,0.05)

Definition at line 67 of file anaplaneclust.C.

TH1D* hisDiffLocalZ = new TH1D("hisdiffLocalz","",100,-0.05,0.05)

Definition at line 69 of file anaplaneclust.C.

TH2D* hisDiffRZ = new TH2D("hisdiffrz","",100,-0.05,0.05,100,-0.00,0.10)

Definition at line 55 of file anaplaneclust.C.

TH1D* hisDiffX = new TH1D("hisdiffx","",100,-0.05,0.05)

Definition at line 58 of file anaplaneclust.C.

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

Definition at line 53 of file anaplaneclust.C.

TH1D* hisDiffY = new TH1D("hisdiffy","",100,-0.05,0.05)

Definition at line 60 of file anaplaneclust.C.

TH1D* hisDiffZ = new TH1D("hisdiffz","",100,-0.05,0.05)

Definition at line 62 of file anaplaneclust.C.

TH2D* hisLocalRZ = new TH2D("hisLocalrz","",100,-5.,5.,100,-0.,10.)

Definition at line 78 of file anaplaneclust.C.

TH2D* hisLocalXY = new TH2D("hisLocalxy","",100,-5.,5.,100,-5.,5.)

Definition at line 76 of file anaplaneclust.C.

Referenced by anaclust().

TH2D* hisLocalXZ = new TH2D("hisLocalxz","",100,-5.,5.,100,-5.,5.)

Definition at line 82 of file anaplaneclust.C.

TH2D* hisLocalYZ = new TH2D("hisLocalyz","",100,-5.,5.,100,-5.,5.)

Definition at line 80 of file anaplaneclust.C.

TH1D* hisPLUV = new TH1D("hispluv","sensor unit vectors",100,0.9999,1.0001)

Definition at line 85 of file anaplaneclust.C.

TH2D* hisrz = new TH2D("hisrz","",400,-20.,20.,400,-15.,25.)

Definition at line 46 of file anaplaneclust.C.

TH2D* hisSensrz = new TH2D("hisSensrz","",400,-20.,20.,400,-15.,25.)

Definition at line 50 of file anaplaneclust.C.

TH2D* hisSensxy = new TH2D("hisSensxy","",400,-15.,15.,400,-15.,15.)

Definition at line 48 of file anaplaneclust.C.

TH2D* hisxy = new TH2D("hisxy","",400,-15.,15.,400,-15.,15.)

Definition at line 44 of file anaplaneclust.C.

TClonesArray* hit_array =new TClonesArray("MvdCluster")

Definition at line 34 of file anaplaneclust.C.

TVector2 localdiff

Definition at line 89 of file anaplaneclust.C.

TVector2 localmc

Definition at line 89 of file anaplaneclust.C.

TVector2 locals

Definition at line 89 of file anaplaneclust.C.

TClonesArray* mc_array =new TClonesArray("MvdMCPoint")

Definition at line 39 of file anaplaneclust.C.

g2 GetParameters&[3] par

Definition at line 169 of file anaplaneclust.C.

Double_t rtime = timer.RealTime()

Definition at line 248 of file anaplaneclust.C.

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

Definition at line 33 of file anaplaneclust.C.

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

Definition at line 38 of file anaplaneclust.C.

TStopwatch timer

Definition at line 27 of file anaplaneclust.C.

Double_t tmpx

Definition at line 88 of file anaplaneclust.C.

Double_t tmpy

Definition at line 88 of file anaplaneclust.C.

Double_t tmpz

Definition at line 88 of file anaplaneclust.C.

TF1* total = new TF1("total","gaus(0)+gaus(3)",-0.03,0.03)

Definition at line 173 of file anaplaneclust.C.

TVector3 vecdiff

Definition at line 87 of file anaplaneclust.C.

TVector3 vecmc

Definition at line 87 of file anaplaneclust.C.

TVector3 vecs

Definition at line 87 of file anaplaneclust.C.

Referenced by for().

bool verbose = false

Definition at line 4 of file anaplaneclust.C.

Functions

Variables

Function Documentation

Variable Documentation