NeutrinoFluxReweight::MakeReweight Class Reference

A class to make the reweight event by event. More...

#include <MakeReweight.h>

Public Member Functions
void	SetOptions (std::string fileIn)
void	calculateWeights (nu_g4numi *nu, const char *tgtcfg, const char *horncfg)
void	calculateWeights (bsim::Dk2Nu *nu, bsim::DkMeta *meta)
	create an interaction chain from the new dk2nu(dkmeta) format More...
std::vector< double >	GetTotalWeights ()
	total weights More...
int	GetNumberOfUniversesUsed ()
	number of universes used in this run More...
std::vector< double >	GetWeights (std::string nameReweighter)
	get the vector of the weights for a given reweighter More...
double	GetCVWeight ()
	get the cv weights More...

Static Public Member Functions
static MakeReweight *	getInstance ()
static void	resetInstance ()

Public Attributes
std::vector< ReweightDriver * >	vec_rws
	vector of Reweighter Drivers, one per universe More...
ReweightDriver *	cv_rw
	Reweighter Drivers for the central value. More...

Private Member Functions
	MakeReweight ()
	~MakeReweight ()
void	Initialize ()
void	ParseOptions ()
void	doTheJob (InteractionChainData *icd)

Private Attributes
std::vector< ParameterTable >	cvPars
std::vector< ParameterTable >	univPars
int	Nuniverses
std::string	mippCorrOption
std::string	fileOptions
std::vector< double >	vec_wgts
std::map< std::string, std::vector< double > >	map_rew_wgts
double	cv_wgt

Static Private Attributes
static MakeReweight *	instance = 0

Detailed Description

A class to make the reweight event by event.

Definition at line 32 of file MakeReweight.h.

Constructor & Destructor Documentation

NeutrinoFluxReweight::MakeReweight::MakeReweight ( )

private

Definition at line 9 of file MakeReweight.cpp.

                            {
  }

NeutrinoFluxReweight::MakeReweight::~MakeReweight ( )

private

Definition at line 160 of file MakeReweight.cpp.

                             {
    
  }

Member Function Documentation

void NeutrinoFluxReweight::MakeReweight::calculateWeights	(	nu_g4numi *	nu,
		const char *	tgtcfg,
		const char *	horncfg
	)

calculate the weights for this event using the old flux ntuple format (filled in nu_g4numi object), the tgtcfg (the target configuration, example: le010z) and the horncfg (horn configuration, example: 185i)

Definition at line 99 of file MakeReweight.cpp.

                                                                                           {
    InteractionChainData inter_chain(nu,tgtcfg,horncfg);
    doTheJob(&inter_chain);
  }

void NeutrinoFluxReweight::MakeReweight::calculateWeights	(	bsim::Dk2Nu *	nu,
		bsim::DkMeta *	meta
	)

create an interaction chain from the new dk2nu(dkmeta) format

Definition at line 104 of file MakeReweight.cpp.

                                                                    {  
    InteractionChainData inter_chain(nu,meta);
    doTheJob(&inter_chain);
  }

void NeutrinoFluxReweight::MakeReweight::doTheJob ( InteractionChainData *  icd )

private

Parse seeting options

Definition at line 109 of file MakeReweight.cpp.

                                                      {
  
    //universe calculation:
    map_rew_wgts.clear();
    for(int ii=0;ii<Nuniverses;ii++){
      vec_wgts[ii] = vec_rws[ii]->calculateWeight(*icd);

      map_rew_wgts["MIPPNumiPionYields"].push_back(vec_rws[ii]->mipp_pion_wgt);
      map_rew_wgts["MIPPNumiKaonYields"].push_back(vec_rws[ii]->mipp_kaon_wgt);
      
      map_rew_wgts["TargetAttenuation"].push_back(vec_rws[ii]->att_wgt);            
      map_rew_wgts["AbsorptionIC"].push_back(vec_rws[ii]->abs_ic_wgt);
      map_rew_wgts["AbsorptionDPIP"].push_back(vec_rws[ii]->abs_dpip_wgt);
      map_rew_wgts["AbsorptionDVOL"].push_back(vec_rws[ii]->abs_dvol_wgt);
      map_rew_wgts["AbsorptionNucleon"].push_back(vec_rws[ii]->abs_nucleon_wgt);
      map_rew_wgts["AbsorptionOther"].push_back(vec_rws[ii]->abs_other_wgt);
      map_rew_wgts["TotalAbsorption"].push_back(vec_rws[ii]->tot_abs_wgt);
      
      map_rew_wgts["ThinTargetpCPion"].push_back(vec_rws[ii]->pC_pi_wgt);
      map_rew_wgts["ThinTargetpCKaon"].push_back(vec_rws[ii]->pC_k_wgt);
      map_rew_wgts["ThinTargetnCPion"].push_back(vec_rws[ii]->nC_pi_wgt);
      map_rew_wgts["ThinTargetpCNucleon"].push_back(vec_rws[ii]->pC_nu_wgt);
      map_rew_wgts["ThinTargetMesonIncident"].push_back(vec_rws[ii]->meson_inc_wgt);
      map_rew_wgts["ThinTargetnucleonA"].push_back(vec_rws[ii]->nuA_wgt);
      map_rew_wgts["Other"].push_back(vec_rws[ii]->other_wgt);

    }
    
    //cv calculation:
    cv_wgt = cv_rw->calculateWeight(*icd);
    
  }

