Page MenuHomeHEPForge
Files F8725607

No OneTemporary
Actions

View Options

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) {
_cR = _cAGQQ_R;
_cL = _cAGQQ_L;
} else {
_cR = _cAGTT_R;
_cL = _cAGTT_L;
}
}
- if(_models!=2) { _cL = 0; _cR = 0; }
+ 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) {
_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) {
_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 = 0; }
+ 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) {
_cR = _cWPTD_R;
_cL = _cWPTD_L;
}
- if(_models!=0) { _cL = 0; _cR = 0; }
+ 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) {
_cR = _cZPTU_R;
_cL = _cZPTU_L;
} else {
if( fabs(aa->id()) != 4 && fabs(bb->id()) != 4 && fabs(cc->id()) != 4) {
_cR = _cZPUU_R;
_cL = _cZPUU_L;
}
if( fabs(aa->id()) == 4 || fabs(bb->id()) == 4 || fabs(cc->id()) == 4) {
_cR = _cZPCC_R;
_cL = _cZPCC_L;
}
}
- if(_models!=1) { _cL = 0; _cR = 0; }
+ if(_models!=1) { _cL = 1E-10; _cR = 1E-10; }
right(_cR);
left(_cL);
norm(1.0);
}

File Metadata

Mime Type
text/x-diff
Expires
Tue, Jan 21, 2:06 AM (1 d, 14 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4237515
Default Alt Text
(13 KB)

Event Timeline