herwig is hosted by Hepforge, IPPP Durham
Herwig  7.2.1
VectorMesonVectorPScalarDecayer.h
1 // -*- C++ -*-
2 //
3 // VectorMesonVectorPScalarDecayer.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_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 
159 private:
160 
164  vector<InvEnergy> _coupling;
165 
169  vector<int> _incoming;
170 
174  vector<int> _outgoingV;
175 
179  vector<int> _outgoingP;
180 
184  vector<double> _maxweight;
185 
189  unsigned int _initsize;
190 
194  mutable RhoDMatrix _rho;
195 
199  mutable vector<Helicity::LorentzPolarizationVector> _vectors[2];
200 
201 };
202 
203 }
204 
205 
206 #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.
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.