Reweight a chain of interactions that are covered by the NuMI target pi+ and pi- yields measured by MIPP. More...

#include <MIPPNumiPionYieldsReweighter.h>

Inheritance diagram for NeutrinoFluxReweight::MIPPNumiPionYieldsReweighter:

Public Member Functions
	MIPPNumiPionYieldsReweighter (int iuniv, const ParameterTable &cv_pars, const ParameterTable &univ_pars)

virtual	~MIPPNumiPionYieldsReweighter ()

virtual std::vector< bool >	canReweight (const InteractionChainData &aa)
	Look through the InteractionChainData input and identify those Interactions that can be reweighted as part of a chain. We return a vector indicating which elements will be assigned a weight by calculateWeight. More...

virtual double	calculateWeight (const InteractionChainData &aa)
	calculate a weight for this interaction chain given the central value parameters and the parameters for this universe. The weight is something like: f(cv)/f(MC) * f(univ)/f(cv) where cv in this case corresponds to the best value of the parameter, given the data. If univ_pars=cv_pars then we are calculating a central value weight. Note, canReweight() should be called to determine which elements of the chain are covered by the weight returned by calculateWeight() More...

Public Member Functions inherited from NeutrinoFluxReweight::IInteractionChainReweighting
virtual	~IInteractionChainReweighting ()

Public Attributes
const ParameterTable &	cvPars

const ParameterTable &	univPars

Private Attributes
int	iUniv

float	prt_no_inter

std::vector< float >	vbin_data_pip

std::vector< float >	vbin_data_pim

Detailed Description

Reweight a chain of interactions that are covered by the NuMI target pi+ and pi- yields measured by MIPP.

Definition at line 12 of file MIPPNumiPionYieldsReweighter.h.

Constructor & Destructor Documentation

NeutrinoFluxReweight::MIPPNumiPionYieldsReweighter::MIPPNumiPionYieldsReweighter	(	int	iuniv,
		const ParameterTable &	cv_pars,
		const ParameterTable &	univ_pars
	)

The constructor. Note, we pass central value and single universe parameters in this constructor only. There is thus a 1 to 1 correspondence between an instance of this class and a given universe.

Definition at line 12 of file MIPPNumiPionYieldsReweighter.cpp.

                                                                                                                                      :cvPars(cv_pars),univPars(univ_pars),iUniv(iuniv){ 
 
      MIPPNumiYieldsBins* MIPPbins =  MIPPNumiYieldsBins::getInstance();
      vbin_data_pip.reserve(MIPPbins->GetNbins_pip_MIPPNuMI());
      vbin_data_pim.reserve(MIPPbins->GetNbins_pim_MIPPNuMI());
      
      //const boost::interprocess::flat_map<std::string, double>& cv_table = cvPars.getMap();
      //const boost::interprocess::flat_map<std::string, double>& univ_table = univPars.getMap();
     prt_no_inter = univPars.getParameterValue("prt_no_interacting");
     char namepar[100];
     for(int ii=0;ii<MIPPbins->GetNbins_pip_MIPPNuMI();ii++){
       sprintf(namepar,"MIPP_NuMI_%s_sys_%d","pip",ii);
       double data_cv  = cvPars.getParameterValue(std::string(namepar));
       double data_sys = univPars.getParameterValue(std::string(namepar));
       sprintf(namepar,"MIPP_NuMI_%s_stats_%d","pip",ii);
       double data_sta = univPars.getParameterValue(std::string(namepar));
       data_sys /= (1.0-prt_no_inter);
       data_sta /= (1.0-prt_no_inter);
       data_cv  /= (1.0-prt_no_inter);
       vbin_data_pip.push_back(data_sys + data_sta - data_cv);
     }
     for(int ii=0;ii<MIPPbins->GetNbins_pim_MIPPNuMI();ii++){
       sprintf(namepar,"MIPP_NuMI_%s_sys_%d","pim",ii);
       double data_cv  = cvPars.getParameterValue(std::string(namepar));
       double data_sys = univPars.getParameterValue(std::string(namepar));
       sprintf(namepar,"MIPP_NuMI_%s_stats_%d","pim",ii);
       double data_sta = univPars.getParameterValue(std::string(namepar));
       data_sys /= (1.0-prt_no_inter);
       data_sta /= (1.0-prt_no_inter);
       data_cv  /= (1.0-prt_no_inter);
       vbin_data_pim.push_back(data_sys + data_sta - data_cv);
     }
 
     
 }

