StTofpMatchMaker.h

/*  -------------------------------------------------------------------------
 * $Log: StTofpMatchMaker.h,v $
 * Revision 1.9  2014/08/06 11:43:49  jeromel
 * Suffix on literals need to be space (later gcc compiler makes it an error) - first wave of fixes
 *
 * Revision 1.8  2005/07/06 22:12:09  fisyak
 * Use templated StThreeVectorD
 *
 * Revision 1.7  2005/04/12 17:32:45  dongx
 * update for year 5 data. Tofp removed, so do nothing in this maker from now on.
 *
 * Revision 1.6  2004/06/09 21:28:05  dongx
 * update matching : checking before projecting track, improve the speed by around an order of magnitude
 *
 * Revision 1.5  2004/03/11 22:29:32  dongx
 * -remove assert()
 * -add member mYear4
 * -use m_Mode to control the output root file
 *
 * Revision 1.4  2003/09/15 22:38:11  geurts
 * dBase updates:
 *  - removed initLocalDb option
 *  - introduced dBase parameters for strobe event definitions
 *
 * Revision 1.3  2003/09/13 19:15:52  geurts
 * Changed passing of StSPtrVecTofData for strobeEvent (see bugtracker ticket #172)
 *
 * Revision 1.2  2003/09/10 19:47:38  perev
 * ansi corrs
 *
 * Revision 1.1  2003/08/08 18:31:26  geurts
 * first release
 *
 *
 * -------------------------------------------------------------------------
 */
#ifndef STTOFPMATCHMAKER_HH     
#define STTOFPMATCHMAKER_HH
#include "StMaker.h"
#include <string>

class StEvent;
class StTrack;
class StHelix;
#include "StThreeVectorD.hh"
class StTrackGeometry;
class StTofGeometry;
class StTofCollection;
class StTofSlatCollection;
class StSPtrVecTofData;
class TOrdCollection;
class TH1D;
class TH2D;



class StTofpMatchMaker : public StMaker {
  static const Int_t NTOFP = 41;
  static const Int_t NPVPD = 6;
public:
  StTofpMatchMaker(const Char_t *name="tofpMatch");
  ~StTofpMatchMaker();

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

  void createHistograms(Bool_t histos=kTRUE);
  void setOuterTrackGeometry();
  void setStandardTrackGeometry();
  void setValidAdcRange(Int_t, Int_t);
  void setValidTdcRange(Int_t, Int_t);
  void setMinHitsPerTrack(Int_t);
  void setMinFitPointsPerTrack(Int_t);
  void setMaxDCA(Float_t);
  void setHistoFileName(const Char_t*);

  Bool_t doPrintMemoryInfo; 
  Bool_t doPrintCpuInfo; 

protected:
  StTrackGeometry* trackGeometry(StTrack*);
  Float_t slatPropagationTime(StThreeVectorD*); // calculate hit position correction
  Float_t startTime(const float); // calculated pvpd startTime
  Int_t getTofData(StTofCollection*); // check, remap and fill local arrays with tof and pvpd data
  Int_t storeMatchData(StTofSlatCollection*, StTofCollection*); 
  Float_t mTofpAdc[NTOFP];
  Float_t mTofpTdc[NTOFP];
  Float_t mPvpdAdc[NPVPD];
  Float_t mPvpdAdcLoRes[NPVPD];
  Float_t mPvpdTdc[NPVPD];


private:
  Bool_t strobeEvent(StSPtrVecTofData&);// check pVPD data for strobe event
  Int_t mStrobeTdcMin[NPVPD]; 
  Int_t mStrobeTdcMax[NPVPD]; 

  void bookHistograms();
  void writeHistogramsToFile();
  StTofGeometry *mTofGeom; 

  static const Int_t mTofpTrayId = 93;   // tofp tray id number in year 4
  Bool_t mHisto; 
  //Bool_t mInitLocalDb; //! initialize from local (true) dbase or STAR dbase (false)
  Bool_t mYear2; 
  Bool_t mYear3; 
  Bool_t mYear4; 
  Bool_t mYear5; 
  Bool_t mOuterTrackGeometry; 
  string mHistoFileName; 

  Bool_t validAdc(float);
  Bool_t validTdc(float);
  Bool_t validEvent(StEvent *);
  Bool_t validTrack(StTrack*);
  Bool_t validTofTrack(StTrack*);

