onlinemonitor/rocmonitor/TFiberHodProc.h (r4864/r3062)

Go to the documentation of this file.

#ifndef TFIBERHODPROC_H
#define TFIBERHODPROC_H

#include "TCBMBeamtimeProc.h"

#include "TFiberHodEvent.h"

class TFiberHodParam;
class TRocEvent;

class TFiberHodProc : public TCBMBeamtimeProc {
   public:
      TFiberHodProc(); // required for ROOT streamer
      TFiberHodProc(const char* name);
      virtual ~TFiberHodProc();

      virtual void InitEvent(TGo4EventElement*);
      virtual void FinalizeEvent();

   protected:

      void InitializeDecoding(Int_t altflag=0);

      void PrintAllHits();

      void PrintAllClusters();

      // Structures for defining decoding
      // -------------------------------------------------------------

      // struct to build decoding table, stored s vector fDecoding
      struct decodingentry_t {
         Int_t FEBchannel;      //readout channel [0..127]
         Int_t fiber;           //fiber number [1..64]
         Int_t plane;           //[1..2],  1=X, 2=y
         Int_t pixel;           //[1..64]
         Int_t lcn;             //[1..16]
         Int_t cable;           //[1..4] linear cable nr
         Double_t position;     // in [mm], 0-point in middle of array, between fibers 32 and 33
      };

      // struct to build hit bank, stored as fHitbank
      struct hitbankentry_t {
         Int_t FEBchannel;
         Int_t fiber;
         Int_t plane;
         Int_t pixel;
         Double_t ADC;
         Double_t ADCraw;
         Double_t time;
#ifdef HODOSCOPE_FULLTIMESTAMP
         uint64_t fulltime; // absolute ts
#endif
         Double_t position;
         Int_t status;         // 0:not clustered, 1:already clustered, 2:crosstalk
      };

      // comparison operator to be able to sort hitbank according to increasing plane and decreasing adc value
      struct adccompare {
        bool operator() ( const hitbankentry_t a, const hitbankentry_t b)
        {
           if (a.plane != b.plane) return a.plane < b.plane;
           return b.ADC < a.ADC;
        }
      };

      // struct to build clusterbank, stored as fClusterbank
      struct clusterbankentry_t {
         Int_t   NrHits;
         Int_t   plane;
         Double_t meanposition;
         Double_t meanfiber;
         Double_t sumADC;
         Double_t meantime;
#ifdef HODOSCOPE_FULLTIMESTAMP
         uint64_t meanfulltime; // absolute ts
#endif
      };

      // define data structures: decoding, hitbank, clusterbank using static vectors
      // ---------------------------------------------------------------------------

      TFiberHodParam* fParam; 


      // define geometry decoding table
      decodingentry_t  fDecoding[128]; 

      TRocEvent*       fRocEvent;  
      TFiberHodEvent*  fHodEvent;  

      std::vector<hitbankentry_t> fHitbank;   

      std::vector <clusterbankentry_t> fClusterbank; 

      Int_t            fSingleEventSlice;

      // Histogram pointers
      TH2  *hSingleEventHit2D; 
      TH2  *hSingleEventCluster2D[5];

      TH1  *hMaxCh;
      TH1  *hMaxFiberall;
      TH1  *hMaxFiberpl1;
      TH1  *hMaxFiberpl2;
      TH1  *hADC;  
      TH2  *hMaxADC;
      TH2  *hMaxADCraw;
      TH2  *hMaxFiber2D;
      TH1  *hNrHits;
      TH1  *hNrClusters;
      TH2  *hNrClusters2D;
      TH1  *hMaxClusterTimediff;
      TH1  *hPosition2D;

      // Conditions
      TGo4WinCond   *fADC_cond;    

   ClassDef(TFiberHodProc,1)
};

#endif