In This Package:

GtInverseBeta.h

Go to the documentation of this file.

#ifndef GTINVERSEBETA_H
#define GTINVERSEBETA_H


#include "GenTools/IHepMCEventMutator.h"
#include "GaudiAlg/GaudiTool.h"
#include "HepMC/GenParticle.h"
#include <CLHEP/Vector/ThreeVector.h>
#include "TF2.h"

#include <string>
using namespace CLHEP;

class GtInverseBeta : public GaudiTool,
  virtual public IHepMCEventMutator
{
 public:
  
  GtInverseBeta(const std::string& type,
                const std::string& name,
                const IInterface* parent);
  virtual ~GtInverseBeta();
  
  virtual StatusCode initialize();
  virtual StatusCode finalize();
  
  // HepMCEventMutator interface
  virtual StatusCode mutate(HepMC::GenEvent& event);
  
  void GenerateInverseBetaPrimaries(TF2*, const double,
                                    Hep3Vector&, Hep3Vector&, Hep3Vector&);
  
 private:
  TF2 *funcInvBetaProb;

  double m_neutrino_angle_in_deg;
  bool m_positron_only;
  bool m_neutron_only;
  
  HepMC::GenParticle* make_particle(Hep3Vector&, double, int);
  
};


#endif  // GTINVERSEBETA_H

---