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

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

/* $Id$ */

///////////////////////////////////////////////////////
//                                                   //
//  Multi Chip Module Simulation Class               //
//                                                   //
///////////////////////////////////////////////////////

#include <TObject.h>

class AliTRDfeeParam;
class AliTRDSimParam;
class AliTRDcalibDB;
class AliTRDgeometry;
class AliTRDtrapAlu;
class AliTRDpadPlane;
class AliTRDarrayADC;

class AliTRDmcmSim : public TObject {

 public:

                    AliTRDmcmSim();
                    AliTRDmcmSim(const AliTRDmcmSim &m);
  virtual          ~AliTRDmcmSim();
  AliTRDmcmSim      &operator=(const AliTRDmcmSim &m);

  virtual void      Copy(TObject &m) const;

          void      Init( Int_t cha, Int_t rob, Int_t mcm, Bool_t newEvent );   // Initialize MCM by the position parameters
  //void      Init( Int_t cha, Int_t rob, Int_t mcm );   // Initialize MCM by the position parameters
          void      SetData(Int_t iadc, Int_t *adc);           // Set ADC data with array 
          void      SetData(Int_t iadc, Int_t it, Int_t adc ); // Set ADC data
          void      SetDataPedestal(Int_t iadc );              // Fill ADC data with pedestal values

          Int_t     GetChaId() const  { return fChaId;  };     // Returns Chamber ID (0-539)
          Int_t     GetRobPos() const { return fRobPos; };     // Returns ROB position (0-7)
          Int_t     GetMcmPos() const { return fMcmPos; };     // Returns MCM position (0-17) (16,17 are mergers)
          Int_t     GetRow() const    { return fRow;    };     // Returns Row number on chamber where the MCM is sitting
	  Int_t     GetCol( Int_t iadc );                      // Get corresponding column (0-143) from for ADC channel iadc = [0:20]
	  // for the ADC/Col mapping, see: http://wiki.kip.uni-heidelberg.de/ti/TRD/index.php/Image:ROB_MCM_numbering.pdf

	  Int_t*    GetPosLUT();

	  Int_t     ProduceRawStream( UInt_t *buf, Int_t bufsize );   // Produce raw data stream from this MCM - old
	  Int_t     ProduceRawStreamV2( UInt_t *buf, Int_t bufsize ); // Produce raw data stream - Read data format
	  Int_t     ProduceTrackletStream( UInt_t *buf, Int_t bufsize ); // produce the tracklet stream for this MCM
	  void      Filter();                                  // Apply digital filters for existing data
	  void      ZSMapping();                               // Do ZS mapping for existing data
	  void      DumpData( char *f, char *target );         // Dump data stored (only for debugging)
	  void      Tracklet();
	  void      SetPosLUT();
	  void      CopyArrays();                              // Copy arrays between containers, internal consistency
	  void      GeneratefZSM1Dim();                        // Generate the ZSM1Dim based on digits array info
	  void      StartfastZS(Int_t pads, Int_t timebis);                    // For ZS in the digitizer
	  void      FlagDigitsArray(AliTRDarrayADC *tempdigs, Int_t valrow);   //Set flags on the digits array

 protected:

          Bool_t    fInitialized;                       // Status whether the class is initialized or not
	  Int_t    fNextEvent;                         // 0, if new event; 1 if not
	  Int_t     fMaxTracklets;                      // maximum number of tracklet-words submitted per mcm **new**
	  Int_t     fChaId;                             // Chamber ID
	  Int_t     fSector;                            // Sector number of the supermodule
	  Int_t     fStack;                             // Chamber stack ID
	  Int_t     fLayer;                             // Chamber layer ID
          Int_t     fRobPos;                            // ROB Position on chamber
          Int_t     fMcmPos;                            // MCM Position on chamber
	  Int_t     fNADC;                              // Number of ADC (usually 21)
	  Int_t     fNTimeBin;                          // Number of Timebins (variable)
          Int_t     fRow;                               // Pad row number (0-11 or 0-15) of the MCM on chamber
	  Int_t   **fADCR;                              // Array with MCM ADC values (Raw)
          Int_t   **fADCF;                              // Array with MCM ADC values (Filtered)
	  Int_t   **fADCT;                              // Array with MCM ADC values (filtered) for tracklet (12bits) //***NEW***
	  Int_t    *fPosLUT;                            // position lookup table **new**
	  UInt_t   *fMCMT;                              // tracklet word for one mcm/trap-chip **new**
	  Int_t   **fZSM;                               // Zero suppression map
          Int_t    *fZSM1Dim;                           // Zero suppression map (1 dimensional projection)

	
	  // Parameter classes
	  AliTRDfeeParam *fFeeParam;                    // FEE parameters
	  AliTRDSimParam *fSimParam;                    // Simulation parameters
	  AliTRDcalibDB  *fCal;                         // Calibration interface
	  AliTRDgeometry *fGeo;                         // Geometry

