Use of PID info (ESD combined PID and/or MC PID) for Inv Mass histos (A.Maire)

Priors for Xi (0,0,2,0,1) and for Omega(0,0,1,1,1).

diff --git a/PWG2/SPECTRA/AliAnalysisTaskCheckCascade.cxx b/PWG2/SPECTRA/AliAnalysisTaskCheckCascade.cxx
index fd9296ef88ebcd06531bb8c302c7b5805e6d740f..455180adc14dd656e442e4c6aebf66a8778b4ebf 100644 (file)
--- a/PWG2/SPECTRA/AliAnalysisTaskCheckCascade.cxx
+++ b/PWG2/SPECTRA/AliAnalysisTaskCheckCascade.cxx
@@ -35,7 +35,7 @@ class AliAODVertex;
 class AliESDv0;
 class AliAODv0;
 
-#include <iostream>
+#include <Riostream.h>
 
 #include "TList.h"
 #include "TH1.h"
@@ -84,6 +84,8 @@ AliAnalysisTaskCheckCascade::AliAnalysisTaskCheckCascade()
 
     fHistMassXiMinus(0), fHistMassXiPlus(0),
     fHistMassOmegaMinus(0), fHistMassOmegaPlus(0),
+    fHistMassWithCombPIDXiMinus(0), fHistMassWithCombPIDXiPlus(0),
+    fHistMassWithCombPIDOmegaMinus(0), fHistMassWithCombPIDOmegaPlus(0),
 
     fHistXiTransvMom(0),    fHistXiTotMom(0),
     fHistBachTransvMom(0),   fHistBachTotMom(0),
@@ -140,6 +142,8 @@ AliAnalysisTaskCheckCascade::AliAnalysisTaskCheckCascade(const char *name)
 
     fHistMassXiMinus(0), fHistMassXiPlus(0),
     fHistMassOmegaMinus(0), fHistMassOmegaPlus(0),
+    fHistMassWithCombPIDXiMinus(0), fHistMassWithCombPIDXiPlus(0),
+    fHistMassWithCombPIDOmegaMinus(0), fHistMassWithCombPIDOmegaPlus(0),
 
     fHistXiTransvMom(0),    fHistXiTotMom(0),
     fHistBachTransvMom(0),   fHistBachTotMom(0),
