AliMUONRecoParam.h

Go to the documentation of this file.

#ifndef ALIMUONRECOPARAM_H
#define ALIMUONRECOPARAM_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice                               */

// $Id$

//  Author: Philippe Pillot

#include "AliDetectorRecoParam.h"
#include "TString.h"
#include <TMath.h>

class AliMUONRecoParam : public AliDetectorRecoParam
{
 public: 
  AliMUONRecoParam();
  virtual ~AliMUONRecoParam();
  
  static AliMUONRecoParam *GetLowFluxParam();
  static AliMUONRecoParam *GetHighFluxParam();
  static AliMUONRecoParam *GetCosmicParam();
  static AliMUONRecoParam *GetCalibrationParam();
  
  void      SetClusteringMode(Option_t* mode) {fClusteringMode = mode; fClusteringMode.ToUpper();}
  Option_t* GetClusteringMode() const {return fClusteringMode.Data();}
  
  Double_t AverageNoisePadCharge() const { return fAverageNoisePadCharge; }
  void AverageNoisePadCharge(Double_t noise) { fAverageNoisePadCharge = noise; }
  
  Double_t LowestPadCharge() const { return fChargeSigmaCut*fAverageNoisePadCharge; }

  Double_t ClusterChargeCut() const { return fClusterChargeCut; }
  void ClusterChargeCut(Double_t n) { fClusterChargeCut=n; }
  
  Double_t LowestClusterCharge() const { return ClusterChargeCut()*LowestPadCharge(); }
     
  void      SetTrackingMode(Option_t* mode) {fTrackingMode = mode; fTrackingMode.ToUpper();}
  Option_t* GetTrackingMode() const {return fTrackingMode.Data();}
  
  void      CombineClusterTrackReco(Bool_t flag) {fCombinedClusterTrackReco = flag;}
  Bool_t    CombineClusterTrackReco() const {return fCombinedClusterTrackReco;}
  
  void      SaveFullClusterInESD(Bool_t flag, Double_t percentOfEvent = 100.) {fSaveFullClusterInESD = flag;
                                 fPercentOfFullClusterInESD = (fSaveFullClusterInESD) ? percentOfEvent : 0.;}
  Bool_t    SaveFullClusterInESD() const {return fSaveFullClusterInESD;}
  Double_t  GetPercentOfFullClusterInESD() const {return fPercentOfFullClusterInESD;}
  
  void     SetMinBendingMomentum(Double_t val) {fMinBendingMomentum = val;}
  Double_t GetMinBendingMomentum() const {return fMinBendingMomentum;}
  void     SetMaxBendingMomentum(Double_t val) {fMaxBendingMomentum = val;}
  Double_t GetMaxBendingMomentum() const {return fMaxBendingMomentum;}
  
  void     SetMaxNonBendingSlope(Double_t val) {fMaxNonBendingSlope = val;}
  Double_t GetMaxNonBendingSlope() const {return fMaxNonBendingSlope;}
  void     SetMaxBendingSlope(Double_t val) {fMaxBendingSlope = val;}
  Double_t GetMaxBendingSlope() const {return fMaxBendingSlope;}
  
  void     SelectOnTrackSlope(Bool_t flag) {fSelectTrackOnSlope = flag;}
  Bool_t   SelectOnTrackSlope() const {return fSelectTrackOnSlope;}
  
  void     SetNonBendingVertexDispersion(Double_t val) {fNonBendingVertexDispersion = val;} 
  Double_t GetNonBendingVertexDispersion() const {return fNonBendingVertexDispersion;}
  void     SetBendingVertexDispersion(Double_t val) {fBendingVertexDispersion = val;} 
  Double_t GetBendingVertexDispersion() const {return fBendingVertexDispersion;}
  
  void     SetMaxNonBendingDistanceToTrack(Double_t val) {fMaxNonBendingDistanceToTrack = val;} 
  Double_t GetMaxNonBendingDistanceToTrack() const {return fMaxNonBendingDistanceToTrack;}
  void     SetMaxBendingDistanceToTrack(Double_t val) {fMaxBendingDistanceToTrack = val;} 
  Double_t GetMaxBendingDistanceToTrack() const {return fMaxBendingDistanceToTrack;}
  
