Adding setter and getter for systematic uncertainty of the track parameters
authormarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 1 Jun 2008 13:41:50 +0000 (13:41 +0000)
committermarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 1 Jun 2008 13:41:50 +0000 (13:41 +0000)
  // 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.,1/(0.3*150*150.)}
  param->SetSystematicError(sysError);

(Marian)

TPC/AliTPCRecoParam.cxx
TPC/AliTPCRecoParam.h

index f2fa94b..6b1ac7d 100644 (file)
   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.,1/(0.3*150*150.)}
+  param->SetSystematicError(sysError);
+
 */
                                                                            //
 ///////////////////////////////////////////////////////////////////////////////
@@ -74,6 +86,7 @@ AliTPCRecoParam::AliTPCRecoParam():
   //
   SetName("TPC");
   SetTitle("TPC");
+  for (Int_t i=0;i<5;i++) fSystematicErrors[i]=0;
 }
 
 //_____________________________________________________________________________
@@ -122,7 +135,8 @@ AliTPCRecoParam *AliTPCRecoParam::GetHLTParam(){
   param->fFirstBin = 80;
   param->fLastBin  = 1000;  
   param->fMaxSnpTracker = 0.9; 
-  param->fMaxC          = 0.06;
+  param->fMaxC          = 0.06; 
+  param->fBYMirror      = kFALSE;
   //
   param->SetName("Hlt Param");
   param->SetTitle("Hlt Param"); 
index 2d3b945..327c6a9 100644 (file)
@@ -55,6 +55,10 @@ class AliTPCRecoParam : public AliDetectorRecoParam
   Bool_t   GetSpecialSeeding() const { return fBSpecialSeeding;}
   Bool_t   GetBYMirror() const { return fBYMirror;}
   void     SetBYMirror(Bool_t mirror)  { fBYMirror = mirror;} //
+  //
+  void     SetSystematicError(Double_t *systematic){ for (Int_t i=0; i<5;i++) fSystematicErrors[i]=systematic[i];}
+  const Double_t * GetSystematicError() const { return fSystematicErrors;}
+
   static   AliTPCRecoParam *GetLowFluxParam();        // make reco parameters for low  flux env.
   static   AliTPCRecoParam *GetHighFluxParam();       // make reco parameters for high flux env. 
   static   AliTPCRecoParam *GetHLTParam(); // special setting for HLT  
@@ -89,7 +93,11 @@ class AliTPCRecoParam : public AliDetectorRecoParam
   Bool_t   fBSpecialSeeding; // special seeding with big inclination angles allowed (for Cosmic and laser)
   Bool_t   fBKinkFinder;     // do kink finder reconstruction
   Int_t    fLastSeedRowSec;     // Most Inner Row to make seeding for secondaries
-  ClassDef(AliTPCRecoParam, 2)
+  //
+  //  misscalibration
+  //
+  Double_t fSystematicErrors[5];  //systematic errors in the track parameters - to be added to TPC covariance matrix    
+  ClassDef(AliTPCRecoParam, 3)
 };