PhokharaInputModule.h

 
/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2014  Belle II Collaboration                              *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Torben Ferber                                            *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#ifndef PHOKHARAINPUTMODULE_H
#define PHOKHARAINPUTMODULE_H
 
#include <framework/core/Module.h>
#include <framework/logging/Logger.h>
 
#include <generators/phokhara/Phokhara.h>
#include <generators/utilities/InitialParticleGeneration.h>
 
#include <mdst/dataobjects/MCParticleGraph.h>
 
#include <vector>
 
namespace Belle2 {
  class PhokharaInputModule : public Module {
 
  public:
 
    PhokharaInputModule();
 
    virtual ~PhokharaInputModule();
 
    virtual void initialize() override;
 
    virtual void event() override;
 
    virtual void terminate() override;
 
  protected:
 
    int m_finalState; 
    bool m_replaceMuonsByVirtualPhoton; 
    int m_nMaxTrials; 
    int m_nSearchMax; 
    int m_LO;  
    int m_NLO;  
    int m_QED;  
    int m_NLOIFI;  
    int m_alpha;  
    int m_pionff;  
    int m_pionstructure;  
    int m_kaonff;  
    int m_narres;  
    int m_protonff;  
    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; 
//     int m_finalState; /**< Final state: mu+mu-(0), pi+pi-(1), 2pi0pi+pi-(2), 2pi+2pi-(3), ppbar(4), nnbar(5), K+K-(6), K0K0bar(7), pi+pi-pi0(8), lamb(->pi-p)lambbar(->pi+pbar)(9), eta pi+ pi- (10) */
    double m_epsilon; 
    std::string m_ParameterFile; 
//     std::string m_InputFile; /**< file that holds all input parameters - NOT USED IN BELLE II IMPLEMENTATION!!!*/
    bool m_BeamEnergySpread; 
    Phokhara m_generator;   
    MCParticleGraph m_mcGraph; 
//     TFile* m_fileExtraInfo;
//     TH1D* m_th1dSDif;
 
    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
 
#endif /* PHOKHARAINPUTMODULE_H */