PhokharaInputModule.h

 
/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <framework/logging/Logger.h>
 
#include <generators/modules/GeneratorBaseModule.h>
#include <generators/phokhara/Phokhara.h>
#include <generators/utilities/InitialParticleGeneration.h>
 
#include <mdst/dataobjects/MCParticleGraph.h>
 
#include <vector>
 
namespace Belle2 {
  class PhokharaInputModule : public GeneratorBaseModule {
 
  public:
 
    PhokharaInputModule();
 
    virtual ~PhokharaInputModule();
 
    virtual void generatorInitialize() override;
 
    virtual void generatorEvent() override;
 
    virtual void terminate() override;
 
  protected:
 
    int m_finalState; 
    bool m_replaceMuonsByVirtualPhoton; 
    int m_nMaxTrials; 
    int m_nSearchMax; 
    int m_weighted;  
    int m_LO;  
    int m_NLO;  
    int m_fullNLO;  
    int m_QED;  
    int m_IFSNLO;  
    int m_alpha;  
    int m_pionff;  
    int m_pionstructure;  
    int m_kaonff;  
    int m_narres;  
    int m_protonff;  
    int m_chi_sw;  
    int m_be_r;  
    double m_cmsEnergy = -1.0; 
    std::vector<double> m_ScatteringAngleRangePhoton; 
    std::vector<double> m_ScatteringAngleRangeFinalStates; 
    double m_MinInvMassHadronsGamma; 
    double m_MinInvMassHadrons; 
    bool m_ForceMinInvMassHadronsCut; 
    double m_MaxInvMassHadrons; 
    double m_MinEnergyGamma; 
    double m_epsilon; 
    double m_beamres; 
    std::string m_ParameterFile; 
    bool m_BeamEnergySpread; 
    Phokhara m_generator;   
    MCParticleGraph m_mcGraph; 
    template <typename T>
    std::vector<T> make_vector(T const& t1, T const& t2);  
    template <typename T>
    std::pair<T, T> vectorToPair(std::vector<T>& vec, const std::string& name = "");
 
  private:
    void initializeGenerator();
 
    bool m_initialized{false}; 
    DBObjPtr<BeamParameters> m_beamParams; 
    InitialParticleGeneration m_initial; 
  }; 
  template <typename T>
  inline std::vector<T> PhokharaInputModule::make_vector(T const& t1, T const& t2)
  {
    std::vector<T> v;
    v.push_back(t1);
    v.push_back(t2);
    return v; 
  }
 
  template <typename T>
  inline std::pair<T, T> PhokharaInputModule::vectorToPair(std::vector<T>& vec, const std::string& name)
  {
    std::pair<T, T> p;
    if (vec.size() != 2) {
      B2ERROR("The parameter " + name + " has to have exactly 2 elements !");
    } else {
      p.first  = vec[0];
      p.second = vec[1];
    }
    return p;  
  }
 
} // end namespace Belle2