[u/mrichter/AliRoot.git] / MUON / AliMUONConstants.h

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

/* $Id$*/
// Revision of includes 07/05/2004
//
/// \ingroup base
/// \class AliMUONConstants
/// \brief MUON global constants

#include <TObject.h>

class AliMUONConstants : public TObject {
 public:
    /// Return number of chambers
  static Int_t    NCh();
    /// Return number of tracking chambers
  static Int_t    NTrackingCh();
    /// Return number of tracking stations
    static Int_t    NTrackingSt() {return fgNTrackingSt;}
    /// Return number of trigger chambers
    static Int_t    NTriggerCh() {return NCh()-NTrackingCh();}
    /// Return number of trigger circuits
    static Int_t    NTriggerCircuit() {return fgNTriggerCircuit;}
    /// Return number of detection element
    static Int_t    NDetElem() {return fgNDetElem;}
    /// Return number of geometry modules
    static Int_t    NGeomModules() {return fgNGeomModules;}
    /// Return position of chamber i
    static Float_t  DefaultChamberZ(Int_t i) {return fgDefaultChamberZ[i];}
    /// Return ratio between trigger chambers
    static Float_t  DefaultRatioTriggerChamber(Int_t i) {return fgDefaultRatioTriggerChamber[i];}
    /// Return Inclination with respect the vertical axis of stations 345
    static Float_t  St345Inclination()       {return fgSt345inclination;}
    /// Return pointer to array of positions
    static Float_t* DefaultChamberZ()        {return fgDefaultChamberZ;}
    /// Return chamber i inner diameter
    static Float_t  Dmin(Int_t i)            {return fgDmin[i];}
    /// Return chamber i inner radius
    static Float_t  Rmin(Int_t i)            {return Dmin(i)/2.0;}
    /// Return chamber i outer diameter
    static Float_t  Dmax(Int_t i)            {return fgDmax[i];}
    /// Return chamber i outer radius
    static Float_t  Rmax(Int_t i)            {return Dmax(i)/2.0;}
    /// Return maximum zoom for event display
    static Int_t    MaxZoom()                {return fgMaxZoom;}
    /// Return half-distance between two half-chambers
    static Float_t  DzCh()                   {return fgDzCh;}
    /// Return half-distance between two slats
    static Float_t  DzSlat()                 {return fgDzSlat;}
    /// Return chamber number according z position of hit.
    static  Int_t ChamberNumber(Float_t z, bool warn = false); 
    /// Return SqrtKx3 for Slat
    static Float_t SqrtKx3()                 {return fgSqrtKx3;}
    /// Return SqrtKy3 for Slat
    static Float_t SqrtKy3()                 {return fgSqrtKy3;}
    /// Return SqrtKx3 for Station 1 & 2
    static Float_t SqrtKx3St1()              {return fgSqrtKx3St1;}
    /// Return SqrtKy3 for Station 1 & 2
    static Float_t SqrtKy3St1()              {return fgSqrtKy3St1;}
    /// Return charge correlation (needed for response and for cluster finder !?)
    static Float_t ChargeCorrel()            {return fgChargeCorrel;}
    /// Return charge correlation for Station 1 & 2 (needed for response and for cluster finder !?)
    static Float_t ChargeCorrelSt1()         {return fgChargeCorrelSt1;}
    /// Return wire pitch
    static Float_t Pitch()    {return fgPitch;}
    /// Return wire pitch for Station 1 & 2
    static Float_t PitchSt1() {return fgPitchSt1;}
    /// Return coil z-position
    static Double_t CoilZ() {return fgCoilZ;}
    /// Return coil lenght
    static Double_t CoilL() {return fgCoilL;}
    /// Return yoke z-position
    static Double_t YokeZ() {return fgYokeZ;}
    /// Return yoke lenght
    static Double_t YokeL() {return fgYokeL;}
    /// Return z-position of absorber begining
    static Double_t AbsZBeg() {return fgkAbsZBeg;}
    /// Return z-position of absorber end
    static Double_t AbsZEnd() {return fgkAbsZEnd;}
    /// Return chamber thickness in X0
    static Double_t ChamberThicknessInX0(Int_t chId) {return (chId >= 0 && chId < 10) ? fgChamberThicknessInX0[chId] : 0.;}
    /// Return Trigger ToF Limit (75 ns)
    static Float_t TriggerTofLimit() {return fgkTriggerTofLimit;}
    /// Return default trigger chamber resolution DUE TO ALIGNMENT ONLY in non bending direction
    static Double_t TriggerNonBendingReso() {return fgkTriggerNonBendingReso;}
    /// Return default trigger chamber resolution DUE TO ALIGNMENT ONLY in bending direction
    static Double_t TriggerBendingReso() {return fgkTriggerBendingReso;}
    /// Return z-position of muon filter begining
    static Double_t MuonFilterZBeg() {return fgkMuonFilterZBeg;}
    /// Return z-position of muon filter end
    static Double_t MuonFilterZEnd() {return fgkMuonFilterZEnd;}
    /// Return radiation length (in cm) of muon filter
    static Double_t MuonFilterX0() {return fgkMuonFilterX0;}

