Page MenuHomeHEPForge
Files F8723523

ATLAS_2012_CONF_2012_103.cc
No OneTemporary
Actions

ATLAS_2012_CONF_2012_103.cc
View Options

// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Tools/BinnedHistogram.hh"
#include "Rivet/RivetAIDA.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/VisibleFinalState.hh"
#include "Rivet/Projections/VetoedFinalState.hh"
#include "Rivet/Projections/IdentifiedFinalState.hh"
#include "Rivet/Projections/FastJets.hh"
#include "Rivet/Tools/RivetMT2.hh"
namespace Rivet {
class ATLAS_2012_CONF_2012_103 : public Analysis {
public:
/// @name Constructors etc.
//@{
/// Constructor
ATLAS_2012_CONF_2012_103()
: Analysis("ATLAS_2012_CONF_2012_103")
{ }
//@}
public:
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
// projection to find the electrons
std::vector<std::pair<double, double> > eta_e;
eta_e.push_back(make_pair(-2.47,2.47));
IdentifiedFinalState elecs(eta_e, 20.0*GeV);
elecs.acceptIdPair(ELECTRON);
addProjection(elecs, "elecs");
// projection to find the muons
std::vector<std::pair<double, double> > eta_m;
eta_m.push_back(make_pair(-2.4,2.4));
IdentifiedFinalState muons(eta_m, 10.0*GeV);
muons.acceptIdPair(MUON);
addProjection(muons, "muons");
// for pTmiss
addProjection(VisibleFinalState(-4.9,4.9),"vfs");
VetoedFinalState vfs;
vfs.addVetoPairId(MUON);
/// Jet finder
addProjection(FastJets(vfs, FastJets::ANTIKT, 0.4),
"AntiKtJets04");
/// Book histograms
_etmiss_HT_7j55 = bookHistogram1D("etmiss_HT_7j55", 8, 0., 16.);
_etmiss_HT_8j55 = bookHistogram1D("etmiss_HT_8j55", 8, 0., 16.);
_etmiss_HT_9j55 = bookHistogram1D("etmiss_HT_9j55", 8, 0., 16.);
_etmiss_HT_6j80 = bookHistogram1D("etmiss_HT_6j80", 8, 0., 16.);
_etmiss_HT_7j80 = bookHistogram1D("etmiss_HT_7j80", 8, 0., 16.);
_etmiss_HT_8j80 = bookHistogram1D("etmiss_HT_8j80", 8, 0., 16.);
_hist_njet55 = bookHistogram1D("hist_njet55", 4, 5.5, 9.5);
_hist_njet80 = bookHistogram1D("hist_njet80", 4, 4.5, 8.5);
_count_7j55 = bookHistogram1D("count_7j55", 1, 0., 1.);
_count_8j55 = bookHistogram1D("count_8j55", 1, 0., 1.);
_count_9j55 = bookHistogram1D("count_9j55", 1, 0., 1.);
_count_6j80 = bookHistogram1D("count_6j80", 1, 0., 1.);
_count_7j80 = bookHistogram1D("count_7j80", 1, 0., 1.);
_count_8j80 = bookHistogram1D("count_8j80", 1, 0., 1.);
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const double weight = event.weight();
// get the jet candidates
Jets cand_jets;
foreach (const Jet& jet,
applyProjection<FastJets>(event, "AntiKtJets04").jetsByPt(20.0*GeV) ) {
if ( fabs( jet.momentum().eta() ) < 2.8 ) {
cand_jets.push_back(jet);
}
}
// candidate muons
ParticleVector cand_mu =
applyProjection<IdentifiedFinalState>(event, "muons").particlesByPt();
// candidate electrons
ParticleVector cand_e =
applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt();
// resolve jet/lepton ambiguity
Jets recon_jets;
foreach ( const Jet& jet, cand_jets ) {
// candidates after |eta| < 2.8
if ( fabs( jet.momentum().eta() ) >= 2.8 ) continue;
bool away_from_e = true;
foreach ( const Particle & e, cand_e ) {
if ( deltaR(e.momentum(),jet.momentum()) <= 0.2 ) {
away_from_e = false;
break;
}
}
if ( away_from_e ) recon_jets.push_back( jet );
}
// only keep electrons more than R=0.4 from jets
ParticleVector recon_e;
foreach ( const Particle & e, cand_e ) {
bool away = true;
foreach ( const Jet& jet, recon_jets ) {
if ( deltaR(e.momentum(),jet.momentum()) < 0.4 ) {
away = false;
break;
}
}
if ( away )
recon_e.push_back( e );
}
// only keep muons more than R=0.4 from jets
ParticleVector recon_mu;
foreach ( const Particle & mu, cand_mu ) {
bool away = true;
foreach ( const Jet& jet, recon_jets ) {
if ( deltaR(mu.momentum(),jet.momentum()) < 0.4 ) {
away = false;
break;
}
}
if ( away )
recon_mu.push_back( mu );
}
// pTmiss
ParticleVector vfs_particles =
applyProjection<VisibleFinalState>(event, "vfs").particles();
FourMomentum pTmiss;
foreach ( const Particle & p, vfs_particles ) {
pTmiss -= p.momentum();
}
double eTmiss = pTmiss.pT();
// now only use recon_jets, recon_mu, recon_e
// reject events with electrons and muons
if ( ! ( recon_mu.empty() && recon_e.empty() ) ) {
MSG_DEBUG("Charged leptons left after selection");
vetoEvent;
}
// calculate H_T
double HT=0;
foreach ( const Jet& jet, recon_jets ) {
if ( jet.momentum().pT() > 40 * GeV )
HT += jet.momentum().pT() ;
}
// number of jets
unsigned int njet55=0, njet80=0;
for (unsigned int ix=0;ix<recon_jets.size();++ix) {
if(recon_jets[ix].momentum().pT()>80.*GeV) ++njet80;
if(recon_jets[ix].momentum().pT()>55.*GeV) ++njet55;
}
double ratio = eTmiss/sqrt(HT);
if(ratio>4.) {
if(njet55>9) njet55 = 9;
if(njet80>8) njet80 = 8;
_hist_njet55->fill(njet55,weight);
_hist_njet80->fill(njet80,weight);
// 7j55
if(njet55>=7)
_count_7j55->fill( 0.5, weight);
// 8j55
if(njet55>=8)
_count_8j55->fill( 0.5, weight) ;
// 8j55
if(njet55==9)
_count_9j55->fill( 0.5, weight) ;
// 6j80
if(njet80>=6)
_count_6j80->fill( 0.5, weight) ;
// 7j80
if(njet80>=7)
_count_7j80->fill( 0.5, weight) ;
// 8j80
if(njet80==8)
_count_8j80->fill( 0.5, weight) ;
}
if(njet55>=7)
_etmiss_HT_7j55->fill( ratio, weight);
// 8j55
if(njet55>=8)
_etmiss_HT_8j55->fill( ratio, weight) ;
// 8j55
if(njet55>=9)
_etmiss_HT_9j55->fill( ratio, weight) ;
// 6j80
if(njet80>=6)
_etmiss_HT_6j80->fill( ratio, weight) ;
// 7j80
if(njet80>=7)
_etmiss_HT_7j80->fill( ratio, weight) ;
// 8j80
if(njet80>=8)
_etmiss_HT_8j80->fill( ratio, weight) ;
}
//@}
void finalize() {
double norm = crossSection()/femtobarn*5.8/sumOfWeights();
scale(_etmiss_HT_7j55,2.*norm);
scale(_etmiss_HT_8j55,2.*norm);
scale(_etmiss_HT_9j55,2.*norm);
scale(_etmiss_HT_6j80,2.*norm);
scale(_etmiss_HT_7j80,2.*norm);
scale(_etmiss_HT_8j80,2.*norm);
scale(_hist_njet55,norm);
scale(_hist_njet80,norm);
scale(_count_7j55,norm);
scale(_count_8j55,norm);
scale(_count_9j55,norm);
scale(_count_6j80,norm);
scale(_count_7j80,norm);
scale(_count_8j80,norm);
}
private:
/// @name Histograms
//@{
AIDA::IHistogram1D* _etmiss_HT_7j55;
AIDA::IHistogram1D* _etmiss_HT_8j55;
AIDA::IHistogram1D* _etmiss_HT_9j55;
AIDA::IHistogram1D* _etmiss_HT_6j80;
AIDA::IHistogram1D* _etmiss_HT_7j80;
AIDA::IHistogram1D* _etmiss_HT_8j80;
AIDA::IHistogram1D* _hist_njet55;
AIDA::IHistogram1D* _hist_njet80;
AIDA::IHistogram1D* _count_7j55;
AIDA::IHistogram1D* _count_8j55;
AIDA::IHistogram1D* _count_9j55;
AIDA::IHistogram1D* _count_6j80;
AIDA::IHistogram1D* _count_7j80;
AIDA::IHistogram1D* _count_8j80;
//@}
};
// The hook for the plugin system
DECLARE_RIVET_PLUGIN(ATLAS_2012_CONF_2012_103);
}

File Metadata

Mime Type
text/x-c
Expires
Mon, Jan 20, 8:39 PM (12 h, 15 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4242301
Default Alt Text
ATLAS_2012_CONF_2012_103.cc (7 KB)

Event Timeline