herwig is hosted by Hepforge, IPPP Durham
Herwig  7.1.5
VectorMesonVectorPScalarDecayer.h
1 // -*- C++ -*-
2 //
3 // VectorMesonVectorPScalarDecayer.h is a part of Herwig - A multi-purpose Monte Carlo event generator
4 // Copyright (C) 2002-2017 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_VectorMesonVectorPScalarDecayer_H
10 #define HERWIG_VectorMesonVectorPScalarDecayer_H
11 // This is the declaration of the VectorMesonVectorPScalarDecayer class.
12 
13 #include "Herwig/Decay/DecayIntegrator.h"
14 #include "Herwig/Decay/DecayPhaseSpaceMode.h"
15 #include "ThePEG/Helicity/LorentzPolarizationVector.h"
16 
17 namespace Herwig {
18 
50 
51 public:
52 
57 
64  virtual int modeNumber(bool & cc, tcPDPtr parent,
65  const tPDVector & children) const;
66 
75  double me2(const int ichan,const Particle & part,
76  const ParticleVector & decay, MEOption meopt) const;
77 
86  bool twoBodyMEcode(const DecayMode & dm, int & mecode, double & coupling) const;
87 
93  virtual void dataBaseOutput(ofstream & os,bool header) const;
94 
95 public:
96 
103  void persistentOutput(PersistentOStream & os) const;
104 
110  void persistentInput(PersistentIStream & is, int version);
112 
116  static void Init();
117 
118 protected:
119 
126  virtual IBPtr clone() const {return new_ptr(*this);}
127 
132  virtual IBPtr fullclone() const {return new_ptr(*this);}
134 
135 protected:
136 
144  virtual void doinit();
145 
149  virtual void doinitrun();
151 
152 private:
153 
158 
163 
164 private:
165 
169  vector<InvEnergy> _coupling;
170 
174  vector<int> _incoming;
175 
179  vector<int> _outgoingV;
180 
184  vector<int> _outgoingP;
185 
189  vector<double> _maxweight;
190 
194  unsigned int _initsize;
195 
199  mutable RhoDMatrix _rho;
200 
204  mutable vector<Helicity::LorentzPolarizationVector> _vectors[2];
205 
206 };
207 
208 }
209 
210 
211 #include "ThePEG/Utilities/ClassTraits.h"
212 
213 namespace ThePEG {
214 
221 template <>
224  typedef Herwig::DecayIntegrator NthBase;
225 };
226 
231 template <>
233  : public ClassTraitsBase<Herwig::VectorMesonVectorPScalarDecayer> {
235  static string className() { return "Herwig::VectorMesonVectorPScalarDecayer"; }
241  static string library() { return "HwVMDecay.so"; }
242 
243 };
244 
247 }
248 #endif /* HERWIG_VectorMesonVectorPScalarDecayer_H */
double me2(const int ichan, const Particle &part, const ParticleVector &decay, MEOption meopt) const
Return the matrix element squared for a given mode and phase-space channel.
virtual void doinitrun()
Initialize this object to the begining of the run phase.
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
vector< int > _outgoingV
the PDG codes for the outgoing vector
vector< InvEnergy > _coupling
coupling for a decay
ThePEG::Ptr< InterfacedBase >::pointer IBPtr
VectorMesonVectorPScalarDecayer()
Default constructor.
VectorMesonVectorPScalarDecayer & operator=(const VectorMesonVectorPScalarDecayer &)=delete
Private and non-existent assignment operator.
static void Init()
Standard Init function used to initialize the interfaces.
bool twoBodyMEcode(const DecayMode &dm, int &mecode, double &coupling) const
Specify the matrix element to be used in the running width calculation.
static ClassDescription< VectorMesonVectorPScalarDecayer > initVectorMesonVectorPScalarDecayer
Describe a concrete class with persistent data.
vector< int > _outgoingP
the PDG codes for the outgoing pseudoscalar
MEOption
Enum for the matrix element option.
void persistentInput(PersistentIStream &is, int version)
Function used to read in object persistently.
virtual int modeNumber(bool &cc, tcPDPtr parent, const tPDVector &children) const
Which of the possible decays is required.
virtual IBPtr clone() const
Make a simple clone of this object.
This class is designed for the decay of a vector meson to another spin-1 particle, either another vector meson or a photon, and a pseduoscalar meson.
vector< tPDPtr > tPDVector
unsigned int _initsize
Initial size of the vectors.
Main class for Decayers implementing multi-channel phase space integration.
ThePEG::Ptr< ParticleData >::transient_const_pointer tcPDPtr
virtual void dataBaseOutput(ofstream &os, bool header) const
Output the setup information for the particle database.
virtual ParticleVector decay(const Particle &parent, const tPDVector &children) const
For a given decay mode and a given particle instance, perform the decay and return the decay products...
void persistentOutput(PersistentOStream &os) const
Function used to write out object persistently.
vector< int > _incoming
the PDG codes for the incoming particle
vector< double > _maxweight
maximum weight for a decay
-*- C++ -*-
vector< Helicity::LorentzPolarizationVector > _vectors[2]
Polarization vectors.
vector< PPtr > ParticleVector
virtual void doinit()
Initialize this object after the setup phase before saving and EventGenerator to disk.