A class to manage and drive the weight calculation procedure. More...

#include <ReweightDriver.h>

Public Member Functions
	ReweightDriver (int iuniv, const ParameterTable &cv_pars, const ParameterTable &univ_pars, std::string fileIn)

	~ReweightDriver ()

double	calculateWeight (const InteractionChainData &icd)

Public Attributes
double	mipp_pion_wgt

double	mipp_kaon_wgt

double	att_wgt

double	abs_ic_wgt

double	abs_dpip_wgt

double	abs_dvol_wgt

double	abs_nucleon_wgt

double	abs_other_wgt

double	tot_abs_wgt

double	pC_pi_wgt

double	pC_k_wgt

double	nC_pi_wgt

double	pC_nu_wgt

double	nuA_wgt

double	meson_inc_wgt

double	other_wgt

MIPPNumiPionYieldsReweighter *	MIPP_NUMI_PION_Universe

MIPPNumiKaonYieldsReweighter *	MIPP_NUMI_KAON_Universe

TargetAttenuationReweighter *	TARG_ATT_Universe

AbsorptionICReweighter *	VOL_ABS_IC_Universe

AbsorptionDPIPReweighter *	VOL_ABS_DPIP_Universe

AbsorptionDVOLReweighter *	VOL_ABS_DVOL_Universe

NucleonAbsorptionOutOfTargetReweighter *	VOL_ABS_NUCLEON_Universe

OtherAbsorptionOutOfTargetReweighter *	VOL_ABS_OTHER_Universe

ThinTargetpCPionReweighter *	THINTARGET_PC_PION_Universe

ThinTargetpCKaonReweighter *	THINTARGET_PC_KAON_Universe

ThinTargetnCPionReweighter *	THINTARGET_NC_PION_Universe

ThinTargetpCNucleonReweighter *	THINTARGET_PC_NUCLEON_Universe

ThinTargetMesonIncidentReweighter *	THINTARGET_MESON_INCIDENT_Universe

ThinTargetnucleonAReweighter *	THINTARGET_NUCLEON_A_Universe

OtherReweighter *	OTHER_Universe

bool	doMIPPNumi

Private Member Functions
void	Configure ()

void	ParseOptions ()

Private Attributes
int	iUniv

const ParameterTable &	cvPars

const ParameterTable &	univPars

std::string	fileOptions

Detailed Description

A class to manage and drive the weight calculation procedure.

There is a one-to-one correspondence between ReweightDrivers and universes.

Definition at line 33 of file ReweightDriver.h.

Constructor & Destructor Documentation

NeutrinoFluxReweight::ReweightDriver::ReweightDriver	(	int	iuniv,
		const ParameterTable &	cv_pars,
		const ParameterTable &	univ_pars,
		std::string	fileIn
	)

the constructor

Definition at line 9 of file ReweightDriver.cpp.

     : iUniv(iuniv), cvPars(cv_pars), univPars(univ_pars), fileOptions(fileIn)
   {
     ParseOptions();
     Configure();
     
   }

NeutrinoFluxReweight::ReweightDriver::~ReweightDriver ( )

Definition at line 263 of file ReweightDriver.cpp.

                                  {
  
     if(doMIPPNumi){
       delete MIPP_NUMI_PION_Universe;
       delete MIPP_NUMI_KAON_Universe;
     }    
     delete TARG_ATT_Universe;
     delete VOL_ABS_IC_Universe;
     delete VOL_ABS_DPIP_Universe;
     delete VOL_ABS_DVOL_Universe;
     delete VOL_ABS_NUCLEON_Universe;
     delete VOL_ABS_OTHER_Universe;
     delete THINTARGET_PC_PION_Universe;
     delete THINTARGET_PC_KAON_Universe;
     delete THINTARGET_NC_PION_Universe;
     delete THINTARGET_PC_NUCLEON_Universe;
     delete THINTARGET_MESON_INCIDENT_Universe;
     delete THINTARGET_NUCLEON_A_Universe;
     delete OTHER_Universe;
   }

Member Function Documentation

double NeutrinoFluxReweight::ReweightDriver::calculateWeight ( const InteractionChainData & icd )

Calculate weight for this event, given a set of central value parameters and a set of randomly varied parameters corresponding to the universe that this ReweightDriver is responsible for.

