Page MenuHomeHEPForge
Files F7877359

No OneTemporary
Actions

View Options

diff --git a/Models/Leptoquarks/LeptoquarkModelSLQFFVertex.cc b/Models/Leptoquarks/LeptoquarkModelSLQFFVertex.cc
--- a/Models/Leptoquarks/LeptoquarkModelSLQFFVertex.cc
+++ b/Models/Leptoquarks/LeptoquarkModelSLQFFVertex.cc
@@ -1,344 +1,344 @@
// -*- C++ -*-
//
// LeptoquarkModelSLQFFVertex.cc is a part of Herwig++ - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2011 The Herwig Collaboration
//
// Herwig++ is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
//
// This is the implementation of the non-inlined, non-templated member
// functions of the LeptoquarkModelSLQFFVertex class.
//
#include "LeptoquarkModelSLQFFVertex.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#include "ThePEG/Config/Constants.h"
using namespace Herwig;
IBPtr LeptoquarkModelSLQFFVertex::clone() const {
return new_ptr(*this);
}
IBPtr LeptoquarkModelSLQFFVertex::fullclone() const {
return new_ptr(*this);
}
LeptoquarkModelSLQFFVertex::LeptoquarkModelSLQFFVertex() {
orderInGem(1);
orderInGs(0);
}
void LeptoquarkModelSLQFFVertex::doinit() {
//S0
addToList( 15, 6,-9911561);
addToList(-15,-6, 9911561);
addToList(-16,-5, 9911561);
addToList( 16, 5,-9911561);
//~S0
addToList(-15,-5, 9921551);
addToList( 15, 5,-9921551);
//S1 triplet
//S1p
addToList(-15,-5, 9931551);
addToList( 15, 5,-9931551);
//S1z
addToList(-15,-6, 9931561);
addToList( 15, 6,-9931561);
addToList(-16,-5, 9931561);
addToList( 16, 5,-9931561);
//S1m
addToList(-16,-6, 9931661);
addToList( 16, 6,-9931661);
//S1/2 doublet
addToList( 15,-6, 9941561);
addToList(-15, 6,-9941561);
addToList(-15, 5,-9941551);
addToList(-16, 6,-9941551);
addToList( 15,-5, 9941551);
addToList( 16,-6, 9941551);
//S1/2 tilde doublet
addToList( 5,-16,-9951651);
addToList(-5, 16, 9951651);
addToList(-5, 15, 9951551);
addToList( 5,-15,-9951551);
//dS0
addToList( 15,-5, 9961551);
addToList(-15, 5,-9961551);
addToList( 16,-6, 9961551);
addToList(-16, 6,-9961551);
//~dS0
addToList( 15,-6, 9971561);
addToList(-15, 6,-9971561);
//dS1 triplet
//dS1p
addToList( 15,-6, 9981561);
addToList(-15, 6,-9981561);
//dS1z
addToList( 16,-6, 9981551);
addToList(-16, 6,-9981551);
addToList( 15,-5, 9981551);
addToList(-15, 5,-9981551);
//dS1m
addToList( 16,-5, 9981651);
addToList(-16, 5,-9981651);
//dS1/2 doublet
addToList(-15,-5, 9991551);
addToList( 15, 5,-9991551);
addToList(-15,-6, 9991561);
addToList( 15, 6,-9991561);
addToList(-16,-5, 9991561);
addToList( 16, 5,-9991561);
//dS1/2 tilde doublet
addToList(-15,-6, 9901561);
addToList( 15, 6,-9901561);
addToList(-16,-6, 9901661);
addToList( 16, 6,-9901661);
_theModel = generator()->standardModel();
tcHwLeptoquarkPtr hwLeptoquark=dynamic_ptr_cast<tcHwLeptoquarkPtr>(_theModel);
if(hwLeptoquark){
_CFF=hwLeptoquark->cfermion();
_cL0 =hwLeptoquark->cleft();
_cR0 =hwLeptoquark->cright();
_cR0t = hwLeptoquark->crighttilde();
_cL1 =hwLeptoquark->cleft1();
_cL12 =hwLeptoquark->cleft12();
_cR12 =hwLeptoquark->cright12();
_cL12t =hwLeptoquark->cleft12tilde();
_derivscale = hwLeptoquark->fscale();
_dcL0 =hwLeptoquark->dcleft();
_dcR0 =hwLeptoquark->dcright();
_dcR0t = hwLeptoquark->dcrighttilde();
_dcL1 =hwLeptoquark->dcleft1();
_dcL12 =hwLeptoquark->dcleft12();
_dcR12 =hwLeptoquark->dcright12();
_dcL12t =hwLeptoquark->dcleft12tilde();
}
FFSVertex::doinit();
}
void LeptoquarkModelSLQFFVertex::persistentOutput(PersistentOStream & os) const {
os << _CFF << _cL0 << _cR0 << _cR0t
<< _cL1 << _cL12 << _cR12 << _cL12t
<< _dcL0 << _dcR0 << _dcR0t
<< _dcL1 << _dcL12 << _dcR12 << _dcL12t
<< ounit(_derivscale,GeV);
}
void LeptoquarkModelSLQFFVertex::persistentInput(PersistentIStream & is, int) {
is >> _CFF >> _cL0 >> _cR0 >> _cR0t
>> _cL1 >> _cL12 >> _cR12 >> _cL12t
>>_dcL0 >> _dcR0 >> _dcR0t
>> _dcL1 >> _dcL12 >> _dcR12 >> _dcL12t
>> iunit(_derivscale,GeV);
}
ClassDescription<LeptoquarkModelSLQFFVertex>
LeptoquarkModelSLQFFVertex::initLeptoquarkModelSLQFFVertex;
// Definition of the static class description member.
void LeptoquarkModelSLQFFVertex::Init() {
static ClassDocumentation<LeptoquarkModelSLQFFVertex> documentation
("The LeptoquarkModelSLQFFVertex class is the implementation"
" of the helicity amplitude calculation of the Leptoquark"
" quark-lepton vertex.");
}
void LeptoquarkModelSLQFFVertex::setCoupling(Energy2,tcPDPtr aa ,tcPDPtr bb, tcPDPtr cc) {
long isc(cc->id()), ism(aa->id()),
ichg(bb->id());
long lqid = isc;
long smid_1 = ism;
long smid_2 = ichg;
- if(fabs(lqid) < 9900000) {
+ if(abs(lqid) < 9900000) {
lqid = ism;
smid_1 = ichg;
smid_2 = isc;
}
- if(fabs(lqid) < 9900000) {
+ if(abs(lqid) < 9900000) {
smid_1 = ism;
smid_2 = isc;
}
- if( fabs(smid_1) > fabs(smid_2) ) { swap(smid_1, smid_2); }
+ if( abs(smid_1) > abs(smid_2) ) { swap(smid_1, smid_2); }
const Energy denom = sqrt(2.) * _derivscale;
const double mtop = getParticleData(ParticleID::t)->mass() / denom;
const double mbot = getParticleData(ParticleID::b)->mass() / denom;
const double mtau = getParticleData(ParticleID::tauminus)->mass() / denom;
//set the couplings to left and right
//S0
- if( fabs(isc) == 9911561 || fabs(ism) == 9911561 || fabs(ichg) == 9911561 ) {
- if(fabs(isc) == 5 || fabs(ism) == 5 || fabs(ichg) == 5) {
+ if( abs(isc) == 9911561 || abs(ism) == 9911561 || abs(ichg) == 9911561 ) {
+ if(abs(isc) == 5 || abs(ism) == 5 || abs(ichg) == 5) {
_cL = -_cL0; _cR = Complex(0.);
}
- if(fabs(isc) == 6 || fabs(ism) == 6 || fabs(ichg) == 6) {
+ if(abs(isc) == 6 || abs(ism) == 6 || abs(ichg) == 6) {
_cL = _cL0;
_cR = _cR0;
}
}
//~S0
- if( fabs(isc) == 9921551 || fabs(ism) == 9921551 || fabs(ichg) == 9921551 ) {
+ if( abs(isc) == 9921551 || abs(ism) == 9921551 || abs(ichg) == 9921551 ) {
_cL = Complex(0.); _cR = _cR0t;
}
//S1 triplet
//Q = + 4/3
- if( fabs(isc) == 9931551 || fabs(ism) == 9931551 || fabs(ichg) == 9931551 ) {
+ if( abs(isc) == 9931551 || abs(ism) == 9931551 || abs(ichg) == 9931551 ) {
_cL = sqrt(2.)* _cL1; _cR = Complex(0.);
}
//Q = + 1/3
- if( fabs(isc) == 9931561 || fabs(ism) == 9931561 || fabs(ichg) == 9931561 ) {
+ if( abs(isc) == 9931561 || abs(ism) == 9931561 || abs(ichg) == 9931561 ) {
_cL = - _cL1; _cR = Complex(0.);
}
//Q = - 2/3
- if( fabs(isc) == 9931661 || fabs(ism) == 9931661 || fabs(ichg) == 9931661 ) {
+ if( abs(isc) == 9931661 || abs(ism) == 9931661 || abs(ichg) == 9931661 ) {
_cL = sqrt(2.) * _cL1; _cR = Complex(0.);
}
//S1/2 doublet
//Q = + 5/3
- if( fabs(isc) == 9941561 || fabs(ism) == 9941561 || fabs(ichg) == 9941561 ) {
+ if( abs(isc) == 9941561 || abs(ism) == 9941561 || abs(ichg) == 9941561 ) {
_cR = _cL12; _cL = _cR12;
}
//Q = + 2/3
- if( fabs(isc) == 9941551 || fabs(ism) == 9941551 || fabs(ichg) == 9941551 ) {
- if(fabs(isc) == 5 || fabs(ism) == 5 || fabs(ichg) == 5) {
+ if( abs(isc) == 9941551 || abs(ism) == 9941551 || abs(ichg) == 9941551 ) {
+ if(abs(isc) == 5 || abs(ism) == 5 || abs(ichg) == 5) {
_cR = Complex(0.); _cL = - _cR12;
}
- if(fabs(isc) == 6 || fabs(ism) == 6 || fabs(ichg) == 6) {
+ if(abs(isc) == 6 || abs(ism) == 6 || abs(ichg) == 6) {
_cL = Complex(0.); _cR = _cL12;
}
}
//S1/2 tilde doublet
//Q = + 2/3
- if( fabs(isc) == 9951551 || fabs(ism) == 9951551 || fabs(ichg) == 9951551 ) {
+ if( abs(isc) == 9951551 || abs(ism) == 9951551 || abs(ichg) == 9951551 ) {
_cR = _cL12t; _cL = Complex(0.);
}
//Q = - 1/3
- if( fabs(isc) == 9951651 || fabs(ism) == 9951651 || fabs(ichg) == 9951651 ) {
+ if( abs(isc) == 9951651 || abs(ism) == 9951651 || abs(ichg) == 9951651 ) {
_cR = _cL12t; _cL = Complex(0.);
}
//dS0
- if( fabs(isc) == 9961551 || fabs(ism) == 9961551 || fabs(ichg) == 9961551) {
- if(fabs(isc) == 5 || fabs(ism) == 5 || fabs(ichg) == 5) {
+ if( abs(isc) == 9961551 || abs(ism) == 9961551 || abs(ichg) == 9961551) {
+ if(abs(isc) == 5 || abs(ism) == 5 || abs(ichg) == 5) {
_cR = _dcL0 * mbot +_dcR0 * mtau;
_cL = _dcR0 * mbot + _dcL0 * mtau;
}
- if(fabs(isc) == 6 || fabs(ism) == 6 || fabs(ichg) == 6) {
+ if(abs(isc) == 6 || abs(ism) == 6 || abs(ichg) == 6) {
_cR = _dcL0 * mtop;
_cL = Complex(0.);
}
}
//d~S0
- if( fabs(isc) == 9971561 || fabs(ism) == 9971561 || fabs(ichg) == 9971561) {
+ if( abs(isc) == 9971561 || abs(ism) == 9971561 || abs(ichg) == 9971561) {
_cR = _dcR0t * mtau;
_cL = _dcR0t * mtop;
}
//dS1 triplet
- if( fabs(isc) == 9981561 || fabs(ism) == 9981561 || fabs(ichg) == 9981561) {
+ if( abs(isc) == 9981561 || abs(ism) == 9981561 || abs(ichg) == 9981561) {
_cR = sqrt(2.) * _dcL1 * mtop;
_cL = sqrt(2.) * _dcL1 * mtau;
}
- if( fabs(isc) == 9981551 || fabs(ism) == 9981551 || fabs(ichg) == 9981551) {
- if(fabs(isc) == 5 || fabs(ism) == 5 || fabs(ichg) == 5) {
+ if( abs(isc) == 9981551 || abs(ism) == 9981551 || abs(ichg) == 9981551) {
+ if(abs(isc) == 5 || abs(ism) == 5 || abs(ichg) == 5) {
_cR = -_dcL1 * mbot;
_cL = -_dcL1 * mtau;
}
- if(fabs(isc) == 6 || fabs(ism) == 6 || fabs(ichg) == 6) {
+ if(abs(isc) == 6 || abs(ism) == 6 || abs(ichg) == 6) {
_cR = _dcL1 * mtop;
_cL = Complex(0.);
}
}
- if( fabs(isc) == 9981651 || fabs(ism) == 9981651 || fabs(ichg) == 9981651) {
+ if( abs(isc) == 9981651 || abs(ism) == 9981651 || abs(ichg) == 9981651) {
_cL = sqrt(2.) * _dcL1 * mbot;
_cR = Complex(0.);
}
//dS1/2 doublet
- if( fabs(isc) == 9991551 || fabs(ism) == 9991551 || fabs(ichg) == 9991551 ) {
+ if( abs(isc) == 9991551 || abs(ism) == 9991551 || abs(ichg) == 9991551 ) {
_cL = _dcL12 * mbot + _dcR12 * mtau;
_cR = _dcR12 * mbot + _dcL12 * mtau;
}
- if( fabs(isc) == 9991561 || fabs(ism) == 9991561 || fabs(ichg) == 9991561 ) {
- if(fabs(isc) == 6 || fabs(ism) == 6 || fabs(ichg) == 6) {
+ if( abs(isc) == 9991561 || abs(ism) == 9991561 || abs(ichg) == 9991561 ) {
+ if(abs(isc) == 6 || abs(ism) == 6 || abs(ichg) == 6) {
_cL = _dcR12 * mtau;
_cR = _dcR12 * mtop;
}
- if(fabs(isc) == 5 || fabs(ism) == 5 || fabs(ichg) == 5) {
+ if(abs(isc) == 5 || abs(ism) == 5 || abs(ichg) == 5) {
_cL = _dcL12 * mbot;
}
}
//dS1/2 tilde doublet
- if( fabs(isc) == 9901561 || fabs(ism) == 9901561 || fabs(ichg) == 9901561 ) {
+ if( abs(isc) == 9901561 || abs(ism) == 9901561 || abs(ichg) == 9901561 ) {
_cL = _dcL12t * mtop;
_cR = _dcL12t * mtau;
}
- if( fabs(isc) == 9901661 || fabs(ism) == 9901661 || fabs(ichg) == 9901661 ) {
+ if( abs(isc) == 9901661 || abs(ism) == 9901661 || abs(ichg) == 9901661 ) {
_cL = _dcL12t * mtop;
_cR = Complex(0.);
}
if(smid_1 > 0) {
left(conj(_cR));
right(conj(_cL));
}
else {
left(_cL);
right(_cR);
}
norm(_CFF);
}
diff --git a/Models/TTbAsymm/TTbAModelAGQQVertex.cc b/Models/TTbAsymm/TTbAModelAGQQVertex.cc
--- a/Models/TTbAsymm/TTbAModelAGQQVertex.cc
+++ b/Models/TTbAsymm/TTbAModelAGQQVertex.cc
@@ -1,101 +1,101 @@
// -*- C++ -*-
//
// TTbAModelAGQQVertex.cc is a part of Herwig++ - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2011 The Herwig Collaboration
//
// Herwig++ is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
//
// This is the implementation of the non-inlined, non-templated member
// functions of the TTbAModelAGQQVertex class.
//
#include "TTbAModelAGQQVertex.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#include "ThePEG/Config/Constants.h"
using namespace Herwig;
IBPtr TTbAModelAGQQVertex::clone() const {
return new_ptr(*this);
}
IBPtr TTbAModelAGQQVertex::fullclone() const {
return new_ptr(*this);
}
TTbAModelAGQQVertex::TTbAModelAGQQVertex() {
orderInGem(1);
orderInGs(1);
addToList(-1,1,63);
addToList(-2,2,63);
addToList(-3,3,63);
addToList(-4,4,63);
addToList(-5,5,63);
addToList(-6,6,63);
}
void TTbAModelAGQQVertex::doinit() {
_theModel = generator()->standardModel();
tcHwTTbAPtr hwTTbA=dynamic_ptr_cast<tcHwTTbAPtr>(_theModel);
if(hwTTbA) {
_cAGQQ_R =hwTTbA->_cAGQQ_right();
_cAGQQ_L =hwTTbA->_cAGQQ_left();
_cAGTT_R =hwTTbA->_cAGTT_right();
_cAGTT_L =hwTTbA->_cAGTT_left();
_models = hwTTbA->_model();
}
FFVVertex::doinit();
}
void TTbAModelAGQQVertex::persistentOutput(PersistentOStream & os) const {
os << _cAGQQ_R << _cAGQQ_L << _cAGTT_R << _cAGTT_L << _models;
}
void TTbAModelAGQQVertex::persistentInput(PersistentIStream & is, int) {
is >> _cAGQQ_R >> _cAGQQ_L >>_cAGTT_R >> _cAGTT_L >> _models;
}
ClassDescription<TTbAModelAGQQVertex>
TTbAModelAGQQVertex::initTTbAModelAGQQVertex;
// Definition of the static class description member.
void TTbAModelAGQQVertex::Init() {
static ClassDocumentation<TTbAModelAGQQVertex> documentation
("The TTbAModelAGQQVertex class is the implementation"
" of the helicity amplitude calculation of the TTbA"
" quark-lepton vertex.");
}
void TTbAModelAGQQVertex::setCoupling(Energy2 q2,tcPDPtr aa ,tcPDPtr bb, tcPDPtr cc) {
double _cL = 0, _cR = 0;
double gstrong = 1.0;
gstrong = strongCoupling(q2);
- if(fabs(aa->id()) == 63 || fabs(bb->id()) == 63 || fabs(cc->id()) == 63) {
- if(fabs(aa->id()) !=6 && fabs(bb->id()) !=6 && fabs(cc->id()) != 6) {
+ if(abs(aa->id()) == 63 || abs(bb->id()) == 63 || abs(cc->id()) == 63) {
+ if(abs(aa->id()) !=6 && abs(bb->id()) !=6 && abs(cc->id()) != 6) {
_cR = _cAGQQ_R;
_cL = _cAGQQ_L;
} else {
_cR = _cAGTT_R;
_cL = _cAGTT_L;
}
}
if(_models!=2) { _cL = 1E-10; _cR = 1E-10; }
left(_cL);
right(_cR);
norm(gstrong);
}
diff --git a/Models/TTbAsymm/TTbAModelSU2XVertex.cc b/Models/TTbAsymm/TTbAModelSU2XVertex.cc
--- a/Models/TTbAsymm/TTbAModelSU2XVertex.cc
+++ b/Models/TTbAsymm/TTbAModelSU2XVertex.cc
@@ -1,184 +1,184 @@
// -*- C++ -*-
//
// TTbAModelSU2XVertex.cc is a part of Herwig++ - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2011 The Herwig Collaboration
//
// Herwig++ is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
//
// This is the implementation of the non-inlined, non-templated member
// functions of the TTbAModelSU2XVertex class.
//
#include "TTbAModelSU2XVertex.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#include "ThePEG/Config/Constants.h"
using namespace Herwig;
IBPtr TTbAModelSU2XVertex::clone() const {
return new_ptr(*this);
}
IBPtr TTbAModelSU2XVertex::fullclone() const {
return new_ptr(*this);
}
TTbAModelSU2XVertex::TTbAModelSU2XVertex() {
orderInGem(1);
orderInGs(1);
addToList(-6,6,70);
addToList(-2,2,70);
addToList(-2,6,70);
addToList(2,-6,70);
addToList(-6,6,71);
addToList(-2,2,71);
addToList(-2,6,71);
addToList(2,-6,71);
addToList(-6,6,-71);
addToList(-2,2,-71);
addToList(-2,6,-71);
addToList(2,-6,-71);
}
void TTbAModelSU2XVertex::doinit() {
_theModel = generator()->standardModel();
tcHwTTbAPtr hwTTbA=dynamic_ptr_cast<tcHwTTbAPtr>(_theModel);
if(hwTTbA) {
_alphaX =hwTTbA->_alphaX_value();
_costhetaX =hwTTbA->_costhetaX_value();
_models =hwTTbA->_model();
}
FFVVertex::doinit();
}
void TTbAModelSU2XVertex::persistentOutput(PersistentOStream & os) const {
os << _alphaX << _costhetaX << _models;
}
void TTbAModelSU2XVertex::persistentInput(PersistentIStream & is, int) {
is >> _alphaX >> _costhetaX >> _models;
}
ClassDescription<TTbAModelSU2XVertex>
TTbAModelSU2XVertex::initTTbAModelSU2XVertex;
// Definition of the static class description member.
void TTbAModelSU2XVertex::Init() {
static ClassDocumentation<TTbAModelSU2XVertex> documentation
("The TTbAModelSU2XVertex class is the implementation"
" of the helicity amplitude calculation of the TTbA"
" SU(2)_X vertex.");
}
void TTbAModelSU2XVertex::setCoupling(Energy2,tcPDPtr aa ,tcPDPtr bb, tcPDPtr cc) {
double _cR = 0, _fac = 0;
_gX = sqrt( 4 * Constants::pi * _alphaX );
double ct = _costhetaX;
double st = sqrt(1 - pow(ct,2));
//Vz
- if( fabs(aa->id()) == 70 || fabs(bb->id()) == 70 || fabs(cc->id()) == 70) {
+ if( abs(aa->id()) == 70 || abs(bb->id()) == 70 || abs(cc->id()) == 70) {
_fac = _gX / 2.0;
if( aa->id() == 6 || bb->id() == 6 || cc->id() == 6) {
if( aa->id() == -6 || bb->id() == -6 || cc->id() == -6) {
_cR = pow(ct,2) - pow(st,2);
}
}
if( aa->id() == 2 || bb->id() == 2 || cc->id() == 2) {
if( aa->id() == -2 || bb->id() == -2 || cc->id() == -2) {
_cR = pow(st,2) - pow(ct,2);
}
}
- if( fabs(aa->id()) == 2 || fabs(bb->id()) == 2 || fabs(cc->id()) == 2) {
- if( fabs(aa->id()) == 6 || fabs(bb->id()) == 6 || fabs(cc->id()) == 6) {
+ if( abs(aa->id()) == 2 || abs(bb->id()) == 2 || abs(cc->id()) == 2) {
+ if( abs(aa->id()) == 6 || abs(bb->id()) == 6 || abs(cc->id()) == 6) {
_cR = 2 * ct * st;
}
}
}
//Ym
if( aa->id() == -71 || bb->id() == -71 || cc->id() == -71 ) {
_fac = _gX / sqrt(2.0);
if( aa->id() == 6 || bb->id() == 6 || cc->id() == 6) {
if( aa->id() == -6 || bb->id() == -6 || cc->id() == -6) {
_cR = - ct * st;
}
}
if( aa->id() == 2 || bb->id() == 2 || cc->id() == 2) {
if( aa->id() == -2 || bb->id() == -2 || cc->id() == -2) {
_cR = ct * st;
}
}
if( aa->id() == 2 || bb->id() == 2 || cc->id() == 2) {
if( aa->id() == -6 || bb->id() == -6 || cc->id() == -6) {
_cR = - pow(st,2);
}
}
if( aa->id() == -2 || bb->id() == -2 || cc->id() == -2) {
if( aa->id() == 6 || bb->id() == 6 || cc->id() == 6) {
_cR = pow(ct,2);
}
}
}
//Yp
if( aa->id() == 71 || bb->id() == 71 || cc->id() == 71 ) {
_fac = _gX / sqrt(2.0);
if( aa->id() == 6 || bb->id() == 6 || cc->id() == 6) {
if( aa->id() == -6 || bb->id() == -6 || cc->id() == -6) {
_cR = - ct * st;
}
}
if( aa->id() == 2 || bb->id() == 2 || cc->id() == 2) {
if( aa->id() == -2 || bb->id() == -2 || cc->id() == -2) {
_cR = ct * st;
}
}
if( aa->id() == 2 || bb->id() == 2 || cc->id() == 2) {
if( aa->id() == -6 || bb->id() == -6 || cc->id() == -6) {
_cR = pow(ct,2);
}
}
if( aa->id() == -2 || bb->id() == -2 || cc->id() == -2) {
if( aa->id() == 6 || bb->id() == 6 || cc->id() == 6) {
_cR = - pow(st,2);
}
}
}
//normalise according to Lagrangian factor
_cR *= _fac;
//If this model is not selected set coupling to zero.
if(_models!=3) { _cR = 1E-10; }
right(_cR);
left(0.);
norm(1.0);
}
diff --git a/Models/TTbAsymm/TTbAModelWPTDVertex.cc b/Models/TTbAsymm/TTbAModelWPTDVertex.cc
--- a/Models/TTbAsymm/TTbAModelWPTDVertex.cc
+++ b/Models/TTbAsymm/TTbAModelWPTDVertex.cc
@@ -1,83 +1,83 @@
// -*- C++ -*-
//
// TTbAModelWPTDVertex.cc is a part of Herwig++ - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2011 The Herwig Collaboration
//
// Herwig++ is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
//
// This is the implementation of the non-inlined, non-templated member
// functions of the TTbAModelWPTDVertex class.
//
#include "TTbAModelWPTDVertex.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#include "ThePEG/Config/Constants.h"
using namespace Herwig;
IBPtr TTbAModelWPTDVertex::clone() const {
return new_ptr(*this);
}
IBPtr TTbAModelWPTDVertex::fullclone() const {
return new_ptr(*this);
}
TTbAModelWPTDVertex::TTbAModelWPTDVertex() {
addToList(-1,6,-34);
addToList(-6,1,34);
orderInGem(1);
orderInGs(1);
}
void TTbAModelWPTDVertex::doinit() {
_theModel = generator()->standardModel();
tcHwTTbAPtr hwTTbA=dynamic_ptr_cast<tcHwTTbAPtr>(_theModel);
if(hwTTbA) {
_cWPTD_R =hwTTbA->_cWPTD_right();
_cWPTD_L =hwTTbA->_cWPTD_left();
_models =hwTTbA->_model();
}
FFVVertex::doinit();
}
void TTbAModelWPTDVertex::persistentOutput(PersistentOStream & os) const {
os << _cWPTD_R << _cWPTD_L << _models;
}
void TTbAModelWPTDVertex::persistentInput(PersistentIStream & is, int) {
is >> _cWPTD_R >> _cWPTD_L >> _models;
}
ClassDescription<TTbAModelWPTDVertex>
TTbAModelWPTDVertex::initTTbAModelWPTDVertex;
// Definition of the static class description member.
void TTbAModelWPTDVertex::Init() {
static ClassDocumentation<TTbAModelWPTDVertex> documentation
("The TTbAModelWPTDVertex class is the implementation"
" of the helicity amplitude calculation of the TTbA"
" quark-lepton vertex.");
}
void TTbAModelWPTDVertex::setCoupling(Energy2,tcPDPtr aa ,tcPDPtr bb, tcPDPtr cc) {
double _cL = 0, _cR = 0;
- if(fabs(aa->id()) == 34 || fabs(bb->id()) == 34 || fabs(cc->id()) == 34) {
+ if(abs(aa->id()) == 34 || abs(bb->id()) == 34 || abs(cc->id()) == 34) {
_cR = _cWPTD_R;
_cL = _cWPTD_L;
}
if(_models!=0) { _cL = 1E-10; _cR = 1E-10; }
left(_cL);
right(_cR);
norm(1.0);
}
diff --git a/Models/TTbAsymm/TTbAModelZPQQVertex.cc b/Models/TTbAsymm/TTbAModelZPQQVertex.cc
--- a/Models/TTbAsymm/TTbAModelZPQQVertex.cc
+++ b/Models/TTbAsymm/TTbAModelZPQQVertex.cc
@@ -1,99 +1,99 @@
// -*- C++ -*-
//
// TTbAModelZPQQVertex.cc is a part of Herwig++ - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2011 The Herwig Collaboration
//
// Herwig++ is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
//
// This is the implementation of the non-inlined, non-templated member
// functions of the TTbAModelZPQQVertex class.
//
#include "TTbAModelZPQQVertex.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#include "ThePEG/Config/Constants.h"
using namespace Herwig;
IBPtr TTbAModelZPQQVertex::clone() const {
return new_ptr(*this);
}
IBPtr TTbAModelZPQQVertex::fullclone() const {
return new_ptr(*this);
}
TTbAModelZPQQVertex::TTbAModelZPQQVertex() {
addToList(-2,6,32);
addToList(-6,2,32);
addToList(-2,2,32);
addToList(-4,4,32);
orderInGem(1);
orderInGs(1);
}
void TTbAModelZPQQVertex::doinit() {
_theModel = generator()->standardModel();
tcHwTTbAPtr hwTTbA=dynamic_ptr_cast<tcHwTTbAPtr>(_theModel);
if(hwTTbA) {
_cZPTU_R =hwTTbA->_cZPTU_right();
_cZPTU_L =hwTTbA->_cZPTU_left();
_cZPUU_R =hwTTbA->_cZPUU_right();
_cZPUU_L =hwTTbA->_cZPUU_left();
_cZPCC_R =hwTTbA->_cZPCC_right();
_cZPCC_L =hwTTbA->_cZPCC_left();
_models =hwTTbA->_model();
}
FFVVertex::doinit();
}
void TTbAModelZPQQVertex::persistentOutput(PersistentOStream & os) const {
os << _cZPTU_R << _cZPTU_L << _cZPUU_R << _cZPUU_L << _cZPCC_R << _cZPCC_L << _models;
}
void TTbAModelZPQQVertex::persistentInput(PersistentIStream & is, int) {
is >> _cZPTU_R >> _cZPTU_L >> _cZPUU_R >> _cZPUU_L >> _cZPCC_R >> _cZPCC_L >> _models;
}
ClassDescription<TTbAModelZPQQVertex>
TTbAModelZPQQVertex::initTTbAModelZPQQVertex;
// Definition of the static class description member.
void TTbAModelZPQQVertex::Init() {
static ClassDocumentation<TTbAModelZPQQVertex> documentation
("The TTbAModelZPQQVertex class is the implementation"
" of the helicity amplitude calculation of the TTbA"
" Z prime Quark-antiQuark vertex.");
}
void TTbAModelZPQQVertex::setCoupling(Energy2,tcPDPtr aa ,tcPDPtr bb, tcPDPtr cc) {
double _cR = 0, _cL = 0;
- if( fabs(aa->id()) == 6 || fabs(bb->id()) == 6 || fabs(cc->id()) == 6) {
+ if( abs(aa->id()) == 6 || abs(bb->id()) == 6 || abs(cc->id()) == 6) {
_cR = _cZPTU_R;
_cL = _cZPTU_L;
} else {
- if( fabs(aa->id()) != 4 && fabs(bb->id()) != 4 && fabs(cc->id()) != 4) {
+ if( abs(aa->id()) != 4 && abs(bb->id()) != 4 && abs(cc->id()) != 4) {
_cR = _cZPUU_R;
_cL = _cZPUU_L;
}
- if( fabs(aa->id()) == 4 || fabs(bb->id()) == 4 || fabs(cc->id()) == 4) {
+ if( abs(aa->id()) == 4 || abs(bb->id()) == 4 || abs(cc->id()) == 4) {
_cR = _cZPCC_R;
_cL = _cZPCC_L;
}
}
if(_models!=1) { _cL = 1E-10; _cR = 1E-10; }
right(_cR);
left(_cL);
norm(1.0);
}
diff --git a/Models/Zprime/ZprimeModelZPQQVertex.cc b/Models/Zprime/ZprimeModelZPQQVertex.cc
--- a/Models/Zprime/ZprimeModelZPQQVertex.cc
+++ b/Models/Zprime/ZprimeModelZPQQVertex.cc
@@ -1,184 +1,184 @@
// -*- C++ -*-
//
// ZprimeModelZPQQVertex.cc is a part of Herwig++ - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2011 The Herwig Collaboration
//
// Herwig++ is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
//
// This is the implementation of the non-inlined, non-templated member
// functions of the ZprimeModelZPQQVertex class.
//
#include "ZprimeModelZPQQVertex.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#include "ThePEG/Config/Constants.h"
using namespace Herwig;
IBPtr ZprimeModelZPQQVertex::clone() const {
return new_ptr(*this);
}
IBPtr ZprimeModelZPQQVertex::fullclone() const {
return new_ptr(*this);
}
ZprimeModelZPQQVertex::ZprimeModelZPQQVertex() {
addToList(-2,6,32);
addToList(-6,2,32);
addToList(-6,6,32);
addToList(-5,5,32);
addToList(-4,4,32);
addToList(-3,3,32);
addToList(-2,2,32);
addToList(-1,1,32);
addToList(-11,11,32);
addToList(-13,13,32);
addToList(-15,15,32);
addToList(-12,12,32);
addToList(-14,14,32);
addToList(-16,16,32);
orderInGem(1);
orderInGs(0);
}
void ZprimeModelZPQQVertex::doinit() {
_theModel = generator()->standardModel();
tcHwZprimePtr hwZprime=dynamic_ptr_cast<tcHwZprimePtr>(_theModel);
if(hwZprime) {
_cZPTU_R =hwZprime->_cZPTU_right();
_cZPTU_L =hwZprime->_cZPTU_left();
_cZPTT_R =hwZprime->_cZPTT_right();
_cZPTT_L =hwZprime->_cZPTT_left();
_cZPUU_R =hwZprime->_cZPUU_right();
_cZPUU_L =hwZprime->_cZPUU_left();
_cZPCC_R =hwZprime->_cZPCC_right();
_cZPCC_L =hwZprime->_cZPCC_left();
_cZPDD_R =hwZprime->_cZPDD_right();
_cZPDD_L =hwZprime->_cZPDD_left();
_cZPBB_R =hwZprime->_cZPBB_right();
_cZPBB_L =hwZprime->_cZPBB_left();
_cZPSS_R =hwZprime->_cZPSS_right();
_cZPSS_L =hwZprime->_cZPSS_left();
_cZPee_R =hwZprime->_cZPee_right();
_cZPee_L =hwZprime->_cZPee_left();
_cZPmm_R =hwZprime->_cZPmm_right();
_cZPmm_L =hwZprime->_cZPmm_left();
_cZPtt_R =hwZprime->_cZPtt_right();
_cZPtt_L =hwZprime->_cZPtt_left();
_cZPnuenue_R =hwZprime->_cZPnuenue_right();
_cZPnuenue_L =hwZprime->_cZPnuenue_left();
_cZPnuenue_R =hwZprime->_cZPnumnum_right();
_cZPnumnum_L =hwZprime->_cZPnumnum_left();
_cZPnutnut_R =hwZprime->_cZPnutnut_right();
_cZPnutnut_L =hwZprime->_cZPnutnut_left();
_cZP_o =hwZprime->_cZPoverallCoup();
}
FFVVertex::doinit();
}
void ZprimeModelZPQQVertex::persistentOutput(PersistentOStream & os) const {
os << _cZPTU_R << _cZPTU_L << _cZPTT_R << _cZPTT_L << _cZPUU_R << _cZPUU_L << _cZPCC_R << _cZPCC_L << _cZPDD_R << _cZPDD_L << _cZPSS_R << _cZPSS_L << _cZPBB_R << _cZPBB_L << _cZPee_R << _cZPee_L << _cZPmm_R << _cZPmm_L << _cZPtt_R << _cZPtt_L << _cZPnuenue_R << _cZPnuenue_L << _cZPnumnum_R << _cZPnumnum_L << _cZPnutnut_R << _cZPnutnut_L << _cZP_o;
}
void ZprimeModelZPQQVertex::persistentInput(PersistentIStream & is, int) {
is >> _cZPTU_R >> _cZPTU_L >> _cZPTT_R >> _cZPTT_L >> _cZPUU_R >> _cZPUU_L >> _cZPCC_R >> _cZPCC_L >> _cZPDD_R >> _cZPDD_L >> _cZPSS_R >> _cZPSS_L >> _cZPBB_R >> _cZPBB_L >> _cZPee_R >> _cZPee_L >> _cZPmm_R >> _cZPmm_L >> _cZPtt_R >> _cZPtt_L >> _cZPnuenue_R >> _cZPnuenue_L >> _cZPnumnum_R >> _cZPnumnum_L >> _cZPnutnut_R >> _cZPnutnut_L >> _cZP_o;
}
ClassDescription<ZprimeModelZPQQVertex>
ZprimeModelZPQQVertex::initZprimeModelZPQQVertex;
// Definition of the static class description member.
void ZprimeModelZPQQVertex::Init() {
static ClassDocumentation<ZprimeModelZPQQVertex> documentation
("The ZprimeModelZPQQVertex class is the implementation"
" of the helicity amplitude calculation of the Zprime"
" Z prime Quark-antiQuark vertex.");
}
void ZprimeModelZPQQVertex::setCoupling(Energy2,tcPDPtr aa ,tcPDPtr bb, tcPDPtr cc) {
double _cR = 1.0, _cL = 1.0;
long ccc(cc->id()), aaa(aa->id()), bbb(bb->id());
- if( fabs(aaa) == 6 || fabs(bbb) == 6 || fabs(ccc) == 6 ) {
- if( fabs(aaa) !=2 && fabs(bbb) !=2 && fabs(ccc) != 2 ) {
+ if( abs(aaa) == 6 || abs(bbb) == 6 || abs(ccc) == 6 ) {
+ if( abs(aaa) !=2 && abs(bbb) !=2 && abs(ccc) != 2 ) {
_cL = _cZPTT_L; _cR = _cZPTT_R;
- } else if( fabs(aaa) ==2 || fabs(bbb) ==2 || fabs(ccc) == 2 ) {
+ } else if( abs(aaa) ==2 || abs(bbb) ==2 || abs(ccc) == 2 ) {
_cL = _cZPTU_L; _cR = _cZPTU_R;
}
}
- if( fabs(aaa) == 5 || fabs(bbb) == 5 || fabs(ccc) == 5 ) {
+ if( abs(aaa) == 5 || abs(bbb) == 5 || abs(ccc) == 5 ) {
_cL = _cZPBB_L; _cR = _cZPBB_R;
}
- if( fabs(aaa) == 4 || fabs(bbb) == 4 || fabs(ccc) == 4 ) {
+ if( abs(aaa) == 4 || abs(bbb) == 4 || abs(ccc) == 4 ) {
_cL = _cZPCC_L; _cR = _cZPCC_R;
}
- if( fabs(aaa) == 3 || fabs(bbb) == 3 || fabs(ccc) == 3 ) {
+ if( abs(aaa) == 3 || abs(bbb) == 3 || abs(ccc) == 3 ) {
_cL = _cZPSS_L; _cR = _cZPSS_R;
}
- if( (fabs(aaa) == 2 || fabs(bbb) == 2 || fabs(ccc) == 2) && (fabs(aaa) !=6 && fabs(bbb) !=6 && fabs(ccc) != 6)) {
+ if( (abs(aaa) == 2 || abs(bbb) == 2 || abs(ccc) == 2) && (abs(aaa) !=6 && abs(bbb) !=6 && abs(ccc) != 6)) {
_cL = _cZPUU_L; _cR = _cZPUU_R;
}
- if( fabs(aaa) == 1 || fabs(bbb) == 1 || fabs(ccc) == 1 ) {
+ if( abs(aaa) == 1 || abs(bbb) == 1 || abs(ccc) == 1 ) {
_cL = _cZPDD_L; _cR = _cZPDD_R;
}
- if( fabs(aaa) == 11 || fabs(bbb) == 11 || fabs(ccc) == 11 ) {
+ if( abs(aaa) == 11 || abs(bbb) == 11 || abs(ccc) == 11 ) {
_cL = _cZPee_L; _cR = _cZPee_R;
}
- if( fabs(aaa) == 13 || fabs(bbb) == 13 || fabs(ccc) == 13 ) {
+ if( abs(aaa) == 13 || abs(bbb) == 13 || abs(ccc) == 13 ) {
_cL = _cZPmm_L; _cR = _cZPmm_R;
}
- if( fabs(aaa) == 15 || fabs(bbb) == 15 || fabs(ccc) == 15 ) {
+ if( abs(aaa) == 15 || abs(bbb) == 15 || abs(ccc) == 15 ) {
_cL = _cZPtt_L; _cR = _cZPtt_R;
}
- if( fabs(aaa) == 12 || fabs(bbb) == 12 || fabs(ccc) == 12 ) {
+ if( abs(aaa) == 12 || abs(bbb) == 12 || abs(ccc) == 12 ) {
_cL = _cZPnuenue_L; _cR = _cZPnuenue_R;
}
- if( fabs(aaa) == 14 || fabs(bbb) == 14 || fabs(ccc) == 14 ) {
+ if( abs(aaa) == 14 || abs(bbb) == 14 || abs(ccc) == 14 ) {
_cL = _cZPnumnum_L; _cR = _cZPnumnum_R;
}
- if( fabs(aaa) == 16 || fabs(bbb) == 16 || fabs(ccc) == 16 ) {
+ if( abs(aaa) == 16 || abs(bbb) == 16 || abs(ccc) == 16 ) {
_cL = _cZPnutnut_L; _cR = _cZPnutnut_R;
}
right(_cR);
left(_cL);
norm(_cZP_o);
}

File Metadata

Mime Type
text/x-diff
Expires
Tue, Nov 19, 3:25 PM (1 d, 16 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3804944
Default Alt Text
(32 KB)

Event Timeline