]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDdigitizer.h
New TRF and PRF. Speedup of the code. Digits from amplification region included
[u/mrichter/AliRoot.git] / TRD / AliTRDdigitizer.h
index 59eb90cb4a95e2edbb0d075c028cd07486143288..53bee41d63d58df3b0ad4c963df071a8a2c20736 100644 (file)
@@ -6,15 +6,13 @@
 /* $Id$ */
 
 #include <TNamed.h>
-#include <TFile.h>
-#include <TF1.h>
 
-#include "AliHit.h" 
-#include "AliTRDdigit.h"
-#include "AliTRDconst.h"
-#include "AliTRDgeometry.h"
+class TFile;
+class TF1;
 
+class AliTRD;
 class AliTRDdigitsManager;
+class AliTRDgeometry;
 
 ///////////////////////////////////////////////////////
 //  Produces digits from the hits information        //
@@ -34,37 +32,56 @@ class AliTRDdigitizer : public TNamed {
   virtual void         Init();
   virtual Bool_t       Open(const Char_t *name, Int_t nEvent = 0);
   virtual Bool_t       MakeDigits();
+  virtual void         ReInit();
+  virtual Bool_t       SumSDigits();
   virtual Bool_t       WriteDigits();
-
-  virtual void         SetGasGain(Float_t gasgain)     { fGasGain       = gasgain;  };
-  virtual void         SetNoise(Float_t noise)         { fNoise         = noise;    };
-  virtual void         SetChipGain(Float_t chipgain)   { fChipGain      = chipgain; };
-  virtual void         SetADCoutRange(Float_t range)   { fADCoutRange   = range;    };
-  virtual void         SetADCinRange(Float_t range)    { fADCinRange    = range;    };
-  virtual void         SetADCthreshold(Int_t thresh)   { fADCthreshold  = thresh;   };
-  virtual void         SetDiffusion(Int_t diffOn = 1)  { fDiffusionOn   = diffOn;   };
-  virtual void         SetDiffusionT(Float_t diff)     { fDiffusionT    = diff;     };
-  virtual void         SetDiffusionL(Float_t diff)     { fDiffusionL    = diff;     };
-  virtual void         SetElAttach(Int_t elOn = 1)     { fElAttachOn    = elOn;     };
-  virtual void         SetElAttachProp(Float_t prop)   { fElAttachProp  = prop;     };
-  virtual void         SetExB(Int_t exbOn = 1)         { fExBOn         = exbOn;    };
-  virtual void         SetLorentzAngle(Float_t angle)  { fLorentzAngle  = angle;    };
-  virtual void         SetPadResponse(TF1 *PRF)        { if (fPRF) delete fPRF;
-                                                         fPRF           = PRF;      };
-
-  AliTRDdigitsManager *Digits()                        { return fDigits;        };
-
-  virtual Float_t      GetGasGain()                    { return fGasGain;       };
-  virtual Float_t      GetNoise()                      { return fNoise;         };
-  virtual Float_t      GetChipGain()                   { return fChipGain;      };
-  virtual Float_t      GetADCoutRange()                { return fADCoutRange;   };
-  virtual Float_t      GetADCinRange()                 { return fADCinRange;    };
-  virtual Int_t        GetADCthreshold()               { return fADCthreshold;  };
-  virtual Float_t      GetDiffusionT()                 { return fDiffusionT;    };
-  virtual Float_t      GetDiffusionL()                 { return fDiffusionL;    };
-  virtual Float_t      GetElAttachProp()               { return fElAttachProp;  };
-  virtual Float_t      GetLorentzAngle()               { return fLorentzAngle;  };
-  virtual TF1         *GetPadResponse()                { return fPRF;           };
+  virtual Bool_t       InitDetector();
+
+  virtual void         SetGasGain(Float_t gasgain)      { fGasGain        = gasgain;  };
+  virtual void         SetNoise(Float_t noise)          { fNoise          = noise;    };
+  virtual void         SetChipGain(Float_t chipgain)    { fChipGain       = chipgain; };
+  virtual void         SetADCoutRange(Float_t range)    { fADCoutRange    = range;    };
+  virtual void         SetADCinRange(Float_t range)     { fADCinRange     = range;    };
+  virtual void         SetADCthreshold(Int_t thresh)    { fADCthreshold   = thresh;   };
+  virtual void         SetDiffusion(Int_t diffOn = 1)   { fDiffusionOn    = diffOn;   };
+  virtual void         SetDiffusionT(Float_t diff)      { fDiffusionT     = diff;     };
+  virtual void         SetDiffusionL(Float_t diff)      { fDiffusionL     = diff;     };
+  virtual void         SetElAttach(Int_t elOn = 1)      { fElAttachOn     = elOn;     };
+  virtual void         SetElAttachProp(Float_t prop)    { fElAttachProp   = prop;     };
+  virtual void         SetExB(Int_t exbOn = 1)          { fExBOn          = exbOn;    };
+  virtual void         SetOmegaTau(Float_t ot)          { fOmegaTau       = ot;      
+                                                          ReInit();                   };
+  virtual void         SetPadResponse(Int_t prfOn = 1)  { fPRFOn          = prfOn;    };
+  virtual void         SetPRF(TF1 *prf);
+  virtual void         SetTimeResponse(Int_t trfOn = 1) { fTRFOn          = trfOn;   
+                                                          ReInit();                   };
+  virtual void         SetDriftVelocity(Float_t v)      { fDriftVelocity  = v;       
+                                                          ReInit();                   };
+  virtual void         SetPadCoupling(Float_t v)        { fPadCoupling    = v;        };
+  virtual void         SetTimeCoupling(Float_t v)       { fTimeCoupling   = v;        };
+  virtual void         SetCompress(Int_t c = 1)         { fCompress       = c;        };
+  virtual void         SetVerbose(Int_t v = 1)          { fVerbose        = v;        };
+  virtual void         SetSDigits(Int_t v = 1)          { fSDigits        = v;        };
+
+  AliTRDdigitsManager *Digits() const                   { return fDigits;             };
+
+  virtual Float_t      GetGasGain() const               { return fGasGain;            };
+  virtual Float_t      GetNoise() const                 { return fNoise;              };
+  virtual Float_t      GetChipGain() const              { return fChipGain;           };
+  virtual Float_t      GetADCoutRange() const           { return fADCoutRange;        };
+  virtual Float_t      GetADCinRange() const            { return fADCinRange;         };
+  virtual Int_t        GetADCthreshold() const          { return fADCthreshold;       };
+  virtual Float_t      GetDiffusionT() const            { return fDiffusionT;         };
+  virtual Float_t      GetDiffusionL() const            { return fDiffusionL;         };
+  virtual Float_t      GetElAttachProp() const          { return fElAttachProp;       };
+  virtual Float_t      GetOmegaTau() const              { return fOmegaTau;           };
+  virtual TF1         *GetPadResponse() const           { return fPRF;                };
+  virtual Float_t      GetDriftVelocity() const         { return fDriftVelocity;      };
+  virtual Float_t      GetPadCoupling() const           { return fPadCoupling;        };
+  virtual Float_t      GetTimeCoupling() const          { return fTimeCoupling;       };
+  virtual Bool_t       GetCompress() const              { return fCompress;           };
+  virtual Bool_t       GetSDigits() const               { return fSDigits;            };
+  virtual Float_t      GetTimeBinWidth() const          { return fTimeBinWidth;       };
 
  protected:
 
@@ -80,6 +97,8 @@ class AliTRDdigitizer : public TNamed {
   Float_t              fChipGain;        // Electronics gain
   Float_t              fADCoutRange;     // ADC output range (number of channels)
   Float_t              fADCinRange;      // ADC input range (input charge)
+  Float_t              fSinRange;        // Input range for summable digits 
+  Float_t              fSoutRange;       // Output range for summable digits 
   Int_t                fADCthreshold;    // ADC threshold in ADC channel
   Int_t                fDiffusionOn;     // Switch for the diffusion
   Float_t              fDiffusionT;      // Diffusion in transverse direction
@@ -87,16 +106,41 @@ class AliTRDdigitizer : public TNamed {
   Int_t                fElAttachOn;      // Switch for the electron attachment
   Float_t              fElAttachProp;    // Propability for electron attachment (for 1m)
   Int_t                fExBOn;           // Switch for the ExB effects
-  Float_t              fLorentzAngle;    // Lorentz angle 
+  Float_t              fOmegaTau;        // Tangens of the Lorentz angle 
   Float_t              fLorentzFactor;   // Factor due to Lorentz force
+  Int_t                fPRFOn;           // Switch for the pad response
   TF1                 *fPRF;             // Pad response function
+  Float_t             *fPRFsmp;          //!Sampled pad response
+  Int_t                fPRFbin;          // Number of bins for the PRF
+  Float_t              fPRFlo;           // Lower boundary of the PRF
+  Float_t              fPRFhi;           // Higher boundary of the PRF
+  Float_t              fPRFwid;          // Bin width of the sampled PRF
+  Int_t                fPRFpad;          // Distance to next pad in PRF
+  Int_t                fTRFOn;           // Switch for the time response
+  Float_t             *fTRFsmp;          //!Integrated time response
+  Int_t                fTRFbin;          // Number of bins for the TRF
+  Float_t              fTRFlo;           // Lower boundary of the TRF
+  Float_t              fTRFhi;           // Higher boundary of the TRF
+  Float_t              fTRFwid;          // Bin width of the integrated TRF
+  Float_t              fDriftVelocity;   // Drift velocity (cm / mus)
+  Float_t              fTimeBinWidth;    // Time bin width in ns
+  Float_t              fPadCoupling;     // Pad coupling factor
+  Float_t              fTimeCoupling;    // Time coupling factor (image charge of moving ions)
+  Bool_t               fCompress;        // Switch to keep only compressed data in memory
+  Int_t                fVerbose;         // Sets the verbose level
+  Bool_t               fSDigits;         // Switch for the summable digits
 
  private:
 
-  virtual Int_t       Diffusion(Float_t driftlength, Float_t *xyz);
-  virtual Int_t       ExB(Float_t driftlength, Float_t *xyz);  
-  
-  ClassDef(AliTRDdigitizer,1)            // Produces TRD-Digits
+  virtual Int_t        Diffusion(Float_t driftlength, Float_t *xyz);
+  virtual Int_t        ExB(Float_t driftlength, Float_t *xyz);  
+  virtual Int_t        PadResponse(Float_t signal, Float_t dist, Float_t *pad);
+  virtual Float_t      TimeResponse(Float_t time);  
+  virtual Bool_t       CheckDetector(Int_t plane, Int_t chamber, Int_t sector);
+  virtual void         SamplePRF();
+  virtual void         SampleTRF();
+
+  ClassDef(AliTRDdigitizer,4)            // Produces TRD-Digits
 
 };