@@ -359,7 +363,7 @@ if (! fHistDcaNegToPrimVertexXi) {
 
 
        // - Effective mass histos for cascades.
-  
+// By cascade hyp  
 if (! fHistMassXiMinus) {
     fHistMassXiMinus = new TH1F("fHistMassXiMinus","#Xi^{-} candidates;M( #Lambda , #pi^{-} ) (GeV/c^{2});Counts", 200,1.2,2.0);
     fListHistCascade->Add(fHistMassXiMinus);
@@ -380,6 +384,27 @@ if (! fHistMassOmegaPlus) {
     fListHistCascade->Add(fHistMassOmegaPlus);
 }
 
+// By cascade hyp + bachelor PID
+if (! fHistMassWithCombPIDXiMinus) {
+    fHistMassWithCombPIDXiMinus = new TH1F("fHistMassWithCombPIDXiMinus","#Xi^{-} candidates, with Bach. comb. PID;M( #Lambda , #pi^{-} ) (GeV/c^{2});Counts", 200,1.2,2.0);
+    fListHistCascade->Add(fHistMassWithCombPIDXiMinus);
+}
+  
+if (! fHistMassWithCombPIDXiPlus) {
+    fHistMassWithCombPIDXiPlus = new TH1F("fHistMassWithCombPIDXiPlus","#Xi^{+} candidates, with Bach. comb. PID;M( #bar{#Lambda}^{0} , #pi^{+} ) (GeV/c^{2});Counts",200,1.2,2.0);
+    fListHistCascade->Add(fHistMassWithCombPIDXiPlus);
+}
+
+if (! fHistMassWithCombPIDOmegaMinus) {
+       fHistMassWithCombPIDOmegaMinus = new TH1F("fHistMassWithCombPIDOmegaMinus","#Omega^{-} candidates, with Bach. comb. PID;M( #Lambda , K^{-} ) (GeV/c^{2});Counts", 250,1.5,2.5);
+    fListHistCascade->Add(fHistMassWithCombPIDOmegaMinus);
+}
+ 
+if (! fHistMassWithCombPIDOmegaPlus) {
+       fHistMassWithCombPIDOmegaPlus = new TH1F("fHistMassWithCombPIDOmegaPlus","#Omega^{+} candidates, with Bach. comb. PID;M( #bar{#Lambda}^{0} , K^{+} ) (GeV/c^{2});Counts", 250,1.5,2.5);
+    fListHistCascade->Add(fHistMassWithCombPIDOmegaPlus);
+}
+
 
 
        // - Complements for QA
@@ -598,9 +623,13 @@ void AliAnalysisTaskCheckCascade::UserExec(Option_t *)
   Double_t lInvMassXiPlus     = 0.;
   Double_t lInvMassOmegaMinus = 0.;
   Double_t lInvMassOmegaPlus  = 0.;
+  
+       // - 4th part of initialisation : PID treatment of the bachelor track
+  Bool_t   lIsBachelorKaon    = kFALSE;
+  Bool_t   lIsBachelorPion    = kFALSE;
 
 
-       // - 4th part of initialisation : extra info for QA
+       // - 5th part of initialisation : extra info for QA
   
   Double_t lXiMomX       = 0., lXiMomY = 0., lXiMomZ = 0.;
   Double_t lXiTransvMom  = 0. ;
@@ -795,7 +824,6 @@ void AliAnalysisTaskCheckCascade::UserExec(Option_t *)
                //-------------
 
        
-       
        if( bachTrackXi->Charge() < 0 ) {
                lV0quality = 0.;
                xi->ChangeMassHypothesis(lV0quality , 3312);    
@@ -811,8 +839,8 @@ void AliAnalysisTaskCheckCascade::UserExec(Option_t *)
                                        
                lV0quality = 0.;
                xi->ChangeMassHypothesis(lV0quality , 3312);    // Back to default hyp.
-               
-       }
+       }// end if negative bachelor
+       
        
        if( bachTrackXi->Charge() >  0 ){
                lV0quality = 0.;
@@ -829,12 +857,41 @@ void AliAnalysisTaskCheckCascade::UserExec(Option_t *)
                
                lV0quality = 0.;
                xi->ChangeMassHypothesis(lV0quality , -3312);   // Back to "default" hyp.
-       }
+       }// end if positive bachelor
+       
+       
+       
+               // - II.Step 5 : PID on the bachelor
+               //-------------
+       
+       // Reasonable guess for the priors for the cascade track sample
+       Double_t lPriorsGuessXi[5]    = {0.0, 0.0, 2, 0, 1};
+       Double_t lPriorsGuessOmega[5] = {0.0, 0.0, 1, 1, 1};
+       AliPID pidXi;           pidXi.SetPriors(    lPriorsGuessXi    );
+       AliPID pidOmega;        pidOmega.SetPriors( lPriorsGuessOmega );
+       
+       if( bachTrackXi->IsOn(AliESDtrack::kESDpid) ){  // Combined PID exists
+               Double_t r[10]; bachTrackXi->GetESDpid(r);
+               pidXi.SetProbabilities(r);
+               pidOmega.SetProbabilities(r);
+               // Check if the bachelor track is a pion
+               Double_t ppion = pidXi.GetProbability(AliPID::kPion);
+               if (ppion > pidXi.GetProbability(AliPID::kElectron) &&
+                   ppion > pidXi.GetProbability(AliPID::kMuon)     &&
+                   ppion > pidXi.GetProbability(AliPID::kKaon)     &&
+                   ppion > pidXi.GetProbability(AliPID::kProton)   )     lIsBachelorPion = kTRUE;
+               // Check if the bachelor track is a kaon
+               Double_t pkaon = pidOmega.GetProbability(AliPID::kKaon);
+               if (pkaon > pidOmega.GetProbability(AliPID::kElectron) &&
+                   pkaon > pidOmega.GetProbability(AliPID::kMuon)     &&
+                   pkaon > pidOmega.GetProbability(AliPID::kPion)     &&
+                   pkaon > pidOmega.GetProbability(AliPID::kProton)   )  lIsBachelorKaon = kTRUE;
                
+       }// end if bachelor track with existing combined PID
        
