CRYHelper.h

Go to the documentation of this file.

#ifndef EVGB_CRYHELPER_H
#define EVGB_CRYHELPER_H
#include <string>
#include <vector>
#include "CLHEP/Random/RandEngine.h"

namespace simb { class MCTruth;  }

class CRYSetup;
class CRYGenerator;
class CRYParticle;

namespace evgb {
  class CRYHelper {
  public:
    CRYHelper();
    explicit CRYHelper(fhicl::ParameterSet     const& pset, 
                       CLHEP::HepRandomEngine&        engine,
                       std::string             const& worldVol="vWorld");
    ~CRYHelper();

    double Sample(simb::MCTruth& mctruth, 
                double       const& surfaceY,
                double       const& detectorLength,
                double*             w,
                double              rantime=0);
    
  private:

    void WorldBox(double* xlo_cm,
                  double* xhi_cm,
                  double* ylo_cm,
                  double* yhi_cm,
                  double* zlo_cm,
                  double* zhi_cm) const;

    void ProjectToBoxEdge(const double xyz[],
                          const double dxyz[],
                          double &xlo, double &xhi,
                          double &ylo, double &yhi,
                          double &zlo, double &zhi,
                          double xyzout[]);

    CRYSetup*      fSetup;           
    CRYGenerator*  fGen;             
    double         fSampleTime;      
    double         fToffset;         
    double         fEthresh;         
    std::string    fWorldVolume;     
    std::string    fLatitude;        
    std::string    fAltitude;        
    std::string    fSubBoxL;         
    double         fBoxDelta;        
    bool           fSingleEventMode; 
  };

  // The following stuff is for the random number gererator
  template<class T> class RNGWrapper {
  public:
    static void set(T* object, double (T::*func)(void));
    static double rng(void);
  private:
    static T* m_obj;
    static double (T::*m_func)(void);
  };// end of RNGWrapper class
  
  template<class T> T* RNGWrapper<T>::m_obj;
  
  template<class T> double (T::*RNGWrapper<T>::m_func)(void);
  
  template<class T> void RNGWrapper<T>::set(T* object, double (T::*func)(void)) {
    m_obj = object; m_func = func;
  }
  
  template<class T> double RNGWrapper<T>::rng(void) { return (m_obj->*m_func)(); }

}
#endif // EVGB_CRYHELPER_H