#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; }
SetVertexer3DDeltaPhiCuts();
SetVertexer3DDCACut();
}
+ void SetSPDVertexerPileupAlgoOff(){fVtxr3DPileupAlgo=3;}
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]; }
Double_t GetMaxRoad() const { return fMaxRoad; }
Double_t GetMaxNormChi2ForGolden(Int_t i) const { return 3.+0.5*i; }
+ void SetSearchForExtraClusters(Bool_t opt=kTRUE){ fSearchForExtras=opt; }
+ Double_t GetSearchForExtraClusters() const { return fSearchForExtras; }
+
Double_t GetXVdef() const { return fXV; }
Double_t GetYVdef() const { return fYV; }
Double_t GetZVdef() const { return fZV; }
void SetNSigZFromBoundaryPlaneEff(Double_t nsigz=1.) {fNSigZFromBoundaryPlaneEff=nsigz;}
Double_t GetNSigZFromBoundaryPlaneEff() const {return fNSigZFromBoundaryPlaneEff;}
//
+ void SetImproveWithVertex(Bool_t impr=kFALSE) { fImproveWithVertex=impr; return; }
+ Bool_t GetImproveWithVertex() const { return fImproveWithVertex; }
void SetExtendedEtaAcceptance(Bool_t ext=kTRUE) { fExtendedEtaAcceptance=ext; return; }
Bool_t GetExtendedEtaAcceptance() const { return fExtendedEtaAcceptance; }
void SetAllowProlongationWithEmptyRoad(Bool_t allow=kTRUE) { fAllowProlongationWithEmptyRoad=allow; return; }
void SetSAUseAllClusters(Bool_t opt=kTRUE) { fSAUseAllClusters=opt; return; }
Bool_t GetSAUseAllClusters() const { return fSAUseAllClusters; }
+ 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; }
void SetTrackleterZetaOverlapCut(Float_t w=0.05) {fTrackleterZetaOverlapCut=w;}
Float_t GetTrackleterPhiOverlapCut() const {return fTrackleterPhiOverlapCut;}
Float_t GetTrackleterZetaOverlapCut() const {return fTrackleterZetaOverlapCut;}
-
+ void SetTrackleterPhiRotationAngle(Float_t w=0.0) {fTrackleterPhiRotationAngle=w;}
+ Float_t GetTrackleterPhiRotationAngle() const {return fTrackleterPhiRotationAngle;}
+ //
+ void SetTrackleterNStdDevCut(Float_t f=1.) {fTrackleterNStdDev = f<0.01 ? 0.01 : f;}
+ Float_t GetTrackleterNStdDevCut() const {return fTrackleterNStdDev;}
+ void SetTrackleterScaleDThetaBySin2T(Bool_t v=kFALSE) {fScaleDTBySin2T = v;}
+ Bool_t GetTrackleterScaleDThetaBySin2T() const {return fScaleDTBySin2T;}
//
void SetSPDRemoveNoisyFlag(Bool_t value) {fSPDRemoveNoisyFlag = value;}
Bool_t GetSPDRemoveNoisyFlag() const {return fSPDRemoveNoisyFlag;}
//
AliESDV0Params *GetESDV0Params() const {return fESDV0Params;}
//
- 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
+ // Lorentz angle
+ Bool_t GetCorrectLorentzAngleSPD() const {return fCorrectLorentzAngleSPD;}
+ Float_t GetTanLorentzAngleHolesSPD() const {return fTanLorentzAngleHolesSPD;}
+ Bool_t GetCorrectLorentzAngleSSD() const {return fCorrectLorentzAngleSSD;}
+ Float_t GetTanLorentzAngleHolesSSD() const {return fTanLorentzAngleHolesSSD;}
+ Float_t GetTanLorentzAngleElectronsSSD() const {return fTanLorentzAngleElectronsSSD;}
+
+ void SetCorrectLorentzAngleSPD(Bool_t flag) {fCorrectLorentzAngleSPD=flag;}
+ void SetTanLorentzAngleHolesSPD(Float_t la) {fTanLorentzAngleHolesSPD=la;}
+ void SetCorrectLorentzAngleSSD(Bool_t flag) {fCorrectLorentzAngleSSD=flag;}
+ void SetTanLorentzAngleHolesSSD(Float_t la) {fTanLorentzAngleHolesSSD=la;}
+ void SetTanLorentzAngleElectronsSSD(Float_t la) {fTanLorentzAngleElectronsSSD=la;}
+
+ //
+ 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)
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)
Double_t fMaxChi2In; // (NOT USED)
Double_t fMaxChi2sR[AliITSgeomTGeo::kNLayers]; // (NOT USED)
Double_t fChi2PerCluster; // (NOT USED)
+ // search for extra clusters
+ Bool_t fSearchForExtras; // swicth yes/no for the search of extra-clusters in RefitInward step
//
// default primary vertex (MI,V2)
Double_t fXV; // x
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
+ 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 fUseBadZonesFromOCDB; // enable using OCDB info on dead modules and chips (MI)
Bool_t fUseSingleBadChannelsFromOCDB; // enable using OCDB info on bad single SPD pixels and SDD anodes (MI)
Float_t fMinClusterChargeSA; // minimum SDD,SSD cluster charge for SA tarcker
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)
Bool_t fTrackleterRemoveClustersFromOverlaps; // Option to skip clusters in the overlaps
Float_t fTrackleterPhiOverlapCut; // Fiducial window in phi for overlap cut
Float_t fTrackleterZetaOverlapCut; // Fiducial window in eta for overlap cut
+ 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 fUseCosmicRunShiftsSSD; // SSD time shifts for cosmic run 2007/2008 (use for data taken up to 18 sept 2008)
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
+
private:
AliESDV0Params * fESDV0Params; // declare the AliESDV0Params to be able to used in AliITSV0Finder
AliITSRecoParam(const AliITSRecoParam & param);
AliITSRecoParam & operator=(const AliITSRecoParam ¶m);
- ClassDef(AliITSRecoParam,29) // ITS reco parameters
+ ClassDef(AliITSRecoParam,39) // ITS reco parameters
};
#endif
-