BeamEnergyCalculator.h

Go to the documentation of this file.

#ifndef BEAM_ENERGY_CALCULATOR_H
#define BEAM_ENERGY_CALCULATOR_H

// BeamEnergyCalculator
// 
// Calculates weights for the different beam energy configurations
// as well as a 1/E flux.  
//
//
// Created: M. Kordosky  -- June 2, 2005
//
// $Author: rhatcher $ 
//
// $Revision: 1.2 $
// 
// $Name:  $
//
// $Id: BeamEnergyCalculator.h,v 1.2 2007/01/30 03:13:14 rhatcher Exp $
//
#ifndef __CINT__
#include <ext/hash_map>
#endif
#include <string>


#include "MCReweight/WeightCalculator.h"
#include "Conventions/BeamType.h"
#include "Conventions/Detector.h"

class Registry;
class MCEventInfo;
class NuParent;
class TH1;
class TDirectory;



class BeamEnergyCalculator : public WeightCalculator {
 public:
  BeamEnergyCalculator(Registry* stdconfig=0);
  //  virtual ~BeamEnergyCalculator();
  virtual void Config();
  virtual void ReweightConfigReset();
  virtual double GetWeight(Registry *eventinfo);
  virtual double GetWeight(MCEventInfo *event,NuParent *parent);
  virtual double GetWeight(BeamType::BeamType_t bt, 
                           Detector::Detector_t det,
                           int nu_idhep, double nu_energy, 
                           double high_energy_limit=120.0,
                           double low_energy_limit=low_energy_cut); 

 private:
  static const double low_energy_cut;
  double DeriveWeight(TH1*,double,double,double);
  double DeriveInverseEWeight(double,double,double);
  static void ConstructName(Detector::Detector_t det, 
                            int nu_idhep, std::string& s);
  TH1* GetHist(BeamType::BeamType_t bt, 
               Detector::Detector_t det, int nu_idhep);

  TDirectory* fFluxDir;
  std::string fFluxFileName;
  bool fFluxFileChanged;

  Detector::Detector_t fDetector;
  BeamType::BeamType_t fBeam;
  double fLowRange;
  double fHighRange;


#ifndef __CINT__
// there is no default hashing function for std::string
// but one exists for char*
//
  struct myhash
  {
    size_t operator()( const std::string& x ) const
    {
      return __gnu_cxx::hash< const char* >()( x.c_str() );
    }
  };

  typedef __gnu_cxx::hash_map<std::string, TH1*,myhash> CacheMap;
  typedef __gnu_cxx::hash_map<std::string, TH1*,myhash>::iterator CacheMapIter;
  CacheMap fCache;
#endif




  //  Int_t fNuID; // stdhep id of neutrino
  //  double fNuEnergy;



};

#endif

---