Page MenuHomeHEPForge
Files F10881744

No OneTemporary
Actions

View Options

diff --git a/MatrixElement/Hadron/MEMinBias.cc b/MatrixElement/Hadron/MEMinBias.cc
--- a/MatrixElement/Hadron/MEMinBias.cc
+++ b/MatrixElement/Hadron/MEMinBias.cc
@@ -1,158 +1,155 @@
// -*- C++ -*-
//
// This is the implementation of the non-inlined, non-templated member
// functions of the MEMinBias class.
//
#include "MEMinBias.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Utilities/SimplePhaseSpace.h"
//#include "ThePEG/Repository/EventGenerator.h"
#include "ThePEG/Handlers/StandardXComb.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
using namespace Herwig;
#include "ThePEG/PDT/EnumParticles.h"
#include "ThePEG/MatrixElement/Tree2toNDiagram.h"
void MEMinBias::getDiagrams() const {
- int maxflav(5);
+ int maxflav(2);
tcPDPtr ph = getParticleData(ParticleID::gamma);
for ( int i = 1; i <= maxflav; ++i ) {
for( int j=1; j <= i; ++j){
tcPDPtr q1 = getParticleData(i);
tcPDPtr q1b = q1->CC();
tcPDPtr q2 = getParticleData(j);
tcPDPtr q2b = q2->CC();
// For each flavour we add:
//qq -> qq
add(new_ptr((Tree2toNDiagram(3), q1, ph, q2, 1, q1, 2, q2, -1)));
//qqb -> qqb
add(new_ptr((Tree2toNDiagram(3), q1, ph, q2b, 1, q1, 2, q2b, -2)));
//qbq -> qbq
- add(new_ptr((Tree2toNDiagram(3), q1b, ph, q2, 1, q1b, 2, q2, -3)));
+ // add(new_ptr((Tree2toNDiagram(3), q1b, ph, q2, 1, q1b, 2, q2, -3)));
//qbqb -> qbqb
- add(new_ptr((Tree2toNDiagram(3), q1b, ph, q2b, 1, q1b, 2, q2b, -4)));
+ add(new_ptr((Tree2toNDiagram(3), q1b, ph, q2b, 1, q1b, 2, q2b, -3)));
}
}
}
Energy2 MEMinBias::scale() const {
- return sqr(2*GeV);
+ //set the scale high as long as I have no pure valence pdf
+ return sqr(100*GeV);
}
int MEMinBias::nDim() const {
return 1;
}
void MEMinBias::setKinematics() {
MEBase::setKinematics(); // Always call the base class method first.
}
bool MEMinBias::generateKinematics(const double *) {
// generate the masses of the particles
for ( int i = 2, N = meMomenta().size(); i < N; ++i ) {
meMomenta()[i] = Lorentz5Momentum(mePartonData()[i]->generateMass());
}
Energy q = 0.0*GeV;
try {
q = SimplePhaseSpace::
getMagnitude(sHat(), meMomenta()[2].mass(), meMomenta()[3].mass());
} catch ( ImpossibleKinematics ) {
return false;
}
Energy pt = 0*GeV;
meMomenta()[2].setVect(Momentum3( pt, pt, q));
meMomenta()[3].setVect(Momentum3(-pt, -pt, -q));
meMomenta()[2].rescaleEnergy();
meMomenta()[3].rescaleEnergy();
jacobian(1.0);
return true;
}
double MEMinBias::me2() const {
return 1.0;
}
CrossSection MEMinBias::dSigHatDR() const {
return me2()*jacobian()/sHat()*sqr(hbarc);
}
unsigned int MEMinBias::orderInAlphaS() const {
return 0;
}
unsigned int MEMinBias::orderInAlphaEW() const {
return 2;
}
Selector<MEBase::DiagramIndex>
MEMinBias::diagrams(const DiagramVector & diags) const {
- // This example corresponds to the diagrams specified in the example
- // in the getDiagrams() function.
-
Selector<DiagramIndex> sel;
for ( DiagramIndex i = 0; i < diags.size(); ++i )
sel.insert(1.0, i);
return sel;
}
Selector<const ColourLines *>
MEMinBias::colourGeometries(tcDiagPtr diag) const {
static ColourLines qq("1 4, 3 5");
static ColourLines qqb("1 4, -3 -5");
- static ColourLines qbq("-1 -4, 3 5");
+ // static ColourLines qbq("-1 -4, 3 5");
static ColourLines qbqb("-1 -4, -3 -5");
Selector<const ColourLines *> sel;
switch(diag->id()){
case -1:
sel.insert(1.0, &qq);
break;
case -2:
sel.insert(1.0, &qqb);
break;
+ // case -3:
+ // sel.insert(1.0, &qbq);
+ //break;
case -3:
- sel.insert(1.0, &qbq);
- break;
- case -4:
sel.insert(1.0, &qbqb);
break;
}
-
return sel;
}
IBPtr MEMinBias::clone() const {
return new_ptr(*this);
}
IBPtr MEMinBias::fullclone() const {
return new_ptr(*this);
}
ClassDescription<MEMinBias> MEMinBias::initMEMinBias;
// Definition of the static class description member.
void MEMinBias::Init() {
static ClassDocumentation<MEMinBias> documentation
("There is no documentation for the MEMinBias class");
}

File Metadata

Mime Type
text/x-diff
Expires
Sat, May 3, 6:47 AM (1 d, 8 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4983137
Default Alt Text
(4 KB)

Event Timeline