Page MenuHomeHEPForge
Files F8723982

No OneTemporary
Actions

View Options

Index: contrib/contribs/ValenciaPlugin/trunk/NEWS
===================================================================
--- contrib/contribs/ValenciaPlugin/trunk/NEWS (revision 1020)
+++ contrib/contribs/ValenciaPlugin/trunk/NEWS (revision 1021)
@@ -1,5 +1,6 @@
+15 June 2017: bug-fix version 2.0.1
24 February 2015: version 2.0.0
21 February 2015: released 2.0-devel
27 January 2015: creating version 2.0-devel
13 March 2014: first development version 1.0-devel
Index: contrib/contribs/ValenciaPlugin/trunk/AUTHORS
===================================================================
--- contrib/contribs/ValenciaPlugin/trunk/AUTHORS (revision 1020)
+++ contrib/contribs/ValenciaPlugin/trunk/AUTHORS (revision 1021)
@@ -1,7 +1,9 @@
The ValenciaJetAlgorithm FastJet contrib was written and is maintained and developed by: Ignacio Garcia Garcia, Marcel Vos, IFIC (UVEG/CSIC) Valencia, Spain
Acknowledging suggestions by Gavin P. Salam (CERN) and J. Thaler (MIT)
-A preprint motivating the algorithm and describing its performance is available on: "A new jet algorithm for lepton colliders", arXiv:1404.4294 (submitted to PLB)
+A paper motivating the algorithm and describing its performance is available on: "A robust jet algorithm for lepton colliders", arXiv:1404.4294, PLB750, (2015) 95-99
+
+More detailed studies are presented in "Jet reconstruction at high-energy lepton colliders", arXiv:1607.05039 (submitted to EPJC).
Please contact: marcel.vos@ific.uv.es for questions, queries and comments
Index: contrib/contribs/ValenciaPlugin/trunk/VERSION
===================================================================
--- contrib/contribs/ValenciaPlugin/trunk/VERSION (revision 1020)
+++ contrib/contribs/ValenciaPlugin/trunk/VERSION (revision 1021)
@@ -1 +1 @@
-2.0.0
+2.0.1
Index: contrib/contribs/ValenciaPlugin/trunk/ValenciaPlugin.cc
===================================================================
--- contrib/contribs/ValenciaPlugin/trunk/ValenciaPlugin.cc (revision 1020)
+++ contrib/contribs/ValenciaPlugin/trunk/ValenciaPlugin.cc (revision 1021)
@@ -1,133 +1,133 @@
// $Id$
//
// 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 "ValenciaPlugin.hh"
#include "fastjet/NNH.hh"
// strings and streams
#include <sstream>
#include <limits>
FASTJET_BEGIN_NAMESPACE // defined in fastjet/internal/base.hh
namespace contrib{
//----------------------------------------------------------------------
/// class that contains the algorithm parameters R and beta
// needed to run the Valencia algorithm
class ValenciaInfo {
public:
ValenciaInfo(double Ri, double betai, double gammai)
{ R_ = Ri; beta_ = betai; gamma_ = gammai;}
double beta() { return beta_; }
double gamma() { return gamma_; }
double R() { return R_; }
private:
double R_, beta_, gamma_;
};
class ValenciaBriefJet {
public:
void init(const PseudoJet & jet, ValenciaInfo * info) {
double norm = 1.0/sqrt(jet.modp2());
// pseudo-jet information needed to calculate distance
nx = jet.px() * norm;
ny = jet.py() * norm;
nz = jet.pz() * norm;
E = jet.E();
// algorithm parameters needed in distance calculation
R = info->R();
beta = info->beta();
// beam distance
// E-to-the-2*beta times sin(polar angle)-to-the-2*gamma
if (E==0. || jet.perp()==0.) diB=0.;
diB = pow(E,2*beta) * pow(jet.perp()/E,2*info->gamma());
}
double distance(const ValenciaBriefJet * jet) const {
// inter-particle distance is similar to Durham and e+e- kt
// two differences:
// R is redefined for greater flexibility
// the energy is elevated to 2*beta
double dij = 1 - nx*jet->nx
- ny*jet->ny
- nz*jet->nz;
- if (jet->E < E)
+ if (pow(jet->E,2*beta) < pow(E,2*beta))
dij *= 2 * pow(jet->E,2*beta);
else
dij *= 2 * pow(E,2*beta);
dij/=pow(R,2);
return dij;
}
double beam_distance() const {
return diB;
}
double E, nx, ny, nz;
double diB;
double R, beta;
};
std::string ValenciaPlugin::description () const {
std::ostringstream desc;
desc << "Valencia plugin with R = " << R() << ", beta = " << beta() << " and gamma = " << gamma();
return desc.str();
}
void ValenciaPlugin::run_clustering(fastjet::ClusterSequence & cs) const {
int njets = cs.jets().size();
ValenciaInfo vinfo(R(),beta(),gamma());
NNH<ValenciaBriefJet,ValenciaInfo> nnh(cs.jets(),&vinfo);
while (njets > 0) {
int i, j, k;
double dij = nnh.dij_min(i, j);
if (j >= 0) {
cs.plugin_record_ij_recombination(i, j, dij, k);
nnh.merge_jets(i, j, cs.jets()[k], k);
} else {
cs.plugin_record_iB_recombination(i, dij);
nnh.remove_jet(i);
}
njets--;
}
}
} // namespace contrib
FASTJET_END_NAMESPACE

File Metadata

Mime Type
text/x-diff
Expires
Mon, Jan 20, 10:10 PM (1 d, 10 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4242591
Default Alt Text
(5 KB)

Event Timeline