]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/AliITSRecoParam.h
cleanup
[u/mrichter/AliRoot.git] / ITS / AliITSRecoParam.h
index 29fc5b4027af5344a94e795fa5a8967bf12c2c76..21e06409cad7260670ef8a1a78558f3b5f61f1f6 100644 (file)
@@ -242,6 +242,9 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Float_t GetTrackleterZetaWindowL1() const {return fTrackleterZetaWindowL1;}
   void    SetTrackleterZetaWindowL2(Float_t w=0.40) {fTrackleterZetaWindowL2=w; return;}
   Float_t GetTrackleterZetaWindowL2() const {return fTrackleterZetaWindowL2;}
+  void    SetTrackleterBuildCl2TrkRefs(Bool_t v=kTRUE) {fTrackleterBuildCl2TrkRefs = v;}
+  Bool_t  GetTrackleterBuildCl2TrkRefs() const { return fTrackleterBuildCl2TrkRefs;}
+  //
   void    SetUpdateOncePerEventPlaneEff(Bool_t use=kTRUE) {fUpdateOncePerEventPlaneEff=use; return;}
   Bool_t  GetUpdateOncePerEventPlaneEff() const {return fUpdateOncePerEventPlaneEff;}
   void    SetMinContVtxPlaneEff(Int_t n=3) {fMinContVtxPlaneEff=n; return;}
@@ -266,6 +269,11 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Double_t GetNSigXFromBoundaryPlaneEff() const {return fNSigXFromBoundaryPlaneEff;}
   void SetNSigZFromBoundaryPlaneEff(Double_t nsigz=1.) {fNSigZFromBoundaryPlaneEff=nsigz;}
   Double_t GetNSigZFromBoundaryPlaneEff() const {return fNSigZFromBoundaryPlaneEff;}
+  void SetDCACutPlaneEff(Double_t dcacpe=999.) {fDCACutPlaneEff=dcacpe;}
+  Double_t GetDCACutPlaneEff() const {return fDCACutPlaneEff;}
+  void SetVertexChi2CutPlaneEff(Double_t vtxchipe=999999999.) {fVertexChi2CutPlaneEff=vtxchipe;}
+  Double_t GetVertexChi2CutPlaneEff() const {return fVertexChi2CutPlaneEff;}
+
   //
   void   SetImproveWithVertex(Bool_t impr=kFALSE) { fImproveWithVertex=impr; return; }
   Bool_t GetImproveWithVertex() const { return fImproveWithVertex; }
@@ -372,6 +380,8 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Float_t GetTrackleterZetaOverlapCut() const {return fTrackleterZetaOverlapCut;}
   void    SetTrackleterPhiRotationAngle(Float_t w=0.0) {fTrackleterPhiRotationAngle=w;}
   Float_t GetTrackleterPhiRotationAngle() const {return fTrackleterPhiRotationAngle;}
+  Bool_t  GetTrackleterStoreSPD2SingleCl()      const {return fTrackleterStoreSPD2SingleCl;}
+  void    SetTrackleterStoreSPD2SingleCl(Bool_t v=kTRUE) {fTrackleterStoreSPD2SingleCl = v;}
   //
   void    SetTrackleterNStdDevCut(Float_t f=1.)          {fTrackleterNStdDev = f<0.01 ? 0.01 : f;}
   Float_t GetTrackleterNStdDevCut()               const  {return fTrackleterNStdDev;}
@@ -469,6 +479,11 @@ class AliITSRecoParam : public AliDetectorRecoParam
   void SetTanLorentzAngleHolesSSD(Float_t la) {fTanLorentzAngleHolesSSD=la;}
   void SetTanLorentzAngleElectronsSSD(Float_t la) {fTanLorentzAngleElectronsSSD=la;}
 
+  // Option for local reconstruction
+  Bool_t SetOptReco(TString r);
+  void ReconstructOnlySPD(){fOptReco="SPD";}
+  TString GetOptReco() const {return fOptReco;}
+
   //
   enum {kMaxClusterPerLayer=70000}; //7000*10;   // max clusters per layer
   enum {kMaxClusterPerLayer5=28000};//7000*10*2/5;  // max clusters per layer
@@ -480,34 +495,34 @@ class AliITSRecoParam : public AliDetectorRecoParam
   static const Int_t fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]; // array with layers not to skip
   static const Int_t fgkLastLayerToTrackTo=0;     // innermost layer
   static const Int_t fgkMaxDetectorPerLayer=1000; // max clusters per layer
