herwig is hosted by Hepforge, IPPP Durham
Herwig  7.1.4
Herwig::ThreeMesonDefaultCurrent Class Reference

The ThreeMesonDefaultCurrent class implements the currents from Z.Phys.C58:445 (1992), this paper uses the form from Z.Phys.C48:445 (1990) for the $a_1$ width and is the default model in TAUOLA. More...

#include <ThreeMesonDefaultCurrent.h>

Inheritance diagram for Herwig::ThreeMesonDefaultCurrent:

Public Member Functions

 ThreeMesonDefaultCurrent ()
 Default constructor.
 
virtual void dataBaseOutput (ofstream &os, bool header, bool create) const
 Output the setup information for the particle database. More...
 
double threeBodyMatrixElement (const int imode, const Energy2 q2, const Energy2 s3, const Energy2 s2, const Energy2 s1, const Energy m1, const Energy m2, const Energy m3) const
 the matrix element for the $a_1$ decay to calculate the running width More...
 
Functions used by the persistent I/O system.
void persistentOutput (PersistentOStream &os) const
 Function used to write out object persistently. More...
 
void persistentInput (PersistentIStream &is, int version)
 Function used to read in object persistently. More...
 
Methods for the construction of the phase space integrator.
virtual bool createMode (int icharge, unsigned int imode, DecayPhaseSpaceModePtr mode, unsigned int iloc, unsigned int ires, DecayPhaseSpaceChannelPtr phase, Energy upp)
 Complete the construction of the decay mode for integration. More...
 
- Public Member Functions inherited from Herwig::ThreeMesonCurrentBase
 ThreeMesonCurrentBase ()
 Default constructor.
 
virtual vector< LorentzPolarizationVectorEcurrent (const int imode, const int ichan, Energy &scale, const ParticleVector &decay, DecayIntegrator::MEOption meopt) const
 Hadronic current. More...
 
virtual bool accept (vector< int > id)
 Accept the decay. More...
 
virtual unsigned int decayMode (vector< int > id)
 Return the decay mode number for a given set of particles in the current. More...
 
virtual tPDVector particles (int icharge, unsigned int imode, int iq, int ia)
 The particles produced by the current. More...
 
- Public Member Functions inherited from Herwig::WeakDecayCurrent
 WeakDecayCurrent ()
 Default constructor.
 
unsigned int numberOfModes () const
 Return the number of modes handled by this current.
 
void decayModeInfo (unsigned int imode, int &iq, int &ia) const
 Information on a decay mode. More...
 
void persistentOutput (PersistentOStream &os) const
 Function used to write out object persistently. More...
 
void persistentInput (PersistentIStream &is, int version)
 Function used to read in object persistently. More...
 
- Public Member Functions inherited from ThePEG::Interfaced
virtual bool defaultInit ()
 
PPtr getParticle (PID) const
 
PDPtr getParticleData (PID) const
 
bool used () const
 
void useMe () const
 
tEGPtr generator () const
 
void persistentOutput (PersistentOStream &os) const
 
void persistentInput (PersistentIStream &is, int version)
 
PPtr getParticle (PID) const
 
PDPtr getParticleData (PID) const
 
bool used () const
 
void useMe () const
 
tEGPtr generator () const
 
void persistentOutput (PersistentOStream &os) const
 
void persistentInput (PersistentIStream &is, int version)
 
- Public Member Functions inherited from ThePEG::InterfacedBase
string fullName () const
 
string name () const
 
string path () const
 
string comment () const
 
void setup (istream &is)
 
virtual void debugme () const
 
void update ()
 
void init ()
 
virtual bool preInitialize () const
 
void initrun ()
 
void finish ()
 
void touch ()
 
void reset ()
 
void clear ()
 
InitState state () const
 
bool locked () const
 
bool touched () const
 
void persistentOutput (PersistentOStream &os) const
 
void persistentInput (PersistentIStream &is, int version)
 
void update ()
 
void init ()
 
virtual bool preInitialize () const
 
void initrun ()
 
void finish ()
 
void touch ()
 
void reset ()
 
void clear ()
 
InitState state () const
 
bool locked () const
 
bool touched () const
 
void persistentOutput (PersistentOStream &os) const
 
void persistentInput (PersistentIStream &is, int version)
 
- Public Member Functions inherited from ThePEG::Base
void debug () const
 
- Public Member Functions inherited from ThePEG::Pointer::ReferenceCounted
CounterType referenceCount () const
 
- Public Member Functions inherited from ThePEG::Named
 Named (const string &newName=string())
 
const string & name () const
 
bool operator== (const Named &other) const
 
