Phokhara.h

Go to the documentation of this file.

//*****************************************************************************
//
// Generator/Phokhara/Phokhara.h
//
// Algorithm runs event generator Phokhara (hep-ph/0710.4227v1)
// and stores output to transient store
//
// Nov 2007 Original BES3 code by Alexey Zhemchugov
//
//*****************************************************************************
 
#ifndef GENERATORMODULESPHOKHARA_H
#define GENERATORMODULESPHOKHARA_H
 
#include "HepMC/GenEvent.h"
#include "HepMC/GenParticle.h"
 
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/ISvcLocator.h"
 
#include <vector>
 
using namespace std;
 
class Phokhara : public Algorithm {
public:
  Phokhara( const string& name, ISvcLocator* pSvcLocator );
 
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();
  StatusCode storeParticles();
 
private:
  // jobOption params
 
  int    m_nm;     // # of events to determine the maximum
  int    m_nlo;    // Born(0), NLO(1)
  double m_w;      // soft photon cutoff
  int    m_pion;   // mu+mu-(0),pi+pi-(1),2pi0pi+pi-(2),2pi+2pi-(3),ppbar(4),nnbar(5),
                   // K+K-(6),K0K0bar(7),pi+pi-pi0(8), Lamb Lambbar->pi-pi+ppbar(9)
  int m_fsr;       // ISR only(0), ISR+FSR(1), ISR+INT+FSR(2)
  int m_fsrnlo;    // yes(1), no(0)
  int m_ivac;      // no(0), yes - Jegerlehner(1), yes - Teubner(2)
  int m_tagged;    // tagged photons(0), untagged photons(1)
  int m_NarrowRes; // none(0), JPsi (1), Psip(2)
  int m_FF_Kaon; // KaonFormFactor: constrained (0), unconstrained (1), Kuhn-Khodjamirian-Bruch
                 // (2)
  int    m_FF_Pion; // KS Pionformfactor(0), GS Pionformfactor old(1), GS Pionformfactor new(2)
  int    m_f0_model; // f0+f0(600): KK model(0), no structure(1), no f0+f0(600)(2), f0 KLOE(3)
  double m_E;        // CMS-energy
  double m_q2min;    // minimal  hadrons(muons)-gamma-inv mass squared
  double m_q2_min_c; // minimal inv. mass squared of the hadrons(muons)
  double m_q2_max_c; // maximal inv. mass squared of the hadrons(muons)
  double m_gmin;     // minimal photon energy/missing energy
  double m_phot1cut; // minimal photon angle/missing momentum angle
  double m_phot2cut; // maximal photon angle/missing momentum angle
  double m_pi1cut;   // minimal hadrons(muons) angle
  double m_pi2cut;   // maximal hadrons(muons) angle
 
  double qqmin, qqmax;
  double cos1min, cos1max, cos2min, cos2max, cos3min, cos3max;
  double dsigm1, dsigm2, sigma1, sigma2, sigma, dsigm, Ar[14], Ar_r[14];
 
  int ievent;
 
  // Initial Seed
  long int m_initSeed;
};
 
#endif