X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;ds=sidebyside;f=TRD%2FAliTRDdigitizer.h;h=304fa8b677f6431c1eb8ae7e91ee605a2bf164c6;hb=fb51265cefc3db171991b44e50d95058cd433a4d;hp=21307d4440e84c574a9908c3291b20010c5ffdb5;hpb=cf31434d07aea087e173751d046261682e01d691;p=u%2Fmrichter%2FAliRoot.git diff --git a/TRD/AliTRDdigitizer.h b/TRD/AliTRDdigitizer.h index 21307d4440e..304fa8b677f 100644 --- a/TRD/AliTRDdigitizer.h +++ b/TRD/AliTRDdigitizer.h @@ -5,174 +5,102 @@ /* $Id$ */ +//////////////////////////////////////////////////////////////////////////// +// // +// Produces digits from the hits information // +// // +//////////////////////////////////////////////////////////////////////////// + #include "AliDigitizer.h" class TFile; class TF1; -class AliRunDigitizer; +class AliDigitizationInput; +class AliRunLoader; class AliTRD; class AliTRDdigitsManager; class AliTRDgeometry; - -/////////////////////////////////////////////////////// -// Produces digits from the hits information // -/////////////////////////////////////////////////////// +class AliTRDarraySignal; +class AliTRDarrayADC; +class AliTRDmcmSim; class AliTRDdigitizer : public AliDigitizer { public: AliTRDdigitizer(); - AliTRDdigitizer(const Text_t* name, const Text_t* title); - AliTRDdigitizer(AliRunDigitizer *manager, const Text_t* name, const Text_t* title); - AliTRDdigitizer(AliRunDigitizer *manager); + 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); - virtual Bool_t Init(); - virtual Bool_t InitDetector(); - virtual Bool_t ReInit(); - virtual void Exec(Option_t* option = 0); + 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 MakeBranch(const Char_t *file = 0); - virtual Bool_t MakeDigits(); + 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 ConvertSDigits(); - virtual Bool_t MergeSDigits(); - virtual Bool_t SDigits2Digits(); - 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 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 SetPadResponse(Int_t prfOn = 1) { fPRFOn = prfOn; }; - virtual void SetTimeResponse(Int_t trfOn = 1) { fTRFOn = trfOn; - ReInit(); }; - virtual void SetCrossTalk(Int_t ctOn = 1) { fCTOn = ctOn; - ReInit(); }; - virtual void SetTailCancelation(Int_t tcOn = 1) { fTCOn = tcOn; }; - virtual void SetNexponential(Int_t nexp) { fTCnexp = nexp; }; - 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 SetTiltingAngle(Float_t v); - virtual void SetCompress(Int_t c = 1) { fCompress = c; }; - virtual void SetDebug(Int_t v = 1) { fDebug = v; }; - virtual void SetSDigits(Int_t v = 1) { fSDigits = v; }; - virtual void SetSDigitsScale(Float_t s) { fSDigitsScale = s; }; - virtual void SetEvent(Int_t v = 0) { fEvent = v; }; - virtual void SetManager(AliTRDdigitsManager *man) { fDigitsManager = man; }; - virtual void SetHitTypeStandard() { if (fTRD) fTRD->SetHitTypeStandard(); }; - virtual void SetHitTypeCompressed() { if (fTRD) fTRD->SetHitTypeCompressed(); }; - - AliTRDdigitsManager *Digits() const { return fDigitsManager; }; - - Float_t GetGasGain() const { return fGasGain; }; - Float_t GetNoise() const { return fNoise; }; - Float_t GetChipGain() const { return fChipGain; }; - Float_t GetADCoutRange() const { return fADCoutRange; }; - Float_t GetADCinRange() const { return fADCinRange; }; - Int_t GetADCthreshold() const { return fADCthreshold; }; - Float_t GetDiffusionT() const { return fDiffusionT; }; - Float_t GetDiffusionL() const { return fDiffusionL; }; - Float_t GetElAttachProp() const { return fElAttachProp; }; - Int_t GetExB() const { return fExBOn; }; - Float_t GetOmegaTau() const { return fOmegaTau; }; - Float_t GetDriftVelocity() const { return fDriftVelocity; }; - Float_t GetPadCoupling() const { return fPadCoupling; }; - Float_t GetTimeCoupling() const { return fTimeCoupling; }; - Bool_t GetCompress() const { return fCompress; }; - Bool_t GetSDigits() const { return fSDigits; }; - Float_t GetSDigitsScale() const { return fSDigitsScale; }; - Float_t GetTimeBinWidth() const { return fTimeBinWidth; }; - Float_t GetTiltingAngle() const; - virtual Float_t GetDiffusionL(Float_t vd, Float_t b); - virtual Float_t GetDiffusionT(Float_t vd, Float_t b); - virtual Float_t GetOmegaTau(Float_t vd, Float_t b); + + 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-file + 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 - - Int_t fEvent; //! Event number - - Int_t *fMasks; //! Masks for the merging - - Float_t fField; // Magnetic field - 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 fOmegaTau; // Tangens of the Lorentz angle - Float_t fLorentzFactor; // Factor due to Lorentz force - Int_t fPRFOn; // Switch for the pad response - 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 - Int_t fCTOn; // Switch for cross talk - Float_t *fCTsmp; //! Integrated cross talk - Int_t fTCOn; // Switch for the tail cancelation - Int_t fTCnexp; // Number of exponential of the digital filter - 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) - Float_t fTiltingAngle; // Tilting angle of the readout pads - Bool_t fCompress; // Switch to keep only compressed data in memory - Int_t fDebug; // Sets the debug level - Bool_t fSDigits; // Switch for the summable digits - Float_t fSDigitsScale; // Scale factor for the summable digits - - private: - - virtual Float_t Col0Tilted(Float_t col0, Float_t rowOffset, Int_t plane); - virtual Float_t CrossTalk(Float_t time); - 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, Int_t plane, Float_t *pad); - virtual Float_t TimeResponse(Float_t time); - virtual void DeConvExp(Double_t *source, Double_t *target, Int_t n, Int_t nexp); - virtual Bool_t CheckDetector(Int_t plane, Int_t chamber, Int_t sector); - virtual void SamplePRF(); - virtual void SampleTRF(); - - ClassDef(AliTRDdigitizer,6) // Produces TRD-Digits -}; + AliTRDmcmSim *fMcmSim; //! MCM simulation for digital processing + + 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 + ClassDef(AliTRDdigitizer,20) // Produces TRD-Digits + +}; #endif