TPC BB param for LHC11a MCs (Francesco)
authordainese <dainese@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 9 May 2011 14:06:53 +0000 (14:06 +0000)
committerdainese <dainese@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 9 May 2011 14:06:53 +0000 (14:06 +0000)
PWG3/vertexingHF/AliAODPidHF.cxx
PWG3/vertexingHF/AliAODPidHF.h
PWG3/vertexingHF/AliRDHFCuts.cxx

index fbd7328..55193c2 100644 (file)
@@ -55,6 +55,7 @@ AliAODPidHF::AliAODPidHF():
   fnSigmaCompat(),
   fMC(kFALSE),
   fOnePad(kFALSE),
+  fMCLowEn2011(kFALSE),
   fPbPb(kFALSE)
 {
  //
@@ -109,6 +110,7 @@ AliAODPidHF::AliAODPidHF(const AliAODPidHF& pid) :
   fnSigmaCompat(pid.fnSigmaCompat),
   fMC(pid.fMC),
   fOnePad(pid.fOnePad),
+  fMCLowEn2011(pid.fMCLowEn2011),
   fPbPb(pid.fPbPb)
   {
   
@@ -705,13 +707,19 @@ void AliAODPidHF::SetBetheBloch(AliTPCPIDResponse &tpcResp) const{
      alephParameters[4] = 4.62254e+00;
 
    } else {  // pp MC
-
-     alephParameters[0] = 2.15898/50.;
-     alephParameters[1] = 1.75295e+01;
-     alephParameters[2] = 3.40030e-09;
-     alephParameters[3] = 1.96178e+00;
-     alephParameters[4] = 3.91720e+00;
-
+     if(fMCLowEn2011){
+       alephParameters[0]=0.0207667;
+       alephParameters[1]=29.9936;
+       alephParameters[2]=3.87866e-11;
+       alephParameters[3]=2.17291;
+       alephParameters[4]=7.1623;
+     }else{
+       alephParameters[0] = 2.15898/50.;
+       alephParameters[1] = 1.75295e+01;
+       alephParameters[2] = 3.40030e-09;
+       alephParameters[3] = 1.96178e+00;
+       alephParameters[4] = 3.91720e+00;
+     }
    }
 
  } else { // Real Data
index d9d5f8f..c87f683 100644 (file)
@@ -47,6 +47,7 @@ class AliAODPidHF : public AliAODPid{
  void SetMatch(Int_t match){fMatch=match;return;}
  void SetCompat(Bool_t comp){fCompat=comp;return;}
  void SetMC(Bool_t mc){fMC=mc;return;}
+ void SetMClowenpp2011(Bool_t mc){fMCLowEn2011=mc;return;}
  void SetOnePad(Bool_t onepad){fOnePad=onepad;return;}
  void SetPbPb(Bool_t pbpb){fPbPb=pbpb;return;}
  void SetPCompatTOF(Double_t pTOF){fPCompatTOF=pTOF;return;}
@@ -67,6 +68,7 @@ class AliAODPidHF : public AliAODPid{
  Bool_t GetCompat() const{return fCompat;}
  Bool_t GetMC() const{return fMC;}
  Bool_t GetOnePad() const{return fOnePad;}
+ Bool_t GetMCLowEn2011() const {return fMCLowEn2011;}
  Bool_t GetPbPb() const{return fPbPb;}
  Double_t GetPCompatTOF() const{return fPCompatTOF;}
  Double_t GetnSigmaCompatTPC() const{return fnSigmaCompat[0];}
@@ -120,11 +122,12 @@ class AliAODPidHF : public AliAODPid{
  Double_t *fnSigmaCompat; //[fnNSigmaCompat]  0: n sigma for TPC compatibility band, 1: for TOF  
  Bool_t fMC; // MC(kTRUE) or real data (kFALSE, default option)
  Bool_t fOnePad; //  real data with one pad clusters 
+ Bool_t fMCLowEn2011; //  MC for low energy MC
  Bool_t fPbPb; //  real data PbPb 
  
 
 
- ClassDef(AliAODPidHF,9) // AliAODPid for heavy flavor PID
+ ClassDef(AliAODPidHF,10) // AliAODPid for heavy flavor PID
 
 };
 
index 4dd0309..c55f814 100644 (file)
@@ -247,6 +247,8 @@ Bool_t AliRDHFCuts::IsEventSelected(AliVEvent *event) {
     if(event->GetRunNumber()>=136851 && event->GetRunNumber()<=139517) fPidHF->SetPbPb(kTRUE);
     // MC
     if(isMC) fPidHF->SetMC(kTRUE);
+    if(isMC && (event->GetRunNumber()>=146686 && event->GetRunNumber()<=146860))
+      fPidHF->SetMClowenpp2011(kTRUE);
   }