  void     SetSigmaCutForTracking(Double_t val) {fSigmaCutForTracking = val;} 
  Double_t GetSigmaCutForTracking() const {return fSigmaCutForTracking;}
  
  void     ImproveTracks(Bool_t flag) {fImproveTracks = flag;} 
  void     ImproveTracks(Bool_t flag, Double_t sigmaCut) {fImproveTracks = flag; fSigmaCutForImprovement = sigmaCut;} 
  Bool_t   ImproveTracks() const {return fImproveTracks;}
  Double_t GetSigmaCutForImprovement() const {return fSigmaCutForImprovement;}
  
  void     SetSigmaCutForTrigger(Double_t val) {fSigmaCutForTrigger = val; fMaxNormChi2MatchTrigger = val*val;} 
  Double_t GetSigmaCutForTrigger() const {return fSigmaCutForTrigger;}
  void     SetStripCutForTrigger(Double_t val) {fStripCutForTrigger = val;} 
  Double_t GetStripCutForTrigger() const {return fStripCutForTrigger;}
  void     SetMaxStripAreaForTrigger(Double_t val) {fMaxStripAreaForTrigger = val;} 
  Double_t GetMaxStripAreaForTrigger() const {return fMaxStripAreaForTrigger;}
  
  Double_t GetMaxNormChi2MatchTrigger() const {return fMaxNormChi2MatchTrigger;}
  
  void     TrackAllTracks(Bool_t flag) {fTrackAllTracks = flag;} 
  Bool_t   TrackAllTracks() const {return fTrackAllTracks;}
  
  void     RecoverTracks(Bool_t flag) {fRecoverTracks = flag;} 
  Bool_t   RecoverTracks() const {return fRecoverTracks;}
  
  void     MakeTrackCandidatesFast(Bool_t flag) {fMakeTrackCandidatesFast = flag;} 
  Bool_t   MakeTrackCandidatesFast() const {return fMakeTrackCandidatesFast;}
  
  void     MakeMoreTrackCandidates(Bool_t flag) {fMakeMoreTrackCandidates = flag;} 
  Bool_t   MakeMoreTrackCandidates() const {return fMakeMoreTrackCandidates;}
  
  void     ComplementTracks(Bool_t flag) {fComplementTracks = flag;} 
  Bool_t   ComplementTracks() const {return fComplementTracks;}
  
  void     RemoveConnectedTracksInSt12(Bool_t flag) {fRemoveConnectedTracksInSt12 = flag;} 
  Bool_t   RemoveConnectedTracksInSt12() const {return fRemoveConnectedTracksInSt12;}
  
  void     UseSmoother(Bool_t flag) {fUseSmoother = flag;} 
  Bool_t   UseSmoother() const {return fUseSmoother;}
  
  void     UseChamber(Int_t iCh, Bool_t flag) {if (iCh >= 0 && iCh < 10) fUseChamber[iCh] = flag;}
  Bool_t   UseChamber(Int_t iCh) const {return (iCh >= 0 && iCh < 10) ? fUseChamber[iCh] : kFALSE;}
  
  void     RequestStation(Int_t iSt, Bool_t flag) {if (iSt >= 0 && iSt < 5) fRequestStation[iSt] = flag;}
  Bool_t   RequestStation(Int_t iSt) const {return (iSt >= 0 && iSt < 5) ? fRequestStation[iSt] : kFALSE;}
  UInt_t   RequestedStationMask() const;
  
  void   BypassSt45(Bool_t st4, Bool_t st5);
  
  Bool_t BypassSt45() const { return fBypassSt45==45; }
  
  Bool_t BypassSt4() const { return BypassSt45() || fBypassSt45==4 ; }
  
  Bool_t BypassSt5() const { return BypassSt45() || fBypassSt45==5 ; }
  
  void    SetHVLimit(Int_t chamberId, Double_t ichamber);
  Double_t HVLimit(Int_t chamberId) const;
  
