Fix bug in tracklet reconstruction and add option to AliTRDfeeParam
[u/mrichter/AliRoot.git] / TRD / AliTRDmcmSim.h
1 #ifndef ALITRDMCMSIM_H
2 #define ALITRDMCMSIM_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice                               */
5
6 /* $Id$ */
7
8 ///////////////////////////////////////////////////////
9 //                                                   //
10 //  Multi Chip Module Simulation Class               //
11 //                                                   //
12 ///////////////////////////////////////////////////////
13
14 #include <TObject.h>
15
16 class AliTRDfeeParam;
17 class AliTRDSimParam;
18 class AliTRDcalibDB;
19 class AliTRDgeometry;
20 class AliTRDtrapAlu;
21 class AliTRDpadPlane;
22
23 class AliTRDmcmSim : public TObject {
24
25  public:
26
27                     AliTRDmcmSim();
28                     AliTRDmcmSim(const AliTRDmcmSim &m);
29   virtual          ~AliTRDmcmSim();
30   AliTRDmcmSim      &operator=(const AliTRDmcmSim &m);
31
32   virtual void      Copy(TObject &m) const;
33
34           void      Init( Int_t cha, Int_t rob, Int_t mcm, Bool_t newEvent );   // Initialize MCM by the position parameters
35   //void      Init( Int_t cha, Int_t rob, Int_t mcm );   // Initialize MCM by the position parameters
36           void      SetData(Int_t iadc, Int_t *adc);           // Set ADC data with array 
37           void      SetData(Int_t iadc, Int_t it, Int_t adc ); // Set ADC data
38           void      SetDataPedestal(Int_t iadc );              // Fill ADC data with pedestal values
39
40           Int_t     GetChaId() const  { return fChaId;  };     // Returns Chamber ID (0-539)
41           Int_t     GetRobPos() const { return fRobPos; };     // Returns ROB position (0-7)
42           Int_t     GetMcmPos() const { return fMcmPos; };     // Returns MCM position (0-17) (16,17 are mergers)
43           Int_t     GetRow() const    { return fRow;    };     // Returns Row number on chamber where the MCM is sitting
44           Int_t     GetCol( Int_t iadc );                      // Get corresponding column (0-143) from for ADC channel iadc = [0:20]
45           // for the ADC/Col mapping, see: http://wiki.kip.uni-heidelberg.de/ti/TRD/index.php/Image:ROB_MCM_numbering.pdf
46
47           Int_t*    GetPosLUT();
48
49           Int_t     ProduceRawStream( UInt_t *buf, Int_t bufsize ); // Produce raw data stream from this MCM
50           Int_t     ProduceTrackletStream( UInt_t *buf, Int_t bufsize ); // produce the tracklet stream for this MCM
51           void      Filter();                                  // Apply digital filters for existing data
52           void      ZSMapping();                               // Do ZS mapping for existing data
53           void      DumpData( char *f, char *target );         // Dump data stored (only for debugging)
54           void      Tracklet();
55           void      SetPosLUT();
56
57  protected:
58
59           Bool_t    fInitialized;                       // Status whether the class is initialized or not
60           Int_t    fNextEvent;                         // 0, if new event; 1 if not
61           Int_t     fMaxTracklets;                      // maximum number of tracklet-words submitted per mcm **new**
62           Int_t     fChaId;                             // Chamber ID
63           Int_t     fSector;                            // Sector number of the supermodule
64           Int_t     fStack;                             // Chamber stack ID
65           Int_t     fLayer;                             // Chamber layer ID
66           Int_t     fRobPos;                            // ROB Position on chamber
67           Int_t     fMcmPos;                            // MCM Position on chamber
68           Int_t     fNADC;                              // Number of ADC (usually 21)
69           Int_t     fNTimeBin;                          // Number of Timebins (variable)
70           Int_t     fRow;                               // Pad row number (0-11 or 0-15) of the MCM on chamber
71           Int_t   **fADCR;                              // Array with MCM ADC values (Raw)
72           Int_t   **fADCF;                              // Array with MCM ADC values (Filtered)
73           Int_t   **fADCT;                              // Array with MCM ADC values (filtered) for tracklet (12bits) //***NEW***
74           Int_t    *fPosLUT;                            // position lookup table **new**
75           UInt_t   *fMCMT;                              // tracklet word for one mcm/trap-chip **new**
76           Int_t   **fZSM;                               // Zero suppression map
77           Int_t    *fZSM1Dim;                           // Zero suppression map (1 dimensional projection)
78
79         
80         
81
82
83           // Parameter classes
84           AliTRDfeeParam *fFeeParam;                    // FEE parameters
85           AliTRDSimParam *fSimParam;                    // Simulation parameters
86           AliTRDcalibDB  *fCal;                         // Calibration interface
87           AliTRDgeometry *fGeo;                         // Geometry
88
89           Bool_t    CheckInitialized();                 // Check whether the class is initialized
90           void      FilterPedestal();                   // Apply pedestal filter
91           void      FilterGain();                       // Apply gain filter
92           void      FilterTail();                       // Apply tail filter
93
94           // Here are the several sub functions used only in FilterTail
95           void      FilterSimDeConvExpA(Int_t *source, Double_t *target, Int_t n, Int_t nexp);
96           void      FilterSimDeConvExpD(Int_t *source, Int_t *target, Int_t n, Int_t nexp);
97           void      FilterSimDeConvExpMI(Int_t *source, Double_t *target, Int_t n);
98           void      FilterSimTailMakerSpline(Double_t *ampin, Double_t *ampout, Double_t lambda, Int_t n);
99           void      FilterSimTailCancelationMI(Double_t *ampin, Double_t *ampout, Double_t norm, Double_t lambda, Int_t n);
100           void      FilterSimDeConvExpEl(Int_t *source, Int_t *target, Int_t n, Int_t nexp);
101
102   ClassDef(AliTRDmcmSim,3)
103
104 };
105
106 #endif