]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/AliITSRecoParam.h
New algorithm for high multiplicity events (trace algorithm besides the existing...
[u/mrichter/AliRoot.git] / ITS / AliITSRecoParam.h
index 334092f175fad709318abd87926408dd6684a06f..29fc5b4027af5344a94e795fa5a8967bf12c2c76 100644 (file)
@@ -15,7 +15,9 @@
 
 #include "AliDetectorRecoParam.h"
 #include "AliITSgeomTGeo.h"
-#include "AliESDV0Params.h"
+//#include "AliESDV0Params.h"
+
+class AliESDV0Params;
 
 class AliITSRecoParam : public AliDetectorRecoParam
 {
@@ -30,10 +32,10 @@ class AliITSRecoParam : public AliDetectorRecoParam
 
   static Int_t GetLayersNotToSkip(Int_t i) { return fgkLayersNotToSkip[i]; }
   static Int_t GetLastLayerToTrackTo() { return fgkLastLayerToTrackTo; }
-  static Int_t GetMaxClusterPerLayer() { return fgkMaxClusterPerLayer; }
-  static Int_t GetMaxClusterPerLayer5() { return fgkMaxClusterPerLayer5; }
-  static Int_t GetMaxClusterPerLayer10() { return fgkMaxClusterPerLayer10; }
-  static Int_t GetMaxClusterPerLayer20() { return fgkMaxClusterPerLayer20; }
+  static Int_t GetMaxClusterPerLayer() { return kMaxClusterPerLayer; }
+  static Int_t GetMaxClusterPerLayer5() { return kMaxClusterPerLayer5; }
+  static Int_t GetMaxClusterPerLayer10() { return kMaxClusterPerLayer10; }
+  static Int_t GetMaxClusterPerLayer20() { return kMaxClusterPerLayer20; }
   static Int_t GetMaxDetectorPerLayer() { return fgkMaxDetectorPerLayer; }
   static Double_t Getriw() { return fgkriw; }
   static Double_t Getdiw() { return fgkdiw; }
@@ -123,6 +125,8 @@ class AliITSRecoParam : public AliDetectorRecoParam
   void SetSPDVertexerPileupAlgoZ(){fVtxr3DPileupAlgo=0;}
   void SetSPDVertexerPileupAlgo3DTwoSteps(){fVtxr3DPileupAlgo=1;}
   void SetSPDVertexerPileupAlgo3DOneShot(){fVtxr3DPileupAlgo=2;}
+  void SetSPDVertexerHighMultAlgoDownscale(){fVtxr3DHighMultAlgo=0;}
+  void SetSPDVertexerHighMultAlgoTraces(){fVtxr3DHighMultAlgo=1;}
   //
   Bool_t   GetSelectBestMIP03()                 const {return fSelectBestMIP03;}
   Bool_t   GetFlagFakes()                       const {return fFlagFakes;}
@@ -139,6 +143,7 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Float_t  GetVertexer3DTightDeltaPhiCut() const {return fVtxr3DPhiCutTight;}
   Float_t  GetVertexer3DDCACut() const {return fVtxr3DDCACut;}
   Int_t    GetSPDVertexerPileupAlgo() const {return fVtxr3DPileupAlgo;}
+  UChar_t  GetSPDVertexerHighMultAlgo() const {return fVtxr3DHighMultAlgo;}
 
   Double_t GetSigmaY2(Int_t i) const { return fSigmaY2[i]; }
   Double_t GetSigmaZ2(Int_t i) const { return fSigmaZ2[i]; }
@@ -315,6 +320,9 @@ class AliITSRecoParam : public AliDetectorRecoParam
   void SetMaxSPDcontrForSAToUseAllClusters(Int_t contr=50) { fMaxSPDcontrForSAToUseAllClusters=contr; return; }
   Int_t GetMaxSPDcontrForSAToUseAllClusters() const { return fMaxSPDcontrForSAToUseAllClusters; }
 
+  void   SetSAUsedEdxInfo(Bool_t opt=kTRUE) { fSAUsedEdxInfo=opt; return; }
+  Bool_t GetSAUsedEdxInfo() const { return fSAUsedEdxInfo; }
+
   void   SetFindV0s(Bool_t find=kTRUE) { fFindV0s=find; return; }
   Bool_t GetFindV0s() const { return fFindV0s; }
 
@@ -462,44 +470,44 @@ class AliITSRecoParam : public AliDetectorRecoParam
   void SetTanLorentzAngleElectronsSSD(Float_t la) {fTanLorentzAngleElectronsSSD=la;}
 
   //
-  enum {fgkMaxClusterPerLayer=70000}; //7000*10;   // max clusters per layer
-  enum {fgkMaxClusterPerLayer5=28000};//7000*10*2/5;  // max clusters per layer
-  enum {fgkMaxClusterPerLayer10=14000};//7000*10*2/10; // max clusters per layer
-  enum {fgkMaxClusterPerLayer20=7000};//7000*10*2/20; // max clusters per layer
+  enum {kMaxClusterPerLayer=70000}; //7000*10;   // max clusters per layer
+  enum {kMaxClusterPerLayer5=28000};//7000*10*2/5;  // max clusters per layer
+  enum {kMaxClusterPerLayer10=14000};//7000*10*2/10; // max clusters per layer
+  enum {kMaxClusterPerLayer20=7000};//7000*10*2/20; // max clusters per layer
 
  protected:
   //
   static const Int_t fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]; // array with layers not to skip
-  static const Int_t fgkLastLayerToTrackTo;  // innermost layer
-  static const Int_t fgkMaxDetectorPerLayer; // max clusters per layer
-  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 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;             // 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
-  static const Double_t fgkSPDdetxlength;     // SPD ladder length in x
+  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)
 
 
   Int_t  fTracker;  // ITS tracker to be used (see AliITSReconstructor)
