Page MenuHomeHEPForge
Files F7877856

No OneTemporary
Actions

View Options

diff --git a/analyses/pluginATLAS/ATLAS_2014_I1298023.cc b/analyses/pluginATLAS/ATLAS_2014_I1298023.cc
--- a/analyses/pluginATLAS/ATLAS_2014_I1298023.cc
+++ b/analyses/pluginATLAS/ATLAS_2014_I1298023.cc
@@ -1,128 +1,128 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/ChargedLeptons.hh"
#include "Rivet/Projections/VetoedFinalState.hh"
#include "Rivet/Projections/MissingMomentum.hh"
#include "Rivet/Projections/FastJets.hh"
#include "Rivet/Projections/DressedLeptons.hh"
namespace Rivet {
/// ATLAS same-sign WW production at 8 TeV (PRL)
class ATLAS_2014_I1298023 : public Analysis {
public:
/// @name Constructors etc.
//@{
ATLAS_2014_I1298023()
: Analysis("ATLAS_2014_I1298023")
{
setNeedsCrossSection(true);
}
//@}
public:
/// @name Analysis methods
//@{
void init() {
const FinalState fs;
// bare leptons
ChargedLeptons bare_leptons(fs);
// dressed leptons
Cut cuts = (Cuts::abseta < 2.5) & (Cuts::pT > 25*GeV);
DressedLeptons leptons(fs, bare_leptons, 0.1, cuts);
declare(leptons, "leptons");
// MET
declare(MissingMomentum(fs), "MissingET");
// jets
VetoedFinalState vfs(fs);
vfs.addVetoPairId(PID::MUON);
vfs.vetoNeutrinos();
declare(FastJets(vfs, FastJets::ANTIKT, 0.4), "jets");
// book histogram
_hist = bookHisto1D(1, 1, 1);
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const double weight = event.weight();
const vector<DressedLepton>& leptons = apply<DressedLeptons>(event, "leptons").dressedLeptons();
if ( leptons.size() < 2 ) vetoEvent;
double minDR_ll = MAXDOUBLE, mll = -1.0;
for (unsigned int i = 0; i < leptons.size(); ++i) {
DressedLepton lep1 = leptons[i];
for (unsigned int j = i + 1; j < leptons.size(); ++j) {
DressedLepton lep2 = leptons[j];
double dr = deltaR(lep1, lep2);
if ( dr < minDR_ll ) minDR_ll = dr;
double m = (lep1.momentum() + lep2.momentum()).mass();
if ( mll < 0. || m < mll ) mll = m;
}
}
if ( minDR_ll <= 0.3 || mll <= 20*GeV ) vetoEvent;
if ( leptons[0].charge() * leptons[1].charge() < 0.0 ) vetoEvent;
const MissingMomentum& met = apply<MissingMomentum>(event, "MissingET");
if ( met.vectorEt().mod() <= 40*GeV ) vetoEvent;
const Jets& all_jets = apply<FastJets>(event, "jets").jetsByPt( (Cuts::abseta < 4.5) && (Cuts::pT > 30*GeV) );
Jets jets;
double minDR_overall = MAXDOUBLE;
foreach (const Jet& jet, all_jets) {
double minDR_jet = MAXDOUBLE, minDR_electrons = MAXDOUBLE;
foreach( DressedLepton lep, leptons ) {
double dr = deltaR(jet, lep);
if ( dr < minDR_jet ) minDR_jet = dr;
if ( lep.abspid() == 11 && dr < minDR_electrons ) minDR_electrons = dr;
}
if ( minDR_electrons < 0.05 ) continue; // veto jet if it overlaps with electron
if ( minDR_jet < minDR_overall ) minDR_overall = minDR_jet;
jets += jet;
}
if ( jets.size() < 2 || minDR_overall <= 0.3 ) vetoEvent;
FourMomentum dijet = jets[0].momentum() + jets[1].momentum();
if ( dijet.mass() <= 500*GeV ) vetoEvent;
// inclusive region
- _hist->fill(0.5, weight);
+ _hist->fill(1.0, weight);
// VBS region
- if ( deltaRap(jets[0], jets[1]) > 2.4 ) _hist->fill(1.5, weight);
+ if ( deltaRap(jets[0], jets[1]) > 2.4 ) _hist->fill(2.0, weight);
}
/// Normalise histograms etc., after the run
void finalize() {
const double scalefactor( crossSection() / sumOfWeights() );
scale(_hist, scalefactor);
}
//@}
private:
/// @name Histograms
//@{
Histo1DPtr _hist;
//@}
};
// The hook for the plugin system
DECLARE_RIVET_PLUGIN(ATLAS_2014_I1298023);
}
diff --git a/analyses/pluginATLAS/ATLAS_2014_I1298023.plot b/analyses/pluginATLAS/ATLAS_2014_I1298023.plot
--- a/analyses/pluginATLAS/ATLAS_2014_I1298023.plot
+++ b/analyses/pluginATLAS/ATLAS_2014_I1298023.plot
@@ -1,7 +1,7 @@
# BEGIN PLOT /ATLAS_2014_I1298023/d.*
-XCustomMajorTicks=0.5 inclusive 1.5 VBS
+XCustomMajorTicks=1.0 inclusive 2.0 VBS
XLabel=phase space
YLabel=$\sigma^\text{fid}$ [pb]
LogY=0
# END PLOT
diff --git a/analyses/pluginATLAS/ATLAS_2014_I1298023.yoda b/analyses/pluginATLAS/ATLAS_2014_I1298023.yoda
--- a/analyses/pluginATLAS/ATLAS_2014_I1298023.yoda
+++ b/analyses/pluginATLAS/ATLAS_2014_I1298023.yoda
@@ -1,10 +1,10 @@
BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2014_I1298023/d01-x01-y01
IsRef: 1
Path: /REF/ATLAS_2014_I1298023/d01-x01-y01
-Title: ~
+Title: doi:10.17182/hepdata.64182.v1/t1
Type: Scatter2D
---
# xval xerr- xerr+ yval yerr- yerr+
-5.000000e-01 5.000000e-01 5.000000e-01 2.100000e+00 5.830952e-01 5.830952e-01
-1.500000e+00 5.000000e-01 5.000000e-01 1.300000e+00 4.472136e-01 4.472136e-01
+1.000000e+00 5.000000e-01 5.000000e-01 2.100000e+00 5.830952e-01 5.830952e-01
+2.000000e+00 5.000000e-01 5.000000e-01 1.300000e+00 4.472136e-01 4.472136e-01
END YODA_SCATTER2D_V2

File Metadata

Mime Type
text/x-diff
Expires
Tue, Nov 19, 4:29 PM (1 d, 13 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3805180
Default Alt Text
(5 KB)

Event Timeline