double NeutrinoFluxReweight::MakeReweight::GetCVWeight

(

)

get the cv weights

Definition at line 153 of file MakeReweight.cpp.

                                  {   
    return cv_wgt;
  }

MakeReweight * NeutrinoFluxReweight::MakeReweight::getInstance ( )

static

Definition at line 164 of file MakeReweight.cpp.

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

int NeutrinoFluxReweight::MakeReweight::GetNumberOfUniversesUsed

(

)

number of universes used in this run

Definition at line 157 of file MakeReweight.cpp.

                                             {
    return Nuniverses;
  }

std::vector< double > NeutrinoFluxReweight::MakeReweight::GetTotalWeights

(

)

total weights

Definition at line 95 of file MakeReweight.cpp.

                                                 {
    return vec_wgts;
  }

std::vector< double > NeutrinoFluxReweight::MakeReweight::GetWeights

(

std::string

nameReweighter

)

get the vector of the weights for a given reweighter

Definition at line 142 of file MakeReweight.cpp.

                                                                    {

    std::vector<double> tmp_vec;
    std::map<std::string,std::vector<double> >::iterator it = map_rew_wgts.begin();
    it = map_rew_wgts.find(nameReweighter);
    if(it!=map_rew_wgts.end())tmp_vec = it->second;
    
    return tmp_vec;
   
  }

void NeutrinoFluxReweight::MakeReweight::Initialize ( void  )

private

Inititalize the job and configurethe ReweighterDrivier

Definition at line 28 of file MakeReweight.cpp.

                               {
    
    //Getting MIPP binning and the correlation parameters:
    CentralValuesAndUncertainties* cvu = CentralValuesAndUncertainties::getInstance();;
    MIPPNumiYieldsBins*  myb =  MIPPNumiYieldsBins::getInstance(); 
    ThinTargetBins*  thinbin =  ThinTargetBins::getInstance(); 
    MIPPNumiMC*  mymc =  MIPPNumiMC::getInstance(); 
    const char* ppfxDir = getenv("PPFX_DIR");
       
    std::cout<<"Initializing correlation parameters"<<std::endl;
    cvu->readFromXML(Form("%s/uncertainties/Parameters_%s.xml",ppfxDir,mippCorrOption.c_str()));
    
    std::cout<<"Initializing bin data conventions"<<std::endl;
    myb->pip_data_from_xml(Form("%s/data/BINS/MIPPNumiData_PIP_Bins.xml",ppfxDir));
    myb->pim_data_from_xml(Form("%s/data/BINS/MIPPNumiData_PIM_Bins.xml",ppfxDir));
    myb->k_pi_data_from_xml(Form("%s/data/BINS/MIPPNumiData_K_PI_Bins.xml",ppfxDir));
    thinbin->pC_pi_from_xml(Form("%s/data/BINS/ThinTarget_pC_pi_Bins.xml",ppfxDir));
    thinbin->barton_pC_pi_from_xml(Form("%s/data/BINS/ThinTargetBarton_pC_pi_Bins.xml",ppfxDir));
    thinbin->pC_k_from_xml(Form("%s/data/BINS/ThinTargetLowxF_pC_k_Bins.xml",ppfxDir));
    thinbin->mipp_pC_k_pi_from_xml(Form("%s/data/BINS/ThinTarget_K_PI_Bins.xml",ppfxDir));
    thinbin->meson_incident_from_xml(Form("%s/data/BINS/ThinTarget_MesonIncident.xml",ppfxDir));
    thinbin->pC_p_from_xml(Form("%s/data/BINS/ThinTarget_pC_p_Bins.xml",ppfxDir));
    thinbin->pC_n_from_xml(Form("%s/data/BINS/ThinTarget_pC_n_Bins.xml",ppfxDir));
    thinbin->material_scaling_from_xml(Form("%s/data/BINS/ThinTarget_material_scaling_Bins.xml",ppfxDir));
    
    std::cout<<"Initializing MC values"<<std::endl;
    mymc->pip_mc_from_xml(Form("%s/data/MIPP/MIPPNuMI_MC_PIP.xml",ppfxDir));
    mymc->pim_mc_from_xml(Form("%s/data/MIPP/MIPPNuMI_MC_PIM.xml",ppfxDir));
    mymc->kap_mc_from_xml(Form("%s/data/MIPP/MIPPNuMI_MC_KAP.xml",ppfxDir));
    mymc->kam_mc_from_xml(Form("%s/data/MIPP/MIPPNuMI_MC_KAM.xml",ppfxDir));
    mymc->k0l_mc_from_xml(Form("%s/data/MIPP/MIPPNuMI_MC_K0L.xml",ppfxDir));
    mymc->k0s_mc_from_xml(Form("%s/data/MIPP/MIPPNuMI_MC_K0S.xml",ppfxDir));

    //Reweighter drivers:
    vec_rws.reserve(Nuniverses);
    std::cout<<"Initializing reweight drivers for "<<Nuniverses<<" universes"<<std::endl;
    
    const int base_universe=1000000;
    // cvPars.reserve(Nuniverses+1);
    univPars.reserve(Nuniverses+1);

    ParameterTable cvPars=cvu->getCVPars();

    std::cout<<"Loading parameters for universe: ";
    for(int ii=0;ii<Nuniverses;ii++){
      std::cout << ii << " " << std::flush;
      // cvPars.push_back(cvu->getCVPars());
      univPars.push_back(cvu->calculateParsForUniverse(ii+base_universe));
      vec_rws.push_back(new ReweightDriver(ii,cvPars,univPars[ii],fileOptions));
      vec_wgts.push_back(1.0);
    }
    std::cout << std::endl;
    
    //Central Value driver:
    //by convention, we use universe -1 to hold the cv. It is in agreement with CentralValueAndUncertainties  
    const int cv_id = -1;
    std::cout<<"Loading parameters for cv"<<std::endl;
    // cvPars.push_back(cvu->getCVPars());
    univPars.push_back(cvu->calculateParsForUniverse(cv_id));
    cv_rw = new ReweightDriver(cv_id,cvPars,univPars[Nuniverses],fileOptions); 

    cv_wgt = 1.0;

    std::cout<<"Done configuring universes"<<std::endl;

  }

