herwig is hosted by Hepforge, IPPP Durham
Herwig  7.2.1
ScalarVectorVectorDecayer.h
1 // -*- C++ -*-
2 //
3 // ScalarVectorVectorDecayer.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_ScalarVectorVectorDecayer_H
10 #define HERWIG_ScalarVectorVectorDecayer_H
11 //
12 // This is the declaration of the ScalarVectorVectorDecayer class.
13 //
14 
15 #include "Herwig/Decay/DecayIntegrator.h"
16 #include "Herwig/Decay/DecayPhaseSpaceMode.h"
17 #include "ThePEG/Helicity/LorentzPolarizationVector.h"
18 
19 namespace Herwig {
20 
21 using namespace ThePEG;
22 
38 
39 public:
40 
45 
52  virtual int modeNumber(bool & cc, tcPDPtr parent,
53  const tPDVector & children) const;
54 
63  double me2( const int ichan,const Particle & part,
64  const ParticleVector & decay, MEOption meopt) const;
65 
74  bool twoBodyMEcode(const DecayMode & dm, int & mecode, double & coupling) const;
75 
81  virtual void dataBaseOutput(ofstream & os,bool header) const;
82 
83 public:
84 
91  void persistentOutput(PersistentOStream & os) const;
92 
98  void persistentInput(PersistentIStream & is, int version);
100 
107  static void Init();
108 
109 protected:
110 
117  virtual IBPtr clone() const {return new_ptr(*this);}
118 
123  virtual IBPtr fullclone() const {return new_ptr(*this);}
125 
126 protected:
127 
135  virtual void doinit();
136 
141  virtual void doinitrun();
143 
144 private:
145 
150  ScalarVectorVectorDecayer & operator=(const ScalarVectorVectorDecayer &) = delete;
151 
152 private:
153 
157  vector<int> _incoming;
158 
162  vector<int> _outgoing1;
163 
167  vector<int> _outgoing2;
168 
172  vector<InvEnergy> _coupling;
173 
177  vector<double> _maxweight;
178 
182  unsigned int _initsize;
183 
187  mutable RhoDMatrix _rho;
188 
192  mutable vector<Helicity::LorentzPolarizationVector> _vectors[2];
193 };
194 
195 }
196 
197 #endif /* HERWIG_ScalarVectorVectorDecayer_H */
RhoDMatrix _rho
Spin density matrix.
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
unsigned int _initsize
initial number of modes
vector< int > _outgoing2
the PDG code for the second outgoing particle
ThePEG::Ptr< InterfacedBase >::pointer IBPtr
vector< double > _maxweight
the maximum weight for the decay
MEOption
Enum for the matrix element option.
vector< int > _outgoing1
the PDG code for the first outgoing particle
vector< tPDPtr > tPDVector
Main class for Decayers implementing multi-channel phase space integration.
ThePEG::Ptr< ParticleData >::transient_const_pointer tcPDPtr
vector< int > _incoming
the PDG code for the incoming particle
virtual IBPtr clone() const
Make a simple clone of this object.
-*- C++ -*-
vector< PPtr > ParticleVector
vector< InvEnergy > _coupling
the coupling for the decay, .
The ScalarVectorVectorDecayer class is designed to perform the decay of a scalar meson to two spin-1 ...