Call mipp_yields on input chain. record interactions covered with a weight in interaction_covered.
Call attenuation on input chain
Call thin target reweighters.
1. Call na49 for interactions not yet covered. record in interaction_covered
2. ditto. call MIPP, then theory

--— PROCESS INTERACTION NODES --— ///

Definition at line 56 of file ReweightDriver.cpp.

                                                                        {
 
     double tot_wgt = 1.0;
     
     //Boolean flags: 
     const int nnodes=icd.interaction_chain.size();
     std::vector<bool> interaction_nodes(nnodes,false);
     std::vector<bool> attenuation_nodes(nnodes,false);
     std::vector<bool> absorption_nodes(nnodes,false);
 
     /// ----- PROCESS INTERACTION NODES ----- ///
     
     //MIPP NuMI Pions:
     bool has_mipp = false;
     mipp_pion_wgt = 1.0;
     if(doMIPPNumi){
       interaction_nodes = MIPP_NUMI_PION_Universe->canReweight(icd);
       for(size_t ii=0;ii<interaction_nodes.size();ii++){
         if(interaction_nodes[ii]==true){
           has_mipp = true;
           mipp_pion_wgt = MIPP_NUMI_PION_Universe->calculateWeight(icd);
           break; 
         }
       }
       tot_wgt *= mipp_pion_wgt;
     }
     
     //MIPP NuMI Kaons:
     mipp_kaon_wgt = 1.0;
     if(!has_mipp && doMIPPNumi){
       interaction_nodes = MIPP_NUMI_KAON_Universe->canReweight(icd);
 
       for(size_t ii=0;ii<interaction_nodes.size();ii++){
         if(interaction_nodes[ii]==true){
           has_mipp = true;
           mipp_kaon_wgt = MIPP_NUMI_KAON_Universe->calculateWeight(icd);
           break; 
         }
       }
       tot_wgt *= mipp_kaon_wgt;
     }
     
     //Thin Target pC->piX:
     pC_pi_wgt = 1.0;
     for(int ii=(interaction_nodes.size()-1);ii>=0;ii--){        
       if(interaction_nodes[ii]==false){
         bool is_rew = THINTARGET_PC_PION_Universe->canReweight((icd.interaction_chain)[ii]);
         if(is_rew){
           double rewval = THINTARGET_PC_PION_Universe->calculateWeight((icd.interaction_chain)[ii]);
           pC_pi_wgt *= rewval;
           interaction_nodes[ii]=true;
         }
       }
     }
     tot_wgt *= pC_pi_wgt;
        
     //Thin Target pC->KX:
     pC_k_wgt = 1.0;    
     for(int ii=(interaction_nodes.size()-1);ii>=0;ii--){        
       if(interaction_nodes[ii]==false){
         bool is_rew = THINTARGET_PC_KAON_Universe->canReweight((icd.interaction_chain)[ii]);
         if(is_rew){
           double rewval = THINTARGET_PC_KAON_Universe->calculateWeight((icd.interaction_chain)[ii]);
           pC_k_wgt *= rewval;
           interaction_nodes[ii]=true;
         }
       }
     }
     tot_wgt *= pC_k_wgt;
     
     //Thin Target nC->piX:
     nC_pi_wgt = 1.0;
     for(int ii=(interaction_nodes.size()-1);ii>=0;ii--){        
       if(interaction_nodes[ii]==false){
         bool is_rew = THINTARGET_NC_PION_Universe->canReweight((icd.interaction_chain)[ii]);
         if(is_rew){
           double rewval = THINTARGET_NC_PION_Universe->calculateWeight((icd.interaction_chain)[ii]);
           nC_pi_wgt *= rewval;
           interaction_nodes[ii]=true;
         }
       }
     }
     tot_wgt *= nC_pi_wgt;
     
     //Thin Target pC->nucleonX:
     pC_nu_wgt = 1.0;
     for(int ii=(interaction_nodes.size()-1);ii>=0;ii--){        
       if(interaction_nodes[ii]==false){
         bool is_rew = THINTARGET_PC_NUCLEON_Universe->canReweight((icd.interaction_chain)[ii]);
         if(is_rew){
           double rewval = THINTARGET_PC_NUCLEON_Universe->calculateWeight((icd.interaction_chain)[ii]);
           pC_nu_wgt *= rewval;
           interaction_nodes[ii]=true;
         }
       }
     }
     tot_wgt *= pC_nu_wgt;
     
     //Thin Target Meson Incident:
     meson_inc_wgt = 1.0; 
     for(int ii=(interaction_nodes.size()-1);ii>=0;ii--){        
       if(interaction_nodes[ii]==false){
         bool is_rew = THINTARGET_MESON_INCIDENT_Universe->canReweight((icd.interaction_chain)[ii]);
         if(is_rew){
           double rewval = THINTARGET_MESON_INCIDENT_Universe->calculateWeight((icd.interaction_chain)[ii]);
           meson_inc_wgt *= rewval;
           interaction_nodes[ii]=true;
         }
       }
     }
     tot_wgt *= meson_inc_wgt;    
     
     //Thin Target Nucleon Incident not hanldle NA49 or Barton:
     nuA_wgt = 1.0; 
     for(int ii=(interaction_nodes.size()-1);ii>=0;ii--){        
       if(interaction_nodes[ii]==false){
         bool is_rew = THINTARGET_NUCLEON_A_Universe->canReweight((icd.interaction_chain)[ii]);
         if(is_rew){
           double rewval = THINTARGET_NUCLEON_A_Universe->calculateWeight((icd.interaction_chain)[ii]);
           nuA_wgt *= rewval;
           interaction_nodes[ii]=true;
         }
       }
     }
     tot_wgt *= nuA_wgt;
     
     //Any other interaction not handled yet:
     other_wgt = 1.0;
     for(int ii=(interaction_nodes.size()-1);ii>=0;ii--){        
       if(interaction_nodes[ii]==false){
         bool is_rew = OTHER_Universe->canReweight((icd.interaction_chain)[ii]);
         if(is_rew){
           double rewval = OTHER_Universe->calculateWeight((icd.interaction_chain)[ii]);
           other_wgt *= rewval;
           interaction_nodes[ii]=true;
         }
       }
     }
     tot_wgt *= other_wgt;
     
     //Target attenuation correction:
     att_wgt = 1.0;
     attenuation_nodes = TARG_ATT_Universe->canReweight(icd);   
     //we just see for the first position (prmary proton)
     if(attenuation_nodes.size()>0 && attenuation_nodes[0]==true){
       att_wgt *= TARG_ATT_Universe->calculateWeight(icd);
     }
     tot_wgt *= att_wgt;
         
     //Absorption correction:
     tot_abs_wgt = 1.0;
     
     // Correction of the pi & K absorption in volumes (Al)    
     abs_ic_wgt = 1.0;    
     absorption_nodes = VOL_ABS_IC_Universe->canReweight(icd);
     //std::cout<<"size of absorption_nodes is "<<absorption_nodes.size()<<std::endl;
 
     if(absorption_nodes.size()>0 && absorption_nodes[0]==true){
       abs_ic_wgt *= VOL_ABS_IC_Universe->calculateWeight(icd);
     }
     tot_wgt     *= abs_ic_wgt;
     tot_abs_wgt *= abs_ic_wgt;
  
     //Correction of the pi & K absorption in volumes (Fe)
     abs_dpip_wgt = 1.0;    
     absorption_nodes = VOL_ABS_DPIP_Universe->canReweight(icd);
     if(absorption_nodes.size()>0 && absorption_nodes[0]==true){
       abs_dpip_wgt *= VOL_ABS_DPIP_Universe->calculateWeight(icd);
     }
     tot_wgt     *= abs_dpip_wgt;
     tot_abs_wgt *= abs_dpip_wgt;
     
     
     //Correction of the pi & K absorption in volumes (He)
     abs_dvol_wgt = 1.0;    
     absorption_nodes = VOL_ABS_DVOL_Universe->canReweight(icd);
     if(absorption_nodes.size()>0 && absorption_nodes[0]==true){
       abs_dvol_wgt *= VOL_ABS_DVOL_Universe->calculateWeight(icd);
     }
     tot_wgt     *= abs_dvol_wgt;
     tot_abs_wgt *= abs_dvol_wgt;
     
     //Correction of nucleons on Al, Fe and He.
     abs_nucleon_wgt = 1.0;    
     absorption_nodes = VOL_ABS_NUCLEON_Universe->canReweight(icd);
     if(absorption_nodes.size()>0 && absorption_nodes[0]==true){
       abs_nucleon_wgt *= VOL_ABS_NUCLEON_Universe->calculateWeight(icd);
     }
     tot_wgt     *= abs_nucleon_wgt;
     tot_abs_wgt *= abs_nucleon_wgt;
     
     //Correction of any other particle on Al, Fe and He.
     abs_other_wgt = 1.0;    
     absorption_nodes = VOL_ABS_OTHER_Universe->canReweight(icd);
     if(absorption_nodes.size()>0 && absorption_nodes[0]==true){
       abs_other_wgt *= VOL_ABS_OTHER_Universe->calculateWeight(icd);
     }
     tot_wgt     *= abs_other_wgt;
     tot_abs_wgt *= abs_other_wgt;
     
     if(tot_wgt!=tot_wgt){
       std::cout<<"Alert nan total wgt... check!!!"<<std::endl;
       return 1.0;
     }
     return tot_wgt;
   }