    /// Return the most probable bending momentum (GeV/c) (used when B = 0)
    static Double_t GetMostProbBendingMomentum() { return fgkMostProbBendingMomentum; }
      
    /// Return average arrival time to chamber i
    static Float_t  AverageChamberT(Int_t i) {return fgAverageChamberT[i];}
    /// Return a reduced total charge
    static Float_t  ReducedQTot(Float_t qtot, Float_t timeDif);
  
 protected:
    /// Default constructor
    AliMUONConstants() : TObject() {}
    /// Destructor
    virtual ~AliMUONConstants(){}

 private:
    static Int_t  fgNTrackingSt;        ///<  Number of Tracking Stations
    static Int_t  fgNTriggerCircuit;    ///<  Number of Trigger Circuits
    static Int_t  fgNDetElem;           ///<  Number of Detection Elements.
    static Int_t  fgNGeomModules;       ///<  Number of Geometry modules   

    static Float_t  fgDefaultChamberZ[14];		//!< Z-positions of chambers
    static Float_t  fgDefaultChamberMinZ[14];		//!< Minimum z-positions of chambers
    static Float_t  fgDefaultChamberMaxZ[14];		//!< Maximum z-positions of chambers
    static Float_t  fgDefaultRatioTriggerChamber[4];	///< Ratio between trigger chambers
    static Float_t  fgSt345inclination;			//!< Inclination with respect the vertical axis of stations 345
    static Float_t  fgDmin[7];				//!< Inner diameter
    static Float_t  fgDmax[7];				//!< Outer diameter

    static Float_t  fgDzCh;             ///< Half-distance between two half-chambers 
    static Float_t  fgDzSlat;           ///< Half-distance between two slat on the same chamber
    static Float_t  fgSqrtKx3;          ///< SqrtKx3 for St2 & Slat
    static Float_t  fgSqrtKy3;          ///< SqrtKy3 for St2 & Slat
    static Float_t  fgSqrtKx3St1;       ///< SqrtKx3 for Station 1 
    static Float_t  fgSqrtKy3St1;       ///< SqrtKy3 for Station 1 
 
    static Float_t  fgChargeCorrel;      ///< Charge correlation for St2 & Slats
    static Float_t  fgChargeCorrelSt1;   ///< Charge correlation for Station 1

    static Float_t  fgPitch;             ///< Wire pitch for St2 & Slats
    static Float_t  fgPitchSt1;          ///< Wire pitch for Station 1

    static Double_t  fgChamberThicknessInX0[10]; ///< default chamber thickness in X0 for reconstruction
    
    static Double_t fgCoilZ; ///< Coil z-position
    static Double_t fgCoilL; ///< Coil lenght
    static Double_t fgYokeZ; ///< Yoke z-position
    static Double_t fgYokeL; ///< Yoke lenght

    static const Double_t fgkAbsZBeg; ///< z-position of the begining of the absorber
    static const Double_t fgkAbsZEnd; ///< z-position of the end of the absorber
    
    static Int_t    fgMaxZoom;          ///< Maximum Zoom for event display
    static Float_t  fgkTriggerTofLimit; ///< Particle above this threshold are discarded in trigger algorithm
    
    static const Double_t fgkTriggerNonBendingReso; ///< default trigger chamber resolution in the non-bending direction, DUE TO ALIGNMENT ONLY
    static const Double_t fgkTriggerBendingReso; ///< default trigger chamber resolution in the bending direction, DUE TO ALIGNMENT ONLY 
    
    static const Double_t fgkMuonFilterZBeg; ///< z-position of the begining of the muon filter
    static const Double_t fgkMuonFilterZEnd; ///< z-position of the end of the muon filter
    static const Double_t fgkMuonFilterX0; ///< x0 of the muon filter
    
    static const Double_t fgkMostProbBendingMomentum; ///< most probable value (GeV/c) of muon momentum in bending plane (used when B = 0) needed to get some "reasonable" corrections for MCS and E loss even if B = 0

    static Float_t fgAverageChamberT[14]; ///<average arrival time to chamber
  
    ClassDef(AliMUONConstants, 0) // MUON global constants 
};
	
