diff --git a/analyses/pluginATLAS/ATLAS_2014_I1282441.cc b/analyses/pluginATLAS/ATLAS_2014_I1282441.cc --- a/analyses/pluginATLAS/ATLAS_2014_I1282441.cc +++ b/analyses/pluginATLAS/ATLAS_2014_I1282441.cc @@ -1,94 +1,91 @@ // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Projections/FinalState.hh" #include "Rivet/Projections/UnstableFinalState.hh" #include "Rivet/Projections/IdentifiedFinalState.hh" namespace Rivet { class ATLAS_2014_I1282441 : public Analysis { public: ATLAS_2014_I1282441() : Analysis("ATLAS_2014_I1282441") { } void init() { // Use a large eta range such that we can discriminate on y /// @todo Convert to use a y-cut directly UnstableFinalState ufs(Cuts::abseta < 10 && Cuts::pT > 500*MeV); IdentifiedFinalState phis(ufs); phis.acceptIdPair(PID::PHI); declare(phis, "Phis"); IdentifiedFinalState kpms(Cuts::abseta < 2.0 && Cuts::pT > 230*MeV); kpms.acceptIdPair(PID::KPLUS); declare(kpms, "Kpms"); - _h_phi_rapidity = bookHisto1D(1,1,1); - _h_phi_pT = bookHisto1D(2,1,1); + _h_phi_rapidity = bookHisto1D(2,1,1); + _h_phi_pT = bookHisto1D(1,1,1); } void analyze(const Event& event) { const Particles& ks_all = apply(event, "Kpms").particles(); Particles kp, km; foreach (const Particle& p, ks_all) { if (!p.hasAncestor(PID::PHI)) { MSG_DEBUG("-- K not from phi."); continue; } if (p.p3().mod() > 800*MeV) { MSG_DEBUG("-- p K too high."); continue; } (p.charge() > 0 ? kp : km).push_back(p); } const Particles& phis_all = apply(event, "Phis").particles(); Particles phis; /// @todo Use particles(Cuts&) instead foreach (const Particle& p, phis_all) { if ( p.absrap() > 0.8 ) { MSG_DEBUG("-- phi Y too high."); continue; } if ( p.pT() > 1.2*GeV ) { MSG_DEBUG("-- phi pT too high."); continue; } phis.push_back(p); } // Find Phi -> KK decays through matching of the kinematics if (!kp.empty() && !km.empty() && !phis.empty()) { const double w = event.weight(); MSG_DEBUG("Numbers of particles: #phi=" << phis.size() << ", #K+=" << kp.size() << ", #K-=" << km.size()); for (size_t ip = 0; ip < phis.size(); ++ip) { const Particle& phi = phis[ip]; for (size_t ikm = 0; ikm < km.size(); ++ikm) { for (size_t ikp = 0; ikp < kp.size(); ++ikp) { const FourMomentum mom = kp[ikp].mom() + km[ikm].mom(); if ( fuzzyEquals(mom.mass(), phi.mass(), 1e-5) ) { MSG_DEBUG("Accepted combinatoric: phi#:" << ip << " K+#:" << ikp << " K-#:" << ikm); _h_phi_rapidity->fill(phi.absrap(), w); _h_phi_pT->fill(phi.pT()/MeV, w); } else { MSG_DEBUG("Rejected combinatoric: phi#:" << ip << " K+#:" << ikp << " K-#:" << ikm << " Mass difference is " << mom.mass()-phi.mass()); } } } } } } void finalize() { scale(_h_phi_rapidity, crossSection()/microbarn/sumOfWeights()); scale(_h_phi_pT, crossSection()/microbarn/sumOfWeights()); } private: Histo1DPtr _h_phi_rapidity, _h_phi_pT; }; - - DECLARE_RIVET_PLUGIN(ATLAS_2014_I1282441); - } diff --git a/analyses/pluginATLAS/ATLAS_2014_I1282441.info b/analyses/pluginATLAS/ATLAS_2014_I1282441.info --- a/analyses/pluginATLAS/ATLAS_2014_I1282441.info +++ b/analyses/pluginATLAS/ATLAS_2014_I1282441.info @@ -1,40 +1,41 @@ Name: ATLAS_2014_I1282441 Year: 2014 Summary: The differential production cross section of the $\phi(1020)$ meson in $\sqrt{s}=7$ TeV $pp$ collisions measured with the ATLAS detector Experiment: ATLAS Collider: LHC InspireID: 1282441 Status: VALIDATED Authors: - Tim Martin - Kiran Joshi References: - arXiv:1402.6162 [hep-ex] + - Eur.Phys.J. C74 (2014) no.7, 2895 RunInfo: Run minimum bias events NumEvents: 1000k Beams: [p+, p+] Energies: [7000] PtCuts: [0.5, 1.2] NeedsCrossSection: yes Description: 'A measurement is presented of the $\phi \rightarrow K^+K^-$ production cross section at $\sqrt{s}$ = 7 TeV using pp collision data corresponding to an integrated luminosity of 383 $\mu$b$^{-1}$ collected with the ATLAS experiment at the LHC. Selection of $\phi(1020)$ mesons is based on the identification of charged kaons by their energy loss in the pixel detector. The differential cross section is measured as a function of the transverse momentum, $pT,\phi$, and rapidity, $y_\phi$, of the $\phi(1020)$ meson in the fiducial region 500 < $pT,\phi$ < 1200 MeV, |$y_\phi$| < 0.8, kaon $pT,K$ > 230 MeV and kaon momentum $p_K$ < 800 MeV. The integrated $\phi(1020)$ production cross section in this fiducial range is measured to be $\sigma_{\phi \rightarrow K^+K^-}$ = 570 \pm 8 (stat) \pm 66 (syst) \pm 20 (lumi) $\mu$b.' BibKey: Aad:2014rca BibTeX: '@article{Aad:2014rca, author = "Aad, Georges and others", title = "{The differential production cross section of the $\phi(1020)$ meson in $\sqrt{s}$ = 7 TeV $pp$ collisions measured with the ATLAS detector}", collaboration = "ATLAS Collaboration", year = "2014", eprint = "1402.6162", archivePrefix = "arXiv", primaryClass = "hep-ex", reportNumber = "CERN-PH-EP-2012-269", SLACcitation = "%%CITATION = ARXIV:1402.6162;%%", }' diff --git a/analyses/pluginATLAS/ATLAS_2014_I1282441.plot b/analyses/pluginATLAS/ATLAS_2014_I1282441.plot --- a/analyses/pluginATLAS/ATLAS_2014_I1282441.plot +++ b/analyses/pluginATLAS/ATLAS_2014_I1282441.plot @@ -1,19 +1,19 @@ -# BEGIN PLOT /ATLAS_2014_I1282441/d01-x01-y01 +# BEGIN PLOT /ATLAS_2014_I1282441/d02-x01-y01 Title=Fiducial $\sigma($\phi\rightarrow K^+K^-)$ differential in $|y_\phi|$ XLabel=$|y_\phi|$ YLabel=$\mathrm{d}\sigma/\mathrm{d}y_\phi$~[$\mu$b] LogY=0 RatioPlotMode=datamc RatioPlotYMin=0.8 RatioPlotYMax=1.6 # END PLOT -# BEGIN PLOT /ATLAS_2014_I1282441/d02-x01-y01 +# BEGIN PLOT /ATLAS_2014_I1282441/d01-x01-y01 Title=Fiducial $\sigma($\phi\rightarrow K^+K^-)$ differential in $p_{T,\phi}$ XLabel=p_{T,\phi}$~[MeV] YLabel=$\mathrm{d}\sigma/\mathrm{d}p_{T,\phi}$~[$\mu$b/MeV] LogY=0 RatioPlotMode=datamc RatioPlotYMin=0.8 RatioPlotYMax=1.6 # END PLOT diff --git a/analyses/pluginATLAS/ATLAS_2014_I1282441.yoda b/analyses/pluginATLAS/ATLAS_2014_I1282441.yoda --- a/analyses/pluginATLAS/ATLAS_2014_I1282441.yoda +++ b/analyses/pluginATLAS/ATLAS_2014_I1282441.yoda @@ -1,35 +1,43 @@ BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2014_I1282441/d01-x01-y01 IsRef: 1 Path: /REF/ATLAS_2014_I1282441/d01-x01-y01 -Title: ~ +Title: doi:10.17182/hepdata.64728.v1/t1 Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ -5.000000e-02 5.000000e-02 5.000000e-02 9.000000e+02 1.000412e+02 1.000412e+02 -1.500000e-01 5.000000e-02 5.000000e-02 8.800000e+02 1.056297e+02 1.056297e+02 -2.500000e-01 5.000000e-02 5.000000e-02 8.400000e+02 9.728494e+01 9.728494e+01 -3.500000e-01 5.000000e-02 5.000000e-02 8.200000e+02 9.656443e+01 9.656443e+01 -4.500000e-01 5.000000e-02 5.000000e-02 8.800000e+02 1.164373e+02 1.164373e+02 -5.500000e-01 5.000000e-02 5.000000e-02 6.600000e+02 9.073924e+01 9.073924e+01 -6.500000e-01 5.000000e-02 5.000000e-02 5.100000e+02 7.358412e+01 7.358412e+01 -7.500000e-01 5.000000e-02 5.000000e-02 3.000000e+02 5.432541e+01 5.432541e+01 +5.350000e+02 3.500000e+01 3.500000e+01 4.400000e-01 8.717798e-02 8.717798e-02 +6.050000e+02 3.500000e+01 3.500000e+01 8.700000e-01 1.652271e-01 1.652271e-01 +6.750000e+02 3.500000e+01 3.500000e+01 1.010000e+00 1.854724e-01 1.854724e-01 +7.450000e+02 3.500000e+01 3.500000e+01 1.190000e+00 2.251666e-01 2.251666e-01 +8.150000e+02 3.500000e+01 3.500000e+01 1.180000e+00 2.264950e-01 2.264950e-01 +8.850000e+02 3.500000e+01 3.500000e+01 1.050000e+00 2.027313e-01 2.027313e-01 +9.550000e+02 3.500000e+01 3.500000e+01 7.900000e-01 1.772005e-01 1.772005e-01 +1.025000e+03 3.500000e+01 3.500000e+01 6.700000e-01 1.558846e-01 1.558846e-01 +1.095000e+03 3.500000e+01 3.500000e+01 4.800000e-01 1.122497e-01 1.122497e-01 +1.165000e+03 3.500000e+01 3.500000e+01 4.600000e-01 1.153256e-01 1.153256e-01 END YODA_SCATTER2D_V2 - BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2014_I1282441/d02-x01-y01 IsRef: 1 Path: /REF/ATLAS_2014_I1282441/d02-x01-y01 -Title: ~ +Title: doi:10.17182/hepdata.64728.v1/t2 Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ -5.350000e+02 3.500000e+01 3.500000e+01 4.400000e-01 6.010957e-02 6.010957e-02 -6.050000e+02 3.500000e+01 3.500000e+01 8.700000e-01 1.064481e-01 1.064481e-01 -6.750000e+02 3.500000e+01 3.500000e+01 1.010000e+00 1.125639e-01 1.125639e-01 -7.450000e+02 3.500000e+01 3.500000e+01 1.190000e+00 1.462044e-01 1.462044e-01 -8.150000e+02 3.500000e+01 3.500000e+01 1.180000e+00 1.481442e-01 1.481442e-01 -8.850000e+02 3.500000e+01 3.500000e+01 1.050000e+00 1.323615e-01 1.323615e-01 -9.550000e+02 3.500000e+01 3.500000e+01 7.900000e-01 1.126478e-01 1.126478e-01 -1.025000e+03 3.500000e+01 3.500000e+01 6.700000e-01 1.007914e-01 1.007914e-01 -1.095000e+03 3.500000e+01 3.500000e+01 4.800000e-01 7.705998e-02 7.705998e-02 -1.165000e+03 3.500000e+01 3.500000e+01 4.600000e-01 8.296511e-02 8.296511e-02 +5.000000e-02 5.000000e-02 5.000000e-02 9.000000e-01 2.240536e-01 2.240536e-01 +1.500000e-01 5.000000e-02 5.000000e-02 8.800000e-01 2.177154e-01 2.177154e-01 +2.500000e-01 5.000000e-02 5.000000e-02 8.400000e-01 2.054264e-01 2.054264e-01 +3.500000e-01 5.000000e-02 5.000000e-02 8.200000e-01 2.054264e-01 2.054264e-01 +4.500000e-01 5.000000e-02 5.000000e-02 8.800000e-01 2.231591e-01 2.231591e-01 +5.500000e-01 5.000000e-02 5.000000e-02 6.600000e-01 1.962142e-01 1.962142e-01 +6.500000e-01 5.000000e-02 5.000000e-02 5.100000e-01 1.529706e-01 1.529706e-01 +7.500000e-01 5.000000e-02 5.000000e-02 3.000000e-01 9.380832e-02 9.380832e-02 END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2014_I1282441/d03-x01-y01 +IsRef: 1 +Path: /REF/ATLAS_2014_I1282441/d03-x01-y01 +Title: doi:10.17182/hepdata.64728.v1/t3 +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +7.000000e+03 0.000000e+00 0.000000e+00 5.700000e+02 7.133022e+01 7.133022e+01 +END YODA_SCATTER2D_V2