#include "AliDetectorRecoParam.h"
#include "AliITSgeomTGeo.h"
-#include "AliESDV0Params.h"
+//#include "AliESDV0Params.h"
+
+class AliESDV0Params;
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; }
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;}
+ Bool_t GetUseImproveKalman() const {return fUseImproveKalman;}
+ void SetSelectBestMIP03(Bool_t v=kTRUE) {fSelectBestMIP03 = v;}
+ void SetFlagFakes(Bool_t v=kTRUE) {fFlagFakes = v;}
+ void SetUseImproveKalman(Bool_t v=kTRUE) {fUseImproveKalman = v;}
+ //
Float_t GetVertexer3DWideFiducialRegionZ() const {return fVtxr3DZCutWide;}
Float_t GetVertexer3DWideFiducialRegionR() const {return fVtxr3DRCutWide;}
Float_t GetVertexer3DNarrowFiducialRegionZ() const {return fVtxr3DZCutNarrow;}
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]; }
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;}
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; }
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; }
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;}
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 {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 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 Int_t fgkLastLayerToTrackTo=0; // innermost layer
+ static const Int_t fgkMaxDetectorPerLayer=1000; // 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 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)
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)
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)
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
+ Bool_t fUseImproveKalman; // (MI) Use ImproveKalman version of AliITSTrackV2 instead of Improve
Bool_t fFindV0s; // flag to enable V0 finder (MI)
Bool_t fStoreLikeSignV0s; // flag to store like-sign V0s (MI)
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)
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
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
AliITSRecoParam(const AliITSRecoParam & param);
AliITSRecoParam & operator=(const AliITSRecoParam ¶m);
- ClassDef(AliITSRecoParam,36) // ITS reco parameters
+ ClassDef(AliITSRecoParam,52) // ITS reco parameters
};
#endif
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff