-#ifndef TRDdigitizer_h
-#define TRDdigitizer_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>
+////////////////////////////////////////////////////////////////////////////
+// //
+// Produces digits from the hits information //
+// //
+////////////////////////////////////////////////////////////////////////////
-#include "AliHit.h"
-#include "AliDigit.h"
-#include "AliTRDconst.h"
-#include "AliTRDgeometry.h"
-#include "AliTRDsegmentArray.h"
+#include "AliDigitizer.h"
-///////////////////////////////////////////////////////
-// Produces digits from the hits information //
-///////////////////////////////////////////////////////
+class TFile;
+class TF1;
-const Int_t kNDict = 3;
+class AliDigitizationInput;
+class AliRunLoader;
-class AliTRDdigitizer : public TNamed {
+class AliTRD;
+class AliTRDdigitsManager;
+class AliTRDgeometry;
+class AliTRDarraySignal;
+class AliTRDarrayADC;
+class AliTRDmcmSim;
+
+class AliTRDdigitizer : public AliDigitizer {
public:
AliTRDdigitizer();
- AliTRDdigitizer(const Text_t* name, const Text_t* title);
- ~AliTRDdigitizer();
-
- virtual void Init();
- virtual Bool_t Open(const Char_t *name, Int_t nEvent = 0);
- virtual Bool_t MakeBranch();
- virtual Bool_t MakeDigits();
- 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 diff_on = 1) { fDiffusionOn = diff_on; };
- virtual void SetDiffusionT(Float_t diff) { fDiffusionT = diff; };
- virtual void SetDiffusionL(Float_t diff) { fDiffusionL = diff; };
- virtual void SetElAttach(Int_t el_on = 1) { fElAttachOn = el_on; };
- virtual void SetElAttachProp(Float_t prop) { fElAttachProp = prop; };
- virtual void SetExB(Int_t exb_on = 1) { fExBOn = exb_on; };
- virtual void SetLorentzAngle(Float_t angle) { fLorentzAngle = angle; };
-
- AliTRDsegmentArray *DigitsArray() { return fDigitsArray; };
- AliTRDsegmentArray *Dictionary(Int_t i) { return fDictionary[i]; };
-
- 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; };
+ AliTRDdigitizer(const Text_t *name, const Text_t *title);
+ AliTRDdigitizer(AliDigitizationInput* digInput, const Text_t *name, const Text_t *title);
+ AliTRDdigitizer(AliDigitizationInput* digInput);
+ AliTRDdigitizer(const AliTRDdigitizer &d);
+ virtual ~AliTRDdigitizer();
+ AliTRDdigitizer &operator=(const AliTRDdigitizer &d);
+
+ virtual void Copy(TObject &d) const;
+ Bool_t InitDetector();
+ void InitOutput(Int_t iEvent);
+ virtual void Digitize(const Option_t * option = 0);
+
+ virtual Bool_t Open(const Char_t *file, Int_t nEvent = 0);
+ virtual Bool_t Open(AliRunLoader * const runLoader, Int_t nEvent = 0);
+ virtual Bool_t MakeBranch(TTree *tree) const;
+ virtual Bool_t WriteDigits() const;
+
+ virtual void AddSDigitsManager(AliTRDdigitsManager *manager);
+ virtual void DeleteSDigitsManager();
+
+ virtual Bool_t MakeDigits();
+
+ Bool_t SortHits(Float_t **hits, Int_t *nhit);
+ Bool_t ConvertHits(Int_t det, const Float_t * const hits, Int_t nhit, AliTRDarraySignal *signals);
+ Bool_t ConvertSignals(Int_t det, AliTRDarraySignal *signals);
+
+ Bool_t Digits2SDigits(AliTRDdigitsManager * const manDig, AliTRDdigitsManager * const manSDig);
+ Bool_t SDigits2Digits();
+ Bool_t MergeSDigits();
+ Bool_t ConvertSDigits();
+
+ Bool_t Signal2ADC(Int_t det, AliTRDarraySignal *signals);
+ Bool_t Signal2SDigits(Int_t det, AliTRDarraySignal *signals);
+ Bool_t CopyDictionary(Int_t det);
+ void CompressOutputArrays(Int_t det);
+
+ void SetCompress(Int_t c = 1) { fCompress = c; }
+ void SetSDigits(Int_t v = 1) { fSDigits = v; }
+ void SetEvent(Int_t v = 0) { fEvent = v; }
+ void SetManager(AliTRDdigitsManager * const man) { fDigitsManager = man; }
+ void SetGeometry(AliTRDgeometry * const geo) { fGeo = geo; }
+ void SetMergeSignalOnly(Bool_t m = kTRUE) { fMergeSignalOnly = m; }
+
+ AliTRDdigitsManager *Digits() const { return fDigitsManager; }
+
+ Bool_t GetCompress() const { return fCompress; }
+ Bool_t GetSDigits() const { return fSDigits; }
+
+ Int_t Diffusion(Float_t vdrift, Double_t absdriftlength, Double_t exbvalue
+ , Double_t &lRow, Double_t &lCol, Double_t &lTime);
+ void RunDigitalProcessing(Int_t det = 0);
protected:
- TFile *fInputFile; //! ALIROOT-filename
-
- AliTRDsegmentArray *fDigitsArray; //! Array of detector segments containing the digits
- AliTRDsegmentArray *fDictionary[3]; //! Dictionary array, connecting MC tracks to the digits
-
- Int_t fEvent; //! Event number
-
- Float_t fGasGain; // Gas gain
- Float_t fNoise; // Electronics noise
- Float_t fChipGain; // Electronics gain
- Float_t fADCoutRange; // ADC output range (number of channels)
- Float_t fADCinRange; // ADC input range (input charge)
- Int_t fADCthreshold; // ADC threshold in ADC channel
- Int_t fDiffusionOn; // Switch for the diffusion
- Float_t fDiffusionT; // Diffusion in transverse direction
- Float_t fDiffusionL; // Diffusion in longitudinal 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 fLorentzFactor; // Factor due to Lorentz force
-
- private:
-
- virtual Int_t Diffusion(Float_t driftlength, Float_t *xyz);
- virtual Int_t ExB(Float_t driftlength, Float_t *xyz);
- virtual Float_t PadResponse(Float_t x);
-
- ClassDef(AliTRDdigitizer,1) // TRD-Digits manager
-
-};
-
-//_____________________________________________________________________________
-class AliTRDdigit : public AliDigitNew {
-
- public:
-
- AliTRDdigit() {};
- AliTRDdigit(Int_t *digits);
- virtual ~AliTRDdigit() {};
+ AliRunLoader *fRunLoader; //! Local pointer
+ AliTRDdigitsManager *fDigitsManager; //! Manager for the output digits
+ AliTRDdigitsManager *fSDigitsManager; //! Manager for the summed input s-digits
+ TList *fSDigitsManagerList; //! List of managers of input s-digits
+ AliTRD *fTRD; //! TRD detector class
+ AliTRDgeometry *fGeo; //! TRD geometry
- virtual Int_t GetAmp() { return fAmplitude; };
- virtual Int_t GetDetector() { return fDetector; };
- virtual Int_t GetRow() { return fRow; };
- virtual Int_t GetCol() { return fCol; };
- virtual Int_t GetTime() { return fTime; };
+ AliTRDmcmSim *fMcmSim; //! MCM simulation for digital processing
- protected:
+ Int_t fEvent; //! Event number
+ Int_t *fMasks; //! Masks for the merging
+ Bool_t fCompress; // Switch to keep only compressed data in memory
+ Bool_t fSDigits; // Switch for the summable digits
+ Bool_t fMergeSignalOnly; // Merge only detectors that contain a signal
- Int_t fDetector; // TRD detector number
- Int_t fRow; // Pad row number
- Int_t fCol; // Pad col number
- Int_t fTime; // Time bucket
- Int_t fAmplitude; // Signal amplitude
-
- ClassDef(AliTRDdigit,1) // Digits for Transition Radiation Detector
+ ClassDef(AliTRDdigitizer,20) // Produces TRD-Digits
};
-
#endif