Fixed a bug in the Digit reader, moved clusterfinder to initialize from OCDB, Added...

[u/mrichter/AliRoot.git] / ITS / AliITSRecoParam.h

diff --git a/ITS/AliITSRecoParam.h b/ITS/AliITSRecoParam.h
index c33218d6141e6568c45586954c4fc9d822a6e250..93b18cc88933810eddc021796c02ca3aaa332a12 100644 (file)
--- a/ITS/AliITSRecoParam.h
+++ b/ITS/AliITSRecoParam.h
@@ -15,6 +15,7 @@
 
 #include "AliDetectorRecoParam.h"
 #include "AliITSgeomTGeo.h"
+#include "AliESDV0Params.h"
 
 class AliITSRecoParam : public AliDetectorRecoParam
 {
@@ -197,8 +198,11 @@ class AliITSRecoParam : public AliDetectorRecoParam
   void   SetClusterMisalErrorY(Float_t e0,Float_t e1,Float_t e2,Float_t e3,Float_t e4,Float_t e5) { fClusterMisalErrorY[0]=e0; fClusterMisalErrorY[1]=e1; fClusterMisalErrorY[2]=e2; fClusterMisalErrorY[3]=e3; fClusterMisalErrorY[4]=e4; fClusterMisalErrorY[5]=e5; return; }
   void   SetClusterMisalErrorZ(Float_t e0,Float_t e1,Float_t e2,Float_t e3,Float_t e4,Float_t e5) { fClusterMisalErrorZ[0]=e0; fClusterMisalErrorZ[1]=e1; fClusterMisalErrorZ[2]=e2; fClusterMisalErrorZ[3]=e3; fClusterMisalErrorZ[4]=e4; fClusterMisalErrorZ[5]=e5; return; }
   void   SetClusterMisalError(Float_t err=0.) { SetClusterMisalErrorY(err,err,err,err,err,err); SetClusterMisalErrorZ(err,err,err,err,err,err); }
-  Float_t GetClusterMisalErrorY(Int_t i) const { return fClusterMisalErrorY[i]; }
-  Float_t GetClusterMisalErrorZ(Int_t i) const { return fClusterMisalErrorZ[i]; }
+  void   SetClusterMisalErrorYBOn(Float_t e0,Float_t e1,Float_t e2,Float_t e3,Float_t e4,Float_t e5) { fClusterMisalErrorYBOn[0]=e0; fClusterMisalErrorYBOn[1]=e1; fClusterMisalErrorYBOn[2]=e2; fClusterMisalErrorYBOn[3]=e3; fClusterMisalErrorYBOn[4]=e4; fClusterMisalErrorYBOn[5]=e5; return; }
+  void   SetClusterMisalErrorZBOn(Float_t e0,Float_t e1,Float_t e2,Float_t e3,Float_t e4,Float_t e5) { fClusterMisalErrorZBOn[0]=e0; fClusterMisalErrorZBOn[1]=e1; fClusterMisalErrorZBOn[2]=e2; fClusterMisalErrorZBOn[3]=e3; fClusterMisalErrorZBOn[4]=e4; fClusterMisalErrorZBOn[5]=e5; return; }
+  void   SetClusterMisalErrorBOn(Float_t err=0.) { SetClusterMisalErrorYBOn(err,err,err,err,err,err); SetClusterMisalErrorZBOn(err,err,err,err,err,err); }
+  Float_t GetClusterMisalErrorY(Int_t i,Double_t b=0.) const { return (TMath::Abs(b)<0.0001 ? fClusterMisalErrorY[i] : fClusterMisalErrorYBOn[i]); }
+  Float_t GetClusterMisalErrorZ(Int_t i,Double_t b=0.) const { return (TMath::Abs(b)<0.0001 ? fClusterMisalErrorZ[i] : fClusterMisalErrorZBOn[i]); }
 
   void   SetUseAmplitudeInfo(Bool_t use=kTRUE) { for(Int_t i=0;i<AliITSgeomTGeo::kNLayers;i++) fUseAmplitudeInfo[i]=use; return; }
   void   SetUseAmplitudeInfo(Int_t ilay,Bool_t use) { fUseAmplitudeInfo[ilay]=use; return; }
@@ -216,8 +220,12 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Bool_t  GetBkgTrackletsPlaneEff() const {return fBkgTrackletsPlaneEff;}
   void    SetTrackleterPhiWindowL1(Float_t w=0.10) {fTrackleterPhiWindowL1=w; return;}
   Float_t GetTrackleterPhiWindowL1() const {return fTrackleterPhiWindowL1;}
+  void    SetTrackleterPhiWindowL2(Float_t w=0.07) {fTrackleterPhiWindowL2=w; return;}
+  Float_t GetTrackleterPhiWindowL2() const {return fTrackleterPhiWindowL2;}
   void    SetTrackleterZetaWindowL1(Float_t w=0.6) {fTrackleterZetaWindowL1=w; return;}
   Float_t GetTrackleterZetaWindowL1() const {return fTrackleterZetaWindowL1;}
+  void    SetTrackleterZetaWindowL2(Float_t w=0.40) {fTrackleterZetaWindowL2=w; return;}
+  Float_t GetTrackleterZetaWindowL2() const {return fTrackleterZetaWindowL2;}
   void    SetUpdateOncePerEventPlaneEff(Bool_t use=kTRUE) {fUpdateOncePerEventPlaneEff=use; return;}
   Bool_t  GetUpdateOncePerEventPlaneEff() const {return fUpdateOncePerEventPlaneEff;}
   void    SetMinContVtxPlaneEff(Int_t n=3) {fMinContVtxPlaneEff=n; return;}
@@ -320,12 +328,12 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Bool_t GetUseCosmicRunShiftsSSD() const { return fUseCosmicRunShiftsSSD; }
 
   // SPD Tracklets (D. Elia)
-  void    SetTrackleterOnlyOneTrackletPerC2(Bool_t use= kTRUE) {fTrackleterOnlyOneTrackletPerC2=use; return; }
-  Bool_t  GetTrackleterOnlyOneTrackletPerC2() const { return fTrackleterOnlyOneTrackletPerC2; }
   void    SetTrackleterPhiWindow(Float_t w=0.08) {fTrackleterPhiWindow=w;}
-  void    SetTrackleterZetaWindow(Float_t w=1.) {fTrackleterZetaWindow=w;}
+  void    SetTrackleterThetaWindow(Float_t w=0.025) {fTrackleterThetaWindow=w;}
+  void    SetTrackleterPhiShift(Float_t w=0.0045) {fTrackleterPhiShift=w;}
   Float_t GetTrackleterPhiWindow() const {return fTrackleterPhiWindow;}
-  Float_t GetTrackleterZetaWindow() const {return fTrackleterZetaWindow;}
+  Float_t GetTrackleterThetaWindow() const {return fTrackleterThetaWindow;}
+  Float_t GetTrackleterPhiShift() const {return fTrackleterPhiShift;}
   void    SetTrackleterRemoveClustersFromOverlaps(Bool_t use=kTRUE) { fTrackleterRemoveClustersFromOverlaps=use; return; }
   Bool_t  GetTrackleterRemoveClustersFromOverlaps() const { return fTrackleterRemoveClustersFromOverlaps; }
   void    SetTrackleterPhiOverlapCut(Float_t w=0.005) {fTrackleterPhiOverlapCut=w;}
@@ -367,6 +375,7 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Float_t GetAlignFilterMaxPt() const {return fAlignFilterMaxPt;}          
   Bool_t  GetAlignFilterFillQANtuples() const {return fAlignFilterFillQANtuples;}     
 
+  AliESDV0Params *GetESDV0Params() const {return fESDV0Params;}
 
 
   enum {fgkMaxClusterPerLayer=70000}; //7000*10;   // max clusters per layer
@@ -494,6 +503,8 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Int_t fClusterErrorsParam; // parametrization for cluster errors (MI), see AliITSRecoParam::GetError()
   Float_t fClusterMisalErrorY[AliITSgeomTGeo::kNLayers]; // [cm] additional error on cluster Y pos. due to misalignment (MI,SA)
   Float_t fClusterMisalErrorZ[AliITSgeomTGeo::kNLayers]; // [cm] additional error on cluster Z pos. due to misalignment (MI,SA)
+  Float_t fClusterMisalErrorYBOn[AliITSgeomTGeo::kNLayers]; // [cm] additional error on cluster Y pos. due to misalignment (MI,SA)
+  Float_t fClusterMisalErrorZBOn[AliITSgeomTGeo::kNLayers]; // [cm] additional error on cluster Z pos. due to misalignment (MI,SA)
 
   Bool_t fUseAmplitudeInfo[AliITSgeomTGeo::kNLayers]; // use cluster charge in cluster-track matching (SDD,SSD) (MI)
 
@@ -504,7 +515,9 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Bool_t fMCTrackletsPlaneEff; // flag to enable the use of MC info for corrections (SPD PlaneEff using tracklets)
   Bool_t fBkgTrackletsPlaneEff; // flag to evaluate background instead of normal use (SPD PlaneEff using tracklets)
   Float_t fTrackleterPhiWindowL1; // Search window in phi for inner layer (1) (SPD PlaneEff using tracklets)
+  Float_t fTrackleterPhiWindowL2; // Search window in phi for outer layer (2) (SPD PlaneEff using tracklets)
   Float_t fTrackleterZetaWindowL1; // Search window in zeta for inner layer (1) (SPD PlaneEff using tracklets)
+  Float_t fTrackleterZetaWindowL2; // Search window in zeta for outer layer (2) (SPD PlaneEff using tracklets)
   Bool_t fUpdateOncePerEventPlaneEff; // option to update chip efficiency once/event (to avoid doubles)
   Int_t  fMinContVtxPlaneEff; // min number of contributors to ESD vtx for SPD PlaneEff using tracklets
   Int_t  fIPlanePlaneEff; // index of the plane (in the range [-1,5]) to study the efficiency (-1 ->Tracklets)
@@ -552,9 +565,9 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Bool_t fUseChargeMatchingInClusterFinderSSD; // SSD
 
   // SPD Tracklets (D. Elia)
-  Bool_t  fTrackleterOnlyOneTrackletPerC2;         // Allow only one tracklet per cluster in the outer layer
   Float_t fTrackleterPhiWindow;                    // Search window in phi
-  Float_t fTrackleterZetaWindow;                   // Search window in eta
+  Float_t fTrackleterThetaWindow;                   // Search window in theta
+  Float_t fTrackleterPhiShift;                     // Phi shift reference value (at 0.5 T) 
   Bool_t  fTrackleterRemoveClustersFromOverlaps;   // Option to skip clusters in the overlaps
   Float_t fTrackleterPhiOverlapCut;                // Fiducial window in phi for overlap cut
   Float_t fTrackleterZetaOverlapCut;               // Fiducial window in eta for overlap cut
@@ -586,8 +599,13 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Float_t fAlignFilterMaxPt;              // max pt
   Bool_t  fAlignFilterFillQANtuples;      // fill QA ntuples  
 
+ private:
+  AliESDV0Params * fESDV0Params;  // declare the AliESDV0Params to be able to used in AliITSV0Finder
 
-  ClassDef(AliITSRecoParam,23) // ITS reco parameters
+  AliITSRecoParam(const AliITSRecoParam & param);
+  AliITSRecoParam & operator=(const AliITSRecoParam &param);
+
+  ClassDef(AliITSRecoParam,26) // ITS reco parameters
 };
 
 #endif