	  Bool_t    CheckInitialized();                 // Check whether the class is initialized
	  void      FilterPedestal();                   // Apply pedestal filter
	  void      FilterGain();                       // Apply gain filter
	  void      FilterTail();                       // Apply tail filter

	  // Here are the several sub functions used only in FilterTail
	  void      FilterSimDeConvExpA(Int_t *source, Double_t *target, Int_t n, Int_t nexp);
	  void      FilterSimDeConvExpD(Int_t *source, Int_t *target, Int_t n, Int_t nexp);
	  void      FilterSimDeConvExpMI(Int_t *source, Double_t *target, Int_t n);
	  void      FilterSimTailMakerSpline(Double_t *ampin, Double_t *ampout, Double_t lambda, Int_t n);
	  void      FilterSimTailCancelationMI(Double_t *ampin, Double_t *ampout, Double_t norm, Double_t lambda, Int_t n);
	  void      FilterSimDeConvExpEl(Int_t *source, Int_t *target, Int_t n, Int_t nexp);

  ClassDef(AliTRDmcmSim,3)

};

#endif
Commit	Line	Data
dfd03fc3	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;
ecf39416	17	class AliTRDSimParam;
ecf39416	18	class AliTRDcalibDB;
dfd03fc3	19	class AliTRDgeometry;
1d93b218	20	class AliTRDtrapAlu;
1d93b218	21	class AliTRDpadPlane;
b65e5048	22	class AliTRDarrayADC;
dfd03fc3	23
	24	class AliTRDmcmSim : public TObject {
	25
	26	public:
	27
ecf39416	28	AliTRDmcmSim();
6e5d4cb2	29	AliTRDmcmSim(const AliTRDmcmSim &m);
ecf39416	30	virtual ~AliTRDmcmSim();
6e5d4cb2	31	AliTRDmcmSim &operator=(const AliTRDmcmSim &m);
dfd03fc3	32
6e5d4cb2	33	virtual void Copy(TObject &m) const;
dfd03fc3	34
96e6312d	35	void Init( Int_t cha, Int_t rob, Int_t mcm, Bool_t newEvent ); // Initialize MCM by the position parameters
96e6312d	36	//void Init( Int_t cha, Int_t rob, Int_t mcm ); // Initialize MCM by the position parameters
dfd03fc3	37	void SetData(Int_t iadc, Int_t *adc); // Set ADC data with array
ecf39416	38	void SetData(Int_t iadc, Int_t it, Int_t adc ); // Set ADC data
dfd03fc3	39	void SetDataPedestal(Int_t iadc ); // Fill ADC data with pedestal values
dfd03fc3	40
0c349049	41	Int_t GetChaId() const { return fChaId; }; // Returns Chamber ID (0-539)
	42	Int_t GetRobPos() const { return fRobPos; }; // Returns ROB position (0-7)
	43	Int_t GetMcmPos() const { return fMcmPos; }; // Returns MCM position (0-17) (16,17 are mergers)
	44	Int_t GetRow() const { return fRow; }; // Returns Row number on chamber where the MCM is sitting
ecf39416	45	Int_t GetCol( Int_t iadc ); // Get corresponding column (0-143) from for ADC channel iadc = [0:20]
ecf39416	46	// for the ADC/Col mapping, see: http://wiki.kip.uni-heidelberg.de/ti/TRD/index.php/Image:ROB_MCM_numbering.pdf
dfd03fc3	47
1d93b218	48	Int_t* GetPosLUT();
1d93b218	49
987ba9a3	50	Int_t ProduceRawStream( UInt_t *buf, Int_t bufsize ); // Produce raw data stream from this MCM - old
987ba9a3	51	Int_t ProduceRawStreamV2( UInt_t *buf, Int_t bufsize ); // Produce raw data stream - Read data format
1d93b218	52	Int_t ProduceTrackletStream( UInt_t *buf, Int_t bufsize ); // produce the tracklet stream for this MCM
ecf39416	53	void Filter(); // Apply digital filters for existing data
	54	void ZSMapping(); // Do ZS mapping for existing data
	55	void DumpData( char f, char target ); // Dump data stored (only for debugging)
1d93b218	56	void Tracklet();
1d93b218	57	void SetPosLUT();
b65e5048	58	void CopyArrays(); // Copy arrays between containers, internal consistency
	59	void GeneratefZSM1Dim(); // Generate the ZSM1Dim based on digits array info
	60	void StartfastZS(Int_t pads, Int_t timebis); // For ZS in the digitizer
	61	void FlagDigitsArray(AliTRDarrayADC *tempdigs, Int_t valrow); //Set flags on the digits array
dfd03fc3	62
	63	protected:
	64
ecf39416	65	Bool_t fInitialized; // Status whether the class is initialized or not
96e6312d	66	Int_t fNextEvent; // 0, if new event; 1 if not
1d93b218	67	Int_t fMaxTracklets; // maximum number of tracklet-words submitted per mcm new
dfd03fc3	68	Int_t fChaId; // Chamber ID
	69	Int_t fSector; // Sector number of the supermodule
	70	Int_t fStack; // Chamber stack ID
	71	Int_t fLayer; // Chamber layer ID
ecf39416	72	Int_t fRobPos; // ROB Position on chamber
ecf39416	73	Int_t fMcmPos; // MCM Position on chamber
dfd03fc3	74	Int_t fNADC; // Number of ADC (usually 21)
ecf39416	75	Int_t fNTimeBin; // Number of Timebins (variable)
ecf39416	76	Int_t fRow; // Pad row number (0-11 or 0-15) of the MCM on chamber
1d93b218	77	Int_t **fADCR; // Array with MCM ADC values (Raw)
dfd03fc3	78	Int_t **fADCF; // Array with MCM ADC values (Filtered)
1d93b218	79	Int_t fADCT; // Array with MCM ADC values (filtered) for tracklet (12bits) //NEW**
	80	Int_t fPosLUT; // position lookup table new*
	81	UInt_t fMCMT; // tracklet word for one mcm/trap-chip new*
	82	Int_t **fZSM; // Zero suppression map
ecf39416	83	Int_t *fZSM1Dim; // Zero suppression map (1 dimensional projection)
ecf39416	84
1d93b218	85
	86
	87
	88
0c349049	89	// Parameter classes
	90	AliTRDfeeParam *fFeeParam; // FEE parameters
	91	AliTRDSimParam *fSimParam; // Simulation parameters
	92	AliTRDcalibDB *fCal; // Calibration interface
	93	AliTRDgeometry *fGeo; // Geometry
dfd03fc3	94
ecf39416	95	Bool_t CheckInitialized(); // Check whether the class is initialized
dfd03fc3	96	void FilterPedestal(); // Apply pedestal filter
	97	void FilterGain(); // Apply gain filter
	98	void FilterTail(); // Apply tail filter
	99
ecf39416	100	// Here are the several sub functions used only in FilterTail
dfd03fc3	101	void FilterSimDeConvExpA(Int_t source, Double_t target, Int_t n, Int_t nexp);
	102	void FilterSimDeConvExpD(Int_t source, Int_t target, Int_t n, Int_t nexp);
	103	void FilterSimDeConvExpMI(Int_t source, Double_t target, Int_t n);
	104	void FilterSimTailMakerSpline(Double_t ampin, Double_t ampout, Double_t lambda, Int_t n);
	105	void FilterSimTailCancelationMI(Double_t ampin, Double_t ampout, Double_t norm, Double_t lambda, Int_t n);
1d93b218	106	void FilterSimDeConvExpEl(Int_t source, Int_t target, Int_t n, Int_t nexp);
dfd03fc3	107
	108	ClassDef(AliTRDmcmSim,3)
	109
	110	};
	111
	112	#endif