[u/mrichter/AliRoot.git] / TRD / AliTRDmcm.h

#ifndef ALITRDMCM_H
#define ALITRDMCM_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */

///////////////////////////////////////////////////////
//                                                   //
//  Multi Chip Module object                         //
//                                                   //
///////////////////////////////////////////////////////

#include <TObject.h>

class AliTRDmcm : public TObject {

 public:

  enum { kMaxTrackletsPerMCM = 4, kMcmCol = 21, kMcmTBmax = 60, kSelClus = 6, kMaxClus = 4 };

  AliTRDmcm();
  AliTRDmcm(const AliTRDmcm &m);
  AliTRDmcm(Int_t id);
  virtual         ~AliTRDmcm();
  AliTRDmcm       &operator=(const AliTRDmcm &m);

  virtual void     Copy(TObject &m) const;

          void     SetRobId(Int_t id)                          { fRobId = id; };
          void     SetChaId(Int_t id)                          { fChaId = id; };
          void     SetRow(Int_t row)                           { fRow = row;  };
          void     SetColRange(Int_t colf, Int_t coll)         { fColFirst = colf; fColLast = coll; };
          void     SetADC(Int_t icol, Int_t itb, Float_t adc)  { fADC[icol][itb] = adc;       };
          void     SetCluster(Int_t icol, Int_t itb)           { fIsClus[icol][itb] = kTRUE;  };
          void     UnSetCluster(Int_t icol, Int_t itb)         { fIsClus[icol][itb] = kFALSE; };

          Int_t   *GetTrkIndex()                               { return &fTrkIndex[0];        };
          Int_t    GetRobId() const { return fRobId; };
          Int_t    GetChaId() const { return fChaId; };
          Int_t    GetRow() const { return fRow; };
          void     GetColRange(Int_t &colf, Int_t &coll) const { colf = fColFirst; coll = fColLast; };
          Float_t  GetADC(Int_t icol, Int_t itb) const         { return fADC[icol][itb];     };
          Int_t    GetNtrkSeeds() const                        { return fNtrkSeeds;          };
          Int_t   *GetSeedCol()                                { return &fSeedCol[0];        };
          Int_t   *GetPadHits()                                { return &fPadHits[0];        };
          Bool_t   IsCluster(Int_t icol, Int_t itim) const     { return fIsClus[icol][itim]; };
          Int_t    Ntrk() const { return fNtrk; };
          Int_t    Id() const { return fId; };

          void     AddTrk(Int_t id);
          void     Reset();
          Bool_t   Run();
          Bool_t   IsCluster(Float_t amp[3]) const;
          void     AddTimeBin(Int_t itime);
          Int_t    CreateSeeds();
          void     Sort(Int_t nel, Int_t *x1, Int_t *x2, Int_t dir);

          void     Filter(Int_t nexp, Int_t ftype = 0);
          void     DeConvExpA(Double_t *source, Double_t *target, Int_t n, Int_t nexp);
          void     DeConvExpD(Double_t *source, Int_t    *target, Int_t n, Int_t nexp);
          void     DeConvExpMI(Double_t *source, Double_t *target, Int_t n);
          void     TailMakerSpline(Double_t *ampin, Double_t *ampout, Double_t lambda, Int_t n);
          void     TailCancelationMI(Double_t *ampin, Double_t *ampout, Double_t norm, Double_t lambda, Int_t n);

 protected:

          Int_t    fNtrk;                              //  Number of found tracklets
          Int_t    fTrkIndex[kMaxTrackletsPerMCM];     //  Index of found tracklets
          Int_t    fRobId;                             //  ROB id
          Int_t    fChaId;                             //  Chamber id
          Int_t    fRow;                               //  Pad row number (0-11 or 0-15)
          Int_t    fColFirst;                          //  First pad column
          Int_t    fColLast;                           //  Last pad column (<)
          Float_t  fADC[kMcmCol][kMcmTBmax];           //! Array with MCM ADC values
          Bool_t   fIsClus[kMcmCol][kMcmTBmax];        //! Flag of a cluster maximum
          Int_t    fTime1;                             //  First time bin for tracking (incl.)
          Int_t    fTime2;                             //  Last  time bin for tracking (incl.)
          Float_t  fClusThr;                           //  Cluster threshold
          Float_t  fPadThr;                            //  Pad threshold
          Int_t    fPadHits[kMcmCol];                  //  Hit counter in pads
          Int_t    fNtrkSeeds;                         //  Number of found seeds
          Int_t    fSeedCol[kMaxTrackletsPerMCM];      //  Column number of found tracklet seeds
                                                       //  Filter parameters (1 = long, 2 = short component)
          Float_t  fR1;                                //  Time constant [microseconds] 
          Float_t  fR2;                                //  Time constant [microseconds]
          Float_t  fC1;                                //  Weight
          Float_t  fC2;                                //  Weight
          Float_t  fPedestal;                          //  ADC baseline (pedestal)

          Int_t    fId;                                //  Dummy id
 
  ClassDef(AliTRDmcm,3)                                //  TRD MCM class

};

#endif
Commit	Line	Data
0ee00e25	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	///////////////////////////////////////////////////////
e3b2b5e5	9	// //
0ee00e25	10	// Multi Chip Module object //
e3b2b5e5	11	// //
0ee00e25	12	///////////////////////////////////////////////////////
	13
	14	#include <TObject.h>
	15
0ee00e25	16	class AliTRDmcm : public TObject {
	17
	18	public:
	19
e3b2b5e5	20	enum { kMaxTrackletsPerMCM = 4, kMcmCol = 21, kMcmTBmax = 60, kSelClus = 6, kMaxClus = 4 };
e3b2b5e5	21
0ee00e25	22	AliTRDmcm();
6d50f529	23	AliTRDmcm(const AliTRDmcm &m);
f162af62	24	AliTRDmcm(Int_t id);
6d50f529	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);
0ee00e25	65
	66	protected:
	67
6d50f529	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
f162af62	93	ClassDef(AliTRDmcm,3) // TRD MCM class
0ee00e25	94
	95	};
	96
	97	#endif