NeutrinoFluxReweight::MIPPNumiMC Class Reference

A class to manage the MC value for MIPP NuMI. More...

#include <MIPPNumiMC.h>

Public Member Functions
void	pip_mc_from_xml (const char *filename)
	Read a xml file name to get the mc value for pip. More...
void	pim_mc_from_xml (const char *filename)
	Read a xml file name to get the mc value for pim. More...
void	kap_mc_from_xml (const char *filename)
	Read a xml file name to get the mc value for kap. More...
void	kam_mc_from_xml (const char *filename)
	Read a xml file name to get the mc value for kam. More...
void	k0l_mc_from_xml (const char *filename)
	Read a xml file name to get the mc value for k0l. More...
void	k0s_mc_from_xml (const char *filename)
	Read a xml file name to get the mc value for k0s. More...
double	getMCval (double pz, double pt, int pdgcode)
	MC value for this HP production. More...

Static Public Member Functions
static MIPPNumiMC *	getInstance ()

Private Member Functions
	MIPPNumiMC ()

Private Attributes
std::vector< double >	pip_cv
std::vector< double >	pim_cv
std::vector< double >	kap_cv
std::vector< double >	kam_cv
std::vector< double >	k0l_cv
std::vector< double >	k0s_cv
std::vector< double >	v_pzmin
std::vector< double >	v_pzmax
std::vector< double >	v_ptmin
std::vector< double >	v_ptmax
bool	ranges_already_filled
double	proton_no_interacting

Static Private Attributes
static MIPPNumiMC *	instance = 0

Detailed Description

A class to manage the MC value for MIPP NuMI.

Definition at line 16 of file MIPPNumiMC.h.

Constructor & Destructor Documentation

NeutrinoFluxReweight::MIPPNumiMC::MIPPNumiMC ( )

private

Definition at line 10 of file MIPPNumiMC.cpp.

                        {
    ranges_already_filled = false;

    //FRaction of protons not interacting in the target or Budal Monitor for
    //LE NuMI mode using FTFP.
    proton_no_interacting = 0.13288294;
  }

Member Function Documentation

MIPPNumiMC * NeutrinoFluxReweight::MIPPNumiMC::getInstance ( )

static

Definition at line 295 of file MIPPNumiMC.cpp.

                                     {
    if (instance == 0) instance = new MIPPNumiMC;
    return instance;
  }

double NeutrinoFluxReweight::MIPPNumiMC::getMCval	(	double	pz,
		double	pt,
		int	pdgcode
	)

MC value for this HP production.

Definition at line 222 of file MIPPNumiMC.cpp.

                                                             {
 
    double cvmc = -1;
    if(abs(pdgcode)!=211 && abs(pdgcode)!=321 && pdgcode!=130 && pdgcode!=310)return cvmc;
    int size = 0;
   
    //pip:
    if(pdgcode==211){
      size = pip_cv.size();
      for(int ii=0;ii<size;++ii){
        if(pz>v_pzmin[ii] && pz<v_pzmax[ii] && pt>v_ptmin[ii] && pt<v_ptmax[ii]){
          cvmc = pip_cv[ii];
        }
      }      
    }

    //pim:
    if(pdgcode==-211){
       size = pim_cv.size();
      for(int ii=0;ii<size;++ii){
        if(pz>v_pzmin[ii] && pz<v_pzmax[ii] && pt>v_ptmin[ii] && pt<v_ptmax[ii]){
          cvmc = pim_cv[ii];
        }
      }  
    }

    //kap
    if(pdgcode==321){
       size = kap_cv.size();
      for(int ii=0;ii<size;++ii){
        if(pz>v_pzmin[ii] && pz<v_pzmax[ii] && pt>v_ptmin[ii] && pt<v_ptmax[ii]){
          cvmc = kap_cv[ii];
        }
      }  
    }

    //kam:
    if(pdgcode==-321){    
       size = kam_cv.size();
       for(int ii=0;ii<size;++ii){
        if(pz>v_pzmin[ii] && pz<v_pzmax[ii] && pt>v_ptmin[ii] && pt<v_ptmax[ii]){
          cvmc = kam_cv[ii];
        }
      }  
    }
    //k0l:
    if(pdgcode== 130){    
       size = k0l_cv.size();
       for(int ii=0;ii<size;++ii){
        if(pz>v_pzmin[ii] && pz<v_pzmax[ii] && pt>v_ptmin[ii] && pt<v_ptmax[ii]){
          cvmc = k0l_cv[ii];
        }
      }  
    }
    
    //k0s:
    if(pdgcode== 310){    
       size = k0s_cv.size();
       for(int ii=0;ii<size;++ii){
        if(pz>v_pzmin[ii] && pz<v_pzmax[ii] && pt>v_ptmin[ii] && pt<v_ptmax[ii]){
          cvmc = k0s_cv[ii];
        }
      }  
    }
    
    //The  values store in the files correspond to the Number of hadron per POT. 
    // But we are going to trasform to Number of hadron per interaction. 
    
   cvmc /=  (1.0-proton_no_interacting);
   return cvmc;
    
  }

