Example of analysis code:
public void process(EventHeader event)
{
// Get reconstructed particle 3Vector's.
List<Hep3Vector> particles = new ArrayList();
List<List<ReconstructedParticle>> rpLists = event.get(ReconstructedParticle.class);
for (List<ReconstructedParticle> collection : rpLists) {
String name = event.getMetaData(collection).getName();
if (!name.equals("ReconCheater")) continue;
for (ReconstructedParticle rp : collection) {
particles.add(rp.getMomentum());
}
}
// Set up Jet finder.
double ycut = 0.0005, dycut = 0.0005;
finder.setYCut(ycut);
finder.setEvent(particles);
NJets = finder.njets();
while (NMJets>MaxNumberJetsExpected) {
ycut += dycut;
finder.setYCut(ycut);
NMJets = finder.njets();
}
aida.cloud1D("# jets").fill(NJets);
}
|