TRD/AliTRDmcm.h

   1 #ifndef ALITRDMCM_H
   2 #define ALITRDMCM_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 object                         //
  11 //                                                   //
  12 ///////////////////////////////////////////////////////
  13
  14 #include <TObject.h>
  15
  16 class AliTRDmcm : public TObject {
  17
  18  public:
  19
  20   enum { kMaxTrackletsPerMCM = 4, kMcmCol = 21, kMcmTBmax = 60, kSelClus = 6, kMaxClus = 4 };
  21
  22   AliTRDmcm();
  23   AliTRDmcm(const AliTRDmcm &m);
  24   AliTRDmcm(Int_t id);
  25   virtual         ~AliTRDmcm();
  26   AliTRDmcm       &operator=(const AliTRDmcm &m);
  27
  28   virtual void     Copy(TObject &m) const;
  29
  30           void     SetRobId(Int_t id)                          { fRobId = id; };
  31           void     SetChaId(Int_t id)                          { fChaId = id; };
  32           void     SetRow(Int_t row)                           { fRow = row;  };
  33           void     SetColRange(Int_t colf, Int_t coll)         { fColFirst = colf; fColLast = coll; };
  34           void     SetADC(Int_t icol, Int_t itb, Float_t adc)  { fADC[icol][itb] = adc;       };
  35           void     SetCluster(Int_t icol, Int_t itb)           { fIsClus[icol][itb] = kTRUE;  };
  36           void     UnSetCluster(Int_t icol, Int_t itb)         { fIsClus[icol][itb] = kFALSE; };
  37
  38           Int_t   *GetTrkIndex()                               { return &fTrkIndex[0];        };
  39           Int_t    GetRobId() const { return fRobId; };
  40           Int_t    GetChaId() const { return fChaId; };
  41           Int_t    GetRow() const { return fRow; };
  42           void     GetColRange(Int_t &colf, Int_t &coll) const { colf = fColFirst; coll = fColLast; };
  43           Float_t  GetADC(Int_t icol, Int_t itb) const         { return fADC[icol][itb];     };
  44           Int_t    GetNtrkSeeds() const                        { return fNtrkSeeds;          };
  45           Int_t   *GetSeedCol()                                { return &fSeedCol[0];        };
  46           Int_t   *GetPadHits()                                { return &fPadHits[0];        };
  47           Bool_t   IsCluster(Int_t icol, Int_t itim) const     { return fIsClus[icol][itim]; };
  48           Int_t    Ntrk() const { return fNtrk; };
  49           Int_t    Id() const { return fId; };
  50
  51           void     AddTrk(Int_t id);
  52           void     Reset();
  53           Bool_t   Run();
  54           Bool_t   IsCluster(Float_t amp[3]) const;
  55           void     AddTimeBin(Int_t itime);
  56           Int_t    CreateSeeds();
  57           void     Sort(Int_t nel, Int_t *x1, Int_t *x2, Int_t dir);
  58
  59           void     Filter(Int_t nexp, Int_t ftype = 0);
  60           void     DeConvExpA(Double_t *source, Double_t *target, Int_t n, Int_t nexp);
  61           void     DeConvExpD(Double_t *source, Int_t    *target, Int_t n, Int_t nexp);
  62           void     DeConvExpMI(Double_t *source, Double_t *target, Int_t n);
  63           void     TailMakerSpline(Double_t *ampin, Double_t *ampout, Double_t lambda, Int_t n);
  64           void     TailCancelationMI(Double_t *ampin, Double_t *ampout, Double_t norm, Double_t lambda, Int_t n);
  65
  66  protected:
  67
  68           Int_t    fNtrk;                              //  Number of found tracklets
  69           Int_t    fTrkIndex[kMaxTrackletsPerMCM];     //  Index of found tracklets
  70           Int_t    fRobId;                             //  ROB id
  71           Int_t    fChaId;                             //  Chamber id
  72           Int_t    fRow;                               //  Pad row number (0-11 or 0-15)
  73           Int_t    fColFirst;                          //  First pad column
  74           Int_t    fColLast;                           //  Last pad column (<)
  75           Float_t  fADC[kMcmCol][kMcmTBmax];           //! Array with MCM ADC values
  76           Bool_t   fIsClus[kMcmCol][kMcmTBmax];        //! Flag of a cluster maximum
  77           Int_t    fTime1;                             //  First time bin for tracking (incl.)
  78           Int_t    fTime2;                             //  Last  time bin for tracking (incl.)
  79           Float_t  fClusThr;                           //  Cluster threshold
  80           Float_t  fPadThr;                            //  Pad threshold
  81           Int_t    fPadHits[kMcmCol];                  //  Hit counter in pads
  82           Int_t    fNtrkSeeds;                         //  Number of found seeds
  83           Int_t    fSeedCol[kMaxTrackletsPerMCM];      //  Column number of found tracklet seeds
  84                                                        //  Filter parameters (1 = long, 2 = short component)
  85           Float_t  fR1;                                //  Time constant [microseconds]
  86           Float_t  fR2;                                //  Time constant [microseconds]
  87           Float_t  fC1;                                //  Weight
  88           Float_t  fC2;                                //  Weight
  89           Float_t  fPedestal;                          //  ADC baseline (pedestal)
  90
  91           Int_t    fId;                                //  Dummy id
  92
  93   ClassDef(AliTRDmcm,3)                                //  TRD MCM class
  94
  95 };
  96
  97 #endif