StMtdSimMaker.h

/***************************************************************************
 *
 * $Id: StMtdSimMaker.h,v 1.11 2015/07/29 01:11:24 smirnovd Exp $
 *
 * Author:  Frank Geurts
 ***************************************************************************
 *
 * Description: StMtdSimMaker virtual base class for Barrel TOF Simulations
 *
 ***************************************************************************
 *
 * $Log: StMtdSimMaker.h,v $
 * Revision 1.11  2015/07/29 01:11:24  smirnovd
 * Initialize static constants outside of class definition
 *
 * C++ forbids initialization of non-integral static const members within the class
 * definition. The syntax is allowed only for integral type variables.
 *
 * Revision 1.10  2015/07/01 17:51:54  marr
 * Add data member mWriteHisto to control the output of histograms
 *
 * Revision 1.9  2014/12/10 17:29:19  marr
 * 1. Use localz and time-of-flight information from GEANT to determine the
 * leading times on east and west sides for MC MTD hits. With this change, the
 * localz of the MC hits can be calcualted exactly the same as the regular
 * hits.
 * 2. Cell Id runs from 0 to 11 as in real data
 * 3. Clean up lines that are commented out.
 *
 * Revision 1.8  2014/08/06 11:43:27  jeromel
 * Suffix on literals need to be space (later gcc compiler makes it an error) - first wave of fixes
 *
 * Revision 1.7  2014/07/16 20:09:11  marr
 * Move the initialization of the GEANT-to-Backleg map using the database to InitRun()
 *
 * Revision 1.5  2014/01/20 18:01:05  geurts
 * bug update: changed default Maker name from TofSim to MtdSim in order to prevent name clashes with StBTofSimMaker [Bill LLope]
 *
 * Revision 1.4  2013/11/14 16:17:08  geurts
 * Correct mapping based on GEANT volume_id [Alex Jentsch]
 *
 * Revision 1.3  2013/07/05 21:57:34  geurts
 * Bug fix and improved mapping [Alex Jentsch]
 *
 * Revision 1.2  2012/03/02 02:18:34  perev
 * New provisional version rewritten by VP
 *
 *
 **************************************************************************/
#ifndef STBMTDSIMMAKER_HH
#define STBMTDSIMMAKER_HH
#include "StMaker.h"

class StEvent;
struct g2t_mtd_hit_st;
class StMtdCollection;
class TH2F;
class TH2I;

// g2t tables
#include "tables/St_g2t_mtd_hit_Table.h"
#include "tables/St_g2t_track_Table.h"
#include "tables/St_g2t_tpc_hit_Table.h"

#include <vector>

class StMtdSimMaker : public StMaker
{
 private:
  Int_t mModuleChannel[5][24];
  

 protected:
  St_DataSet        *mGeantData;        
  StEvent           *mEvent;            
  StMtdCollection   *mMtdCollection;   
  int                mNMtdHits;
  g2t_mtd_hit_st    *mMtdHitsFromGeant;

  //define some constants
  enum {
   kNBackleg = 99,
   kNModule  = 99,
   kNCell    = 12,     
   kAMP      = 50000,     
   kADCBINWIDTH = 25,
   kTDCBINWIDTH = 50
 };

  static const float kMtdPadWidth;      
  Bool_t mBookHisto;
  Bool_t mWriteHisto;
  Bool_t mWriteStEvent;  



  string mHistFile;    //for QA histograms
  TH1F* mBetaHist;    
  TH1F* mPathLHist;    
  TH1F* mTofHist;    
  TH1F* mRecMass;    

  TH2F* mCellGeant;    
  TH1F* mVpdGeant;     
  TH2F* mNCellGeant;   
  TH2F* mNVpdGeant;    
  TH1F* mDeGeant;      
  TH1F* mTofGeant;     

  TH2F* mCellSeen;     
  TH1F* mVpdSeen;      
  TH2F* mNCellSeen;    
  TH2F* mNVpdSeen;     
  TH1F* mDeSeen;       
  TH1F* mT0Seen;       
  TH1F* mTofSeen;      
  TH1F* mTofResSeen;   
  TH1F* mVpdResSeen;   

  TH2F* mCellReco;     
  TH1F* mVpdReco;      
  TH2F* mNCellReco;    
  TH2F* mNVpdReco;     
  TH1F* mTDCReco;      
  TH1F* mADCReco;      
  TH1F* mT0Reco;       
  TH1F* mTofResReco;   
  TH1F* mVpdResReco;   
  TH2F* mTACorr;       
  TH1F* mModHist;       
  TH2I* QABacklegChannel; 

    TH1F* mdE;           
    TH1F* mdS;           
    TH1F* mNumberOfPhotoelectrons;  
    TH1F* mT;            
    TH1F* mTime;         
    TH1F* mTime1;        
    TH1F* mPMlength;     
    TH1F* mAdc;          
    TH1F* mTdc;          

    TVolume *starHall;

    Int_t        bookHistograms();
    Int_t        writeHistograms();


 public:
    StMtdSimMaker(const char *name="MtdSim");
    virtual ~StMtdSimMaker();

    void           Reset();
    virtual Int_t  Init();
    Int_t          InitRun(Int_t);
    Int_t          FinishRun(Int_t);
    virtual Int_t  Make();
    virtual Int_t  Finish();
    Int_t          FastCellResponse();
    Int_t          CalcCellId(Int_t volume_id, Float_t ylocal,
                              Int_t &ibackleg,Int_t &imodule,Int_t &icell);

    virtual const char *GetCVS() const
    {static const char cvs[]="Tag $Name:  $ $Id: StMtdSimMaker.h,v 1.11 2015/07/29 01:11:24 smirnovd Exp $ built " __DATE__ " " __TIME__ ; return cvs;}

    ClassDef(StMtdSimMaker,2)
};
#endif