evtgen is hosted by Hepforge, IPPP Durham
EvtGen  2.0.0
Monte Carlo generator of particle decays, in particular the weak decays of heavy flavour particles such as B mesons.
EvtTaulnunu.cpp
Go to the documentation of this file.
1 
2 /***********************************************************************
3 * Copyright 1998-2020 CERN for the benefit of the EvtGen authors *
4 * *
5 * This file is part of EvtGen. *
6 * *
7 * EvtGen is free software: you can redistribute it and/or modify *
8 * it under the terms of the GNU General Public License as published by *
9 * the Free Software Foundation, either version 3 of the License, or *
10 * (at your option) any later version. *
11 * *
12 * EvtGen is distributed in the hope that it will be useful, *
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of *
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
15 * GNU General Public License for more details. *
16 * *
17 * You should have received a copy of the GNU General Public License *
18 * along with EvtGen. If not, see <https://www.gnu.org/licenses/>. *
19 ***********************************************************************/
20 
22 
24 #include "EvtGenBase/EvtGenKine.hh"
25 #include "EvtGenBase/EvtPDL.hh"
27 #include "EvtGenBase/EvtPatches.hh"
28 #include "EvtGenBase/EvtReport.hh"
30 
31 #include <iostream>
32 #include <stdlib.h>
33 #include <string>
34 
35 std::string EvtTaulnunu::getName()
36 {
37  return "TAULNUNU";
38 }
39 
41 {
42  return new EvtTaulnunu;
43 }
44 
46 {
47  // check that there are 0 arguments
48  checkNArg( 0 );
49  checkNDaug( 3 );
50 
52 
56 }
57 
59 {
60  setProbMax( 650.0 );
61 }
62 
64 {
65  static EvtId TAUM = EvtPDL::getId( "tau-" );
66 
68 
69  EvtParticle *l, *nul, *nut;
70 
71  l = p->getDaug( 0 );
72  nul = p->getDaug( 1 );
73  nut = p->getDaug( 2 );
74 
75  EvtVector4C l1, l2, tau1, tau2;
76 
77  if ( p->getId() == TAUM ) {
78  tau1 = EvtLeptonVACurrent( nut->spParentNeutrino(), p->sp( 0 ) );
79  tau2 = EvtLeptonVACurrent( nut->spParentNeutrino(), p->sp( 1 ) );
80  l1 = EvtLeptonVACurrent( l->spParent( 0 ), nul->spParentNeutrino() );
81  l2 = EvtLeptonVACurrent( l->spParent( 1 ), nul->spParentNeutrino() );
82 
83  } else {
84  tau1 = EvtLeptonVACurrent( p->sp( 0 ), nut->spParentNeutrino() );
85  tau2 = EvtLeptonVACurrent( p->sp( 1 ), nut->spParentNeutrino() );
86  l1 = EvtLeptonVACurrent( nul->spParentNeutrino(), l->spParent( 0 ) );
87  l2 = EvtLeptonVACurrent( nul->spParentNeutrino(), l->spParent( 1 ) );
88  }
89 
90  vertex( 0, 0, tau1 * l1 );
91  vertex( 0, 1, tau1 * l2 );
92  vertex( 1, 0, tau2 * l1 );
93  vertex( 1, 1, tau2 * l2 );
94 
95  return;
96 }
void initProbMax() override
Definition: EvtTaulnunu.cpp:58
void decay(EvtParticle *p) override
Definition: EvtTaulnunu.cpp:63
virtual EvtDiracSpinor sp(int) const
std::string getName() override
Definition: EvtTaulnunu.cpp:35
EvtId getId() const
virtual EvtDiracSpinor spParent(int) const
EvtId * getDaugs()
Definition: EvtDecayBase.hh:66
void setProbMax(double prbmx)
Definition: EvtId.hh:27
void vertex(const EvtComplex &amp)
Definition: EvtDecayAmp.hh:37
double initializePhaseSpace(unsigned int numdaughter, EvtId *daughters, bool forceResetMasses=false, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
void checkNDaug(int d1, int d2=-1)
void checkSpinParent(EvtSpinType::spintype sp)
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
static EvtId getId(const std::string &name)
Definition: EvtPDL.cpp:287
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
virtual EvtDiracSpinor spParentNeutrino() const
int getNDaug() const
Definition: EvtDecayBase.hh:65
EvtVector4C EvtLeptonVACurrent(const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
EvtParticle * getDaug(int i)
Definition: EvtParticle.cpp:91
EvtDecayBase * clone() override
Definition: EvtTaulnunu.cpp:40
void init() override
Definition: EvtTaulnunu.cpp:45