  // event counters
  Int_t mEventCounter;          
  Int_t mAcceptedEventCounter;  
  Int_t mTofEventCounter;       
  Int_t mTofStrobeEventCounter; 
  Int_t mAcceptAndStrobe;       
  Int_t mAcceptAndBeam;         

  // various cut-offs and ranges
  Float_t mMinValidTdc; 
  Float_t mMaxValidTdc; 
  Float_t mMinValidAdc; 
  Float_t mMaxValidAdc; 
  UInt_t mMinHitsPerTrack; 
  UInt_t mMinFitPointsPerTrack; 
  Float_t mMaxDCA; 


  // TOFp histogram collections
  TOrdCollection* mHitPosHistNames;
  TOrdCollection* mTrackHistNames; 
  TOrdCollection* mOccupancyHistNames; 

  TH2D *hTofpHitMap1; 
  TH2D *hTofpHitMap2; 
  TH2D *hTofpHitMap3; 
  TH2D *hTofpHitMap4; 
  TH1D *hTofpSlatIdA0; 
  TH1D *hTofpSlatIdA1; 
  TH1D *hTofpSlatIdB1; 
  TH1D *hTofpSlatIdD1; 
  TH1D *hTofpSlatIdD2; 
  TH1D *hTofpSlatIdE1; 
  TH1D *hTofpSlatIdE2; 
  TH1D *hTofpSlatIdE3; 
  TH1D *hTofpSlatIdE4; 
  TH1D *hTofpSlatIdE5; 
  TH1D *hTofpSlatIdF1; 
  TH1D *hTofpSlatHitVecSize;     

  TH1D *hTofpNumberOfTrackHits;  
  TH1D *hTofpPtTrack;            
  TH1D *hTofpDCATrackprimVertex; 


  // Occupancy plots
  TH1D *hTofpNumberOfValidAdc;        
  TH1D *hTofpNumberOfValidTdc;        
  TH1D *hTofpNumberOfValidSlats;      
  TH1D *hTofpNumberOfGlobalTracks;    
  TH1D *hTofpNumberOfHitSlats;        
  TH1D *hTofpNumberOfSingleHitTracks; 
  TH1D *hTofpNumberOfSingleValidHitTracks;

  // track matching plots
  TOrdCollection* mMatchHistNames;
  TH2D *hTofpMatchHit[NTOFP]; 
  TH2D *hTofpMatchNoHit[NTOFP]; 

  virtual const char *GetCVS() const 
    {static const char cvs[]="Tag $Name:  $ $Id: StTofpMatchMaker.h,v 1.9 2014/08/06 11:43:49 jeromel Exp $ built " __DATE__ " " __TIME__ ; return cvs;}

  ClassDef(StTofpMatchMaker,0)
};

inline void StTofpMatchMaker::setValidAdcRange(Int_t min, Int_t max){mMinValidAdc=min; mMaxValidAdc=max;}
inline void StTofpMatchMaker::setValidTdcRange(Int_t min, Int_t max){mMinValidTdc=min; mMaxValidTdc=max;}
inline void StTofpMatchMaker::setOuterTrackGeometry(){mOuterTrackGeometry=true;}
inline void StTofpMatchMaker::setStandardTrackGeometry(){mOuterTrackGeometry=false;}
inline void StTofpMatchMaker::setMinHitsPerTrack(Int_t nhits){mMinHitsPerTrack=nhits;}
inline void StTofpMatchMaker::setMinFitPointsPerTrack(Int_t nfitpnts){mMinFitPointsPerTrack=nfitpnts;}
inline void StTofpMatchMaker::setMaxDCA(Float_t maxdca){mMaxDCA=maxdca;}
inline void StTofpMatchMaker::setHistoFileName(const Char_t* filename){mHistoFileName=filename;}
inline void StTofpMatchMaker::createHistograms(Bool_t histos){mHisto = histos;}
inline bool StTofpMatchMaker::validAdc(float adc){return((adc>=mMinValidAdc) && (adc<=mMaxValidAdc));}
inline bool StTofpMatchMaker::validTdc(float tdc){return((tdc>=mMinValidTdc) && (tdc<=mMaxValidTdc));}

#endif