void NeutrinoFluxReweight::MIPPNumiMC::k0l_mc_from_xml

(

const char *

filename

)

Read a xml file name to get the mc value for k0l.

Definition at line 156 of file MIPPNumiMC.cpp.

                                                      {
    using boost::property_tree::ptree;
    ptree top;
    read_xml(filename,top,2); 
    ptree bins;
    bins = top.get_child("mcbin.MIPPNuMI_MC_k0l");
    ptree::iterator it;    
    // we know that pip, pim, kap, kam, k0l and k0s have the same binning
    double cv,pzmin,pzmax,ptmin,ptmax;
    
    for(it = bins.begin(); it!=bins.end(); ++it){
   
      std::string cv_string=it->second.get<std::string>("cvmc");
      std::string pz_string=it->second.get<std::string>("pzrange");
      std::string pt_string=it->second.get<std::string>("ptrange");
      std::stringstream ss1(cv_string);
      std::stringstream ss2(pz_string);
      std::stringstream ss3(pt_string);
      ss1 >> cv;
      ss2 >> pzmin >> pzmax;
      ss3 >> ptmin >> ptmax;
   
      k0l_cv.push_back(cv);
      if(ranges_already_filled==false){
        v_pzmin.push_back(pzmin);
        v_pzmax.push_back(pzmax);
        v_ptmin.push_back(ptmin);
        v_ptmax.push_back(ptmax);
      }
    }
    ranges_already_filled=true;
  }

void NeutrinoFluxReweight::MIPPNumiMC::k0s_mc_from_xml

(

const char *

filename

)

Read a xml file name to get the mc value for k0s.

Definition at line 189 of file MIPPNumiMC.cpp.

                                                      {
    using boost::property_tree::ptree;
    ptree top;
    read_xml(filename,top,2); 
    ptree bins;
    bins = top.get_child("mcbin.MIPPNuMI_MC_k0s");
    ptree::iterator it;    
    // we know that pip, pim, kap, kam, k0l and k0s have the same binning
    double cv,pzmin,pzmax,ptmin,ptmax;
    
    for(it = bins.begin(); it!=bins.end(); ++it){
   
      std::string cv_string=it->second.get<std::string>("cvmc");
      std::string pz_string=it->second.get<std::string>("pzrange");
      std::string pt_string=it->second.get<std::string>("ptrange");
      std::stringstream ss1(cv_string);
      std::stringstream ss2(pz_string);
      std::stringstream ss3(pt_string);
      ss1 >> cv;
      ss2 >> pzmin >> pzmax;
      ss3 >> ptmin >> ptmax;
   
      k0s_cv.push_back(cv);
      if(ranges_already_filled==false){
        v_pzmin.push_back(pzmin);
        v_pzmax.push_back(pzmax);
        v_ptmin.push_back(ptmin);
        v_ptmax.push_back(ptmax);
      }
    }
    ranges_already_filled=true;
  }

