Index: contrib/contribs/Centauro/tags/1.0/Centauro.cc =================================================================== --- contrib/contribs/Centauro/tags/1.0/Centauro.cc (revision 1274) +++ contrib/contribs/Centauro/tags/1.0/Centauro.cc (revision 1275) @@ -1,134 +0,0 @@ -// 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 Index: contrib/contribs/Centauro/tags/1.0/AUTHORS =================================================================== --- contrib/contribs/Centauro/tags/1.0/AUTHORS (revision 1274) +++ contrib/contribs/Centauro/tags/1.0/AUTHORS (revision 1275) @@ -1,5 +0,0 @@ -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/VERSION =================================================================== --- contrib/contribs/Centauro/tags/1.0/VERSION (revision 1274) +++ contrib/contribs/Centauro/tags/1.0/VERSION (revision 1275) @@ -1 +0,0 @@ -1.0 Index: contrib/contribs/Centauro/tags/1.0/example.ref =================================================================== --- contrib/contribs/Centauro/tags/1.0/example.ref (revision 1274) +++ contrib/contribs/Centauro/tags/1.0/example.ref (revision 1275) @@ -1,471 +0,0 @@ -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/COPYING =================================================================== --- contrib/contribs/Centauro/tags/1.0/COPYING (revision 1274) +++ contrib/contribs/Centauro/tags/1.0/COPYING (revision 1275) @@ -1,343 +0,0 @@ -====================================================================== -====================================================================== -====================================================================== - 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/NEWS =================================================================== --- contrib/contribs/Centauro/tags/1.0/NEWS (revision 1274) +++ contrib/contribs/Centauro/tags/1.0/NEWS (revision 1275) @@ -1 +0,0 @@ -08/03/20. First commit \ No newline at end of file Index: contrib/contribs/Centauro/tags/1.0/Makefile =================================================================== --- contrib/contribs/Centauro/tags/1.0/Makefile (revision 1274) +++ contrib/contribs/Centauro/tags/1.0/Makefile (revision 1275) @@ -1,77 +0,0 @@ -# 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/README =================================================================== --- contrib/contribs/Centauro/tags/1.0/README (revision 1274) +++ contrib/contribs/Centauro/tags/1.0/README (revision 1275) @@ -1,3 +0,0 @@ -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/Centauro.hh =================================================================== --- contrib/contribs/Centauro/tags/1.0/Centauro.hh (revision 1274) +++ contrib/contribs/Centauro/tags/1.0/Centauro.hh (revision 1275) @@ -1,71 +0,0 @@ -// 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__ Index: contrib/contribs/Centauro/tags/1.0/example.cc =================================================================== --- contrib/contribs/Centauro/tags/1.0/example.cc (revision 1274) +++ contrib/contribs/Centauro/tags/1.0/example.cc (revision 1275) @@ -1,83 +0,0 @@ -// 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 <