1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
|
// -*- C++ -*-
//
// MultiCutBase.cc is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2011 Leif Lonnblad
//
// ThePEG 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 MultiCutBase class.
//
#include "MultiCutBase.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/EventRecord/Particle.h"
#include "ThePEG/Repository/CurrentGenerator.h"
using namespace ThePEG;
void MultiCutBase::describe() const {
CurrentGenerator::log() << fullName() << " has no description.\n\n";
}
Energy2 MultiCutBase::minS(const tcPDVector &) const {
return ZERO;
}
Energy2 MultiCutBase::maxS(const tcPDVector &) const {
return Constants::MaxEnergy2;
}
bool MultiCutBase::passCuts(tcCutsPtr, const tcPDVector & ptype,
const vector<LorentzMomentum> & p) const {
long NN = (1 << ptype.size());
// NN is the number of different combinations that can be made.
for ( long ii = 1; ii < NN; ++ii ) {
long mask = ii;
tcPDVector pt;
LorentzMomentum sum;
int i = -1;
while ( mask ) {
++i;
int sel = mask%2;
mask /= 2;
if ( !sel ) continue;
pt.push_back(ptype[i]);
sum += p[i];
}
if ( pt.size() < 2 ) continue;
if ( sum.m2() < minS(pt) ) return false;
if ( sum.m2() >= maxS(pt) ) return false;
}
return true;
}
bool MultiCutBase::passCuts(tcCutsPtr parent, const tcPVector & p) const {
tcPDVector ptype(p.size());
vector<LorentzMomentum> mom(p.size());
for ( int i = 0, N = p.size(); i < N; ++i ) {
ptype[i] = p[i]->dataPtr();
mom[i] = p[i]->momentum();
}
return passCuts(parent, ptype, mom);
}
AbstractNoPIOClassDescription<MultiCutBase> MultiCutBase::initMultiCutBase;
// Definition of the static class description member.
void MultiCutBase::Init() {
static ClassDocumentation<MultiCutBase> documentation
("This class corresponds to a kinematical cut to be made on a set of "
"outgoing particles from a hard sub-process.");
}
|