-  static const Double_t fgkriw=80.0;              // TPC inner wall radius
-  static const Double_t fgkdiw=0.0053;            // TPC inner wall x/X0
-  static const Double_t fgkX0iw=30.0;             // TPC inner wall X0 
-  static const Double_t fgkrcd=61.0;              // TPC central drum radius
-  static const Double_t fgkdcd=0.0053;            // TPC central drum x/X0
-  static const Double_t fgkX0cd=30.0;             // TPC central drum X0
-  static const Double_t fgkyr=12.8;               // TPC rods y (tracking c.s.)
-  static const Double_t fgkdr=0.03;               // TPC rods x/X0
-  static const Double_t fgkzm=0.2;                // TPC membrane z
-  static const Double_t fgkdm=0.40;               // TPC membrane x/X0
-  static const Double_t fgkrs=50.0;               // ITS screen radius
-  static const Double_t fgkds=0.001;              // ITS screed x/X0
-  static const Double_t fgkrInsideITSscreen=49.0; // inside ITS screen radius
-  static const Double_t fgkrInsideSPD1=3.5;       // inside SPD1 radius
-  static const Double_t fgkrPipe=3.;              // pipe radius
-  static const Double_t fgkrInsidePipe=2.7;       // inside pipe radius
-  static const Double_t fgkrOutsidePipe=3.3;      // outside pipe radius
-  static const Double_t fgkdPipe=0.0028;          // pipe x/X0
-  static const Double_t fgkrInsideShield[2]; // inside SPD (0) SDD (1) shield radius
-  static const Double_t fgkrOutsideShield[2]; // outside SPD (0) SDD (1) shield radius
-  static const Double_t fgkdshield[2];        // SPD (0) SDD (1) shield x/X0
-  static const Double_t fgkX0shield[2];       // SPD (0) SDD (1) shield X0
-  static const Double_t fgkX0Air=21.82;       // air X0
-  static const Double_t fgkX0Be=65.19;        // Berillium X0
-  static const Double_t fgkBoundaryWidth=0.2; // to define track at detector boundary
-  static const Double_t fgkDeltaXNeighbDets=0.5; // max difference in radius between neighbouring detectors 
-  static const Double_t fgkSPDdetzlength=6.960;     // SPD ladder length in z (=7.072-2*0.056)
-  static const Double_t fgkSPDdetxlength=1.298;     // SPD ladder length in x (=1.410-2*0.056)
+  static const Double_t fgkriw;                   // TPC inner wall radius
+  static const Double_t fgkdiw;                   // TPC inner wall x/X0
+  static const Double_t fgkX0iw;                  // TPC inner wall X0 
+  static const Double_t fgkrcd;                   // TPC central drum radius
+  static const Double_t fgkdcd;                   // TPC central drum x/X0
+  static const Double_t fgkX0cd;                  // TPC central drum X0
+  static const Double_t fgkyr;                    // TPC rods y (tracking c.s.)
+  static const Double_t fgkdr;                    // TPC rods x/X0
+  static const Double_t fgkzm;                    // TPC membrane z
+  static const Double_t fgkdm;                    // TPC membrane x/X0
+  static const Double_t fgkrs;                    // ITS screen radius
+  static const Double_t fgkds;                    // ITS screed x/X0
+  static const Double_t fgkrInsideITSscreen;      // inside ITS screen radius
+  static const Double_t fgkrInsideSPD1;           // inside SPD1 radius
+  static const Double_t fgkrPipe;                 // pipe radius
+  static const Double_t fgkrInsidePipe;           // inside pipe radius
+  static const Double_t fgkrOutsidePipe;          // outside pipe radius
+  static const Double_t fgkdPipe;                 // pipe x/X0
+  static const Double_t fgkrInsideShield[2];      // inside SPD (0) SDD (1) shield radius
+  static const Double_t fgkrOutsideShield[2];     // outside SPD (0) SDD (1) shield radius
+  static const Double_t fgkdshield[2];            // SPD (0) SDD (1) shield x/X0
+  static const Double_t fgkX0shield[2];           // SPD (0) SDD (1) shield X0
+  static const Double_t fgkX0Air;                 // air X0
+  static const Double_t fgkX0Be;                  // Berillium X0
+  static const Double_t fgkBoundaryWidth;         // to define track at detector boundary
+  static const Double_t fgkDeltaXNeighbDets;      // max difference in radius between neighbouring detectors 
+  static const Double_t fgkSPDdetzlength;         // SPD ladder length in z (=7.072-2*0.056)
+  static const Double_t fgkSPDdetxlength;         // SPD ladder length in x (=1.410-2*0.056)
 
 
   Int_t  fTracker;  // ITS tracker to be used (see AliITSReconstructor)
