Page MenuHomeHEPForge
Files F10881894

example.cc
No OneTemporary
Actions

example.cc
View Options

#include <cstddef>
#include <fstream>
#include <iostream>
#include <iomanip>
#include <vector>
#include <fastjet/ClusterSequence.hh>
#include "KTClusCXX.hh"
using namespace fastjet;
class Input {
public:
Input() {};
~Input() {};
const double PIMASS = 0.139;
std::vector <PseudoJet> input_particles;
std::size_t N;
void fill(const std::string& inp="")
{
input_particles.clear();
std::size_t i = 0;
double px, py, pz, E;
if (inp.size()) {
std::ifstream myfile;
myfile.open (inp.c_str());
while (myfile >> px >> py >> pz >> E) {
input_particles.push_back(PseudoJet(px, py, pz, E));
i++;
}
myfile.close();
} else {
while (std::cin >> px >> py >> pz >> E) {
input_particles.push_back(PseudoJet(px, py, pz, E));
i++;
}
}
N = i;
}
};
inline std::ostream& operator<<(std::ostream& os, const PseudoJet &j) {
os << j.px() << " " << j.py() << " " << j.pz() << " " << j.e();
return os;
}
int main( int argc, char** argv)
{
///2111 1.0 0.0001 1.0
if (argc<2) std::cout <<"#Usage: "<<argv[0]<<" <input file=standard input> <mode=2111> <Radius=1.0> <YCUT=0.0001> <ECUT=1.0>" <<std::endl;
std::string smode = argc>2?std::string(argv[2]):"2111";
bool inclusive = (smode.length() > 4 && smode[4] == 'i');
int mode = std::atoi(smode.substr(0, 4).c_str());
double Radius = argc>3?std::atof(argv[3]):1.0;
double YCUT = argc>4?std::atof(argv[4]):0.0001;
double ECUT =argc>5?std::atof(argv[5]):1.0;
Input IN;
IN.fill(argc>1?argv[1]:"");
std::vector < PseudoJet > fastjets;
fastjets.reserve(IN.N);
ClusterSequence cs(IN.input_particles, contrib::KTClusCXXJetDefinition(mode, Radius));
if (inclusive) {
fastjets = cs.inclusive_jets();
} else {
fastjets = cs.exclusive_jets(ECUT*ECUT*YCUT);
}
std::vector < double > fastjets_dmerge;
std::sort(fastjets.begin(), fastjets.end(), [](const PseudoJet & a, const PseudoJet & b) -> bool { return a.E() > b.E(); });
for (std::size_t i = 0; i < IN.N; i++) {
fastjets_dmerge.push_back(cs.exclusive_dmerge_max(i));
}
std::cout << std::fixed << std::setw(11) << std::setprecision(6);
std::cout<<"Input particles:"<<std::endl;
for (auto& i: IN.input_particles) {
std::cout<<i<<std::endl;
}
std::cout<<"Parameters: "<< smode << " " << Radius << " " << YCUT<< " " << ECUT <<std::endl;
std::cout<<"Pseudojets:"<<std::endl;
for (auto& j: fastjets) {
std::cout<<j<<std::endl;
}
std::cout<<"exclusive_dmerge:"<<std::endl;
for (auto& d: fastjets_dmerge) {
std::cout<<std::min(1000000.0,d)<<std::endl;
}
return 0;
}

File Metadata

Mime Type
text/x-c++
Expires
Sat, May 3, 6:55 AM (18 h, 52 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4983194
Default Alt Text
example.cc (2 KB)

Event Timeline