MC strange particles selection added
authorjotwinow <jotwinow@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 11 Nov 2010 18:10:01 +0000 (18:10 +0000)
committerjotwinow <jotwinow@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 11 Nov 2010 18:10:01 +0000 (18:10 +0000)
PWG0/dNdPt/AlidNdPtCutAnalysisPbPb.cxx

index 4521028..89b6d5d 100644 (file)
@@ -148,7 +148,7 @@ void AlidNdPtCutAnalysisPbPb::Init(){
   // 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
@@ -164,7 +164,7 @@ void AlidNdPtCutAnalysisPbPb::Init(){
   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
@@ -377,6 +377,7 @@ void AlidNdPtCutAnalysisPbPb::FillHistograms(AliESDtrack *const esdTrack, AliSta
   // 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
@@ -387,13 +388,23 @@ void AlidNdPtCutAnalysisPbPb::FillHistograms(AliESDtrack *const esdTrack, AliSta
     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