//\r
// THnSparse track histograms\r
//\r
+ //nClust:chi2PerClust:nClust/nFindableClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromMaterial:isPrim:charge\r
+ Int_t binsRecMCTrackHist[12]={160,20,40,100,100,30,90,ptNbins, 2, 2, 2, 3};\r
+ Double_t minRecMCTrackHist[12]={0., 0., 0., -10.,-10.,-1.5, 0., ptMin, 0., 0., 0., -1.};\r
+ Double_t maxRecMCTrackHist[12]={160.,10.,1.2, 10.,10.,1.5, 2.*TMath::Pi(), ptMax, 2., 2., 2., 2.};\r
\r
- //nClust:chi2PerClust:nClust/nFindableClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isPrim:polarity\r
- Int_t binsRecMCTrackHist[11]={160,80,80,100,100,90,90,ptNbins, 3, 2, 2};\r
- Double_t minRecMCTrackHist[11]={0., 0., 0., -5.,-5.,-1.5, 0., ptMin, -1., 0., 0.};\r
- Double_t maxRecMCTrackHist[11]={160.,10.,1.2, 5.,5.,1.5, 2.*TMath::Pi(), ptMax, 2., 2., 2.};\r
-\r
- fRecMCTrackHist = new THnSparseF("fRecMCTrackHist","nClust:chi2PerClust:nClust/nFindableClust:DCAy:DCAz:eta:phi:pt:kinkIdx:isPrim:polarity",11,binsRecMCTrackHist,minRecMCTrackHist,maxRecMCTrackHist);\r
+ fRecMCTrackHist = new THnSparseF("fRecMCTrackHist","nClust:chi2PerClust:nClust/nFindableClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromMaterial:isPrim:charge",12,binsRecMCTrackHist,minRecMCTrackHist,maxRecMCTrackHist);\r
fRecMCTrackHist->SetBinEdges(7,binsPt);\r
\r
fRecMCTrackHist->GetAxis(0)->SetTitle("nClust");\r
fRecMCTrackHist->GetAxis(6)->SetTitle("#phi (rad)");\r
fRecMCTrackHist->GetAxis(7)->SetTitle("p_{T} (GeV/c)");\r
fRecMCTrackHist->GetAxis(8)->SetTitle("hasStrangeMother");\r
- fRecMCTrackHist->GetAxis(9)->SetTitle("isPrim");\r
- fRecMCTrackHist->GetAxis(10)->SetTitle("polarity");\r
+ fRecMCTrackHist->GetAxis(9)->SetTitle("isFromMaterial");\r
+ fRecMCTrackHist->GetAxis(10)->SetTitle("isPrim");\r
+ fRecMCTrackHist->GetAxis(11)->SetTitle("charge");\r
fRecMCTrackHist->Sumw2();\r
\r
// init output folder\r
vtxESD = esdEvent->GetPrimaryVertexTPC();\r
}\r
}\r
-\r
Bool_t isEventOK = evtCuts->AcceptEvent(esdEvent,mcEvent,vtxESD) && isRecVertex; \r
//printf("isEventOK %d, isEventTriggered %d \n",isEventOK, isEventTriggered);\r
//printf("GetAnalysisMode() %d \n",GetAnalysisMode());\r
//\r
Bool_t isPrim = kTRUE;\r
Bool_t hasStrangeMother = kFALSE;\r
+ Bool_t isFromMaterial = kFALSE;\r
\r
if(IsUseMCInfo()) {\r
if(!stack) return;\r
if(particle->GetPDG() && particle->GetPDG()->Charge()==0.) return;\r
isPrim = stack->IsPhysicalPrimary(label);\r
\r
- //if(isPrim && pt > 1.5 && kinkIdx == -1) printf("nClust %d \n", nClust);\r
- \r
// check whether has stange mother\r
//\r
Int_t motherPdg = -1; \r
if( (motherPdg == 3122) || (motherPdg == -3122) || (motherPdg == 310)) // lambda, antilambda, k0s\r
{\r
if( (mech == 4) || (mech == 5) ) hasStrangeMother = kTRUE;\r
+ } \r
+ else {\r
+ //if(isPrim==0 && mech == 13) \r
+ //printf("mech %d \n", mech);\r
+ if(!isPrim) isFromMaterial = kTRUE; \r
}\r
}\r
+\r
\r
// fill histo\r
- Int_t polarity = -2;\r
- if (esdTrack->Charge() < 0.) polarity = 0; \r
- else polarity = 1; \r
- Double_t vRecMCTrackHist[11] = { nClust,chi2PerCluster,clustPerFindClust,b[0],b[1],eta,phi,pt,hasStrangeMother,isPrim, polarity }; \r
+ Int_t charge = esdTrack->Charge();\r
+ Double_t vRecMCTrackHist[12] = { nClust,chi2PerCluster,clustPerFindClust,b[0],b[1],eta,phi,pt,hasStrangeMother,isFromMaterial,isPrim,charge }; \r
fRecMCTrackHist->Fill(vRecMCTrackHist);\r
}\r
\r