bool operator< (const Named &other) const
 

Static Public Member Functions

static void Init ()
 Standard Init function used to initialize the interfaces.
 
- Static Public Member Functions inherited from Herwig::ThreeMesonCurrentBase
static void Init ()
 Standard Init function used to initialize the interfaces.
 
- Static Public Member Functions inherited from Herwig::WeakDecayCurrent
static void Init ()
 Standard Init function used to initialize the interfaces.
 
- Static Public Member Functions inherited from ThePEG::Interfaced
static void Init ()
 
- Static Public Member Functions inherited from ThePEG::InterfacedBase
static void Init ()
 
- Static Public Member Functions inherited from ThePEG::Base
static void Init ()
 

Protected Member Functions

virtual bool acceptMode (int) const
 Can the current handle a particular set of mesons. More...
 
virtual FormFactors calculateFormFactors (const int ichan, const int imode, Energy2 q2, Energy2 s1, Energy2 s2, Energy2 s3) const
 Calculate the form factor for the current. More...
 
Clone Methods.
virtual IBPtr clone () const
 Make a simple clone of this object. More...
 
virtual IBPtr fullclone () const
 Make a clone of this object, possibly modifying the cloned object to make it sane. More...
 
Standard Interfaced functions.
virtual void doinit ()
 Initialize this object after the setup phase before saving and EventGenerator to disk. More...
 
virtual void doinitrun ()
 Initialize this object to the begining of the run phase.
 
virtual void doupdate ()
 Check sanity of the object during the setup phase.
 
- Protected Member Functions inherited from Herwig::WeakDecayCurrent
void addDecayMode (int iq, int ia)
 Add a decay mode to the list. More...
 
void setInitialModes (unsigned int nmodes)
 Set initial number of modes. More...
 
- Protected Member Functions inherited from ThePEG::Interfaced
void reporeg (IBPtr object, string name) const
 
bool setDefaultReference (PtrT &ptr, string classname, string objectname)
 
 Interfaced (const string &newName)
 
 Interfaced (const Interfaced &i)
 
void setGenerator (tEGPtr generator)
 
- Protected Member Functions inherited from ThePEG::InterfacedBase
 InterfacedBase (string newName)
 
 InterfacedBase (const InterfacedBase &i)
 
virtual void readSetup (istream &is)
 
virtual void dofinish ()
 
virtual IVector getReferences ()
 
virtual void rebind (const TranslationMap &)
 
virtual void readSetup (istream &is)
 
virtual void dofinish ()
 
virtual IVector getReferences ()
 
virtual void rebind (const TranslationMap &)
 
- Protected Member Functions inherited from ThePEG::Pointer::ReferenceCounted
 ReferenceCounted (const ReferenceCounted &)
 
ReferenceCountedoperator= (const ReferenceCounted &)
 
 ReferenceCounted (const ReferenceCounted &)
 
ReferenceCountedoperator= (const ReferenceCounted &)
 
- Protected Member Functions inherited from ThePEG::Named
const Namedoperator= (const Named &other)
 
const string & name (const string &newName)
 

Private Member Functions

ThreeMesonDefaultCurrentoperator= (const ThreeMesonDefaultCurrent &)
 Private and non-existent assignment operator.
 
Complex BrhoF123 (Energy2 q2, int ires) const
 The $\rho$ Breit-Wigner for the $F_{1,2,3}$ form factors. More...
 
Complex BrhoF5 (Energy2 q2, int ires) const
 The $\rho$ Breit-Wigner for the $F_5$ form factors. More...
 
Complex BKstarF123 (Energy2 q2, int ires) const
 The $K^*$ Breit-Wigner for the $F_{1,2,3}$ form factors. More...
 
Complex BKstarF5 (Energy2 q2, int ires) const
 The $K^*$ Breit-Wigner for the $F_5$ form factors. More...
 
Complex FKrho (Energy2 si, Energy2 sj, int ires) const
 Mixed Breit Wigner for the $F_5$ form factor. More...
 
Complex a1BreitWigner (Energy2 q2) const
 $a_1$ Breit-Wigner More...
 
Complex K1BreitWigner (Energy2 q2) const
 The $K_1$ Breit-Wigner. More...
 
Energy a1Width (Energy2 q2) const
 The $a_1$ running width. More...
 
double g (Energy2 q2) const
 The $g(Q^2)$ function of Kuhn and Santamaria.
 
void inita1Width (int iopt)
 Initialize the $a_1$ running width. More...
 
Complex rhoKBreitWigner (Energy2 q2, unsigned int itype, unsigned int ires) const
 Breit-Wigners for the $\rho$ and $K^*$. More...
 

