/// get the clustering (pre-clustering) mode
Option_t* GetClusteringMode() const {return fClusteringMode.Data();}
+ /// Get the (truncated) average of sigmas of pedestal measurements, i.e. noise, of pads
+ Double_t AverageNoisePadCharge() const { return fAverageNoisePadCharge; }
+ /// Set the average of sigmas of pedestal measurements, i.e. noise, of pads
+ void AverageNoisePadCharge(Double_t noise) { fAverageNoisePadCharge = noise; }
+
+ /// Get the lowest charge we allow for pads
+ Double_t LowestPadCharge() const { return fChargeSigmaCut*fAverageNoisePadCharge; }
+
+ /// Get the cut applied to cut on cluster charge (the charge is cut if below fClusterChargeCut*LowestPadCharge())
+ Double_t ClusterChargeCut() const { return fClusterChargeCut; }
+ /// Set the cut applied to cut on cluster charge (the charge is cut if below fClusterChargeCut*LowestPadCharge())
+ void ClusterChargeCut(Double_t n) { fClusterChargeCut=n; }
+
+ /// Get the lowest possible cluster charge
+ Double_t LowestClusterCharge() const { return ClusterChargeCut()*LowestPadCharge(); }
+
/// set the tracking mode
void SetTrackingMode(Option_t* mode) {fTrackingMode = mode; fTrackingMode.ToUpper();}
/// get the tracking mode
/// return kTRUE/kFALSE if tracks are selected according to their slope/impact parameter
Bool_t SelectOnTrackSlope() const {return fSelectTrackOnSlope;}
- /// set the vertex dispersion (cm) in non bending plane (used for original tracking only)
+ /// set the vertex dispersion (cm) in non bending plane
void SetNonBendingVertexDispersion(Double_t val) {fNonBendingVertexDispersion = val;}
- /// return the vertex dispersion (cm) in non bending plane (used for original tracking only)
+ /// return the vertex dispersion (cm) in non bending plane
Double_t GetNonBendingVertexDispersion() const {return fNonBendingVertexDispersion;}
- /// set the vertex dispersion (cm) in bending plane (used for original tracking only)
+ /// set the vertex dispersion (cm) in bending plane
void SetBendingVertexDispersion(Double_t val) {fBendingVertexDispersion = val;}
- /// return the vertex dispersion (cm) in bending plane (used for original tracking only)
+ /// return the vertex dispersion (cm) in bending plane
Double_t GetBendingVertexDispersion() const {return fBendingVertexDispersion;}
/// set the maximum distance to the track to search for compatible cluster(s) in non bending direction
Double_t GetSigmaCutForImprovement() const {return fSigmaCutForImprovement;}
/// set the cut in sigma to apply on track during trigger hit pattern search
- void SetSigmaCutForTrigger(Double_t val) {fSigmaCutForTrigger = val;}
+ void SetSigmaCutForTrigger(Double_t val) {fSigmaCutForTrigger = val; fMaxNormChi2MatchTrigger = val*val;}
/// return the cut in sigma to apply on track during trigger hit pattern search
Double_t GetSigmaCutForTrigger() const {return fSigmaCutForTrigger;}
/// set the cut in strips to apply on trigger track during trigger chamber efficiency
/// return the maximum search area in strips to apply on trigger track during trigger chamber efficiency
Double_t GetMaxStripAreaForTrigger() const {return fMaxStripAreaForTrigger;}
- /// set the maximum normalized chi2 of tracking/trigger track matching
- void SetMaxNormChi2MatchTrigger(Double_t val) {fMaxNormChi2MatchTrigger = val;}
/// return the maximum normalized chi2 of tracking/trigger track matching
Double_t GetMaxNormChi2MatchTrigger() const {return fMaxNormChi2MatchTrigger;}
/// return kTRUE if we should replace clusters in St 5 by generated clusters from trigger tracks
Bool_t BypassSt5() const { return BypassSt45() || fBypassSt45==5 ; }
- /// Set Low and High threshold for St12 HV
- void SetHVSt12Limits(float low, float high) { fHVSt12Limits[0]=low; fHVSt12Limits[1]=high; }
- /// Retrieve low limit for St12's HV
- Float_t HVSt12LowLimit() const { return fHVSt12Limits[0]; }
- /// Retrieve high limit for St12's HV
- Float_t HVSt12HighLimit() const { return fHVSt12Limits[1]; }
-
- /// Set Low and High threshold for St345 HV
- void SetHVSt345Limits(float low, float high) { fHVSt345Limits[0]=low; fHVSt345Limits[1]=high; }
- /// Retrieve low limit for St345's HV
- Float_t HVSt345LowLimit() const { return fHVSt345Limits[0]; }
- /// Retrieve high limit for St345's HV
- Float_t HVSt345HighLimit() const { return fHVSt345Limits[1]; }
+ /// Set HV threshold for chambers (chamberId=0..9, use -1 to set all chambers equal)
+ void SetHVLimit(Int_t chamberId, Double_t ichamber);
+ /// Retrieve HV limit for chamber (chamberId=0..9)
+ Double_t HVLimit(Int_t chamberId) const;
/// Set Low and High threshold for pedestal mean
void SetPedMeanLimits(float low, float high) { fPedMeanLimits[0]=low; fPedMeanLimits[1]=high; }
/// Retrieve high value of DE occupancy limit
Float_t DEOccupancyHighLimit() const { return fDEOccupancyLimits[1]; }
- /// Set the missing pad fraction limit
- void SetMissingPadFractionLimit(float v) { fMissingPadFractionLimit = v; }
- /// Get the missing pad fraction limit
- Float_t MissingPadFractionLimit() const { return fMissingPadFractionLimit; }
-
/// Set the fraction of buspatches outside the occupancy limits
void SetFractionOfBuspatchOutsideOccupancyLimit(float v) { fFractionOfBuspatchOutsideOccupancyLimit = v; }
/// Get the fraction of buspatches outside the occupancy limits
virtual void Print(Option_t *option = "") const;
- private:
+ /// Get the max event size (soft limit)
+ virtual Double_t EventSizeSoftLimit() const { return fEventSizeSoftLimit; }
+
+ /// Get the max event size (hard limit)
+ virtual Double_t EventSizeHardLimit() const { return fEventSizeHardLimit; }
+
+ /// Set the max event size limits
+ virtual void SetEventSizeLimits(Double_t soft, Double_t hard) { fEventSizeSoftLimit=soft; fEventSizeHardLimit=hard; }
+
+ /// Get the percentage of token lost error we allow
+ virtual Double_t TokenLostLimit() const { return fTokenLostLimit; }
+
+ /// Set the percentage of token lost error we allow
+ virtual void SetTokenLostLimit(Double_t limit) { fTokenLostLimit = limit; }
+
+ /// Whether or not we try to recover corrupted raw data
+ virtual Bool_t TryRecover() const { return fTryRecover; }
+
+ /// Set the try recover corrupted raw data (use kTRUE only if you know what you are doing. Should be left to kFALSE by default)
+ virtual void TryRecover(Bool_t flag) { fTryRecover = flag; }
+
+ /// Discard or not the mono-cathod clusters by assigning to them different resolutions (use default values)
+ void DiscardMonoCathodClusters(Bool_t flag) { fDiscardMonoCathodClusters = flag; }
+ /// Discard or not the mono-cathod clusters by assigning to them different resolutions (use given values)
+ void DiscardMonoCathodClusters(Bool_t flag, Double_t resNB, Double_t resB) { fDiscardMonoCathodClusters = flag;
+ fMonoCathodClNonBendingRes = resNB; fMonoCathodClBendingRes = resB; }
+ /// Check whether to discard or not the mono-cathod clusters
+ Bool_t DiscardMonoCathodClusters() const { return fDiscardMonoCathodClusters; }
+ /// Get the non-bending resolution of mono-cathod clusters when the non-bending plane is missing
+ Double_t GetMonoCathodClNonBendingRes() const { return fMonoCathodClNonBendingRes; }
+ /// Get the bending resolution of mono-cathod clusters when the bending plane is missing
+ Double_t GetMonoCathodClBendingRes() const { return fMonoCathodClBendingRes; }
+
+ /// Create object ready to be put in OCDB
+ static TObjArray* Create(const char* settings);
+
+ /// Show what is the OCDB for that run
+ static void Show(Int_t runNumber, const char* ocdbPath="raw://");
+
+private:
void SetDefaultLimits();
Double32_t fGainA1Limits[2]; ///< Low and High threshold for gain a0 parameter
Double32_t fGainA2Limits[2]; ///< Low and High threshold for gain a1 parameter
Double32_t fGainThresLimits[2]; ///< Low and High threshold for gain threshold parameter
- Double32_t fHVSt12Limits[2]; ///< Low and High threshold for St12 HV
- Double32_t fHVSt345Limits[2]; ///< Low and High threshold for St345 HV
+ Double32_t fHVSt12Limits[2]; ///< DEPRECATED. See fHVLimits
+ Double32_t fHVSt345Limits[2]; ///< DEPRECATED. See fHVLimits
Double32_t fPedMeanLimits[2]; ///< Low and High threshold for pedestal mean
Double32_t fPedSigmaLimits[2]; ///< Low and High threshold for pedestal sigma
Double32_t fBuspatchOccupancyLimits[2]; ///< low and high thresholds for bus patch occupancy cut
Double32_t fDEOccupancyLimits[2]; ///< low and high thresholds for DE occupancy cut
- Double32_t fMissingPadFractionLimit; ///< above this fraction, we consider we have too few pads alive...
+ Double32_t fMissingPadFractionLimit; ///< DEPRECATED
Double32_t fFractionOfBuspatchOutsideOccupancyLimit; ///< above this limit, we consider we have too many buspatches out of the allowed occupancy range
+ Double32_t fAverageNoisePadCharge; ///< the (truncated, typically at 10%) mean of the sigma of the pedestals, in femto-coulomb
+ Double32_t fClusterChargeCut; ///< the cluster is cut if its charge is below fClusterChargeCut*LowestPadCharge()
+
+ Double32_t fEventSizeSoftLimit; ///< (soft) limit on mean event size per event (KB)
+ Double32_t fEventSizeHardLimit; ///< (hard) limit on mean event size per event (KB)
+
+ Double32_t fTokenLostLimit; ///< limit on the fraction of token lost error per event we allow
+
+ Bool_t fTryRecover; ///< try to recover corrupted raw data
+
+ Double32_t fHVLimit[10]; // HV limit (below which we consider that chamber efficiency is to be considered zero)
+
+ Double32_t fDiscardMonoCathodClusters; // assign a different resolution to mono-cathod clusters in the direction of the missing plane
+ Double32_t fMonoCathodClNonBendingRes; // resolution of mono-cathod clusters in the non-bending direction when the non-bending plane is missing
+ Double32_t fMonoCathodClBendingRes; // resolution of mono-cathod clusters in the bending direction when the bending plane is missing
+
// functions
void SetLowFluxParam();
void SetHighFluxParam();
void SetCosmicParam();
void SetCalibrationParam();
- ClassDef(AliMUONRecoParam,15) // MUON reco parameters
+ ClassDef(AliMUONRecoParam,170) // MUON reco parameters
+ // we're at 167 not because we had that many versions, but because at some point (version 15->16)
+ // 166 was committed by error, and we did not to go reverse afterwards...
};
#endif