Reweight a MC survival probabiity when the particles through volumes. More...

#include <NucleonAbsorptionOutOfTargetReweighter.h>

Inheritance diagram for NeutrinoFluxReweight::NucleonAbsorptionOutOfTargetReweighter:

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

virtual	~NucleonAbsorptionOutOfTargetReweighter ()

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	inel_piAl_xsec

Detailed Description

Reweight a MC survival probabiity when the particles through volumes.

Definition at line 16 of file NucleonAbsorptionOutOfTargetReweighter.h.

Constructor & Destructor Documentation

NeutrinoFluxReweight::NucleonAbsorptionOutOfTargetReweighter::NucleonAbsorptionOutOfTargetReweighter	(	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 NucleonAbsorptionOutOfTargetReweighter.cpp.

                                                                                                                                                          :cvPars(cv_pars),univPars(univ_pars),iUniv(iuniv){ 
 
     // const boost::interprocess::flat_map<std::string, double>& dsig_table = univPars.getMap();
     inel_piAl_xsec = univPars.getParameterValue("inel_piAl_xsec");
 
   }

NeutrinoFluxReweight::NucleonAbsorptionOutOfTargetReweighter::~NucleonAbsorptionOutOfTargetReweighter ( )

virtual

Definition at line 18 of file NucleonAbsorptionOutOfTargetReweighter.cpp.

                                                                                  {
     
   }

Member Function Documentation

double NeutrinoFluxReweight::NucleonAbsorptionOutOfTargetReweighter::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 40 of file NucleonAbsorptionOutOfTargetReweighter.cpp.

                                                                                               {
     
     std::vector<ParticlesThroughVolumesData>  vec_ptv = aa.ptv_info;
     double shift = inel_piAl_xsec;
     
     double NA_mb    = 6.02E-4;
     double wgt      = 1.0;
     double tot_dist = 0.0;
     double low_val  = 1.E-20;   
     
     for(int index_vol =0;index_vol<3;index_vol++){
       for(int ii=0;ii<3;ii++){  
         tot_dist = vec_ptv[index_vol].AmountMat[ii];
         if(tot_dist<low_val)continue;
         if(vec_ptv[index_vol].Pdgs[ii]!=2212 && vec_ptv[index_vol].Pdgs[ii]!=2112) continue;    
         tot_dist *= NA_mb;
         tot_dist *= shift;
         wgt *= exp(-1.0*tot_dist);   
       }
     }
     
     return wgt;
     
   }

std::vector< bool > NeutrinoFluxReweight::NucleonAbsorptionOutOfTargetReweighter::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 21 of file NucleonAbsorptionOutOfTargetReweighter.cpp.

                                                                                                    {
  
     std::vector<bool> this_nodes;
     //int index_vol = 0;
     double low_val = 1.E-20;    
   
     std::vector<ParticlesThroughVolumesData>  vec_ptv = aa.ptv_info;
 
     bool passVOL = false;
     //Cheking at least one ancestor with amount of materail value:
     for(int index_vol =0;index_vol<3;index_vol++){
       for(int ii=0;ii<3;ii++){
         passVOL = passVOL || (vec_ptv[index_vol].AmountMat[ii] >low_val && (vec_ptv[index_vol].Pdgs[ii]==2212 || vec_ptv[index_vol].Pdgs[ii]==2112));
       }
     }
     this_nodes.push_back(passVOL);
     return this_nodes;
     
   }