Private Attributes

vector< double > _rhoF123wgts
 Parameters for the $\rho$ Breit-Wigner in the $F_{1,2,3}$ form factors.
 
vector< double > _kstarF123wgts
 Parameters for the $K^*$ Breit-Wigner in the $F_{1,2,3}$ form factors.
 
vector< double > _rhoF5wgts
 Parameters for the $\rho$ Breit-Wigner in the $F_5$ form factors.
 
vector< double > _kstarF5wgts
 Parameters for the $K^*$ Breit-Wigner in the $F_5$ form factors.
 
double _rhoKstarwgt
 The relative weight of the $\rho$ and $K^*$ where needed.
 
vector< Energy > _a1runwidth
 The $a_1$ width for the running $a_1$ width calculation.
 
vector< Energy2 > _a1runq2
 The $q^2$ for the running $a_1$ width calculation.
 
Interpolator< Energy, Energy2 >::Ptr _a1runinter
 The interpolator for the running $a_1$ width calculation.
 
bool _initializea1
 Initialize the running $a_1$ width.
 
Energy _a1mass
 The mass of the $a_1$ resonances.
 
Energy _a1width
 The width of the $a_1$ resonances.
 
Energy _k1mass
 The mass of the $aK1$ resonances.
 
Energy _k1width
 The width of the $K_1$ resonances.
 
Energy _fpi
 The pion decay constant, $f_\pi$.
 
Energy _mpi
 The pion mass.
 
Energy _mK
 The kaon mass.
 
bool _rhoparameters
 use local values of the $\rho$ masses and widths
 
vector< Energy > _rhoF123masses
 The $\rho$ masses for the $F_{1,2,3}$ form factors.
 
vector< Energy > _rhoF5masses
 The $\rho$ masses for the $F_5$ form factors.
 
vector< Energy > _rhoF123widths
 The $\rho$ widths for the $F_{1,2,3}$ form factors.
 
vector< Energy > _rhoF5widths
 The $\rho$ widths for the $F_5$ form factors.
 
bool _kstarparameters
 use local values of the $K^*$ resonances masses and widths
 
vector< Energy > _kstarF123masses
 The $K^*$ masses for the $F_{1,2,3}$ form factors.
 
vector< Energy > _kstarF5masses
 The $K^*$ masses for the $F_5$ form factors.
 
vector< Energy > _kstarF123widths
 The $K^*$ widths for the $F_{1,2,3}$ form factors.
 
vector< Energy > _kstarF5widths
 The $K^*$ widths for the $F_5$ form factors.
 
bool _a1parameters
 Use local values of the $a_1$ parameters.
 
bool _k1parameters
 Use local values of the $K_1$ parameters.
 
bool _a1opt
 Option for the $a_1$ width.
 
Energy _maxmass
 The maximum mass of the hadronic system.
 
Energy _maxcalc
 The maximum mass when the running width was calculated.
 

Friends

class Defaulta1MatrixElement
 The matrix element for the running $a_1$ width is a friend to keep some members private.
 

Additional Inherited Members

- Public Types inherited from ThePEG::InterfacedBase
enum  InitState
 
- Public Types inherited from ThePEG::Pointer::ReferenceCounted
typedef unsigned int CounterType
 
- Public Attributes inherited from ThePEG::InterfacedBase
 initializing
 
 uninitialized
 
 initialized
 
 runready
 
- Public Attributes inherited from ThePEG::Pointer::ReferenceCounted
const unsigned long uniqueId
 
- Static Protected Member Functions inherited from ThePEG::Interfaced
static void registerRepository (IBPtr)
 
static void registerRepository (IBPtr, string newName)
 

Detailed Description

The ThreeMesonDefaultCurrent class implements the currents from Z.Phys.C58:445 (1992), this paper uses the form from Z.Phys.C48:445 (1990) for the $a_1$ width and is the default model in TAUOLA.

The following three meson modes are implemented.

  • $ \pi^- \pi^- \pi^+ $, (imode=0)
  • $ \pi^0 \pi^0 \pi^- $, (imode=1)
  • $ K^- \pi^- K^+ $, (imode=2)
  • $ K^0 \pi^- \bar{K}^0$, (imode=3)
  • $ K^- \pi^0 K^0 $, (imode=4)
  • $ \pi^0 \pi^0 K^- $, (imode=5)
  • $ K^- \pi^- \pi^+ $, (imode=6)
  • $ \pi^- \bar{K}^0 \pi^0 $, (imode=7)
  • $ \pi^- \pi^0 \eta $, (imode=8)

    using the currents from TAUOLA

