TRD/AliTRDgtuParam.h

   1 #ifndef ALITRDGTUPARAM_H
   2 #define ALITRDGTUPARAM_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id: AliTRDgtuParam.h 27496 2008-07-22 08:35:45Z cblume $ */
   7
   8 // --------------------------------------------------------
   9 //
  10 // Singleton class to hold the parameters steering the GTU
  11 // tracking
  12 //
  13 // --------------------------------------------------------
  14
  15 #include "TObject.h"
  16
  17 class AliTRDgeometry;
  18
  19 class AliTRDgtuParam : public TObject {
  20  public:
  21   virtual ~AliTRDgtuParam();
  22
  23   static AliTRDgtuParam *Instance(); // Singleton
  24   static void Terminate();
  25
  26   static Int_t GetNLinks() { return fgkNLinks; }
  27   static Int_t GetNLayers() { return fgkNLinks/2; }
  28   static Int_t GetNZChannels() { return fgkNZChannels; }
  29   static Int_t GetNRefLayers() { return fgkNRefLayers; }
  30
  31   static Float_t GetChamberThickness() { return 3.0; }
  32
  33   // ----- Bin widths (granularity) -----
  34   static Float_t GetBinWidthY() { return fgkBinWidthY; }
  35   static Float_t GetBinWidthdY() { return fgkBinWidthdY; }
  36
  37   // ----- Bit Widths (used for internal representation) -----
  38   static Int_t GetBitWidthY() { return fgkBitWidthY; }
  39   static Int_t GetBitWidthdY() { return fgkBitWidthdY; }
  40   static Int_t GetBitWidthYProj() { return fgkBitWidthYProj; }
  41   static Int_t GetBitExcessY() { return fgkBitExcessY; }
  42   static Int_t GetBitExcessAlpha() { return fgkBitExcessAlpha; }
  43   static Int_t GetBitExcessYProj() { return fgkBitExcessYProj; }
  44
  45   AliTRDgeometry* GetGeo() const { return fGeo; }
  46   Float_t GetVertexSize() const { return fVertexSize; }
  47   Int_t GetCiAlpha(Int_t layer) const;
  48   Int_t GetCiYProj(Int_t layer) const;
  49   Int_t GetYt(Int_t stack, Int_t layer, Int_t zrow) const;
  50   Int_t GetDeltaY() const { return fgkDeltaY; }
  51   Int_t GetDeltaAlpha() const { return fgkDeltaAlpha; }
  52   Int_t GetZSubchannel(Int_t stack, Int_t layer, Int_t zchannel, Int_t zpos) const;
  53   Int_t GetRefLayer(Int_t refLayerIdx) const;
  54 //  Bool_t GetFitParams(TVectorD &rhs, Int_t k); // const
  55   Bool_t GetIntersectionPoints(Int_t k, Float_t &x1, Float_t &x2); // const
  56   Float_t GetPt(Int_t a, Float_t b, Float_t x1, Float_t x2) const;
  57
  58   Bool_t IsInZChannel(Int_t stack, Int_t layer, Int_t zchannel, Int_t zpos) const;
  59
  60   void SetVertexSize(Float_t vertexsize) { fVertexSize = vertexsize; }
  61
  62   // z-channel map
  63   Int_t GenerateZChannelMap(); // could have different modes (for beam-beam, cosmics, ...)
  64   Bool_t DisplayZChannelMap(Int_t zchannel = -1, Int_t subch = 0) const;
  65
  66   // variables for pt-reconstruction (not used at the moment)
  67   Bool_t GenerateRecoCoefficients(Int_t trackletMask);
  68   Float_t GetAki(Int_t k, Int_t i);
  69   Float_t GetBki(Int_t k, Int_t i);
  70   Float_t GetCki(Int_t k, Int_t i);
  71 //  Float_t GetD(Int_t k) const;
  72
  73   // B-field
  74   void SetMagField(Float_t field) { fMagField = field; }
  75   Float_t GetMagField() const { return fMagField; }
  76
  77  protected:
  78   static const Int_t fgkNZChannels = 3; // No. of z-channels
  79   static const Int_t fgkNLinks = 12;    // No. of links
  80   static const Int_t fgkFixLayer = 2;   // which layer is fixed for the generation of the z-channel map
  81   static const Int_t fgkDeltaY = 19;    // accepted deviation in y_proj, default: 9
  82   static const Int_t fgkDeltaAlpha = 21; // accepted deviation in alpha, default: 11
  83   static const Int_t fgkNRefLayers = 3;  // no. of reference layers
  84
  85   static const Float_t fgkBinWidthY; // bin width for y-position
  86   static const Float_t fgkBinWidthdY; // bin width for deflection length
  87
  88   static const Int_t fgkBitWidthY; // bit width for y-position
  89   static const Int_t fgkBitWidthdY; // bit width for deflection length
  90   static const Int_t fgkBitWidthYProj; // bit width for projected y-position
  91   static const Int_t fgkBitExcessY; // excess bits for y-position
  92   static const Int_t fgkBitExcessAlpha; // excess bits for alpha
  93   static const Int_t fgkBitExcessYProj; // excess bits for projected y-position
  94
  95   Float_t fVertexSize;          // assumed vertex size (z-dir.) for the z-channel map
  96
  97   Int_t fZChannelMap[5][16][6][16];               // must be changed
  98   Int_t fZSubChannel[5][fgkNZChannels][6][16];    // must be changed
  99
 100   Int_t fCurrTrackletMask; // current tracklet mask for which the coefficients have been calculated
 101   Float_t fAki[6]; // coefficients used for the fit, calculated for the current tracklet mask
 102   Float_t fBki[6]; // coefficients used for the fit, calculated for the current tracklet mask
 103   Float_t fCki[6]; // coefficients used for the fit, calculated for the current tracklet mask
 104
 105   Int_t *fRefLayers;            //[fgkNRefLayers] reference layers for track finding
 106
 107   Float_t fMagField;            // magnetic field in T
 108
 109   AliTRDgeometry *fGeo;         //! pointer to the TRD geometry
 110
 111   static AliTRDgtuParam *fgInstance; // instance pointer
 112
 113  private:
 114   AliTRDgtuParam();                          // instance only via Instance()
 115   AliTRDgtuParam(const AliTRDgtuParam &rhs); // not implemented
 116   AliTRDgtuParam& operator=(const AliTRDgtuParam &rhs); // not implemented
 117
 118   ClassDef(AliTRDgtuParam, 1);
 119 };
 120
 121 #endif