void SetClusterFinder(Int_t cf=0) { fClusterFinder=cf; }
void SetClusterFinderV2() { SetClusterFinder(0); }
- void SetClusterFinderOrig() { SetClusterFinder(1); }
+ void SetClusterFinderSDDfast() { SetClusterFinder(1); }
Int_t GetClusterFinder() const { return fClusterFinder; }
void SetPID(Int_t pid=0) {fPID=pid;}
void SetDefaultPID() {SetPID(0);}
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;}
Bool_t GetRequireClusterInInnerLayerPlaneEff() const {return fRequireClusterInInnerLayerPlaneEff;}
void SetOnlyConstraintPlaneEff(Bool_t con=kFALSE) { fOnlyConstraintPlaneEff=con; }
Bool_t GetOnlyConstraintPlaneEff() const { return fOnlyConstraintPlaneEff; }
- void SetNSigXFromBoundaryPlaneEff(Double_t nsigx=1.) {fNSigXFromBoundaryPlaneEff=nsigx;}
+ void SetNSigXFromBoundaryPlaneEff(Double_t nsigx=0.) {if(nsigx<0.)fNSigXFromBoundaryPlaneEff=TMath::Abs(nsigx);else fNSigXFromBoundaryPlaneEff=nsigx;}
Double_t GetNSigXFromBoundaryPlaneEff() const {return fNSigXFromBoundaryPlaneEff;}
- void SetNSigZFromBoundaryPlaneEff(Double_t nsigz=1.) {fNSigZFromBoundaryPlaneEff=nsigz;}
+ void SetNSigZFromBoundaryPlaneEff(Double_t nsigz=0.) {if(nsigz<0.)fNSigZFromBoundaryPlaneEff=TMath::Abs(nsigz);else fNSigZFromBoundaryPlaneEff=nsigz;}
Double_t GetNSigZFromBoundaryPlaneEff() const {return fNSigZFromBoundaryPlaneEff;}
+ void SetDistXFromBoundaryPlaneEff(Double_t distx=0.) {if(distx<0.)fDistXFromBoundaryPlaneEff=TMath::Abs(distx);else fDistXFromBoundaryPlaneEff=distx;}
+ Double_t GetDistXFromBoundaryPlaneEff() const {return fDistXFromBoundaryPlaneEff;}
+ void SetDistZFromBoundaryPlaneEff(Double_t distz=0.) {if(distz<0.)fDistZFromBoundaryPlaneEff=TMath::Abs(distz);else fDistZFromBoundaryPlaneEff=distz;}
+ Double_t GetDistZFromBoundaryPlaneEff() const {return fDistZFromBoundaryPlaneEff;}
+ void SetSwitchOnMaxDistNSigFrmBndPlaneEff(Bool_t flagbnd=kFALSE) {fSwitchOnMaxDistNSigFrmBndPlaneEff=flagbnd;}
+ Bool_t GetSwitchOnMaxDistNSigFrmBndPlaneEff() const {return fSwitchOnMaxDistNSigFrmBndPlaneEff;}
+ void SetSwitchOffStdSearchClusPlaneEff(Bool_t flagstdclus=kFALSE) {fSwitchOffStdSearchClusPlaneEff=flagstdclus;}
+ Bool_t GetSwitchOffStdSearchClusPlaneEff() const {return fSwitchOffStdSearchClusPlaneEff;}
+ void SetNSigXSearchClusterPlaneEff(Double_t nsigclx=0.) {if(nsigclx<0.)fNSigXSearchClusterPlaneEff=TMath::Abs(nsigclx);else fNSigXSearchClusterPlaneEff=nsigclx;}
+ Double_t GetNSigXSearchClusterPlaneEff() const {return fNSigXSearchClusterPlaneEff;}
+ void SetNSigZSearchClusterPlaneEff(Double_t nsigclz=0.) {if(nsigclz<0.)fNSigZSearchClusterPlaneEff=TMath::Abs(nsigclz);else fNSigZSearchClusterPlaneEff=nsigclz;}
+ Double_t GetNSigZSearchClusterPlaneEff() const {return fNSigZSearchClusterPlaneEff;}
+ void SetDistXSearchClusterPlaneEff(Double_t distclx=0.) {if(distclx<0.)fDistXSearchClusterPlaneEff=TMath::Abs(distclx);else fDistXSearchClusterPlaneEff=distclx;}
+ Double_t GetDistXSearchClusterPlaneEff() const {return fDistXSearchClusterPlaneEff;}
+ void SetDistZSearchClusterPlaneEff(Double_t distclz=0.) {if(distclz<0.)fDistZSearchClusterPlaneEff=TMath::Abs(distclz);else fDistZSearchClusterPlaneEff=distclz;}
+ Double_t GetDistZSearchClusterPlaneEff() const {return fDistZSearchClusterPlaneEff;}
+ void SetSwitchOnMaxDistNSigSrhClusPlaneEff(Bool_t flagbndcl=kFALSE) {fSwitchOnMaxDistNSigSrhClusPlaneEff=flagbndcl;}
+ Bool_t GetSwitchOnMaxDistNSigSrhClusPlaneEff() const {return fSwitchOnMaxDistNSigSrhClusPlaneEff;}
+ 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; }
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 {kMaxClusterPerLayer=70000}; //7000*10; // max clusters per layer
enum {kMaxClusterPerLayer5=28000};//7000*10*2/5; // max clusters per layer
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)
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 fDistXFromBoundaryPlaneEff; // accept one track for PlaneEff if distance from border (in loc x or z)
+ Double_t fDistZFromBoundaryPlaneEff; // is greater than fDistXFromBoundaryPlaneEff centimeters
+ Bool_t fSwitchOnMaxDistNSigFrmBndPlaneEff; //if kTRUE,use max(fDistXFromBoundaryPlaneEff,fNSigXFromBoundaryPlaneEff) to accept tracks
+ Bool_t fSwitchOffStdSearchClusPlaneEff; //if kTRUE,use fNSigXSearchClusterPlaneEff and fDistXSearchClusterPlaneEff
+ Double_t fNSigXSearchClusterPlaneEff; // cluster search in distance from track impact point (in loc x or z)
+ Double_t fNSigZSearchClusterPlaneEff; // less than fNSigXSearchClusterPlaneEff * Track_precision
+ Double_t fDistXSearchClusterPlaneEff; // cluster found in distance from track impact point (in loc x or z)
+ Double_t fDistZSearchClusterPlaneEff; // is greater than fDistXSearchClusterPlaneEff centimeters
+ Bool_t fSwitchOnMaxDistNSigSrhClusPlaneEff; //if kTRUE,use max(fDistXSearchClusterPlaneEff,fNSigXSearchClusterPlaneEff) to accept tracks
+ 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)
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,39) // ITS reco parameters
+ ClassDef(AliITSRecoParam,53) // ITS reco parameters
};
#endif
+
+