Page MenuHomeHEPForge
Files F9501698

No OneTemporary
Actions

View Options

diff --git a/Models/General/HPDiagram.h b/Models/General/HPDiagram.h
--- a/Models/General/HPDiagram.h
+++ b/Models/General/HPDiagram.h
@@ -1,216 +1,219 @@
// -*- C++ -*-
//
// HPDiagram.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2017 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_HPDiagram_H
#define HERWIG_HPDiagram_H
//
// This is the declaration of the HPDiagram struct.
//
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#include "ThePEG/Helicity/Vertex/VertexBase.h"
namespace Herwig {
using namespace ThePEG;
using Helicity::VertexBasePtr;
/** Pair of particle ids. */
typedef pair<long, long> IDPair;
/** Pair of bool's*/
typedef pair<bool, bool> BPair;
/** Convenient typedef of VertexBasePtr */
typedef VertexBasePtr VBPtr;
/** Pair of VertexBasePtrs */
typedef pair<VBPtr, VBPtr> VBPair;
/** Pair of int,double */
typedef pair<unsigned int, double> CFPair;
/**
* The HPDiagram struct contains information about a 2->2 hard-process that
* has been automatically generated by HardProcessConstructor.
*/
struct HPDiagram {
/** Enumeration for channel type */
enum Channel {UNDEFINED = -1, sChannel, tChannel, fourPoint};
/** Standard Constructor */
HPDiagram() : incoming(make_pair(0, 0)), outgoing(make_pair(0, 0)),
ordered(make_pair(true,true)), channelType(UNDEFINED),
colourFlow(0), ids(4, 0) {}
/** Constructor taking ids as arguments.*/
HPDiagram(IDPair a, IDPair b)
: incoming(a), outgoing(b), ordered(make_pair(true,true)),
channelType(UNDEFINED), colourFlow(0), ids(4, 0){
ids[0] = incoming.first;
ids[1] = incoming.second;
ids[2] = outgoing.first;
ids[3] = outgoing.second;
}
/** Incoming particle id's */
IDPair incoming;
/** Outgoing particle id's */
IDPair outgoing;
/** Particle ordering for t-channel diagrams */
BPair ordered;
/** ParticleData pointer to intermediate, null for 4-point vertices */
PDPtr intermediate;
/** The two vertices for the diagram */
VBPair vertices;
/** Enum of channel type */
Channel channelType;
/** Store colour flow information */
vector<CFPair> colourFlow;
/** Store the ids in a vector for easy use of comparison operator. */
vector<long> ids;
/**
* Test whether this and x are the same process
* @param x The other process to check
*/
bool sameProcess(const HPDiagram & x) const {
return ( x.incoming == incoming && x.outgoing == outgoing );
}
};
/**
* Test whether two diagrams are identical.
*/
inline bool operator==(const HPDiagram & x, const HPDiagram & y) {
// check incoming
if( x.incoming != y.incoming) return false;
// check type of diagram
if( x.channelType != y.channelType ) return false;
if( x.outgoing == y.outgoing &&
x.ordered == y.ordered ) {
// 4 point
if(x.channelType==HPDiagram::fourPoint) {
// check vertex
return x.vertices.first==y.vertices.first;
}
// all others
else {
// check vertex
if(x.vertices.first != y.vertices.first ||
x.vertices.second != y.vertices.second )
return false;
// check intermediate
return x.intermediate == y.intermediate;
}
}
//diagram is also the same if the outgoing particles are
//swapped and the ordering is opposite
else if( x.outgoing.first != x.outgoing.second &&
y.outgoing.first != y.outgoing.second &&
x.outgoing.first == y.outgoing.second &&
x.outgoing.second == y.outgoing.first &&
x.channelType == y.channelType) {
// 4 point
if(x.channelType==HPDiagram::fourPoint) {
// check vertex
return x.vertices.first==y.vertices.first;
}
// t/u channel
else if(x.channelType==HPDiagram::tChannel) {
- // check vertex
- if(x.vertices.first != y.vertices.first ||
- x.vertices.second != y.vertices.second)
- return false;
+ // check vertex and intermediate
+ if(x.vertices.first == y.vertices.first &&
+ x.vertices.second == y.vertices.second &&
+ x.intermediate == y.intermediate )
+ return true;
tPDPtr inter = x.intermediate;
if(inter->CC()) inter = inter->CC();
- // check intermediate
- return inter == y.intermediate;
+ if(x.vertices.first == y.vertices.second &&
+ x.vertices.second == y.vertices.first &&
+ inter == y.intermediate )
+ return true;
}
// s-channel
else {
// check vertex
if(x.vertices.first != y.vertices.first ||
x.vertices.second != y.vertices.second )
return false;
// check intermediate
return x.intermediate == y.intermediate;
}
}
return false;
}
/**
* Test whether one diagram is 'less' than another. Does a
* lexicographic comparison of the external states.
*/
inline bool operator<(const HPDiagram & x, const HPDiagram & y) {
for( size_t i = 0; i < 4; ++i ) {
if( x.ids[i] < y.ids[i] )
return true;
else if( x.ids[i] > y.ids[i] )
return false;
}
return false;
}
/**
* Output to a stream
*/
inline ostream & operator<<(ostream & os, const HPDiagram & diag) {
os << diag.incoming.first << " " << diag.incoming.second << " -> ";
if(diag.intermediate)
os << diag.intermediate->id() << " -> ";
os << diag.outgoing.first << " " << diag.outgoing.second
<< " channel: " << diag.channelType << " ";
if(diag.channelType == HPDiagram::tChannel)
os << "ordering " << diag.ordered.first << " "
<< diag.ordered.second << " ";
for(size_t cf = 0; cf < diag.colourFlow.size(); ++cf)
os << "(" << diag.colourFlow[cf].first << ","
<<diag.colourFlow[cf].second << ")";
os << '\n';
return os;
}
/**
* Output operator to allow the structure to be persistently written
* @param os The output stream
* @param x The HPDiagram
*/
inline PersistentOStream & operator<<(PersistentOStream & os,
const HPDiagram & x) {
os << x.incoming << x.outgoing << x.ordered << x.intermediate
<< x.vertices << x.channelType << x.colourFlow << x.ids;
return os;
}
/**
* Input operator to allow persistently written data to be read in
* @param is The input stream
* @param x The HPDiagram
*/
inline PersistentIStream & operator>>(PersistentIStream & is,
HPDiagram & x) {
int chan = -1;
is >> x.incoming >> x.outgoing >> x.ordered >> x.intermediate
>> x.vertices >> chan >> x.colourFlow >> x.ids;
x.channelType = HPDiagram::Channel(chan);
return is;
}
}
#endif /* HERWIG_HPDiagram_H */

File Metadata

Mime Type
text/x-diff
Expires
Sun, Feb 23, 2:59 PM (5 h, 14 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4486766
Default Alt Text
(6 KB)

Event Timeline