  void    SetPedMeanLimits(float low, float high) { fPedMeanLimits[0]=low; fPedMeanLimits[1]=high; }
  Float_t PedMeanLowLimit() const { return fPedMeanLimits[0]; }
  Float_t PedMeanHighLimit() const { return fPedMeanLimits[1]; }
  
  void    SetPedSigmaLimits(float low, float high) { fPedSigmaLimits[0]=low; fPedSigmaLimits[1]=high; }
  Float_t PedSigmaLowLimit() const { return fPedSigmaLimits[0]; }
  Float_t PedSigmaHighLimit() const { return fPedSigmaLimits[1]; }
    
  void   SetPadGoodnessMask(UInt_t mask) { fPadGoodnessMask=mask; }
  UInt_t PadGoodnessMask() const { return fPadGoodnessMask; }
  
  Double_t ChargeSigmaCut() const { return fChargeSigmaCut; }
  
  void ChargeSigmaCut(Double_t value) { fChargeSigmaCut=value; }
  
  void     SetDefaultNonBendingReso(Int_t iCh, Double_t val) {if (iCh >= 0 && iCh < 10) fDefaultNonBendingReso[iCh] = val;}
  Double_t GetDefaultNonBendingReso(Int_t iCh) const {return (iCh >= 0 && iCh < 10) ? fDefaultNonBendingReso[iCh] : FLT_MAX;}
  void     SetDefaultBendingReso(Int_t iCh, Double_t val) {if (iCh >= 0 && iCh < 10) fDefaultBendingReso[iCh] = val;}
  Double_t GetDefaultBendingReso(Int_t iCh) const {return (iCh >= 0 && iCh < 10) ? fDefaultBendingReso[iCh] : FLT_MAX;}
  
  void SetMaxTriggerTracks(Int_t maxTriggerTracks) {fMaxTriggerTracks = maxTriggerTracks;}
  Int_t GetMaxTriggerTracks() const {return fMaxTriggerTracks;}
  
  void SetMaxTrackCandidates(Int_t maxTrackCandidates) {fMaxTrackCandidates = maxTrackCandidates;}
  Int_t GetMaxTrackCandidates() const {return fMaxTrackCandidates;}
  
  void SetManuOccupancyLimits(float low, float high) { fManuOccupancyLimits[0]=low; fManuOccupancyLimits[1]=high; }
  Float_t ManuOccupancyLowLimit() const { return fManuOccupancyLimits[0]; }
  Float_t ManuOccupancyHighLimit() const { return fManuOccupancyLimits[1]; }

  void SetBuspatchOccupancyLimits(float low, float high) { fBuspatchOccupancyLimits[0]=low; fBuspatchOccupancyLimits[1]=high; }
  Float_t BuspatchOccupancyLowLimit() const { return fBuspatchOccupancyLimits[0]; }
  Float_t BuspatchOccupancyHighLimit() const { return fBuspatchOccupancyLimits[1]; }

  void SetDEOccupancyLimits(float low, float high) { fDEOccupancyLimits[0]=low; fDEOccupancyLimits[1]=high; }
  Float_t DEOccupancyLowLimit() const { return fDEOccupancyLimits[0]; }
  Float_t DEOccupancyHighLimit() const { return fDEOccupancyLimits[1]; }
  
  void SetFractionOfBuspatchOutsideOccupancyLimit(float v) { fFractionOfBuspatchOutsideOccupancyLimit = v; }
  Float_t FractionOfBuspatchOutsideOccupancyLimit() const { return fFractionOfBuspatchOutsideOccupancyLimit; }

  virtual void Print(Option_t *option = "") const;
  
  virtual Double_t EventSizeSoftLimit() const { return fEventSizeSoftLimit; }
  
  virtual Double_t EventSizeHardLimit() const { return fEventSizeHardLimit; }

  virtual void SetEventSizeLimits(Double_t soft, Double_t hard) { fEventSizeSoftLimit=soft; fEventSizeHardLimit=hard; }
  
  virtual Double_t TokenLostLimit() const { return fTokenLostLimit; }

