StMtdHitMaker.h

#ifndef STAR_StMtdHitMaker_H
#define STAR_StMtdHitMaker_H

/***************************************************************************
 *
 * $Id: StMtdHitMaker.h,v 1.13 2015/01/22 22:10:58 marr Exp $ 
 * StMtdHitMaker - class to fill StEvent from DAQ reader
 *--------------------------------------------------------------------------
 *
 ***************************************************************************/
#include "StRTSBaseMaker.h"
#include "TH3D.h"
#include "StMtdUtil/StMtdConstants.h"
struct mtd_t;

class StEvent;
class StMtdCollection;
class StMtdRawHit;
class StMtdHit;
class StBTofINLCorr;

// 
#if !defined(ST_NO_TEMPLATE_DEF_ARGS) || defined(__CINT__)
typedef vector<Int_t>  IntVec;
typedef vector<UInt_t>  UIntVec;
#else
typedef vector<Int_t, allocator<Int_t>>  IntVec;
typedef vector<UInt_t, allocator<UInt_t>>  UIntVec;
#endif

struct MtdRawHit {
  UInt_t   tdc;               
  UInt_t   dataword;          
  UChar_t  fiberid;           
  UChar_t  backlegid;            
  UChar_t  globaltdcchan;     
};

struct MTDSingleHit{
        Int_t                   fiberId;
        Int_t           backleg; // 1-30
        Int_t           tray;  // 1-5
        Int_t           channel; // 0-23
        IntVec                  leadingEdgeTime;     
        IntVec                  trailingEdgeTime;     
};

struct MTDOneSideHit{
        Int_t           fiberId;
        Int_t           backleg; // 1-30
        Int_t           tray;  // 1-5
        Int_t           channel; // 0-23
        Double_t        leadingEdgeTime;     
        Double_t                trailingEdgeTime;     
};

#ifndef ST_NO_TEMPLATE_DEF_ARGS
  typedef vector<MTDSingleHit> mtdSingleHitVector;
#else
  typedef vector<MTDSingleHit,allocator<MTDSingleHit>> mtdSingleHitVector;
#endif
  typedef vector<MTDSingleHit>::iterator mtdSingleHitVectorIter;

#define VHRBIN2PS 24.4140625  // Very High resolution mode, pico-second per bin
                              // 1000*25/1024 (ps/chn)
#define HRBIN2PS 97.65625     // High resolution mode, pico-second per bin
                              // 97.65625= 1000*100/1024  (ps/chn)

class StMtdHitMaker:public StRTSBaseMaker {
 private:
  StEvent                 *mStEvent;
  mtd_t                   *fMtd;
  Int_t                   mUseMuDst;
  Bool_t                  mCosmicFlag;
  TString                 mCosmicTrigTimeWinFile;
  Bool_t                  mTriggerWndSelection;
  Int_t                   mSwapBacklegInRun13;                         // 0 - do not swap; 1 - for first part of Run13
                                                                       // 2 - for second part of Run13
  Bool_t                  mReverseStripInRun14;                        // Reverse strips in backleg 7, module 5 in Run14
  Int_t                   mYear;                                       
  Int_t                   mNValidTrays;                                

  StMtdCollection*        mMtdCollection;                              
  StBTofINLCorr*          mINLCorr;                                    

  static const Int_t nTHUB     = 2;
  static const Int_t mNFIBER   = 1;
  UInt_t                  mTriggerTimeStamp[nTHUB];                    
  Int_t                   mTray2TdigMap[gMtdNBacklegs][gMtdNModules];  
  Int_t                   mTrayId[gMtdNBacklegs][gMtdNModules];        
  Int_t                   mTdigId[gMtdNModulesAll];                    
  Int_t                   mtdStrip[gMtdNChannels];                     
  Double_t                mTriggerTimeWindow[gMtdNModulesAll][2];      

  vector<MtdRawHit> MtdLeadingHits;
  vector<MtdRawHit> MtdTrailingHits;
  mtdSingleHitVector mSingleHitVec[gMtdNModulesAll];

  TH3D *hxhyhz;


 protected:
  StMtdCollection *GetMtdCollection();
  Int_t UnpackMtdRawData();
  StRtsTable *GetNextRaw();
  StRtsTable *GetNextRaw(Int_t sec);
  Int_t tdcChan2globalStrip11(Int_t, Int_t, Int_t, Int_t);
  Int_t tdcChan2globalStrip(Int_t, Int_t, Int_t, Int_t);
  Int_t getTdigBoardId(Int_t backlegid, Int_t tray, Int_t chn);
  Int_t getLocalTdcChan(Int_t backlegid, Int_t tray, Int_t chn);
  Int_t getTdigLocalChan(Int_t backlegid, Int_t itdigid);
  void fillMtdHeader();
  void fillMtdRawHitCollection();
  void fillMtdHitCollection();
  void fillMtdSingleHits();
  void fillStEvent();     
  IntVec  GetValidTrays();
  IntVec  GetValidChannel(int backleg, int tray, int &fiber);
  UIntVec GetLeadingTdc(int backleg, int tray, int channel);
  UIntVec GetTrailingTdc(int backleg, int tray, int channel);

 public:

  StMtdHitMaker(const char *name="mtd_raw");     
  ~StMtdHitMaker() ;
  void setUseMuDst(Int_t val);
  void setTriggerWndSelection(Bool_t val);
  void setSwapBacklegInRun13(Int_t swap);
  void setReverseStripInRun14(Bool_t re);
  void setCosmicEvent(Bool_t val);
  void setCosmicTrigTimeWinFile(const char *file);

  void   Clear(Option_t* option="");
  Int_t  Init();
  Int_t  InitRun(Int_t);
  Int_t  FinishRun(Int_t);
  Int_t  Finish();
  Int_t  Make();

  vector<MtdRawHit> getLeadingHits();
  vector<MtdRawHit> getTrailingHits();
  
  virtual const char *GetCVS() const {
    static const char cvs[]="Tag $Name:  $Id: built " __DATE__ " " __TIME__ ; return cvs;
  }
  
  ClassDef(StMtdHitMaker, 4)    
};

inline vector<MtdRawHit> StMtdHitMaker::getLeadingHits()      { return MtdLeadingHits; }
inline vector<MtdRawHit> StMtdHitMaker::getTrailingHits()     { return MtdTrailingHits;}
inline void StMtdHitMaker::setUseMuDst(Int_t val)             { mUseMuDst = val;}
inline void StMtdHitMaker::setTriggerWndSelection(Bool_t val) { mTriggerWndSelection = val;}
inline void StMtdHitMaker::setSwapBacklegInRun13(Int_t swap)  { mSwapBacklegInRun13 = swap; }
inline void StMtdHitMaker::setReverseStripInRun14(Bool_t re)  { mReverseStripInRun14 = re; }
inline void StMtdHitMaker::setCosmicEvent(Bool_t val)         { mCosmicFlag = val; }
inline void StMtdHitMaker::setCosmicTrigTimeWinFile(const char *file)
{ mCosmicTrigTimeWinFile = file; }

#endif