New task to check TPC-ITS track prolongation eff with cosmics

[u/mrichter/AliRoot.git] / TPC / AliTPCRecoParam.cxx

diff --git a/TPC/AliTPCRecoParam.cxx b/TPC/AliTPCRecoParam.cxx
index 7bedd22c03b95a172c854a02b2b3e8c348989bab..6a077623ba1ab8abdec16e8eb1cdaa5a51842ade 100644 (file)
--- a/TPC/AliTPCRecoParam.cxx
+++ b/TPC/AliTPCRecoParam.cxx
@@ -18,7 +18,33 @@
 //                                                                           //
 // Class with TPC reconstruction parameters                                  //
 //                                                                           //  
-//                                                                           //
+//
+/*
+  The reconstruction parameters are used in the AliTPCclustererMI and AliTPCtrackerMI
+  
+  They are retrieved:
+  0. User speciefied it in reconstruction macro
+  1. if (not 0) from OCDB  - AliTPCcalibDB::GetRecoParam(eventtype)
+  2. if (not 0 or 1) default parameter - High flux enevironment used  
+
+  FIXME:
+  In the future  reconstruction parameters should be changed on event basis
+  But for the moment, event types are still not defined 
+
+
+  // Setting for systematic errors addition
+  [0] - systematic RMSY
+  [1] - systematic RMSZ
+  [2] - systematic RMSSNP
+  [3] - systematic RMSTheta
+  [4] - systematic RMSCuravture -  systematic error in 1/cm not in 1/pt
+  //
+  //  How to add it example - 3 mm systematic error y, 3 cm systematic error z (drift)
+  Double_t sysError[5]={0.3,3, 0.3/150., 3./150.,0.3/(150*150.)}
+  param->SetSystematicError(sysError);
+
+*/
+                                                                           //
 ///////////////////////////////////////////////////////////////////////////////
 
 
@@ -31,20 +57,51 @@ ClassImp(AliTPCRecoParam)
 
 //_____________________________________________________________________________
 AliTPCRecoParam::AliTPCRecoParam():
+  AliDetectorRecoParam(),
+  fBClusterSharing(kTRUE),
   fCtgRange(1.05),       
   fMaxSnpTracker(0.95),
   fMaxSnpTrack(0.999),
+  fDumpSignal(kFALSE),
   fFirstBin(0),
   fLastBin(-1),
   fBCalcPedestal(kFALSE),
   fBDoUnfold(kTRUE),
   fDumpAmplitudeMin(100),
-  fMaxNoise(3.),
-  fBKinkFinder(kTRUE)
+  fMaxNoise(2.),
+  //
+  fMinMaxCutAbs(4.),
+  fMinLeftRightCutAbs(6.),
+  fMinUpDownCutAbs(6.),
+  //
+  fMinMaxCutSigma(4.),
+  fMinLeftRightCutSigma(7.),
+  fMinUpDownCutSigma(8.),
+  fMaxC(0.3),
+  fBSpecialSeeding(kFALSE),
+  fBKinkFinder(kTRUE),
+  fLastSeedRowSec(120),
+  fUseFieldCorrection(0),      // use field correction
+  fUseRPHICorrection(0),      // use rphi correction
+  fUseRadialCorrection(0),    // use radial correction
+  fUseQuadrantAlignment(0),   // use quadrant alignment
+  fUseSectorAlignment(0),     // use sector alignment
+  fUseDriftCorrectionTime(0), // use drift correction time
+  fUseDriftCorrectionGY(0),   // use drif correction global y
+  fUseGainCorrectionTime(0),  // use gain correction time
+  fUseExBCorrection(1),  // use ExB correction
+  //
+  fUseTotCharge(kTRUE),          // switch use total or max charge
+  fMinFraction(0.01),           // truncated mean - lower threshold
+  fMaxFaction(0.7),            // truncated mean - upper threshold
+  fUseTOFCorrection(kTRUE)
 {
   //
   // constructor
   //
+  SetName("TPC");
+  SetTitle("TPC");
+  for (Int_t i=0;i<5;i++) fSystematicErrors[i]=0;
 }
 
 //_____________________________________________________________________________
@@ -66,6 +123,8 @@ AliTPCRecoParam *AliTPCRecoParam::GetLowFluxParam(){
   param->fCtgRange = 10;
   param->fFirstBin = 0;
   param->fLastBin  = 1000;
+  param->SetName("Low Flux");
+  param->SetTitle("Low Flux");
   return param;
 }
 
@@ -76,7 +135,27 @@ AliTPCRecoParam *AliTPCRecoParam::GetHighFluxParam(){
   AliTPCRecoParam *param = new AliTPCRecoParam;
   param->fCtgRange = 1.05;
   param->fFirstBin = 0;
-  param->fLastBin  = 1000;
+  param->fLastBin  = 1000;  
+  param->fUseTotCharge=kFALSE;
+  param->SetName("High Flux");
+  param->SetTitle("High Flux");
+  return param;
+}
+
+AliTPCRecoParam *AliTPCRecoParam::GetHLTParam(){
+  //
+  // make reco parameters for high flux env.
+  //
+  AliTPCRecoParam *param = new AliTPCRecoParam;
+  param->fCtgRange = 1.05;
+  param->fFirstBin = 80;
+  param->fLastBin  = 1000;  
+  param->fMaxSnpTracker = 0.9; 
+  param->fMaxC          = 0.06; 
+  //
+  param->SetName("Hlt Param");
+  param->SetTitle("Hlt Param"); 
+  param->fBKinkFinder   = kFALSE;
   return param;
 }
 
@@ -85,12 +164,22 @@ AliTPCRecoParam *AliTPCRecoParam::GetLaserTestParam(Bool_t bPedestal){
   // special setting for laser
   //
   AliTPCRecoParam *param = new AliTPCRecoParam;
+  param->fDumpSignal=kTRUE;
   param->fCtgRange = 10.05;
   param->fFirstBin = 0;
   param->fLastBin  = 1000;
   param->fBCalcPedestal = bPedestal;
   param->fBDoUnfold     = kFALSE;
+  param->fDumpAmplitudeMin = 150;
   param->fBKinkFinder   = kFALSE;
+  param->fMaxSnpTracker = 0.98;
+  param->fMaxC          = 0.02;
+  param->fBSpecialSeeding = kTRUE;
+  param->fUseTOFCorrection=kFALSE;
+  //
+  //
+  param->SetName("Laser Flux");
+  param->SetTitle("Laser Flux");
   return param;
 }
 
@@ -99,12 +188,19 @@ AliTPCRecoParam *AliTPCRecoParam::GetCosmicTestParam(Bool_t bPedestal){
   // special setting for cosmic 
   // 
   AliTPCRecoParam *param = new AliTPCRecoParam;
+  param->fDumpSignal=kTRUE;
   param->fCtgRange = 10.05;    // full TPC
   param->fFirstBin = 60;
   param->fLastBin  = 1000;
   param->fBCalcPedestal = bPedestal;
   param->fBDoUnfold     = kFALSE;
+  param->fBSpecialSeeding = kTRUE;
+  param->fMaxC          = 0.07;
   param->fBKinkFinder   = kFALSE;
+  param->fUseTOFCorrection =kFALSE;
+  param->SetName("Cosmic Flux");
+  param->SetTitle("Cosmic Flux");
+
   return param;
 }