@@ -518,6 +526,7 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Float_t fVtxr3DPhiCutTight; // tight deltaPhi cut to define tracklets in vertexer 3D
   Float_t fVtxr3DDCACut;      // cut on tracklet-to-tracklet DCA in vertexer3D
   Int_t   fVtxr3DPileupAlgo;  // pileup algorithm (0 = VtxZ, 1 = 3D - 2 step, 2 = 3D all in once)
+  UChar_t fVtxr3DHighMultAlgo; // downscaling if 0 - traces if 1
 
   Int_t fLayersToSkip[AliITSgeomTGeo::kNLayers]; // array with layers to skip (MI,SA)
 
@@ -639,6 +648,7 @@ class AliITSRecoParam : public AliDetectorRecoParam
   Bool_t fSAOnePointTracks; // one-cluster tracks in SA (only for cosmics!)
   Bool_t fSAUseAllClusters; // do not skip clusters used by MI (same track twice in AliESDEvent!)
   Int_t fMaxSPDcontrForSAToUseAllClusters; // maximum nContr of SPD vertex for which trackerSA will reuse all ITS clusters
+  Bool_t fSAUsedEdxInfo;   // use/not use dE/dx in ITS for assign mass hypothesis
 
   Bool_t fSelectBestMIP03;          // (MI) Multiply norm chi2 by interpolated one in hypthesis analysis
   Bool_t fFlagFakes;                // (MI) preform shared cluster analysis and flag candidates for fakes
@@ -738,7 +748,7 @@ class AliITSRecoParam : public AliDetectorRecoParam
   AliITSRecoParam(const AliITSRecoParam & param);
   AliITSRecoParam & operator=(const AliITSRecoParam &param);
 
-  ClassDef(AliITSRecoParam,37) // ITS reco parameters
+  ClassDef(AliITSRecoParam,39) // ITS reco parameters
 };
 
 #endif