See also
ThreeMesonCurrentBase,
WeakDecayCurrent.
Defaulta1MatrixElement

Definition at line 48 of file ThreeMesonDefaultCurrent.h.

Member Function Documentation

◆ a1BreitWigner()

Complex Herwig::ThreeMesonDefaultCurrent::a1BreitWigner ( Energy2  q2) const
inlineprivate

$a_1$ Breit-Wigner

Parameters
q2The scale $q^2$ for the Breit-Wigner
Returns
The Breit-Wigner

Definition at line 325 of file ThreeMesonDefaultCurrent.h.

References sqrt().

◆ a1Width()

Energy Herwig::ThreeMesonDefaultCurrent::a1Width ( Energy2  q2) const
inlineprivate

The $a_1$ running width.

Parameters
q2The scale $q^2$ for the Breit-Wigner
Returns
The $a_1$ running width.

Definition at line 349 of file ThreeMesonDefaultCurrent.h.

References sqr(), and sqrt().

◆ acceptMode()

virtual bool Herwig::ThreeMesonDefaultCurrent::acceptMode ( int  ) const
protectedvirtual

Can the current handle a particular set of mesons.

As this current includes all the allowed modes this is always true.

Implements Herwig::ThreeMesonCurrentBase.

◆ BKstarF123()

Complex Herwig::ThreeMesonDefaultCurrent::BKstarF123 ( Energy2  q2,
int  ires 
) const
inlineprivate

The $K^*$ Breit-Wigner for the $F_{1,2,3}$ form factors.

Parameters
q2The scale $q^2$ for the Breit-Wigner
iresWhich $\rho$ multiplet
Returns
The Breit-Wigner

Definition at line 260 of file ThreeMesonDefaultCurrent.h.

◆ BKstarF5()

Complex Herwig::ThreeMesonDefaultCurrent::BKstarF5 ( Energy2  q2,
int  ires 
) const
inlineprivate

The $K^*$ Breit-Wigner for the $F_5$ form factors.

Parameters
q2The scale $q^2$ for the Breit-Wigner
iresWhich $\rho$ multiplet
Returns
The Breit-Wigner

Definition at line 285 of file ThreeMesonDefaultCurrent.h.

◆ BrhoF123()

Complex Herwig::ThreeMesonDefaultCurrent::BrhoF123 ( Energy2  q2,
int  ires 
) const
inlineprivate

The $\rho$ Breit-Wigner for the $F_{1,2,3}$ form factors.

Parameters
q2The scale $q^2$ for the Breit-Wigner
iresWhich $\rho$ multiplet
Returns
The Breit-Wigner

Definition at line 209 of file ThreeMesonDefaultCurrent.h.

◆ BrhoF5()

Complex Herwig::ThreeMesonDefaultCurrent::BrhoF5 ( Energy2  q2,
int  ires 
) const
inlineprivate

The $\rho$ Breit-Wigner for the $F_5$ form factors.

Parameters
q2The scale $q^2$ for the Breit-Wigner
iresWhich $\rho$ multiplet
Returns
The Breit-Wigner

Definition at line 235 of file ThreeMesonDefaultCurrent.h.

◆ calculateFormFactors()

virtual FormFactors Herwig::ThreeMesonDefaultCurrent::calculateFormFactors ( const int  ichan,
const int  imode,
Energy2  q2,
Energy2  s1,
Energy2  s2,
Energy2  s3 
) const
protectedvirtual

Calculate the form factor for the current.

Implements the form factors described above.

Parameters
ichanThe phase space channel
imodeThe mode
q2The scale $q^2$ for the current.
s1The invariant mass squared of particles 2 and 3, $s_1=m^2_{23}$.
s2The invariant mass squared of particles 1 and 3, $s_2=m^2_{13}$.
s3The invariant mass squared of particles 1 and 2, $s_3=m^2_{12}$.

Implements Herwig::ThreeMesonCurrentBase.

◆ clone()

virtual IBPtr Herwig::ThreeMesonDefaultCurrent::clone ( ) const
inlineprotectedvirtual

Make a simple clone of this object.

Returns
a pointer to the new object.

Implements ThePEG::InterfacedBase.

Definition at line 163 of file ThreeMesonDefaultCurrent.h.

◆ createMode()

virtual bool Herwig::ThreeMesonDefaultCurrent::createMode ( int  icharge,
unsigned int  imode,
DecayPhaseSpaceModePtr  mode,
unsigned int  iloc,
unsigned int  ires,
DecayPhaseSpaceChannelPtr  phase,
Energy  upp 
)
virtual

Complete the construction of the decay mode for integration.

