ThinTargetBins.h

Go to the documentation of this file.

#ifndef THINTARGETBINS_H
#define THINTARGETBINS_H

#include <utility>
#include <map>
#include <vector>

namespace NeutrinoFluxReweight{
  
  /*! \class ThinTargetBins
   *  \brief A class to manage the bin definitions for MIPP Numi Yields
   */
  class ThinTargetBins{


  private:    
    ThinTargetBins();

  public:
    static ThinTargetBins* getInstance();
    
    //! Read a NA49 data pip xml file name to parse the bins
    void pC_pi_from_xml(const char* filename);

    //!Barton:
    void barton_pC_pi_from_xml(const char* filename);

    //! Read a NA49 data K xml file name to parse the bins
    void pC_k_from_xml(const char* filename);
    
    //! Read a NA49 data prt xml file name to parse the bins
    void pC_p_from_xml(const char* filename);
    
    //! Read a NA49 data neutron xml file name to parse the bins
    void pC_n_from_xml(const char* filename);

    //!MIPP k/pi:
    void mipp_pC_k_pi_from_xml(const char* filename);
    
    //! Read a pion incident 
    void meson_incident_from_xml(const char* filename);
    
    //! Read a pion incident 
    void material_scaling_from_xml(const char* filename);
    
    //! Return the Bin ID for this data
    int BinID_pC_pi(double xf, double pt,int pdgcode);
    
    //! Return the Bin ID for this data
    int barton_BinID_pC_pi(double xf, double pt,int pdgcode);
    
    //! Return the Bin ID for this data
    int BinID_pC_p(double xf, double pt,int pdgcode);

    //! Return the Bin ID for this data
    int BinID_pC_n(double xf,int pdgcode);

    //! Return the Bin ID for this data
    int BinID_pC_k(double xf, double pt,int pdgcode);
    
    //! Return the MIPP Thin Target Bin ID for this data
    int mipp_BinID_pC_k(double pz, double pt,int pdgcode);
    
    //! Return Pion incident bin
    int meson_inc_BinID(double xf, double pt,int pdgcode);
    
    //! Return material scaling bin
    int material_scaling_BinID(double xf, double pt,int pdgcode);

    //Pions production:
    std::vector<double> pC_pi_xfmin, pC_pi_xfmax, pC_pi_ptmin, pC_pi_ptmax;
    std::vector<double> b_pC_pi_xfmin, b_pC_pi_xfmax, b_pC_pi_ptmin, b_pC_pi_ptmax;
    
    //Proton production:
    std::vector<double> pC_p_xfmin, pC_p_xfmax, pC_p_ptmin, pC_p_ptmax;

    //Neutron production:
    std::vector<double> pC_n_xfmin, pC_n_xfmax;

    //Kaons production:
    std::vector<double> pC_k_xfmin, pC_k_xfmax, pC_k_ptmin, pC_k_ptmax;
    std::vector<double> mipp_pC_k_pzmin, mipp_pC_k_pzmax, mipp_pC_k_ptmin, mipp_pC_k_ptmax;
    
    //Pion incident
    std::vector<double> meson_inc_xfmin, meson_inc_xfmax, meson_inc_ptmin, meson_inc_ptmax;
    
    //Material scaling
    std::vector<double> mat_scal_xfmin, mat_scal_xfmax, mat_scal_ptmin, mat_scal_ptmax;

    //Number of bins:
    int GetNbins_pC_piX_NA49();
    int GetNbins_pC_piX_Barton();
    
    int GetNbins_pC_pX_NA49();
    int GetNbins_pC_nX_NA49();
    
    int GetNbins_pC_KX_NA49();
    int GetNbins_pC_KX_MIPP();
    
    int GetNbins_meson_incident();
    int GetNbins_material_scaling();

  private:
    static ThinTargetBins* instance;
    
  };

  
  
}
#endif