\r
//\r
if(fAnalysisFolder) delete fAnalysisFolder; fAnalysisFolder=0;\r
+ \r
+ if (fBinsMult) delete[] fBinsMult; fBinsMult=0;\r
+ if (fBinsPt) delete[] fBinsPt; fBinsPt=0;\r
+ if (fBinsPtCorr) delete[] fBinsPtCorr; fBinsPtCorr=0;\r
+ if (fBinsEta) delete[] fBinsEta; fBinsEta=0;\r
+ if (fBinsZv) delete[] fBinsMult; fBinsZv=0;\r
}\r
\r
//_____________________________________________________________________________\r
-void AlidNdPtAnalysis::Init(){\r
- \r
- Int_t multNbins = 27;\r
- Int_t ptNbinsTrackEventCorr = 36;\r
- Int_t ptNbins = 68;\r
- Int_t etaNbins = 30;\r
- Int_t zvNbins = 12;\r
-\r
- Double_t binsMultDefault[28] = {-0.5, 0.5 , 1.5 , 2.5 , 3.5 , 4.5 , 5.5 , 6.5 , 7.5 , 8.5,\r
+void AlidNdPtAnalysis::Init()\r
+{\r
+ //define default binning\r
+ Double_t binsMultDefault[28] = {-0.5, 0.5 , 1.5 , 2.5 , 3.5 , 4.5 , 5.5 , 6.5 , 7.5 , 8.5,\r
9.5, 10.5, 11.5, 12.5, 13.5, 14.5, 15.5, 16.5, 17.5, 18.5,\r
19.5,20.5, 21.5, 22.5, 23.5, 24.5, 29.5, 149.5};\r
-\r
- Double_t binsPtTrackEventCorrDefault[37] = {0.,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45,0.5,0.55,0.6,0.65,0.7,0.75,0.8,0.85,0.9,0.95,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.2,2.4,2.6,3.0,4.0,50.0};\r
-\r
- Double_t binsPtDefault[69] = {0.,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45,0.5,0.55,0.6,0.65,0.7,0.75,0.8,0.85,0.9,0.95,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.2,2.4,2.6,2.8,3.0,3.2,3.4,3.6,3.8,4.0,4.5,5.0,5.5,6.0,6.5,7.0,8.0,9.0,10.0,11.0,12.0,13.0,14.0,15.0,16.0,18.0,20.0,22.0,24.0,26.0,28.0,30.0,32.0,34.0,36.0,40.0,45.0,50.0};\r
-\r
- Double_t binsEtaDefault[31] = {-1.5,-1.4,-1.3,-1.2,-1.1,-1.0,-0.9,-0.8,-0.7,-0.6,-0.5,-0.4,-0.3,-0.2,-0.1,0.,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5};\r
-\r
- Double_t binsZvDefault[13] = {-30.,-25.,-20.,-15.,-10.,-5.,0.,5.,10.,15.,20.,25.,30.};\r
+ Double_t binsPtDefault[69] = {0.,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45,0.5,0.55,0.6,0.65,0.7,0.75,0.8,0.85,0.9,0.95,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.2,2.4,2.6,2.8,3.0,3.2,3.4,3.6,3.8,4.0,4.5,5.0,5.5,6.0,6.5,7.0,8.0,9.0,10.0,11.0,12.0,13.0,14.0,15.0,16.0,18.0,20.0,22.0,24.0,26.0,28.0,30.0,32.0,34.0,36.0,40.0,45.0,50.0};\r
+ Double_t binsPtCorrDefault[37] = {0.,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45,0.5,0.55,0.6,0.65,0.7,0.75,0.8,0.85,0.9,0.95,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.2,2.4,2.6,3.0,4.0,50.0}; \r
+ Double_t binsEtaDefault[31] = {-1.5,-1.4,-1.3,-1.2,-1.1,-1.0,-0.9,-0.8,-0.7,-0.6,-0.5,-0.4,-0.3,-0.2,-0.1,0.,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5};\r
+ Double_t binsZvDefault[13] = {-30.,-25.,-20.,-15.,-10.,-5.,0.,5.,10.,15.,20.,25.,30.};\r
+\r
+ // if no binning is set, use the default\r
+ if (!fBinsMult) { SetBinsMult(27,binsMultDefault); }\r
+ if (!fBinsPt) { SetBinsPt(68,binsPtDefault); }\r
+ if (!fBinsPtCorr) { SetBinsPtCorr(36,binsPtCorrDefault); }\r
+ if (!fBinsEta) { SetBinsEta(30,binsEtaDefault); }\r
+ if (!fBinsZv) { SetBinsZv(12,binsZvDefault); } \r
\r
- Double_t* binsMult = binsMultDefault;\r
- Double_t* binsPtTrackEventCorr = binsPtTrackEventCorrDefault;\r
- Double_t* binsPt = binsPtDefault;\r
- Double_t* binsEta = binsEtaDefault;\r
- Double_t* binsZv = binsZvDefault; \r
-\r
- if (fMultNbins > 0) { multNbins = fMultNbins; binsMult = fBinsMult; }\r
- if (fPtNbins > 0) { ptNbins = fPtNbins; binsPt = fBinsPt; }\r
- if (fPtCorrNbins > 0) { ptNbinsTrackEventCorr = fPtCorrNbins; binsPtTrackEventCorr = fBinsPtCorr; }\r
- if (fEtaNbins > 0) { etaNbins = fEtaNbins; binsEta = fBinsEta; }\r
- if (fZvNbins > 0) { zvNbins = fZvNbins; binsZv = fBinsZv; } \r
-\r
//Int_t binsTrackMatrix[3]={zvNbins,ptNbins,etaNbins};\r
- Int_t binsTrackEventCorrMatrix[3]={zvNbins,ptNbinsTrackEventCorr,etaNbins};\r
-\r
+ Int_t binsTrackEventCorrMatrix[3]={fZvNbins,fPtCorrNbins,fEtaNbins};\r
\r
//\r
// Generic histograms to be corrected\r
//\r
- Int_t binsEventHist[2]={100,150};\r
- Double_t minEventHist[2]={-25.,-0.5}; \r
- Double_t maxEventHist[2]={25.,149.5}; \r
+ Int_t binsEventHist[2]={fZvNbins,fMultNbins};\r
+ //Double_t minEventHist[2]={-fBinsZv[0],fBinsMult[0]}; \r
+ //Double_t maxEventHist[2]={fBinsZv[fZvNbins],fBinsMult[fMultNbins]}; \r
\r
- fRecEventHist = new THnSparseF("fRecEventHist","Zv:multMB",2,binsEventHist,minEventHist,maxEventHist);\r
+ fRecEventHist = new THnSparseF("fRecEventHist","Zv:multMB",2,binsEventHist); //,minEventHist,maxEventHist);\r
+ fRecEventHist->SetBinEdges(0,fBinsZv);\r
+ fRecEventHist->SetBinEdges(1,fBinsMult);\r
fRecEventHist->GetAxis(0)->SetTitle("Zv (cm)");\r
fRecEventHist->GetAxis(1)->SetTitle("multiplicity MB");\r
fRecEventHist->Sumw2();\r
\r
//\r
- Int_t binsTrackHist[4]={100,ptNbins,etaNbins,150};\r
- Double_t minTrackHist[4]={-25.,0.,-1.5,-0.5}; \r
- Double_t maxTrackHist[4]={25.,50.,1.5,149.5}; \r
+ Int_t binsTrackHist[4]={fZvNbins,fPtNbins,fEtaNbins,fMultNbins};\r
+ // Double_t minTrackHist[4]={-25.,0.,-1.5,-0.5}; \r
+ // Double_t maxTrackHist[4]={25.,50.,1.5,149.5}; \r
\r
- fRecTrackHist = new THnSparseF("fRecTrackHist","Zv:pT:eta:multRec",4,binsTrackHist,minTrackHist,maxTrackHist);\r
- fRecTrackHist->SetBinEdges(1,binsPt);\r
- fRecTrackHist->SetBinEdges(2,binsEta);\r
+ fRecTrackHist = new THnSparseF("fRecTrackHist","Zv:pT:eta:multRec",4,binsTrackHist); //,minTrackHist,maxTrackHist);\r
+ fRecTrackHist->SetBinEdges(0,fBinsZv);\r
+ fRecTrackHist->SetBinEdges(1,fBinsPt);\r
+ fRecTrackHist->SetBinEdges(2,fBinsEta);\r
+ fRecTrackHist->SetBinEdges(3,fBinsMult);\r
fRecTrackHist->GetAxis(0)->SetTitle("Zv (cm)");\r
fRecTrackHist->GetAxis(1)->SetTitle("p_{T} (GeV/c)");\r
fRecTrackHist->GetAxis(2)->SetTitle("#eta");\r
//\r
// rec. vs MC correlation matrices\r
//\r
- Int_t binsMultTrueEventMatrix[3]={150,150,150};\r
- Double_t minMultTrueEventMatrix[3]={-0.5,-0.5,-0.5}; \r
- Double_t maxMultTrueEventMatrix[3]={149.5,149.5,149.5}; \r
- fEventMultCorrelationMatrix = new THnSparseF("fEventMultCorrelationMatrix","mult:true_mult:multMB",3,binsMultTrueEventMatrix,minMultTrueEventMatrix,maxMultTrueEventMatrix);\r
+ Int_t binsMultTrueEventMatrix[3]={fMultNbins,fMultNbins,fMultNbins};\r
+// Double_t minMultTrueEventMatrix[3]={-0.5,-0.5,-0.5}; \r
+// Double_t maxMultTrueEventMatrix[3]={149.5,149.5,149.5}; \r
+ fEventMultCorrelationMatrix = new THnSparseF("fEventMultCorrelationMatrix","mult:true_mult:multMB",3,binsMultTrueEventMatrix); //,minMultTrueEventMatrix,maxMultTrueEventMatrix);\r
+ fEventMultCorrelationMatrix->SetBinEdges(0,fBinsMult);\r
+ fEventMultCorrelationMatrix->SetBinEdges(1,fBinsMult);\r
+ fEventMultCorrelationMatrix->SetBinEdges(2,fBinsMult);\r
fEventMultCorrelationMatrix->GetAxis(0)->SetTitle("track multiplicity");\r
fEventMultCorrelationMatrix->GetAxis(1)->SetTitle("true multiplicity");\r
fEventMultCorrelationMatrix->GetAxis(2)->SetTitle("MB multiplicity");\r
fEventMultCorrelationMatrix->Sumw2();\r
\r
- Int_t binsTrackPtCorrelationMatrix[3]={ptNbins,ptNbins,etaNbins};\r
+ Int_t binsTrackPtCorrelationMatrix[3]={fPtCorrNbins,fPtCorrNbins,fEtaNbins};\r
fTrackPtCorrelationMatrix = new THnSparseF("fTrackPtCorrelationMatrix","Pt:mcPt:mcEta",3,binsTrackPtCorrelationMatrix);\r
- fTrackPtCorrelationMatrix->SetBinEdges(0,binsPt);\r
- fTrackPtCorrelationMatrix->SetBinEdges(1,binsPt);\r
- fTrackPtCorrelationMatrix->SetBinEdges(2,binsEta);\r
+ fTrackPtCorrelationMatrix->SetBinEdges(0,fBinsPtCorr);\r
+ fTrackPtCorrelationMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fTrackPtCorrelationMatrix->SetBinEdges(2,fBinsEta);\r
fTrackPtCorrelationMatrix->GetAxis(0)->SetTitle("Pt (GeV/c)");\r
fTrackPtCorrelationMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fTrackPtCorrelationMatrix->GetAxis(2)->SetTitle("mcEta");\r
//\r
// Efficiency and contamination correction matrices\r
//\r
- Int_t binsEventMatrix[2]={zvNbins,multNbins};\r
- Double_t minEventMatrix[2]={-25.,-0.5}; \r
- Double_t maxEventMatrix[2]={25.,149.5}; \r
+ Int_t binsEventMatrix[2]={fZvNbins,fMultNbins};\r
+// Double_t minEventMatrix[2]={-25.,-0.5}; \r
+// Double_t maxEventMatrix[2]={25.,149.5}; \r
\r
- fGenEventMatrix = new THnSparseF("fGenEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fGenEventMatrix->SetBinEdges(0,binsZv);\r
- fGenEventMatrix->SetBinEdges(1,binsMult);\r
+ fGenEventMatrix = new THnSparseF("fGenEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fGenEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fGenEventMatrix->SetBinEdges(1,fBinsMult);\r
fGenEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fGenEventMatrix->Sumw2();\r
\r
- fGenSDEventMatrix = new THnSparseF("fGenSDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fGenSDEventMatrix->SetBinEdges(0,binsZv);\r
- fGenSDEventMatrix->SetBinEdges(1,binsMult);\r
+ fGenSDEventMatrix = new THnSparseF("fGenSDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fGenSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fGenSDEventMatrix->SetBinEdges(1,fBinsMult);\r
fGenSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenSDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fGenSDEventMatrix->Sumw2();\r
\r
- fGenDDEventMatrix = new THnSparseF("fGenDDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fGenDDEventMatrix->SetBinEdges(0,binsZv);\r
- fGenDDEventMatrix->SetBinEdges(1,binsMult);\r
+ fGenDDEventMatrix = new THnSparseF("fGenDDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fGenDDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fGenDDEventMatrix->SetBinEdges(1,fBinsMult);\r
fGenDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenDDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fGenDDEventMatrix->Sumw2();\r
\r
- fGenNDEventMatrix = new THnSparseF("fGenNDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fGenNDEventMatrix->SetBinEdges(0,binsZv);\r
- fGenNDEventMatrix->SetBinEdges(1,binsMult);\r
+ fGenNDEventMatrix = new THnSparseF("fGenNDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fGenNDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fGenNDEventMatrix->SetBinEdges(1,fBinsMult);\r
fGenNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenNDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fGenNDEventMatrix->Sumw2();\r
\r
- fGenNSDEventMatrix = new THnSparseF("fGenNSDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fGenNSDEventMatrix->SetBinEdges(0,binsZv);\r
- fGenNSDEventMatrix->SetBinEdges(1,binsMult);\r
+ fGenNSDEventMatrix = new THnSparseF("fGenNSDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fGenNSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fGenNSDEventMatrix->SetBinEdges(1,fBinsMult);\r
fGenNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenNSDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fGenNSDEventMatrix->Sumw2();\r
\r
//\r
- fTriggerEventMatrix = new THnSparseF("fTriggerEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fTriggerEventMatrix->SetBinEdges(0,binsZv);\r
- fTriggerEventMatrix->SetBinEdges(1,binsMult);\r
+ fTriggerEventMatrix = new THnSparseF("fTriggerEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fTriggerEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fTriggerEventMatrix->SetBinEdges(1,fBinsMult);\r
fTriggerEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fTriggerEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fTriggerEventMatrix->Sumw2();\r
\r
- fTriggerSDEventMatrix = new THnSparseF("fTriggerSDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fTriggerSDEventMatrix->SetBinEdges(0,binsZv);\r
- fTriggerSDEventMatrix->SetBinEdges(1,binsMult);\r
+ fTriggerSDEventMatrix = new THnSparseF("fTriggerSDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fTriggerSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fTriggerSDEventMatrix->SetBinEdges(1,fBinsMult);\r
fTriggerSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fTriggerSDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fTriggerSDEventMatrix->Sumw2();\r
\r
- fTriggerDDEventMatrix = new THnSparseF("fTriggerDDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fTriggerDDEventMatrix->SetBinEdges(0,binsZv);\r
- fTriggerDDEventMatrix->SetBinEdges(1,binsMult);\r
+ fTriggerDDEventMatrix = new THnSparseF("fTriggerDDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fTriggerDDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fTriggerDDEventMatrix->SetBinEdges(1,fBinsMult);\r
fTriggerDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fTriggerDDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fTriggerDDEventMatrix->Sumw2();\r
\r
- fTriggerNDEventMatrix = new THnSparseF("fTriggerNDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fTriggerNDEventMatrix->SetBinEdges(0,binsZv);\r
- fTriggerNDEventMatrix->SetBinEdges(1,binsMult);\r
+ fTriggerNDEventMatrix = new THnSparseF("fTriggerNDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fTriggerNDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fTriggerNDEventMatrix->SetBinEdges(1,fBinsMult);\r
fTriggerNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fTriggerNDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fTriggerNDEventMatrix->Sumw2();\r
\r
- fTriggerNSDEventMatrix = new THnSparseF("fTriggerNSDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fTriggerNSDEventMatrix->SetBinEdges(0,binsZv);\r
- fTriggerNSDEventMatrix->SetBinEdges(1,binsMult);\r
+ fTriggerNSDEventMatrix = new THnSparseF("fTriggerNSDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fTriggerNSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fTriggerNSDEventMatrix->SetBinEdges(1,fBinsMult);\r
fTriggerNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fTriggerNSDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fTriggerNSDEventMatrix->Sumw2();\r
\r
//\r
- fRecEventMatrix = new THnSparseF("fRecEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fRecEventMatrix->SetBinEdges(0,binsZv);\r
- fRecEventMatrix->SetBinEdges(1,binsMult);\r
+ fRecEventMatrix = new THnSparseF("fRecEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fRecEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fRecEventMatrix->SetBinEdges(1,fBinsMult);\r
fRecEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fRecEventMatrix->Sumw2();\r
\r
- fRecSDEventMatrix = new THnSparseF("fRecSDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fRecSDEventMatrix->SetBinEdges(0,binsZv);\r
- fRecSDEventMatrix->SetBinEdges(1,binsMult);\r
+ fRecSDEventMatrix = new THnSparseF("fRecSDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fRecSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fRecSDEventMatrix->SetBinEdges(1,fBinsMult);\r
fRecSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecSDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fRecSDEventMatrix->Sumw2();\r
\r
- fRecDDEventMatrix = new THnSparseF("fRecDDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fRecDDEventMatrix->SetBinEdges(0,binsZv);\r
- fRecDDEventMatrix->SetBinEdges(1,binsMult);\r
+ fRecDDEventMatrix = new THnSparseF("fRecDDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fRecDDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fRecDDEventMatrix->SetBinEdges(1,fBinsMult);\r
fRecDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecDDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fRecDDEventMatrix->Sumw2();\r
\r
- fRecNDEventMatrix = new THnSparseF("fRecNDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fRecNDEventMatrix->SetBinEdges(0,binsZv);\r
- fRecNDEventMatrix->SetBinEdges(1,binsMult);\r
+ fRecNDEventMatrix = new THnSparseF("fRecNDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fRecNDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fRecNDEventMatrix->SetBinEdges(1,fBinsMult);\r
fRecNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecNDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fRecNDEventMatrix->Sumw2();\r
\r
- fRecNSDEventMatrix = new THnSparseF("fRecNSDEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fRecNSDEventMatrix->SetBinEdges(0,binsZv);\r
- fRecNSDEventMatrix->SetBinEdges(1,binsMult);\r
+ fRecNSDEventMatrix = new THnSparseF("fRecNSDEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fRecNSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fRecNSDEventMatrix->SetBinEdges(1,fBinsMult);\r
fRecNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecNSDEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fRecNSDEventMatrix->Sumw2();\r
\r
- fRecCandleEventMatrix = new THnSparseF("fRecCandleEventMatrix","mcZv:multMB",2,binsEventMatrix,minEventMatrix,maxEventMatrix);\r
- fRecCandleEventMatrix->SetBinEdges(0,binsZv);\r
- fRecCandleEventMatrix->SetBinEdges(1,binsMult);\r
+ fRecCandleEventMatrix = new THnSparseF("fRecCandleEventMatrix","mcZv:multMB",2,binsEventMatrix); //,minEventMatrix,maxEventMatrix);\r
+ fRecCandleEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fRecCandleEventMatrix->SetBinEdges(1,fBinsMult);\r
fRecCandleEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecCandleEventMatrix->GetAxis(1)->SetTitle("multiplicity MB");\r
fRecCandleEventMatrix->Sumw2();\r
//\r
\r
fGenTrackEventMatrix = new THnSparseF("fGenTrackEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fGenTrackEventMatrix->SetBinEdges(0,binsZv);\r
- fGenTrackEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fGenTrackEventMatrix->SetBinEdges(2,binsEta);\r
+ fGenTrackEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fGenTrackEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fGenTrackEventMatrix->SetBinEdges(2,fBinsEta);\r
fGenTrackEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenTrackEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fGenTrackEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fGenTrackEventMatrix->Sumw2();\r
\r
fGenTrackSDEventMatrix = new THnSparseF("fGenTrackSDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fGenTrackSDEventMatrix->SetBinEdges(0,binsZv);\r
- fGenTrackSDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fGenTrackSDEventMatrix->SetBinEdges(2,binsEta);\r
+ fGenTrackSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fGenTrackSDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fGenTrackSDEventMatrix->SetBinEdges(2,fBinsEta);\r
fGenTrackSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenTrackSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fGenTrackSDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fGenTrackSDEventMatrix->Sumw2();\r
\r
fGenTrackDDEventMatrix = new THnSparseF("fGenTrackDDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fGenTrackDDEventMatrix->SetBinEdges(0,binsZv);\r
- fGenTrackDDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fGenTrackDDEventMatrix->SetBinEdges(2,binsEta);\r
+ fGenTrackDDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fGenTrackDDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fGenTrackDDEventMatrix->SetBinEdges(2,fBinsEta);\r
fGenTrackDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenTrackDDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fGenTrackDDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fGenTrackDDEventMatrix->Sumw2();\r
\r
fGenTrackNDEventMatrix = new THnSparseF("fGenTrackNDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fGenTrackNDEventMatrix->SetBinEdges(0,binsZv);\r
- fGenTrackNDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fGenTrackNDEventMatrix->SetBinEdges(2,binsEta);\r
+ fGenTrackNDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fGenTrackNDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fGenTrackNDEventMatrix->SetBinEdges(2,fBinsEta);\r
fGenTrackNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenTrackNDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fGenTrackNDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fGenTrackNDEventMatrix->Sumw2();\r
\r
fGenTrackNSDEventMatrix = new THnSparseF("fGenTrackNSDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fGenTrackNSDEventMatrix->SetBinEdges(0,binsZv);\r
- fGenTrackNSDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fGenTrackNSDEventMatrix->SetBinEdges(2,binsEta);\r
+ fGenTrackNSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fGenTrackNSDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fGenTrackNSDEventMatrix->SetBinEdges(2,fBinsEta);\r
fGenTrackNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenTrackNSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fGenTrackNSDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
\r
//\r
fTriggerTrackEventMatrix = new THnSparseF("fTriggerTrackEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fTriggerTrackEventMatrix->SetBinEdges(0,binsZv);\r
- fTriggerTrackEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fTriggerTrackEventMatrix->SetBinEdges(2,binsEta);\r
+ fTriggerTrackEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fTriggerTrackEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fTriggerTrackEventMatrix->SetBinEdges(2,fBinsEta);\r
fTriggerTrackEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fTriggerTrackEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fTriggerTrackEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fTriggerTrackEventMatrix->Sumw2();\r
\r
fTriggerTrackSDEventMatrix = new THnSparseF("fTriggerTrackSDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fTriggerTrackSDEventMatrix->SetBinEdges(0,binsZv);\r
- fTriggerTrackSDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fTriggerTrackSDEventMatrix->SetBinEdges(2,binsEta);\r
+ fTriggerTrackSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fTriggerTrackSDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fTriggerTrackSDEventMatrix->SetBinEdges(2,fBinsEta);\r
fTriggerTrackSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fTriggerTrackSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fTriggerTrackSDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fTriggerTrackSDEventMatrix->Sumw2();\r
\r
fTriggerTrackDDEventMatrix = new THnSparseF("fTriggerTrackDDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fTriggerTrackDDEventMatrix->SetBinEdges(0,binsZv);\r
- fTriggerTrackDDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fTriggerTrackDDEventMatrix->SetBinEdges(2,binsEta);\r
+ fTriggerTrackDDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fTriggerTrackDDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fTriggerTrackDDEventMatrix->SetBinEdges(2,fBinsEta);\r
fTriggerTrackDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fTriggerTrackDDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fTriggerTrackDDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fTriggerTrackDDEventMatrix->Sumw2();\r
\r
fTriggerTrackNDEventMatrix = new THnSparseF("fTriggerTrackNDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fTriggerTrackNDEventMatrix->SetBinEdges(0,binsZv);\r
- fTriggerTrackNDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fTriggerTrackNDEventMatrix->SetBinEdges(2,binsEta);\r
+ fTriggerTrackNDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fTriggerTrackNDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fTriggerTrackNDEventMatrix->SetBinEdges(2,fBinsEta);\r
fTriggerTrackNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fTriggerTrackNDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fTriggerTrackNDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fTriggerTrackNDEventMatrix->Sumw2();\r
\r
fTriggerTrackNSDEventMatrix = new THnSparseF("fTriggerTrackNSDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fTriggerTrackNSDEventMatrix->SetBinEdges(0,binsZv);\r
- fTriggerTrackNSDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fTriggerTrackNSDEventMatrix->SetBinEdges(2,binsEta);\r
+ fTriggerTrackNSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fTriggerTrackNSDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fTriggerTrackNSDEventMatrix->SetBinEdges(2,fBinsEta);\r
fTriggerTrackNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fTriggerTrackNSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fTriggerTrackNSDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
\r
//\r
fRecTrackEventMatrix = new THnSparseF("fRecTrackEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fRecTrackEventMatrix->SetBinEdges(0,binsZv);\r
- fRecTrackEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fRecTrackEventMatrix->SetBinEdges(2,binsEta);\r
+ fRecTrackEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fRecTrackEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fRecTrackEventMatrix->SetBinEdges(2,fBinsEta);\r
fRecTrackEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecTrackEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fRecTrackEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fRecTrackEventMatrix->Sumw2();\r
\r
fRecTrackSDEventMatrix = new THnSparseF("fRecTrackSDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fRecTrackSDEventMatrix->SetBinEdges(0,binsZv);\r
- fRecTrackSDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fRecTrackSDEventMatrix->SetBinEdges(2,binsEta);\r
+ fRecTrackSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fRecTrackSDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fRecTrackSDEventMatrix->SetBinEdges(2,fBinsEta);\r
fRecTrackSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecTrackSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fRecTrackSDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fRecTrackSDEventMatrix->Sumw2();\r
\r
fRecTrackDDEventMatrix = new THnSparseF("fRecTrackDDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fRecTrackDDEventMatrix->SetBinEdges(0,binsZv);\r
- fRecTrackDDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fRecTrackDDEventMatrix->SetBinEdges(2,binsEta);\r
+ fRecTrackDDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fRecTrackDDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fRecTrackDDEventMatrix->SetBinEdges(2,fBinsEta);\r
fRecTrackDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecTrackDDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fRecTrackDDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fRecTrackDDEventMatrix->Sumw2();\r
\r
fRecTrackNDEventMatrix = new THnSparseF("fRecTrackNDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fRecTrackNDEventMatrix->SetBinEdges(0,binsZv);\r
- fRecTrackNDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fRecTrackNDEventMatrix->SetBinEdges(2,binsEta);\r
+ fRecTrackNDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fRecTrackNDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fRecTrackNDEventMatrix->SetBinEdges(2,fBinsEta);\r
fRecTrackNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecTrackNDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fRecTrackNDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
fRecTrackNDEventMatrix->Sumw2();\r
\r
fRecTrackNSDEventMatrix = new THnSparseF("fRecTrackNSDEventMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fRecTrackNSDEventMatrix->SetBinEdges(0,binsZv);\r
- fRecTrackNSDEventMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fRecTrackNSDEventMatrix->SetBinEdges(2,binsEta);\r
+ fRecTrackNSDEventMatrix->SetBinEdges(0,fBinsZv);\r
+ fRecTrackNSDEventMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fRecTrackNSDEventMatrix->SetBinEdges(2,fBinsEta);\r
fRecTrackNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecTrackNSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fRecTrackNSDEventMatrix->GetAxis(2)->SetTitle("mcEta");\r
//\r
//fGenTrackMatrix = new THnSparseF("fGenTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackMatrix);\r
fGenTrackMatrix = new THnSparseF("fGenTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fGenTrackMatrix->SetBinEdges(0,binsZv);\r
+ fGenTrackMatrix->SetBinEdges(0,fBinsZv);\r
//fGenTrackMatrix->SetBinEdges(1,binsPt);\r
- fGenTrackMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fGenTrackMatrix->SetBinEdges(2,binsEta);\r
+ fGenTrackMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fGenTrackMatrix->SetBinEdges(2,fBinsEta);\r
fGenTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenTrackMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fGenTrackMatrix->GetAxis(2)->SetTitle("mcEta");\r
\r
//fGenPrimTrackMatrix = new THnSparseF("fGenPrimTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackMatrix);\r
fGenPrimTrackMatrix = new THnSparseF("fGenPrimTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fGenPrimTrackMatrix->SetBinEdges(0,binsZv);\r
+ fGenPrimTrackMatrix->SetBinEdges(0,fBinsZv);\r
//fGenPrimTrackMatrix->SetBinEdges(1,binsPt);\r
- fGenPrimTrackMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fGenPrimTrackMatrix->SetBinEdges(2,binsEta);\r
+ fGenPrimTrackMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fGenPrimTrackMatrix->SetBinEdges(2,fBinsEta);\r
fGenPrimTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fGenPrimTrackMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fGenPrimTrackMatrix->GetAxis(2)->SetTitle("mcEta");\r
\r
//fRecPrimTrackMatrix = new THnSparseF("fRecPrimTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackMatrix);\r
fRecPrimTrackMatrix = new THnSparseF("fRecPrimTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fRecPrimTrackMatrix->SetBinEdges(0,binsZv);\r
+ fRecPrimTrackMatrix->SetBinEdges(0,fBinsZv);\r
//fRecPrimTrackMatrix->SetBinEdges(1,binsPt);\r
- fRecPrimTrackMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fRecPrimTrackMatrix->SetBinEdges(2,binsEta);\r
+ fRecPrimTrackMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fRecPrimTrackMatrix->SetBinEdges(2,fBinsEta);\r
fRecPrimTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecPrimTrackMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fRecPrimTrackMatrix->GetAxis(2)->SetTitle("mcEta");\r
//\r
//fRecTrackMatrix = new THnSparseF("fRecTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackMatrix);\r
fRecTrackMatrix = new THnSparseF("fRecTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fRecTrackMatrix->SetBinEdges(0,binsZv);\r
+ fRecTrackMatrix->SetBinEdges(0,fBinsZv);\r
//fRecTrackMatrix->SetBinEdges(1,binsPt);\r
- fRecTrackMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fRecTrackMatrix->SetBinEdges(2,binsEta);\r
+ fRecTrackMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fRecTrackMatrix->SetBinEdges(2,fBinsEta);\r
fRecTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecTrackMatrix->GetAxis(1)->SetTitle("Pt (GeV/c)");\r
fRecTrackMatrix->GetAxis(2)->SetTitle("Eta");\r
\r
//fRecSecTrackMatrix = new THnSparseF("fRecSecTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackMatrix);\r
fRecSecTrackMatrix = new THnSparseF("fRecSecTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fRecSecTrackMatrix->SetBinEdges(0,binsZv);\r
+ fRecSecTrackMatrix->SetBinEdges(0,fBinsZv);\r
//fRecSecTrackMatrix->SetBinEdges(1,binsPt);\r
- fRecSecTrackMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fRecSecTrackMatrix->SetBinEdges(2,binsEta);\r
+ fRecSecTrackMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fRecSecTrackMatrix->SetBinEdges(2,fBinsEta);\r
fRecSecTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecSecTrackMatrix->GetAxis(1)->SetTitle("Pt (GeV/c)");\r
fRecSecTrackMatrix->GetAxis(2)->SetTitle("Eta");\r
//\r
//fRecMultTrackMatrix = new THnSparseF("fRecMultTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackMatrix);\r
fRecMultTrackMatrix = new THnSparseF("fRecMultTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);\r
- fRecMultTrackMatrix->SetBinEdges(0,binsZv);\r
+ fRecMultTrackMatrix->SetBinEdges(0,fBinsZv);\r
//fRecMultTrackMatrix->SetBinEdges(1,binsPt);\r
- fRecMultTrackMatrix->SetBinEdges(1,binsPtTrackEventCorr);\r
- fRecMultTrackMatrix->SetBinEdges(2,binsEta);\r
+ fRecMultTrackMatrix->SetBinEdges(1,fBinsPtCorr);\r
+ fRecMultTrackMatrix->SetBinEdges(2,fBinsEta);\r
fRecMultTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");\r
fRecMultTrackMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");\r
fRecMultTrackMatrix->GetAxis(2)->SetTitle("mcEta");\r
fRecEventHist1->Sumw2();\r
\r
//\r
- Int_t binsRecEventHist2[3]={zvNbins,150,150};\r
- Double_t minRecEventHist2[3]={-25.,-0.5,-0.5}; \r
- Double_t maxRecEventHist2[3]={25.,149.5,149.5}; \r
+ Int_t binsRecEventHist2[3]={fZvNbins,fMultNbins,fMultNbins};\r
+// Double_t minRecEventHist2[3]={-25.,-0.5,-0.5}; \r
+// Double_t maxRecEventHist2[3]={25.,149.5,149.5}; \r
\r
- fRecEventHist2 = new THnSparseF("fRecEventHist2","Zv:multMB:mult",3,binsRecEventHist2,minRecEventHist2,maxRecEventHist2);\r
- fRecEventHist2->SetBinEdges(0,binsZv);\r
+ fRecEventHist2 = new THnSparseF("fRecEventHist2","Zv:multMB:mult",3,binsRecEventHist2); //,minRecEventHist2,maxRecEventHist2);\r
+ fRecEventHist2->SetBinEdges(0,fBinsZv);\r
+ fRecEventHist2->SetBinEdges(1,fBinsMult);\r
+ fRecEventHist2->SetBinEdges(2,fBinsMult);\r
fRecEventHist2->GetAxis(0)->SetTitle("Zv (cm)");\r
fRecEventHist2->GetAxis(1)->SetTitle("multiplicity MB");\r
fRecEventHist2->GetAxis(2)->SetTitle("multiplicity");\r
fRecMCEventHist1->Sumw2();\r
\r
//\r
- Int_t binsRecMCEventHist2[3]={100,100,150};\r
+ Int_t binsRecMCEventHist2[3]={100,100,fMultNbins};\r
Double_t minRecMCEventHist2[3]={-10.0*kFact,-10.0*kFact,0.0}; \r
Double_t maxRecMCEventHist2[3]={10.0*kFact,10.0*kFact,149.50}; \r
\r
fRecMCEventHist2 = new THnSparseF("fRecMCEventHist2","Xv-mcXv:Zv-mcZv:mult",3,binsRecMCEventHist2,minRecMCEventHist2,maxRecMCEventHist2);\r
+ fRecMCEventHist2->SetBinEdges(2,fBinsMult); \r
fRecMCEventHist2->GetAxis(0)->SetTitle("Xv-mcXv (cm)");\r
fRecMCEventHist2->GetAxis(1)->SetTitle("Zv-mcZv (cm)");\r
fRecMCEventHist2->GetAxis(2)->SetTitle("multiplicity");\r
fRecMCEventHist2->Sumw2();\r
\r
- Int_t binsRecMCEventHist3[2]={150,5};\r
+ Int_t binsRecMCEventHist3[2]={fMultNbins,5};\r
Double_t minRecMCEventHist3[2]={-0.5,0.0}; \r
Double_t maxRecMCEventHist3[2]={149.50,5.0}; \r
fRecMCEventHist3 = new THnSparseF("fRecMCEventHist3","mult:EventType (ND, DD, SD)",2,binsRecMCEventHist3,minRecMCEventHist3,maxRecMCEventHist3);\r
+ fRecMCEventHist3->SetBinEdges(0,fBinsMult); \r
fRecMCEventHist3->GetAxis(0)->SetTitle("multiplicity");\r
fRecMCEventHist3->GetAxis(1)->SetTitle("EventType");\r
fRecMCEventHist3->Sumw2();\r
{\r
// THnSparse track histograms\r
\r
- Int_t binsMCTrackHist1[3]= {ptNbins, etaNbins, 90};\r
+ Int_t binsMCTrackHist1[3]= {fPtCorrNbins, fEtaNbins, 90};\r
Double_t minMCTrackHist1[3]={0.,-1.,0.}; \r
Double_t maxMCTrackHist1[3]={10.,1.,2.*TMath::Pi()}; \r
sprintf(name,"fMCTrackHist1_%d",i);\r
sprintf(title,"mcPt:mcEta:mcPhi");\r
\r
fMCTrackHist1[i] = new THnSparseF(name,title,3,binsMCTrackHist1,minMCTrackHist1,maxMCTrackHist1);\r
- fMCTrackHist1[i]->SetBinEdges(0,binsPt);\r
- fMCTrackHist1[i]->SetBinEdges(1,binsEta);\r
+ fMCTrackHist1[i]->SetBinEdges(0,fBinsPtCorr);\r
+ fMCTrackHist1[i]->SetBinEdges(1,fBinsEta);\r
fMCTrackHist1[i]->GetAxis(0)->SetTitle("mcPt (GeV/c)");\r
fMCTrackHist1[i]->GetAxis(1)->SetTitle("mcEta");\r
fMCTrackHist1[i]->GetAxis(2)->SetTitle("mcPhi (rad)");\r
fMCTrackHist1[i]->Sumw2();\r
\r
- Int_t binsMCPrimTrackHist1[5]= {ptNbins,etaNbins,6,20,4000};\r
+ Int_t binsMCPrimTrackHist1[5]= {fPtCorrNbins,fEtaNbins,6,20,4000};\r
Double_t minMCPrimTrackHist1[5]={0.,-1.,0.,0.,0.}; \r
Double_t maxMCPrimTrackHist1[5]={10.,1.,6.,20.,4000.}; \r
sprintf(name,"fMCPrimTrackHist1_%d",i);\r
sprintf(title,"mcPt:mcEta:pid:mech:mother");\r
\r
fMCPrimTrackHist1[i] = new THnSparseF(name,title,5,binsMCPrimTrackHist1,minMCPrimTrackHist1,maxMCPrimTrackHist1);\r
- fMCPrimTrackHist1[i]->SetBinEdges(0,binsPt);\r
- fMCPrimTrackHist1[i]->SetBinEdges(1,binsEta);\r
+ fMCPrimTrackHist1[i]->SetBinEdges(0,fBinsPtCorr);\r
+ fMCPrimTrackHist1[i]->SetBinEdges(1,fBinsEta);\r
fMCPrimTrackHist1[i]->GetAxis(0)->SetTitle("mcPt (GeV/c)");\r
fMCPrimTrackHist1[i]->GetAxis(1)->SetTitle("mcEta");\r
fMCPrimTrackHist1[i]->GetAxis(2)->SetTitle("pid");\r
fMCPrimTrackHist2[i]->GetAxis(2)->SetTitle("mother");\r
fMCPrimTrackHist2[i]->Sumw2();\r
\r
- Int_t binsMCSecTrackHist1[5]= {ptNbins,etaNbins,6,20,4000};\r
+ Int_t binsMCSecTrackHist1[5]= {fPtCorrNbins,fEtaNbins,6,20,4000};\r
Double_t minMCSecTrackHist1[5]={0.,-1.,0.,0.,0.}; \r
- Double_t maxMCSecTrackHist1[5]={10.,1.,6.,20.,4000.}; \r
+ Double_t maxMCSecTrackHist1[5]={10.,1.,6.,20.,4000.}; \r
sprintf(name,"fMCSecTrackHist1_%d",i);\r
sprintf(title,"mcPt:mcEta:pid:mech:mother");\r
\r
fMCSecTrackHist1[i] = new THnSparseF(name,title,5,binsMCSecTrackHist1,minMCSecTrackHist1,maxMCSecTrackHist1);\r
- fMCSecTrackHist1[i]->SetBinEdges(0,binsPt);\r
- fMCSecTrackHist1[i]->SetBinEdges(1,binsEta);\r
+ fMCSecTrackHist1[i]->SetBinEdges(0,fBinsPtCorr);\r
+ fMCSecTrackHist1[i]->SetBinEdges(1,fBinsEta);\r
fMCSecTrackHist1[i]->GetAxis(0)->SetTitle("mcPt (GeV/c)");\r
fMCSecTrackHist1[i]->GetAxis(1)->SetTitle("mcEta");\r
fMCSecTrackHist1[i]->GetAxis(2)->SetTitle("pid");\r
\r
// \r
\r
- Int_t binsRecTrackHist1[3]={ptNbins,etaNbins,90};\r
+ Int_t binsRecTrackHist1[3]={fPtNbins,fEtaNbins,90};\r
Double_t minRecTrackHist1[3]={0.,-1.,0.}; \r
Double_t maxRecTrackHist1[3]={10.,1.,2.*TMath::Pi()};\r
sprintf(name,"fRecTrackHist1_%d",i);\r
sprintf(title,"Pt:Eta:Phi");\r
fRecTrackHist1[i] = new THnSparseF(name,title,3,binsRecTrackHist1,minRecTrackHist1,maxRecTrackHist1);\r
- fRecTrackHist1[i]->SetBinEdges(0,binsPt);\r
- fRecTrackHist1[i]->SetBinEdges(1,binsEta);\r
+ fRecTrackHist1[i]->SetBinEdges(0,fBinsPt);\r
+ fRecTrackHist1[i]->SetBinEdges(1,fBinsEta);\r
fRecTrackHist1[i]->GetAxis(0)->SetTitle("p_{T} (GeV/c)");\r
fRecTrackHist1[i]->GetAxis(1)->SetTitle("#eta");\r
fRecTrackHist1[i]->GetAxis(2)->SetTitle("#phi (rad)");\r
fRecTrackHist1[i]->Sumw2();\r
\r
// \r
- Int_t binsRecTrackMultHist1[2]={ptNbins,150};\r
- Double_t minRecTrackMultHist1[2]={0.,-0.5}; \r
- Double_t maxRecTrackMultHist1[2]={10.,149.5};\r
+ Int_t binsRecTrackMultHist1[2]={fPtNbins,fMultNbins};\r
+// Double_t minRecTrackMultHist1[2]={0.,-0.5}; \r
+// Double_t maxRecTrackMultHist1[2]={10.,149.5};\r
sprintf(name,"fRecTrackMultHist_%d",i);\r
sprintf(title,"Pt:Mult");\r
- fRecTrackMultHist1[i] = new THnSparseF(name,title,2,binsRecTrackMultHist1,minRecTrackMultHist1,maxRecTrackMultHist1);\r
- fRecTrackMultHist1[i]->SetBinEdges(0,binsPt);\r
+ fRecTrackMultHist1[i] = new THnSparseF(name,title,2,binsRecTrackMultHist1); //,minRecTrackMultHist1,maxRecTrackMultHist1);\r
+ fRecTrackMultHist1[i]->SetBinEdges(0,fBinsPt);\r
+ fRecTrackMultHist1[i]->SetBinEdges(1,fBinsMult);\r
fRecTrackMultHist1[i]->GetAxis(0)->SetTitle("Pt (GeV/c)");\r
fRecTrackMultHist1[i]->GetAxis(1)->SetTitle("multiplicity");\r
fRecTrackMultHist1[i]->Sumw2();\r
}\r
\r
- Int_t binsRecMCTrackHist1[4] = {ptNbins,etaNbins,100,100};\r
+ Int_t binsRecMCTrackHist1[4] = {fPtCorrNbins,fEtaNbins,100,100};\r
Double_t minRecMCTrackHist1[4]={0.,-1.,-0.5,-0.5}; \r
Double_t maxRecMCTrackHist1[4]={20.,1.,0.5,0.5}; \r
sprintf(name,"fRecMCTrackHist1");\r
sprintf(title,"mcPt:mcEta:(Pt-mcPt)/mcPt:(Eta-mcEta)");\r
fRecMCTrackHist1 = new THnSparseF(name,title,4,binsRecMCTrackHist1,minRecMCTrackHist1,maxRecMCTrackHist1);\r
- fRecMCTrackHist1->SetBinEdges(0,binsPt);\r
- fRecMCTrackHist1->SetBinEdges(1,binsEta);\r
+ fRecMCTrackHist1->SetBinEdges(0,fBinsPtCorr);\r
+ fRecMCTrackHist1->SetBinEdges(1,fBinsEta);\r
fRecMCTrackHist1->GetAxis(0)->SetTitle("mcPt (GeV/c)");\r
fRecMCTrackHist1->GetAxis(1)->SetTitle("mcEta");\r
fRecMCTrackHist1->GetAxis(2)->SetTitle("(Pt-mcPt)/mcPt");\r
fRecMCTrackHist1->GetAxis(3)->SetTitle("Eta-mcEta");\r
\r
- Int_t binsMCMultRecTrackHist1[3] = {ptNbins,etaNbins,6};\r
+ Int_t binsMCMultRecTrackHist1[3] = {fPtCorrNbins,fEtaNbins,6};\r
Double_t minMCMultRecTrackHist1[3]={0.,-1.,0.}; \r
Double_t maxMCMultRecTrackHist1[3]={20.,1.,6.}; \r
sprintf(name,"fMCMultRecTrackHist1");\r
sprintf(title,"mcPt:mcEta:pid");\r
fMCMultRecTrackHist1 = new THnSparseF(name,title,3,binsMCMultRecTrackHist1,minMCMultRecTrackHist1,maxMCMultRecTrackHist1);\r
- fMCMultRecTrackHist1->SetBinEdges(0,binsPt);\r
- fMCMultRecTrackHist1->SetBinEdges(1,binsEta);\r
+ fMCMultRecTrackHist1->SetBinEdges(0,fBinsPtCorr);\r
+ fMCMultRecTrackHist1->SetBinEdges(1,fBinsEta);\r
fMCMultRecTrackHist1->GetAxis(0)->SetTitle("mcPt (GeV/c)");\r
fMCMultRecTrackHist1->GetAxis(1)->SetTitle("mcEta");\r
fMCMultRecTrackHist1->GetAxis(2)->SetTitle("pid");\r
\r
//nClust:chi2PerClust:pt:eta:phi\r
- Int_t binsRecTrackHist2[5]={160,100,ptNbins,etaNbins,90};\r
+ Int_t binsRecTrackHist2[5]={160,100,fPtNbins,fEtaNbins,90};\r
Double_t minRecTrackHist2[5]={0., 0., 0, -1.5, 0.};\r
Double_t maxRecRecTrackHist2[5]={160.,10., 16, 1.5, 2.*TMath::Pi()};\r
\r
fRecTrackHist2 = new THnSparseF("fRecTrackHist2","nClust:chi2PerClust:pt:eta:phi",5,binsRecTrackHist2,minRecTrackHist2,maxRecRecTrackHist2);\r
- fRecTrackHist2->SetBinEdges(2,binsPt);\r
- fRecTrackHist2->SetBinEdges(3,binsEta);\r
+ fRecTrackHist2->SetBinEdges(2,fBinsPt);\r
+ fRecTrackHist2->SetBinEdges(3,fBinsEta);\r
fRecTrackHist2->GetAxis(0)->SetTitle("nClust");\r
fRecTrackHist2->GetAxis(1)->SetTitle("chi2PerClust");\r
fRecTrackHist2->GetAxis(2)->SetTitle("p_{T} (GeV/c)");\r
//\r
// 0-multiplicity bin for LHC background correction\r
//\r
+ /* bin0 done in the train\r
if( GetAnalysisMode() == AlidNdPtHelper::kTPCITS || \r
GetAnalysisMode() == AlidNdPtHelper::kTPCTrackSPDvtx || \r
GetAnalysisMode() == AlidNdPtHelper::kTPCTrackSPDvtxUpdate || \r
GetAnalysisMode() == AlidNdPtHelper::kTPCITSHybridTrackSPDvtx || \r
GetAnalysisMode() == AlidNdPtHelper::kTPCITSHybridTrackSPDvtxDCArPt ) \r
{\r
- physicsSelection->SetBin0CallbackViaPointer(&AlidNdPtAnalysis::IsBinZeroTrackSPDvtx);\r
+ physicsSelection->SetBin0CallbackViaPointer(&AlidNdPtAnalysis::IsBinZeroTrackSPDvtx);\r
} else {\r
- physicsSelection->SetBin0CallbackViaPointer(&AlidNdPtAnalysis::IsBinZeroSPDvtx);\r
+ physicsSelection->SetBin0CallbackViaPointer(&AlidNdPtAnalysis::IsBinZeroSPDvtx);\r
}\r
+ */\r
}\r
}\r
\r
if(!stack) return;\r
\r
Int_t label = TMath::Abs(esdTrack->GetLabel()); \r
- if(label > stack->GetNtrack()) \r
- {\r
- printf("esd track label: %d, stack->GetNtrack(): %d\n", label, stack->GetNtrack());\r
- return;\r
- }\r
+ //if(label == 0) return;\r
\r
TParticle* particle = stack->Particle(label);\r
if(!particle) return;\r