FastJet: example/88-HI-pt-asymmetry.cc Source File

00001 #include "fastjet/ClusterSequeneArea.hh"
00002 #include "fastjet/tools/Selector.hh"
00003 #include "fastjet/tools/BackgroundEstimator"
00004 #include<iostream> // needed for io
00005 
00006 using namespace std;
00007 using namespace fastjet;
00008 
00009 
00010 int main() {
00011 
00012   Selector jet_acceptance = SelectorPtMin(25.0) && SelectorAbsRapMax(2.8);
00013   Selector background_acceptance = SelectorStripRange(1.0) * (!SelectorNHardest(2));
00014   JetDefinition jet_def_akt04 = JetDefinition(antikt_algorithm, 0.4);
00015   JetDefinition jet_def_bkgd  = JetDefinition(kt_algorithm, 0.4);
00016   AreaDefinition area_def; // can the default work here?
00017 
00018   // read in an event
00019   vector<PseudoJet> input_particles;
00020   double px, py , pz, E;
00021   while (cin >> px >> py >> pz >> E) {
00022     input_particles.push_back(fastjet::PseudoJet(px,py,pz,E)); 
00023   }
00024 
00025   if (jets.size() > 2 && jets[0].perp() > 100) {
00026     double A_J = (jets[0].perp()-jets[1].perp())/(jets[0].perp()+jets[1].perp());
00027   }
00028 }