void NeutrinoFluxReweight::ReweightDriver::Configure ( void )

private

Configures each of the reweighing tools.

Definition at line 17 of file ReweightDriver.cpp.

                                 {
     
     //Creating the vector of reweighters:
     
     if(doMIPPNumi){
       MIPP_NUMI_PION_Universe = new MIPPNumiPionYieldsReweighter(iUniv,cvPars,univPars);
       MIPP_NUMI_KAON_Universe = new MIPPNumiKaonYieldsReweighter(iUniv,cvPars,univPars);
     }
     
     TARG_ATT_Universe = new TargetAttenuationReweighter(iUniv,cvPars,univPars);    
     VOL_ABS_IC_Universe = new AbsorptionICReweighter(iUniv,cvPars,univPars);
     VOL_ABS_DPIP_Universe = new AbsorptionDPIPReweighter(iUniv,cvPars,univPars);
     VOL_ABS_DVOL_Universe = new AbsorptionDVOLReweighter(iUniv,cvPars,univPars);
     VOL_ABS_NUCLEON_Universe = new NucleonAbsorptionOutOfTargetReweighter(iUniv,cvPars,univPars);
     VOL_ABS_OTHER_Universe = new OtherAbsorptionOutOfTargetReweighter(iUniv,cvPars,univPars);
 
     THINTARGET_PC_PION_Universe = new ThinTargetpCPionReweighter(iUniv,cvPars,univPars);
     THINTARGET_PC_KAON_Universe = new ThinTargetpCKaonReweighter(iUniv,cvPars,univPars);
     THINTARGET_NC_PION_Universe = new ThinTargetnCPionReweighter(iUniv,cvPars,univPars);
     THINTARGET_PC_NUCLEON_Universe = new ThinTargetpCNucleonReweighter(iUniv,cvPars,univPars);
     THINTARGET_MESON_INCIDENT_Universe = new ThinTargetMesonIncidentReweighter(iUniv,cvPars,univPars);
     THINTARGET_NUCLEON_A_Universe = new ThinTargetnucleonAReweighter(iUniv,cvPars,univPars);
     OTHER_Universe = new OtherReweighter(iUniv,cvPars,univPars);
     
   }

