80 |
|
// the main routine: |
81 |
|
// jet finding using FastJet |
82 |
|
|
83 |
– |
cout << "Main finding routine..." << endl; |
84 |
– |
|
83 |
|
// clear the list if it is filled |
84 |
|
if (JetsOut.size()!=0){ |
85 |
|
JetsOut.clear(); |
126 |
|
for (unsigned int ijet = 0; ijet < SortedJets.size(); ijet++) { |
127 |
|
Jet* jet = new Jet(); |
128 |
|
JetsOut.push_back(jet); |
129 |
< |
jet->pt = SortedJets[ijet].pt(); |
130 |
< |
jet->eta = SortedJets[ijet].eta(); |
131 |
< |
jet->phi = SortedJets[ijet].phi(); |
132 |
< |
jet->energy = SortedJets[ijet].E(); |
129 |
> |
jet->set_pt ( SortedJets[ijet].pt()); |
130 |
> |
jet->set_eta ( SortedJets[ijet].eta()); |
131 |
> |
jet->set_phi ( SortedJets[ijet].phi()); |
132 |
> |
jet->set_energy ( SortedJets[ijet].E()); |
133 |
|
} |
134 |
|
|
135 |
|
// get the association of jetparticles to jets |