herwig is hosted by Hepforge, IPPP Durham
Herwig 7.3.0
EtaPiPiFermionsDecayer.h
1// -*- C++ -*-
2//
3// EtaPiPiFermionsDecayer.h is a part of Herwig - A multi-purpose Monte Carlo event generator
4// Copyright (C) 2002-2019 The Herwig Collaboration
5//
6// Herwig is licenced under version 3 of the GPL, see COPYING for details.
7// Please respect the MCnet academic guidelines, see GUIDELINES for details.
8//
9#ifndef HERWIG_EtaPiPiFermionsDecayer_H
10#define HERWIG_EtaPiPiFermionsDecayer_H
11// This is the declaration of the EtaPiPiFermionsDecayer class.
12
13#include "Herwig/Utilities/Kinematics.h"
14#include "Herwig/Decay/DecayIntegrator.h"
15#include "Herwig/Decay/PhaseSpaceMode.h"
16#include "Herwig/Decay/FormFactors/OmnesFunction.fh"
17#include "ThePEG/Helicity/LorentzSpinorBar.h"
18
19namespace Herwig {
20using namespace ThePEG;
21
57
58public:
59
64
71 virtual int modeNumber(bool & cc, tcPDPtr parent,
72 const tPDVector & children) const;
73
83 double me2(const int ichan,const Particle & part,
84 const tPDVector & outgoing,
85 const vector<Lorentz5Momentum> & momenta,
86 MEOption meopt) const;
87
91 virtual void constructSpinInfo(const Particle & part,
93
99 virtual void dataBaseOutput(ofstream & os,bool header) const;
100
101public:
102
110
116 void persistentInput(PersistentIStream & is, int version);
118
122 static void Init();
123
124protected:
125
132 virtual IBPtr clone() const {return new_ptr(*this);}
133
138 virtual IBPtr fullclone() const {return new_ptr(*this);}
140
141protected:
142
145
151 virtual void doinit();
152
156 virtual void doinitrun();
158
159public:
160
164 string setUpDecayMode(string arg);
165
166private:
167
172
173private:
174
178 Energy fPi_;
179
183 vector<int> incoming_;
184
188 vector<int> lepton_;
189
193 vector<double> coupling_;
194
198 vector<double> maxWeight_;
199
203 vector<int> option_;
204
208 InvEnergy2 aConst_;
209
213 double cConst_;
214
218 Energy mRho_;
219
223 Energy rhoWidth_;
224
228 double rhoConst_;
229
233 Energy mPi_;
234
239
243 OmnesFunctionPtr omnesFunction_;
244
249
253 mutable vector<Helicity::LorentzSpinor <SqrtEnergy> > wave_;
254
258 mutable vector<Helicity::LorentzSpinorBar<SqrtEnergy> > wavebar_;
259};
260}
261
262#endif /* HERWIG_EtaPiPiFermionsDecayer_H */
Main class for Decayers implementing multi-channel phase space integration.
MEOption
Enum for the matrix element option.
The EtaPiPiFermionsDecayer class implements the decay of the or to using either a VMD type model o...
virtual void dataBaseOutput(ofstream &os, bool header) const
Output the setup information for the particle database.
string setUpDecayMode(string arg)
Set the parameters for a decay mode.
void persistentInput(PersistentIStream &is, int version)
Function used to read in object persistently.
void persistentOutput(PersistentOStream &os) const
Function used to write out object persistently.
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
virtual IBPtr clone() const
Make a simple clone of this object.
vector< double > coupling_
Coupling for the decay, .
virtual void doinitrun()
Initialize this object to the begining of the run phase.
virtual int modeNumber(bool &cc, tcPDPtr parent, const tPDVector &children) const
Which of the possible decays is required.
double me2(const int ichan, const Particle &part, const tPDVector &outgoing, const vector< Lorentz5Momentum > &momenta, MEOption meopt) const
Return the matrix element squared for a given mode and phase-space channel.
vector< int > option_
The option for the energy dependence of the prefactor.
RhoDMatrix rho_
Spin densit matrix.
Energy fPi_
the pion decay constant, .
vector< int > incoming_
the PDG code for the incoming particle
virtual void constructSpinInfo(const Particle &part, ParticleVector outgoing) const
Construct the SpinInfos for the particles produced in the decay.
vector< Helicity::LorentzSpinorBar< SqrtEnergy > > wavebar_
Barred spinors for the fermions.
EtaPiPiFermionsDecayer()
Default constructor.
EtaPiPiFermionsDecayer & operator=(const EtaPiPiFermionsDecayer &)=delete
Private and non-existent assignment operator.
virtual void doinit()
Initialize this object after the setup phase before saving and EventGenerator to disk.
double rhoConst_
Constant for the running width.
static void Init()
Standard Init function used to initialize the interfaces.
bool localParameters_
Use local values of the parameters.
vector< Helicity::LorentzSpinor< SqrtEnergy > > wave_
Spinors for the fermions.
vector< double > maxWeight_
The maximum weight.
double cConst_
The constants for the Omnes function form.
InvEnergy2 aConst_
The constants for the omnes function form.
vector< int > lepton_
The PDG code of the leptons.
OmnesFunctionPtr omnesFunction_
Object calculating the Omnes function.
-*- C++ -*-
ThePEG::Ptr< InterfacedBase >::pointer IBPtr
vector< tPDPtr > tPDVector
vector< PPtr > ParticleVector
ThePEG::Ptr< ParticleData >::transient_const_pointer tcPDPtr