ITS/UPGRADE/AliITSURecoParam.h

   1 #ifndef ALIITSURECOPARAM_H
   2 #define ALIITSURECOPARAM_H
   3 /* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id: AliITSURecoParam.h 57215 2012-06-17 14:47:08Z masera $ */
   7
   8 ///////////////////////////////////////////////////////////////////////////////
   9 //                                                                           //
  10 // Class with ITS reconstruction parameters                                  //
  11 //                                                                           //
  12 ///////////////////////////////////////////////////////////////////////////////
  13
  14 class AliITSUTrackCond;
  15
  16 #include <TObjArray.h>
  17 #include "AliDetectorRecoParam.h"
  18
  19 class AliITSURecoParam : public AliDetectorRecoParam
  20 {
  21  public:
  22   enum {kClus2Tracks,kPropBack,kRefitInw,kNTrackingPhases};   // tracking phases
  23   //
  24   AliITSURecoParam();
  25   AliITSURecoParam(Int_t nLr);
  26   virtual ~AliITSURecoParam();
  27
  28   static AliITSURecoParam *GetLowFluxParam();// make reco parameters for low flux env.
  29   static AliITSURecoParam *GetHighFluxParam();// make reco parameters for high flux env.
  30   static AliITSURecoParam *GetCosmicTestParam();// special setting for cosmic
  31
  32   Double_t    GetMaxDforV0dghtrForProlongation() const {return fMaxDforV0dghtrForProlongation;}
  33   Double_t    GetMaxDForProlongation()           const {return fMaxDForProlongation;}
  34   Double_t    GetMaxDZForProlongation()          const {return fMaxDZForProlongation;}
  35   Double_t    GetMinPtForProlongation()          const {return fMinPtForProlongation;}
  36   Double_t    GetTanLorentzAngle(Int_t lr)       const;
  37   Double_t    GetSigmaY2(Int_t lr)               const;
  38   Double_t    GetSigmaZ2(Int_t lr)               const;
  39   Bool_t      GetAllowDiagonalClusterization(Int_t lr) const;
  40   //
  41   Double_t    GetTPCITSWallRMin()                const {return fTPCITSWallRMin;}
  42   Double_t    GetTPCITSWallRMax()                const {return fTPCITSWallRMax;}
  43   Double_t    GetTPCITSWallZSpanH()              const {return fTPCITSWallZSpanH;}
  44   Double_t    GetTPCITSWallMaxStep()             const {return fTPCITSWallMaxStep;}
  45   TObjArray*  GetTrackingConditions()            const {return (TObjArray*)&fTrackingConditions;}
  46   Int_t       GetNTrackingConditions()           const {return fTrackingConditions.GetEntriesFast();}
  47   AliITSUTrackCond* GetTrackingCondition(Int_t i) const {return (AliITSUTrackCond*)fTrackingConditions[i];}
  48   Bool_t      GetUseMatLUT(Int_t step)           const {return (step<0||step>=kNTrackingPhases) ? kFALSE:fUseMatLUT[step];}
  49   Int_t GetTracker() const {return fTracker;}
  50   Bool_t GetSAonly() const {return fSAonly;}
  51   Int_t GetMaxROCycle() const {return fMaxROCycle;}
  52   //
  53   void        SetNLayers(Int_t n);
  54   void        SetTanLorentzAngle(Int_t lr, Double_t v);
  55   void        SetSigmaY2(Int_t lr, Double_t v);
  56   void        SetSigmaZ2(Int_t lr, Double_t v);
  57   void        SetAllowDiagonalClusterization(Int_t lr, Bool_t v);
  58   //
  59   void        SetMaxDforV0dghtrForProlongation(Double_t v)            {fMaxDforV0dghtrForProlongation = v;}
  60   void        SetMaxDForProlongation(Double_t v)                      {fMaxDForProlongation = v;}
  61   void        SetMaxDZForProlongation(Double_t v)                     {fMaxDZForProlongation = v;}
  62   void        SetMinPtForProlongation(Double_t v)                     {fMinPtForProlongation = v;}
  63   //
  64   void        SetTPCITSWallRMin(double v)                             {fTPCITSWallRMin = v;}
  65   void        SetTPCITSWallRMax(double v)                             {fTPCITSWallRMax = v;}
  66   void        SetTPCITSWallZSpanH(double v)                           {fTPCITSWallZSpanH = v;}
  67   void        SetTPCITSWallMaxStep(double v)                          {fTPCITSWallMaxStep = v;}
  68   //
  69   void        SetUseMatLUT(Int_t step, Bool_t v)                      {if (step>=0&&step<kNTrackingPhases) fUseMatLUT[step]=v;}
  70   void        AddTrackingCondition(AliITSUTrackCond* cond);
  71   virtual void Print(Option_t *opt="")  const;
  72   void SetTracker(Int_t opt) {fTracker=opt;}
  73   void SetSAonly(Bool_t sa=kTRUE) {fSAonly=sa;}
  74   void SetMaxROCycle(Int_t ro) {fMaxROCycle=ro;}
  75   //
  76  protected:
  77   Int_t          fNLayers;          // number of layers
  78   //
  79   Double_t       fMaxDforV0dghtrForProlongation; // max. rphi imp. par. cut for V0 daughter
  80   Double_t       fMaxDForProlongation; // max. rphi imp. par. cut
  81   Double_t       fMaxDZForProlongation; // max. 3D imp. par. cut
  82   Double_t       fMinPtForProlongation; // min. pt cut
  83   //
  84   Double_t       fTPCITSWallRMin;       // minR
  85   Double_t       fTPCITSWallRMax;       // maxR
  86   Double_t       fTPCITSWallZSpanH;     // half Z span
  87   Double_t       fTPCITSWallMaxStep;    // max tracking step
  88   //
  89   Bool_t         fUseMatLUT[kNTrackingPhases]; // allow usage of material lookup table at given tracking step
  90   Bool_t*        fAllowDiagonalClusterization; //[fNLayers] allow clusters of pixels with common corners only
  91   Double_t*      fTanLorentzAngle;  //[fNLayers] Optional Lorentz angle for each layer
  92   Double_t*      fSigmaY2;          //[fNLayers] addition to road width^2
  93   Double_t*      fSigmaZ2;          //[fNLayers] addition to road width^2
  94   //
  95   TObjArray      fTrackingConditions; // array of tracking conditions for different iterations
  96   //
  97   static const Double_t fgkMaxDforV0dghtrForProlongation;      // default
  98   static const Double_t fgkMaxDForProlongation;                // default
  99   static const Double_t fgkMaxDZForProlongation;               // default
 100   static const Double_t fgkMinPtForProlongation;               // default
 101   static const Double_t fgkTanLorentzAngle;                    // default
 102   static const Double_t fgkSigmaRoadY;                         // default
 103   static const Double_t fgkSigmaRoadZ;                         // default
 104   //
 105   // hardwired params for TPC-ITS border layer
 106   static const Double_t fgkTPCITSWallRMin;                     // fiducial R min
 107   static const Double_t fgkTPCITSWallRMax;                     // fiducial R max
 108   static const Double_t fgkTPCITSWallZSpanH;                   // half Z span
 109   static const Double_t fgkTPCITSWallMaxStep;                  // max tracking step
 110   //
 111   static const Bool_t   fgkUseMatLUT[kNTrackingPhases];        // def. mat. LUT usage
 112   static const Bool_t   fgkAllowDiagonalClusterization;        // clusters of pixels with common corners
 113   //
 114   Int_t fTracker;  // 0 - global, 1 - cooked, 2 - CA
 115   Bool_t fSAonly;  // kTRUE, if ITSU standalone tracking only
 116   Int_t fMaxROCycle; // max RO cycle accepted by the cluster finder
 117  private:
 118   AliITSURecoParam(const AliITSURecoParam & param);
 119   AliITSURecoParam & operator=(const AliITSURecoParam &param);
 120
 121   ClassDef(AliITSURecoParam,6) // ITS reco parameters
 122 };
 123
 124 //_____________________________________________________________________________
 125 inline Double_t AliITSURecoParam::GetTanLorentzAngle(Int_t lr) const
 126 {
 127   // get tg of Lorentz Angle for the layer
 128   return (lr<fNLayers)&&fTanLorentzAngle ? fTanLorentzAngle[lr]:0;
 129 }
 130
 131 //_____________________________________________________________________________
 132 inline Double_t AliITSURecoParam::GetSigmaY2(Int_t lr) const
 133 {
 134   // get tg of Lorentz Angle for the layer
 135   return (lr<fNLayers)&&fSigmaY2 ? fSigmaY2[lr]:fgkSigmaRoadY*fgkSigmaRoadY; //0;
 136 }
 137
 138 //_____________________________________________________________________________
 139 inline Double_t AliITSURecoParam::GetSigmaZ2(Int_t lr) const
 140 {
 141   // get tg of Lorentz Angle for the layer
 142   return (lr<fNLayers)&&fSigmaZ2 ? fSigmaZ2[lr]:fgkSigmaRoadZ*fgkSigmaRoadZ;//0;
 143 }
 144
 145 //_____________________________________________________________________________
 146 inline Bool_t AliITSURecoParam::GetAllowDiagonalClusterization(Int_t lr) const
 147 {
 148   // are diagonal clusters permitted
 149   return (lr<fNLayers)&&fAllowDiagonalClusterization ? fAllowDiagonalClusterization[lr]:fgkAllowDiagonalClusterization;
 150 }
 151
 152 #endif
 153
 154