Page MenuHomeHEPForge
Files F8308608

example.cc
No OneTemporary
Actions

example.cc
View Options

//----------------------------------------------------------------------
// Example on how to use this contribution
//
// run it with
// ./example < ../data/single-event.dat
//----------------------------------------------------------------------
//
// Copyright (c) 2014-, Marcel Vos and Ignacio Garcia
//
//----------------------------------------------------------------------
// This file is part of FastJet contrib.
//
// It is free software; you can redistribute it and/or modify it under
// the terms of the GNU General Public License as published by the
// Free Software Foundation; either version 2 of the License, or (at
// your option) any later version.
//
// It is distributed in the hope that it will be useful, but WITHOUT
// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
// or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public
// License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this code. If not, see <http://www.gnu.org/licenses/>.
//----------------------------------------------------------------------
#include <iostream>
#include <iomanip>
#include "fastjet/JetDefinition.hh"
#include "fastjet/PseudoJet.hh"
#include "fastjet/Selector.hh"
#include "ValenciaPlugin.hh"
#include "fastjet/EECambridgePlugin.hh"
using namespace std;
using namespace fastjet;
void read_event(vector<PseudoJet> &full_event);
//----------------------------------------------------------------------
int main(){
// example for the creation and execution of the fastjet plugin
// for the Valencia jet algorithm
//
// the algorithm is intended for e+e- collisions, but can also be run
// on the simulated LHC data that ship with the fjcontrib area
// ./example ../data/Pythia-Zp2jets-lhc-pileup-1ev.dat
//
// read in input particles
//----------------------------------------------------------
vector<PseudoJet> full_event;
read_event(full_event);
full_event = SelectorAbsRapMax(4.0)(full_event);
// create the jet definition using the plugin mechanism
//----------------------------------------------------------
fastjet::contrib::ValenciaPlugin * vlc_plugin = new fastjet::contrib::ValenciaPlugin(1.2, 0.8);
fastjet::JetDefinition jet_def(vlc_plugin);
// do the clustering and get the jets
// only exclusive clustering is supported
//-----------------------------------------------------------
ClusterSequence clust_seq(full_event, jet_def);
vector<PseudoJet> jets = clust_seq.inclusive_jets(20.);
cout << "hard jets (pT > 20 GeV) in inclusive clustering " << endl;
for (unsigned int i=0; i<jets.size(); i++){
const PseudoJet &jet = jets[i];
cout << "pt = " << jet.pt()
<< ", rap = " << jet.rap() << endl;
}
cout << endl;
jets = clust_seq.exclusive_jets(4);
// print the result
//----------------------------------------------------------
cout << "hard jets in exclusive N=4 clustering " << endl;
for (unsigned int i=0; i<jets.size(); i++){
const PseudoJet &jet = jets[i];
cout << "pt = " << jet.pt()
<< ", rap = " << jet.rap() << endl;
}
cout << endl;
jets = clust_seq.exclusive_jets(500.);
// print the result
//----------------------------------------------------------
cout << "hard jets in exclusive clustering up to d = 500" << endl;
for (unsigned int i=0; i<jets.size(); i++){
const PseudoJet &jet = jets[i];
cout << "pt = " << jet.pt()
<< ", rap = " << jet.rap() << endl;
}
cout << endl;
return 0;
}
//------------------------------------------------------------------------
// read the event
void read_event(vector<PseudoJet> &full_event){
string line;
while (getline(cin, line)) {
istringstream linestream(line);
// take substrings to avoid problems when there are extra "pollution"
// characters (e.g. line-feed).
double px,py,pz,E;
linestream >> px >> py >> pz >> E;
PseudoJet particle(px,py,pz,E);
// push event onto back of full_event vector
full_event.push_back(particle);
}
}

File Metadata

Mime Type
text/x-c
Expires
Sat, Dec 21, 12:35 PM (1 d, 20 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4022792
Default Alt Text
example.cc (3 KB)

Event Timeline