herwig is hosted by Hepforge, IPPP Durham
Herwig 7.3.0
GeneralCurrentDecayer.h
1// -*- C++ -*-
2//
3// GeneralCurrentDecayer.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_GeneralCurrentDecayer_H
10#define HERWIG_GeneralCurrentDecayer_H
11//
12// This is the declaration of the GeneralCurrentDecayer class.
13//
14
15#include "Herwig/Decay/DecayIntegrator.h"
16#include "Herwig/Decay/WeakCurrents/WeakCurrent.h"
17#include "Herwig/Decay/PhaseSpaceMode.h"
18#include "ThePEG/Helicity/Vertex/VertexBase.h"
19#include "GeneralCurrentDecayer.fh"
20
21namespace Herwig {
22using namespace ThePEG;
23using Helicity::VertexBasePtr;
24
32
33public:
34
39 maxmass_(5.*GeV), wgtmax_(0.) {}
40
49 virtual int modeNumber(bool & cc, tcPDPtr parent,const tPDVector & children) const;
50
57 virtual Energy partialWidth(tPDPtr inpart, tPDPtr outa,
58 vector<tPDPtr> currout) = 0;
60
64 void setDecayInfo(PDPtr in, PDPtr out, const vector<tPDPtr> & outCurrent,
65 VertexBasePtr vertex, WeakCurrentPtr current,
66 Energy maxmass);
67
68public:
69
77
83 void persistentInput(PersistentIStream & is, int version);
85
92 static void Init();
93
94protected:
95
103 virtual void doinit();
104
109 virtual void doinitrun();
111
112protected:
113
119 int modeNumber(bool & cc, vector<long> id) const;
120
125 unsigned int mode() const { return mode_; }
126
130 WeakCurrentPtr weakCurrent() const { return current_; }
131
136 VertexBasePtr vertex() const { return theVertex_; }
137
138private:
139
145
146private:
147
151 VertexBasePtr theVertex_;
152
157
162
166 vector<tPDPtr> currentOut_;
167
171 WeakCurrentPtr current_;
172
176 Energy maxmass_;
177
181 unsigned int mode_;
182
187
191 double wgtmax_;
192
196 vector<double> weights_;
197};
198
199}
200
201#endif /* HERWIG_GeneralCurrentDecayer_H */
Main class for Decayers implementing multi-channel phase space integration.
Here is the documentation of the GeneralCurrentDecayer class.
GeneralCurrentDecayer & operator=(const GeneralCurrentDecayer &)=delete
The assignment operator is private and must never be called.
VertexBasePtr vertex() const
Get vertex pointer.
WeakCurrentPtr current_
Pointer to the current.
unsigned int mode() const
Access to the map between the number of the mode and the modes in the current.
virtual int modeNumber(bool &cc, tcPDPtr parent, const tPDVector &children) const
Which of the possible decays is required.
unsigned int mode_
mapping of the modes to the currents
Energy maxmass_
Maximum mass difference.
WeakCurrentPtr weakCurrent() const
Access to the weak current.
static void Init()
The standard Init function used to initialize the interfaces.
int modeNumber(bool &cc, vector< long > id) const
The number of the mode.
void setDecayInfo(PDPtr in, PDPtr out, const vector< tPDPtr > &outCurrent, VertexBasePtr vertex, WeakCurrentPtr current, Energy maxmass)
set up the decay
void persistentOutput(PersistentOStream &os) const
Function used to write out object persistently.
GeneralCurrentDecayer()
The default constructor.
int wgtloc_
location of the weights
PDPtr outpart_
First outgoing particle.
virtual Energy partialWidth(tPDPtr inpart, tPDPtr outa, vector< tPDPtr > currout)=0
Function to return partial Width.
double wgtmax_
the maximum weight
void persistentInput(PersistentIStream &is, int version)
Function used to read in object persistently.
VertexBasePtr theVertex_
Pointer to vertex set in inheriting class.
vector< double > weights_
The weights for the different channels.
vector< tPDPtr > currentOut_
Outgoing particles from the current.
virtual void doinit()
Initialize this object after the setup phase before saving an EventGenerator to disk.
virtual void doinitrun()
Initialize this object.
-*- C++ -*-
ThePEG::Ptr< ParticleData >::transient_pointer tPDPtr
vector< tPDPtr > tPDVector
ThePEG::Ptr< ParticleData >::transient_const_pointer tcPDPtr
ThePEG::Ptr< ParticleData >::pointer PDPtr