@@ -627,6 +642,8 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Bool_t fOnlyConstraintPlaneEff;  // if kTRUE, use only constrained tracks at primary vertex for Plane Eff.
   Double_t fNSigXFromBoundaryPlaneEff;  // accept one track for PlaneEff if distance from border (in loc x or z)
   Double_t fNSigZFromBoundaryPlaneEff;  // is greater than fNSigXFromBoundaryPlaneEff * Track_precision
+  Double_t fDCACutPlaneEff;  // this set the cut on DCA in rphi plane when evaluating PlaneEff(SPD0)
+  Double_t fVertexChi2CutPlaneEff; // and also with a cut on the chi2
 
   Bool_t fImproveWithVertex;    // use the method AliITStrackV2::Improve() to point to the vertex during prolongation
   Bool_t fExtendedEtaAcceptance;  // enable jumping from TPC to SPD at large eta (MI)
@@ -680,7 +697,10 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Float_t fTrackleterPhiRotationAngle;             // Angle to rotate cluster in the SPD inner layer for combinatorial reco only
   Float_t fTrackleterNStdDev;      // cut on the number of standard deviations
   Bool_t  fScaleDTBySin2T;         // scale Dtheta by 1/sin^2(theta)
-
+  Bool_t  fTrackleterStoreSPD2SingleCl;            // request storing of L2 singles
+  //
+  Bool_t  fTrackleterBuildCl2TrkRefs;           // build cluster to track references in AliMultiplicity
+  //
   Bool_t fUseCosmicRunShiftsSSD; // SSD time shifts for cosmic run 2007/2008 (use for data taken up to 18 sept 2008)
 
 
@@ -694,7 +714,7 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Float_t fVertexerFastSmearY;  // gaussian sigma for y MC vertex smearing
   Float_t fVertexerFastSmearZ;  // gaussian sigma for z MC vertex smearing
 
-  // PWG1/AliAlignmentDataFilterITS configuration
+  // PWGPP/AliAlignmentDataFilterITS configuration
   Bool_t  fAlignFilterCosmics;            // flag for cosmics case
   Bool_t  fAlignFilterCosmicMergeTracks;  // merge cosmic tracks
   Int_t   fAlignFilterMinITSPoints;       // min points per track
@@ -734,13 +754,16 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Float_t fMultCutGammaSFromDecay;        // min path*P for gammas
   Float_t fMultCutK0SFromDecay;           // min path*P for K0s
   Float_t fMultCutMaxDCA;                 // max DCA for V0 at ESD vertex
-  //
   // Lorentz angle
   Bool_t fCorrectLorentzAngleSPD;         // flag to enable correction
   Float_t fTanLorentzAngleHolesSPD;       // angle for holes in SPD
   Bool_t fCorrectLorentzAngleSSD;         // flag to enable correction
   Float_t fTanLorentzAngleHolesSSD;       // tan(angle) for holes in SSD @ B = 0.5 T
   Float_t fTanLorentzAngleElectronsSSD;   // tan(angle) for electrons in SSD @ B = 0.5 T
+  //
+  // Possibility of reconstructing only part of the ITS
+  TString fOptReco;                      // "All" by default. It can be any 
+                                         // combination of "SPD" "SDD" and "SSD"
 
  private:
   AliESDV0Params * fESDV0Params;  // declare the AliESDV0Params to be able to used in AliITSV0Finder
@@ -748,7 +771,7 @@ class AliITSRecoParam : public AliDetectorRecoParam
   AliITSRecoParam(const AliITSRecoParam & param);
   AliITSRecoParam & operator=(const AliITSRecoParam &param);
 
-  ClassDef(AliITSRecoParam,39) // ITS reco parameters
+  ClassDef(AliITSRecoParam,52) // ITS reco parameters
 };
 
 #endif