herwig is hosted by Hepforge, IPPP Durham
Herwig 7.3.0
Spin3MesonTensorVectorDecayer.h
1// -*- C++ -*-
2//
3// Spin3MesonTensorVectorDecayer.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_Spin3MesonTensorVectorDecayer_H
10#define HERWIG_Spin3MesonTensorVectorDecayer_H
11//
12// This is the declaration of the Spin3MesonTensorVectorDecayer class.
13//
14#include "Herwig/Decay/DecayIntegrator.h"
15#include "ThePEG/Helicity/LorentzPolarizationVector.h"
16#include "Herwig/Decay/PhaseSpaceMode.h"
17
18namespace Herwig {
19using namespace ThePEG;
20
47
48public:
49
56 virtual int modeNumber(bool & cc, tcPDPtr parent,
57 const tPDVector & children) const;
58
68 double me2(const int ichan,const Particle & part,
69 const tPDVector & outgoing,
70 const vector<Lorentz5Momentum> & momenta,
71 MEOption meopt) const;
72
76 virtual void constructSpinInfo(const Particle & part,
78
87 bool twoBodyMEcode(const DecayMode & dm, int & mecode, double & coupling) const;
88
94 virtual void dataBaseOutput(ofstream & os,bool header) const;
95
96public:
97
105
111 void persistentInput(PersistentIStream & is, int version);
113
117 static void Init();
118
119protected:
120
127 virtual IBPtr clone() const {return new_ptr(*this);}
128
133 virtual IBPtr fullclone() const {return new_ptr(*this);}
135
136protected:
137
145 virtual void doinit();
146
150 virtual void doinitrun();
152
153private:
154
159
160public:
161
165 string setUpDecayMode(string arg);
166
167private:
168
172 vector<int> incoming_;
173
177 vector<pair<long,long> > outgoing_;
178
182 vector<Energy> coupling_;
183
187 vector<double> maxWeight_;
188
193 mutable vector<Helicity::LorentzRank3Tensor<double> > rank3_;
194
198 mutable vector<Helicity::LorentzTensor<double> > tensors_;
199
203 mutable vector<Helicity::LorentzPolarizationVector> vectors_;
204
209
210};
211
212}
213
214#endif /* HERWIG_Spin3MesonTensorVectorDecayer_H */
Main class for Decayers implementing multi-channel phase space integration.
MEOption
Enum for the matrix element option.
The Spin3MesonTensorVectorDecayer class is designed for the decay of a rank 3 tensor meson to a tenso...
static void Init()
Standard Init function used to initialize the interfaces.
virtual IBPtr clone() const
Make a simple clone of this object.
bool twoBodyMEcode(const DecayMode &dm, int &mecode, double &coupling) const
Specify the matrix element to be used in the running width calculation.
vector< pair< long, long > > outgoing_
the PDG codes the outgoing particles (vector,scalar)
double me2(const int ichan, const Particle &part, const tPDVector &outgoing, const vector< Lorentz5Momentum > &momenta, MEOption meopt) const
Return the matrix element squared for a given mode and phase-space channel.
vector< int > incoming_
the PDG codes for the incoming particles
virtual int modeNumber(bool &cc, tcPDPtr parent, const tPDVector &children) const
Which of the possible decays is required.
void persistentInput(PersistentIStream &is, int version)
Function used to read in object persistently.
vector< double > maxWeight_
the maximum weight for the decay
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
vector< Helicity::LorentzPolarizationVector > vectors_
Storage of the polarization vectors.
string setUpDecayMode(string arg)
Set the parameters for a decay mode.
virtual void doinit()
Initialize this object after the setup phase before saving and EventGenerator to disk.
vector< Energy > coupling_
the coupling for the decay
vector< Helicity::LorentzTensor< double > > tensors_
Polarization tensors for the decay product.
Spin3MesonTensorVectorDecayer & operator=(const Spin3MesonTensorVectorDecayer &)=delete
Private and non-existent assignment operator.
virtual void doinitrun()
Initialize this object to the begining of the run phase.
virtual void dataBaseOutput(ofstream &os, bool header) const
Output the setup information for the particle database.
virtual void constructSpinInfo(const Particle &part, ParticleVector outgoing) const
Construct the SpinInfos for the particles produced in the decay.
void persistentOutput(PersistentOStream &os) const
Function used to write out object persistently.
vector< Helicity::LorentzRank3Tensor< double > > rank3_
Storage of polarization tensors to try and increase speed.
-*- C++ -*-
ThePEG::Ptr< InterfacedBase >::pointer IBPtr
vector< tPDPtr > tPDVector
vector< PPtr > ParticleVector
ThePEG::Ptr< ParticleData >::transient_const_pointer tcPDPtr