CRYHelper.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
/// \file CRYHelper.h
/// \brief Interface to the CRY cosmic ray generator
///
/// For documentation on CRY, see: http://nuclear.llnl.gov/simulation/
/// and http://nuclear.llnl.gov/simulations/additional_bsd.html
///
/// This class assumes that the user has a ROOT TGeoManager defined
/// 
/// \author  messier@indiana.edu
////////////////////////////////////////////////////////////////////////
#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 {
    /// Interface to the CRY cosmic-ray generator
  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;           ///< CRY configuration                              
    CRYGenerator*  fGen;             ///< The CRY generator                              
    double         fSampleTime;      ///< Amount of time to sample (seconds)             
    double         fToffset;         ///< Shift in time of particles (s)                 
    double         fEthresh;         ///< Cut on kinetic energy (GeV)            
    std::string    fWorldVolume;     ///< Name of the world volume                       
    std::string    fLatitude;        ///< Latitude of detector need space after value 
    std::string    fAltitude;        ///< Altitude of detector need space after value 
    std::string    fSubBoxL;         ///< Length of subbox (m) need space after value 
    double         fBoxDelta;        ///< Adjustment to the size of the world box in  
                                     ///< each dimension to avoid G4 rounding errors  
    bool           fSingleEventMode; ///< flag to turn on producing only a single cosmic ray
  };

  // 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
////////////////////////////////////////////////////////////////////////