herwig is hosted by Hepforge, IPPP Durham
Herwig 7.3.0
SemiLeptonicScalarDecayer.h
1// -*- C++ -*-
2//
3// SemiLeptonicScalarDecayer.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_SemiLeptonicScalarDecayer_H
10#define HERWIG_SemiLeptonicScalarDecayer_H
11//
12// This is the declaration of the SemiLeptonicScalarDecayer class.
13//
14#include "Herwig/Decay/DecayIntegrator.h"
15#include "Herwig/Decay/FormFactors/ScalarFormFactor.h"
16#include "Herwig/Decay/WeakCurrents/LeptonNeutrinoCurrent.h"
17#include "ThePEG/Helicity/LorentzTensor.h"
18
19namespace Herwig {
20using namespace ThePEG;
21
38
39public:
40
45
52 virtual bool accept(tcPDPtr parent, const tPDVector & children) const;
53
60 virtual int modeNumber(bool & cc, tcPDPtr parent,
61 const tPDVector & children) const ;
62
72 double me2(const int ichan,const Particle & part,
73 const tPDVector & outgoing,
74 const vector<Lorentz5Momentum> & momenta,
75 MEOption meopt) const;
76
80 virtual void constructSpinInfo(const Particle & part,
82
88 virtual void dataBaseOutput(ofstream & os,bool header) const;
89
90public:
91
99
105 void persistentInput(PersistentIStream & is, int version);
107
111 static void Init();
112
113protected:
114
121 virtual IBPtr clone() const {return new_ptr(*this);}
122
127 virtual IBPtr fullclone() const {return new_ptr(*this);}
129
130protected:
131
134
140 virtual void doinit();
141
145 virtual void doinitrun();
147
148private:
149
154
155private:
156
160 LeptonNeutrinoCurrentPtr _current;
161
165 ScalarFormFactorPtr _form;
166
170 vector<double> _maxwgt;
171
175 vector<int> _modemap;
176
181
185 mutable vector<Helicity::LorentzPolarizationVector> _vectors;
186
190 mutable vector<Helicity::LorentzTensor<double> > _tensors;
191
195 mutable vector<unsigned int> _constants;
196
200 mutable vector<PDT::Spin> _ispin;
201
205 mutable unsigned int _imes;
206};
207
208}
209
210
211#endif /* HERWIG_SemiLeptonicScalarDecayer_H */
Main class for Decayers implementing multi-channel phase space integration.
MEOption
Enum for the matrix element option.
The SemiLeptonicScalarDecayer class is designed for the semileptonic decay of a (pesudo)scalar meson ...
vector< Helicity::LorentzTensor< double > > _tensors
Polarization vectors for the decay products.
SemiLeptonicScalarDecayer & operator=(const SemiLeptonicScalarDecayer &)=delete
Private and non-existent assignment operator.
ScalarFormFactorPtr _form
The form factor.
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
virtual bool accept(tcPDPtr parent, const tPDVector &children) const
Check if this decayer can perfom the decay for a particular mode.
vector< PDT::Spin > _ispin
Spins of the particles.
void persistentOutput(PersistentOStream &os) const
Function used to write out object persistently.
virtual void doinitrun()
Initialize this object to the begining of the run phase.
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.
LeptonNeutrinoCurrentPtr _current
The leptonic current.
virtual IBPtr clone() const
Make a simple clone of this object.
virtual void doinit()
Initialize this object after the setup phase before saving and EventGenerator to disk.
vector< Helicity::LorentzPolarizationVector > _vectors
Polarization vectors for the decay products.
virtual int modeNumber(bool &cc, tcPDPtr parent, const tPDVector &children) const
Which of the possible decays is required.
SemiLeptonicScalarDecayer()
Default constructor.
unsigned int _imes
Location of the outgoing meson.
void persistentInput(PersistentIStream &is, int version)
Function used to read in object persistently.
static void Init()
Standard Init function used to initialize the interfaces.
virtual void constructSpinInfo(const Particle &part, ParticleVector outgoing) const
Construct the SpinInfos for the particles produced in the decay.
virtual void dataBaseOutput(ofstream &os, bool header) const
Output the setup information for the particle database.
RhoDMatrix _rho
Spin density matrix.
vector< int > _modemap
mapping of the mode to the form factor
vector< unsigned int > _constants
Constants for the mapping of the leptonic current.
vector< double > _maxwgt
the maximum weight for the integration of a given decay
-*- C++ -*-
ThePEG::Ptr< InterfacedBase >::pointer IBPtr
vector< tPDPtr > tPDVector
vector< PPtr > ParticleVector
ThePEG::Ptr< ParticleData >::transient_const_pointer tcPDPtr