-                                                       
-
-               // - II.Step 5 : extra info for QA (ESD)
+               
+               
+               // - II.Step 6 : extra info for QA (ESD)
                // miscellaneous pieces onf info that may help regarding data quality assessment.
                //-------------
 
@@ -952,7 +1009,13 @@ void AliAnalysisTaskCheckCascade::UserExec(Option_t *)
        if( lChargeXi > 0 )             lInvMassOmegaPlus       = xi->MassOmega();
 
        
-               // - II.Step 5 : extra info for QA (AOD)
+               // - II.Step 5 : PID on the bachelor
+               //-------------
+       
+       
+       // To be developed
+       
+               // - II.Step 6 : extra info for QA (AOD)
                // miscellaneous pieces onf info that may help regarding data quality assessment.
                //-------------
 
@@ -1027,11 +1090,21 @@ void AliAnalysisTaskCheckCascade::UserExec(Option_t *)
        fHistDcaNegToPrimVertexXi       ->Fill( lDcaNegToPrimVertexXi      );
                
        
-       // - III.Step 4
-       if( lChargeXi < 0 )     fHistMassXiMinus        ->Fill( lInvMassXiMinus );
-       if( lChargeXi > 0 )     fHistMassXiPlus         ->Fill( lInvMassXiPlus );
-       if( lChargeXi < 0 )     fHistMassOmegaMinus     ->Fill( lInvMassOmegaMinus );
-       if( lChargeXi > 0 )     fHistMassOmegaPlus      ->Fill( lInvMassOmegaPlus );
+       // - III.Step 4+5
+       if( lChargeXi < 0 ){
+                                       fHistMassXiMinus               ->Fill( lInvMassXiMinus    );
+                                       fHistMassOmegaMinus            ->Fill( lInvMassOmegaMinus );
+               if(lIsBachelorPion)     fHistMassWithCombPIDXiMinus    ->Fill( lInvMassXiMinus    );
+               if(lIsBachelorKaon)     fHistMassWithCombPIDOmegaMinus ->Fill( lInvMassOmegaMinus );
+       }
+       
+       if( lChargeXi > 0 ){
+                                       fHistMassXiPlus                ->Fill( lInvMassXiPlus     );
+                                       fHistMassOmegaPlus             ->Fill( lInvMassOmegaPlus  );
+               if(lIsBachelorPion)     fHistMassWithCombPIDXiPlus     ->Fill( lInvMassXiPlus     );
+               if(lIsBachelorKaon)     fHistMassWithCombPIDOmegaPlus  ->Fill( lInvMassOmegaPlus  );
+       }
+       
 
 //     if( bachTrackXi->Charge() < 0 ) fHistMassXiMinus        ->Fill( lInvMassXiMinus    );
 //     if( bachTrackXi->Charge() > 0 ) fHistMassXiPlus         ->Fill( lInvMassXiPlus     );