void NeutrinoFluxReweight::ReweightDriver::ParseOptions ( )

private

Parse the option to select which reweighter use

Definition at line 43 of file ReweightDriver.cpp.

                                    {
     //Parsing the file input:
     using boost::property_tree::ptree;
     ptree top;
     std::string val = "";
     read_xml(fileOptions.c_str(),top,2); // option 2 removes comment strings
     ptree& options = top.get_child("inputs.Settings");
     
     val = options.get<std::string>("Reweighters");
     if(val=="MIPPNuMIOn")doMIPPNumi = true;
     else doMIPPNumi = false;
     
   }

Member Data Documentation

double NeutrinoFluxReweight::ReweightDriver::abs_dpip_wgt

Definition at line 65 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::abs_dvol_wgt

Definition at line 66 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::abs_ic_wgt

Absorption weight

Definition at line 64 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::abs_nucleon_wgt

Definition at line 67 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::abs_other_wgt

Definition at line 68 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::att_wgt

Target attenuation weight

Definition at line 61 of file ReweightDriver.h.

const ParameterTable& NeutrinoFluxReweight::ReweightDriver::cvPars

private

Definition at line 123 of file ReweightDriver.h.

bool NeutrinoFluxReweight::ReweightDriver::doMIPPNumi

Definition at line 111 of file ReweightDriver.h.

