herwig is hosted by Hepforge, IPPP Durham
Herwig  7.2.1
ScalarScalarScalarDecayer.h
1 // -*- C++ -*-
2 //
3 // ScalarScalarScalarDecayer.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_ScalarScalarScalarDecayer_H
10 #define HERWIG_ScalarScalarScalarDecayer_H
11 // This is the declaration of the ScalarScalarScalarDecayer class.
12 
13 #include "Herwig/Decay/DecayIntegrator.h"
14 #include "Herwig/Decay/DecayPhaseSpaceMode.h"
15 
16 namespace Herwig {
17 using namespace ThePEG;
18 
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 
143  ScalarScalarScalarDecayer & operator=(const ScalarScalarScalarDecayer &) = delete;
144 
145 private:
146 
150  vector<int> _incoming;
151 
155  vector<int> _outgoing1;
156 
160  vector<int> _outgoing2;
161 
165  vector<Energy> _coupling;
166 
170  vector<double> _maxweight;
171 
175  unsigned int _initsize;
176 
180  mutable RhoDMatrix _rho;
181 };
182 
183 }
184 
185 #endif /* HERWIG_ScalarScalarScalarDecayer_H */
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
vector< double > _maxweight
the maximum weight for the decay
vector< int > _outgoing1
the PDG code for the first outgoing particle
ThePEG::Ptr< InterfacedBase >::pointer IBPtr
virtual IBPtr clone() const
Make a simple clone of this object.
unsigned int _initsize
initial number of modes
RhoDMatrix _rho
Spin density matrix.
MEOption
Enum for the matrix element option.
vector< tPDPtr > tPDVector
Main class for Decayers implementing multi-channel phase space integration.
ThePEG::Ptr< ParticleData >::transient_const_pointer tcPDPtr
vector< int > _outgoing2
the PDG code for the second incoming particle
-*- C++ -*-
vector< int > _incoming
the PDG code for the incoming particle
The ScalarScalarScalarDecayer class is designed for the decay of a scalar meson to two scalar mesons...
vector< PPtr > ParticleVector
vector< Energy > _coupling
the coupling for the decay