herwig is hosted by Hepforge, IPPP Durham
Herwig 7.3.0
SFFDecayer.h
1// -*- C++ -*-
2//
3// SFFDecayer.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_SFFDecayer_H
10#define HERWIG_SFFDecayer_H
11//
12// This is the declaration of the SFFDecayer class.
13//
14
15#include "GeneralTwoBodyDecayer.h"
16#include "ThePEG/Repository/EventGenerator.h"
17#include "ThePEG/Helicity/Vertex/Scalar/FFSVertex.h"
18#include "ThePEG/Helicity/Vertex/Scalar/VSSVertex.h"
19#include "ThePEG/Helicity/Vertex/Vector/FFVVertex.h"
20
21namespace Herwig {
22using namespace ThePEG;
23using Helicity::FFSVertexPtr;
24
35
36public:
37
47 double me2(const int ichan,const Particle & part,
48 const tPDVector & outgoing,
49 const vector<Lorentz5Momentum> & momenta,
50 MEOption meopt) const;
51
55 virtual void constructSpinInfo(const Particle & part,
57
64 virtual Energy partialWidth(PMPair inpart, PMPair outa,
65 PMPair outb) const;
66
71 POWHEGType output = FSR;
72 for(auto vertex : vertex_) {
73 if(vertex->orderInAllCouplings()!=1) {
74 output = No;
75 break;
76 }
77 }
78 return output;
79 }
80
84 virtual double threeBodyME(const int , const Particle & inpart,
85 const ParticleVector & decay,
87 MEOption meopt);
88
92 void identifyVertices(const int iferm, const int ianti,
93 const Particle & inpart, const ParticleVector & decay,
94 AbstractFFVVertexPtr & outgoingVertexF,
95 AbstractFFVVertexPtr & outgoingVertexA,
96 ShowerInteraction inter);
97
102 vector<VertexBasePtr>,
103 map<ShowerInteraction,VertexBasePtr> &,
104 const vector<map<ShowerInteraction,VertexBasePtr> > &,
105 map<ShowerInteraction,VertexBasePtr>);
107
108public:
109
117
123 void persistentInput(PersistentIStream & is, int version);
125
132 static void Init();
133
134protected:
135
142 virtual IBPtr clone() const;
143
148 virtual IBPtr fullclone() const;
150
151private:
152
157 SFFDecayer & operator=(const SFFDecayer &) = delete;
158
159private:
160
164 vector<AbstractFFSVertexPtr> vertex_;
165
169 vector<FFSVertexPtr> perturbativeVertex_;
170
174 map<ShowerInteraction,AbstractVSSVertexPtr> incomingVertex_;
175
179 map<ShowerInteraction,AbstractFFVVertexPtr> outgoingVertex1_;
180
184 map<ShowerInteraction,AbstractFFVVertexPtr> outgoingVertex2_;
185
190
195
199 mutable vector<SpinorWaveFunction> wave_;
200
204 mutable vector<SpinorBarWaveFunction> wavebar_;
205
210
215
219 mutable vector<SpinorWaveFunction> wave3_;
220
224 mutable vector<SpinorBarWaveFunction> wavebar3_;
225
229 mutable vector<VectorWaveFunction> gluon_;
230
231
232};
233
234}
235
236#endif /* HERWIG_SFFDecayer_H */
MEOption
Enum for the matrix element option.
The GeneralTwoBodyDecayer class is designed to be the base class for 2 body decays for some general m...
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...
pair< tcPDPtr, Energy > PMPair
A ParticleData ptr and (possible) mass pair.
POWHEGType
Virtual members to be overridden by inheriting classes which implement hard corrections.
Definition: HwDecayerBase.h:76
The SFFDecayer class implements the decay of a scalar to 2 fermions in a general model.
Definition: SFFDecayer.h:34
map< ShowerInteraction, AbstractFFVVertexPtr > outgoingVertex2_
Abstract pointer to AbstractFFVVertex for QCD radiation from outgoing (anti)fermion.
Definition: SFFDecayer.h:184
vector< SpinorWaveFunction > wave3_
Spinor wavefunction for 3 body decay.
Definition: SFFDecayer.h:219
void identifyVertices(const int iferm, const int ianti, const Particle &inpart, const ParticleVector &decay, AbstractFFVVertexPtr &outgoingVertexF, AbstractFFVVertexPtr &outgoingVertexA, ShowerInteraction inter)
Indentify outgoing vertices for the fermion and antifermion.
vector< SpinorBarWaveFunction > wavebar3_
Barred spinor wavefunction for 3 body decay.
Definition: SFFDecayer.h:224
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
static void Init()
The standard Init function used to initialize the interfaces.
RhoDMatrix rho_
Spin density matrix.
Definition: SFFDecayer.h:189
ScalarWaveFunction swave3_
Scalar wavefunction for 3 body decay.
Definition: SFFDecayer.h:214
vector< AbstractFFSVertexPtr > vertex_
Abstract pointer to AbstractFFSVertex.
Definition: SFFDecayer.h:164
virtual Energy partialWidth(PMPair inpart, PMPair outa, PMPair outb) const
Function to return partial Width.
virtual IBPtr clone() const
Make a simple clone of this object.
vector< FFSVertexPtr > perturbativeVertex_
Pointer to the perturbative vertex.
Definition: SFFDecayer.h:169
virtual void constructSpinInfo(const Particle &part, ParticleVector outgoing) const
Construct the SpinInfos for the particles produced in the decay.
virtual POWHEGType hasPOWHEGCorrection()
Has a POWHEG style correction.
Definition: SFFDecayer.h:70
map< ShowerInteraction, AbstractVSSVertexPtr > incomingVertex_
Abstract pointer to AbstractVSSVertex for QCD radiation from incoming scalar.
Definition: SFFDecayer.h:174
SFFDecayer & operator=(const SFFDecayer &)=delete
The assignment operator is private and must never be called.
void persistentOutput(PersistentOStream &os) const
Function used to write out object persistently.
map< ShowerInteraction, AbstractFFVVertexPtr > outgoingVertex1_
Abstract pointer to AbstractFFVVertex for QCD radiation from outgoing (anti)fermion.
Definition: SFFDecayer.h:179
vector< SpinorWaveFunction > wave_
Spinor wavefunction.
Definition: SFFDecayer.h:199
void persistentInput(PersistentIStream &is, int version)
Function used to read in object persistently.
RhoDMatrix rho3_
Spin density matrix for 3 body decay.
Definition: SFFDecayer.h:209
ScalarWaveFunction swave_
Scalar wavefunction.
Definition: SFFDecayer.h:194
virtual double threeBodyME(const int, const Particle &inpart, const ParticleVector &decay, ShowerInteraction inter, MEOption meopt)
Three-body matrix element including additional QCD radiation.
vector< VectorWaveFunction > gluon_
Vector wavefunction for 3 body decay.
Definition: SFFDecayer.h:229
vector< SpinorBarWaveFunction > wavebar_
Barred spinor wavefunction.
Definition: SFFDecayer.h:204
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.
virtual void setDecayInfo(PDPtr incoming, PDPair outgoing, vector< VertexBasePtr >, map< ShowerInteraction, VertexBasePtr > &, const vector< map< ShowerInteraction, VertexBasePtr > > &, map< ShowerInteraction, VertexBasePtr >)
Set the information on the decay.
ShowerInteraction
Handy header file to be included in all Shower classes.
-*- C++ -*-
ThePEG::Ptr< InterfacedBase >::pointer IBPtr
vector< tPDPtr > tPDVector
vector< PPtr > ParticleVector
ThePEG::Ptr< ParticleData >::pointer PDPtr
pair< PDPtr, PDPtr > PDPair