StVpdSimMaker.h

/***************************************************************************
 *
 * $Id: StVpdSimMaker.h,v 1.2 2021/01/27 04:33:30 geurts Exp $
 *
 * Author: Nickolas Luttrell (Rice University)
 ***************************************************************************
 *
 * Description: StVpdSimMaker.h   - Virtual base class for Vertex Position
 *  detector simulations
 *
 ***************************************************************************
 *
 * $Log: StVpdSimMaker.h,v $
 * Revision 1.2  2021/01/27 04:33:30  geurts
 * move StVpdSimConfig.h to StBTofUtil area for shared used between StBTofCalibMaker and StVpdSimMaker
 *
 * Revision 1.1  2017/03/02 18:38:17  jeromel
 * First version of Vpd simulations
 *
 *
 ***************************************************************************/

#ifndef StVpdSimMaker_HH
#define StVpdSimMaker_HH

#include "StMaker.h"
#include "St_DataSet.h"
class TH1F;
class TH2F;
class TH3F;
class TNtuple;
class StEvent;
class StBTofCollection;
class StVpdSimConfig;

#include "tables/St_g2t_vpd_hit_Table.h"
#include "StMcEvent/StMcEvent.hh"
#include "StMcEvent/StMcBTofHitCollection.hh"
#include "StBTofUtil/StVpdSimConfig.h"

#include <vector>
#ifndef ST_NO_NAMESPACES

#endif

class StVpdSimMaker : public StMaker {

public:
        StVpdSimMaker(const char *name = "VpdSim");
        virtual ~StVpdSimMaker();

        void           Reset();
        virtual int  Init();
        int          InitRun(int);
        int          FinishRun(int);
        virtual int  Make();
        virtual int  Finish();

    // Define get and set functions

        StBTofCollection*  GetVpdCollection()  const { return mVpdCollection; }
        StMcBTofHitCollection* GetMcBTofHitCollection() const { return mMcBTofHitCollection; }

        string   pullHistFileName();
    string   getParamsFileName() { return mParamsFileName; }
    void    setParamsFile(string fileName = "db/vpdSimParams/vpdSimParams.dat") { mParamsFileName = fileName; }
    void    setBookHisto(bool bookHist) { mBookHisto = bookHist; }

        virtual const char *GetCVS() const{
                static const char cvs[] = "Tag $Name:  $ $Id: StVpdSimMaker.h,v 1.2 2021/01/27 04:33:30 geurts Exp $ built " __DATE__ " " __TIME__; return cvs;
        }

protected:


    StMcBTofHitCollection * mMcBTofHitCollection = nullptr;     
    St_DataSet            * mGeantData           = nullptr;     
    StEvent               * mEvent               = nullptr;     
    StMcEvent             * mMcEvent             = nullptr;     
    StBTofCollection      * mVpdCollection      = nullptr;      
    std::map<int, StVpdSimConfig::SingleTubeParams>     mSimParams; 

        struct VpdSingleHit {
        int tray;       
        int tubeId;     
        double tof;     
        double t0;      
        double de;      
        double pathL;   
        double q;       
        };

        //  Various general use variables

    StVpdSimConfig* mSimConfig; 
        bool mBookHisto;            
    bool mUseFileParameters;    
    string mParamsFileName;     

        int       mNHits = 0;
        double    mVx = 0;
        double    mVy = 0;
        double    mVz = 0;
        double    mSumTubeTime = 0;      
        double    mTubeTAvg = 0;         

    double    mTStart = 0;           
        double    mTubeTAvgWest = 0;     
        double    mTubeTAvgEast = 0;     
        float     mVpdVertex = 0;        

        // Histogram variables
        string   mHistoFileName = "";    
        TH1F* mNRawHitsWest;    
    TH1F* mNRawHitsEast;    
        TH1F* mTubeHitsWest;    
        TH1F* mTubeHitsEast;    
        TH1F* mNHitsWest;       
        TH1F* mNHitsEast;       
        TH1F* mLeTimeWest;      
        TH1F* mLeTimeEast;      
        TH1F* mTStartHist;      
        TH2F* mLeTubeHitsWest;  
    TH2F* mLeTubeHitsEast;  
        TH1F* mZVertexHist;     
        TH1F* mVpdVertexHist;   
        TH1F* mVpdResHist;      
    TH3F* mResVsNumHits;    

                                                        // Functions

        int        vpdResponse(VpdSingleHit &Hit, g2t_vpd_hit_st* vpd_hit, int vId);
        double     thresholdCut(std::vector<VpdSingleHit> Hits, std::vector<int> Tube_Counts, TH1F* TubeHits, TH1F* NHits);
        int        storeMcVpdHit(StMcBTofHit* mcVpdHit);
        int        fillEvent();
        int        bookHistograms();

        ClassDef(StVpdSimMaker, 2)
};

#endif /* StVpdSimMaker_h */



// end of StVpdSimMaker.h