14 using namespace JetReco;
23 list<InputItem>
input;
24 for (InputCollection::const_iterator towerIter = fInput.begin();
25 towerIter != fInput.end(); ++towerIter) {
26 input.push_back(*towerIter);
29 input.sort(compCandidate);
32 while( !input.empty() && input.front()->et() > theSeedThreshold ) {
34 double eta0=input.front()->eta();
35 double phi0=input.front()->phi();
41 list< list<InputItem>::iterator> cone;
47 for(list<InputItem>::iterator inp=input.begin();
48 inp!=input.end();inp++){
51 theConeRadius*theConeRadius) {
53 eta+= tower->
et()*tower->
eta();
54 double dphi=tower->
phi()-phi0;
57 phi+=tower->
et()*dphi;
66 if(fabs(eta-eta0)<.001 && fabs(phi-phi0)<.001)
break;
73 list< list<InputItem>::iterator>::const_iterator inp;
74 for(inp=cone.begin();inp!=cone.end();inp++) {
75 jetConstituents.push_back(**inp);
78 fOutput->push_back (
ProtoJet (jetConstituents));
82 sort (fOutput->begin (), fOutput->end (), compJets);
std::vector< ProtoJet > OutputCollection
virtual double et() const =0
transverse energy
std::vector< InputItem > InputCollection
Transient Jet class used by the reconstruction algorithms.
double deltaR2(const Vector1 &v1, const Vector2 &v2)
void run(const JetReco::InputCollection &fInput, JetReco::OutputCollection *fOutput) const
Find the ProtoJets from the collection of input Candidates.
virtual double phi() const =0
momentum azimuthal angle
virtual double eta() const =0
momentum pseudorapidity