GetMedium replaced by CurrentMedium (I.Hrivnacova)
[u/mrichter/AliRoot.git] / TOF / AliTOFGeometry.h
index 5c50996..51c12e8 100644 (file)
 ///////////////////////////////////////////////////////////////////////////////
 
 
+#include "TObject.h"
+#include "TGeoManager.h"
 
-class AliTOFGeometry{
+class AliTOFGeometry: public TObject{
 
  public:
   AliTOFGeometry();
   virtual ~AliTOFGeometry();
 
-  static  Int_t NStripA()     { return fgkNStripA;};
-  static  Int_t NStripB()     { return fgkNStripB;};
-  static  Int_t NStripC()     { return fgkNStripC;};
-  static  Int_t NpadX()       { return fgkNpadX;};
-  static  Int_t NpadZ()       { return fgkNpadZ;};
-  static  Int_t NSectors()    { return fgkNSectors;};
-  static  Int_t NPlates()     { return fgkNPlates;};
-  static  Int_t NPadXSector() { return (fgkNStripA + 2*fgkNStripB +
-                                      2*fgkNStripC)*fgkNpadX*fgkNpadZ;};
-  static  Int_t TimeDiff()    { return fgkTimeDiff;};
-  static  Int_t MaxTOFTree()  { return fgkMaxTOFTree;};
+  static  Int_t NStripA()     { return kNStripA;};
+  static  Int_t NStripB()     { return kNStripB;};
+  virtual Int_t NStripC()     { return kNStripC;};
+  static  Int_t NMaxNstrip()  { return kMaxNstrip;};
+  static  Int_t NpadX()       { return kNpadX;};
+  static  Int_t NpadZ()       { return kNpadZ;};
+  static  Int_t NpadXStrip()  { return kNpadX*kNpadZ;};
+  static  Int_t NSectors()    { return kNSectors;};
+  static  Int_t NPlates()     { return kNPlates;};
+  virtual Int_t NPadXSector() { return (kNStripA + 2*kNStripB +
+                                       2*kNStripC)*kNpadX*kNpadZ;};
+
+  virtual Float_t RinTOF()    { return fgkxTOF;};
+  virtual Float_t Rmin()      { return fgkRmin;};
+  virtual Float_t Rmax()      { return fgkRmax;};
+
+  static  Float_t XPad()      { return fgkXPad;};
+  static  Float_t ZPad()      { return fgkZPad;};
 
+  static  Float_t StripLength() { return fgkStripLength;};
+
+  static  Int_t TimeDiff()    { return fgkTimeDiff;};
+  static  Int_t MaxTOFTree()  { return kMaxTOFTree;};
 
-  static  Float_t Rmin()     { return fgkRmin;};
-  static  Float_t Rmax()     { return fgkRmax;};
-  static  Float_t ZlenA()    { return fgkZlenA;};
-  static  Float_t ZlenB()    { return fgkZlenB;};
-  static  Float_t ZlenC()    { return fgkZlenC;};
-  static  Float_t XPad()     { return fgkXPad;};
-  static  Float_t ZPad()     { return fgkZPad;};
-  static  Float_t MaxhZtof() { return fgkMaxhZtof;};
+  static  Int_t NDDL()        { return kNDDL;};
+  static  Int_t NTRM()        { return kNTRM;}
+  static  Int_t NTdc()        { return kNTdc;};
+  static  Int_t NCh()         { return kNCh;};
+  static  Int_t NPadXTRM()    { return kNCh*kNTdc;};
 
+  virtual  Float_t ZlenA()       { return kZlenA;};
+  virtual  Float_t ZlenB()       { return kZlenB;};
+  virtual  Float_t ZlenC()       { return kZlenC;};
+  virtual  Float_t MaxhZtof()    { return kMaxhZtof;};
 
   static  Float_t SigmaForTail1() { return fgkSigmaForTail1;};
   static  Float_t SigmaForTail2() { return fgkSigmaForTail2;};
-  static  Float_t SpeedOfLight()  { return fgkSpeedOfLight;};
-  static  Float_t PionMass()      { return fgkPionMass;};
-  static  Float_t KaonMass()      { return fgkKaonMass;};
-  static  Float_t ProtonMass()    { return fgkProtonMass;};
-  static  Float_t ElectronMass()  { return fgkElectronMass;};
-  static  Float_t MuonMass()      { return fgkMuonMass;};
  
+  static  Double_t GetAlpha()  { return 2 * 3.14159265358979323846 / kNSectors; }; 
  
+  static Float_t TdcBinWidth() {return fgkTdcBin;};
+
 
   virtual void    Init();
-  virtual void    GetPos(Int_t *det,Float_t *pos);
-  virtual void    GetDetID(Float_t *pos,Int_t *det);
-  virtual Int_t   GetPlate(Float_t *pos);
-  virtual Int_t   GetStrip(Float_t *pos);
-  virtual Int_t   GetSector(Float_t *pos);
-  virtual Int_t   GetPadX(Float_t *pos);
-  virtual Int_t   GetPadZ(Float_t *pos);
-  virtual Float_t GetX(Int_t *det);
-  virtual Float_t GetY(Int_t *det);
-  virtual Float_t GetZ(Int_t *det);
-  virtual Float_t GetMinPlateTheta(Int_t iPlate);
-  virtual Float_t GetMaxPlateTheta(Int_t iPlate);
-  virtual Float_t GetMinStripTheta(Int_t iPlate, Int_t iStrip);
-  virtual Float_t GetMaxStripTheta(Int_t iPlate, Int_t iStrip);
-  virtual Float_t GetStripTheta(Int_t iPlate, Int_t iStrip);
-  virtual Float_t GetAngles(Int_t iplate, Int_t istrip)  const {return fAngles[iplate][istrip];};
-  virtual Float_t GetHeights(Int_t iplate, Int_t istrip) const {return fHeights[iplate][istrip];};
-
-  private:
-
-  static const Int_t fgkNStripA    = 15; // number of strips in A type module 
-  static const Int_t fgkNStripB    = 19; // number of strips in B type module 
-  static const Int_t fgkNStripC    = 20; // number of strips in C type module 
-  static const Int_t fgkNpadX      = 48; // Number of pads along X 
-  static const Int_t fgkNpadZ      = 2;  // Number of pads along Z
-  static const Int_t fgkNSectors   = 18; // Number of Sectors
-  static const Int_t fgkNPlates    = 5;  // Number of Plates
-  static const Int_t fgkMaxNstrip  = 20; // Max. number of strips
-  static const Int_t fgkMaxTOFTree  = 5;    // numer of geom. levels: 
-  static const Int_t fgkTimeDiff   = 25000;// Min signal separation (ps)
-
-  static const Float_t fgkRmin     = 370.; // Inner radius of the TOF (cm)
-  static const Float_t fgkRmax     = 399;  // Outer radius of the TOF (cm)
-  static const Float_t fgkZlenA    = 106.0;// length (cm) of the A module
-  static const Float_t fgkZlenB    = 141.0;// length (cm) of the B module
-  static const Float_t fgkZlenC    = 177.5;// length (cm) of the C module
-  static const Float_t fgkXPad     = 2.5;  // Pad size in the x direction (cm)
-  static const Float_t fgkZPad     = 3.5;  // Pad size in the z direction (cm)
-  static const Float_t fgkMaxhZtof = 371.5;// Max half z-size of TOF (cm)
-
-
-  static const Float_t fgkSigmaForTail1= 2.;//Sig1 for simulation of TDC tails 
-  static const Float_t fgkSigmaForTail2= 0.5;//Sig2 for simulation of TDC tails
-  static const Float_t fgkSpeedOfLight = 0.299792458;// c (10^9 m/s)
-  static const Float_t fgkPionMass     = 0.13957;// pion mass (Gev/c^2)
-  static const Float_t fgkKaonMass     = 0.49368;// kaon mass (Gev/c^2)
-  static const Float_t fgkProtonMass   = 0.93827;// proton mass (Gev/c^2)
-  static const Float_t fgkElectronMass = 0.00051;// electron mass (Gev/c^2)
-  static const Float_t fgkMuonMass     = 0.10566;// muon mass (Gev/c^2)
-
-
-  static const Float_t fgkDprecMin = 0.0000075;//num.prec.tolerance on Thmin 
-  static const Float_t fgkDprecMax = 0.0000100;//num.prec.tolerance on Thma 
-  static const Float_t fgkDprecCen = 0.0000005;//num.prec.tolerance on <Theta> 
-  Float_t fAngles[fgkNPlates][fgkMaxNstrip]; //Strip Tilt Angles
-  Float_t fHeights[fgkNPlates][fgkMaxNstrip];//Strip heights
+  virtual void    ImportGeometry() {};
+  virtual void    SetHoles(Bool_t holes) {fHoles = holes;};
+  virtual Bool_t  GetHoles() const {return fHoles;};
+  virtual Bool_t  IsInsideThePadPar(Int_t */*det*/, Float_t */*pos*/) {return kFALSE;};
+  virtual Float_t DistanceToPadPar(Int_t */*det*/, Float_t */*pos*/, Float_t *dist3d=0) {return dist3d[0];};
+  virtual Bool_t  IsInsideThePad(Int_t */*det*/,TGeoHMatrix /*mat*/, Float_t */*pos*/){return kFALSE;};
+  virtual Float_t DistanceToPad(Int_t */*det*/,TGeoHMatrix /*mat*/, Float_t */*pos*/, Float_t *dist3d=0){return dist3d[0];};
+  virtual void    GetVolumePath(Int_t */*ind*/, Char_t */*path*/ ){};
+  virtual void    GetPos(Int_t */*det*/,Float_t */*pos*/){};
+  virtual void    GetPosPar(Int_t */*det*/,Float_t */*pos*/);
+  virtual void    GetDetID(Float_t */*pos*/,Int_t */*det*/);
+  virtual Int_t   GetPlate(Float_t */*pos*/) {return -1;};
+  virtual Int_t   GetStrip(Float_t */*pos*/) {return -1;};
+  virtual Int_t   GetSector(Float_t */*pos*/) {return -1;};
+  virtual Int_t   GetPadX(Float_t */*pos*/) {return -1;};
+  virtual Int_t   GetPadZ(Float_t */*pos*/) {return -1;};
+  virtual Float_t GetX(Int_t */*det*/) {return -500.;};
+  virtual Float_t GetY(Int_t */*det*/) {return -500.;};
+  virtual Float_t GetZ(Int_t */*det*/) {return -500.;};
+
+  Float_t GetAngles(Int_t iplate, Int_t istrip)  const {return fAngles[iplate][istrip];};
+  Float_t GetHeights(Int_t iplate, Int_t istrip) const {return fHeights[iplate][istrip];};
+  Float_t GetDistances(Int_t iplate, Int_t istrip) const {return fDistances[iplate][istrip];};
+
+  //private:
+  protected:
+
+  enum {
+    kNStripA    = 15, // number of strips in A type module 
+    kNStripB    = 19, // number of strips in B type module 
+    kNpadX      = 48, // Number of pads along X 
+    kNpadZ      = 2,  // Number of pads along Z
+    kNSectors   = 18, // Number of Sectors
+    kNPlates    = 5,  // Number of Plates
+    kMaxTOFTree = 5,  // numer of geom. levels: 
+    kMaxNstrip  = 20  // Max. number of strips
+  };
+
+  // DAQ characteristics
+  // cfr. TOF-TDR pag. 105 for Glossary
+  // TARODA : TOF-ALICE Read Out and Data Acquisition system
+  enum {
+    kNDDL        =    4, // Number of DDL (Detector Data Link) per sector
+    kNTRM        =   10, // Number of TRM ( Readout Module) per DDL
+    kNTdc        =   30, // Number of Tdc (Time to Digital Converter) per TRM
+    kNCh         =    8  // Number of channels per Tdc
+  };
+
+  static const Int_t fgkTimeDiff;      // Min signal separation (ps)
+
+  mutable Int_t kNStripC;       // number of strips in C type module 
+
+  mutable Float_t kZlenA;       // length (cm) of the A module
+  mutable Float_t kZlenB;       // length (cm) of the B module
+  mutable Float_t kZlenC;       // length (cm) of the C module
+  mutable Float_t kMaxhZtof;    // Max half z-size of TOF (cm)
+
+  mutable Float_t fgkRmin;     // Inner radius of the TOF (cm)
+  mutable Float_t fgkRmax;     // Outer radius of the TOF (cm)
+  mutable Float_t fgkxTOF;     // Inner TOF Radius used in Reconstruction (cm)
+
+  static const Float_t fgkStripLength; // Strip Length (rho X phi direction) (cm)
+
+  static const Float_t fgkXPad;     // Pad size in the x direction (cm)
+  static const Float_t fgkZPad;     // Pad size in the z direction (cm)
+
+  static const Float_t fgkSigmaForTail1;//Sig1 for simulation of TDC tails 
+  static const Float_t fgkSigmaForTail2;//Sig2 for simulation of TDC tails
+
+  Bool_t fHoles; //logical for geometry version (w/wo holes)
+
+  Float_t fAngles[kNPlates][kMaxNstrip];   //Strip Tilt Angles
+  Float_t fHeights[kNPlates][kMaxNstrip];  //Strip heights
+  Float_t fDistances[kNPlates][kMaxNstrip];//Strip distances
+
   Float_t fPhiSec; //sector Phi width (deg)
 
-  ClassDef(AliTOFGeometry,0) // TOF Geometry base class
+  static const Float_t fgkTdcBin;   // time-window for the TDC bins [ps]
+
+  ClassDef(AliTOFGeometry,3) // TOF Geometry base class
 };
 
 #endif