std::string NeutrinoFluxReweight::ReweightDriver::fileOptions

private

Definition at line 126 of file ReweightDriver.h.

int NeutrinoFluxReweight::ReweightDriver::iUniv

private

Definition at line 122 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::meson_inc_wgt

Meson incident weights

Definition at line 87 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::mipp_kaon_wgt

MIPP NuMI kaons yield weight

Definition at line 58 of file ReweightDriver.h.

MIPPNumiKaonYieldsReweighter* NeutrinoFluxReweight::ReweightDriver::MIPP_NUMI_KAON_Universe

Definition at line 93 of file ReweightDriver.h.

MIPPNumiPionYieldsReweighter* NeutrinoFluxReweight::ReweightDriver::MIPP_NUMI_PION_Universe

Definition at line 92 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::mipp_pion_wgt

MIPP NuMI yield weight

Definition at line 55 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::nC_pi_wgt

Thin target neutron on carbon producing pions weights

Definition at line 78 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::nuA_wgt

nuA

Definition at line 84 of file ReweightDriver.h.

OtherReweighter* NeutrinoFluxReweight::ReweightDriver::OTHER_Universe

Definition at line 108 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::other_wgt

Any other hadronic interaction not corrected yet

Definition at line 90 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::pC_k_wgt

Thin target proton on carbon producing kaons weights

Definition at line 75 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::pC_nu_wgt

Thin target proton on carbon producing nucleons weights

Definition at line 81 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::pC_pi_wgt

Thin target proton on carbon producing pions weights

Definition at line 72 of file ReweightDriver.h.

TargetAttenuationReweighter* NeutrinoFluxReweight::ReweightDriver::TARG_ATT_Universe

Definition at line 94 of file ReweightDriver.h.

ThinTargetMesonIncidentReweighter* NeutrinoFluxReweight::ReweightDriver::THINTARGET_MESON_INCIDENT_Universe

Definition at line 104 of file ReweightDriver.h.

ThinTargetnCPionReweighter* NeutrinoFluxReweight::ReweightDriver::THINTARGET_NC_PION_Universe

Definition at line 102 of file ReweightDriver.h.

ThinTargetnucleonAReweighter* NeutrinoFluxReweight::ReweightDriver::THINTARGET_NUCLEON_A_Universe

Definition at line 106 of file ReweightDriver.h.

ThinTargetpCKaonReweighter* NeutrinoFluxReweight::ReweightDriver::THINTARGET_PC_KAON_Universe

Definition at line 101 of file ReweightDriver.h.

ThinTargetpCNucleonReweighter* NeutrinoFluxReweight::ReweightDriver::THINTARGET_PC_NUCLEON_Universe

Definition at line 103 of file ReweightDriver.h.

ThinTargetpCPionReweighter* NeutrinoFluxReweight::ReweightDriver::THINTARGET_PC_PION_Universe

Definition at line 100 of file ReweightDriver.h.

double NeutrinoFluxReweight::ReweightDriver::tot_abs_wgt

Definition at line 69 of file ReweightDriver.h.

const ParameterTable& NeutrinoFluxReweight::ReweightDriver::univPars

private

Definition at line 124 of file ReweightDriver.h.

AbsorptionDPIPReweighter* NeutrinoFluxReweight::ReweightDriver::VOL_ABS_DPIP_Universe

Definition at line 96 of file ReweightDriver.h.

AbsorptionDVOLReweighter* NeutrinoFluxReweight::ReweightDriver::VOL_ABS_DVOL_Universe

Definition at line 97 of file ReweightDriver.h.

AbsorptionICReweighter* NeutrinoFluxReweight::ReweightDriver::VOL_ABS_IC_Universe

Definition at line 95 of file ReweightDriver.h.

NucleonAbsorptionOutOfTargetReweighter* NeutrinoFluxReweight::ReweightDriver::VOL_ABS_NUCLEON_Universe

Definition at line 98 of file ReweightDriver.h.

OtherAbsorptionOutOfTargetReweighter* NeutrinoFluxReweight::ReweightDriver::VOL_ABS_OTHER_Universe

Definition at line 99 of file ReweightDriver.h.

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

ppfx/include/ReweightDriver.h
ppfx/src/ReweightDriver.cpp

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation