-#ifndef TRD_H
-#define TRD_H
+#ifndef ALITRD_H
+#define ALITRD_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice */
+
+/* $Id$ */
+
////////////////////////////////////////////////
// Manager and hits classes for set: TRD //
////////////////////////////////////////////////
-
+
+#include <TLorentzVector.h>
#include "AliDetector.h"
-#include "AliHit.h"
-
+#include <TVirtualMC.h>
+
+class TFile;
+class TLorentzVector;
+
+class AliRun;
+class AliDigit;
+
+class AliTRDhit;
+class AliTRDsim;
+class AliTRDgeometry;
+class AliTRDtrackHits;
+
+//_____________________________________________________________________________
class AliTRD : public AliDetector {
-
-protected:
- Int_t fGasMix; // Gas mixture. 0: Xe/Isobutane 1: Xe/CO2
- Int_t fSensSelect; // Switch to select only parts of the detector
- Int_t fSensPlane; // Sensitive detector plane
- Int_t fSensChamber; // Sensitive detector chamber
- Int_t fSensSector; // Sensitive detector sector
+ public:
-public:
AliTRD();
AliTRD(const char *name, const char *title);
- virtual ~AliTRD() {}
- virtual void AddHit(Int_t, Int_t*, Float_t*);
- virtual void BuildGeometry();
- virtual void CreateMaterials();
- Int_t DistancetoPrimitive(Int_t px, Int_t py);
- virtual void Init();
- virtual Int_t IsVersion() const =0;
- virtual void StepManager()=0;
- virtual void SetGasMix(Int_t imix = 0);
- virtual void SetSensPlane(Int_t iplane = 0);
- virtual void SetSensChamber(Int_t ichamber = 0);
- virtual void SetSensSector(Int_t isector = 0);
-
- ClassDef(AliTRD,1) // Transition Radiation Detector base class
-};
+ AliTRD(const AliTRD &trd);
+ virtual ~AliTRD();
+
+ AliTRD &operator=(const AliTRD &trd);
+
+ virtual void AddHit(Int_t track, Int_t det, Float_t *hits, Int_t q, Bool_t inDrift);
+ virtual void BuildGeometry();
+ virtual void Copy(TObject &trd);
+ virtual void CreateGeometry();
+ virtual void CreateMaterials();
+ virtual void DrawModule() const;
+ Int_t DistancetoPrimitive(Int_t px, Int_t py) const;
+ virtual void LoadPoints(Int_t track);
+ virtual void Init();
+ virtual Int_t IsVersion() const = 0;
+ virtual void MakeBranch(Option_t* option);
+ virtual void ResetDigits();
+ virtual void StepManager() = 0;
+ virtual void SetTreeAddress();
+
+ virtual void SetGasMix(Int_t imix = 0);
+ virtual void SetHits() {};
+ virtual void SetPHOShole();
+ virtual void SetRICHhole();
+ virtual void SetDrawTR(Int_t idraw = 1) { fDrawTR = idraw; };
+ virtual void SetDisplayType(Int_t type = 0) { fDisplayType = type; };
+
+ AliTRDgeometry *GetGeometry() const { return fGeometry; };
+
+ virtual void SetSensChamber(Int_t ichamber) = 0;
+ virtual void SetSensPlane(Int_t iplane) = 0;
+ virtual void SetSensSector(Int_t isector) = 0;
+ virtual void SetSensSector(Int_t isector, Int_t nsector) = 0;
+
+ virtual void SetHitTypeStandard() { fHitType = 1; };
+ virtual void SetHitTypeCompressed() { fHitType = 2; };
+ virtual void SetHitTypeBoth() { fHitType = 3; };
+
+ virtual Int_t GetSensChamber() const = 0;
+ virtual Int_t GetSensPlane() const = 0;
+ virtual Int_t GetSensSector() const = 0;
+ virtual Int_t GetSensSectorRange() const = 0;
+
+ virtual void Hits2Digits();
+ virtual void Hits2SDigits();
+ virtual AliDigitizer* CreateDigitizer(AliRunDigitizer* manager) const;
+ virtual void SDigits2Digits();
+ virtual void Digits2Raw();
-////////////////////////////////////////////////
-// Geometry parameter
-////////////////////////////////////////////////
+ virtual AliTRDsim *CreateTR() = 0;
+ virtual AliTRDsim *GetTR() const = 0;
-const Int_t nsect = 18; // Number of sectors in the full detector
-const Int_t nmodul = 6; // Number of modules in each sector
-const Int_t ncham = 6; // Number of different chambers
-const Int_t narmsec = 5; // Number of sectors in one arm (geometry 1)
+ virtual void MakeBranch2(Option_t *opt=" ", const char *file=0);
+ virtual void SetTreeAddress2();
+ virtual void AddHit2(Int_t track, Int_t det, Float_t *hits, Int_t q, Bool_t inDrift);
+ virtual void ResetHits();
-const Float_t rmin = 294.0; // r-dimensions of the TRD-frame
-const Float_t rmax = 368.0;
+ virtual AliHit* FirstHit(Int_t track);
+ virtual AliHit* NextHit();
+ virtual AliHit* FirstHit2(Int_t track);
+ virtual AliHit* NextHit2();
-const Float_t zmax1 = 378.35; // z-dimensions of the TRD-frame
-const Float_t zmax2 = 302.0;
+ virtual void FinishPrimary();
+ virtual void RemapTrackHitIDs(Int_t *map);
-const Float_t zleni = 110.0; // z-dimension of the inner chamber
-const Float_t zlenn = 156.0; // z-dimension of the neighbouring chambers
-const Float_t zleno = 156.0; // z-dimension of the outer chambers
+ protected:
-const Float_t widpl1 = 99.6; // rphi-dimension of the innermost plane
+ Int_t fGasMix; // Gas mixture. 0: Xe/Isobutane 1: Xe/CO2
-const Float_t alframe = 1.0; // Thickness of the aluminium of the support frame
-
-const Float_t ccframe = 1.0; // Thickness of the carbon chamber frame
-
-// Thicknesses of the the material layers
-const Float_t sethick = 0.02; // Radiator seal
-const Float_t rathick = 4.2; // Radiator
-const Float_t pethick = 0.20; // PE-layer in the radiator
-const Float_t mythick = 0.005; // Mylar-layer
-const Float_t xethick = 3.5; // Gas mixture
-const Float_t cuthick = 0.001; // Pad plane
-const Float_t suthick = 0.06; // HEXCEL+G10 support structure (= 0.31% X0)
-const Float_t fethick = 0.0044; // FEE + signal lines (= 0.31% X0)
-const Float_t cothick = 0.02; // PE of the cooling device
-const Float_t wathick = 0.01; // Cooling water
-
-// Position of the material layers
-const Float_t sezpos = -5.657; // Radiator seal
-const Float_t razpos = -3.557; // Radiator
-const Float_t pezpos = 0.0; // PE-layer in the radiator
-const Float_t myzpos = -1.4545; // Mylar-layer
-const Float_t xezpos = 0.298; // Gas mixture
-const Float_t cuzpos = 2.047; // Pad plane
-const Float_t suzpos = 3.046; // HEXCEL+G10 support structure
-const Float_t fezpos = 4.0482; // FEE + signal lines
-const Float_t cozpos = 4.1504; // PE of the cooling devices
-const Float_t wazpos = 4.3004; // Cooling water
+ AliTRDgeometry *fGeometry; // The TRD geometry
-////////////////////////////////////////////////
-// Parameter for the energy loss calculation
-////////////////////////////////////////////////
+ Float_t fGasDensity; // The density of the drift gas
+ Float_t fFoilDensity; // The density of the entrance window foil
-const Float_t kPoti = 12.1; // First ionization potential (eV)
-const Float_t kEend = 50000.0; // Maximum energy (50 keV);
+ Int_t fDrawTR; // Switches marking the TR photons in the display
+ Int_t fDisplayType; // Display type (0: normal, 1: detailed)
-//_____________________________________________________________________________
-class AliTRDhit : public AliHit {
-public:
- Int_t fSector; // TRD sector number
- Int_t fChamber; // TRD chamber number
- Int_t fPlane; // TRD plane number
- Float_t fQ; // Charge created by a hit (geometry 2)
-
-public:
- AliTRDhit() {}
- AliTRDhit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits);
- virtual ~AliTRDhit() {}
-
- ClassDef(AliTRDhit,1) // Hits for Transition Radiation Detector
+ AliTRDtrackHits *fTrackHits; //! Compressed hits for a given track
+ Int_t fHitType; // Switch for the different hit data structures
+
+ ClassDef(AliTRD,6) // Transition Radiation Detector base class
+
};
#endif