Fixing bug #59311
authorfca <fca@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 25 Nov 2009 09:40:44 +0000 (09:40 +0000)
committerfca <fca@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 25 Nov 2009 09:40:44 +0000 (09:40 +0000)
ITS/AliITSRecoParam.cxx
ITS/AliITSV0Finder.cxx
STEER/AliESDV0Params.cxx
STEER/AliESDV0Params.h

index 9866a1f..51180b1 100644 (file)
@@ -231,6 +231,8 @@ fESDV0Params(NULL)
   fESDV0Params->SetMinNormDistForb3(1.0);
   fESDV0Params->SetMinNormDistForb4(4.0);
   fESDV0Params->SetMinNormDistForb5(5.0);
+  fESDV0Params->SetMinNormDistForbProt(2.0);
+  fESDV0Params->SetMaxPidProbPionForb(0.5);
 
   fESDV0Params->SetMinRTPCdensity(40.);
   fESDV0Params->SetMaxRTPCdensity0(110.);
@@ -484,6 +486,8 @@ AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
   param->GetESDV0Params()->SetMinPABestConst(0.99);
   param->GetESDV0Params()->SetMinNormDistForbTgl0(1.);
   param->GetESDV0Params()->SetMinNormDistForb1(2.);
+  param->GetESDV0Params()->SetMinNormDistForbProt(1.);
+  param->GetESDV0Params()->SetMaxPidProbPionForb(0.7);
   param->GetESDV0Params()->SetLikelihood01Cut(0.3);
   param->GetESDV0Params()->SetLikelihood1Cut(0.35);
   param->GetESDV0Params()->SetCombinedCut(0.4);
index 113200d..ecaba15 100644 (file)
@@ -204,6 +204,8 @@ void AliITSV0Finder::FindV02(AliESDEvent *event,
   const Float_t kMinNormDistForb3= AliITSReconstructor::GetRecoParam()->GetESDV0Params()->GetMinNormDistForb3();\r
   const Float_t kMinNormDistForb4= AliITSReconstructor::GetRecoParam()->GetESDV0Params()->GetMinNormDistForb4();\r
   const Float_t kMinNormDistForb5= AliITSReconstructor::GetRecoParam()->GetESDV0Params()->GetMinNormDistForb5();\r
+  const Float_t kMinNormDistForbProt= AliITSReconstructor::GetRecoParam()->GetESDV0Params()->GetMinNormDistForbProt();\r
+  const Float_t kMaxPidProbPionForb= AliITSReconstructor::GetRecoParam()->GetESDV0Params()->GetMaxPidProbPionForb();\r
 \r
   const Float_t kMinRTPCdensity= AliITSReconstructor::GetRecoParam()->GetESDV0Params()->GetMinRTPCdensity();\r
   const Float_t kMaxRTPCdensity0= AliITSReconstructor::GetRecoParam()->GetESDV0Params()->GetMaxRTPCdensity0();\r
@@ -466,10 +468,15 @@ void AliITSV0Finder::FindV02(AliESDEvent *event,
       normdist[itsindex]*=-1;\r
     }\r
     if (isProton){\r
-      if (normdist[itsindex]>2) forbidden[itsindex]=kFALSE;    \r
+      if (normdist[itsindex]>kMinNormDistForbProt) forbidden[itsindex]=kFALSE; \r
       normdist[itsindex]*=-1;\r
     }\r
 \r
+    // We allow all tracks that are not pions\r
+    if( (pid[1]+pid[2])< kMaxPidProbPionForb ){\r
+      forbidden[itsindex]=kFALSE;\r
+    }\r
+\r
     //\r
     // Causality cuts in TPC volume\r
     //\r
index 9517a4d..54130dd 100644 (file)
@@ -91,7 +91,8 @@ AliESDV0Params::AliESDV0Params() :
   fkMinNormDistForb3(1.0),
   fkMinNormDistForb4(4.0),
   fkMinNormDistForb5(5.0),
-  
+  fkMinNormDistForbProt(2.0),
+  fkMaxPidProbPionForb(0.5),
   fkMinRTPCdensity(40.),
   fkMaxRTPCdensity0(110.),
   fkMaxRTPCdensity10(120.),
index da84582..f497fcb 100644 (file)
@@ -43,6 +43,9 @@ class AliESDV0Params : public TObject{
   void SetMinNormDistForb3(Float_t kMinNormDistForb3=1.0){fkMinNormDistForb3=kMinNormDistForb3;}
   void SetMinNormDistForb4(Float_t kMinNormDistForb4=4.0){fkMinNormDistForb4=kMinNormDistForb4;}
   void SetMinNormDistForb5(Float_t kMinNormDistForb5=5.0){fkMinNormDistForb5=kMinNormDistForb5;}
+  void SetMinNormDistForbProt(Float_t kMinNormDistForbProt=2.0){fkMinNormDistForbProt=kMinNormDistForbProt;}
+  void SetMaxPidProbPionForb(Float_t kMaxPidProbPionForb=0.5){fkMaxPidProbPionForb=kMaxPidProbPionForb;}
+
 
   void SetMinRTPCdensity(Float_t kMinRTPCdensity=40.0){fkMinRTPCdensity=kMinRTPCdensity;}
   void SetMaxRTPCdensity0( Float_t kMaxRTPCdensity0=110.0){fkMaxRTPCdensity0=kMaxRTPCdensity0;}
@@ -97,6 +100,8 @@ class AliESDV0Params : public TObject{
   Float_t GetMinNormDistForb3() const {return fkMinNormDistForb3;}
   Float_t GetMinNormDistForb4() const {return fkMinNormDistForb4;}
   Float_t GetMinNormDistForb5() const {return fkMinNormDistForb5;}
+  Float_t GetMinNormDistForbProt() const {return fkMinNormDistForbProt;}
+  Float_t GetMaxPidProbPionForb() const {return fkMaxPidProbPionForb;}
 
   Float_t GetMinRTPCdensity() const  {return fkMinRTPCdensity;}
   Float_t GetMaxRTPCdensity0() const {return fkMaxRTPCdensity0;}
@@ -181,6 +186,9 @@ class AliESDV0Params : public TObject{
   Float_t fkMinNormDistForb3;       // Minimum  normalize distance Forbid cond 3  
   Float_t fkMinNormDistForb4;       // Minimum  normalize distance Forbid cond 4 
   Float_t fkMinNormDistForb5;       // Minimum  normalize distance Forbid cond 5
+  Float_t fkMinNormDistForbProt;    // Minimum  normalize distance to not Forbid Proton
+  Float_t fkMaxPidProbPionForb;     // Max pid prob to decleare not pion  
+
 
   Float_t fkMinRTPCdensity;         // Minimum R TPC density cond  
   Float_t fkMaxRTPCdensity0;        // Maximum R TPC density cond  
@@ -214,7 +222,7 @@ class AliESDV0Params : public TObject{
 
 
 
-  ClassDef(AliESDV0Params,3)      // ESD V0 vertex - error and likelihood parameterization constant
+  ClassDef(AliESDV0Params,4)      // ESD V0 vertex - error and likelihood parameterization constant
 };