TPC/Rec/AliTPCRecoParam.h

   1 #ifndef ALITPCRECOPARAM_H
   2 #define ALITPCRECOPARAM_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 ///////////////////////////////////////////////////////////////////////////////
   7 //                                                                           //
   8 // Class with TPC reconstruction parameters                                  //
   9 //                                                                           //
  10 ///////////////////////////////////////////////////////////////////////////////
  11
  12
  13 #include "AliDetectorRecoParam.h"
  14
  15 class AliTPCRecoParam : public AliDetectorRecoParam
  16 {
  17  public:
  18   AliTPCRecoParam();
  19   virtual ~AliTPCRecoParam();
  20   virtual void Print(const Option_t* option="") const;
  21   static   Bool_t  GetUseTimeCalibration();
  22   static   void    SetUseTimeCalibration(Bool_t useTimeCalibration);
  23
  24   void     SetUseHLTClusters(Int_t useHLTClusters){fUseHLTClusters=useHLTClusters;}
  25   Int_t    GetUseHLTClusters() const {return fUseHLTClusters;}
  26   void     SetUseHLTPreSeeding(Int_t useHLTPreSeeding){fUseHLTPreSeeding=useHLTPreSeeding;}
  27   Int_t    GetUseHLTPreSeeding() const {return fUseHLTPreSeeding;}
  28   void     SetClusterSharing(Bool_t sharing){fBClusterSharing=sharing;}
  29   Bool_t   GetClusterSharing() const {return fBClusterSharing;}
  30   Double_t GetCtgRange() const     { return fCtgRange;}
  31   Double_t GetMaxSnpTracker() const{ return fMaxSnpTracker;}
  32   Double_t GetMaxSnpTrack() const  { return fMaxSnpTrack;}
  33   Bool_t   GetUseOuterDetectors() const { return fUseOuterDetectors;}
  34   void     SetUseOuterDetectors(Bool_t flag)  { fUseOuterDetectors=flag;}
  35   void     SetMaxChi2TPCTRD(Double_t maxChi2){fMaxChi2TPCTRD=maxChi2;}
  36   Double_t GetMaxChi2TPCTRD() const {return fMaxChi2TPCTRD;}
  37   void     SetMaxChi2TPCITS(Double_t maxChi2){fMaxChi2TPCITS=maxChi2;}
  38   Double_t GetMaxChi2TPCITS() const {return fMaxChi2TPCITS;}
  39   Double_t GetCutSharedClusters(Int_t index)const { return fCutSharedClusters[index];}
  40   void  SetCutSharedClusters(Int_t index, Float_t value){ fCutSharedClusters[index]=value;}
  41   Int_t GetClusterMaxRange(Int_t index)const { return fClusterMaxRange[index];}
  42   void     SetClusterMaxRange(Int_t index, Int_t value){ fClusterMaxRange[index]=value;}
  43   //
  44   Bool_t   DumpSignal()     const  { return fDumpSignal;}
  45   void     SetTimeInterval(Int_t first, Int_t last) { fFirstBin=first, fLastBin =last;}
  46   Int_t    GetFirstBin() const     { return fFirstBin;}
  47   Int_t    GetLastBin() const      { return fLastBin;}
  48   void     SetTimeBinRange(Int_t first, Int_t last){ fFirstBin = first; fLastBin = last;}
  49   Bool_t   GetCalcPedestal()       const  { return fBCalcPedestal;}
  50   Bool_t   GetDoUnfold()           const  { return fBDoUnfold;}
  51   void     SetDoUnfold(Bool_t unfold)     { fBDoUnfold = unfold;}
  52   Float_t  GetDumpAmplitudeMin()   const  { return fDumpAmplitudeMin;}
  53   Float_t  GetMaxNoise()           const  { return fMaxNoise;}
  54   //
  55   Int_t    GetUseOnePadCluster()   const  { return fUseOnePadCluster;}
  56   Bool_t   GetUseHLTOnePadCluster()const  { return fUseHLTOnePadCluster;}
  57   Float_t  GetMinMaxCutAbs()       const  { return fMinMaxCutAbs; }
  58   Float_t  GetMinLeftRightCutAbs() const  { return fMinLeftRightCutAbs;}  // minimal amplitude left right - PRF
  59   Float_t  GetMinUpDownCutAbs()    const  { return fMinUpDownCutAbs;}  // minimal amplitude up-down - TRF
  60   Float_t  GetMinMaxCutSigma()       const  { return fMinMaxCutSigma; }
  61   Float_t  GetMinLeftRightCutSigma() const  { return fMinLeftRightCutSigma;}  // minimal amplitude left right - PRF
  62   Float_t  GetMinUpDownCutSigma()    const  { return fMinUpDownCutSigma;}  // minimal amplitude up-down - TRF
  63   //
  64   void SetUseOnePadCluster(Int_t use)      {   fUseOnePadCluster = use;}
  65   void SetUseHLTOnePadCluster(Bool_t use)  {   fUseHLTOnePadCluster = use;}
  66   void SetMinMaxCutAbs(Float_t th)         {   fMinMaxCutAbs=th; }
  67   void SetMinLeftRightCutAbs(Float_t th)   {   fMinLeftRightCutAbs=th;}  // minimal amplitude left right - PRF
  68   void SetMinUpDownCutAbs(Float_t th)      {   fMinUpDownCutAbs=th;}  // minimal amplitude up-down - TRF
  69   void SetMinMaxCutSigma(Float_t th)       {   fMinMaxCutSigma=th; }
  70   void SetMinLeftRightCutSigma(Float_t th) {   fMinLeftRightCutSigma=th;}  // minimal amplitude left right - PRF
  71   void SetMinUpDownCutSigma(Float_t th)    {   fMinUpDownCutSigma=th;}  // minimal amplitude up-down - TRF
  72   void  SetUseTotCharge(Bool_t flag) {fUseTotCharge = flag;}
  73   void  SetCtgRange(Double_t ctgRange) {fCtgRange = ctgRange;}
  74   void  SetUseMultiplicityCorrectionDedx(Bool_t flag) {fUseMultiplicityCorrectionDedx = flag;}
  75
  76   void  SetUseAlignmentTime(Bool_t flag) {fUseAlignmentTime = flag;}
  77   void  SetNeighborRowsDedx(Int_t nRows) {fNeighborRowsDedx = nRows;}
  78   void SetCorrectionHVandPTMode(Int_t value){ fGainCorrectionHVandPTMode =value;}
  79   void SetSkipTimeBins(Double_t value) {fSkipTimeBins=value;}
  80   //
  81   Int_t    GetLastSeedRowSec()       const  { return fLastSeedRowSec;}
  82   Int_t    GetSeedGapPrim()        const  { return fSeedGapPrim;}
  83   Int_t    GetSeedGapSec()         const  { return fSeedGapSec;}
  84   void     SetDoKinks(Bool_t on)   { fBKinkFinder=on; }
  85   Bool_t   GetDoKinks() const      { return fBKinkFinder;}
  86   Double_t GetKinkAngleCutChi2(Int_t index) const {return fKinkAngleCutChi2[index];}
  87   void     SetKinkAngleCutChi2(Int_t index,Double_t value) {fKinkAngleCutChi2[index]=value;}
  88   void     SetSeedGapPrim(Int_t seedGapPrim)         { fSeedGapPrim = seedGapPrim;}
  89   void     SetSeedGapSec(Int_t seedGapSec)          { fSeedGapSec  = seedGapSec;}
  90   Float_t  GetMaxC()    const      { return fMaxC;}
  91   Bool_t   GetSpecialSeeding() const { return fBSpecialSeeding;}
  92   //
  93   //
  94
  95   //
  96   // Correction setup
  97   //
  98   void  SetUseFieldCorrection(Int_t flag){fUseFieldCorrection=flag;}
  99   void  SetUseComposedCorrection(Bool_t flag){fUseComposedCorrection=flag;}
 100   void  SetUseRPHICorrection(Int_t flag){fUseRPHICorrection=flag;}
 101   void  SetUseRadialCorrection(Int_t flag){fUseRadialCorrection=flag;}
 102   void  SetUseQuadrantAlignment(Int_t flag){fUseQuadrantAlignment=flag;}
 103   void  SetUseSectorAlignment(Int_t flag){fUseSectorAlignment=flag;}
 104   void  SetUseDriftCorrectionTime(Int_t flag){fUseDriftCorrectionTime=flag;}
 105   void  SetUseDriftCorrectionGY(Int_t flag){fUseDriftCorrectionGY=flag;}
 106   void  SetUseGainCorrectionTime(Int_t flag){fUseGainCorrectionTime=flag;}
 107   void  SetUseExBCorrection(Int_t flag){fUseExBCorrection=flag;}
 108   void  SetUseTOFCorrection(Bool_t flag) {fUseTOFCorrection = flag;}
 109   void  SetUseIonTailCorrection(Int_t flag) {fUseIonTailCorrection = flag;}
 110   void  SetCrosstalkCorrection(Float_t crosstalkCorrection) {fCrosstalkCorrection= crosstalkCorrection; }
 111   void  SetCrosstalkCorrectionMissingCharge(Float_t crosstalkCorrection) {fCrosstalkCorrectionMissingCharge= crosstalkCorrection; }
 112  //
 113   Int_t GetUseFieldCorrection() const {return fUseFieldCorrection;}
 114   Int_t GetUseComposedCorrection() const {return fUseComposedCorrection;}
 115   Int_t GetUseRPHICorrection() const {return fUseRPHICorrection;}
 116   Int_t GetUseRadialCorrection() const {return fUseRadialCorrection;}
 117   Int_t GetUseQuadrantAlignment() const {return fUseQuadrantAlignment;}
 118   Int_t GetUseSectorAlignment() const {return fUseSectorAlignment;}
 119   Int_t GetUseDriftCorrectionTime() const {return fUseDriftCorrectionTime;}
 120   Int_t GetUseDriftCorrectionGY() const {return fUseDriftCorrectionGY;}
 121   Int_t GetUseGainCorrectionTime() const {return fUseGainCorrectionTime;}
 122   Int_t GetUseExBCorrection() const {return fUseExBCorrection;}
 123   Bool_t GetUseTOFCorrection() {return fUseTOFCorrection;}
 124   Int_t GetUseIonTailCorrection() const {return fUseIonTailCorrection;}
 125   Double_t GetCrosstalkCorrection() const {return fCrosstalkCorrection;}
 126  Double_t GetCrosstalkCorrectionMissingCharge() const {return fCrosstalkCorrectionMissingCharge;}
 127
 128   Bool_t GetUseMultiplicityCorrectionDedx() const {return fUseMultiplicityCorrectionDedx;}
 129   Int_t  GetGainCorrectionHVandPTMode() const  { return   fGainCorrectionHVandPTMode;}
 130   Double_t  GetSkipTimeBins() const {return fSkipTimeBins;}
 131
 132   Bool_t GetUseAlignmentTime() const {return fUseAlignmentTime;}
 133   //
 134   Bool_t   GetUseTotCharge() const {return fUseTotCharge;}          // switch use total or max charge
 135   Float_t  GetMinFraction() const {return fMinFraction;}           // truncated mean - lower threshold
 136   Float_t  GetMaxFraction() const {return fMaxFaction;}            // truncated mean - upper threshold
 137   Int_t    GetNeighborRowsDedx() const {return fNeighborRowsDedx;}
 138
 139   //
 140   void     SetSystematicError(Double_t *systematic){ for (Int_t i=0; i<5;i++) fSystematicErrors[i]=systematic[i];}
 141   const Double_t * GetSystematicError() const { return fSystematicErrors;}
 142   const Double_t * GetSystematicErrorClusterInner() const { return fSystematicErrorClusterInner;}
 143
 144   void    SetUseSystematicCorrelation(Bool_t useCorrelation)  {fUseSystematicCorrelation=useCorrelation;}
 145   Bool_t  GetUseSystematicCorrelation() const { return fUseSystematicCorrelation;}
 146
 147   static   AliTPCRecoParam *GetLowFluxParam();        // make reco parameters for low  flux env.
 148   static   AliTPCRecoParam *GetHighFluxParam();       // make reco parameters for high flux env.
 149   static   AliTPCRecoParam *GetHLTParam(); // special setting for HLT
 150   static   AliTPCRecoParam *GetLaserTestParam(Bool_t bPedestal);  // special setting for laser
 151   static   AliTPCRecoParam *GetCosmicTestParam(Bool_t bPedestal); // special setting for cosmic
 152   //
 153  protected:
 154
 155   Int_t    fUseHLTClusters;  // allows usage of HLT clusters instead of RAW data
 156   Int_t    fUseHLTPreSeeding; // Usage of HLT pre-seeding
 157   Bool_t   fBClusterSharing; // allows or disable cluster sharing during tracking
 158   Double_t fCtgRange;        // +-fCtgRange is the ctg(Theta) window used for clusterization and tracking (MI)
 159   Double_t fMaxSnpTracker;   // max sin of local angle  - for TPC tracker
 160   Double_t fMaxSnpTrack;     // max sin of local angle  - for track
 161   Bool_t   fUseOuterDetectors; // switch - to use the outer detectors
 162   Double_t fMaxChi2TPCTRD;     // maximal allowed chi2 between the TRD in and TPC out to be accepted for refit
 163   Double_t fMaxChi2TPCITS;     // maximal allowed chi2 between the ITS in and TPC out to be accepted for backpropagation
 164   //
 165   //
 166   Double_t fCutSharedClusters[2]; // cut value - maximal amount  of shared clusters
 167   Int_t fClusterMaxRange[2];   // neighborhood  - to define local maxima for cluster
 168   //
 169   //   clusterer parameters
 170   //
 171   Bool_t   fDumpSignal;      // Dump Signal information flag
 172   Int_t    fFirstBin;        // first time bin used by cluster finder
 173   Int_t    fLastBin;         // last time bin  used by cluster finder
 174   Bool_t   fBCalcPedestal;   // calculate Pedestal
 175   Bool_t   fBDoUnfold;       // do unfolding of clusters
 176   Float_t  fDumpAmplitudeMin; // minimal amplitude of signal to be dumped
 177   Float_t  fMaxNoise;        // maximal noise sigma on pad to be used in cluster finder
 178   Int_t    fUseOnePadCluster; // flag - use one pad cluster -0 not use >0 use
 179   Bool_t   fUseHLTOnePadCluster; // flag - use one HLT pad cluster for tracking
 180   Float_t  fMinMaxCutAbs;    // minimal amplitude at cluster maxima
 181   Float_t  fMinLeftRightCutAbs;  // minimal amplitude left right - PRF
 182   Float_t  fMinUpDownCutAbs;  // minimal amplitude up-down - TRF
 183   Float_t  fMinMaxCutSigma;    // minimal amplitude at cluster maxima
 184   Float_t  fMinLeftRightCutSigma;  // minimal amplitude left right - PRF
 185   Float_t  fMinUpDownCutSigma;  // minimal amplitude up-down - TRF
 186   //
 187   //
 188   Float_t  fMaxC;            // maximal curvature for tracking
 189   Bool_t   fBSpecialSeeding; // special seeding with big inclination angles allowed (for Cosmic and laser)
 190   Bool_t   fBKinkFinder;       // do kink finder reconstruction
 191   Double_t fKinkAngleCutChi2[2];   // angular cut for kinks
 192   Int_t    fLastSeedRowSec;     // Most Inner Row to make seeding for secondaries
 193   Int_t    fSeedGapPrim;   // seeding gap for primary tracks
 194   Int_t    fSeedGapSec;   // seeding gap for secondary tracks
 195
 196   //
 197   // Correction switches
 198   //
 199   Int_t fUseFieldCorrection;     // use field correction
 200   Bool_t fUseComposedCorrection; // flag to use composed correction
 201   Int_t fUseRPHICorrection;      // use rphi correction
 202   Int_t fUseRadialCorrection;    // use radial correction
 203   Int_t fUseQuadrantAlignment;   // use quadrant alignment
 204   Int_t fUseSectorAlignment;     // use sector alignment
 205   Int_t fUseDriftCorrectionTime; // use drift correction time
 206   Int_t fUseDriftCorrectionGY;   // use drif correction global y
 207   Int_t fUseGainCorrectionTime;  // use gain correction time
 208   Int_t fUseExBCorrection;       // use ExB correction
 209   Bool_t fUseMultiplicityCorrectionDedx; // use Dedx multiplicity correction
 210   Bool_t fUseAlignmentTime;              // use time dependent alignment correction
 211   Int_t fUseIonTailCorrection;   // use ion tail correction
 212   Double_t fCrosstalkCorrection;   // crosstalk correction factor (fro each signal substracted by (mean signal in wite patch)xfCrosstalkCorrection) - Effect important only after removing oc capacitors in 2012
 213   Double_t fCrosstalkCorrectionMissingCharge;   // crosstalk correction factor - missing charge factor (from each signal substracted by (mean signal in wite patch)xfCrosstalkCorrection) - Effect important only after removing  capacitors in 2012
 214  //
 215   // dEdx switches
 216   //
 217   Bool_t   fUseTotCharge;          // switch use total or max charge
 218   Float_t fMinFraction;           // truncated mean - lower threshold
 219   Float_t fMaxFaction;            // truncated mean - upper threshold
 220   Int_t   fNeighborRowsDedx;      // number of neighboring rows to identify cluster below thres in dEdx calculation 0 -> switch off
 221   Int_t   fGainCorrectionHVandPTMode; // switch for the usage of GainCorrectionHVandPT (see AliTPCcalibDB::GetGainCorrectionHVandPT
 222   Double_t fSkipTimeBins;        // number of time bins to be skiiped (corrupted signal druing gating opening)
 223
 224   Bool_t fUseTOFCorrection;  // switch - kTRUE use TOF correction kFALSE - do not use
 225   //
 226   //  misscalibration
 227   //
 228   Double_t fSystematicErrors[5];  //systematic errors in the track parameters - to be added to TPC covariance matrix
 229   Double_t fSystematicErrorClusterInner[2];  // systematic error of the cluster - used to downscale the information
 230   Bool_t fUseSystematicCorrelation;         // switch to use the correlation for the sys
 231 public:
 232   static Bool_t fgUseTimeCalibration; // flag usage the time dependent calibration
 233                                       // to be switched off for pass 0 reconstruction
 234                                       // Use static function, other option will be to use
 235                                       // additional specific storage ?
 236   ClassDef(AliTPCRecoParam, 21)
 237 };
 238
 239
 240 #endif