STEER/TFndRdt.h

// @(#)fROOT/STEER:$Name:  $:$Id: TFndRdt.h,v 1.23 2007/09/05 11:09:42 Diego_Faso Exp $
// Revision Author: Diego Faso <mailto:faso@to.infn.it>, 2005/06/24
// Original Author: P. Cerello <mailto:cerello@to.infn.it>, 2000

#ifndef FROOT_TFndRdt
#define FROOT_TFndRdt

//                                                                      //
//                               TFndRdt                                //
//                                                                      //
// This class gives the Description of the Raw Data Tape event.         //
// The private part includes some data members                          //
// (event length, pointers to the whole buffer and the subdetector      //
// buffers) and a set of member functions retrieving some information   //
// needed to build the Hit Data Tape.                                   //
// The public part allows to dump the event and invoke the member       //
// functions which transform the RDT sub-events into HDT sub-events.    //
//                                                                      //

#include "TObject.h"
#include "TFndHdt.h"
#include "TFndFeeMap.h"

#include "FROOT.h"
using namespace FROOT;

typedef struct {
  Int_t Layer; // Layer 1 = ISIM 2 = OSIM 
  Int_t Module; // Module number ISIM (1-8) OSIM (1-10)
  Int_t Side;  // (See FROOT::E_Fnd_Sil_Side)
  Int_t Strip; // Strip number
  Int_t Spy;   // 1=Spy
  Int_t Noise; // Noise level: 0=OK ... 6=dead
  Int_t AdcCh; // cram channel number
  Int_t AdcVal;// adc value
} SILADCHIT;

typedef struct {
  Int_t    NumStripCl; // Number of Strip of Cluster
  Int_t    AdcChCl;  // Center of Mass of Cluster (adc channel) 
  Double_t StripChCl; // Center of Mass of Cluster
  Double_t NoiseCl;  // Noise level of Cluster
  Double_t AdcValCl; // Sum Adc Values of Cluster
} SILCLRAW;

// Ascending adc val cluster sorting function
struct SAscendingAdcValClSort
{
     bool operator()(const SILCLRAW& rpStart,const SILCLRAW& rpEnd)
     {
          return rpStart.AdcValCl < rpEnd.AdcValCl;
     }
};

// Descending adc val cluster sorting function
struct SDescendingAdcValClSort
{
     bool operator()(const SILCLRAW& rpStart,const SILCLRAW& rpEnd)
     {
          return rpStart.AdcValCl > rpEnd.AdcValCl;
     }
};


class TFndRdt : public TObject {

 public:
  enum E_FIN_RDT_DetID { FIN_RDT_Eqp_Gts = 0,
                         FIN_RDT_Eqp_Tof = 1, 
                         FIN_RDT_Eqp_Isi = 2, 
                         FIN_RDT_Eqp_Osi = 3, 
                         FIN_RDT_Eqp_Lmd = 4,
                         FIN_RDT_Eqp_Stb = 5 
  };

  enum E_FIN_RDT_EqpID { 
    FIN_EQ_CORBO =         0,  // corbo (CES RCB8047)
    FIN_EQ_CAEN_V513 =     1,  // vme i/o register (CAEN V513)
    FIN_EQ_TRIG_SEL =      2,  // vme trigger selector
    FIN_EQ_CAEN_V529 =     3,  // vme pattern unit (CAEN V259)
    FIN_EQ_CAEN_V560 =     4,  // vme scaler (CAEN V560)
    FIN_EQ_CAEN_VN1485 =   5,  // vme adc 9U (CAEN VN1485)
    FIN_EQ_CAEN_VN1488 =   6,  // vme tdc 9U (CAEN VN1488)
    FIN_EQ_LECROY_1277 =   7,  // vme tdc 9U (LeCroy 1277)
    FIN_EQ_ADC_6U =        8,  // vme adc 6U(CAEN xxxxx )
    FIN_EQ_CAEN_673A =     9,  // vme tdc 6U64 (CAEN 673A )
    FIN_EQ_LECROY_4532 =  20, // malu (LeCroy 4532)
    FIN_EQ_CAMAC_IO_REG = 21, // camac i/o register (CAEN )
    FIN_EQ_CAMAC_SCALER = 22, // camac scaler (CAEN C285 or LeCroy 2251)
    FIN_EQ_LeCroy_2249 =  23, // camac adc (LeCroy 2249)
    FIN_EQ_LeCroy_2277 =  24, // camac tdc (LeCroy 2277)
    FIN_EQ_LeCroy_2229 =  25 // camac tdc (LeCroy 2229)
  };
  
 private:
  
  Int_t fDebugLev; // taken from TFndRun
  Int_t fCurRunType; // see TFndRun and FROOT.h
  TFndFeeMap *fFndFeeMap; // maybe it will be replaced by a (coming soon) global pointer to TFndDB
  //  TFndHdt *fCurHdtEv; // current HDT event taken from TFndRun (do not delete it here)