void NeutrinoFluxReweight::MIPPNumiMC::kam_mc_from_xml

(

const char *

filename

)

Read a xml file name to get the mc value for kam.

Definition at line 123 of file MIPPNumiMC.cpp.

                                                    {
    using boost::property_tree::ptree;
    ptree top;
    read_xml(filename,top,2); 
    ptree bins;
    bins = top.get_child("mcbin.MIPPNuMI_MC_kam");
    ptree::iterator it;    
    // we know that pip, pim, kap and kam have the same binning
    double cv,pzmin,pzmax,ptmin,ptmax;
    
    for(it = bins.begin(); it!=bins.end(); ++it){
   
      std::string cv_string=it->second.get<std::string>("cvmc");
      std::string pz_string=it->second.get<std::string>("pzrange");
      std::string pt_string=it->second.get<std::string>("ptrange");
      std::stringstream ss1(cv_string);
      std::stringstream ss2(pz_string);
      std::stringstream ss3(pt_string);
      ss1 >> cv;
      ss2 >> pzmin >> pzmax;
      ss3 >> ptmin >> ptmax;
   
      kam_cv.push_back(cv);
      if(ranges_already_filled==false){
        v_pzmin.push_back(pzmin);
        v_pzmax.push_back(pzmax);
        v_ptmin.push_back(ptmin);
        v_ptmax.push_back(ptmax);
      }
    }
    ranges_already_filled=true;
  }

void NeutrinoFluxReweight::MIPPNumiMC::kap_mc_from_xml

(

const char *

filename

)

Read a xml file name to get the mc value for kap.

Definition at line 89 of file MIPPNumiMC.cpp.

                                                    {
    using boost::property_tree::ptree;
    ptree top;
    read_xml(filename,top,2); 
    ptree bins;
    bins = top.get_child("mcbin.MIPPNuMI_MC_kap");
    ptree::iterator it;    
    // we know that pip, pim, kap and kam have the same binning
    double cv,pzmin,pzmax,ptmin,ptmax;
    
    for(it = bins.begin(); it!=bins.end(); ++it){
   
      std::string cv_string=it->second.get<std::string>("cvmc");
      std::string pz_string=it->second.get<std::string>("pzrange");
      std::string pt_string=it->second.get<std::string>("ptrange");
      
      std::stringstream ss1(cv_string);
      std::stringstream ss2(pz_string);
      std::stringstream ss3(pt_string);
      ss1 >> cv;
      ss2 >> pzmin >> pzmax;
      ss3 >> ptmin >> ptmax;
   
      kap_cv.push_back(cv);
      if(ranges_already_filled==false){
        v_pzmin.push_back(pzmin);
        v_pzmax.push_back(pzmax);
        v_ptmin.push_back(ptmin);
        v_ptmax.push_back(ptmax);
      }
    }
    ranges_already_filled=true;
  }

void NeutrinoFluxReweight::MIPPNumiMC::pim_mc_from_xml

(

const char *

filename

)

Read a xml file name to get the mc value for pim.

Definition at line 54 of file MIPPNumiMC.cpp.

                                                    {
    using boost::property_tree::ptree;
    ptree top;
    read_xml(filename,top,2); 
    ptree bins;
    bins = top.get_child("mcbin.MIPPNuMI_MC_pim");
    ptree::iterator it;    
    // we know that pip, pim, kap and kam have the same binning
    double cv,pzmin,pzmax,ptmin,ptmax;
    
    for(it = bins.begin(); it!=bins.end(); ++it){
   
      std::string cv_string=it->second.get<std::string>("cvmc");
      std::string pz_string=it->second.get<std::string>("pzrange");
      std::string pt_string=it->second.get<std::string>("ptrange");

      std::stringstream ss1(cv_string);
      std::stringstream ss2(pz_string);
      std::stringstream ss3(pt_string);
      ss1 >> cv;
      ss2 >> pzmin >> pzmax;
      ss3 >> ptmin >> ptmax;
  
      pim_cv.push_back(cv);
      if(ranges_already_filled==false){
        v_pzmin.push_back(pzmin);
        v_pzmax.push_back(pzmax);
        v_ptmin.push_back(ptmin);
        v_ptmax.push_back(ptmax);
      }
    }
    ranges_already_filled=true;
}

