diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z.cc b/analyses/pluginATLAS/ATLAS_2013_I1217863.cc rename from analyses/pluginATLAS/ATLAS_2013_I1217863_Z.cc rename to analyses/pluginATLAS/ATLAS_2013_I1217863.cc --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z.cc +++ b/analyses/pluginATLAS/ATLAS_2013_I1217863.cc @@ -1,195 +1,290 @@ // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Projections/FinalState.hh" #include "Rivet/Projections/FastJets.hh" #include "Rivet/Projections/ZFinder.hh" #include "Rivet/Projections/VetoedFinalState.hh" #include "Rivet/Projections/LeadingParticlesFinalState.hh" namespace Rivet { - class ATLAS_2013_I1217863_Z : public Analysis { + class ATLAS_2013_I1217863 : public Analysis { public: + /// Constructor - ATLAS_2013_I1217863_Z(string name="ATLAS_2013_I1217863_Z") - : Analysis(name) - { - // the electron mode is used by default - _mode = 1; - } - + ///@brief: Electroweak Wjj production at 8 TeV + DEFAULT_RIVET_ANALYSIS_CTOR(ATLAS_2013_I1217863); + //@} /// @name Analysis methods //@{ /// Book histograms and initialise projections before the run void init() { + // Get options from the new option system + _mode = 2; + _doZ = true; + _doW = true; + if ( getOption("LMODE") == "EL" ) { _mode = 2;} + if ( getOption("LMODE") == "MU" ) _mode = 3; + if ( getOption("LMODE") == "ZEL" ) { + _mode = 2; + _doW = false; + } + if ( getOption("LMODE") == "ZMU" ) { + _mode = 3; + _doW = false; + } + if ( getOption("LMODE") == "WEL" ) { + _mode = 2; + _doZ = false; + } + if ( getOption("LMODE") == "WMU" ) { + _mode = 3; + _doZ = false; + } + FinalState fs; declare(fs, "FS"); Cut cuts = Cuts::abseta < 2.47 && Cuts::pT > 25*GeV; // Z finder - ZFinder zf(fs, cuts, _mode==3? PID::MUON : PID::ELECTRON, 40.0*GeV, 1000.0*GeV, 0.1, ZFinder::CLUSTERNODECAY, ZFinder::NOTRACK); - declare(zf, "ZF"); + if (_doZ) { + ZFinder zf(fs, cuts, _mode==3? PID::MUON : PID::ELECTRON, 40.0*GeV, 1000.0*GeV, 0.1, ZFinder::CLUSTERNODECAY, ZFinder::NOTRACK); + declare(zf, "ZF"); + } + + if (_doW) { + // W finder for electrons and muons + WFinder wf(fs, cuts, _mode==3? PID::MUON : PID::ELECTRON, 0.0*GeV, 1000.0*GeV, 35.0*GeV, 0.1, + WFinder::CLUSTERNODECAY, WFinder::NOTRACK, WFinder::TRANSMASS); + declare(wf, "WF"); + } // leading photon LeadingParticlesFinalState photonfs(FinalState(Cuts::abseta < 2.37 && Cuts::pT > 15*GeV)); photonfs.addParticleId(PID::PHOTON); declare(photonfs, "LeadingPhoton"); // jets VetoedFinalState jet_fs(fs); - jet_fs.addVetoOnThisFinalState(getProjection("ZF")); + if (_doZ) { jet_fs.addVetoOnThisFinalState(getProjection("ZF")); } + if (_doW) { jet_fs.addVetoOnThisFinalState(getProjection("WF")); } jet_fs.addVetoOnThisFinalState(getProjection("LeadingPhoton")); FastJets jets(jet_fs, FastJets::ANTIKT, 0.4); jets.useInvisibles(true); declare(jets, "Jets"); // FS excluding the leading photon VetoedFinalState vfs(fs); vfs.addVetoOnThisFinalState(photonfs); declare(vfs, "isolatedFS"); // Book histograms - _hist_EgammaT_incl = bookHisto1D(11, 1, _mode); // dSigma / dE^gamma_T for Njet >= 0 - _hist_EgammaT_excl = bookHisto1D(12, 1, _mode); // dSigma / dE^gamma_T for Njet = 0 - _hist_Njet_EgammaT15 = bookHisto1D(17, 1, _mode); // dSigma / dNjet for E^gamma_T >= 15 - _hist_Njet_EgammaT60 = bookHisto1D(18, 1, _mode); // dSigma / dNjet for E^gamma_T >= 60 - _hist_mZgamma = bookHisto1D(20, 1, _mode); // dSigma / dm^{Zgamma} + if (_doZ) { + _hist_EgammaT_inclZ = bookHisto1D(11, 1, _mode); // dSigma / dE^gamma_T for Njet >= 0 + _hist_EgammaT_exclZ = bookHisto1D(12, 1, _mode); // dSigma / dE^gamma_T for Njet = 0 + _hist_Njet_EgammaT15Z = bookHisto1D(17, 1, _mode); // dSigma / dNjet for E^gamma_T >= 15 + _hist_Njet_EgammaT60Z = bookHisto1D(18, 1, _mode); // dSigma / dNjet for E^gamma_T >= 60 + _hist_mZgamma = bookHisto1D(20, 1, _mode); // dSigma / dm^{Zgamma} + } + + if (_doW){ + // Book histograms + _hist_EgammaT_inclW = bookHisto1D( 7, 1, _mode); // dSigma / dE^gamma_T for Njet >= 0 + _hist_EgammaT_exclW = bookHisto1D( 8, 1, _mode); // dSigma / dE^gamma_T for Njet = 0 + _hist_Njet_EgammaT15W = bookHisto1D(15, 1, _mode); // dSigma / dNjet for E^gamma_T > 15 + _hist_Njet_EgammaT60W = bookHisto1D(16, 1, _mode); // dSigma / dNjet for E^gamma_T > 60 + _hist_mWgammaT = bookHisto1D(19, 1, _mode); // dSigma / dm^{Wgamma}_T + } } /// Perform the per-event analysis void analyze(const Event& event) { const double weight = event.weight(); // retrieve leading photon Particles photons = apply(event, "LeadingPhoton").particles(); if (photons.size() != 1) vetoEvent; const Particle& leadingPhoton = photons[0]; if (leadingPhoton.Et() < 15.0*GeV) vetoEvent; if (leadingPhoton.abseta() > 2.37) vetoEvent; // check photon isolation double coneEnergy(0.0); Particles fs = apply(event, "isolatedFS").particles(); for (const Particle& p : fs) { if ( deltaR(leadingPhoton, p) < 0.4 ) coneEnergy += p.E(); } if (coneEnergy / leadingPhoton.E() >= 0.5 ) vetoEvent; - // retrieve Z boson candidate - const ZFinder& zf = apply(event, "ZF"); - if ( zf.bosons().size() != 1 ) vetoEvent; // only one Z boson candidate - const Particle& Zboson = zf.boson(); - if ( !(Zboson.mass() > 40.0*GeV) ) vetoEvent; - - // check charge of constituent leptons - const ParticleVector& leptons = zf.constituents(); - if (leptons.size() != 2 || leptons[0].charge() * leptons[1].charge() > 0.) vetoEvent; - - // check photon-lepton overlap - for (const Particle& p : leptons) { - if ( !(p.pT() > 25.0*GeV && p.abseta() < 2.47 && deltaR(leadingPhoton, p) > 0.7) ) vetoEvent; + if (_doW) { + // retrieve W boson candidate + const WFinder& wf = apply(event, "WF"); + if ( wf.bosons().size() == 1 ) { + + // retrieve constituent neutrino + const Particle& neutrino = wf.constituentNeutrino(); + if ( (neutrino.pT() > 35.0*GeV) ) { + + // retrieve constituent lepton + const Particle& lepton = wf.constituentLepton(); + if ( lepton.pT() > 25.0*GeV && lepton.abseta() < 2.47 ) { + + // check photon-lepton overlap + if ( deltaR(leadingPhoton, lepton) > 0.7 ) { + + // count jets + const FastJets& jetfs = apply(event, "Jets"); + Jets jets = jetfs.jets(cmpMomByEt); + int goodJets = 0; + for (const Jet& j : jets) { + if ( !(j.Et() > 30.0*GeV) ) break; + if ( (j.abseta() < 4.4) && \ + (deltaR(leadingPhoton, j) > 0.3) && \ + (deltaR(lepton, j) > 0.3) ) ++goodJets; + } + + double Njets = double(goodJets) + 0.5; + double photonEt = leadingPhoton.Et()*GeV; + + const FourMomentum& lep_gamma = lepton.momentum() + leadingPhoton.momentum(); + double term1 = sqrt(lep_gamma.mass2() + lep_gamma.pT2()) + neutrino.Et(); + double term2 = (lep_gamma + neutrino.momentum()).pT2(); + double mWgammaT = sqrt(term1 * term1 - term2) * GeV; + + _hist_EgammaT_inclW->fill(photonEt, weight); + + _hist_Njet_EgammaT15W->fill(Njets, weight); + + if ( !goodJets ) _hist_EgammaT_exclW->fill(photonEt, weight); + + if (photonEt > 40.0*GeV) { + _hist_mWgammaT->fill(mWgammaT, weight); + if (photonEt > 60.0*GeV) _hist_Njet_EgammaT60W->fill(Njets, weight); + } + } + } + } + } } - // count jets - const FastJets& jetfs = apply(event, "Jets"); - Jets jets = jetfs.jets(cmpMomByEt); - int goodJets = 0; - foreach (const Jet& j, jets) { - if ( !(j.Et() > 30.0*GeV) ) break; - if ( (j.abseta() < 4.4) && \ - (deltaR(leadingPhoton, j) > 0.3) && \ - (deltaR(leptons[0], j) > 0.3) && \ - (deltaR(leptons[1], j) > 0.3) ) ++goodJets; + if (_doZ ){ + + // retrieve Z boson candidate + const ZFinder& zf = apply(event, "ZF"); + if ( zf.bosons().size() == 1 ) { + const Particle& Zboson = zf.boson(); + if ( (Zboson.mass() > 40.0*GeV) ) { + + // check charge of constituent leptons + const ParticleVector& leptons = zf.constituents(); + if (leptons.size() == 2 && leptons[0].charge() * leptons[1].charge() < 0.) { + + bool lpass = true; + // check photon-lepton overlap + for (const Particle& p : leptons) { + if ( !(p.pT() > 25.0*GeV && p.abseta() < 2.47 && deltaR(leadingPhoton, p) > 0.7) ) lpass = false; + } + if ( lpass ) { + + // count jets + const FastJets& jetfs = apply(event, "Jets"); + Jets jets = jetfs.jets(cmpMomByEt); + int goodJets = 0; + foreach (const Jet& j, jets) { + if ( !(j.Et() > 30.0*GeV) ) break; + if ( (j.abseta() < 4.4) && \ + (deltaR(leadingPhoton, j) > 0.3) && \ + (deltaR(leptons[0], j) > 0.3) && \ + (deltaR(leptons[1], j) > 0.3) ) ++goodJets; + } + + double Njets = double(goodJets) + 0.5; + double photonEt = leadingPhoton.Et()*GeV; + double mZgamma = (Zboson.momentum() + leadingPhoton.momentum()).mass() * GeV; + + _hist_EgammaT_inclZ->fill(photonEt, weight); + + _hist_Njet_EgammaT15Z->fill(Njets, weight); + + if ( !goodJets ) _hist_EgammaT_exclZ->fill(photonEt, weight); + + if (photonEt >= 40.0*GeV) { + _hist_mZgamma->fill(mZgamma, weight); + if (photonEt >= 60.0*GeV) _hist_Njet_EgammaT60Z->fill(Njets, weight); + } + } + } + } + } } - - double Njets = double(goodJets) + 0.5; - double photonEt = leadingPhoton.Et()*GeV; - double mZgamma = (Zboson.momentum() + leadingPhoton.momentum()).mass() * GeV; - - _hist_EgammaT_incl->fill(photonEt, weight); - - _hist_Njet_EgammaT15->fill(Njets, weight); - - if ( !goodJets ) _hist_EgammaT_excl->fill(photonEt, weight); - - if (photonEt >= 40.0*GeV) { - _hist_mZgamma->fill(mZgamma, weight); - if (photonEt >= 60.0*GeV) _hist_Njet_EgammaT60->fill(Njets, weight); - } - } - /// Normalise histograms etc., after the run void finalize() { const double xs_fb = crossSection()/femtobarn; const double sumw = sumOfWeights(); const double sf = xs_fb / sumw; - scale(_hist_EgammaT_excl, sf); - scale(_hist_EgammaT_incl, sf); + if (_doZ) { + scale(_hist_EgammaT_exclZ, sf); + scale(_hist_EgammaT_inclZ, sf); + normalize(_hist_Njet_EgammaT15Z); + normalize(_hist_Njet_EgammaT60Z); + normalize(_hist_mZgamma); + } - normalize(_hist_Njet_EgammaT15); - normalize(_hist_Njet_EgammaT60); - normalize(_hist_mZgamma); + if (_doW) { + scale(_hist_EgammaT_exclW, sf); + scale(_hist_EgammaT_inclW, sf); + normalize(_hist_Njet_EgammaT15W); + normalize(_hist_Njet_EgammaT60W); + normalize(_hist_mWgammaT); + } } //@} protected: size_t _mode; + bool _doW; + bool _doZ; private: /// @name Histograms //@{ - Histo1DPtr _hist_EgammaT_incl; - Histo1DPtr _hist_EgammaT_excl; - Histo1DPtr _hist_Njet_EgammaT15; - Histo1DPtr _hist_Njet_EgammaT60; + Histo1DPtr _hist_EgammaT_inclZ; + Histo1DPtr _hist_EgammaT_exclZ; + Histo1DPtr _hist_Njet_EgammaT15Z; + Histo1DPtr _hist_Njet_EgammaT60Z; Histo1DPtr _hist_mZgamma; + Histo1DPtr _hist_EgammaT_inclW; + Histo1DPtr _hist_EgammaT_exclW; + Histo1DPtr _hist_Njet_EgammaT15W; + Histo1DPtr _hist_Njet_EgammaT60W; + Histo1DPtr _hist_mWgammaT; + //@} }; - - class ATLAS_2013_I1217863_Z_EL : public ATLAS_2013_I1217863_Z { - public: - ATLAS_2013_I1217863_Z_EL() - : ATLAS_2013_I1217863_Z("ATLAS_2013_I1217863_Z_EL") - { - _mode = 2; - } - }; - - - class ATLAS_2013_I1217863_Z_MU : public ATLAS_2013_I1217863_Z { - public: - ATLAS_2013_I1217863_Z_MU() - : ATLAS_2013_I1217863_Z("ATLAS_2013_I1217863_Z_MU") - { - _mode = 3; - } - }; - - - DECLARE_RIVET_PLUGIN(ATLAS_2013_I1217863_Z); - DECLARE_RIVET_PLUGIN(ATLAS_2013_I1217863_Z_EL); - DECLARE_RIVET_PLUGIN(ATLAS_2013_I1217863_Z_MU); + DECLARE_RIVET_PLUGIN(ATLAS_2013_I1217863); } + + diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z.info b/analyses/pluginATLAS/ATLAS_2013_I1217863.info rename from analyses/pluginATLAS/ATLAS_2013_I1217863_Z.info rename to analyses/pluginATLAS/ATLAS_2013_I1217863.info --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z.info +++ b/analyses/pluginATLAS/ATLAS_2013_I1217863.info @@ -1,45 +1,47 @@ -Name: ATLAS_2013_I1217863_Z +Name: ATLAS_2013_I1217863 Year: 2013 Summary: Z + gamma production at 7 TeV Experiment: ATLAS Collider: LHC InspireID: 1217863 Status: VALIDATED Authors: - Chritian Gutschow References: - Phys.Rev. D87 (2013) 112003 - doi:10.1103/PhysRevD.87.112003 - arXiv:1302.1283 [hep-ex] RunInfo: Z+gamma in the electron channel NumEvents: 1000000 Beams: [p+, p+] Energies: [7000] PtCuts: [25,15] NeedCrossSection: True +Options: + - LMODE=EL,MU,ZEL,ZMU,WEL,WMU Description: Measurements of the differential fiducial cross sections for the production of a Z boson in association with a high-energy photon are measured using pp collisions at $\sqrt{s}=7~\TeV$. The analysis uses a data sample with an integrated luminosity of 4.6/fb collected by the ATLAS detector during the 2011 LHC data-taking period. Events are selected using leptonic decays of the Z bosons with the requirement of an associated isolated photon. The default routine will consider the electron decay - channel of the Z boson. Use ATLAS_2013_I1217863_Z_EL and ATLAS_2013_I1217863_Z_MU to specify the decay channel directly. + channel of the Z boson. Use LMODE to specify the decay channel directly. BibKey: Aad:2013iz BibTeX: '@article{Aad:2013izg, author = "Aad, Georges and others", title = "{Measurements of $W\gamma$ and $Z\gamma$ production in $pp$ collisions at $\sqrt{s}$=7  TeV with the ATLAS detector at the LHC}", collaboration = "ATLAS Collaboration", journal = "Phys.Rev.", number = "11", volume = "D87", pages = "112003", doi = "10.1103/PhysRevD.87.112003", year = "2013", eprint = "1302.1283", archivePrefix = "arXiv", primaryClass = "hep-ex", reportNumber = "CERN-PH-EP-2012-345", SLACcitation = "%%CITATION = ARXIV:1302.1283;%%", }' diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863.plot b/analyses/pluginATLAS/ATLAS_2013_I1217863.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginATLAS/ATLAS_2013_I1217863.plot @@ -0,0 +1,139 @@ +# BEGIN PLOT /ATLAS_2013_I1217863/d..-x..-y.. +LogY=1 +LegendYPos=0.90 +LegendXPos=0.90 +LegendAlign=r +RatioPlotYMax=1.5 +RatioPlotYMin=0.5 +XTwosidedTicks=1 +YTwosidedTicks=1 +ErrorBars=1 +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d1[1-2]-x01-y.. +LogX=1 +RatioPlotYMax=2.0 +RatioPlotYMin=0.0 +XLabel=$E^\gamma_\text{T}$ [GeV] +YLabel=$\text{d}\sigma(pp \rightarrow \ell^+ \ell^- \gamma) / \text{d}E^\gamma_\text{T}$ [fb $\text{GeV}^{-1}$] +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d1[7-8]-x01-y.. +XLabel=Jet multiplicity +YLabel=$1/\sigma_{Z\gamma} \times \text{d}\sigma_{Z\gamma} / \text{d} N_\text{jet}$ +# END PLOT + +# Z electron + +# BEGIN PLOT /ATLAS_2013_I1217863/d20-x01-y.. +XLabel=$m^{Z\gamma}$ [GeV] +YLabel=$1/\sigma_{Z\gamma} \times \text{d}\sigma_{Z\gamma} / \text{d} m^{Z\gamma}$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d11-x01-y01 +Title=$N_\text{jet} \geq 0$, using electrons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d12-x01-y01 +Title=$N_\text{jet} = 0$, using electrons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d17-x01-y01 +Title=$E^\gamma_\text{T} \geq 15$ GeV, using electrons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d18-x01-y01 +Title=$E^\gamma_\text{T} \geq 60$ GeV, using electrons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d20-x01-y01 +Title=$E^\gamma_\text{T} > 40$ GeV, using electrons, dressed level +# END PLOT + + +# W electron + +# BEGIN PLOT /ATLAS_2013_I1217863/d0[7-8]-x01-y.. +LogX=1 +RatioPlotYMax=2.0 +RatioPlotYMin=0.0 +XLabel=$E^\gamma_\text{T}$ [GeV] +YLabel=$\text{d}\sigma(pp \rightarrow \ell \nu \gamma) / \text{d}E^\gamma_\text{T}$ [fb $\text{GeV}^{-1}$] +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d1[5-6]-x01-y.. +XLabel=Jet multiplicity +YLabel=$1/\sigma_{W\gamma} \times \text{d}\sigma_{W\gamma} / \text{d} N_\text{jet}$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d07-x01-y02 +Title=$N_\text{jet} \geq 0$, using electrons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d08-x01-y01 +Title=$N_\text{jet} = 0$, using electrons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d15-x01-y01 +Title=$E^\gamma_\text{T} \geq 15$ GeV, using electrons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d16-x01-y01 +Title=$E^\gamma_\text{T} \geq 60$ GeV, using electrons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d19-x01-y01 +Title=$E^\gamma_\text{T} > 40$ GeV, using electrons, dressed level +# END PLOT + + +# Z, muon + +# BEGIN PLOT /ATLAS_2013_I1217863/d11-x01-y01 +Title=$N_\text{jet} \geq 0$, using muons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d12-x01-y01 +Title=$N_\text{jet} = 0$, using muons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d17-x01-y01 +Title=$E^\gamma_\text{T} \geq 15$ GeV, using muons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d18-x01-y01 +Title=$E^\gamma_\text{T} \geq 60$ GeV, using muons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d20-x01-y01 +Title=$E^\gamma_\text{T} > 40$ GeV, using muons, dressed level +# END PLOT + + +# W, muon + +# BEGIN PLOT /ATLAS_2013_I1217863/d19-x01-y.. +XLabel=$m^{W\gamma}_\text{T}$ [GeV] +YLabel=$1/\sigma_{W\gamma} \times \text{d}\sigma_{W\gamma} / \text{d} m^{W\gamma}_\text{T}$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d07-x01-y01 +Title=$N_\text{jet} \geq 0$, using muons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d08-x01-y01 +Title=$N_\text{jet} = 0$, using muons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d15-x01-y01 +Title=$E^\gamma_\text{T} \geq 15$ GeV, using muons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d16-x01-y01 +Title=$E^\gamma_\text{T} \geq 60$ GeV, using muons, dressed level +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1217863/d19-x01-y01 +Title=$E^\gamma_\text{T} > 40$ GeV, using muons, dressed level +# END PLOT + diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863.yoda b/analyses/pluginATLAS/ATLAS_2013_I1217863.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginATLAS/ATLAS_2013_I1217863.yoda @@ -0,0 +1,278 @@ +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d07-x01-y02 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d07-x01-y02 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.750000e+01 2.500000e+00 2.500000e+00 1.947072e+02 3.245120e+01 3.245120e+01 +2.500000e+01 5.000000e+00 5.000000e+00 8.518440e+01 1.115510e+01 1.115510e+01 +3.500000e+01 5.000000e+00 5.000000e+00 4.360630e+01 5.070500e+00 5.070500e+00 +5.000000e+01 1.000000e+01 1.000000e+01 1.409599e+01 1.825380e+00 1.825380e+00 +8.000000e+01 2.000000e+01 2.000000e+01 5.070500e+00 5.070500e-01 5.070500e-01 +5.500000e+02 4.500000e+02 4.500000e+02 9.126900e-02 1.216920e-02 1.216920e-02 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d08-x01-y02 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d08-x01-y02 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.750000e+01 2.500000e+00 2.500000e+00 1.379176e+02 2.231020e+01 2.231020e+01 +2.500000e+01 5.000000e+00 5.000000e+00 5.526845e+01 7.200110e+00 7.200110e+00 +3.500000e+01 5.000000e+00 5.000000e+00 2.393276e+01 3.346530e+00 3.346530e+00 +5.000000e+01 1.000000e+01 1.000000e+01 6.997290e+00 1.216920e+00 1.216920e+00 +8.000000e+01 2.000000e+01 2.000000e+01 2.129610e+00 3.042300e-01 3.042300e-01 +5.500000e+02 4.500000e+02 4.500000e+02 3.042300e-02 6.084600e-03 6.084600e-03 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d15-x01-y02 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d15-x01-y02 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +5.000000e-01 5.000000e-01 5.000000e-01 5.800000e-01 3.000000e-02 3.000000e-02 +1.500000e+00 5.000000e-01 5.000000e-01 3.300000e-01 3.000000e-02 3.000000e-02 +2.500000e+00 5.000000e-01 5.000000e-01 1.200000e-01 1.000000e-02 1.000000e-02 +3.500000e+00 5.000000e-01 5.000000e-01 3.700000e-02 8.000000e-03 8.000000e-03 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d16-x01-y02 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d16-x01-y02 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +5.000000e-01 5.000000e-01 5.000000e-01 3.300000e-01 3.300000e-02 3.300000e-02 +1.500000e+00 5.000000e-01 5.000000e-01 4.000000e-01 2.600000e-02 2.600000e-02 +2.500000e+00 5.000000e-01 5.000000e-01 1.600000e-01 1.800000e-02 1.800000e-02 +3.500000e+00 5.000000e-01 5.000000e-01 1.100000e-01 2.400000e-02 2.400000e-02 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d19-x01-y02 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d19-x01-y02 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +7.500000e+01 2.500000e+01 2.500000e+01 7.400000e-04 1.980000e-04 1.980000e-04 +1.250000e+02 2.500000e+01 2.500000e+01 5.400000e-03 4.400000e-04 4.400000e-04 +1.750000e+02 2.500000e+01 2.500000e+01 6.400000e-03 4.800000e-04 4.800000e-04 +2.500000e+02 5.000000e+01 5.000000e+01 2.400000e-03 1.500000e-04 1.500000e-04 +3.500000e+02 5.000000e+01 5.000000e+01 7.900000e-04 8.700000e-05 8.700000e-05 +7.000000e+02 3.000000e+02 3.000000e+02 8.333333e-05 1.416667e-05 1.416667e-05 +END YODA_SCATTER2D_V2 + + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d07-x01-y03 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d07-x01-y03 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.750000e+01 2.500000e+00 2.500000e+00 1.947072e+02 3.245120e+01 3.245120e+01 +2.500000e+01 5.000000e+00 5.000000e+00 8.518440e+01 1.115510e+01 1.115510e+01 +3.500000e+01 5.000000e+00 5.000000e+00 4.360630e+01 5.070500e+00 5.070500e+00 +5.000000e+01 1.000000e+01 1.000000e+01 1.409599e+01 1.825380e+00 1.825380e+00 +8.000000e+01 2.000000e+01 2.000000e+01 5.070500e+00 5.070500e-01 5.070500e-01 +5.500000e+02 4.500000e+02 4.500000e+02 9.126900e-02 1.216920e-02 1.216920e-02 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d08-x01-y03 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d08-x01-y03 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.750000e+01 2.500000e+00 2.500000e+00 1.379176e+02 2.231020e+01 2.231020e+01 +2.500000e+01 5.000000e+00 5.000000e+00 5.526845e+01 7.200110e+00 7.200110e+00 +3.500000e+01 5.000000e+00 5.000000e+00 2.393276e+01 3.346530e+00 3.346530e+00 +5.000000e+01 1.000000e+01 1.000000e+01 6.997290e+00 1.216920e+00 1.216920e+00 +8.000000e+01 2.000000e+01 2.000000e+01 2.129610e+00 3.042300e-01 3.042300e-01 +5.500000e+02 4.500000e+02 4.500000e+02 3.042300e-02 6.084600e-03 6.084600e-03 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d15-x01-y03 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d15-x01-y03 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +5.000000e-01 5.000000e-01 5.000000e-01 5.800000e-01 3.000000e-02 3.000000e-02 +1.500000e+00 5.000000e-01 5.000000e-01 3.300000e-01 3.000000e-02 3.000000e-02 +2.500000e+00 5.000000e-01 5.000000e-01 1.200000e-01 1.000000e-02 1.000000e-02 +3.500000e+00 5.000000e-01 5.000000e-01 3.700000e-02 8.000000e-03 8.000000e-03 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d16-x01-y03 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d16-x01-y03 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +5.000000e-01 5.000000e-01 5.000000e-01 3.300000e-01 3.300000e-02 3.300000e-02 +1.500000e+00 5.000000e-01 5.000000e-01 4.000000e-01 2.600000e-02 2.600000e-02 +2.500000e+00 5.000000e-01 5.000000e-01 1.600000e-01 1.800000e-02 1.800000e-02 +3.500000e+00 5.000000e-01 5.000000e-01 1.100000e-01 2.400000e-02 2.400000e-02 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d19-x01-y03 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d19-x01-y03 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +7.500000e+01 2.500000e+01 2.500000e+01 7.400000e-04 1.980000e-04 1.980000e-04 +1.250000e+02 2.500000e+01 2.500000e+01 5.400000e-03 4.400000e-04 4.400000e-04 +1.750000e+02 2.500000e+01 2.500000e+01 6.400000e-03 4.800000e-04 4.800000e-04 +2.500000e+02 5.000000e+01 5.000000e+01 2.400000e-03 1.500000e-04 1.500000e-04 +3.500000e+02 5.000000e+01 5.000000e+01 7.900000e-04 8.700000e-05 8.700000e-05 +7.000000e+02 3.000000e+02 3.000000e+02 8.333333e-05 1.416667e-05 1.416667e-05 +END YODA_SCATTER2D_V2 + + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d11-x01-y02 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d11-x01-y02 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.750000e+01 2.500000e+00 2.500000e+00 1.216920e+02 1.216920e+01 1.216920e+01 +2.500000e+01 5.000000e+00 5.000000e+00 4.309925e+01 4.259220e+00 4.259220e+00 +3.500000e+01 5.000000e+00 5.000000e+00 1.318330e+01 1.419740e+00 1.419740e+00 +5.000000e+01 1.000000e+01 1.000000e+01 5.009654e+00 6.186010e-01 6.186010e-01 +8.000000e+01 2.000000e+01 2.000000e+01 1.419740e+00 1.926790e-01 1.926790e-01 +5.500000e+02 4.500000e+02 4.500000e+02 1.825380e-02 7.098700e-03 7.098700e-03 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d12-x01-y02 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d12-x01-y02 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.750000e+01 2.500000e+00 2.500000e+00 1.074946e+02 1.115510e+01 1.115510e+01 +2.500000e+01 5.000000e+00 5.000000e+00 3.478363e+01 4.157810e+00 4.157810e+00 +3.500000e+01 5.000000e+00 5.000000e+00 9.431130e+00 1.115510e+00 1.115510e+00 +5.000000e+01 1.000000e+01 1.000000e+01 3.285684e+00 4.766270e-01 4.766270e-01 +8.000000e+01 2.000000e+01 2.000000e+01 1.024241e+00 1.622560e-01 1.622560e-01 +5.500000e+02 4.500000e+02 4.500000e+02 7.098700e-03 4.056400e-03 4.056400e-03 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d17-x01-y02 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d17-x01-y02 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +5.000000e-01 5.000000e-01 5.000000e-01 8.000000e-01 2.700000e-02 2.700000e-02 +1.500000e+00 5.000000e-01 5.000000e-01 1.500000e-01 1.600000e-02 1.600000e-02 +2.500000e+00 5.000000e-01 5.000000e-01 5.200000e-02 1.100000e-02 1.100000e-02 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d18-x01-y02 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d18-x01-y02 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +5.000000e-01 5.000000e-01 5.000000e-01 6.000000e-01 3.800000e-02 3.800000e-02 +1.500000e+00 5.000000e-01 5.000000e-01 2.800000e-01 2.600000e-02 2.600000e-02 +2.500000e+00 5.000000e-01 5.000000e-01 1.200000e-01 1.900000e-02 1.900000e-02 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d20-x01-y02 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d20-x01-y02 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.000000e+02 5.000000e+01 5.000000e+01 1.400000e-03 1.800000e-04 1.800000e-04 +1.750000e+02 2.500000e+01 2.500000e+01 8.800000e-03 4.800000e-04 4.800000e-04 +2.500000e+02 5.000000e+01 5.000000e+01 3.200000e-03 2.200000e-04 2.200000e-04 +6.500000e+02 3.500000e+02 3.500000e+02 1.571429e-04 2.142857e-05 2.142857e-05 +END YODA_SCATTER2D_V2 + + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d11-x01-y03 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d11-x01-y03 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.750000e+01 2.500000e+00 2.500000e+00 1.216920e+02 1.216920e+01 1.216920e+01 +2.500000e+01 5.000000e+00 5.000000e+00 4.309925e+01 4.259220e+00 4.259220e+00 +3.500000e+01 5.000000e+00 5.000000e+00 1.318330e+01 1.419740e+00 1.419740e+00 +5.000000e+01 1.000000e+01 1.000000e+01 5.009654e+00 6.186010e-01 6.186010e-01 +8.000000e+01 2.000000e+01 2.000000e+01 1.419740e+00 1.926790e-01 1.926790e-01 +5.500000e+02 4.500000e+02 4.500000e+02 1.825380e-02 7.098700e-03 7.098700e-03 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d12-x01-y03 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d12-x01-y03 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.750000e+01 2.500000e+00 2.500000e+00 1.074946e+02 1.115510e+01 1.115510e+01 +2.500000e+01 5.000000e+00 5.000000e+00 3.478363e+01 4.157810e+00 4.157810e+00 +3.500000e+01 5.000000e+00 5.000000e+00 9.431130e+00 1.115510e+00 1.115510e+00 +5.000000e+01 1.000000e+01 1.000000e+01 3.285684e+00 4.766270e-01 4.766270e-01 +8.000000e+01 2.000000e+01 2.000000e+01 1.024241e+00 1.622560e-01 1.622560e-01 +5.500000e+02 4.500000e+02 4.500000e+02 7.098700e-03 4.056400e-03 4.056400e-03 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d17-x01-y03 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d17-x01-y03 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +5.000000e-01 5.000000e-01 5.000000e-01 8.000000e-01 2.700000e-02 2.700000e-02 +1.500000e+00 5.000000e-01 5.000000e-01 1.500000e-01 1.600000e-02 1.600000e-02 +2.500000e+00 5.000000e-01 5.000000e-01 5.200000e-02 1.100000e-02 1.100000e-02 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d18-x01-y03 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d18-x01-y03 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +5.000000e-01 5.000000e-01 5.000000e-01 6.000000e-01 3.800000e-02 3.800000e-02 +1.500000e+00 5.000000e-01 5.000000e-01 2.800000e-01 2.600000e-02 2.600000e-02 +2.500000e+00 5.000000e-01 5.000000e-01 1.200000e-01 1.900000e-02 1.900000e-02 +END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863/d20-x01-y03 +IsRef: 1 +Path: /REF/ATLAS_2013_I1217863/d20-x01-y03 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.000000e+02 5.000000e+01 5.000000e+01 1.400000e-03 1.800000e-04 1.800000e-04 +1.750000e+02 2.500000e+01 2.500000e+01 8.800000e-03 4.800000e-04 4.800000e-04 +2.500000e+02 5.000000e+01 5.000000e+01 3.200000e-03 2.200000e-04 2.200000e-04 +6.500000e+02 3.500000e+02 3.500000e+02 1.571429e-04 2.142857e-05 2.142857e-05 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_W.cc b/analyses/pluginATLAS/ATLAS_2013_I1217863_W.cc deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_W.cc +++ /dev/null @@ -1,200 +0,0 @@ -// -*- C++ -*- -#include "Rivet/Analysis.hh" -#include "Rivet/Projections/FinalState.hh" -#include "Rivet/Projections/FastJets.hh" -#include "Rivet/Projections/WFinder.hh" -#include "Rivet/Projections/VetoedFinalState.hh" -#include "Rivet/Projections/LeadingParticlesFinalState.hh" - -namespace Rivet { - - - class ATLAS_2013_I1217863_W : public Analysis { - public: - - /// Constructor - ATLAS_2013_I1217863_W(string name="ATLAS_2013_I1217863_W") - : Analysis(name) - { - // the electron mode is used by default - _mode = 1; - } - - - /// @name Analysis methods - //@{ - - /// Book histograms and initialise projections before the run - void init() { - - FinalState fs; - declare(fs, "FS"); - - Cut cuts = Cuts::abseta < 2.47 && Cuts::pT > 25*GeV; - - // W finder for electrons and muons - WFinder wf(fs, cuts, _mode==3? PID::MUON : PID::ELECTRON, 0.0*GeV, 1000.0*GeV, 35.0*GeV, 0.1, - WFinder::CLUSTERNODECAY, WFinder::NOTRACK, WFinder::TRANSMASS); - declare(wf, "WF"); - - // leading photon - LeadingParticlesFinalState photonfs(FinalState(Cuts::abseta < 2.37 && Cuts::pT > 15*GeV)); - photonfs.addParticleId(PID::PHOTON); - declare(photonfs, "LeadingPhoton"); - - // jets - VetoedFinalState jet_fs(fs); - jet_fs.addVetoOnThisFinalState(getProjection("WF")); - jet_fs.addVetoOnThisFinalState(getProjection("LeadingPhoton")); - FastJets jets(jet_fs, FastJets::ANTIKT, 0.4); - jets.useInvisibles(true); - declare(jets, "Jets"); - - // FS excluding the leading photon - VetoedFinalState vfs(fs); - vfs.addVetoOnThisFinalState(photonfs); - declare(vfs, "isolatedFS"); - - - // Book histograms - _hist_EgammaT_incl = bookHisto1D( 7, 1, _mode); // dSigma / dE^gamma_T for Njet >= 0 - _hist_EgammaT_excl = bookHisto1D( 8, 1, _mode); // dSigma / dE^gamma_T for Njet = 0 - _hist_Njet_EgammaT15 = bookHisto1D(15, 1, _mode); // dSigma / dNjet for E^gamma_T > 15 - _hist_Njet_EgammaT60 = bookHisto1D(16, 1, _mode); // dSigma / dNjet for E^gamma_T > 60 - _hist_mWgammaT = bookHisto1D(19, 1, _mode); // dSigma / dm^{Wgamma}_T - - } - - - /// Perform the per-event analysis - void analyze(const Event& event) { - - const double weight = event.weight(); - - // retrieve leading photon - Particles photons = apply(event, "LeadingPhoton").particles(); - if (photons.size() != 1) vetoEvent; - const Particle& leadingPhoton = photons[0]; - if (leadingPhoton.Et() < 15.0*GeV) vetoEvent; - if (leadingPhoton.abseta() > 2.37) vetoEvent; - - // check photon isolation - double coneEnergy(0.0); - Particles fs = apply(event, "isolatedFS").particles(); - for (const Particle& p : fs) { - if ( deltaR(leadingPhoton, p) < 0.4 ) coneEnergy += p.E(); - } - if ( coneEnergy / leadingPhoton.E() >= 0.5 ) vetoEvent; - - // retrieve W boson candidate - const WFinder& wf = apply(event, "WF"); - if ( wf.bosons().size() != 1 ) vetoEvent; // only one W boson candidate - //const Particle& Wboson = wf.boson(); - - // retrieve constituent neutrino - const Particle& neutrino = wf.constituentNeutrino(); - if ( !(neutrino.pT() > 35.0*GeV) ) vetoEvent; - - // retrieve constituent lepton - const Particle& lepton = wf.constituentLepton(); - if ( !(lepton.pT() > 25.0*GeV && lepton.abseta() < 2.47) ) vetoEvent; - - // check photon-lepton overlap - if ( !(deltaR(leadingPhoton, lepton) > 0.7) ) vetoEvent; - - // count jets - const FastJets& jetfs = apply(event, "Jets"); - Jets jets = jetfs.jets(cmpMomByEt); - int goodJets = 0; - for (const Jet& j : jets) { - if ( !(j.Et() > 30.0*GeV) ) break; - if ( (j.abseta() < 4.4) && \ - (deltaR(leadingPhoton, j) > 0.3) && \ - (deltaR(lepton, j) > 0.3) ) ++goodJets; - } - - double Njets = double(goodJets) + 0.5; - double photonEt = leadingPhoton.Et()*GeV; - - const FourMomentum& lep_gamma = lepton.momentum() + leadingPhoton.momentum(); - double term1 = sqrt(lep_gamma.mass2() + lep_gamma.pT2()) + neutrino.Et(); - double term2 = (lep_gamma + neutrino.momentum()).pT2(); - double mWgammaT = sqrt(term1 * term1 - term2) * GeV; - - _hist_EgammaT_incl->fill(photonEt, weight); - - _hist_Njet_EgammaT15->fill(Njets, weight); - - if ( !goodJets ) _hist_EgammaT_excl->fill(photonEt, weight); - - if (photonEt > 40.0*GeV) { - _hist_mWgammaT->fill(mWgammaT, weight); - if (photonEt > 60.0*GeV) _hist_Njet_EgammaT60->fill(Njets, weight); - } - - } - - - /// Normalise histograms etc., after the run - void finalize() { - - const double xs_fb = crossSection()/femtobarn; - const double sumw = sumOfWeights(); - const double sf = xs_fb / sumw; - - scale(_hist_EgammaT_excl, sf); - scale(_hist_EgammaT_incl, sf); - - normalize(_hist_Njet_EgammaT15); - normalize(_hist_Njet_EgammaT60); - normalize(_hist_mWgammaT); - - } - - //@} - - protected: - - size_t _mode; - - private: - - /// @name Histograms - //@{ - - Histo1DPtr _hist_EgammaT_incl; - Histo1DPtr _hist_EgammaT_excl; - Histo1DPtr _hist_Njet_EgammaT15; - Histo1DPtr _hist_Njet_EgammaT60; - Histo1DPtr _hist_mWgammaT; - - //@} - - }; - - - class ATLAS_2013_I1217863_W_EL : public ATLAS_2013_I1217863_W { - public: - ATLAS_2013_I1217863_W_EL() - : ATLAS_2013_I1217863_W("ATLAS_2013_I1217863_W_EL") - { - _mode = 2; - } - }; - - - class ATLAS_2013_I1217863_W_MU : public ATLAS_2013_I1217863_W { - public: - ATLAS_2013_I1217863_W_MU() - : ATLAS_2013_I1217863_W("ATLAS_2013_I1217863_W_MU") - { - _mode = 3; - } - }; - - - DECLARE_RIVET_PLUGIN(ATLAS_2013_I1217863_W); - DECLARE_RIVET_PLUGIN(ATLAS_2013_I1217863_W_EL); - DECLARE_RIVET_PLUGIN(ATLAS_2013_I1217863_W_MU); - -} diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_W.info b/analyses/pluginATLAS/ATLAS_2013_I1217863_W.info deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_W.info +++ /dev/null @@ -1,45 +0,0 @@ -Name: ATLAS_2013_I1217863_W -Year: 2013 -Summary: W + gamma production at 7 TeV -Experiment: ATLAS -Collider: LHC -InspireID: 1217863 -Status: VALIDATED -Authors: - - Chritian Gutschow -References: - - Phys.Rev. D87 (2013) 112003 - - doi:10.1103/PhysRevD.87.112003 - - arXiv:1302.1283 [hep-ex] -RunInfo: - W+gamma in the electron channel -NumEvents: 1000000 -Beams: [p+, p+] -Energies: [7000] -PtCuts: [25,15] -NeedCrossSection: True -Description: - Measurements of the differential fiducial cross sections for the production of a W boson in association with a high-energy - photon are measured using pp collisions at $\sqrt{s}=7~\TeV$. The analysis uses a data sample with an integrated luminosity of - 4.6/fb collected by the ATLAS detector during the 2011 LHC data-taking period. Events are selected using leptonic decays - of the W bosons with the requirement of an associated isolated photon. The default routine will consider the electron decay - channel of the W boson. Use ATLAS_2013_I1217863_W_EL and ATLAS_2013_I1217863_W_MU to specify the decay channel directly. -BibKey: Aad:2013iz -BibTeX: '@article{Aad:2013izg, - author = "Aad, Georges and others", - title = "{Measurements of $W\gamma$ and $Z\gamma$ production in $pp$ collisions - at $\sqrt{s}$=7  TeV with the ATLAS detector at the - LHC}", - collaboration = "ATLAS Collaboration", - journal = "Phys.Rev.", - number = "11", - volume = "D87", - pages = "112003", - doi = "10.1103/PhysRevD.87.112003", - year = "2013", - eprint = "1302.1283", - archivePrefix = "arXiv", - primaryClass = "hep-ex", - reportNumber = "CERN-PH-EP-2012-345", - SLACcitation = "%%CITATION = ARXIV:1302.1283;%%", -}' diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_W.plot b/analyses/pluginATLAS/ATLAS_2013_I1217863_W.plot deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_W.plot +++ /dev/null @@ -1,50 +0,0 @@ -# BEGIN PLOT /ATLAS_2013_I1217863_W/d..-x..-y.. -LogY=1 -LegendYPos=0.90 -LegendXPos=0.90 -LegendAlign=r -RatioPlotYMax=1.5 -RatioPlotYMin=0.5 -XTwosidedTicks=1 -YTwosidedTicks=1 -ErrorBars=1 -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W/d0[7-8]-x01-y.. -LogX=1 -RatioPlotYMax=2.0 -RatioPlotYMin=0.0 -XLabel=$E^\gamma_\text{T}$ [GeV] -YLabel=$\text{d}\sigma(pp \rightarrow \ell \nu \gamma) / \text{d}E^\gamma_\text{T}$ [fb $\text{GeV}^{-1}$] -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W/d1[5-6]-x01-y.. -XLabel=Jet multiplicity -YLabel=$1/\sigma_{W\gamma} \times \text{d}\sigma_{W\gamma} / \text{d} N_\text{jet}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W/d19-x01-y.. -XLabel=$m^{W\gamma}_\text{T}$ [GeV] -YLabel=$1/\sigma_{W\gamma} \times \text{d}\sigma_{W\gamma} / \text{d} m^{W\gamma}_\text{T}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W/d07-x01-y01 -Title=$N_\text{jet} \geq 0$, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W/d08-x01-y01 -Title=$N_\text{jet} = 0$, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W/d15-x01-y01 -Title=$E^\gamma_\text{T} \geq 15$ GeV, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W/d16-x01-y01 -Title=$E^\gamma_\text{T} \geq 60$ GeV, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W/d19-x01-y01 -Title=$E^\gamma_\text{T} > 40$ GeV, electron channel, dressed level -# END PLOT - diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_W.yoda b/analyses/pluginATLAS/ATLAS_2013_I1217863_W.yoda deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_W.yoda +++ /dev/null @@ -1,70 +0,0 @@ -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W/d07-x01-y01 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W/d07-x01-y01 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.947072e+02 3.245120e+01 3.245120e+01 -2.500000e+01 5.000000e+00 5.000000e+00 8.518440e+01 1.115510e+01 1.115510e+01 -3.500000e+01 5.000000e+00 5.000000e+00 4.360630e+01 5.070500e+00 5.070500e+00 -5.000000e+01 1.000000e+01 1.000000e+01 1.409599e+01 1.825380e+00 1.825380e+00 -8.000000e+01 2.000000e+01 2.000000e+01 5.070500e+00 5.070500e-01 5.070500e-01 -5.500000e+02 4.500000e+02 4.500000e+02 9.126900e-02 1.216920e-02 1.216920e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W/d08-x01-y01 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W/d08-x01-y01 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.379176e+02 2.231020e+01 2.231020e+01 -2.500000e+01 5.000000e+00 5.000000e+00 5.526845e+01 7.200110e+00 7.200110e+00 -3.500000e+01 5.000000e+00 5.000000e+00 2.393276e+01 3.346530e+00 3.346530e+00 -5.000000e+01 1.000000e+01 1.000000e+01 6.997290e+00 1.216920e+00 1.216920e+00 -8.000000e+01 2.000000e+01 2.000000e+01 2.129610e+00 3.042300e-01 3.042300e-01 -5.500000e+02 4.500000e+02 4.500000e+02 3.042300e-02 6.084600e-03 6.084600e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W/d15-x01-y01 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W/d15-x01-y01 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 5.800000e-01 3.000000e-02 3.000000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 3.300000e-01 3.000000e-02 3.000000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 1.200000e-01 1.000000e-02 1.000000e-02 -3.500000e+00 5.000000e-01 5.000000e-01 3.700000e-02 8.000000e-03 8.000000e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W/d16-x01-y01 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W/d16-x01-y01 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 3.300000e-01 3.300000e-02 3.300000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 4.000000e-01 2.600000e-02 2.600000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 1.600000e-01 1.800000e-02 1.800000e-02 -3.500000e+00 5.000000e-01 5.000000e-01 1.100000e-01 2.400000e-02 2.400000e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W/d19-x01-y01 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W/d19-x01-y01 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -7.500000e+01 2.500000e+01 2.500000e+01 7.400000e-04 1.980000e-04 1.980000e-04 -1.250000e+02 2.500000e+01 2.500000e+01 5.400000e-03 4.400000e-04 4.400000e-04 -1.750000e+02 2.500000e+01 2.500000e+01 6.400000e-03 4.800000e-04 4.800000e-04 -2.500000e+02 5.000000e+01 5.000000e+01 2.400000e-03 1.500000e-04 1.500000e-04 -3.500000e+02 5.000000e+01 5.000000e+01 7.900000e-04 8.700000e-05 8.700000e-05 -7.000000e+02 3.000000e+02 3.000000e+02 8.333333e-05 1.416667e-05 1.416667e-05 -END YODA_SCATTER2D_V2 diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_EL.info b/analyses/pluginATLAS/ATLAS_2013_I1217863_W_EL.info deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_EL.info +++ /dev/null @@ -1,45 +0,0 @@ -Name: ATLAS_2013_I1217863_W_EL -Year: 2013 -Summary: W + gamma production at 7 TeV -Experiment: ATLAS -Collider: LHC -InspireID: 1217863 -Status: VALIDATED -Authors: - - Chritian Gutschow -References: - - Phys.Rev. D87 (2013) 112003 - - doi:10.1103/PhysRevD.87.112003 - - arXiv:1302.1283 [hep-ex] -RunInfo: - W+gamma in the electron channel -NumEvents: 1000000 -Beams: [p+, p+] -Energies: [7000] -PtCuts: [25,15] -NeedCrossSection: True -Description: - Measurements of the differential fiducial cross sections for the production of a W boson in association with a high-energy - photon are measured using pp collisions at $\sqrt{s}=7~\TeV$. The analysis uses a data sample with an integrated luminosity of - 4.6/fb collected by the ATLAS detector during the 2011 LHC data-taking period. Events are selected using leptonic decays - of the W bosons with the requirement of an associated isolated photon. The default routine will consider the electron decay - channel of the W boson. Use ATLAS_2013_I1217863_W_EL and ATLAS_2013_I1217863_W_MU to specify the decay channel directly. -BibKey: Aad:2013iz -BibTeX: '@article{Aad:2013izg, - author = "Aad, Georges and others", - title = "{Measurements of $W\gamma$ and $Z\gamma$ production in $pp$ collisions - at $\sqrt{s}$=7  TeV with the ATLAS detector at the - LHC}", - collaboration = "ATLAS Collaboration", - journal = "Phys.Rev.", - number = "11", - volume = "D87", - pages = "112003", - doi = "10.1103/PhysRevD.87.112003", - year = "2013", - eprint = "1302.1283", - archivePrefix = "arXiv", - primaryClass = "hep-ex", - reportNumber = "CERN-PH-EP-2012-345", - SLACcitation = "%%CITATION = ARXIV:1302.1283;%%", -}' diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_EL.plot b/analyses/pluginATLAS/ATLAS_2013_I1217863_W_EL.plot deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_EL.plot +++ /dev/null @@ -1,50 +0,0 @@ -# BEGIN PLOT /ATLAS_2013_I1217863_W_EL/d..-x..-y.. -LogY=1 -LegendYPos=0.90 -LegendXPos=0.90 -LegendAlign=r -RatioPlotYMax=1.5 -RatioPlotYMin=0.5 -XTwosidedTicks=1 -YTwosidedTicks=1 -ErrorBars=1 -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_EL/d0[7-8]-x01-y.. -LogX=1 -RatioPlotYMax=2.0 -RatioPlotYMin=0.0 -XLabel=$E^\gamma_\text{T}$ [GeV] -YLabel=$\text{d}\sigma(pp \rightarrow \ell \nu \gamma) / \text{d}E^\gamma_\text{T}$ [fb $\text{GeV}^{-1}$] -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_EL/d1[5-6]-x01-y.. -XLabel=Jet multiplicity -YLabel=$1/\sigma_{W\gamma} \times \text{d}\sigma_{W\gamma} / \text{d} N_\text{jet}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_EL/d19-x01-y.. -XLabel=$m^{W\gamma}_\text{T}$ [GeV] -YLabel=$1/\sigma_{W\gamma} \times \text{d}\sigma_{W\gamma} / \text{d} m^{W\gamma}_\text{T}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_EL/d07-x01-y01 -Title=$N_\text{jet} \geq 0$, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_EL/d08-x01-y01 -Title=$N_\text{jet} = 0$, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_EL/d15-x01-y01 -Title=$E^\gamma_\text{T} \geq 15$ GeV, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_EL/d16-x01-y01 -Title=$E^\gamma_\text{T} \geq 60$ GeV, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_EL/d19-x01-y01 -Title=$E^\gamma_\text{T} > 40$ GeV, electron channel, dressed level -# END PLOT - diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_EL.yoda b/analyses/pluginATLAS/ATLAS_2013_I1217863_W_EL.yoda deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_EL.yoda +++ /dev/null @@ -1,70 +0,0 @@ -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W_EL/d07-x01-y02 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W_EL/d07-x01-y02 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.947072e+02 3.245120e+01 3.245120e+01 -2.500000e+01 5.000000e+00 5.000000e+00 8.518440e+01 1.115510e+01 1.115510e+01 -3.500000e+01 5.000000e+00 5.000000e+00 4.360630e+01 5.070500e+00 5.070500e+00 -5.000000e+01 1.000000e+01 1.000000e+01 1.409599e+01 1.825380e+00 1.825380e+00 -8.000000e+01 2.000000e+01 2.000000e+01 5.070500e+00 5.070500e-01 5.070500e-01 -5.500000e+02 4.500000e+02 4.500000e+02 9.126900e-02 1.216920e-02 1.216920e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W_EL/d08-x01-y02 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W_EL/d08-x01-y02 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.379176e+02 2.231020e+01 2.231020e+01 -2.500000e+01 5.000000e+00 5.000000e+00 5.526845e+01 7.200110e+00 7.200110e+00 -3.500000e+01 5.000000e+00 5.000000e+00 2.393276e+01 3.346530e+00 3.346530e+00 -5.000000e+01 1.000000e+01 1.000000e+01 6.997290e+00 1.216920e+00 1.216920e+00 -8.000000e+01 2.000000e+01 2.000000e+01 2.129610e+00 3.042300e-01 3.042300e-01 -5.500000e+02 4.500000e+02 4.500000e+02 3.042300e-02 6.084600e-03 6.084600e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W_EL/d15-x01-y02 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W_EL/d15-x01-y02 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 5.800000e-01 3.000000e-02 3.000000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 3.300000e-01 3.000000e-02 3.000000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 1.200000e-01 1.000000e-02 1.000000e-02 -3.500000e+00 5.000000e-01 5.000000e-01 3.700000e-02 8.000000e-03 8.000000e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W_EL/d16-x01-y02 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W_EL/d16-x01-y02 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 3.300000e-01 3.300000e-02 3.300000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 4.000000e-01 2.600000e-02 2.600000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 1.600000e-01 1.800000e-02 1.800000e-02 -3.500000e+00 5.000000e-01 5.000000e-01 1.100000e-01 2.400000e-02 2.400000e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W_EL/d19-x01-y02 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W_EL/d19-x01-y02 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -7.500000e+01 2.500000e+01 2.500000e+01 7.400000e-04 1.980000e-04 1.980000e-04 -1.250000e+02 2.500000e+01 2.500000e+01 5.400000e-03 4.400000e-04 4.400000e-04 -1.750000e+02 2.500000e+01 2.500000e+01 6.400000e-03 4.800000e-04 4.800000e-04 -2.500000e+02 5.000000e+01 5.000000e+01 2.400000e-03 1.500000e-04 1.500000e-04 -3.500000e+02 5.000000e+01 5.000000e+01 7.900000e-04 8.700000e-05 8.700000e-05 -7.000000e+02 3.000000e+02 3.000000e+02 8.333333e-05 1.416667e-05 1.416667e-05 -END YODA_SCATTER2D_V2 diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_MU.info b/analyses/pluginATLAS/ATLAS_2013_I1217863_W_MU.info deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_MU.info +++ /dev/null @@ -1,45 +0,0 @@ -Name: ATLAS_2013_I1217863_W_MU -Year: 2013 -Summary: W + gamma production at 7 TeV -Experiment: ATLAS -Collider: LHC -InspireID: 1217863 -Status: VALIDATED -Authors: - - Chritian Gutschow -References: - - Phys.Rev. D87 (2013) 112003 - - doi:10.1103/PhysRevD.87.112003 - - arXiv:1302.1283 [hep-ex] -RunInfo: - W+gamma in the muon channel -NumEvents: 1000000 -Beams: [p+, p+] -Energies: [7000] -PtCuts: [25,15] -NeedCrossSection: True -Description: - Measurements of the differential fiducial cross sections for the production of a W boson in association with a high-energy - photon are measured using pp collisions at $\sqrt{s}=7~\TeV$. The analysis uses a data sample with an integrated luminosity of - 4.6/fb collected by the ATLAS detector during the 2011 LHC data-taking period. Events are selected using leptonic decays - of the W bosons with the requirement of an associated isolated photon. The default routine will consider the electron decay - channel of the W boson. Use ATLAS_2013_I1217863_W_EL and ATLAS_2013_I1217863_W_MU to specify the decay channel directly. -BibKey: Aad:2013iz -BibTeX: '@article{Aad:2013izg, - author = "Aad, Georges and others", - title = "{Measurements of $W\gamma$ and $Z\gamma$ production in $pp$ collisions - at $\sqrt{s}$=7  TeV with the ATLAS detector at the - LHC}", - collaboration = "ATLAS Collaboration", - journal = "Phys.Rev.", - number = "11", - volume = "D87", - pages = "112003", - doi = "10.1103/PhysRevD.87.112003", - year = "2013", - eprint = "1302.1283", - archivePrefix = "arXiv", - primaryClass = "hep-ex", - reportNumber = "CERN-PH-EP-2012-345", - SLACcitation = "%%CITATION = ARXIV:1302.1283;%%", -}' diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_MU.plot b/analyses/pluginATLAS/ATLAS_2013_I1217863_W_MU.plot deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_MU.plot +++ /dev/null @@ -1,50 +0,0 @@ -# BEGIN PLOT /ATLAS_2013_I1217863_W_MU/d..-x..-y.. -LogY=1 -LegendYPos=0.90 -LegendXPos=0.90 -LegendAlign=r -RatioPlotYMax=1.5 -RatioPlotYMin=0.5 -XTwosidedTicks=1 -YTwosidedTicks=1 -ErrorBars=1 -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_MU/d0[7-8]-x01-y.. -LogX=1 -RatioPlotYMax=2.0 -RatioPlotYMin=0.0 -XLabel=$E^\gamma_\text{T}$ [GeV] -YLabel=$\text{d}\sigma(pp \rightarrow \ell \nu \gamma) / \text{d}E^\gamma_\text{T}$ [fb $\text{GeV}^{-1}$] -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_MU/d1[5-6]-x01-y.. -XLabel=Jet multiplicity -YLabel=$1/\sigma_{W\gamma} \times \text{d}\sigma_{W\gamma} / \text{d} N_\text{jet}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_MU/d19-x01-y.. -XLabel=$m^{W\gamma}_\text{T}$ [GeV] -YLabel=$1/\sigma_{W\gamma} \times \text{d}\sigma_{W\gamma} / \text{d} m^{W\gamma}_\text{T}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_MU/d07-x01-y01 -Title=$N_\text{jet} \geq 0$, muon channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_MU/d08-x01-y01 -Title=$N_\text{jet} = 0$, muon channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_MU/d15-x01-y01 -Title=$E^\gamma_\text{T} \geq 15$ GeV, muon channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_MU/d16-x01-y01 -Title=$E^\gamma_\text{T} \geq 60$ GeV, muon channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_W_MU/d19-x01-y01 -Title=$E^\gamma_\text{T} > 40$ GeV, muon channel, dressed level -# END PLOT - diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_MU.yoda b/analyses/pluginATLAS/ATLAS_2013_I1217863_W_MU.yoda deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_W_MU.yoda +++ /dev/null @@ -1,70 +0,0 @@ -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W_MU/d07-x01-y03 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W_MU/d07-x01-y03 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.947072e+02 3.245120e+01 3.245120e+01 -2.500000e+01 5.000000e+00 5.000000e+00 8.518440e+01 1.115510e+01 1.115510e+01 -3.500000e+01 5.000000e+00 5.000000e+00 4.360630e+01 5.070500e+00 5.070500e+00 -5.000000e+01 1.000000e+01 1.000000e+01 1.409599e+01 1.825380e+00 1.825380e+00 -8.000000e+01 2.000000e+01 2.000000e+01 5.070500e+00 5.070500e-01 5.070500e-01 -5.500000e+02 4.500000e+02 4.500000e+02 9.126900e-02 1.216920e-02 1.216920e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W_MU/d08-x01-y03 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W_MU/d08-x01-y03 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.379176e+02 2.231020e+01 2.231020e+01 -2.500000e+01 5.000000e+00 5.000000e+00 5.526845e+01 7.200110e+00 7.200110e+00 -3.500000e+01 5.000000e+00 5.000000e+00 2.393276e+01 3.346530e+00 3.346530e+00 -5.000000e+01 1.000000e+01 1.000000e+01 6.997290e+00 1.216920e+00 1.216920e+00 -8.000000e+01 2.000000e+01 2.000000e+01 2.129610e+00 3.042300e-01 3.042300e-01 -5.500000e+02 4.500000e+02 4.500000e+02 3.042300e-02 6.084600e-03 6.084600e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W_MU/d15-x01-y03 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W_MU/d15-x01-y03 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 5.800000e-01 3.000000e-02 3.000000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 3.300000e-01 3.000000e-02 3.000000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 1.200000e-01 1.000000e-02 1.000000e-02 -3.500000e+00 5.000000e-01 5.000000e-01 3.700000e-02 8.000000e-03 8.000000e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W_MU/d16-x01-y03 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W_MU/d16-x01-y03 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 3.300000e-01 3.300000e-02 3.300000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 4.000000e-01 2.600000e-02 2.600000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 1.600000e-01 1.800000e-02 1.800000e-02 -3.500000e+00 5.000000e-01 5.000000e-01 1.100000e-01 2.400000e-02 2.400000e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_W_MU/d19-x01-y03 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_W_MU/d19-x01-y03 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -7.500000e+01 2.500000e+01 2.500000e+01 7.400000e-04 1.980000e-04 1.980000e-04 -1.250000e+02 2.500000e+01 2.500000e+01 5.400000e-03 4.400000e-04 4.400000e-04 -1.750000e+02 2.500000e+01 2.500000e+01 6.400000e-03 4.800000e-04 4.800000e-04 -2.500000e+02 5.000000e+01 5.000000e+01 2.400000e-03 1.500000e-04 1.500000e-04 -3.500000e+02 5.000000e+01 5.000000e+01 7.900000e-04 8.700000e-05 8.700000e-05 -7.000000e+02 3.000000e+02 3.000000e+02 8.333333e-05 1.416667e-05 1.416667e-05 -END YODA_SCATTER2D_V2 diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z.plot b/analyses/pluginATLAS/ATLAS_2013_I1217863_Z.plot deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z.plot +++ /dev/null @@ -1,50 +0,0 @@ -# BEGIN PLOT /ATLAS_2013_I1217863_Z/d..-x..-y.. -LogY=1 -LegendYPos=0.90 -LegendXPos=0.90 -LegendAlign=r -RatioPlotYMax=1.5 -RatioPlotYMin=0.5 -XTwosidedTicks=1 -YTwosidedTicks=1 -ErrorBars=1 -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z/d1[1-2]-x01-y.. -LogX=1 -RatioPlotYMax=2.0 -RatioPlotYMin=0.0 -XLabel=$E^\gamma_\text{T}$ [GeV] -YLabel=$\text{d}\sigma(pp \rightarrow \ell^+ \ell^- \gamma) / \text{d}E^\gamma_\text{T}$ [fb $\text{GeV}^{-1}$] -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z/d1[7-8]-x01-y.. -XLabel=Jet multiplicity -YLabel=$1/\sigma_{Z\gamma} \times \text{d}\sigma_{Z\gamma} / \text{d} N_\text{jet}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z/d20-x01-y.. -XLabel=$m^{Z\gamma}$ [GeV] -YLabel=$1/\sigma_{Z\gamma} \times \text{d}\sigma_{Z\gamma} / \text{d} m^{Z\gamma}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z/d11-x01-y01 -Title=$N_\text{jet} \geq 0$, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z/d12-x01-y01 -Title=$N_\text{jet} = 0$, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z/d17-x01-y01 -Title=$E^\gamma_\text{T} \geq 15$ GeV, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z/d18-x01-y01 -Title=$E^\gamma_\text{T} \geq 60$ GeV, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z/d20-x01-y01 -Title=$E^\gamma_\text{T} > 40$ GeV, electron channel, dressed level -# END PLOT - diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z.yoda b/analyses/pluginATLAS/ATLAS_2013_I1217863_Z.yoda deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z.yoda +++ /dev/null @@ -1,66 +0,0 @@ -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z/d11-x01-y01 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z/d11-x01-y01 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.216920e+02 1.216920e+01 1.216920e+01 -2.500000e+01 5.000000e+00 5.000000e+00 4.309925e+01 4.259220e+00 4.259220e+00 -3.500000e+01 5.000000e+00 5.000000e+00 1.318330e+01 1.419740e+00 1.419740e+00 -5.000000e+01 1.000000e+01 1.000000e+01 5.009654e+00 6.186010e-01 6.186010e-01 -8.000000e+01 2.000000e+01 2.000000e+01 1.419740e+00 1.926790e-01 1.926790e-01 -5.500000e+02 4.500000e+02 4.500000e+02 1.825380e-02 7.098700e-03 7.098700e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z/d12-x01-y01 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z/d12-x01-y01 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.074946e+02 1.115510e+01 1.115510e+01 -2.500000e+01 5.000000e+00 5.000000e+00 3.478363e+01 4.157810e+00 4.157810e+00 -3.500000e+01 5.000000e+00 5.000000e+00 9.431130e+00 1.115510e+00 1.115510e+00 -5.000000e+01 1.000000e+01 1.000000e+01 3.285684e+00 4.766270e-01 4.766270e-01 -8.000000e+01 2.000000e+01 2.000000e+01 1.024241e+00 1.622560e-01 1.622560e-01 -5.500000e+02 4.500000e+02 4.500000e+02 7.098700e-03 4.056400e-03 4.056400e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z/d17-x01-y01 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z/d17-x01-y01 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 8.000000e-01 2.700000e-02 2.700000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 1.500000e-01 1.600000e-02 1.600000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 5.200000e-02 1.100000e-02 1.100000e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z/d18-x01-y01 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z/d18-x01-y01 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 6.000000e-01 3.800000e-02 3.800000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 2.800000e-01 2.600000e-02 2.600000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 1.200000e-01 1.900000e-02 1.900000e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z/d20-x01-y01 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z/d20-x01-y01 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.000000e+02 5.000000e+01 5.000000e+01 1.400000e-03 1.800000e-04 1.800000e-04 -1.750000e+02 2.500000e+01 2.500000e+01 8.800000e-03 4.800000e-04 4.800000e-04 -2.500000e+02 5.000000e+01 5.000000e+01 3.200000e-03 2.200000e-04 2.200000e-04 -6.500000e+02 3.500000e+02 3.500000e+02 1.571429e-04 2.142857e-05 2.142857e-05 -END YODA_SCATTER2D_V2 diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_EL.info b/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_EL.info deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_EL.info +++ /dev/null @@ -1,45 +0,0 @@ -Name: ATLAS_2013_I1217863_Z_EL -Year: 2013 -Summary: Z + gamma production at 7 TeV -Experiment: ATLAS -Collider: LHC -InspireID: 1217863 -Status: VALIDATED -Authors: - - Chritian Gutschow -References: - - Phys.Rev. D87 (2013) 112003 - - doi:10.1103/PhysRevD.87.112003 - - arXiv:1302.1283 [hep-ex] -RunInfo: - Z+gamma in the electron channel -NumEvents: 1000000 -Beams: [p+, p+] -Energies: [7000] -PtCuts: [25,15] -NeedCrossSection: True -Description: - Measurements of the differential fiducial cross sections for the production of a Z boson in association with a high-energy - photon are measured using pp collisions at $\sqrt{s}=7~\TeV$. The analysis uses a data sample with an integrated luminosity of - 4.6 fb-1 collected by the ATLAS detector during the 2011 LHC data-taking period. Events are selected using leptonic decays - of the Z bosons with the requirement of an associated isolated photon. The default routine will consider the electron decay - channel of the Z boson. Use ATLAS_2013_I1217863_Z_EL and ATLAS_2013_I1217863_Z_MU to specify the decay channel directly. -BibKey: Aad:2013iz -BibTeX: '@article{Aad:2013izg, - author = "Aad, Georges and others", - title = "{Measurements of $W\gamma$ and $Z\gamma$ production in $pp$ collisions - at $\sqrt{s}$=7  TeV with the ATLAS detector at the - LHC}", - collaboration = "ATLAS Collaboration", - journal = "Phys.Rev.", - number = "11", - volume = "D87", - pages = "112003", - doi = "10.1103/PhysRevD.87.112003", - year = "2013", - eprint = "1302.1283", - archivePrefix = "arXiv", - primaryClass = "hep-ex", - reportNumber = "CERN-PH-EP-2012-345", - SLACcitation = "%%CITATION = ARXIV:1302.1283;%%", -}' diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_EL.plot b/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_EL.plot deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_EL.plot +++ /dev/null @@ -1,50 +0,0 @@ -# BEGIN PLOT /ATLAS_2013_I1217863_Z_EL/d..-x..-y.. -LogY=1 -LegendYPos=0.90 -LegendXPos=0.90 -LegendAlign=r -RatioPlotYMax=1.5 -RatioPlotYMin=0.5 -XTwosidedTicks=1 -YTwosidedTicks=1 -ErrorBars=1 -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_EL/d1[1-2]-x01-y.. -LogX=1 -RatioPlotYMax=2.0 -RatioPlotYMin=0.0 -XLabel=$E^\gamma_\text{T}$ [GeV] -YLabel=$\text{d}\sigma(pp \rightarrow \ell^+ \ell^- \gamma) / \text{d}E^\gamma_\text{T}$ [fb $\text{GeV}^{-1}$] -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_EL/d1[7-8]-x01-y.. -XLabel=Jet multiplicity -YLabel=$1/\sigma_{Z\gamma} \times \text{d}\sigma_{Z\gamma} / \text{d} N_\text{jet}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_EL/d20-x01-y.. -XLabel=$m^{Z\gamma}$ [GeV] -YLabel=$1/\sigma_{Z\gamma} \times \text{d}\sigma_{Z\gamma} / \text{d} m^{Z\gamma}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_EL/d11-x01-y01 -Title=$N_\text{jet} \geq 0$, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_EL/d12-x01-y01 -Title=$N_\text{jet} = 0$, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_EL/d17-x01-y01 -Title=$E^\gamma_\text{T} \geq 15$ GeV, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_EL/d18-x01-y01 -Title=$E^\gamma_\text{T} \geq 60$ GeV, electron channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_EL/d20-x01-y01 -Title=$E^\gamma_\text{T} > 40$ GeV, electron channel, dressed level -# END PLOT - diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_EL.yoda b/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_EL.yoda deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_EL.yoda +++ /dev/null @@ -1,66 +0,0 @@ -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z_EL/d11-x01-y02 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z_EL/d11-x01-y02 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.216920e+02 1.216920e+01 1.216920e+01 -2.500000e+01 5.000000e+00 5.000000e+00 4.309925e+01 4.259220e+00 4.259220e+00 -3.500000e+01 5.000000e+00 5.000000e+00 1.318330e+01 1.419740e+00 1.419740e+00 -5.000000e+01 1.000000e+01 1.000000e+01 5.009654e+00 6.186010e-01 6.186010e-01 -8.000000e+01 2.000000e+01 2.000000e+01 1.419740e+00 1.926790e-01 1.926790e-01 -5.500000e+02 4.500000e+02 4.500000e+02 1.825380e-02 7.098700e-03 7.098700e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z_EL/d12-x01-y02 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z_EL/d12-x01-y02 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.074946e+02 1.115510e+01 1.115510e+01 -2.500000e+01 5.000000e+00 5.000000e+00 3.478363e+01 4.157810e+00 4.157810e+00 -3.500000e+01 5.000000e+00 5.000000e+00 9.431130e+00 1.115510e+00 1.115510e+00 -5.000000e+01 1.000000e+01 1.000000e+01 3.285684e+00 4.766270e-01 4.766270e-01 -8.000000e+01 2.000000e+01 2.000000e+01 1.024241e+00 1.622560e-01 1.622560e-01 -5.500000e+02 4.500000e+02 4.500000e+02 7.098700e-03 4.056400e-03 4.056400e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z_EL/d17-x01-y02 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z_EL/d17-x01-y02 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 8.000000e-01 2.700000e-02 2.700000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 1.500000e-01 1.600000e-02 1.600000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 5.200000e-02 1.100000e-02 1.100000e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z_EL/d18-x01-y02 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z_EL/d18-x01-y02 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 6.000000e-01 3.800000e-02 3.800000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 2.800000e-01 2.600000e-02 2.600000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 1.200000e-01 1.900000e-02 1.900000e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z_EL/d20-x01-y02 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z_EL/d20-x01-y02 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.000000e+02 5.000000e+01 5.000000e+01 1.400000e-03 1.800000e-04 1.800000e-04 -1.750000e+02 2.500000e+01 2.500000e+01 8.800000e-03 4.800000e-04 4.800000e-04 -2.500000e+02 5.000000e+01 5.000000e+01 3.200000e-03 2.200000e-04 2.200000e-04 -6.500000e+02 3.500000e+02 3.500000e+02 1.571429e-04 2.142857e-05 2.142857e-05 -END YODA_SCATTER2D_V2 diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_MU.info b/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_MU.info deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_MU.info +++ /dev/null @@ -1,45 +0,0 @@ -Name: ATLAS_2013_I1217863_Z_MU -Year: 2013 -Summary: Z + gamma production at 7 TeV -Experiment: ATLAS -Collider: LHC -InspireID: 1217863 -Status: VALIDATED -Authors: - - Chritian Gutschow -References: - - Phys.Rev. D87 (2013) 112003 - - doi:10.1103/PhysRevD.87.112003 - - arXiv:1302.1283 [hep-ex] -RunInfo: - Z+gamma in the muon channel -NumEvents: 1000000 -Beams: [p+, p+] -Energies: [7000] -PtCuts: [25,15] -NeedCrossSection: True -Description: - Measurements of the differential fiducial cross sections for the production of a Z boson in association with a high-energy - photon are measured using pp collisions at $\sqrt{s}=7~\TeV$. The analysis uses a data sample with an integrated luminosity of - 4.6/fb collected by the ATLAS detector during the 2011 LHC data-taking period. Events are selected using leptonic decays - of the Z bosons with the requirement of an associated isolated photon. The default routine will consider the electron decay - channel of the Z boson. Use ATLAS_2013_I1217863_Z_EL and ATLAS_2013_I1217863_Z_MU to specify the decay channel directly. -BibKey: Aad:2013iz -BibTeX: '@article{Aad:2013izg, - author = "Aad, Georges and others", - title = "{Measurements of $W\gamma$ and $Z\gamma$ production in $pp$ collisions - at $\sqrt{s}$=7  TeV with the ATLAS detector at the - LHC}", - collaboration = "ATLAS Collaboration", - journal = "Phys.Rev.", - number = "11", - volume = "D87", - pages = "112003", - doi = "10.1103/PhysRevD.87.112003", - year = "2013", - eprint = "1302.1283", - archivePrefix = "arXiv", - primaryClass = "hep-ex", - reportNumber = "CERN-PH-EP-2012-345", - SLACcitation = "%%CITATION = ARXIV:1302.1283;%%", -}' diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_MU.plot b/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_MU.plot deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_MU.plot +++ /dev/null @@ -1,50 +0,0 @@ -# BEGIN PLOT /ATLAS_2013_I1217863_Z_MU/d..-x..-y.. -LogY=1 -LegendYPos=0.90 -LegendXPos=0.90 -LegendAlign=r -RatioPlotYMax=1.5 -RatioPlotYMin=0.5 -XTwosidedTicks=1 -YTwosidedTicks=1 -ErrorBars=1 -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_MU/d1[1-2]-x01-y.. -LogX=1 -RatioPlotYMax=2.0 -RatioPlotYMin=0.0 -XLabel=$E^\gamma_\text{T}$ [GeV] -YLabel=$\text{d}\sigma(pp \rightarrow \ell^+ \ell^- \gamma) / \text{d}E^\gamma_\text{T}$ [fb $\text{GeV}^{-1}$] -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_MU/d1[7-8]-x01-y.. -XLabel=Jet multiplicity -YLabel=$1/\sigma_{Z\gamma} \times \text{d}\sigma_{Z\gamma} / \text{d} N_\text{jet}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_MU/d20-x01-y.. -XLabel=$m^{Z\gamma}$ [GeV] -YLabel=$1/\sigma_{Z\gamma} \times \text{d}\sigma_{Z\gamma} / \text{d} m^{Z\gamma}$ -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_MU/d11-x01-y01 -Title=$N_\text{jet} \geq 0$, muon channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_MU/d12-x01-y01 -Title=$N_\text{jet} = 0$, muon channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_MU/d17-x01-y01 -Title=$E^\gamma_\text{T} \geq 15$ GeV, muon channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_MU/d18-x01-y01 -Title=$E^\gamma_\text{T} \geq 60$ GeV, muon channel, dressed level -# END PLOT - -# BEGIN PLOT /ATLAS_2013_I1217863_Z_MU/d20-x01-y01 -Title=$E^\gamma_\text{T} > 40$ GeV, muon channel, dressed level -# END PLOT - diff --git a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_MU.yoda b/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_MU.yoda deleted file mode 100644 --- a/analyses/pluginATLAS/ATLAS_2013_I1217863_Z_MU.yoda +++ /dev/null @@ -1,66 +0,0 @@ -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z_MU/d11-x01-y03 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z_MU/d11-x01-y03 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.216920e+02 1.216920e+01 1.216920e+01 -2.500000e+01 5.000000e+00 5.000000e+00 4.309925e+01 4.259220e+00 4.259220e+00 -3.500000e+01 5.000000e+00 5.000000e+00 1.318330e+01 1.419740e+00 1.419740e+00 -5.000000e+01 1.000000e+01 1.000000e+01 5.009654e+00 6.186010e-01 6.186010e-01 -8.000000e+01 2.000000e+01 2.000000e+01 1.419740e+00 1.926790e-01 1.926790e-01 -5.500000e+02 4.500000e+02 4.500000e+02 1.825380e-02 7.098700e-03 7.098700e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z_MU/d12-x01-y03 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z_MU/d12-x01-y03 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.750000e+01 2.500000e+00 2.500000e+00 1.074946e+02 1.115510e+01 1.115510e+01 -2.500000e+01 5.000000e+00 5.000000e+00 3.478363e+01 4.157810e+00 4.157810e+00 -3.500000e+01 5.000000e+00 5.000000e+00 9.431130e+00 1.115510e+00 1.115510e+00 -5.000000e+01 1.000000e+01 1.000000e+01 3.285684e+00 4.766270e-01 4.766270e-01 -8.000000e+01 2.000000e+01 2.000000e+01 1.024241e+00 1.622560e-01 1.622560e-01 -5.500000e+02 4.500000e+02 4.500000e+02 7.098700e-03 4.056400e-03 4.056400e-03 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z_MU/d17-x01-y03 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z_MU/d17-x01-y03 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 8.000000e-01 2.700000e-02 2.700000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 1.500000e-01 1.600000e-02 1.600000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 5.200000e-02 1.100000e-02 1.100000e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z_MU/d18-x01-y03 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z_MU/d18-x01-y03 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -5.000000e-01 5.000000e-01 5.000000e-01 6.000000e-01 3.800000e-02 3.800000e-02 -1.500000e+00 5.000000e-01 5.000000e-01 2.800000e-01 2.600000e-02 2.600000e-02 -2.500000e+00 5.000000e-01 5.000000e-01 1.200000e-01 1.900000e-02 1.900000e-02 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2013_I1217863_Z_MU/d20-x01-y03 -IsRef: 1 -Path: /REF/ATLAS_2013_I1217863_Z_MU/d20-x01-y03 -Title: ~ -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -1.000000e+02 5.000000e+01 5.000000e+01 1.400000e-03 1.800000e-04 1.800000e-04 -1.750000e+02 2.500000e+01 2.500000e+01 8.800000e-03 4.800000e-04 4.800000e-04 -2.500000e+02 5.000000e+01 5.000000e+01 3.200000e-03 2.200000e-04 2.200000e-04 -6.500000e+02 3.500000e+02 3.500000e+02 1.571429e-04 2.142857e-05 2.142857e-05 -END YODA_SCATTER2D_V2 diff --git a/analyses/pluginATLAS/ATLAS_2014_I1306294.cc b/analyses/pluginATLAS/ATLAS_2014_I1306294.cc --- a/analyses/pluginATLAS/ATLAS_2014_I1306294.cc +++ b/analyses/pluginATLAS/ATLAS_2014_I1306294.cc @@ -1,212 +1,212 @@ // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Projections/FinalState.hh" #include "Rivet/Projections/ZFinder.hh" #include "Rivet/Projections/FastJets.hh" #include "Rivet/Projections/HeavyHadrons.hh" #include "Rivet/Projections/VetoedFinalState.hh" namespace Rivet { class ATLAS_2014_I1306294 : public Analysis { public: /// @name Constructors etc. //@{ /// Constructors ///@brief: Electroweak Wjj production at 8 TeV DEFAULT_RIVET_ANALYSIS_CTOR(ATLAS_2014_I1306294); //@} public: /// @name Analysis methods //@{ /// Book histograms and initialise projections before the run void init() { // Get options from the new option system - _mode = 0; + _mode = 1; if ( getOption("LMODE") == "EL" ) _mode = 1; if ( getOption("LMODE") == "MU" ) _mode = 2; FinalState fs; Cut cuts = Cuts::etaIn(-2.5,2.5) & (Cuts::pT > 20.0*GeV); ZFinder zfinder(fs, cuts, _mode==1? PID::ELECTRON : PID::MUON, 76.0*GeV, 106.0*GeV, 0.1, ZFinder::CLUSTERNODECAY, ZFinder::NOTRACK); declare(zfinder, "ZFinder"); //FastJets jetpro1( getProjection("ZFinder").remainingFinalState(), FastJets::ANTIKT, 0.4); VetoedFinalState jet_fs(fs); jet_fs.addVetoOnThisFinalState(getProjection("ZFinder")); FastJets jetpro1(jet_fs, FastJets::ANTIKT, 0.4); jetpro1.useInvisibles(); declare(jetpro1, "AntiKtJets04"); declare(HeavyHadrons(), "BHadrons"); //Histograms with data binning _h_bjet_Pt = bookHisto1D( 3, 1, 1); _h_bjet_Y = bookHisto1D( 5, 1, 1); _h_bjet_Yboost = bookHisto1D( 7, 1, 1); _h_bjet_DY20 = bookHisto1D( 9, 1, 1); _h_bjet_ZdPhi20 = bookHisto1D(11, 1, 1); _h_bjet_ZdR20 = bookHisto1D(13, 1, 1); _h_bjet_ZPt = bookHisto1D(15, 1, 1); _h_bjet_ZY = bookHisto1D(17, 1, 1); _h_2bjet_dR = bookHisto1D(21, 1, 1); _h_2bjet_Mbb = bookHisto1D(23, 1, 1); _h_2bjet_ZPt = bookHisto1D(25, 1, 1); _h_2bjet_ZY = bookHisto1D(27, 1, 1); } //========================================================================================== /// Perform the per-event analysis void analyze(const Event& e) { //--------------------------- const double weight = e.weight(); // -- check we have a Z: const ZFinder& zfinder = apply(e, "ZFinder"); if(zfinder.bosons().size() != 1) vetoEvent; const ParticleVector boson_s = zfinder.bosons(); const Particle boson_f = boson_s[0] ; const ParticleVector zleps = zfinder.constituents(); //--------------------------- //--------------------------- //------------- stop processing the event if no true b-partons or hadrons are found const Particles& allBs = apply(e, "BHadrons").bHadrons(5.0*GeV); Particles stableBs; foreach(Particle p, allBs) { if(p.abseta() < 2.5) stableBs += p; } if( stableBs.empty() ) vetoEvent; //--------------------------- // -- get the b-jets: const Jets& jets = apply(e, "AntiKtJets04").jetsByPt(Cuts::pT >20.0*GeV && Cuts::abseta <2.4); Jets b_jets; foreach(const Jet& jet, jets) { //veto overlaps with Z leptons: bool veto = false; foreach(const Particle& zlep, zleps) { if(deltaR(jet, zlep) < 0.5) veto = true; } if(veto) continue; foreach(const Particle& bhadron, stableBs) { if( deltaR(jet, bhadron) <= 0.3 ) { b_jets.push_back(jet); break; // match } } // end loop on b-hadrons } //and make sure we have at least 1: if(b_jets.empty()) vetoEvent; //--------------------------- // fill the plots: const double ZpT = boson_f.pT()/GeV; const double ZY = boson_f.absrap(); _h_bjet_ZPt->fill(ZpT, weight); _h_bjet_ZY ->fill(ZY, weight); foreach(const Jet& jet, b_jets) { _h_bjet_Pt->fill(jet.pT()/GeV, weight ); _h_bjet_Y ->fill(jet.absrap(), weight ); const double Yboost = 0.5 * fabs(boson_f.rapidity() + jet.rapidity()); _h_bjet_Yboost->fill(Yboost, weight ); if(ZpT > 20.) { const double ZBDY = fabs( boson_f.rapidity() - jet.rapidity() ); const double ZBDPHI = fabs( deltaPhi(jet.phi(), boson_f.phi()) ); const double ZBDR = deltaR(jet, boson_f, RAPIDITY); _h_bjet_DY20->fill( ZBDY, weight); _h_bjet_ZdPhi20->fill(ZBDPHI, weight); _h_bjet_ZdR20->fill( ZBDR, weight); } } //loop over b-jets if (b_jets.size() < 2) return; _h_2bjet_ZPt->fill(ZpT, weight); _h_2bjet_ZY ->fill(ZY, weight); const double BBDR = deltaR(b_jets[0], b_jets[1], RAPIDITY); const double Mbb = (b_jets[0].momentum() + b_jets[1].momentum()).mass(); _h_2bjet_dR ->fill(BBDR, weight); _h_2bjet_Mbb->fill(Mbb, weight); } // end of analysis loop /// Normalise histograms etc., after the run void finalize() { const double normfac = crossSection() / sumOfWeights(); scale( _h_bjet_Pt, normfac); scale( _h_bjet_Y, normfac); scale( _h_bjet_Yboost, normfac); scale( _h_bjet_DY20, normfac); scale( _h_bjet_ZdPhi20, normfac); scale( _h_bjet_ZdR20, normfac); scale( _h_bjet_ZPt, normfac); scale( _h_bjet_ZY, normfac); scale( _h_2bjet_dR, normfac); scale( _h_2bjet_Mbb, normfac); scale( _h_2bjet_ZPt, normfac); scale( _h_2bjet_ZY, normfac); } //@} protected: // Data members like post-cuts event weight counters go here size_t _mode; private: Histo1DPtr _h_bjet_Pt; Histo1DPtr _h_bjet_Y; Histo1DPtr _h_bjet_Yboost; Histo1DPtr _h_bjet_DY20; Histo1DPtr _h_bjet_ZdPhi20; Histo1DPtr _h_bjet_ZdR20; Histo1DPtr _h_bjet_ZPt; Histo1DPtr _h_bjet_ZY; Histo1DPtr _h_2bjet_dR; Histo1DPtr _h_2bjet_Mbb; Histo1DPtr _h_2bjet_ZPt; Histo1DPtr _h_2bjet_ZY; }; // The hook for the plugin system DECLARE_RIVET_PLUGIN(ATLAS_2014_I1306294); }