// THnSparse track histograms\r
//\r
\r
- //nClust:chi2PerClust:nClust/nFindableClust:DCAy:DCAz:eta:phi:pt:kinkIdx:isPrim:polarity\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., -1.,-1.,-1.5, 0., ptMin, -1., 0., 0.};\r
Double_t maxRecMCTrackHist[11]={160.,10.,1.2, 1.,1.,1.5, 2.*TMath::Pi(), ptMax, 2., 2., 2.};\r
fRecMCTrackHist->GetAxis(5)->SetTitle("#eta");\r
fRecMCTrackHist->GetAxis(6)->SetTitle("#phi (rad)");\r
fRecMCTrackHist->GetAxis(7)->SetTitle("p_{T} (GeV/c)");\r
- fRecMCTrackHist->GetAxis(8)->SetTitle("kinkIdx"); // 0 - no kink, -1 - kink mother, 1 - kink daugther \r
+ fRecMCTrackHist->GetAxis(8)->SetTitle("hasStrangeMother");\r
fRecMCTrackHist->GetAxis(9)->SetTitle("isPrim");\r
fRecMCTrackHist->GetAxis(10)->SetTitle("polarity");\r
fRecMCTrackHist->Sumw2();\r
// Fill rec vs MC information\r
//\r
Bool_t isPrim = kTRUE;\r
+ Bool_t hasStrangeMother = kTRUE;\r
\r
if(IsUseMCInfo()) {\r
if(!stack) 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
+ TParticle* mother = 0; \r
+ \r
+ Int_t motherLabel = particle->GetMother(0); \r
+ if(motherLabel>0) mother = stack->Particle(motherLabel); \r
+ if(mother) motherPdg = TMath::Abs(mother->GetPdgCode()); // take abs for visualisation only \r
+ Int_t mech = particle->GetUniqueID(); // production mechanism \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,kinkIdx,isPrim, polarity }; \r
+ Double_t vRecMCTrackHist[11] = { nClust,chi2PerCluster,clustPerFindClust,b[0],b[1],eta,phi,pt,hasStrangeMother,isPrim, polarity }; \r
fRecMCTrackHist->Fill(vRecMCTrackHist);\r
}\r
\r