  UInt_t *fRawEvent; // pointer to the (complete) raw event
  
  std::vector<SILADCHIT> fSilHits;
  Bool_t fMustEvalCommonMode;
  Bool_t fMustFindSilClusters;

  Int_t fSilAdcVal[(Int_t)E_FIN_OUTER_LAYER+1][(Int_t)K_N_Sil_Modules+1][(Int_t)K_N_Sil_Sides-1][(Int_t)K_N_Sil_Strips+1]; 
  Int_t fSilAdcNoi[(Int_t)E_FIN_OUTER_LAYER+1][(Int_t)K_N_Sil_Modules+1][(Int_t)K_N_Sil_Sides-1][(Int_t)K_N_Sil_Strips+1]; 
  Int_t fSilAdcCh[(Int_t)E_FIN_OUTER_LAYER+1][(Int_t)K_N_Sil_Modules+1][(Int_t)K_N_Sil_Sides-1][(Int_t)K_N_Sil_Strips+1]; 

  void SetAdc_Val_Noi_Zero() {
    for(Int_t nlay = 1; nlay < (Int_t) E_FIN_OUTER_LAYER+1; nlay++)
      for(Int_t nmod = 1; nmod < (Int_t) K_N_Sil_Modules+1; nmod++)
        for(Int_t nside = 0; nside < (Int_t) K_N_Sil_Sides - 1; nside++)
          for(Int_t nstrip = 1; nstrip < (Int_t) K_N_Sil_Strips + 1; nstrip++) {
            SetSilAdcVal(nlay,nmod,nside,nstrip,0);
            SetSilAdcNoi(nlay,nmod,nside,nstrip,0);
            SetSilAdcCh(nlay,nmod,nside,nstrip,0);
          }
  }
  //
  void SetSilAdcVal(Int_t Sil_Layer,Int_t Sil_Module,Int_t Sil_Side,Int_t Sil_Nstrip,Int_t val) { fSilAdcVal[Sil_Layer][Sil_Module][Sil_Side][Sil_Nstrip] = val; }
  void SetSilAdcNoi(Int_t Sil_Layer,Int_t Sil_Module,Int_t Sil_Side,Int_t Sil_Nstrip,Int_t val) { fSilAdcNoi[Sil_Layer][Sil_Module][Sil_Side][Sil_Nstrip] = val; }
  void SetSilAdcCh(Int_t Sil_Layer,Int_t Sil_Module,Int_t Sil_Side,Int_t Sil_Nstrip,Int_t val) { fSilAdcCh[Sil_Layer][Sil_Module][Sil_Side][Sil_Nstrip] = val; }
  Int_t GetSilAdcVal(Int_t Sil_Layer,Int_t Sil_Module,Int_t Sil_Side,Int_t Sil_Nstrip) {return fSilAdcVal[Sil_Layer][Sil_Module][Sil_Side][Sil_Nstrip]; }
  Int_t GetSilAdcNoi(Int_t Sil_Layer,Int_t Sil_Module,Int_t Sil_Side,Int_t Sil_Nstrip) {return fSilAdcNoi[Sil_Layer][Sil_Module][Sil_Side][Sil_Nstrip]; }
  Int_t GetSilAdcCh(Int_t Sil_Layer,Int_t Sil_Module,Int_t Sil_Side,Int_t Sil_Nstrip) {return fSilAdcCh[Sil_Layer][Sil_Module][Sil_Side][Sil_Nstrip]; }
  /*   Double_t SilPositionStrip(Int_t side, Double_t barcha); */
  //
  void FillLmdAdcHits();
  void FillLmdTdcHits();
  void FillTofAdcHits(UInt_t *det);
  void FillTofTdcHits(UInt_t *det);
  void EvalCommonModeNoise();
  void FindSilClusters();
  
  UInt_t *GetDetector(E_FIN_RDT_DetID det_id);
  UInt_t *GetEquipment(E_FIN_RDT_DetID det_id,Int_t eqp_num);

  Int_t FillGtsHit(); // return 0 in case of normal buffer
  void  FillTofHits();
  void  FillSilHits(Int_t io);
  void  FillLmdHits();
  void  FillStbHits();
  
 public:
  
  TFndRdt();            // Default constructor
  ~TFndRdt();           // Destructor
  
  void Reset(); // Come back to the just created object
  void Reload(); // Reload needed info from current fndrun (Reset method is called here)
  
  void RdtToHdt(UInt_t *raw,TFndHdt *fev); // Fill the HDT structure for the current event

  void PrintEquipment(Int_t eqp_id, Int_t neqp); // Print the content of Equipment number "neqp"  
  ClassDef(TFndRdt,0)  // Raw event decoding (and HDT filling)

};

#endif // FROOT_TFndRdt

---