diff --git a/PWG2/SPECTRA/AliAnalysisTaskCheckCascade.h b/PWG2/SPECTRA/AliAnalysisTaskCheckCascade.h
index 8e57ff8ee97291e6ceb2a41c57846d32e0282eab..be31d8e778300cf8bebd701202c569a077423e35 100644 (file)
--- a/PWG2/SPECTRA/AliAnalysisTaskCheckCascade.h
+++ b/PWG2/SPECTRA/AliAnalysisTaskCheckCascade.h
@@ -91,6 +91,11 @@ class AliAnalysisTaskCheckCascade : public AliAnalysisTaskSE {
        TH1F    *fHistMassXiPlus;                       //! reconstructed cascade effective mass, under Xi+ hyp.
        TH1F    *fHistMassOmegaMinus;                   //! reconstructed cascade effective mass, under Omega- hyp.
        TH1F    *fHistMassOmegaPlus;                    //! reconstructed cascade effective mass, under Omega+ hyp.
+       
+       TH1F    *fHistMassWithCombPIDXiMinus;           //! reconstructed Xi- effective mass, with bach. comb PID
+       TH1F    *fHistMassWithCombPIDXiPlus;            //! reconstructed Xi+ effective mass, with bach. comb PID
+       TH1F    *fHistMassWithCombPIDOmegaMinus;        //! reconstructed Omega- effective mass, with bach. comb PID
+       TH1F    *fHistMassWithCombPIDOmegaPlus;         //! reconstructed Omega+ effective mass, with bach. comb PID
 
        // - Complements for QA
        TH1F    *fHistXiTransvMom;                      //! Xi transverse momentum 

diff --git a/PWG2/SPECTRA/AliAnalysisTaskCheckPerformanceCascade.cxx b/PWG2/SPECTRA/AliAnalysisTaskCheckPerformanceCascade.cxx
index fea83821534dffe9100635fc7e8413194f5a9148..557d92c743dc10af5ae5129050cba94c279f7533 100644 (file)
--- a/PWG2/SPECTRA/AliAnalysisTaskCheckPerformanceCascade.cxx
+++ b/PWG2/SPECTRA/AliAnalysisTaskCheckPerformanceCascade.cxx
@@ -150,6 +150,14 @@ fListHistCascade(0),
     fHistMassXiPlus(0),
     fHistMassOmegaMinus(0),
     fHistMassOmegaPlus(0),
+    
+       // - Effective mass histos with combined PID
+    fHistMassWithCombPIDXiMinus(0), fHistMassWithCombPIDXiPlus(0),
+    fHistMassWithCombPIDOmegaMinus(0), fHistMassWithCombPIDOmegaPlus(0),
+    
+       // - Effective mass histos with perfect MC PID on the bachelor
+    fHistMassWithMcPIDXiMinus(0), fHistMassWithMcPIDXiPlus(0),
+    fHistMassWithMcPIDOmegaMinus(0), fHistMassWithMcPIDOmegaPlus(0),
        
        // - Effective mass histos for the cascade candidates associated with MC
     fHistAsMCMassXiMinus(0),           
@@ -308,6 +316,14 @@ fListHistCascade(0),
     fHistMassOmegaMinus(0),
     fHistMassOmegaPlus(0),
        
+       // - Effective mass histos with combined PID
+    fHistMassWithCombPIDXiMinus(0), fHistMassWithCombPIDXiPlus(0),
+    fHistMassWithCombPIDOmegaMinus(0), fHistMassWithCombPIDOmegaPlus(0),
+    
+       // - Effective mass histos with perfect MC PID on the bachelor
+    fHistMassWithMcPIDXiMinus(0), fHistMassWithMcPIDXiPlus(0),
+    fHistMassWithMcPIDOmegaMinus(0), fHistMassWithMcPIDOmegaPlus(0),
+    
        // - Effective mass histos for the cascade candidates associated with MC
     fHistAsMCMassXiMinus(0),           
     fHistAsMCMassXiPlus(0),            
@@ -767,6 +783,53 @@ void AliAnalysisTaskCheckPerformanceCascade::UserCreateOutputObjects()
          fListHistCascade->Add(fHistMassOmegaPlus);
   }
   
