FluxWeights.cxx

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// \file    FluxWeights.cxx
// \brief   Store flux weigths for neutrino correction
// \author  Leo Aliaga - laliaga@fnal.gov
////////////////////////////////////////////////////////////////////////

#include "MCReweight/FluxWeights.h"

namespace fxwgt
{
  FluxWeights::FluxWeights(){
    fcv_hp = 1;
    fcv_bf = 1;
    fcv_ad = 1;
    funiv_hp.clear();
    funiv_bf.clear();
    funiv_ad.clear();
    fnuniv = 0;
  }
  /////
  void FluxWeights::SetHadronProductionCentralValue(float cv){
    fcv_hp = cv;
  }
  void FluxWeights::SetBeamFocusingCentralValue(float cv){
    fcv_bf = cv;
  }
  void FluxWeights::SetHadronProductionMultiUniverses(std::vector<float> vuniv){
    funiv_hp = vuniv;
  }
  void FluxWeights::SetBeamFocusingMultiUniverses(std::vector<float> vuniv){
    funiv_bf = vuniv;
  }
  void FluxWeights::SetAdditionalCentralValue(float cv){
    fcv_ad = cv;    
  }
  void FluxWeights::SetAdditionalMultiUniverses(std::vector<float> vuniv){
    funiv_ad = vuniv;    
  }
  ////    
  float FluxWeights::GetHadronProductionCentralValue() const{
    return fcv_hp;
  }
  float FluxWeights::GetBeamFocusingCentralValue() const{
    return fcv_bf;
  }
  float FluxWeights::GetTotalCentralValue() const{
    return fcv_hp*fcv_bf;
  }
  float FluxWeights::GetAdditionalCentralValue() const{
    return fcv_ad;
  }  
  std::vector<float> FluxWeights::GetHadronProductionMultiUniverses() const{
    return funiv_hp;
  }
  std::vector<float> FluxWeights::GetBeamFocusingMultiUniverses() const{
    return funiv_bf;
  }
  std::vector<float> FluxWeights::GetTotalMultiUniverses() const{
    std::vector<float> tmp_vec;
    for(unsigned int ii=0;ii<FluxWeights::GetNumberOfUniverses();ii++){
      tmp_vec.push_back(funiv_hp[ii]*funiv_bf[ii]);
    }
    return tmp_vec;
  }
  std::vector<float> FluxWeights::GetAdditionalMultiUniverses() const{
    return funiv_ad;
  }
  unsigned int FluxWeights::GetNumberOfUniverses() const{
    return funiv_hp.size();
  }
}