#endif
Commit	Line	Data
	1	#ifndef ALIMUONCONSTANTS_H
	2	#define ALIMUONCONSTANTS_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id$*/
	7	// Revision of includes 07/05/2004
	8	//
	9	/// \ingroup base
	10	/// \class AliMUONConstants
	11	/// \brief MUON global constants
	12
	13	#include <TObject.h>
	14
	15	class AliMUONConstants : public TObject {
	16	public:
	17	/// Return number of chambers
	18	static Int_t NCh();
	19	/// Return number of tracking chambers
	20	static Int_t NTrackingCh();
	21	/// Return number of tracking stations
	22	static Int_t NTrackingSt() {return fgNTrackingSt;}
	23	/// Return number of trigger chambers
	24	static Int_t NTriggerCh() {return NCh()-NTrackingCh();}
	25	/// Return number of trigger circuits
	26	static Int_t NTriggerCircuit() {return fgNTriggerCircuit;}
	27	/// Return number of detection element
	28	static Int_t NDetElem() {return fgNDetElem;}
	29	/// Return number of geometry modules
	30	static Int_t NGeomModules() {return fgNGeomModules;}
	31	/// Return position of chamber i
	32	static Float_t DefaultChamberZ(Int_t i) {return fgDefaultChamberZ[i];}
	33	/// Return ratio between trigger chambers
	34	static Float_t DefaultRatioTriggerChamber(Int_t i) {return fgDefaultRatioTriggerChamber[i];}
	35	/// Return Inclination with respect the vertical axis of stations 345
	36	static Float_t St345Inclination() {return fgSt345inclination;}
	37	/// Return pointer to array of positions
	38	static Float_t* DefaultChamberZ() {return fgDefaultChamberZ;}
	39	/// Return chamber i inner diameter
	40	static Float_t Dmin(Int_t i) {return fgDmin[i];}
	41	/// Return chamber i inner radius
	42	static Float_t Rmin(Int_t i) {return Dmin(i)/2.0;}
	43	/// Return chamber i outer diameter
	44	static Float_t Dmax(Int_t i) {return fgDmax[i];}
	45	/// Return chamber i outer radius
	46	static Float_t Rmax(Int_t i) {return Dmax(i)/2.0;}
	47	/// Return maximum zoom for event display
	48	static Int_t MaxZoom() {return fgMaxZoom;}
	49	/// Return half-distance between two half-chambers
	50	static Float_t DzCh() {return fgDzCh;}
	51	/// Return half-distance between two slats
	52	static Float_t DzSlat() {return fgDzSlat;}
	53	/// Return chamber number according z position of hit.
	54	static Int_t ChamberNumber(Float_t z, bool warn = false);
	55	/// Return SqrtKx3 for Slat
	56	static Float_t SqrtKx3() {return fgSqrtKx3;}
	57	/// Return SqrtKy3 for Slat
	58	static Float_t SqrtKy3() {return fgSqrtKy3;}
	59	/// Return SqrtKx3 for Station 1 & 2
	60	static Float_t SqrtKx3St1() {return fgSqrtKx3St1;}
	61	/// Return SqrtKy3 for Station 1 & 2
	62	static Float_t SqrtKy3St1() {return fgSqrtKy3St1;}
	63	/// Return charge correlation (needed for response and for cluster finder !?)
	64	static Float_t ChargeCorrel() {return fgChargeCorrel;}
	65	/// Return charge correlation for Station 1 & 2 (needed for response and for cluster finder !?)
	66	static Float_t ChargeCorrelSt1() {return fgChargeCorrelSt1;}
	67	/// Return wire pitch
	68	static Float_t Pitch() {return fgPitch;}
	69	/// Return wire pitch for Station 1 & 2
	70	static Float_t PitchSt1() {return fgPitchSt1;}
	71	/// Return coil z-position
	72	static Double_t CoilZ() {return fgCoilZ;}
	73	/// Return coil lenght
	74	static Double_t CoilL() {return fgCoilL;}
	75	/// Return yoke z-position
	76	static Double_t YokeZ() {return fgYokeZ;}
	77	/// Return yoke lenght
	78	static Double_t YokeL() {return fgYokeL;}
	79	/// Return z-position of absorber begining
	80	static Double_t AbsZBeg() {return fgkAbsZBeg;}
	81	/// Return z-position of absorber end
	82	static Double_t AbsZEnd() {return fgkAbsZEnd;}
	83	/// Return chamber thickness in X0
	84	static Double_t ChamberThicknessInX0(Int_t chId) {return (chId >= 0 && chId < 10) ? fgChamberThicknessInX0[chId] : 0.;}
	85	/// Return Trigger ToF Limit (75 ns)
	86	static Float_t TriggerTofLimit() {return fgkTriggerTofLimit;}
	87	/// Return default trigger chamber resolution DUE TO ALIGNMENT ONLY in non bending direction
	88	static Double_t TriggerNonBendingReso() {return fgkTriggerNonBendingReso;}
	89	/// Return default trigger chamber resolution DUE TO ALIGNMENT ONLY in bending direction
	90	static Double_t TriggerBendingReso() {return fgkTriggerBendingReso;}
	91	/// Return z-position of muon filter begining
	92	static Double_t MuonFilterZBeg() {return fgkMuonFilterZBeg;}
	93	/// Return z-position of muon filter end
	94	static Double_t MuonFilterZEnd() {return fgkMuonFilterZEnd;}
	95	/// Return radiation length (in cm) of muon filter
	96	static Double_t MuonFilterX0() {return fgkMuonFilterX0;}
	97
	98	/// Return the most probable bending momentum (GeV/c) (used when B = 0)
	99	static Double_t GetMostProbBendingMomentum() { return fgkMostProbBendingMomentum; }
	100
	101	/// Return average arrival time to chamber i
	102	static Float_t AverageChamberT(Int_t i) {return fgAverageChamberT[i];}
	103	/// Return a reduced total charge
	104	static Float_t ReducedQTot(Float_t qtot, Float_t timeDif);
	105
	106	protected:
	107	/// Default constructor
	108	AliMUONConstants() : TObject() {}
	109	/// Destructor
	110	virtual ~AliMUONConstants(){}
	111
	112	private:
	113	static Int_t fgNTrackingSt; ///< Number of Tracking Stations
	114	static Int_t fgNTriggerCircuit; ///< Number of Trigger Circuits
	115	static Int_t fgNDetElem; ///< Number of Detection Elements.
	116	static Int_t fgNGeomModules; ///< Number of Geometry modules
	117
	118	static Float_t fgDefaultChamberZ[14]; //!< Z-positions of chambers
	119	static Float_t fgDefaultChamberMinZ[14]; //!< Minimum z-positions of chambers
	120	static Float_t fgDefaultChamberMaxZ[14]; //!< Maximum z-positions of chambers
	121	static Float_t fgDefaultRatioTriggerChamber[4]; ///< Ratio between trigger chambers
	122	static Float_t fgSt345inclination; //!< Inclination with respect the vertical axis of stations 345
	123	static Float_t fgDmin[7]; //!< Inner diameter
	124	static Float_t fgDmax[7]; //!< Outer diameter
	125
	126	static Float_t fgDzCh; ///< Half-distance between two half-chambers
	127	static Float_t fgDzSlat; ///< Half-distance between two slat on the same chamber
	128	static Float_t fgSqrtKx3; ///< SqrtKx3 for St2 & Slat
	129	static Float_t fgSqrtKy3; ///< SqrtKy3 for St2 & Slat
	130	static Float_t fgSqrtKx3St1; ///< SqrtKx3 for Station 1
	131	static Float_t fgSqrtKy3St1; ///< SqrtKy3 for Station 1
	132
	133	static Float_t fgChargeCorrel; ///< Charge correlation for St2 & Slats
	134	static Float_t fgChargeCorrelSt1; ///< Charge correlation for Station 1
	135
	136	static Float_t fgPitch; ///< Wire pitch for St2 & Slats
	137	static Float_t fgPitchSt1; ///< Wire pitch for Station 1
	138
	139	static Double_t fgChamberThicknessInX0[10]; ///< default chamber thickness in X0 for reconstruction
	140
	141	static Double_t fgCoilZ; ///< Coil z-position
	142	static Double_t fgCoilL; ///< Coil lenght
	143	static Double_t fgYokeZ; ///< Yoke z-position
	144	static Double_t fgYokeL; ///< Yoke lenght
	145
	146	static const Double_t fgkAbsZBeg; ///< z-position of the begining of the absorber
	147	static const Double_t fgkAbsZEnd; ///< z-position of the end of the absorber
	148
	149	static Int_t fgMaxZoom; ///< Maximum Zoom for event display
	150	static Float_t fgkTriggerTofLimit; ///< Particle above this threshold are discarded in trigger algorithm
	151
	152	static const Double_t fgkTriggerNonBendingReso; ///< default trigger chamber resolution in the non-bending direction, DUE TO ALIGNMENT ONLY
	153	static const Double_t fgkTriggerBendingReso; ///< default trigger chamber resolution in the bending direction, DUE TO ALIGNMENT ONLY
	154
	155	static const Double_t fgkMuonFilterZBeg; ///< z-position of the begining of the muon filter
	156	static const Double_t fgkMuonFilterZEnd; ///< z-position of the end of the muon filter
	157	static const Double_t fgkMuonFilterX0; ///< x0 of the muon filter
	158
	159	static const Double_t fgkMostProbBendingMomentum; ///< most probable value (GeV/c) of muon momentum in bending plane (used when B = 0) needed to get some "reasonable" corrections for MCS and E loss even if B = 0
	160
	161	static Float_t fgAverageChamberT[14]; ///<average arrival time to chamber
	162
	163	ClassDef(AliMUONConstants, 0) // MUON global constants
	164	};
	165
	166	#endif
	167