+  
+  
+               // - Effective mass histos with combined PID
+  
+  if (! fHistMassWithCombPIDXiMinus) {
+    fHistMassWithCombPIDXiMinus = new TH1F("fHistMassWithCombPIDXiMinus","#Xi^{-} candidates, with Bach. comb. PID;M( #Lambda , #pi^{-} ) (GeV/c^{2});Counts", 200,1.2,2.0);
+    fListHistCascade->Add(fHistMassWithCombPIDXiMinus);
+  }
+  
+  if (! fHistMassWithCombPIDXiPlus) {
+    fHistMassWithCombPIDXiPlus = new TH1F("fHistMassWithCombPIDXiPlus","#Xi^{+} candidates, with Bach. comb. PID;M( #bar{#Lambda}^{0} , #pi^{+} ) (GeV/c^{2});Counts",200,1.2,2.0);
+    fListHistCascade->Add(fHistMassWithCombPIDXiPlus);
+  }
+
+  if (! fHistMassWithCombPIDOmegaMinus) {
+       fHistMassWithCombPIDOmegaMinus = new TH1F("fHistMassWithCombPIDOmegaMinus","#Omega^{-} candidates, with Bach. comb. PID;M( #Lambda , K^{-} ) (GeV/c^{2});Counts", 250,1.5,2.5);
+    fListHistCascade->Add(fHistMassWithCombPIDOmegaMinus);
+  }
+ 
+  if (! fHistMassWithCombPIDOmegaPlus) {
+       fHistMassWithCombPIDOmegaPlus = new TH1F("fHistMassWithCombPIDOmegaPlus","#Omega^{+} candidates, with Bach. comb. PID;M( #bar{#Lambda}^{0} , K^{+} ) (GeV/c^{2});Counts", 250,1.5,2.5);
+    fListHistCascade->Add(fHistMassWithCombPIDOmegaPlus);
+  }
+  
+               // - Effective mass histos with perfect MC PID on the bachelor
+  
+  if (! fHistMassWithMcPIDXiMinus) {
+    fHistMassWithMcPIDXiMinus = new TH1F("fHistMassWithMcPIDXiMinus","#Xi^{-} candidates, with Bach. MC PID;M( #Lambda , #pi^{-} ) (GeV/c^{2});Counts", 200,1.2,2.0);
+    fListHistCascade->Add(fHistMassWithMcPIDXiMinus);
+  }
+  
+  if (! fHistMassWithMcPIDXiPlus) {
+    fHistMassWithMcPIDXiPlus = new TH1F("fHistMassWithMcPIDXiPlus","#Xi^{+} candidates, with Bach. MC PID;M( #bar{#Lambda}^{0} , #pi^{+} ) (GeV/c^{2});Counts",200,1.2,2.0);
+    fListHistCascade->Add(fHistMassWithMcPIDXiPlus);
+  }
+
+  if (! fHistMassWithMcPIDOmegaMinus) {
+       fHistMassWithMcPIDOmegaMinus = new TH1F("fHistMassWithMcPIDOmegaMinus","#Omega^{-} candidates, with Bach. MC PID;M( #Lambda , K^{-} ) (GeV/c^{2});Counts", 250,1.5,2.5);
+    fListHistCascade->Add(fHistMassWithMcPIDOmegaMinus);
+  }
+ 
+  if (! fHistMassWithMcPIDOmegaPlus) {
+       fHistMassWithMcPIDOmegaPlus = new TH1F("fHistMassWithMcPIDOmegaPlus","#Omega^{+} candidates, with Bach. MC PID;M( #bar{#Lambda}^{0} , K^{+} ) (GeV/c^{2});Counts", 250,1.5,2.5);
+    fListHistCascade->Add(fHistMassWithMcPIDOmegaPlus);
+  }
+  
+  
                // - Effective mass histos for cascades candidates ASSOCIATED with MC.
   
   if (! fHistAsMCMassXiMinus) {
@@ -1441,6 +1504,63 @@ for (Int_t iXi = 0; iXi < ncascades; iXi++)
        if( lChargeXi > 0 )     fHistMassOmegaPlus      ->Fill( lInvMassOmegaPlus );
        
        
+       // - Step 4 : PID info
+       //-------------
+       
+       Bool_t   lIsBachelorKaon    = kFALSE;
+       Bool_t   lIsBachelorPion    = kFALSE;
+       
+       // Reasonable guess for the priors for the cascade track sample
+       Double_t lPriorsGuessXi[5]    = {0.0, 0.0, 2, 0, 1};
+       Double_t lPriorsGuessOmega[5] = {0.0, 0.0, 1, 1, 1};
+       AliPID pidXi;           pidXi.SetPriors(    lPriorsGuessXi    );
+       AliPID pidOmega;        pidOmega.SetPriors( lPriorsGuessOmega );
+       
+       if( bachTrackXi->IsOn(AliESDtrack::kESDpid) ){  // Combined PID exists
+               Double_t r[10]; bachTrackXi->GetESDpid(r);
+               pidXi.SetProbabilities(r);
+               pidOmega.SetProbabilities(r);
+               // Check if the bachelor track is a pion
+               Double_t ppion = pidXi.GetProbability(AliPID::kPion);
+               if (ppion > pidXi.GetProbability(AliPID::kElectron) &&
+                   ppion > pidXi.GetProbability(AliPID::kMuon)     &&
+                   ppion > pidXi.GetProbability(AliPID::kKaon)     &&
+                   ppion > pidXi.GetProbability(AliPID::kProton)   )     lIsBachelorPion = kTRUE;
+               // Check if the bachelor track is a kaon
+               Double_t pkaon = pidOmega.GetProbability(AliPID::kKaon);
+               if (pkaon > pidOmega.GetProbability(AliPID::kElectron) &&
+                   pkaon > pidOmega.GetProbability(AliPID::kMuon)     &&
+                   pkaon > pidOmega.GetProbability(AliPID::kPion)     &&
+                   pkaon > pidOmega.GetProbability(AliPID::kProton)   )  lIsBachelorKaon = kTRUE;
+               
+       }// end if bachelor track with existing combined PID
+       
+                  
+       if( lChargeXi < 0 && lIsBachelorPion )    fHistMassWithCombPIDXiMinus     ->Fill( lInvMassXiMinus );
+       if( lChargeXi > 0 && lIsBachelorPion )    fHistMassWithCombPIDXiPlus      ->Fill( lInvMassXiPlus );
+       if( lChargeXi < 0 && lIsBachelorKaon )    fHistMassWithCombPIDOmegaMinus  ->Fill( lInvMassOmegaMinus );
+       if( lChargeXi > 0 && lIsBachelorKaon )    fHistMassWithCombPIDOmegaPlus   ->Fill( lInvMassOmegaPlus );
+         
+       
+       Bool_t   lIsBachelorMCPiMinus  = kFALSE;
+       Bool_t   lIsBachelorMCPiPlus   = kFALSE;
+       Bool_t   lIsBachelorMCKMinus   = kFALSE;
+       Bool_t   lIsBachelorMCKPlus    = kFALSE;        
+       
+       Int_t      lblBachForPID  = (Int_t) TMath::Abs( bachTrackXi->GetLabel() );
+       TParticle* mcBachForPID   = lMCstack->Particle( lblBachForPID );
+       
+       if( mcBachForPID->GetPdgCode() == -211) lIsBachelorMCPiMinus = kTRUE;
+       if( mcBachForPID->GetPdgCode() ==  211) lIsBachelorMCPiPlus  = kTRUE;
+       if( mcBachForPID->GetPdgCode() == -321) lIsBachelorMCKMinus  = kTRUE;
+       if( mcBachForPID->GetPdgCode() ==  321) lIsBachelorMCKPlus   = kTRUE;
+       
+       if( lChargeXi < 0 && lIsBachelorMCPiMinus )    fHistMassWithMcPIDXiMinus     ->Fill( lInvMassXiMinus );
+       if( lChargeXi > 0 && lIsBachelorMCPiPlus  )    fHistMassWithMcPIDXiPlus      ->Fill( lInvMassXiPlus );
+       if( lChargeXi < 0 && lIsBachelorMCKMinus  )    fHistMassWithMcPIDOmegaMinus  ->Fill( lInvMassOmegaMinus );
+       if( lChargeXi > 0 && lIsBachelorMCKPlus   )    fHistMassWithMcPIDOmegaPlus   ->Fill( lInvMassOmegaPlus );
+       
+
        
        
        // - Step 4 : MC association
@@ -1471,7 +1591,7 @@ for (Int_t iXi = 0; iXi < ncascades; iXi++)
        Int_t lblMotherPosV0Dghter = mcPosV0Dghter->GetFirstMother() ; 
        Int_t lblMotherNegV0Dghter = mcNegV0Dghter->GetFirstMother();
        
-               if(lblMotherPosV0Dghter != lblMotherNegV0Dghter) continue; // same mother
+               if( lblMotherPosV0Dghter != lblMotherNegV0Dghter) continue; // same mother
                if( lblMotherPosV0Dghter < 0 ) continue; // mother != primary (!= -1)
                if( lblMotherNegV0Dghter < 0 ) continue;
                                        

diff --git a/PWG2/SPECTRA/AliAnalysisTaskCheckPerformanceCascade.h b/PWG2/SPECTRA/AliAnalysisTaskCheckPerformanceCascade.h
index 33ef021a4de0385232dee33123ba7b86934a0f6e..af4f92eab463631f9d867c5497fb65b915cef241 100644 (file)
--- a/PWG2/SPECTRA/AliAnalysisTaskCheckPerformanceCascade.h
+++ b/PWG2/SPECTRA/AliAnalysisTaskCheckPerformanceCascade.h
@@ -161,6 +161,19 @@ class AliAnalysisTaskCheckPerformanceCascade : public AliAnalysisTaskSE {
        TH1F    *fHistMassOmegaMinus;                   //! reconstructed cascade effective mass, under Omega- hyp.
        TH1F    *fHistMassOmegaPlus;                    //! reconstructed cascade effective mass, under Omega+ hyp.
        
+       // - Effective mass histos with reconstruction combined PID
+       TH1F    *fHistMassWithCombPIDXiMinus;           //! reconstructed Xi- effective mass, with bach. comb PID
+       TH1F    *fHistMassWithCombPIDXiPlus;            //! reconstructed Xi+ effective mass, with bach. comb PID
+       TH1F    *fHistMassWithCombPIDOmegaMinus;        //! reconstructed Omega- effective mass, with bach. comb PID
+       TH1F    *fHistMassWithCombPIDOmegaPlus;         //! reconstructed Omega+ effective mass, with bach. comb PID
+
+       // - Effective mass histos with perfect MC PID
+       TH1F    *fHistMassWithMcPIDXiMinus;             //! reconstructed Xi- effective mass, with MC bach. PID
+       TH1F    *fHistMassWithMcPIDXiPlus;              //! reconstructed Xi+ effective mass, with MC bach. PID
+       TH1F    *fHistMassWithMcPIDOmegaMinus;          //! reconstructed Omega- effective mass, with MC bach. PID
+       TH1F    *fHistMassWithMcPIDOmegaPlus;           //! reconstructed Omega+ effective mass, with MC bach. PID
+
+       
        // - Effective mass histos for the cascade candidates associated with MC
        TH1F    *fHistAsMCMassXiMinus;                  //! reconstr. cascade effective mass, under Xi- hyp. for Associated cand.
        TH1F    *fHistAsMCMassXiPlus;                   //! reconstr. cascade effective mass, under Xi+ hyp. for Associated cand.