Merge with TRD-develop
[u/mrichter/AliRoot.git] / TRD / AliTRD.h
index 6245efa..f691731 100644 (file)
-#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 "AliDetector.h"
-#include "AliHit.h" 
+
+class AliRun;
+class AliDigit;
+
+class AliTRDhit;
+class AliTRDsim;
+class AliTRDgeometry;
+
+//_____________________________________________________________________________
 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();
 
-////////////////////////////////////////////////
-//  Geometry parameter
-////////////////////////////////////////////////
+          AliTRD    &operator=(const AliTRD &trd);
 
-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)
-
-const Float_t rmin    = 294.0;   // r-dimensions of the TRD-frame
-const Float_t rmax    = 368.0;
-
-const Float_t zmax1   = 378.35;  // z-dimensions of the TRD-frame
-const Float_t zmax2   = 302.0;
-
-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
-
-const Float_t widpl1  = 99.6;    // rphi-dimension of the innermost plane
-
-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
+  virtual void       AddHit(Int_t track, Int_t det, Float_t *hits, Int_t q);
+  virtual void       AddDigit(Int_t *digits, Int_t *amp);    
+  virtual void       AddRecPoint(Float_t *pos, Int_t *digits
+                               , Int_t det, Float_t amp, Int_t *tracks);
+  virtual void       BuildGeometry();
+  virtual void       Copy(TObject &trd);
+  virtual void       CreateGeometry();
+  virtual void       CreateMaterials();
+  virtual void       DrawModule();
+  Int_t              DistancetoPrimitive(Int_t px, Int_t py);
+  virtual void       LoadPoints(Int_t track);    
+  TObjArray         *RecPoints() const              { return fRecPoints;   };
+  virtual void       Init();
+  virtual Int_t      IsVersion() const = 0;
+  virtual void       MakeBranch(Option_t* option);     
+  virtual void       ResetRecPoints();
+  virtual void       StepManager() = 0; 
+  virtual void       SetTreeAddress();
 
-////////////////////////////////////////////////
-// Parameter for the energy loss calculation
-////////////////////////////////////////////////
+  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;  };
 
-const Float_t kPoti = 12.1;      // First ionization potential (eV)
-const Float_t kEend = 50000.0;   // Maximum energy (50 keV);
+  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 Int_t      GetSensChamber() const     = 0;
+  virtual Int_t      GetSensPlane() const       = 0;
+  virtual Int_t      GetSensSector() const      = 0;
+  virtual Int_t      GetSensSectorRange() const = 0; 
+
+  virtual AliTRDsim *CreateTR()     = 0;
+  virtual AliTRDsim *GetTR() const  = 0;
+
+ protected:
+
+  Int_t              fGasMix;            //  Gas mixture. 0: Xe/Isobutane 1: Xe/CO2
+
+  AliTRDgeometry    *fGeometry;          //  The TRD geometry
+
+  TObjArray         *fRecPoints;         //  Array of reconstructed points
+  Int_t              fNRecPoints;        //! Number of reconstructed points
+
+  Float_t            fGasDensity;        //  The density of the drift gas
+  Float_t            fFoilDensity;       //  The density of the entrance window foil
+
+  Int_t              fDrawTR;            //  Switches marking the TR photons in the display
+  Int_t              fDisplayType;       //  Display type (0: normal, 1: detailed) 
+
+  ClassDef(AliTRD,1)                     //  Transition Radiation Detector base class
 
-//_____________________________________________________________________________
-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
 };
 
 #endif