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
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
|
#ifndef JETREC_JETSPLITMERGETOOL_H
#define JETREC_JETSPLITMERGETOOL_H
//----------------------------------------------------------------------
// This file distributed with FastJet has been obtained from SpartyJet
// v2.20.0 by Pierre-Antoine Delsart, Kurtis L. Geerlings, Joey
// Huston, Brian T. Martin and Chris Vermilion
// For details, see http://www.pa.msu.edu/~huston/SpartyJet/
// http://projects.hepforge.org/spartyjet/
//
// Changes from the original file are listed below.
//----------------------------------------------------------------------
//******************************************************************************
// Filename : JetSplitMergeTool
// Author : Ambreesh Gupta
// Created : November, 2001
//
// DESCRIPTION:
//
// Cone algorithms generally produce overlapping jets. These jets are then
// passed through a Split Merge algorithm to create seperated jets. This
// is an implementation of one such algorithm.
//
// The steps of the algorithm are:
// - <to be written>
//
// PROPERTIES (JobOption Parameters):
//
// ProtoJetContainerLoc string "Default" Key for ProtoJet input.
// JetContainerLoc string "Default" Key for Jet list to output.
// OverlapFraction double 0.5
//
// HISTORY
// 05Nov01 agupta. First version
// 15Aug02 agupta. Fix phi-wrapping.
//*****************************************************************************
// History of changes from the original JetSplitMergeTool.hh file in
// SpartyJet v2.20
//
// 2009-01-15 Gregory Soyez <soyez@fastjet.fr>
//
// * put the code in the fastjet::atlas namespace
//
// 2009-02-14 Gregory Soyez <soyez@fastjet.fr>
//
// * imported into FastJet
// * removed the message logs
// * replaced StatusCode by int
//class StoreGateSvc;
//Includes
#include "Jet.hh"
#include <vector>
#include <list>
//#include "JetCore/CustomMessage.hh"
#include "fastjet/internal/base.hh"
FASTJET_BEGIN_NAMESPACE
namespace atlas {
//typedef int StatusCode ;
class JetSplitMergeTool
{
public:
typedef Jet::jet_list_t jetcollection_t;
//Constructors
JetSplitMergeTool();
//Destructor
virtual ~JetSplitMergeTool();
virtual int execute( jetcollection_t* theJets );
//Handle negative energy ProtoJets:
double etaTrue(Jet::constit_vect_t::iterator pj);
double phiTrue(Jet::constit_vect_t::iterator pj);
//private:
void split_merge();
double m_f;
jetcollection_t m_preJet;
jetcollection_t m_jet;
jetcollection_t* m_jetVec;
int m_ctr;
int m_dctr;
// Message m_log;
};
} // namespace atlas
FASTJET_END_NAMESPACE
#endif // JETREC_JETSPLITMERGETOOL_H
|
