MCReweight.cxx

Go to the documentation of this file.

#ifndef mcreweight_cxx
#define mcreweight_cxx
#include "MCReweight/MCReweight.h"
#include "MessageService/MsgService.h"
#include <cassert>

CVSID("$Id: MCReweight.cxx,v 1.5 2005/04/13 19:17:03 cbs Exp $");

MCReweight* MCReweight::fgInstance = 0;

//*****************************************************
MCReweight::MCReweight()
{
  fWeightCalculators.clear();    
}

//*****************************************************
MCReweight::~MCReweight()
{
  MSG("MCReweight",Msg::kInfo) << "Closing down MCReweight object" 
                               << endl;
}

//*****************************************************
MCReweight& MCReweight::Instance()
{
  static Cleaner cleaner;
  
  //all copied from Calibrator package...
  if (!fgInstance) {
    cleaner.UseMe();           // dummy call to quiet compiler warnings
    fgInstance = new MCReweight();
    if(!fgInstance){
      MSG("MCReweight", Msg::kError)
        << "No MCReweight Instance - fatal." << endl;
      assert(fgInstance); // Kill job is there is no instance
    }    
  }
  return *fgInstance;
}

//*****************************************************
double MCReweight::ComputeWeight(Registry *event, Registry *rwtconfig){
  
  std::vector<WeightCalculator*>::iterator beg = fWeightCalculators.begin();
  std::vector<WeightCalculator*>::iterator end = fWeightCalculators.end();
  double totalWeight = 1;
  while(beg!=end) {
    (*beg)->SetReweightConfig(rwtconfig);
    totalWeight *= (*beg)->GetWeight(event);
    beg++;
  }
  return totalWeight;
}

//*****************************************************
double MCReweight::ComputeWeight(MCEventInfo *event, NuParent *parent,
                                 Registry *rwtconfig){
  
  std::vector<WeightCalculator*>::iterator beg = fWeightCalculators.begin();
  std::vector<WeightCalculator*>::iterator end = fWeightCalculators.end();
  double totalWeight = 1;
  while(beg!=end) {
    (*beg)->SetReweightConfig(rwtconfig);
    totalWeight *= (*beg)->GetWeight(event,parent);
    beg++;
  }
  return totalWeight;
}

//*****************************************************
void MCReweight::AddWeightCalculator(WeightCalculator *wc){

  std::vector<WeightCalculator*>::iterator beg = fWeightCalculators.begin();
  std::vector<WeightCalculator*>::iterator end = fWeightCalculators.end();
  while(beg!=end){
    if((*beg)->GetName()==wc->GetName()){
      MSG("MCReweight", Msg::kError) 
        << "Attempting to add " << wc->GetName() << " more than once! " 
        << "This is not a good idea, so not doing it."
        << endl;
      return;
    }
    beg++;
  }
  wc->SetStandardConfig(&fStandardConfig);
  fWeightCalculators.push_back(wc);
  
}

//*****************************************************
void MCReweight::ClearWeightCalculators(){
  
  fWeightCalculators.clear();
  
}

//*****************************************************
void MCReweight::ResetAllReweightConfigs() {
  std::vector<WeightCalculator*>::iterator beg = fWeightCalculators.begin();
  std::vector<WeightCalculator*>::iterator end = fWeightCalculators.end();
  while(beg!=end){
    (*beg)->ReweightConfigReset();
    beg++;
  }
}

//*****************************************************
void MCReweight::SetStandardConfig(Registry *config) {
  if(!config) return;
  fStandardConfig = *config;
  std::vector<WeightCalculator*>::iterator beg = fWeightCalculators.begin();
  std::vector<WeightCalculator*>::iterator end = fWeightCalculators.end();
  while(beg!=end){
    (*beg)->SetStandardConfig(&fStandardConfig);
    beg++;
  } 
}

//*****************************************************
void MCReweight::PrintReweightConfig(ostream & stream) {
  std::vector<WeightCalculator*>::iterator beg = fWeightCalculators.begin();
  std::vector<WeightCalculator*>::iterator end = fWeightCalculators.end();
  while(beg!=end){
    (*beg)->PrintReweightConfig(stream);
    beg++;
  }
}
#endif

---