NeutrinoFluxReweight::MIPPNumiPionYieldsReweighter::~MIPPNumiPionYieldsReweighter ( )

virtual

Definition at line 47 of file MIPPNumiPionYieldsReweighter.cpp.

                                                              {
     
   }

Member Function Documentation

double NeutrinoFluxReweight::MIPPNumiPionYieldsReweighter::calculateWeight ( const InteractionChainData & aa )

virtual

calculate a weight for this interaction chain given the central value parameters and the parameters for this universe. The weight is something like: f(cv)/f(MC) * f(univ)/f(cv) where cv in this case corresponds to the best value of the parameter, given the data. If univ_pars=cv_pars then we are calculating a central value weight. Note, canReweight() should be called to determine which elements of the chain are covered by the weight returned by calculateWeight()

Implements NeutrinoFluxReweight::IInteractionChainReweighting.

Definition at line 86 of file MIPPNumiPionYieldsReweighter.cpp.

                                                                                     {
     
     double wgt = 1.0;
     
     MIPPNumiYieldsBins*  MIPPbins =  MIPPNumiYieldsBins::getInstance();
     MIPPNumiMC*  MCval =  MIPPNumiMC::getInstance();    
     
     TargetData tar = aa.tar_info;
     int binID = MIPPbins->BinID(tar.Pz,tar.Pt,tar.Tar_pdg);
     
     if(binID<0){
       std::cout<<"BINID ZERO:" <<tar.Pz<<" "<<tar.Pt<<" " <<tar.Tar_pdg<<std::endl;
       return 1.0;
     }
     
     //Getting MC:
     double binC = MCval->getMCval(tar.Pz,tar.Pt,tar.Tar_pdg);
     if(binC<1.e-18){
       std::cout<<"LOW MC VAL: "<<binC <<std::endl;
       return 1.0;
     }
     
     if(tar.Tar_pdg == 211)wgt = vbin_data_pip[binID]/binC;
     if(tar.Tar_pdg ==-211)wgt = vbin_data_pim[binID]/binC;
     
     if(wgt<0)std::cout<<"TTMIPPPI check wgt(<0) "<<iUniv<<" "<<tar.Pz<<" "<<tar.Pt<<" "<<tar.Tar_pdg<<std::endl;
     return wgt;
     
   }

std::vector< bool > NeutrinoFluxReweight::MIPPNumiPionYieldsReweighter::canReweight ( const InteractionChainData & aa )

virtual

Look through the InteractionChainData input and identify those Interactions that can be reweighted as part of a chain. We return a vector indicating which elements will be assigned a weight by calculateWeight.

Implements NeutrinoFluxReweight::IInteractionChainReweighting.

Definition at line 50 of file MIPPNumiPionYieldsReweighter.cpp.

                                                                                          {
  
     MIPPNumiYieldsBins*  MIPPbins =  MIPPNumiYieldsBins::getInstance();
     std::vector<bool> this_nodes;
     for(size_t ii=0;ii<(aa.interaction_chain).size();ii++){
       this_nodes.push_back(false);
     }
    
     //Look for MIPP Numi events
     //if the code find a MIPP Numi event, it will look 
     //for how many interaction nodes covers
     //if not, return all nodes false.
     //bool is_there_mipp = false;   
     TargetData tar = aa.tar_info;
     //Cheking if the particle is a pion plus and pion minus:
     if(tar.Tar_pdg != 211 && tar.Tar_pdg != -211)return this_nodes;
     int binID = MIPPbins->BinID(tar.Pz,tar.Pt,tar.Tar_pdg);
     if(binID<0) return this_nodes;    
     
     //Now that we know that we have a MIPP Numi event, 
     //we will see how many nodes are covered.
     std::vector<InteractionData> this_interactions = aa.interaction_chain; 
     
     //Now we have the index of the hadron that exit the target in the 
     //ancesty chain:
     if(tar.Idx_ancestry>=0){
       for(int ii=0;ii<tar.Idx_ancestry;ii++){
         this_nodes[ii] = true;
       }
     }
     else{
       //  std::cout<<"==>>SOMETHING WEIRD WITH MIPP NUMI "<<tar.Idx_ancestry<<" "<<tar.Tar_pdg<<std::endl;
     }
     
     return this_nodes;
   }