Added TPC BetheBloch parameters for PbPb (preliminary) (Rossella)

diff --git a/PWG3/vertexingHF/AliAODPidHF.cxx b/PWG3/vertexingHF/AliAODPidHF.cxx
index 515debe296846ec6c406e071fa3bac0a90e32210..ac1518d255516d0222edebcd4f9bfd8346f7eaa3 100644 (file)
--- a/PWG3/vertexingHF/AliAODPidHF.cxx
+++ b/PWG3/vertexingHF/AliAODPidHF.cxx
@@ -49,7 +49,8 @@ AliAODPidHF::AliAODPidHF():
   fMatch(0),
   fCompat(kFALSE),
   fMC(kFALSE),
-  fOnePad(kFALSE)
+  fOnePad(kFALSE),
+  fPbPb(kFALSE)
 {
  //
  // Default constructor
@@ -95,7 +96,8 @@ AliAODPidHF::AliAODPidHF(const AliAODPidHF& pid) :
   fMatch(pid.fMatch),
   fCompat(pid.fCompat),
   fMC(pid.fMC),
-  fOnePad(pid.fOnePad)
+  fOnePad(pid.fOnePad),
+  fPbPb(pid.fPbPb)
   {
   
   for(Int_t i=0;i<5;i++){
@@ -661,15 +663,23 @@ void AliAODPidHF::SetBetheBloch(AliTPCPIDResponse tpcResp) const{
 
  }else{
 
-  if(fOnePad){
-   
+  if(fOnePad) {
+
    alephParameters[0] =1.34490e+00/50; 
    alephParameters[1] = 2.69455e+01; 
    alephParameters[2] = TMath::Exp(-2.97552+01); 
    alephParameters[3] = 2.35339e+00; 
    alephParameters[4] = 5.98079e+00;
-  }else{
 
+  } else if(fPbPb) {
+
+   alephParameters[0] = 1.25202/50.; 
+   alephParameters[1] = 2.74992e+01; 
+   alephParameters[2] = TMath::Exp(-3.31517e+01); 
+   alephParameters[3] = 2.46246; 
+   alephParameters[4] = 6.78938;
+
+  } else {
     alephParameters[0] = 0.0283086/0.97;
     alephParameters[1] = 2.63394e+01;
     alephParameters[2] = 5.04114e-11;

diff --git a/PWG3/vertexingHF/AliAODPidHF.h b/PWG3/vertexingHF/AliAODPidHF.h
index 7e21d81dee3c6d28dc60eccbff4955704cca463f..850acf4494141e0869fac1a4786c15ff92ba4449 100644 (file)
--- a/PWG3/vertexingHF/AliAODPidHF.h
+++ b/PWG3/vertexingHF/AliAODPidHF.h
@@ -44,6 +44,7 @@ class AliAODPidHF : public AliAODPid{
  void SetCompat(Bool_t comp){fCompat=comp;return;}
  void SetMC(Bool_t mc){fMC=mc;return;}
  void SetOnePad(Bool_t onepad){fOnePad=onepad;return;}
+ void SetPbPb(Bool_t pbpb){fPbPb=pbpb;return;}
  
  //Getters
  Double_t GetSigma(Int_t idet) const{return fnSigma[idet];}
@@ -58,6 +59,8 @@ class AliAODPidHF : public AliAODPid{
  Int_t GetMatch() const{return fMatch;}
  Bool_t GetCompat() const{return fCompat;}
  Bool_t GetMC() const{return fMC;}
+ Bool_t GetOnePad() const{return fOnePad;}
+ Bool_t GetPbPb() const{return fPbPb;}
 
  Int_t RawSignalPID (AliAODTrack *track, TString detector) const;
  Bool_t IsKaonRaw (AliAODTrack *track, TString detector) const;
@@ -103,6 +106,7 @@ class AliAODPidHF : public AliAODPid{
  Bool_t fCompat; // compatibility region : useful only if fMatch=1
  Bool_t fMC; // MC(kTRUE) or real data (kFALSE, default option)
  Bool_t fOnePad; //  real data with one pad clusters 
+ Bool_t fPbPb; //  real data PbPb 
  
 
 

diff --git a/PWG3/vertexingHF/AliRDHFCuts.cxx b/PWG3/vertexingHF/AliRDHFCuts.cxx
index c50e946c8b0ff58326ef4e6e39a6dd459e1c0f29..926b88cdcc04fa549f3d2b776f26ce746e3747f2 100644 (file)
--- a/PWG3/vertexingHF/AliRDHFCuts.cxx
+++ b/PWG3/vertexingHF/AliRDHFCuts.cxx
@@ -185,7 +185,10 @@ Bool_t AliRDHFCuts::IsEventSelected(AliVEvent *event) {
 
   // switch to settings for 1-pad cls in TPC
   if(fPidHF) {
-    if(event->GetRunNumber()>121693) fPidHF->SetOnePad(kTRUE);
+    if(event->GetRunNumber()>121693 && event->GetRunNumber()<136851) 
+      fPidHF->SetOnePad(kTRUE);
+    if(event->GetRunNumber()>=136851 && event->GetRunNumber()<=139517) 
+      fPidHF->SetPbPb(kTRUE);
   }
 
   if(fOptPileup==kRejectPileupEvent){