TPC/AliTPCRecoParam.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16
  17 ///////////////////////////////////////////////////////////////////////////////
  18 //                                                                           //
  19 // Class with TPC reconstruction parameters                                  //
  20 //                                                                           //
  21 //
  22 /*
  23   The reconstruction parameters are used in the AliTPCclustererMI and AliTPCtrackerMI
  24
  25   They are retrieved:
  26   0. User speciefied it in reconstruction macro
  27   1. if (not 0) from OCDB  - AliTPCcalibDB::GetRecoParam(eventtype)
  28   2. if (not 0 or 1) default parameter - High flux enevironment used
  29
  30   FIXME:
  31   In the future  reconstruction parameters should be changed on event basis
  32   But for the moment, event types are still not defined
  33
  34
  35   // Setting for systematic errors addition
  36   [0] - systematic RMSY
  37   [1] - systematic RMSZ
  38   [2] - systematic RMSSNP
  39   [3] - systematic RMSTheta
  40   [4] - systematic RMSCuravture -  systematic error in 1/cm not in 1/pt
  41   //
  42   //  How to add it example - 3 mm systematic error y, 3 cm systematic error z (drift)
  43   Double_t sysError[5]={0.3,3, 0.3/150., 3./150.,0.3/(150*150.)}
  44   param->SetSystematicError(sysError);
  45
  46 */
  47                                                                            //
  48 ///////////////////////////////////////////////////////////////////////////////
  49
  50
  51 #include "AliTPCRecoParam.h"
  52
  53 ClassImp(AliTPCRecoParam)
  54
  55
  56 Bool_t AliTPCRecoParam::fgUseTimeCalibration=kTRUE; // flag usage the time dependent calibration
  57                                       // to be switched off for pass 0 reconstruction
  58                                       // Use static function, other option will be to use
  59                                       // additional specific storage ?
  60
  61 //_____________________________________________________________________________
  62 AliTPCRecoParam::AliTPCRecoParam():
  63   AliDetectorRecoParam(),
  64   fBClusterSharing(kTRUE),
  65   fCtgRange(1.05),
  66   fMaxSnpTracker(0.95),
  67   fMaxSnpTrack(0.999),
  68   fUseOuterDetectors(kFALSE),
  69   fDumpSignal(kFALSE),
  70   fFirstBin(0),
  71   fLastBin(-1),
  72   fBCalcPedestal(kFALSE),
  73   fBDoUnfold(kTRUE),
  74   fDumpAmplitudeMin(100),
  75   fMaxNoise(2.),
  76   //
  77   fUseOnePadCluster(kTRUE),
  78   fMinMaxCutAbs(4.),
  79   fMinLeftRightCutAbs(6.),
  80   fMinUpDownCutAbs(6.),
  81   //
  82   fMinMaxCutSigma(4.),
  83   fMinLeftRightCutSigma(7.),
  84   fMinUpDownCutSigma(8.),
  85   fMaxC(0.3),
  86   fBSpecialSeeding(kFALSE),
  87   fBKinkFinder(kTRUE),
  88   fLastSeedRowSec(120),
  89   fSeedGapPrim(6),
  90   fSeedGapSec(6),
  91   fUseFieldCorrection(2),      // use field correction
  92   fUseComposedCorrection(kFALSE),      // use field correction
  93   fUseRPHICorrection(0),      // use rphi correction
  94   fUseRadialCorrection(0),    // use radial correction
  95   fUseQuadrantAlignment(0),   // use quadrant alignment
  96   fUseSectorAlignment(0),     // use sector alignment
  97   fUseDriftCorrectionTime(1), // use drift correction time
  98   fUseDriftCorrectionGY(1),   // use drif correction global y
  99   fUseGainCorrectionTime(0),  // use gain correction time
 100   fUseExBCorrection(1),  // use ExB correction
 101   //
 102   fUseTotCharge(kTRUE),          // switch use total or max charge
 103   fMinFraction(0.01),           // truncated mean - lower threshold
 104   fMaxFaction(0.7),            // truncated mean - upper threshold
 105   fUseTOFCorrection(kTRUE)
 106 {
 107   //
 108   // constructor
 109   //
 110   SetName("TPC");
 111   SetTitle("TPC");
 112   for (Int_t i=0;i<5;i++) fSystematicErrors[i]=0;
 113   fCutSharedClusters[0]=0.5; // maximal allowed fraction of shared clusters - shorter track
 114   fCutSharedClusters[1]=0.25; // maximal allowed fraction of shared clusters - longer  track
 115   fClusterMaxRange[0]=1;     // y - pad      range
 116   fClusterMaxRange[1]=1;     // z - time bin range
 117   fKinkAngleCutChi2[0]=9;    // angular cut for kink finder - to create a kink
 118                              // ~ about 5 % rate  for high pt kink finder
 119   fKinkAngleCutChi2[1]=12;    // angular cut for kink finder - to use the partial track                             // form kink
 120                              // ~ about 2 % rate  for high pt kink finder
 121 }
 122
 123 //_____________________________________________________________________________
 124 AliTPCRecoParam::~AliTPCRecoParam()
 125 {
 126   //
 127   // destructor
 128   //
 129 }
 130
 131
 132
 133
 134 AliTPCRecoParam *AliTPCRecoParam::GetLowFluxParam(){
 135   //
 136   // make default reconstruction  parameters for low  flux env.
 137   //
 138   AliTPCRecoParam *param = new AliTPCRecoParam;
 139   param->fCtgRange = 10;
 140   param->fFirstBin = 0;
 141   param->fLastBin  = 1000;
 142   param->SetName("Low Flux");
 143   param->SetTitle("Low Flux");
 144   return param;
 145 }
 146
 147 AliTPCRecoParam *AliTPCRecoParam::GetHighFluxParam(){
 148   //
 149   // make reco parameters for high flux env.
 150   //
 151   AliTPCRecoParam *param = new AliTPCRecoParam;
 152   param->fCtgRange = 1.05;
 153   param->fFirstBin = 0;
 154   param->fLastBin  = 1000;
 155   param->fUseTotCharge=kFALSE;
 156   param->SetName("High Flux");
 157   param->SetTitle("High Flux");
 158   return param;
 159 }
 160
 161 AliTPCRecoParam *AliTPCRecoParam::GetHLTParam(){
 162   //
 163   // make reco parameters for high flux env.
 164   //
 165   AliTPCRecoParam *param = new AliTPCRecoParam;
 166   param->fCtgRange = 1.05;
 167   param->fFirstBin = 80;
 168   param->fLastBin  = 1000;
 169   param->fMaxSnpTracker = 0.9;
 170   param->fMaxC          = 0.06;
 171   //
 172   param->SetName("Hlt Param");
 173   param->SetTitle("Hlt Param");
 174   param->fBKinkFinder   = kFALSE;
 175   return param;
 176 }
 177
 178 AliTPCRecoParam *AliTPCRecoParam::GetLaserTestParam(Bool_t bPedestal){
 179   //
 180   // special setting for laser
 181   //
 182   AliTPCRecoParam *param = new AliTPCRecoParam;
 183   param->fDumpSignal=kTRUE;
 184   param->fCtgRange = 10.05;
 185   param->fFirstBin = 0;
 186   param->fLastBin  = 1000;
 187   param->fBCalcPedestal = bPedestal;
 188   param->fBDoUnfold     = kFALSE;
 189   param->fDumpAmplitudeMin = 150;
 190   param->fBKinkFinder   = kFALSE;
 191   param->fMaxSnpTracker = 0.98;
 192   param->fMaxC          = 0.02;
 193   param->fBSpecialSeeding = kTRUE;
 194   param->fUseTOFCorrection=kFALSE;
 195   //
 196   //
 197   param->SetName("Laser Flux");
 198   param->SetTitle("Laser Flux");
 199   return param;
 200 }
 201
 202 AliTPCRecoParam *AliTPCRecoParam::GetCosmicTestParam(Bool_t bPedestal){
 203   //
 204   // special setting for cosmic
 205   //
 206   AliTPCRecoParam *param = new AliTPCRecoParam;
 207   param->fDumpSignal=kTRUE;
 208   param->fCtgRange = 10.05;    // full TPC
 209   param->fFirstBin = 60;
 210   param->fLastBin  = 1000;
 211   param->fBCalcPedestal = bPedestal;
 212   param->fBDoUnfold     = kFALSE;
 213   param->fBSpecialSeeding = kTRUE;
 214   param->fMaxC          = 0.07;
 215   param->fBKinkFinder   = kFALSE;
 216   param->fUseTOFCorrection =kFALSE;
 217   param->SetName("Cosmic Flux");
 218   param->SetTitle("Cosmic Flux");
 219
 220   return param;
 221 }
 222
 223
 224 Bool_t  AliTPCRecoParam::GetUseTimeCalibration(){
 225   //
 226   // get
 227   //
 228   return fgUseTimeCalibration;
 229 }
 230 void    AliTPCRecoParam::SetUseTimeCalibration(Bool_t useTimeCalibration) {
 231   //
 232   // set
 233   //
 234   fgUseTimeCalibration = useTimeCalibration;
 235 }
 236