This version addes the mesons for the current

Parameters
ichargeThe total charge of the outgoing particles in the current.
imodeThe mode in the current being asked for.
modeThe phase space mode for the integration
ilocThe location of the of the first particle from the current in the list of outgoing particles.
iresThe location of the first intermediate for the current.
phaseThe prototype phase space channel for the integration.
uppThe maximum possible mass the particles in the current are allowed to have.
Returns
Whether the current was sucessfully constructed.

Implements Herwig::WeakDecayCurrent.

◆ dataBaseOutput()

virtual void Herwig::ThreeMesonDefaultCurrent::dataBaseOutput ( ofstream &  os,
bool  header,
bool  create 
) const
virtual

Output the setup information for the particle database.

Parameters
osThe stream to output the information to
headerWhether or not to output the information for MySQL
createWhether or not to add a statement creating the object

Reimplemented from Herwig::ThreeMesonCurrentBase.

◆ doinit()

virtual void Herwig::ThreeMesonDefaultCurrent::doinit ( )
protectedvirtual

Initialize this object after the setup phase before saving and EventGenerator to disk.

Exceptions
InitExceptionif object could not be initialized properly.

Reimplemented from ThePEG::InterfacedBase.

◆ FKrho()

Complex Herwig::ThreeMesonDefaultCurrent::FKrho ( Energy2  si,
Energy2  sj,
int  ires 
) const
inlineprivate

Mixed Breit Wigner for the $F_5$ form factor.

Parameters
siThe scale $s_1$.
sjThe scale $s_2$.
iresWhich resonances to use
Returns
The mixed Breit-Wigner

Definition at line 311 of file ThreeMesonDefaultCurrent.h.

◆ fullclone()

virtual IBPtr Herwig::ThreeMesonDefaultCurrent::fullclone ( ) const
inlineprotectedvirtual

Make a clone of this object, possibly modifying the cloned object to make it sane.

Returns
a pointer to the new object.

Reimplemented from ThePEG::InterfacedBase.

Definition at line 169 of file ThreeMesonDefaultCurrent.h.

◆ inita1Width()

void Herwig::ThreeMesonDefaultCurrent::inita1Width ( int  iopt)
private

Initialize the $a_1$ running width.

Parameters
ioptInitialization option (-1 full calculation, 0 set up the interpolation)

◆ K1BreitWigner()

Complex Herwig::ThreeMesonDefaultCurrent::K1BreitWigner ( Energy2  q2) const
inlineprivate

The $K_1$ Breit-Wigner.

Parameters
q2The scale $q^2$ for the Breit-Wigner
Returns
The Breit-Wigner

Definition at line 337 of file ThreeMesonDefaultCurrent.h.

References sqr().

◆ persistentInput()

void Herwig::ThreeMesonDefaultCurrent::persistentInput ( PersistentIStream is,
int  version 
)

Function used to read in object persistently.

Parameters
isthe persistent input stream read from.
versionthe version number of the object when written.

◆ persistentOutput()

void Herwig::ThreeMesonDefaultCurrent::persistentOutput ( PersistentOStream os) const

Function used to write out object persistently.

Parameters
osthe persistent output stream written to.

◆ rhoKBreitWigner()

Complex Herwig::ThreeMesonDefaultCurrent::rhoKBreitWigner ( Energy2  q2,
unsigned int  itype,
unsigned int  ires 
) const
private

Breit-Wigners for the $\rho$ and $K^*$.

Parameters
q2The scale $q^2$ for the Breit-Wigner.
itypeThe type of Breit-Wigner, i.e. which masses and widths to use.x
iresWhich multiplet to use.

◆ threeBodyMatrixElement()

double Herwig::ThreeMesonDefaultCurrent::threeBodyMatrixElement ( const int  imode,
const Energy2  q2,
const Energy2  s3,
const Energy2  s2,
const Energy2  s1,
const Energy  m1,
const Energy  m2,
const Energy  m3 
) const

the matrix element for the $a_1$ decay to calculate the running width

Parameters
imodeThe mode for which the matrix element is needed.
q2The mass of the decaying off-shell $a_1$, $q^2$.
s3The invariant mass squared of particles 1 and 2, $s_3=m^2_{12}$.
s2The invariant mass squared of particles 1 and 3, $s_2=m^2_{13}$.
s1The invariant mass squared of particles 2 and 3, $s_1=m^2_{23}$.
m1The mass of the first outgoing particle.
m2The mass of the second outgoing particle.
m3The mass of the third outgoing particle.
Returns
The matrix element squared summed over spins.

The documentation for this class was generated from the following file: