Possibility to use nSigma in ITS for PID (Annalisa)
authorfprino <fprino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 3 Mar 2013 23:38:11 +0000 (23:38 +0000)
committerfprino <fprino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 3 Mar 2013 23:38:11 +0000 (23:38 +0000)
PWGHF/vertexingHF/AliAODPidHF.cxx
PWGHF/vertexingHF/AliAODPidHF.h

index f3b6987..8386031 100644 (file)
@@ -953,3 +953,39 @@ void AliAODPidHF::SetUpCombinedPID(){
   break;
  }
 }
+//-----------------------------
+
+Int_t AliAODPidHF::GetnSigmaITS(AliAODTrack *track,Int_t species, Double_t &nsigma) const{
+  // get n sigma for ITS
+
+  Double_t nsigmaITS=-999;
+
+  if (!CheckITSPIDStatus(track)) return -1;
+
+  if (fOldPid) {
+    Double_t mom=track->P();
+    AliAODPid *pidObj = track->GetDetPid();
+    Double_t dedx=pidObj->GetITSsignal();
+
+    AliITSPIDResponse itsResponse;
+    AliPID::EParticleType type=AliPID::EParticleType(species);
+    nsigmaITS = TMath::Abs(itsResponse.GetNumberOfSigmas(mom,dedx,type));
+
+  } // old pid
+  else { // new pid
+
+    AliPID::EParticleType type=AliPID::EParticleType(species);
+    nsigmaITS = TMath::Abs(fPidResponse->NumberOfSigmasITS(track,type));
+
+  } //new pid
+
+  nsigma = nsigmaITS;
+
+  return 1;
+
+
+
+
+
+}
+
index d9696d8..58d476d 100644 (file)
@@ -89,6 +89,7 @@ class AliAODPidHF : public AliAODPid{
  
  Int_t GetnSigmaTPC(AliAODTrack *track, Int_t species, Double_t &sigma) const;
  Int_t GetnSigmaTOF(AliAODTrack *track, Int_t species, Double_t &sigma) const;
+ Int_t GetnSigmaITS(AliAODTrack *track, Int_t species, Double_t &sigma) const;
  Double_t GetSigma(Int_t idet) const{return fnSigma[idet];}
  Double_t GetTofSigma() const{return fTOFSigma;}
  //void GetPriors(Double_t *priors) const{priors=fPriors;return;}