-#ifndef ALITRDDIGITIZER_h
-#define ALITRDDIGITIZER_h
+#ifndef ALITRDDIGITIZER_H
+#define ALITRDDIGITIZER_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
/* $Id$ */
#include <TNamed.h>
-#include <TFile.h>
-#include <TF1.h>
-#include "AliHit.h"
-#include "AliTRDdigit.h"
-#include "AliTRDconst.h"
-#include "AliTRDgeometry.h"
+// Time response function of the preamp
+Double_t TRFlandau(Double_t *x, Double_t *par);
+class TFile;
+class TF1;
+
+class AliTRD;
class AliTRDdigitsManager;
+class AliTRDgeometry;
///////////////////////////////////////////////////////
// Produces digits from the hits information //
AliTRDdigitizer();
AliTRDdigitizer(const Text_t* name, const Text_t* title);
- AliTRDdigitizer(AliTRDdigitizer &d);
+ AliTRDdigitizer(const AliTRDdigitizer &d);
virtual ~AliTRDdigitizer();
+ AliTRDdigitizer &operator=(const AliTRDdigitizer &d);
- virtual void Copy(AliTRDdigitizer &d);
+ virtual void Copy(TObject &d);
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; };
-
- inline AliTRDdigitizer &operator=(AliTRDdigitizer &d);
+ 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 SetTRF(TF1 *trf);
+ 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 TF1 *GetTimeResponse() const { return fTRF; };
+ 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:
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
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
+ TF1 *fTRF; // Time response function of the shaper
+ Float_t *fTRFint; //!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,3) // Produces TRD-Digits
-//_____________________________________________________________________________
-AliTRDdigitizer &AliTRDdigitizer::operator=(AliTRDdigitizer &d)
-{
- //
- // Assignment operator
- //
-
- if (this != &d) d.Copy(*this);
- return *this;
-
-}
+};
#endif