add setting of sigma range for PID in tracking

sevdokim · Sep 7, 2016 · e48af86 · e48af86
1 parent b2244f6
commit e48af86
Show file tree

Hide file tree

Showing 4 changed files with 12 additions and 1 deletion.
diff --git a/STEER/STEERBase/AliPIDResponse.h b/STEER/STEERBase/AliPIDResponse.h
@@ -205,6 +205,9 @@ class AliPIDResponse : public TNamed {
 
   void SetNoTOFmism(Bool_t value=kTRUE){fNoTOFmism=value;};
 
+  void    SetProbabilityRangeNsigma(Float_t range) { fRange = range; }
+  Float_t GetProbabilityRangeNsigma() const        { return fRange;  }
+
 protected:
   AliITSPIDResponse   fITSResponse;    //PID response function of the ITS
   AliTPCPIDResponse   fTPCResponse;    //PID response function of the TPC

diff --git a/TPC/TPCbase/AliTPCParam.cxx b/TPC/TPCbase/AliTPCParam.cxx
@@ -133,6 +133,7 @@ AliTPCParam::AliTPCParam()
 	     fTotalNormFac(0.),
 	     fNoiseNormFac(0.),
 	     fUseGlitchFilter(kTRUE),
+	     fSigmaRangePIDinTracking(15.),
 	     fNominalVoltage(),
 	     fMaxVoltageDeviation(40.),
 	     fMaxDipVoltage(2.),

diff --git a/TPC/TPCbase/AliTPCParam.h b/TPC/TPCbase/AliTPCParam.h
@@ -224,6 +224,7 @@ class AliTPCParam : public AliDetectorParam {
   void  SetADCSat(Int_t adcsat)    {  fADCSat  = adcsat;}
   void  SetADCDynRange(Float_t adcdynrange) {fADCDynRange = adcdynrange;}
   void  SetUseGlitchFilter(Bool_t use) {fUseGlitchFilter=use;}
+  void  SetSigmaRangePIDinTracking(Float_t range) { fSigmaRangePIDinTracking = range; }
   //
   // High voltage parameters
   //
@@ -373,6 +374,8 @@ class AliTPCParam : public AliDetectorParam {
   Float_t  GetTotalNormFac() const {return fTotalNormFac;}
   Float_t  GetNoiseNormFac() const {return fNoiseNormFac;}
   Bool_t   GetUseGlitchFilter() const {return fUseGlitchFilter;}
+  Float_t  GetSigmaRangePIDinTracking() const { return fSigmaRangePIDinTracking; }
+
   //
   // High voltage parameters
   //
@@ -519,6 +522,7 @@ protected :
   Float_t fTotalNormFac;    ///< full normalisation factor - calculated
   Float_t fNoiseNormFac;    ///< normalisation factor to transform noise in electron to ADC channel
   Bool_t  fUseGlitchFilter; ///< use the glitch filter to suppress single digits above threshold
+  Float_t fSigmaRangePIDinTracking; ///< maximum sigma range allowed for PID probability range calculation used in tracking
 
   //---------------------------------------------------------------------
   // High voltage parameters
@@ -554,7 +558,7 @@ protected :
   void CleanGeoMatrices();
 
   /// \cond CLASSIMP
-  ClassDef(AliTPCParam,12)  //parameter  object for set:TPC
+  ClassDef(AliTPCParam,13)  //parameter  object for set:TPC
   /// \endcond
 };
 

diff --git a/TPC/TPCrec/AliTPCReconstructor.cxx b/TPC/TPCrec/AliTPCReconstructor.cxx
@@ -298,10 +298,13 @@ void AliTPCReconstructor::GetPidSettings(AliESDpid *esdPID)
     AliTPCParam* param = AliTPCcalibDB::Instance()->GetParameters();
     if (param) {
       TVectorD *paramBB=param->GetBetheBlochParameters();
+      const Float_t maxnSigmaRange = param->GetSigmaRangePIDinTracking();
       if (paramBB){
         esdPID->GetTPCResponse().SetBetheBlochParameters((*paramBB)(0),(*paramBB)(1),(*paramBB)(2),(*paramBB)(3),(*paramBB)(4));
         AliInfo(Form("Setting BB parameters from OCDB (AliTPCParam): %.2g, %.2g, %.2g, %.2g, %.2g",
                      (*paramBB)(0),(*paramBB)(1),(*paramBB)(2),(*paramBB)(3),(*paramBB)(4)));
+        AliInfoF("Setting max sigma PID range to %.2f", maxnSigmaRange);
+        esdPID->SetProbabilityRangeNsigma(maxnSigmaRange);
       } else {
         AliError("Couldn't get BB parameters from OCDB, the old default values will be used instead");
       }