/Users/jmonk/Physics/ForIA/ForIA/AnalysisTools/ReWeightMultiplicity.hh

Go to the documentation of this file.

#ifndef FORIA_REWEIGHT_MULTIPLICITY
#define FORIA_REWEIGHT_MULTIPLICITY

#include "ForIA/AnalysisTools/TruthSelection.hh"

#include <vector>

namespace ForIA{
  
  using std::vector;
  
  class Event;
  
  class MultiplicityFunction{
    
  public:
    virtual double fn(size_t mult)const{return 1.;}
    
    virtual ~MultiplicityFunction(){}
    
  };
  
  class ReWeightMultiplicity{
    
  public:
    
    ReWeightMultiplicity();
    
    void initialise();
    
    void initialise(const MultiplicityFunction *fn);
    
    double weight(const Event &evt);
    
    double weight(size_t nTracks);
    
  private:
    
    // The weight for a given (offset) multiplicity
    vector<double> m_weights;
    
    // the offset multiplicity - the multiplicity of the first bin in the 
    // vector/histogram of multiplicity
    size_t m_startMultiplicity;
    
    // the highest multiplicity for which we have a weight
    size_t m_endMultiplicity;

    // selects the truth particles from the event
    TruthSelection m_truthSelection;

    
  };
}


#endif