herwig is hosted by Hepforge, IPPP Durham
Herwig  7.2.1
ScalarMesonTensorScalarDecayer.h
1 // -*- C++ -*-
2 //
3 // ScalarMesonTensorScalarDecayer.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_ScalarMesonTensorScalarDecayer_H
10 #define HERWIG_ScalarMesonTensorScalarDecayer_H
11 //
12 // This is the declaration of the ScalarMesonTensorScalarDecayer class.
13 //
14 #include "Herwig/Decay/DecayIntegrator.h"
15 #include "Herwig/Decay/DecayPhaseSpaceMode.h"
16 #include "ThePEG/Helicity/LorentzTensor.h"
17 
18 namespace Herwig {
19 using namespace Herwig;
20 
36 
37 public:
38 
43 
50  virtual int modeNumber(bool & cc, tcPDPtr parent,
51  const tPDVector & children) const;
52 
61  double me2( const int ichan,const Particle & part,
62  const ParticleVector & decay, MEOption meopt) const;
63 
72  bool twoBodyMEcode(const DecayMode & dm, int & mecode, double & coupling) const;
73 
79  virtual void dataBaseOutput(ofstream & os,bool header) const;
80 
81 public:
82 
89  void persistentOutput(PersistentOStream & os) const;
90 
96  void persistentInput(PersistentIStream & is, int version);
98 
102  static void Init();
103 
104 protected:
105 
112  virtual IBPtr clone() const {return new_ptr(*this);}
113 
118  virtual IBPtr fullclone() const {return new_ptr(*this);}
120 
121 protected:
122 
130  virtual void doinit();
131 
135  virtual void doinitrun();
137 
138 private:
139 
144 
145  private:
146 
150  vector<int> _incoming;
151 
155  vector<int> _outgoingT;
156 
160  vector<int> _outgoingS;
161 
165  vector<InvEnergy> _coupling;
166 
170  vector<double> _maxweight;
171 
175  unsigned int _initsize;
176 
180  mutable RhoDMatrix _rho;
181 
185  mutable vector<Helicity::LorentzTensor<double> > _tensors;
186 };
187 
188 }
189 
190 
191 #endif /* HERWIG_ScalarMesonTensorScalarDecayer_H */
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.
ScalarMesonTensorScalarDecayer & operator=(const ScalarMesonTensorScalarDecayer &)=delete
Private and non-existent assignment operator.
vector< double > _maxweight
the maximum weight for the decay
virtual int modeNumber(bool &cc, tcPDPtr parent, const tPDVector &children) const
Which of the possible decays is required.
virtual void dataBaseOutput(ofstream &os, bool header) const
Output the setup information for the particle database.
void persistentOutput(PersistentOStream &os) const
Function used to write out object persistently.
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.
void persistentInput(PersistentIStream &is, int version)
Function used to read in object persistently.
ThePEG::Ptr< InterfacedBase >::pointer IBPtr
ScalarMesonTensorScalarDecayer()
Default constructor.
vector< InvEnergy > _coupling
the coupling for the decay
MEOption
Enum for the matrix element option.
vector< Helicity::LorentzTensor< double > > _tensors
Polarization tensors for the decay product.
vector< tPDPtr > tPDVector
vector< int > _outgoingT
the PDG code for the tensor meson
The ScalarMesonTensorScalarDecayer class is designed for the decay of a (pseudo)scalar meson to a ten...
Main class for Decayers implementing multi-channel phase space integration.
ThePEG::Ptr< ParticleData >::transient_const_pointer tcPDPtr
vector< int > _outgoingS
the PDG code for the scalar meson
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...
virtual void doinit()
Initialize this object after the setup phase before saving and EventGenerator to disk.
vector< int > _incoming
the PDG code for the incoming particle
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
-*- C++ -*-
virtual IBPtr clone() const
Make a simple clone of this object.
unsigned int _initsize
initial number of modes
virtual void doinitrun()
Initialize this object to the begining of the run phase.
vector< PPtr > ParticleVector