/
Py8toJetInput.h
51 lines (34 loc) · 1.4 KB
/
Py8toJetInput.h
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
#ifndef gen_Py8toJetInput_h
#define gen_Py8toJetInput_h
#include "Pythia8/FastJet3.h" // Py8 overhead on top of FastJets package
#include "Pythia8/Event.h"
namespace lhef {
class LHEEvent;
}
class Py8toJetInput
{
public:
typedef Pythia8::Event Event;
typedef Pythia8::Particle Particle;
Py8toJetInput(): fJetEtaMax(10.) {}
~Py8toJetInput() {}
virtual const std::vector<fastjet::PseudoJet> fillJetAlgoInput( const Event&, const Event&,
const lhef::LHEEvent* lhee=0,
const std::vector<int>* partonList=0 );
void setJetEtaMax( double max ) { fJetEtaMax=max; return; }
protected:
enum partonTypes { ID_TOP=6, ID_GLUON=21, ID_PHOTON=22 };
double fJetEtaMax;
int getAncestor( int, const Event&, const Event& );
std::vector<fastjet::PseudoJet> fJetInput;
};
class Py8toJetInputHEPEVT : public Py8toJetInput
{
public:
Py8toJetInputHEPEVT() {}
~Py8toJetInputHEPEVT() {}
const std::vector<fastjet::PseudoJet> fillJetAlgoInput( const Event&, const Event&,
const lhef::LHEEvent*,
const std::vector<int>* partonList=0 );
};
#endif