|
[Rivet-svn] r2885 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgSat Jan 22 16:46:23 GMT 2011
Author: fsiegert Date: Sat Jan 22 16:46:22 2011 New Revision: 2885 Log: Add ATLAS_2010_S8919674: W+jets measurement in both electron and muon channel with jet multiplicities and jet_1,2 pT. The aida output from Hepforge didn't work for this analysis, so the aida file was built via a modified yoda output. Might be worth dumping the original aida file again as soon as HepData works. The two multiplicity ratio plots are also not filled yet, but since they are just calculated from the existing plots I don't consider them crucial. Until I have run it using the MC samples used in the analysis I still leave it marked as UNVALIDATED. Added: trunk/data/anainfo/ATLAS_2010_S8919674.info trunk/data/plotinfo/ATLAS_2010_S8919674.plot trunk/data/refdata/ATLAS_2010_S8919674.aida trunk/src/Analyses/ATLAS_2010_S8919674.cc Modified: trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/src/Analyses/Makefile.am Added: trunk/data/anainfo/ATLAS_2010_S8919674.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/ATLAS_2010_S8919674.info Sat Jan 22 16:46:22 2011 (r2885) @@ -0,0 +1,41 @@ +Name: ATLAS_2010_S8919674 +Year: 2010 +Summary: W+jets jet multiplicities and pT +Experiment: ATLAS +Collider: LHC +SpiresID: 8919674 +Status: UNVALIDATED +Authors: + - Frank Siegert <frank.siegert at cern.ch> +References: + - arXiv:1012.5382 [hep-ex] +RunInfo: + W+jets events ideally with matrix element corrections to describe the + higher jet multiplicities correctly. Both channels, electron and muon, + are part of this analysis and should be run simultaneously. +NumEvents: 1000000 +Beams: [p+, p+] +Energies: [7000] +PtCuts: [0] +Description: + 'Cross sections, in both the electron and muon decay modes of the W boson, + are presented as a function of jet multiplicity and of the transverse + momentum of the leading and next-to-leading jets in the event. + Measurements are also presented of the ratio of cross sections for inclusive + jet multiplicities. + The results, based on an integrated luminosity of 1.3 pb-1, have been + corrected for all known detector effects and are quoted in a limited and + well-defined range of jet and lepton kinematics.' +BibKey: Aad:2010pg +BibTeX: '@Article{Aad:2010pg, + author = "Collaboration, The ATLAS", + title = "{Measurement of the production cross section for W-bosons + in association with jets in pp collisions at sqrt(s) = 7 + TeV with the ATLAS detector}", + year = "2010", + eprint = "1012.5382", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = 1012.5382;%%" +}' +ToDo: Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Fri Jan 21 18:55:51 2011 (r2884) +++ trunk/data/anainfo/Makefile.am Sat Jan 22 16:46:22 2011 (r2885) @@ -8,6 +8,7 @@ ALICE_2010_S8706239.info \ ATLAS_2010_S8591806.info \ ATLAS_2010_S8817804.info \ + ATLAS_2010_S8919674.info \ ATLAS_2010_CONF_2010_031.info \ ATLAS_2010_CONF_2010_046.info \ ATLAS_2010_CONF_2010_049.info \ Added: trunk/data/plotinfo/ATLAS_2010_S8919674.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/ATLAS_2010_S8919674.plot Sat Jan 22 16:46:22 2011 (r2885) @@ -0,0 +1,39 @@ +# BEGIN PLOT /ATLAS_2010_S8919674/d01-x01-y01 +Title=Inclusive jet multiplicity (electron channel) +XLabel=$N_{\text{jet}}$ +YLabel=$\sigma(W \, + \geq N_{\text{jet}}\;\text{jets})$ [pb] +XMajorTickMarks=10 +XMinorTickMarks=0 +# END PLOT + +# BEGIN PLOT /ATLAS_2010_S8919674/d02-x01-y01 +Title=Inclusive jet multiplicity (muon channel) +XLabel=$N_{\text{jet}}$ +YLabel=$\sigma(W \, + \geq N_{\text{jet}}\;\text{jets})$ [pb] +XMajorTickMarks=10 +XMinorTickMarks=0 +# END PLOT + +# BEGIN PLOT /ATLAS_2010_S8919674/d05-x01-y01 +Title=$p_\perp$ of 1st jet (electron channel) +XLabel=$p_\perp$(1st jet) [GeV] +YLabel=$\text{d}\sigma/\text{d}p_\perp$ [pb/GeV] +# END PLOT + +# BEGIN PLOT /ATLAS_2010_S8919674/d06-x01-y01 +Title=$p_\perp$ of 1st jet (muon channel) +XLabel=$p_\perp$(1st jet) [GeV] +YLabel=$\text{d}\sigma/\text{d}p_\perp$ [pb/GeV] +# END PLOT + +# BEGIN PLOT /ATLAS_2010_S8919674/d07-x01-y01 +Title=$p_\perp$ of 2nd jet (electron channel) +XLabel=$p_\perp$(2nd jet) [GeV] +YLabel=$\text{d}\sigma/\text{d}p_\perp$ [pb/GeV] +# END PLOT + +# BEGIN PLOT /ATLAS_2010_S8919674/d08-x01-y01 +Title=$p_\perp$ of 2nd jet (muon channel) +XLabel=$p_\perp$(2nd jet) [GeV] +YLabel=$\text{d}\sigma/\text{d}p_\perp$ [pb/GeV] +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Fri Jan 21 18:55:51 2011 (r2884) +++ trunk/data/plotinfo/Makefile.am Sat Jan 22 16:46:22 2011 (r2885) @@ -8,6 +8,7 @@ ALICE_2010_S8624100.plot \ ATLAS_2010_S8591806.plot \ ATLAS_2010_S8817804.plot \ + ATLAS_2010_S8919674.plot \ ATLAS_2010_CONF_2010_031.plot \ ATLAS_2010_CONF_2010_046.plot \ ATLAS_2010_CONF_2010_049.plot \ Added: trunk/data/refdata/ATLAS_2010_S8919674.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/ATLAS_2010_S8919674.aida Sat Jan 22 16:46:22 2011 (r2885) @@ -0,0 +1,185 @@ +<?xml version="1.0" encoding="ISO-8859-1" ?> +<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd"> +<aida version="3.3"> + <implementation version="1.1" package="FreeHEP"/> + <dataPointSet name="d01-x01-y01" dimension="2" + path="/REF/ATLAS_2010_S8919674" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2010_S8919674/d01-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2010_S8919674.aida/REF/ATLAS_2010_S8919674/d01-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="0.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="6.214947e+02" value="4.530000e+03" errorMinus="6.214947e+02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.553783e+02" value="8.400000e+02" errorMinus="1.553783e+02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="2.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="4.420338e+01" value="2.100000e+02" errorMinus="4.420338e+01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="3.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.596555e+01" value="4.700000e+01" errorMinus="1.596555e+01"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" + path="/REF/ATLAS_2010_S8919674" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2010_S8919674/d02-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2010_S8919674.aida/REF/ATLAS_2010_S8919674/d02-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="0.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="6.516000e+02" value="4.580000e+03" errorMinus="6.516000e+02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="1.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.445984e+02" value="8.400000e+02" errorMinus="1.445984e+02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="2.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="4.754135e+01" value="2.300000e+02" errorMinus="4.754135e+01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="3.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="1.924695e+01" value="6.400000e+01" errorMinus="1.924695e+01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e-01" value="4.000000e+00" errorMinus="5.000000e-01"/> + <measurement errorPlus="8.570782e+00" value="1.900000e+01" errorMinus="8.570782e+00"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d05-x01-y01" dimension="2" + path="/REF/ATLAS_2010_S8919674" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2010_S8919674/d05-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2010_S8919674.aida/REF/ATLAS_2010_S8919674/d05-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="2.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="7.243666e+00" value="3.760000e+01" errorMinus="7.243666e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="3.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="3.608877e+00" value="1.990000e+01" errorMinus="3.608877e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="4.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="1.653053e+00" value="8.000000e+00" errorMinus="1.653053e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="6.000000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="9.880853e-01" value="4.500000e+00" errorMinus="9.880853e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="8.000000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="4.943284e-01" value="2.000000e+00" errorMinus="4.943284e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.500000e+01" value="1.050000e+02" errorMinus="1.500000e+01"/> + <measurement errorPlus="3.000000e-01" value="8.000000e-01" errorMinus="3.000000e-01"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d06-x01-y01" dimension="2" + path="/REF/ATLAS_2010_S8919674" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2010_S8919674/d06-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2010_S8919674.aida/REF/ATLAS_2010_S8919674/d06-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="2.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="6.807126e+00" value="3.560000e+01" errorMinus="6.807126e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="3.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="3.293911e+00" value="1.880000e+01" errorMinus="3.293911e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="4.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="1.889936e+00" value="1.020000e+01" errorMinus="1.889936e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="6.000000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="9.880853e-01" value="4.800000e+00" errorMinus="9.880853e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="8.000000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="5.608058e-01" value="2.000000e+00" errorMinus="5.608058e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.500000e+01" value="1.050000e+02" errorMinus="1.500000e+01"/> + <measurement errorPlus="3.232051e-01" value="1.000000e+00" errorMinus="3.232051e-01"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d07-x01-y01" dimension="2" + path="/REF/ATLAS_2010_S8919674" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2010_S8919674/d07-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2010_S8919674.aida/REF/ATLAS_2010_S8919674/d07-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="2.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="2.990574e+00" value="1.232000e+01" errorMinus="2.990574e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="3.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="9.640812e-01" value="3.320000e+00" errorMinus="9.640812e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="4.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="6.484348e-01" value="2.170000e+00" errorMinus="6.484348e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="6.000000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="3.246890e-01" value="9.400000e-01" errorMinus="3.246890e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="8.000000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="1.705892e-01" value="3.300000e-01" errorMinus="1.705892e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.500000e+01" value="1.050000e+02" errorMinus="1.500000e+01"/> + <measurement errorPlus="7.483315e-02" value="1.100000e-01" errorMinus="7.483315e-02"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d08-x01-y01" dimension="2" + path="/REF/ATLAS_2010_S8919674" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2010_S8919674/d08-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2010_S8919674.aida/REF/ATLAS_2010_S8919674/d08-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="2.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="2.935891e+00" value="1.230000e+01" errorMinus="2.935891e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="3.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="1.010885e+00" value="4.210000e+00" errorMinus="1.010885e+00"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.000000e+00" value="4.500000e+01" errorMinus="5.000000e+00"/> + <measurement errorPlus="5.589804e-01" value="1.410000e+00" errorMinus="5.589804e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="6.000000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="4.354881e-01" value="1.520000e+00" errorMinus="4.354881e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="8.000000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="1.652246e-01" value="2.700000e-01" errorMinus="1.652246e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.500000e+01" value="1.050000e+02" errorMinus="1.500000e+01"/> + <measurement errorPlus="9.209308e-02" value="1.100000e-01" errorMinus="9.209308e-02"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Fri Jan 21 18:55:51 2011 (r2884) +++ trunk/data/refdata/Makefile.am Sat Jan 22 16:46:22 2011 (r2885) @@ -10,6 +10,7 @@ ALICE_2010_S8624100.aida \ ATLAS_2010_S8591806.aida \ ATLAS_2010_S8817804.aida \ + ATLAS_2010_S8919674.aida \ ATLAS_2010_CONF_2010_031.aida \ ATLAS_2010_CONF_2010_046.aida \ ATLAS_2010_CONF_2010_049.aida \ Added: trunk/src/Analyses/ATLAS_2010_S8919674.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/ATLAS_2010_S8919674.cc Sat Jan 22 16:46:22 2011 (r2885) @@ -0,0 +1,203 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/IdentifiedFinalState.hh" +#include "Rivet/Projections/MissingMomentum.hh" +#include "Rivet/Projections/FastJets.hh" +#include "Rivet/Projections/ClusteredPhotons.hh" +#include "Rivet/Projections/LeadingParticlesFinalState.hh" + +namespace Rivet { + + + class ATLAS_2010_S8919674 : public Analysis { + public: + + /// @name Constructors etc. + //@{ + + /// Constructor + ATLAS_2010_S8919674() + : Analysis("ATLAS_2010_S8919674") + { + setNeedsCrossSection(true); + } + + //@} + + + public: + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + /// Initialise and register projections (selections on the final state) + // projection to find the electrons + std::vector<std::pair<double, double> > eta_e; + eta_e.push_back(make_pair(-2.47,-1.52)); + eta_e.push_back(make_pair(-1.37,1.37)); + eta_e.push_back(make_pair(1.52,2.47)); + IdentifiedFinalState elecs(eta_e, 20.0*GeV); + elecs.acceptIdPair(ELECTRON); + addProjection(elecs, "elecs"); + // projection for finding the photons which have to be clustered into + // the lepton later + ClusteredPhotons cphotons_e(FinalState(), elecs, 0.1); + addProjection(cphotons_e, "cphotons_e"); + + // 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, 20.0*GeV); + muons.acceptIdPair(MUON); + addProjection(muons, "muons"); + // projection for finding the photons which have to be clustered into + // the lepton later + ClusteredPhotons cphotons_m(FinalState(), muons, 0.1); + addProjection(cphotons_m, "cphotons_m"); + + // Leading neutrinos for Etmiss + FinalState fs; + LeadingParticlesFinalState muon_neutrino(fs); + muon_neutrino.addParticleIdPair(NU_MU); + addProjection(muon_neutrino, "muon_neutrino"); + LeadingParticlesFinalState elec_neutrino(fs); + elec_neutrino.addParticleIdPair(NU_E); + addProjection(elec_neutrino, "elec_neutrino"); + + // Input for the jets: No neutrinos, no muons, and no electron which + // passed the electron cuts ("elecs" finalstate from above) + VetoedFinalState veto; + veto.addVetoOnThisFinalState(elecs); + veto.addVetoPairId(MUON); + veto.vetoNeutrinos(); + FastJets jets(veto, FastJets::ANTIKT, 0.4); + addProjection(jets, "jets"); + + + + /// book histograms + _h_el_njet_inclusive = bookHistogram1D(1,1,1); + _h_mu_njet_inclusive = bookHistogram1D(2,1,1); + + _h_el_pT_jet1 = bookHistogram1D(5,1,1); + _h_mu_pT_jet1 = bookHistogram1D(6,1,1); + + _h_el_pT_jet2 = bookHistogram1D(7,1,1); + _h_mu_pT_jet2 = bookHistogram1D(8,1,1); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const double weight = event.weight(); + + const FinalState& elecs = applyProjection<FinalState>(event, "elecs"); + ParticleVector elec_neutrino=applyProjection<FinalState>(event, "elec_neutrino").particles(); + if (elecs.size()==1 && elec_neutrino.size()>0) { + FourMomentum lepton=elecs.particles()[0].momentum(); + foreach (const Particle& photon, + applyProjection<FinalState>(event, "cphotons_e").particles()) { + lepton+=photon.momentum(); + } + FourMomentum p_miss = elec_neutrino[0].momentum(); + double mT=sqrt(2.0*lepton.pT()*p_miss.Et()*(1.0-cos(lepton.phi()-p_miss.phi()))); + if (p_miss.Et()>25.0*GeV && mT>40.0*GeV) { + Jets jets; + foreach (const Jet& jet, applyProjection<FastJets>(event, "jets").jetsByPt(20.0*GeV)) { + if (fabs(jet.eta())<2.8 && deltaR(lepton, jet.momentum())>0.5) { + jets.push_back(jet); + } + } + + _h_el_njet_inclusive->fill(0, weight); + if (jets.size()>=1) { + _h_el_njet_inclusive->fill(1, weight); + _h_el_pT_jet1->fill(jets[0].momentum().pT(), weight); + } + if (jets.size()>=2) { + _h_el_njet_inclusive->fill(2, weight); + _h_el_pT_jet2->fill(jets[1].momentum().pT(), weight); + } + if (jets.size()>=3) { + _h_el_njet_inclusive->fill(3, weight); + } + } + } + + const FinalState& muons = applyProjection<FinalState>(event, "muons"); + ParticleVector muon_neutrino=applyProjection<FinalState>(event, "muon_neutrino").particles(); + if (muons.size()==1 && muon_neutrino.size()>0) { + FourMomentum lepton=muons.particles()[0].momentum(); + foreach (const Particle& photon, + applyProjection<FinalState>(event, "cphotons_m").particles()) { + lepton+=photon.momentum(); + } + FourMomentum p_miss = muon_neutrino[0].momentum(); + double mT=sqrt(2.0*lepton.pT()*p_miss.Et()*(1.0-cos(lepton.phi()-p_miss.phi()))); + if (p_miss.Et()>25.0*GeV && mT>40.0*GeV) { + Jets jets; + foreach (const Jet& jet, applyProjection<FastJets>(event, "jets").jetsByPt(20.0*GeV)) { + if (fabs(jet.eta())<2.8 && deltaR(lepton, jet.momentum())>0.5) { + jets.push_back(jet); + } + } + + _h_mu_njet_inclusive->fill(0, weight); + if (jets.size()>=1) { + _h_mu_njet_inclusive->fill(1, weight); + _h_mu_pT_jet1->fill(jets[0].momentum().pT(), weight); + } + if (jets.size()>=2) { + _h_mu_njet_inclusive->fill(2, weight); + _h_mu_pT_jet2->fill(jets[1].momentum().pT(), weight); + } + if (jets.size()>=3) { + _h_mu_njet_inclusive->fill(3, weight); + } + } + } + + } + + + /// Normalise histograms etc., after the run + void finalize() { + double normfac=crossSection()/sumOfWeights(); + scale(_h_el_njet_inclusive, normfac); + scale(_h_mu_njet_inclusive, normfac); + scale(_h_el_pT_jet1, normfac); + scale(_h_mu_pT_jet1, normfac); + scale(_h_el_pT_jet2, normfac); + scale(_h_mu_pT_jet2, normfac); + } + + //@} + + + private: + + /// @name Histograms + //@{ + + AIDA::IHistogram1D * _h_el_njet_inclusive; + AIDA::IHistogram1D * _h_mu_njet_inclusive; + AIDA::IHistogram1D * _h_el_pT_jet1; + AIDA::IHistogram1D * _h_mu_pT_jet1; + AIDA::IHistogram1D * _h_el_pT_jet2; + AIDA::IHistogram1D * _h_mu_pT_jet2; + //@} + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<ATLAS_2010_S8919674> plugin_ATLAS_2010_S8919674; + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Fri Jan 21 18:55:51 2011 (r2884) +++ trunk/src/Analyses/Makefile.am Sat Jan 22 16:46:22 2011 (r2885) @@ -46,7 +46,8 @@ RivetALICEAnalyses_la_SOURCES = \ ALICE_2010_S8706239.cc \ ALICE_2010_S8625980.cc \ - ALICE_2010_S8624100.cc + ALICE_2010_S8624100.cc \ + ATLAS_2010_S8919674.cc lib_LTLIBRARIES += RivetATLASAnalyses.la
More information about the Rivet-svn mailing list |