#ifndef ALITOFGEOMETRY_H
#define ALITOFGEOMETRY_H
+
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
// //
///////////////////////////////////////////////////////////////////////////////
+#include "TGeoMatrix.h"
+#include "TObject.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() const { return fNStripC;};
+ 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 NStripXSector() const { return (kNStripA + 2*kNStripB +
+ 2*fNStripC);};
+ virtual Int_t NPadXSector() const { return (kNStripA + 2*kNStripB +
+ 2*fNStripC)*kNpadX*kNpadZ;};
+
+ virtual Float_t RinTOF() const { return fxTOF;};
+ virtual Float_t Rmin() const { return fRmin;};
+ virtual Float_t Rmax() const { return fRmax;};
- 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 Float_t StripLength() { return fgkStripLength;};
+
+ static Int_t TimeDiff() { return fgkTimeDiff;};
+ static Int_t MaxTOFTree() { return kMaxTOFTree;};
+
+ static Int_t NDDL() { return kNDDL;};
+ static Int_t NTRM() { return kNTRM;}
+ static Int_t NTdc() { return kNTdc;};
+ static Int_t NChain() { return kNChain;};
+ static Int_t NCh() { return kNCh;};
+ static Int_t NPadXTRM() { return kNCh*kNTdc*kNChain;};
+
+ virtual Float_t ZlenA() const { return fZlenA;};
+ virtual Float_t ZlenB() const { return fZlenB;};
+ virtual Float_t ZlenC() const { return fZlenC;};
+ virtual Float_t MaxhZtof() const { return fMaxhZtof;};
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;};
+ static Float_t ToTBinWidth() {return fgkToTBin;};
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*/) const {return kFALSE;};
+ virtual Float_t DistanceToPadPar(Int_t */*det*/, Float_t */*pos*/, Float_t *dist3d=0) const {return dist3d[0];};
+ virtual Bool_t IsInsideThePad(Int_t */*det*/,TGeoHMatrix /*mat*/, Float_t */*pos*/) const {return kFALSE;};
+ virtual Float_t DistanceToPad(Int_t */*det*/,TGeoHMatrix /*mat*/, Float_t */*pos*/, Float_t *dist3d=0) const {return dist3d[0];};
+ virtual void GetVolumePath(Int_t */*ind*/, Char_t */*path*/ ){};
+ virtual void GetVolumePath(Int_t /*sector*/, Char_t */*path*/ ){};
+ virtual void GetVolumePath(Int_t /*sector*/, Int_t /*plate*/, Int_t /*strip*/, Char_t */*path*/ ){};
+ virtual void GetPos(Int_t */*det*/,Float_t */*pos*/){};
+ virtual void GetPosPar(Int_t *det,Float_t *pos) const;
+ virtual void GetDetID(Float_t *pos,Int_t *det) const;
+ virtual Int_t GetPlate(Float_t */*pos*/) const {return -1;};
+ virtual Int_t GetStrip(Float_t */*pos*/) const {return -1;};
+ virtual Int_t GetSector(Float_t */*pos*/) const {return -1;};
+ virtual Int_t GetPadX(Float_t */*pos*/) const {return -1;};
+ virtual Int_t GetPadZ(Float_t */*pos*/) const {return -1;};
+ virtual Float_t GetX(Int_t */*det*/) const {return -500.;};
+ virtual Float_t GetY(Int_t */*det*/) const {return -500.;};
+ virtual Float_t GetZ(Int_t */*det*/) const {return -500.;};
+ virtual void DetToStripRF(Int_t nPadX, Int_t nPadZ,
+ Float_t &x, Float_t &z) const;
+ virtual void DetToSectorRF(Int_t /*vol*/[5], Double_t ** /*coord*/) { };
+
+ 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 = 12, // Number of TRM ( Readout Module) per DDL
+ kNTdc = 15, // Number of Tdc (Time to Digital Converter) per TRM
+ kNChain = 2, // Number of chains per TRM
+ kNCh = 8 // Number of channels per Tdc
+ };
+
+ static const Int_t fgkTimeDiff; // Min signal separation (ps)
+
+ mutable Int_t fNStripC; // number of strips in C type module
+
+ mutable Float_t fZlenA; // length (cm) of the A module
+ mutable Float_t fZlenB; // length (cm) of the B module
+ mutable Float_t fZlenC; // length (cm) of the C module
+ mutable Float_t fMaxhZtof; // Max half z-size of TOF (cm)
+
+ mutable Float_t fRmin; // Inner radius of the TOF (cm)
+ mutable Float_t fRmax; // Outer radius of the TOF (cm)
+ mutable Float_t fxTOF; // 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-of-flight bin width [ps]
+ static const Float_t fgkToTBin; // time-over-threshold bin width [ps]
+
+ ClassDef(AliTOFGeometry,4) // TOF Geometry base class
};
#endif