void NeutrinoFluxReweight::MakeReweight::ParseOptions ( )

private

Parse seeting options

Definition at line 17 of file MakeReweight.cpp.

                                 {
    //Parsing the file input:
    using boost::property_tree::ptree;
    ptree top;
    read_xml(fileOptions.c_str(),top,2); // option 2 removes comment strings
    ptree& options = top.get_child("inputs.Settings");

    mippCorrOption = options.get<std::string>("MIPPCorrOption");
    Nuniverses = options.get<int>("NumberOfUniverses");
    
  }

void NeutrinoFluxReweight::MakeReweight::resetInstance ( )

static

Definition at line 169 of file MakeReweight.cpp.

   {
      delete instance; 
      instance = 0; 
   }

void NeutrinoFluxReweight::MakeReweight::SetOptions

(

std::string

fileIn

)

Definition at line 11 of file MakeReweight.cpp.

                                              {
    fileOptions = fileIn;
    ParseOptions();
    Initialize();
  }

Member Data Documentation

ReweightDriver* NeutrinoFluxReweight::MakeReweight::cv_rw

Reweighter Drivers for the central value.

Definition at line 68 of file MakeReweight.h.

double NeutrinoFluxReweight::MakeReweight::cv_wgt

private

Definition at line 86 of file MakeReweight.h.

std::vector<ParameterTable> NeutrinoFluxReweight::MakeReweight::cvPars

private

Definition at line 80 of file MakeReweight.h.

std::string NeutrinoFluxReweight::MakeReweight::fileOptions

private

Definition at line 83 of file MakeReweight.h.

MakeReweight * NeutrinoFluxReweight::MakeReweight::instance = 0

staticprivate

Definition at line 88 of file MakeReweight.h.

std::map<std::string,std::vector<double> > NeutrinoFluxReweight::MakeReweight::map_rew_wgts

private

Definition at line 85 of file MakeReweight.h.

std::string NeutrinoFluxReweight::MakeReweight::mippCorrOption

private

Definition at line 82 of file MakeReweight.h.

int NeutrinoFluxReweight::MakeReweight::Nuniverses

private

Definition at line 81 of file MakeReweight.h.

std::vector<ParameterTable> NeutrinoFluxReweight::MakeReweight::univPars

private

Definition at line 80 of file MakeReweight.h.

std::vector<ReweightDriver*> NeutrinoFluxReweight::MakeReweight::vec_rws

vector of Reweighter Drivers, one per universe

Definition at line 65 of file MakeReweight.h.

std::vector<double> NeutrinoFluxReweight::MakeReweight::vec_wgts

private

Definition at line 84 of file MakeReweight.h.

---

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

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