Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F7877359
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Flag For Later
Size
32 KB
Subscribers
None
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
Details
Attached
Mime Type
text/x-diff
Expires
Tue, Nov 19, 3:25 PM (1 d, 19 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3804944
Default Alt Text
(32 KB)
Attached To
rHERWIGHG herwighg
Event Timeline
Log In to Comment