Index: contrib/contribs/Centauro/tags/1.0.0/Centauro.cc
===================================================================
--- contrib/contribs/Centauro/tags/1.0.0/Centauro.cc (revision 0)
+++ contrib/contribs/Centauro/tags/1.0.0/Centauro.cc (revision 1278)
@@ -0,0 +1,134 @@
+// This code 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 .
+//----------------------------------------------------------------------x
+
+#include "Centauro.hh"
+#include "fastjet/NNH.hh"
+
+// strings and streams
+#include
+#include
+
+FASTJET_BEGIN_NAMESPACE // defined in fastjet/internal/base.hh
+
+namespace contrib{
+
+ //----------------------------------------------------------------------
+ /// class that contains the algorithm parameters R, photon energy and photon longitudinal momentum.
+ class CentauroInfo {
+
+ public:
+ CentauroInfo(double Ri, double gammaEi, double gammaPzi)
+ { R_ = Ri; gammaE_ = gammaEi; gammaPz_ = gammaPzi;}
+
+ double gammaPz() { return gammaPz_; }
+ double gammaE() { return gammaE_; }
+ double R() { return R_; }
+
+ private:
+ double R_, gammaE_, gammaPz_;
+ };
+
+ class CentauroBriefJet {
+ public: //For definitions see https://arxiv.org/abs/2006.10751
+ // n = (1,0,0,1)
+ // nbar = (1,0,0,-1)
+ // P = Q/2x(1,0,0,1) //proton
+ // q = Q(0,0,0,-1) //virtual photon
+ // etabar = -2Q/(nbar*q)*pT/(n*p) , so in the Breit frame: // etabar = +2*pT/(n*p) = 2*pT/(E-pz)
+ // The distance is: (for f=x)
+ // dij = [(etabar_i - etabar_j)^{2} + 2*etabar_i*etabar_j(1-cos(phi_i - phi_j))]/R^{2}
+
+ void init(const PseudoJet & jet, CentauroInfo * info) {
+
+ R = info->R();
+ gammaE = info->gammaE();
+ gammaPz = info->gammaPz();
+
+ // photon 4-momentum is q = (gammaE, 0 , 0 , gammaPz);
+
+ 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;
+ pT = jet.perp();
+ phi = jet.phi();
+ if(gammaE!=0 and gammaPz!=0){ //gammaE and gammaPz passed, so not running in Breit frame
+ Q = sqrt(-1.0*(gammaE*gammaE-gammaPz*gammaPz));
+ etabar = -2.0*(Q/(gammaE+gammaPz))*(pT/(jet.E()-jet.pz()));
+ }
+ else{ //gammaE and gammaPz not passed, so assume that it is running in the Breit frame
+ etabar = +2.0*pT/(jet.E()-jet.pz());
+ }
+ // beam distance
+ diB = 1.0;
+ }
+
+ double distance(const CentauroBriefJet * jet) const {
+
+ double dij = pow(etabar - jet->etabar, 2.0) + 2*etabar*jet->etabar*(1-cos(phi- jet->phi));
+ dij = dij/pow(R,2.0);
+
+ return dij;
+ }
+
+ double beam_distance() const {
+ return diB;
+ }
+
+ double pT, phi, nx, ny, nz;
+ double etabar;
+ double diB;
+ double R, gammaE, gammaPz, Q;
+ };
+
+
+ std::string CentauroPlugin::description () const {
+ std::ostringstream desc;
+ desc << "Centauro plugin with R = " << R();
+ if(gammaE()==0 and gammaPz()==0){
+ desc << " gamma E and gamma Pz parameters were not given --> assume you are giving particles momenta in Breit frame";
+ }
+ return desc.str();
+ }
+
+ void CentauroPlugin::run_clustering(fastjet::ClusterSequence & cs) const {
+ int njets = cs.jets().size();
+ CentauroInfo vinfo(R(), gammaE(), gammaPz());
+
+ NNH 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
Property changes on: contrib/contribs/Centauro/tags/1.0.0/Centauro.cc
___________________________________________________________________
Added: svn:keywords
## -0,0 +1 ##
+Id
\ No newline at end of property
Index: contrib/contribs/Centauro/tags/1.0.0/AUTHORS
===================================================================
--- contrib/contribs/Centauro/tags/1.0.0/AUTHORS (revision 0)
+++ contrib/contribs/Centauro/tags/1.0.0/AUTHORS (revision 1278)
@@ -0,0 +1,5 @@
+The Centauro jet algorithm is described in arXiv:2006.10751, "Asymmetric jet clustering in deep-inelastic scattering",
+
+by Miguel Arratia, Yiannis Makris, Duff Neill, Felix Ringer, Nobuo Sato.
+
+Coded in the format of FastJet contrib classes by M. Arratia.
Index: contrib/contribs/Centauro/tags/1.0.0/VERSION
===================================================================
--- contrib/contribs/Centauro/tags/1.0.0/VERSION (revision 0)
+++ contrib/contribs/Centauro/tags/1.0.0/VERSION (revision 1278)
@@ -0,0 +1 @@
+1.0.0
Index: contrib/contribs/Centauro/tags/1.0.0/example.ref
===================================================================
--- contrib/contribs/Centauro/tags/1.0.0/example.ref (revision 0)
+++ contrib/contribs/Centauro/tags/1.0.0/example.ref (revision 1278)
@@ -0,0 +1,471 @@
+Particles that are input
+-0.880741 -1.23313 -157.431 157.439
+-0.00517126 0.238155 -9.7396 9.74352
+0.0362281 0.269475 -6.92434 6.93108
+-0.220663 -0.14382 -0.683861 0.746004
+1.27168 1.04223 -6.17402 6.39073
+-0.569559 -0.362776 -58.543 58.5545
+0.283999 -0.46682 -49.6979 49.7097
+0.651053 1.39709 -62.7226 62.7486
+0.143456 -0.0312881 -6.93824 6.94119
+0.493156 2.16278 -14.8866 15.0516
+0.239681 -0.0786237 -1.9341 1.95546
+0.335549 0.0516403 -0.834654 0.911804
+-0.785387 -0.781052 -1.53678 1.89949
+0.109489 -0.175467 -0.684313 0.728382
+-1.33953 -1.06775 -6.45274 6.74204
+-0.461043 -0.0168194 -0.835725 0.964759
+-0.858633 -2.02526 -6.11693 6.50194
+0.150685 -0.281284 -0.303887 0.462226
+1.99802 1.13 -3.44596 4.24547
+0.327288 -0.236459 0.155347 0.45458
+1.00764 1.44403 0.852359 1.96126
+0.0379999 -0.36601 0.152872 0.422205
+10.0152 -2.32116 -15.0786 18.2503
+11.3599 -1.23253 -15.8349 19.5277
+1.66877 -1.15858 -2.87099 3.51982
+3.20144 -0.3128 -0.551375 3.3958
+-30.7151 7.18282 -30.3287 43.7617
+-7.59561 2.59734 -7.53885 11.0235
+-0.92136 -0.510777 -0.14196 1.07212
+-1.07946 -0.66682 0.531688 1.38277
+-0.739634 -0.618217 2.51194 2.69418
+0.207418 0.279428 0.996016 1.06425
+-0.153917 -0.775168 2.57381 2.73728
+-1.32156 -0.637192 234.435 234.44
+-1.05753 0.0953188 464.947 464.948
+-0.627641 0.236196 71.0829 71.0862
+-0.208202 -0.437336 0.902859 1.03404
+-0.0179792 0.0902484 0.264846 0.313195
+0.10202 0.109801 1.77874 1.79049
+2.07066 2.40441 1.0424 3.46962
+0.797326 0.534987 57.3601 57.3683
+-0.209682 0.0247006 -23.7891 23.7952
+0.310251 0.131903 -314.235 314.236
+-0.255101 0.141414 -264.421 264.421
+-0.0667217 0.167567 -440.299 440.299
+0.123081 0.0544788 -9.23343 9.23441
+0.0149245 0.0926409 -13.1527 13.1531
+-0.727346 0.810304 -69.3094 69.3181
+1.34449 -0.0354347 -75.6843 75.6963
+-0.267367 -0.0908127 -12.8334 12.8373
+0.347806 0.115655 -5.5306 5.54274
+0.106686 -0.019479 -1.03504 1.0407
+0.0216701 -0.0151984 -0.0808218 0.0850455
+-0.0908971 -0.179593 -2.14713 2.15654
+0.16191 0.109639 -5.15499 5.16059
+0.16588 -0.432303 -10.972 10.9928
+-0.0738683 0.218955 -5.28589 5.29278
+0.603248 0.0236333 -3.20082 3.29444
+-0.0357684 0.185609 -0.623546 0.651567
+-0.353371 0.62957 -2.86075 2.95045
+0.0481173 0.0105768 -0.5957 0.597734
+0.207583 -0.0308763 -0.74871 0.777567
+0.530402 1.47778 -3.93222 4.33708
+0.0413924 0.148665 -0.0938036 0.22824
+0.0836184 1.01756 -5.1672 5.35025
+0.348582 0.740061 -2.32778 2.47128
+0.28408 0.144808 -0.359861 0.50065
+0.0615749 0.122668 -0.256791 0.322894
+0.444708 -0.470324 -41.1153 41.1207
+0.0124634 0.401554 -14.9711 14.9771
+-0.031622 -0.143397 -40.0395 40.0428
+-0.308748 0.180695 -13.8308 13.8354
+-0.335663 0.0993823 -10.08 10.0861
+-0.320188 -0.202498 -41.8447 41.8569
+0.106188 -0.315185 -14.3999 14.4044
+-0.252772 0.906934 -58.5239 58.5316
+-0.179212 -0.0360309 -33.2666 33.2674
+0.96756 0.359378 -47.8071 47.8275
+-0.0714174 0.0908505 -6.71525 6.71769
+-0.0778202 1.40326 -90.8847 90.9004
+0.0244495 0.0548209 -26.7828 26.7832
+-0.166031 -0.0217473 -18.6254 18.6267
+0.015705 0.0168585 -0.320317 0.321145
+0.310761 -0.197165 -2.52162 2.54834
+0.0555181 0.101241 -2.85688 2.85921
+0.492405 0.463807 -11.5989 11.6186
+0.0246279 0.350083 -1.47621 1.52376
+-0.13378 0.0207567 -0.692298 0.705411
+0.142599 0.215925 -1.44223 1.46525
+1.15091 1.39326 -11.9161 12.0524
+-0.281736 0.117135 -3.50441 3.52043
+-0.29383 0.547755 -4.06079 4.11046
+0.232077 0.272558 -2.44178 2.47183
+-0.116685 0.356148 -1.01678 1.09261
+0.430131 0.339061 -3.71537 3.75812
+0.229221 0.534237 -4.02785 4.07198
+-0.952433 0.033435 -3.68614 3.80734
+-0.265694 -0.257273 -2.34127 2.3703
+-1.03114 -1.05862 -4.07097 4.43165
+-0.0950318 0.0651839 -0.470922 0.504507
+0.0664907 -0.228339 -1.87629 1.89645
+-0.561111 -0.695389 -1.96068 2.1592
+0.189189 -0.147209 -0.209635 0.318449
+0.0792705 0.0251845 -0.0983276 0.128788
+0.161174 0.0711834 -0.369058 0.43212
+-0.00547856 -0.177247 -0.188786 0.294222
+1.25928 0.602768 -1.97317 2.59285
+0.202956 0.165086 0.0141931 0.29686
+1.2749 0.41628 0.264907 1.37416
+1.07623 0.331172 0.552579 1.26205
+0.161704 -0.119772 0.0749213 0.256099
+-0.0323955 0.249423 0.156091 0.32727
+0.478365 0.0155339 0.0730681 0.484163
+0.888763 0.0859887 0.317118 0.947554
+0.462039 0.0959242 0.95491 1.07425
+0.492961 0.292406 0.404665 0.715364
+0.809038 0.335743 0.124281 0.895651
+-0.221861 0.0554978 0.114004 0.291168
+1.09472 -0.623643 -0.882151 1.54435
+1.23124 -0.738088 -1.69685 2.22699
+0.982367 -0.198001 -1.1324 1.51857
+0.307264 -0.361637 -0.886308 1.01499
+3.52671 -0.870934 -6.0396 7.04929
+0.030231 -0.109307 -0.145418 0.231275
+8.31951 -1.26965 -13.6925 16.0727
+1.91112 -0.260403 -2.87686 3.46641
+0.530177 -0.247449 -1.16457 1.31073
+1.74203 -0.282291 -2.81689 3.32404
+7.16624 -0.898926 -11.6921 13.7429
+4.02785 -0.75388 -6.13676 7.43872
+2.35206 -0.359454 -3.91487 4.58335
+0.887401 0.111805 0.0320276 0.905807
+0.150072 0.0233403 0.0493208 0.212082
+4.29627 -1.71423 -0.169014 4.63083
+1.84307 -0.435415 0.398654 1.93531
+17.9175 -3.95419 4.19095 18.8211
+3.73811 -1.15883 1.28947 4.12293
+40.1307 -9.94554 9.43724 42.4084
+46.9323 -12.1458 11.2984 49.7778
+49.9481 -13.6951 11.6388 53.0835
+55.4632 -13.9171 13.2275 58.6926
+8.22602 -2.01464 1.97951 8.69739
+552.609 -143.384 127.406 584.951
+-58.5643 13.6721 -60.7292 85.4679
+-20.7056 4.93672 -21.3461 30.1458
+-63.1393 15.4319 -63.1699 90.6389
+-184.709 43.8836 -185.99 265.774
+-132.344 31.1976 -133.812 190.771
+-249.626 59.9807 -252.053 359.78
+-84.6956 20.2531 -84.9406 121.649
+-1.48961 0.539472 -1.16977 1.97428
+-15.691 4.0929 -16.0309 22.8077
+-6.51615 1.45304 -7.04142 9.70426
+-1.68878 0.621736 -1.74539 2.51086
+-2.9066 0.85783 -2.17682 3.73393
+-4.32166 0.487667 -4.01063 5.91771
+-0.131328 0.171002 0.296093 0.366278
+-0.13788 0.0239116 0.133612 0.193481
+0.122857 0.103882 0.18563 0.24565
+0.0880826 -0.0369855 0.0702433 0.118577
+0.300034 -0.0143367 0.819116 0.883548
+0.110318 -0.327697 0.569549 0.680751
+0.024617 -0.0632054 0.446281 0.451407
+0.346265 -0.137049 1.8336 1.87103
+0.107272 0.234714 0.810601 0.983521
+-0.122867 -0.342083 0.426945 0.577822
+-0.0537484 -0.141832 0.165822 0.224726
+-0.192158 -0.205106 0.651038 0.722719
+0.0719282 -0.638981 1.32607 1.48034
+0.0770785 -0.0961695 0.712638 0.736561
+-0.102479 0.25704 1.00071 1.04761
+-0.356965 0.0132487 8.272 8.28088
+0.911974 -0.544036 10.4508 10.5046
+0.213035 0.082474 12.4439 12.4814
+0.0796437 -0.219034 1.895 1.91437
+0.0257845 -0.0205888 0.413465 0.41478
+0.340701 -0.977123 20.8001 20.8258
+-0.0484682 -0.417098 51.2421 51.2524
+0.0897243 -0.147796 33.5942 33.595
+0.330917 0.205356 15.3759 15.3815
+-0.451343 -1.06297 198.648 198.652
+-0.509621 0.280855 526.283 526.284
+0.043921 -0.0909597 62.6116 62.6118
+-0.492147 -0.225927 263.163 263.165
+0.0256742 0.0647083 60.1293 60.1295
+-0.212947 0.335618 10.9531 10.9612
+-0.334505 0.00530422 10.0385 10.0451
+-0.216888 0.0832963 14.9197 14.9221
+0.219073 0.140813 1.5508 1.57869
+0.115749 0.0973728 0.277607 0.345579
+-0.270018 -0.265979 3.81164 3.83044
+-0.050352 -0.0502248 1.59578 1.59736
+0.372829 -0.0813045 0.383384 0.558636
+0.204553 0.553167 1.40673 1.53173
+1.11246 0.636345 0.616124 1.70366
+-0.00457604 -0.0333414 -0.0266276 0.146019
+0.203473 -0.097539 -0.0758289 0.275944
+0.219887 0.595547 0.725965 0.974439
+1.65124 1.32091 1.55523 2.62862
+0.148563 0.330211 2.28693 2.49879
+0.143486 0.174823 1.25361 1.28147
+0.446611 -0.370857 3.34549 3.52308
+0.0645162 -0.324011 0.86536 0.936735
+0.729718 -0.023067 8.75816 8.78964
+0.0818687 -1.05908 6.0449 6.13752
+-0.0165521 0.049697 0.93817 0.939631
+0.185679 -0.265143 3.52967 3.54723
+0.24623 -0.0254016 9.58382 9.58803
+0.456426 0.154724 36.4834 36.49
+0.242001 -0.475077 8.89784 8.91489
+-0.583645 0.271306 97.8985 97.9007
+0.109873 -0.317333 192.152 192.153
+-0.0884675 -0.0343951 37.9163 37.9166
+-0.347987 -1.61636 3219.7 3219.7
+-0.450278 -0.0845034 -15.5305 15.5379
+-0.0136189 0.0749919 -11.8648 11.8658
+-0.235612 0.254626 -464.623 464.623
+-0.336693 0.000126715 -314.791 314.791
+-0.590802 0.271431 -62.057 62.0624
+0.0950294 -0.0979437 -6.50751 6.50894
+-0.0258241 -0.16519 -7.66134 7.66316
+-0.0283747 -0.0128187 -1.18609 1.1865
+-0.244547 -0.036785 -2.3235 2.33662
+1.11342 -0.00849554 -8.17565 8.25113
+0.742307 0.0450352 -4.78039 4.83789
+0.0631615 -0.0147898 -0.337208 0.34339
+-0.0615524 -0.073034 -0.410951 0.421904
+0.272176 0.935645 -1.18993 1.538
+0.0195568 -0.0613884 -0.243814 0.252183
+0.00157538 -0.0127041 -5.17987 5.17988
+-0.23365 0.239859 -61.5028 61.5037
+-0.211158 -1.22643 -54.4484 54.4628
+-0.172415 -0.671014 -46.7073 46.7127
+-0.00835489 0.119193 -1.62732 1.6317
+0.178258 0.23706 -3.45154 3.46426
+-0.015716 -0.0172284 -0.832185 0.832512
+0.205512 -0.227745 -4.74005 4.74996
+0.1434 0.131935 -4.22377 4.22826
+0.250522 0.510955 -14.5866 14.5977
+0.0907355 0.221224 -1.62452 1.64202
+0.0442827 0.68279 -4.27125 4.32571
+-0.02869 0.127284 -0.95867 0.967508
+0.0168298 0.0212463 -1.24988 1.25017
+0.0191627 0.0396505 -0.109601 0.118118
+-0.182616 0.158348 -0.349358 0.424823
+0.0736987 -0.47576 -0.556793 0.736069
+0.0718737 -0.10731 -0.092269 0.158729
+0.382305 0.247318 -0.976494 1.07743
+-0.00996699 0.0144861 -0.0618196 0.0642717
+-0.089312 0.0045092 -0.421947 0.43132
+-0.0173901 0.0980487 -0.228073 0.248864
+0.0470684 0.106722 -0.0220929 0.118714
+0.223635 0.146211 -0.200759 0.334207
+0.230494 -0.178628 0.27383 0.400023
+-0.0034267 -0.0207838 -0.00992506 0.0232856
+-0.0198611 0.0339168 0.0860578 0.0946084
+-0.149509 -0.0992705 0.172169 0.248696
+0.138859 0.0621639 0.18205 0.237252
+0.436405 0.171141 0.252614 0.532497
+0.113578 -0.375531 -0.481926 0.621431
+0.0656587 -0.117427 -0.310663 0.338543
+0.679416 -0.295888 -1.31087 1.50583
+0.210556 -0.0146095 -0.359833 0.417166
+0.889197 -0.285985 -1.28992 1.59259
+0.970963 -0.463859 -1.50072 1.84664
+0.837859 -0.0907723 -1.08119 1.37084
+0.177904 -0.0469476 -0.153921 0.239886
+0.469609 -0.218861 -0.50301 0.722116
+0.451137 -0.114036 -0.591412 0.752527
+1.0715 -0.0504322 -1.25983 1.65464
+0.0668091 -0.0378517 -0.0814591 0.111945
+0.631639 -0.127324 -0.646231 0.912575
+0.268405 0.00839595 -0.188646 0.328175
+0.389287 -0.174337 -0.650956 0.778256
+1.33893 -0.36273 -1.92186 2.3702
+0.00450944 -0.0453765 -0.0434901 0.0630138
+0.313353 -0.100881 -0.365535 0.491917
+1.45193 -0.377662 -2.34366 2.78271
+1.17749 -0.214244 -1.72329 2.09812
+0.235585 -0.220075 0.140179 0.351544
+0.373054 -0.203334 0.0678202 0.430249
+5.70366 -1.38471 1.03047 5.95911
+0.313298 -0.0441957 0.0509067 0.320469
+12.0524 -3.08888 2.78043 12.7488
+31.6297 -7.9897 7.48225 33.4703
+0.129817 -0.0603423 0.0449213 0.150038
+5.60412 -1.78273 1.39294 6.04356
+15.2636 -3.91265 3.52274 16.146
+12.642 -3.31563 2.82037 13.3704
+-6.82806 1.65752 -7.06422 9.96359
+-28.7311 6.68994 -29.6661 41.8367
+-23.0806 5.28047 -23.7594 33.5429
+-9.58673 2.35668 -9.63028 13.7921
+-2.49416 0.760594 -2.34655 3.50794
+-12.8361 4.17247 -12.393 18.3238
+-0.221894 0.0177591 -0.160145 0.274224
+-0.354718 0.16272 -0.437598 0.586341
+-1.50778 0.484668 -1.37224 2.09555
+-1.02491 0.345425 -1.0884 1.5344
+-0.685101 0.138786 -0.693147 0.984417
+-0.134848 0.00278312 -0.0667383 0.150485
+-1.71723 0.427862 -1.54836 2.35146
+-0.118508 -0.00234705 -0.125537 0.172653
+-0.248198 0.158231 -0.146802 0.328923
+-0.996664 0.338963 -0.478314 1.15629
+0.00883286 -0.0568652 0.366604 0.371093
+-0.0824728 0.0441798 0.387952 0.399075
+0.119963 -0.0647633 0.301229 0.330642
+0.00186304 -0.004997 -0.0147955 0.0157273
+-0.0136056 0.00781108 0.00427791 0.0162612
+0.0895633 -0.121489 0.551825 0.572094
+0.309377 0.0801079 4.73961 4.75037
+0.144343 0.0774488 4.14383 4.14707
+0.0800207 -0.283154 44.4007 44.4016
+0.00773677 0.0117046 0.771776 0.771903
+-0.111439 -0.0611335 10.4091 10.4099
+-0.21797 0.0686322 20.7719 20.7732
+-0.0517191 0.14414 13.4281 13.429
+-0.0315432 -0.0288834 4.11543 4.11565
+-0.021049 -0.191329 14.6251 14.6264
+-0.413306 -0.994221 74.9928 75.0006
+0.00171897 0.000308655 0.116277 0.11629
+-0.0197979 0.491867 1.49059 1.56977
+0.149821 0.265803 0.0262281 0.306244
+0.148572 0.0742812 -0.00738766 0.166271
+-0.00508379 0.000228477 -0.00219567 0.00554239
+0.190522 0.072154 1.04008 1.05984
+0.17894 0.0666765 2.90455 2.91082
+0.194374 -0.0680286 3.01187 3.01891
+0.0150257 0.016201 3.84059 3.84066
+0.0864382 -0.177535 10.9517 10.9534
+0.0112701 0.00776215 6.26456 6.26457
+-0.088069 -0.201175 24.3997 24.4007
+-0.0428591 0.0438332 44.4944 44.4945
+0.0444439 -0.049705 94.1332 94.1332
+0.396364 0.067377 -289.884 289.884
+0.0542375 0.0178424 -96.5357 96.5358
+-0.0134027 -0.0410107 -5.27787 5.27805
+-0.00741191 -0.0297613 -58.9037 58.9037
+-0.281117 0.856545 -1345.21 1345.21
+-0.0217591 0.0295681 -1.43166 1.43214
+-0.0411872 -0.143315 -4.19509 4.19773
+1.61243e-05 -0.0197893 -18.6653 18.6653
+-0.00370808 -0.164426 -20.7847 20.7853
+0.186445 0.0253239 0.0515968 0.195103
+0.103622 -0.0670424 -0.034011 0.128019
+0.467247 0.223657 -0.0501274 0.520438
+0.189184 0.0409518 0.0506236 0.200076
+0.0261106 0.0529553 -0.0657671 0.0883818
+0.161813 0.0242326 0.00236871 0.163634
+0.356526 0.012488 9.39339 9.4012
+0.50335 -0.50074 16.8455 16.861
+-0.13982 0.0757052 -118.92 118.92
+0.00182717 -0.0111141 -74.4315 74.4315
+#--------------------------------------------------------------------------
+# FastJet release 3.3.4
+# M. Cacciari, G.P. Salam and G. Soyez
+# A software package for jet finding and analysis at colliders
+# http://fastjet.fr
+#
+# Please cite EPJC72(2012)1896 [arXiv:1111.6097] if you use this package
+# for scientific work and optionally PLB641(2006)57 [hep-ph/0512210].
+#
+# FastJet is provided without warranty under the GNU GPL v2 or higher.
+# It uses T. Chan's closest pair algorithm, S. Fortune's Voronoi code
+# and 3rd party plugin jet algorithms. See COPYING file for details.
+#--------------------------------------------------------------------------
+jets in inclusive clustering
+ rap = -1.6988 e 6851.82 n constituents 208
+ rap = 8.12777 e 3219.7 n constituents 1
+ rap = 0.220458 e 919.866 n constituents 22
+ rap = 6.86 e 526.284 n constituents 1
+ rap = 6.76652 e 464.948 n constituents 1
+ rap = 6.18589 e 263.165 n constituents 1
+ rap = 5.76251 e 234.44 n constituents 1
+ rap = 5.83349 e 198.652 n constituents 1
+ rap = 6.96296 e 192.153 n constituents 1
+ rap = 5.69476 e 97.9007 n constituents 1
+ rap = 7.94575 e 94.1332 n constituents 1
+ rap = 4.93668 e 75.0006 n constituents 1
+ rap = 5.33538 e 71.0862 n constituents 1
+ rap = 6.5887 e 62.6118 n constituents 1
+ rap = 6.64775 e 60.1295 n constituents 1
+ rap = 4.77275 e 57.3683 n constituents 1
+ rap = 4.60225 e 51.2524 n constituents 1
+ rap = 7.28041 e 44.4945 n constituents 1
+ rap = 5.70976 e 44.4016 n constituents 1
+ rap = 6.10765 e 37.9166 n constituents 1
+ rap = 4.65714 e 36.49 n constituents 1
+ rap = 5.71168 e 33.595 n constituents 1
+ rap = 5.40365 e 24.4007 n constituents 1
+ rap = 3.6945 e 20.8258 n constituents 1
+ rap = 5.20291 e 20.7732 n constituents 1
+ rap = 3.84122 e 16.861 n constituents 1
+ rap = 4.30872 e 15.3815 n constituents 1
+ rap = 4.70143 e 14.9221 n constituents 1
+ rap = 5.02367 e 14.6264 n constituents 1
+ rap = 5.16695 e 13.429 n constituents 1
+ rap = 3.2495 e 12.4814 n constituents 1
+ rap = 3.95163 e 10.9612 n constituents 1
+ rap = 4.70894 e 10.9534 n constituents 1
+ rap = 2.98232 e 10.5046 n constituents 1
+ rap = 5.0986 e 10.4099 n constituents 1
+ rap = 4.01465 e 10.0451 n constituents 1
+ rap = 0.24229 e 9.86 n constituents 15
+ rap = 4.21162 e 9.58803 n constituents 1
+ rap = 3.89309 e 9.4012 n constituents 1
+ rap = 0.597202 e 9.33815 n constituents 6
+ rap = 3.4757 e 8.91489 n constituents 1
+ rap = 3.16158 e 8.78964 n constituents 1
+ rap = 3.76495 e 8.28088 n constituents 1
+ rap = 6.81954 e 6.26457 n constituents 1
+ rap = 2.43961 e 6.13752 n constituents 1
+ rap = 3.39098 e 4.75037 n constituents 1
+ rap = 3.92422 e 4.14707 n constituents 1
+ rap = 5.25985 e 4.11565 n constituents 1
+ rap = 5.85113 e 3.84066 n constituents 1
+ rap = 3.00384 e 3.83044 n constituents 1
+ rap = 0.703275 e 3.68956 n constituents 3
+ rap = 2.99954 e 3.54723 n constituents 1
+ rap = 1.82762 e 3.52308 n constituents 1
+ rap = 3.37707 e 3.01891 n constituents 1
+ rap = 3.4162 e 2.91082 n constituents 1
+ rap = 1.74045 e 2.73728 n constituents 1
+ rap = 1.67615 e 2.69418 n constituents 1
+ rap = 2.64067 e 1.91437 n constituents 1
+ rap = 2.29737 e 1.87103 n constituents 1
+ rap = 2.85806 e 1.79049 n constituents 1
+ rap = 3.80441 e 1.59736 n constituents 1
+ rap = 2.36007 e 1.57869 n constituents 1
+ rap = 1.82727 e 1.56977 n constituents 1
+ rap = 1.57865 e 1.53173 n constituents 1
+ rap = 1.45048 e 1.48034 n constituents 1
+ rap = 0.40534 e 1.38277 n constituents 1
+ rap = 2.25537 e 1.28147 n constituents 1
+ rap = 1.4167 e 1.07425 n constituents 1
+ rap = -0.133193 e 1.07212 n constituents 1
+ rap = 1.70382 e 1.06425 n constituents 1
+ rap = 2.33287 e 1.05984 n constituents 1
+ rap = 1.88846 e 1.04761 n constituents 1
+ rap = 1.34612 e 1.03404 n constituents 1
+ rap = 0.961642 e 0.974439 n constituents 1
+ rap = 3.57931 e 0.939631 n constituents 1
+ rap = 1.61438 e 0.936735 n constituents 1
+ rap = 1.63717 e 0.883548 n constituents 1
+ rap = 0.947404 e 0.802549 n constituents 2
+ rap = 4.70069 e 0.771903 n constituents 1
+ rap = 2.05196 e 0.736561 n constituents 1
+ rap = 1.47654 e 0.722719 n constituents 1
+ rap = 1.2099 e 0.680751 n constituents 1
+ rap = 0.810459 e 0.677214 n constituents 2
+ rap = 0.787932 e 0.640465 n constituents 2
+ rap = 2.00773 e 0.572094 n constituents 1
+ rap = 2.58282 e 0.451407 n constituents 1
+ rap = 0.379276 e 0.422205 n constituents 1
+ rap = 3.22292 e 0.41478 n constituents 1
+ rap = 0.837601 e 0.400023 n constituents 1
+ rap = 2.12965 e 0.399075 n constituents 1
+ rap = 2.55093 e 0.371093 n constituents 1
+ rap = 1.12235 e 0.366278 n constituents 1
+ rap = 0.422165 e 0.351544 n constituents 1
+ rap = 1.53362 e 0.330642 n constituents 1
+ rap = 0.405615 e 0.307429 n constituents 2
+ rap = 0.852332 e 0.248696 n constituents 1
+ rap = 0.98604 e 0.24565 n constituents 1
+ rap = 1.0138 e 0.237252 n constituents 1
+ rap = 0.849045 e 0.193481 n constituents 1
+ rap = 4.89158 e 0.11629 n constituents 1
+ rap = 1.52532 e 0.0946084 n constituents 1
+
Index: contrib/contribs/Centauro/tags/1.0.0/COPYING
===================================================================
--- contrib/contribs/Centauro/tags/1.0.0/COPYING (revision 0)
+++ contrib/contribs/Centauro/tags/1.0.0/COPYING (revision 1278)
@@ -0,0 +1,343 @@
+======================================================================
+======================================================================
+======================================================================
+ GNU GENERAL PUBLIC LICENSE
+ Version 2, June 1991
+
+ Copyright (C) 1989, 1991 Free Software Foundation, Inc.
+ 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+ Everyone is permitted to copy and distribute verbatim copies
+ of this license document, but changing it is not allowed.
+
+ Preamble
+
+ The licenses for most software are designed to take away your
+freedom to share and change it. By contrast, the GNU General Public
+License is intended to guarantee your freedom to share and change free
+software--to make sure the software is free for all its users. This
+General Public License applies to most of the Free Software
+Foundation's software and to any other program whose authors commit to
+using it. (Some other Free Software Foundation software is covered by
+the GNU Library General Public License instead.) You can apply it to
+your programs, too.
+
+ When we speak of free software, we are referring to freedom, not
+price. Our General Public Licenses are designed to make sure that you
+have the freedom to distribute copies of free software (and charge for
+this service if you wish), that you receive source code or can get it
+if you want it, that you can change the software or use pieces of it
+in new free programs; and that you know you can do these things.
+
+ To protect your rights, we need to make restrictions that forbid
+anyone to deny you these rights or to ask you to surrender the rights.
+These restrictions translate to certain responsibilities for you if you
+distribute copies of the software, or if you modify it.
+
+ For example, if you distribute copies of such a program, whether
+gratis or for a fee, you must give the recipients all the rights that
+you have. You must make sure that they, too, receive or can get the
+source code. And you must show them these terms so they know their
+rights.
+
+ We protect your rights with two steps: (1) copyright the software, and
+(2) offer you this license which gives you legal permission to copy,
+distribute and/or modify the software.
+
+ Also, for each author's protection and ours, we want to make certain
+that everyone understands that there is no warranty for this free
+software. If the software is modified by someone else and passed on, we
+want its recipients to know that what they have is not the original, so
+that any problems introduced by others will not reflect on the original
+authors' reputations.
+
+ Finally, any free program is threatened constantly by software
+patents. We wish to avoid the danger that redistributors of a free
+program will individually obtain patent licenses, in effect making the
+program proprietary. To prevent this, we have made it clear that any
+patent must be licensed for everyone's free use or not licensed at all.
+
+ The precise terms and conditions for copying, distribution and
+modification follow.
+
+ GNU GENERAL PUBLIC LICENSE
+ TERMS AND CONDITIONS FOR COPYING, DISTRIBUTION AND MODIFICATION
+
+ 0. This License applies to any program or other work which contains
+a notice placed by the copyright holder saying it may be distributed
+under the terms of this General Public License. The "Program", below,
+refers to any such program or work, and a "work based on the Program"
+means either the Program or any derivative work under copyright law:
+that is to say, a work containing the Program or a portion of it,
+either verbatim or with modifications and/or translated into another
+language. (Hereinafter, translation is included without limitation in
+the term "modification".) Each licensee is addressed as "you".
+
+Activities other than copying, distribution and modification are not
+covered by this License; they are outside its scope. The act of
+running the Program is not restricted, and the output from the Program
+is covered only if its contents constitute a work based on the
+Program (independent of having been made by running the Program).
+Whether that is true depends on what the Program does.
+
+ 1. You may copy and distribute verbatim copies of the Program's
+source code as you receive it, in any medium, provided that you
+conspicuously and appropriately publish on each copy an appropriate
+copyright notice and disclaimer of warranty; keep intact all the
+notices that refer to this License and to the absence of any warranty;
+and give any other recipients of the Program a copy of this License
+along with the Program.
+
+You may charge a fee for the physical act of transferring a copy, and
+you may at your option offer warranty protection in exchange for a fee.
+
+ 2. You may modify your copy or copies of the Program or any portion
+of it, thus forming a work based on the Program, and copy and
+distribute such modifications or work under the terms of Section 1
+above, provided that you also meet all of these conditions:
+
+ a) You must cause the modified files to carry prominent notices
+ stating that you changed the files and the date of any change.
+
+ b) You must cause any work that you distribute or publish, that in
+ whole or in part contains or is derived from the Program or any
+ part thereof, to be licensed as a whole at no charge to all third
+ parties under the terms of this License.
+
+ c) If the modified program normally reads commands interactively
+ when run, you must cause it, when started running for such
+ interactive use in the most ordinary way, to print or display an
+ announcement including an appropriate copyright notice and a
+ notice that there is no warranty (or else, saying that you provide
+ a warranty) and that users may redistribute the program under
+ these conditions, and telling the user how to view a copy of this
+ License. (Exception: if the Program itself is interactive but
+ does not normally print such an announcement, your work based on
+ the Program is not required to print an announcement.)
+
+These requirements apply to the modified work as a whole. If
+identifiable sections of that work are not derived from the Program,
+and can be reasonably considered independent and separate works in
+themselves, then this License, and its terms, do not apply to those
+sections when you distribute them as separate works. But when you
+distribute the same sections as part of a whole which is a work based
+on the Program, the distribution of the whole must be on the terms of
+this License, whose permissions for other licensees extend to the
+entire whole, and thus to each and every part regardless of who wrote it.
+
+Thus, it is not the intent of this section to claim rights or contest
+your rights to work written entirely by you; rather, the intent is to
+exercise the right to control the distribution of derivative or
+collective works based on the Program.
+
+In addition, mere aggregation of another work not based on the Program
+with the Program (or with a work based on the Program) on a volume of
+a storage or distribution medium does not bring the other work under
+the scope of this License.
+
+ 3. You may copy and distribute the Program (or a work based on it,
+under Section 2) in object code or executable form under the terms of
+Sections 1 and 2 above provided that you also do one of the following:
+
+ a) Accompany it with the complete corresponding machine-readable
+ source code, which must be distributed under the terms of Sections
+ 1 and 2 above on a medium customarily used for software interchange; or,
+
+ b) Accompany it with a written offer, valid for at least three
+ years, to give any third party, for a charge no more than your
+ cost of physically performing source distribution, a complete
+ machine-readable copy of the corresponding source code, to be
+ distributed under the terms of Sections 1 and 2 above on a medium
+ customarily used for software interchange; or,
+
+ c) Accompany it with the information you received as to the offer
+ to distribute corresponding source code. (This alternative is
+ allowed only for noncommercial distribution and only if you
+ received the program in object code or executable form with such
+ an offer, in accord with Subsection b above.)
+
+The source code for a work means the preferred form of the work for
+making modifications to it. For an executable work, complete source
+code means all the source code for all modules it contains, plus any
+associated interface definition files, plus the scripts used to
+control compilation and installation of the executable. However, as a
+special exception, the source code distributed need not include
+anything that is normally distributed (in either source or binary
+form) with the major components (compiler, kernel, and so on) of the
+operating system on which the executable runs, unless that component
+itself accompanies the executable.
+
+If distribution of executable or object code is made by offering
+access to copy from a designated place, then offering equivalent
+access to copy the source code from the same place counts as
+distribution of the source code, even though third parties are not
+compelled to copy the source along with the object code.
+
+ 4. You may not copy, modify, sublicense, or distribute the Program
+except as expressly provided under this License. Any attempt
+otherwise to copy, modify, sublicense or distribute the Program is
+void, and will automatically terminate your rights under this License.
+However, parties who have received copies, or rights, from you under
+this License will not have their licenses terminated so long as such
+parties remain in full compliance.
+
+ 5. You are not required to accept this License, since you have not
+signed it. However, nothing else grants you permission to modify or
+distribute the Program or its derivative works. These actions are
+prohibited by law if you do not accept this License. Therefore, by
+modifying or distributing the Program (or any work based on the
+Program), you indicate your acceptance of this License to do so, and
+all its terms and conditions for copying, distributing or modifying
+the Program or works based on it.
+
+ 6. Each time you redistribute the Program (or any work based on the
+Program), the recipient automatically receives a license from the
+original licensor to copy, distribute or modify the Program subject to
+these terms and conditions. You may not impose any further
+restrictions on the recipients' exercise of the rights granted herein.
+You are not responsible for enforcing compliance by third parties to
+this License.
+
+ 7. If, as a consequence of a court judgment or allegation of patent
+infringement or for any other reason (not limited to patent issues),
+conditions are imposed on you (whether by court order, agreement or
+otherwise) that contradict the conditions of this License, they do not
+excuse you from the conditions of this License. If you cannot
+distribute so as to satisfy simultaneously your obligations under this
+License and any other pertinent obligations, then as a consequence you
+may not distribute the Program at all. For example, if a patent
+license would not permit royalty-free redistribution of the Program by
+all those who receive copies directly or indirectly through you, then
+the only way you could satisfy both it and this License would be to
+refrain entirely from distribution of the Program.
+
+If any portion of this section is held invalid or unenforceable under
+any particular circumstance, the balance of the section is intended to
+apply and the section as a whole is intended to apply in other
+circumstances.
+
+It is not the purpose of this section to induce you to infringe any
+patents or other property right claims or to contest validity of any
+such claims; this section has the sole purpose of protecting the
+integrity of the free software distribution system, which is
+implemented by public license practices. Many people have made
+generous contributions to the wide range of software distributed
+through that system in reliance on consistent application of that
+system; it is up to the author/donor to decide if he or she is willing
+to distribute software through any other system and a licensee cannot
+impose that choice.
+
+This section is intended to make thoroughly clear what is believed to
+be a consequence of the rest of this License.
+
+ 8. If the distribution and/or use of the Program is restricted in
+certain countries either by patents or by copyrighted interfaces, the
+original copyright holder who places the Program under this License
+may add an explicit geographical distribution limitation excluding
+those countries, so that distribution is permitted only in or among
+countries not thus excluded. In such case, this License incorporates
+the limitation as if written in the body of this License.
+
+ 9. The Free Software Foundation may publish revised and/or new versions
+of the General Public License from time to time. Such new versions will
+be similar in spirit to the present version, but may differ in detail to
+address new problems or concerns.
+
+Each version is given a distinguishing version number. If the Program
+specifies a version number of this License which applies to it and "any
+later version", you have the option of following the terms and conditions
+either of that version or of any later version published by the Free
+Software Foundation. If the Program does not specify a version number of
+this License, you may choose any version ever published by the Free Software
+Foundation.
+
+ 10. If you wish to incorporate parts of the Program into other free
+programs whose distribution conditions are different, write to the author
+to ask for permission. For software which is copyrighted by the Free
+Software Foundation, write to the Free Software Foundation; we sometimes
+make exceptions for this. Our decision will be guided by the two goals
+of preserving the free status of all derivatives of our free software and
+of promoting the sharing and reuse of software generally.
+
+ NO WARRANTY
+
+ 11. BECAUSE THE PROGRAM IS LICENSED FREE OF CHARGE, THERE IS NO WARRANTY
+FOR THE PROGRAM, TO THE EXTENT PERMITTED BY APPLICABLE LAW. EXCEPT WHEN
+OTHERWISE STATED IN WRITING THE COPYRIGHT HOLDERS AND/OR OTHER PARTIES
+PROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY OF ANY KIND, EITHER EXPRESSED
+OR IMPLIED, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF
+MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE. THE ENTIRE RISK AS
+TO THE QUALITY AND PERFORMANCE OF THE PROGRAM IS WITH YOU. SHOULD THE
+PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF ALL NECESSARY SERVICING,
+REPAIR OR CORRECTION.
+
+ 12. IN NO EVENT UNLESS REQUIRED BY APPLICABLE LAW OR AGREED TO IN WRITING
+WILL ANY COPYRIGHT HOLDER, OR ANY OTHER PARTY WHO MAY MODIFY AND/OR
+REDISTRIBUTE THE PROGRAM AS PERMITTED ABOVE, BE LIABLE TO YOU FOR DAMAGES,
+INCLUDING ANY GENERAL, SPECIAL, INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING
+OUT OF THE USE OR INABILITY TO USE THE PROGRAM (INCLUDING BUT NOT LIMITED
+TO LOSS OF DATA OR DATA BEING RENDERED INACCURATE OR LOSSES SUSTAINED BY
+YOU OR THIRD PARTIES OR A FAILURE OF THE PROGRAM TO OPERATE WITH ANY OTHER
+PROGRAMS), EVEN IF SUCH HOLDER OR OTHER PARTY HAS BEEN ADVISED OF THE
+POSSIBILITY OF SUCH DAMAGES.
+
+ END OF TERMS AND CONDITIONS
+
+ How to Apply These Terms to Your New Programs
+
+ If you develop a new program, and you want it to be of the greatest
+possible use to the public, the best way to achieve this is to make it
+free software which everyone can redistribute and change under these terms.
+
+ To do so, attach the following notices to the program. It is safest
+to attach them to the start of each source file to most effectively
+convey the exclusion of warranty; and each file should have at least
+the "copyright" line and a pointer to where the full notice is found.
+
+
+ Copyright (C)
+
+ This program 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.
+
+ This program 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 program; if not, write to the Free Software
+ Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+
+
+Also add information on how to contact you by electronic and paper mail.
+
+If the program is interactive, make it output a short notice like this
+when it starts in an interactive mode:
+
+ Gnomovision version 69, Copyright (C) year name of author
+ Gnomovision comes with ABSOLUTELY NO WARRANTY; for details type `show w'.
+ This is free software, and you are welcome to redistribute it
+ under certain conditions; type `show c' for details.
+
+The hypothetical commands `show w' and `show c' should show the appropriate
+parts of the General Public License. Of course, the commands you use may
+be called something other than `show w' and `show c'; they could even be
+mouse-clicks or menu items--whatever suits your program.
+
+You should also get your employer (if you work as a programmer) or your
+school, if any, to sign a "copyright disclaimer" for the program, if
+necessary. Here is a sample; alter the names:
+
+ Yoyodyne, Inc., hereby disclaims all copyright interest in the program
+ `Gnomovision' (which makes passes at compilers) written by James Hacker.
+
+ , 1 April 1989
+ Ty Coon, President of Vice
+
+This General Public License does not permit incorporating your program into
+proprietary programs. If your program is a subroutine library, you may
+consider it more useful to permit linking proprietary applications with the
+library. If this is what you want to do, use the GNU Library General
+Public License instead of this License.
Index: contrib/contribs/Centauro/tags/1.0.0/NEWS
===================================================================
--- contrib/contribs/Centauro/tags/1.0.0/NEWS (revision 0)
+++ contrib/contribs/Centauro/tags/1.0.0/NEWS (revision 1278)
@@ -0,0 +1 @@
+08/04/20 releasing version 1.0.0
\ No newline at end of file
Property changes on: contrib/contribs/Centauro/tags/1.0.0/NEWS
___________________________________________________________________
Added: svn:executable
## -0,0 +1 ##
+*
\ No newline at end of property
Index: contrib/contribs/Centauro/tags/1.0.0/Makefile
===================================================================
--- contrib/contribs/Centauro/tags/1.0.0/Makefile (revision 0)
+++ contrib/contribs/Centauro/tags/1.0.0/Makefile (revision 1278)
@@ -0,0 +1,77 @@
+# If you are using this Makefile standalone and fastjet-config is not
+# in your path, edit this line to specify the full path
+FASTJETCONFIG=fastjet-config
+PREFIX=`$(FASTJETCONFIG) --prefix`
+CXX=g++
+CXXFLAGS= -O3 -Wall -g
+install_script = $(SHELL) ../utils/install-sh
+check_script = ../utils/check.sh
+
+# global contrib-wide Makefile include may override some of the above
+# variables (leading "-" means don't give an error if you can't find
+# the file)
+-include ../.Makefile.inc
+
+#------------------------------------------------------------------------
+# things that are specific to this contrib
+NAME=Centauro
+SRCS=Centauro.cc
+EXAMPLES=example
+INSTALLED_HEADERS=Centauro.hh
+#------------------------------------------------------------------------
+
+CXXFLAGS+= $(shell $(FASTJETCONFIG) --cxxflags)
+LDFLAGS += -lm $(shell $(FASTJETCONFIG) --libs)
+
+OBJS = $(SRCS:.cc=.o)
+EXAMPLES_SRCS = $(EXAMPLES:=.cc)
+
+install_HEADER = $(install_script) -c -m 644
+install_LIB = $(install_script) -c -m 644
+install_DIR = $(install_script) -d
+install_DATA = $(install_script) -c -m 644
+install_PROGRAM = $(install_script) -c -s
+install_SCRIPT = $(install_script) -c
+
+.PHONY: clean distclean examples check install
+
+# compilation of the code (default target)
+all: lib$(NAME).a
+
+lib$(NAME).a: $(OBJS)
+ ar cru lib$(NAME).a $(OBJS)
+ ranlib lib$(NAME).a
+
+# building the examples
+examples: $(EXAMPLES)
+
+# the following construct makes it possible to automatically build
+# each of the examples listed in $EXAMPLES
+$(EXAMPLES): % : %.o all
+ $(CXX) -o $@ $< -L. -l$(NAME) $(LDFLAGS)
+
+# check that everything went fine
+check: examples
+ @for prog in $(EXAMPLES); do\
+ $(check_script) $${prog} ../data/single-event.dat || exit 1; \
+ done
+ @echo "All tests successful"
+
+# cleaning the directory
+clean:
+ rm -f *~ *.o
+
+distclean: clean
+ rm -f lib$(NAME).a $(EXAMPLES)
+
+# install things in PREFIX/...
+install: all
+ $(install_DIR) $(PREFIX)/include/fastjet/contrib
+ for header in $(INSTALLED_HEADERS); do\
+ $(install_HEADER) $$header $(PREFIX)/include/fastjet/contrib/;\
+ done
+ $(install_DIR) $(PREFIX)/lib
+ $(install_LIB) lib$(NAME).a $(PREFIX)/lib
+
+depend:
+ makedepend -Y -- -- $(SRCS) $(EXAMPLES_SRCS)
Index: contrib/contribs/Centauro/tags/1.0.0/README
===================================================================
--- contrib/contribs/Centauro/tags/1.0.0/README (revision 0)
+++ contrib/contribs/Centauro/tags/1.0.0/README (revision 1278)
@@ -0,0 +1,3 @@
+Implementation of the Centauro jet algorithm (described in "Asymmetric jet clustering in deep-inelastic scattering" by Miguel Arratia, Yiannis Makris, Duff Neill, Felix Ringer, Nobuo Sato, arXiv:2006.10751) using the FastJet contrib classes.
+
+The Centauro jet algorithm is intended to be used in deep-inelastic scattering (DIS) and is suited for the family of frames connected to the Breit frame by longitudinal boosts. Its distance measure accounts for the forward-backward asymmetry in the Breit frame. It is longitudinally invariant and can cluster jets with Born kinematics, which enables novel studies of TMD observables.
\ No newline at end of file
Index: contrib/contribs/Centauro/tags/1.0.0/Centauro.hh
===================================================================
--- contrib/contribs/Centauro/tags/1.0.0/Centauro.hh (revision 0)
+++ contrib/contribs/Centauro/tags/1.0.0/Centauro.hh (revision 1278)
@@ -0,0 +1,71 @@
+// 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 .
+//----------------------------------------------------------------------
+
+#ifndef __FASTJET_CONTRIB_CENTAUROJETALGORITHM_HH__
+#define __FASTJET_CONTRIB_CENTAUROJETALGORITHM_HH__
+
+#include
+#include "fastjet/JetDefinition.hh"
+#include "fastjet/ClusterSequence.hh"
+FASTJET_BEGIN_NAMESPACE // defined in fastjet/internal/base.hh
+
+namespace contrib{
+
+ class CentauroPlugin : public JetDefinition::Plugin {
+ public:
+
+ /// Constructor for the Centauro Plugin class.
+
+ /// Three floating point arguments are specified to set the parameters
+ /// the radius parameter R, and gammaE and gammaPz are the energy and pz of virtual photon
+
+ CentauroPlugin (double R, double gammaE, double gammaPz) : _R(R), _gammaE(gammaE), _gammaPz(gammaPz){}
+
+ /// if only one argument is passed, assume that it runs in Breit frame so gammaE and gammaPz info not required so set
+ /// so they are set to zero
+
+ CentauroPlugin (double R) : _R(R), _gammaE(0), _gammaPz(0){}
+ /// copy constructor
+ CentauroPlugin (const CentauroPlugin & plugin) {
+ *this = plugin;
+ }
+
+ // the things that are required by base class
+ virtual std::string description () const;
+ virtual void run_clustering(ClusterSequence &) const;
+
+
+ virtual double R() const {return _R;}
+ virtual double gammaE() const {return _gammaE;}
+ virtual double gammaPz() const {return _gammaPz;}
+
+
+
+ /// avoid the warning whenever the user requests "exclusive" jets
+ /// from the cluster sequence
+ virtual bool exclusive_sequence_meaningful() const {return true;}
+
+ private:
+ double _R;
+ double _gammaE;
+ double _gammaPz;
+ };
+
+} // namespace contrib
+
+FASTJET_END_NAMESPACE
+
+#endif // __FASTJET_CONTRIB_CENTAUROJETALGORITHM_HH__
Property changes on: contrib/contribs/Centauro/tags/1.0.0/Centauro.hh
___________________________________________________________________
Added: svn:keywords
## -0,0 +1 ##
+Id
\ No newline at end of property
Index: contrib/contribs/Centauro/tags/1.0.0/example.cc
===================================================================
--- contrib/contribs/Centauro/tags/1.0.0/example.cc (revision 0)
+++ contrib/contribs/Centauro/tags/1.0.0/example.cc (revision 1278)
@@ -0,0 +1,83 @@
+// run it with
+// ./example < single-epDIS-event.dat
+//----------------------------------------------------------------------
+//
+//----------------------------------------------------------------------
+// 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 .
+
+#include
+#include
+
+#include "fastjet/JetDefinition.hh"
+#include "fastjet/PseudoJet.hh"
+#include "fastjet/Selector.hh"
+#include "Centauro.hh"
+#include "fastjet/EECambridgePlugin.hh"
+
+
+using namespace std;
+using namespace fastjet;
+
+void read_event(vector &full_event);
+
+//----------------------------------------------------------------------
+int main(){
+ // example for the creation and execution of the fastjet plugin
+ // for the Centauro jet algorithm
+ //
+ // read in input particles
+ //----------------------------------------------------------
+ vector full_event;
+ read_event(full_event);
+
+ // create the jet definition using the plugin mechanism
+ //----------------------------------------------------------
+ //The Centauro jet algorithm is described in arXiv:2006.10751, "Asymmetric jet clustering in deep-inelastic scattering", Miguel Arratia, Yiannis Makris, Duff Neill, Felix Ringer, Nobuo Sato.
+ fastjet::contrib::CentauroPlugin * centauro_plugin = new fastjet::contrib::CentauroPlugin(1.0);
+ fastjet::JetDefinition jet_def(centauro_plugin);
+ ClusterSequence clust_seq(full_event, jet_def);
+
+ vector jets = clust_seq.inclusive_jets(0);
+
+ vector sortedJets = sorted_by_E(jets);
+ cout << "jets in inclusive clustering " << endl;
+ for (unsigned int i=0; i constituents = jet.constituents();
+ cout << " rap = " << jet.rap() << " e " << jet.e() << " n constituents " << constituents.size() << endl;
+ }
+ cout << endl;
+
+ return 0;
+}
+
+
+//------------------------------------------------------------------------
+// read the event
+void read_event(vector &full_event){
+ string line;
+
+ std::cout << "Particles that are input" << std::endl;
+ 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);
+ std::cout << px << " " << py << " " << pz << " " << E <