  virtual void SetTokenLostLimit(Double_t limit) { fTokenLostLimit = limit; }

  virtual Bool_t TryRecover() const { return fTryRecover; }

  virtual void TryRecover(Bool_t flag) { fTryRecover = flag; }

  void     DiscardMonoCathodClusters(Bool_t flag) { fDiscardMonoCathodClusters = flag; }
  void     DiscardMonoCathodClusters(Bool_t flag, Double_t resNB, Double_t resB) { fDiscardMonoCathodClusters = flag;
                                     fMonoCathodClNonBendingRes = resNB; fMonoCathodClBendingRes = resB; }
  Bool_t   DiscardMonoCathodClusters() const { return fDiscardMonoCathodClusters; }
  Double_t GetMonoCathodClNonBendingRes() const { return fMonoCathodClNonBendingRes; }
  Double_t GetMonoCathodClBendingRes() const { return fMonoCathodClBendingRes; }
  
  static TObjArray* Create(const char* settings);
  
  static void Show(Int_t runNumber, const char* ocdbPath="raw://");
  
private:
  
  void SetDefaultLimits();
  
 private:
  
  TString fClusteringMode; 
  
  TString fTrackingMode; 
  
  Double32_t fMinBendingMomentum; 
  Double32_t fMaxBendingMomentum; 
  Double32_t fMaxNonBendingSlope; 
  Double32_t fMaxBendingSlope;    
  
  Double32_t fNonBendingVertexDispersion; 
  Double32_t fBendingVertexDispersion;    
  
  Double32_t fMaxNonBendingDistanceToTrack; 
  Double32_t fMaxBendingDistanceToTrack;    
  
  Double32_t fSigmaCutForTracking; 

  Double32_t fSigmaCutForImprovement; 
  
  Double32_t fSigmaCutForTrigger; 

  Double32_t fStripCutForTrigger; 

  Double32_t fMaxStripAreaForTrigger; 
  
  Double32_t fMaxNormChi2MatchTrigger; 
  
  Double32_t fPercentOfFullClusterInESD; 
  
  Bool_t     fCombinedClusterTrackReco; 
  
  Bool_t     fTrackAllTracks; 
  
  Bool_t     fRecoverTracks; 
  
  Bool_t     fMakeTrackCandidatesFast; 
  
  Bool_t     fMakeMoreTrackCandidates; 
  
  Bool_t     fComplementTracks; 
  
  Bool_t     fImproveTracks; 
  
  Bool_t     fUseSmoother; 
  
  Bool_t     fSaveFullClusterInESD; 
  
  Int_t      fBypassSt45; 
  
  Bool_t     fUseChamber[10]; 
  
  Bool_t     fRequestStation[5]; 
  
  Double32_t fHVSt12Limits[2]; 
  Double32_t fHVSt345Limits[2]; 
  Double32_t fPedMeanLimits[2]; 
  Double32_t fPedSigmaLimits[2]; 
  
  UInt_t     fPadGoodnessMask; 
  
  Double32_t fChargeSigmaCut; 
  
  Double32_t fDefaultNonBendingReso[10]; 
  Double32_t fDefaultBendingReso[10]; 
  
  Bool_t     fRemoveConnectedTracksInSt12; 
  
  Int_t      fMaxTriggerTracks; 
  Int_t      fMaxTrackCandidates; 
  
  Bool_t     fSelectTrackOnSlope; 
  
  Double32_t fManuOccupancyLimits[2]; 
  Double32_t fBuspatchOccupancyLimits[2]; 
  Double32_t fDEOccupancyLimits[2]; 

  Double32_t fMissingPadFractionLimit; 
  Double32_t fFractionOfBuspatchOutsideOccupancyLimit; 

  Double32_t fAverageNoisePadCharge; 
  Double32_t fClusterChargeCut; 
  
  Double32_t fEventSizeSoftLimit; 
  Double32_t fEventSizeHardLimit; 
  
  Double32_t fTokenLostLimit; 
  
  Bool_t     fTryRecover; 

  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,171) // MUON reco parameters
  // we're at 171 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