void NeutrinoFluxReweight::MIPPNumiMC::pip_mc_from_xml

(

const char *

filename

)

Read a xml file name to get the mc value for pip.

Definition at line 18 of file MIPPNumiMC.cpp.

                                                      {
    using boost::property_tree::ptree;
    ptree top;
    read_xml(filename,top,2); 
    ptree bins;
    bins = top.get_child("mcbin.MIPPNuMI_MC_pip");
    ptree::iterator it;    
    // we know that pip, pim, kap and kam have the same binning
    double cv,pzmin,pzmax,ptmin,ptmax;
    
    for(it = bins.begin(); it!=bins.end(); ++it){
   
      std::string cv_string=it->second.get<std::string>("cvmc");
      std::string pz_string=it->second.get<std::string>("pzrange");
      std::string pt_string=it->second.get<std::string>("ptrange");

      std::stringstream ss1(cv_string);
      std::stringstream ss2(pz_string);
      std::stringstream ss3(pt_string);
      ss1 >> cv;
      ss2 >> pzmin >> pzmax;
      ss3 >> ptmin >> ptmax;
   
      pip_cv.push_back(cv);
      if(ranges_already_filled==false){
        v_pzmin.push_back(pzmin);
        v_pzmax.push_back(pzmax);
        v_ptmin.push_back(ptmin);
        v_ptmax.push_back(ptmax);
      }
      
    }
    ranges_already_filled=true;
    
  }

Member Data Documentation

MIPPNumiMC * NeutrinoFluxReweight::MIPPNumiMC::instance = 0

staticprivate

Definition at line 50 of file MIPPNumiMC.h.

std::vector<double> NeutrinoFluxReweight::MIPPNumiMC::k0l_cv

private

Definition at line 46 of file MIPPNumiMC.h.

std::vector<double> NeutrinoFluxReweight::MIPPNumiMC::k0s_cv

private

Definition at line 46 of file MIPPNumiMC.h.

std::vector<double> NeutrinoFluxReweight::MIPPNumiMC::kam_cv

private

Definition at line 46 of file MIPPNumiMC.h.

std::vector<double> NeutrinoFluxReweight::MIPPNumiMC::kap_cv

private

Definition at line 46 of file MIPPNumiMC.h.

std::vector<double> NeutrinoFluxReweight::MIPPNumiMC::pim_cv

private

Definition at line 46 of file MIPPNumiMC.h.

std::vector<double> NeutrinoFluxReweight::MIPPNumiMC::pip_cv

private

Definition at line 46 of file MIPPNumiMC.h.

double NeutrinoFluxReweight::MIPPNumiMC::proton_no_interacting

private

Definition at line 49 of file MIPPNumiMC.h.

bool NeutrinoFluxReweight::MIPPNumiMC::ranges_already_filled

private

Definition at line 48 of file MIPPNumiMC.h.

std::vector<double> NeutrinoFluxReweight::MIPPNumiMC::v_ptmax

private

Definition at line 47 of file MIPPNumiMC.h.

std::vector<double> NeutrinoFluxReweight::MIPPNumiMC::v_ptmin

private

Definition at line 47 of file MIPPNumiMC.h.

std::vector<double> NeutrinoFluxReweight::MIPPNumiMC::v_pzmax

private

Definition at line 47 of file MIPPNumiMC.h.

std::vector<double> NeutrinoFluxReweight::MIPPNumiMC::v_pzmin

private

Definition at line 47 of file MIPPNumiMC.h.

---

The documentation for this class was generated from the following files:

• ppfx/include/MIPPNumiMC.h
• ppfx/src/MIPPNumiMC.cpp