|
[Rivet-svn] r3476 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgWed Nov 9 16:24:46 GMT 2011
Author: buckley Date: Wed Nov 9 16:24:45 2011 New Revision: 3476 Log: Adding ATLAS W pT analysis -- build currently broken because of a missing AIDA file Added: trunk/data/anainfo/ATLAS_2011_I925932.info trunk/data/plotinfo/ATLAS_2011_I925932.plot trunk/data/refdata/ATLAS_2011_I925932.aida trunk/src/Analyses/ATLAS_2011_I925932.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_2011_I925932.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/ATLAS_2011_I925932.info Wed Nov 9 16:24:45 2011 (r3476) @@ -0,0 +1,35 @@ +Name: ATLAS_2011_I925932 +Year: 2011 +Summary: Measurement of the W pT with electrons and muons at 7 TeV +Experiment: ATLAS +Collider: LHC +InspireID: 925932 +Status: VALIDATED +Authors: + - Elena Yatsenko <elena.yatsenko at desy.de> + - Judith Katzy <jkatzy at mail.cern.ch> +References: + - arXiv:1108.6308v1 [hep-ex] +RunInfo: + Run with $W$ decays to $e\nu_e$ and/or $\mu\nu_{\mu}$. +NumEvents: 100k +Beams: [p+, p+] +Energies: [7000] +PtCuts: [0] +Description: + 'The W pT at $\sqrt{s} = 7$\;TeV is measured using $W\to e\nu_e$ and $W\to \mu\nu_{\mu}$ + decay channels. The dressed leptons kinematics calculated from the sum of the post-FSR + lepton momentum and the momenta of all photons radiated in a cone around the lepton, + while the bare uses the lepton kinematics after all QED radiation.' +BibKey: +BibTeX: '@article{Aad:2011fp, + author = "Aad, Georges and others", + title = "Measurement of the transverse momentum distribution of $W$ + bosons in $pp$ Collisions at $\sqrt{s} = 7$~TeV with the {ATLAS + Detector}", + collaboration = "The ATLAS", + year = "2011", + eprint = "1108.6308", + archivePrefix = "arXiv", + primaryClass = "hep-ex", +}' Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Wed Nov 9 13:46:55 2011 (r3475) +++ trunk/data/anainfo/Makefile.am Wed Nov 9 16:24:45 2011 (r3476) @@ -28,9 +28,10 @@ ATLAS_2011_S9126244.info \ ATLAS_2011_S9128077.info \ ATLAS_2011_S9131140.info \ - ATLAS_2011_I919017.info \ ATLAS_2011_S9108483.info \ ATLAS_2011_S9212183.info \ + ATLAS_2011_I919017.info \ + ATLAS_2011_I925932.info \ BELLE_2006_S6265367.info \ CDF_1988_S1865951.info \ CDF_1990_S2089246.info \ Added: trunk/data/plotinfo/ATLAS_2011_I925932.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/ATLAS_2011_I925932.plot Wed Nov 9 16:24:45 2011 (r3476) @@ -0,0 +1,27 @@ +# BEGIN PLOT /ATLAS_2011_I925932/d01-x01-y01 +Title=$W\to e\nu_e$ $p_\perp$ with "dressed" kinematics +XLabel=$p_\perp^W$ [GeV] +YLabel=$(1/\sigma_{fid})(\mathrm{d}\sigma_{fid}/\mathrm{d}p_\perp^W)$ [GeV$^{-1}$] +YMin=0.000001 +# END PLOT + +# BEGIN PLOT /ATLAS_2011_I925932/d01-x01-y02 +Title=$W\to e\nu_e$ $p_\perp$ with "bare" kinematics +XLabel=$p_\perp^W$ [GeV] +YLabel=$(1/\sigma_{fid})(\mathrm{d}\sigma_{fid}/\mathrm{d}p_\perp^W)$ [GeV$^{-1}$] +YMin=0.000001 +# END PLOT + +# BEGIN PLOT /ATLAS_2011_I925932/d02-x01-y01 +Title=$W\to \mu\nu_{\mu}$ $p_\perp$ with "dressed" kinematics +XLabel=$p_\perp^W$ [GeV] +YLabel=$(1/\sigma_{fid})(\mathrm{d}\sigma_{fid}/\mathrm{d}p_\perp^W)$ [GeV$^{-1}$] +YMin=0.000001 +# END PLOT + +# BEGIN PLOT /ATLAS_2011_I925932/d02-x01-y02 +Title=$W\to \mu\nu_{\mu}$ $p_\perp$ with "bare" kinematics +XLabel=$p_\perp^W$ [GeV] +YLabel=$(1/\sigma_{fid})(\mathrm{d}\sigma_{fid}/\mathrm{d}p_\perp^W)$ [GeV$^{-1}$] +YMin=0.000001 +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Wed Nov 9 13:46:55 2011 (r3475) +++ trunk/data/plotinfo/Makefile.am Wed Nov 9 16:24:45 2011 (r3476) @@ -23,7 +23,9 @@ ATLAS_2011_S9126244.plot \ ATLAS_2011_S9128077.plot \ ATLAS_2011_S9131140.plot \ - ATLAS_2011_I919017.plot \ + ATLAS_2011_S9212183.plot \ + ATLAS_2011_I919017.plot \ + ATLAS_2011_I925932.plot \ CDF_1988_S1865951.plot \ CDF_1990_S2089246.plot \ CDF_1993_S2742446.plot \ Added: trunk/data/refdata/ATLAS_2011_I925932.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/ATLAS_2011_I925932.aida Wed Nov 9 16:24:45 2011 (r3476) @@ -0,0 +1,217 @@ +<?xml version="1.0" ?> +<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd"> +<aida version="3.3"> + <implementation version="1.1" package="Rivet"/> + <dataPointSet name="d01-x01-y01" dimension="2" + path="/REF/ATLAS_2011_I925932" title="$W\to e\nu_e$ $p_\perp$ with Dressed kinematics"> + <dimenstion dim="0" title="$p_\perp(W)$ [GeV]"/> + <dimenstion dim="1" title="$1/\sigma_{fid}\mathrm{d}\sigma_{fid}/\mathrm{d}p_\perp$ [GeV$^{-1}$]"/> + <annotation> + <item key="Title" value="$W\to e\nu_e$ $p_\perp$ with Dressed kinematics" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2011_I925932/d01-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="4.000000e+00" value="4.000000e+00" errorMinus="4.000000e+00"/> + <measurement errorPlus="1.579120e-03" value="5.545800e-02" errorMinus="1.579120e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.500000e+00" value="1.550000e+01" errorMinus="7.500000e+00"/> + <measurement errorPlus="7.397108e-04" value="2.515300e-02" errorMinus="7.397108e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.500000e+00" value="3.050000e+01" errorMinus="7.500000e+00"/> + <measurement errorPlus="3.259063e-04" value="6.757800e-03" errorMinus="3.259063e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="8.500000e+00" value="4.650000e+01" errorMinus="8.500000e+00"/> + <measurement errorPlus="1.249200e-04" value="2.461600e-03" errorMinus="1.249200e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="6.500000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="7.206898e-05" value="9.349100e-04" errorMinus="7.206898e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="8.500000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="3.739402e-05" value="3.729700e-04" errorMinus="3.739402e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.250000e+01" value="1.075000e+02" errorMinus="1.250000e+01"/> + <measurement errorPlus="2.070585e-05" value="1.802800e-04" errorMinus="2.070585e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.250000e+01" value="1.325000e+02" errorMinus="1.250000e+01"/> + <measurement errorPlus="1.138824e-05" value="9.488200e-05" errorMinus="1.138824e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.500000e+01" value="1.600000e+02" errorMinus="1.500000e+01"/> + <measurement errorPlus="4.538466e-06" value="3.541300e-05" errorMinus="4.538466e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.750000e+01" value="1.925000e+02" errorMinus="1.750000e+01"/> + <measurement errorPlus="2.325012e-06" value="1.583500e-05" errorMinus="2.325012e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.500000e+01" value="2.550000e+02" errorMinus="4.500000e+01"/> + <measurement errorPlus="9.084369e-07" value="4.669900e-06" errorMinus="9.084369e-07"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y02" dimension="2" + path="/REF/ATLAS_2011_I925932" title="$W\to e\nu_e$ $p_\perp$ with Bare kinematics"> + <dimenstion dim="0" title="$p_\perp(W)$ [GeV]"/> + <dimenstion dim="1" title="$1/\sigma_{fid}\mathrm{d}\sigma_{fid}/\mathrm{d}p_\perp$ [GeV$^{-1}$]"/> + <annotation> + <item key="Title" value="$W\to e\nu_e$ $p_\perp$ with Bare kinematics" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2011_I925932/d01-x01-y02" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="4.000000e+00" value="4.000000e+00" errorMinus="4.000000e+00"/> + <measurement errorPlus="1.542774e-03" value="5.418400e-02" errorMinus="1.542774e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.500000e+00" value="1.550000e+01" errorMinus="7.500000e+00"/> + <measurement errorPlus="7.528459e-04" value="2.559500e-02" errorMinus="7.528459e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.500000e+00" value="3.050000e+01" errorMinus="7.500000e+00"/> + <measurement errorPlus="3.354188e-04" value="6.956100e-03" errorMinus="3.354188e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="8.500000e+00" value="4.650000e+01" errorMinus="8.500000e+00"/> + <measurement errorPlus="1.248247e-04" value="2.460400e-03" errorMinus="1.248247e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="6.500000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="7.087936e-05" value="9.194400e-04" errorMinus="7.087936e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="8.500000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="3.651806e-05" value="3.642600e-04" errorMinus="3.651806e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.250000e+01" value="1.075000e+02" errorMinus="1.250000e+01"/> + <measurement errorPlus="2.013000e-05" value="1.751900e-04" errorMinus="2.013000e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.250000e+01" value="1.325000e+02" errorMinus="1.250000e+01"/> + <measurement errorPlus="1.103404e-05" value="9.192900e-05" errorMinus="1.103404e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.500000e+01" value="1.600000e+02" errorMinus="1.500000e+01"/> + <measurement errorPlus="4.396985e-06" value="3.431000e-05" errorMinus="4.396985e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.750000e+01" value="1.925000e+02" errorMinus="1.750000e+01"/> + <measurement errorPlus="2.237396e-06" value="1.523800e-05" errorMinus="2.237396e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.500000e+01" value="2.550000e+02" errorMinus="4.500000e+01"/> + <measurement errorPlus="8.733905e-07" value="4.489500e-06" errorMinus="8.733905e-07"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" + path="/REF/ATLAS_2011_I925932" title="$W\to \mu\nu_{\mu}$ $p_\perp$ with Dressed kinematics"> + <dimenstion dim="0" title="$p_\perp(W)$ [GeV]"/> + <dimenstion dim="1" title="$1/\sigma_{fid}\mathrm{d}\sigma_{fid}/\mathrm{d}p_\perp$ [GeV$^{-1}$]"/> + <annotation> + <item key="Title" value="$W\to \mu\nu_{\mu}$ $p_\perp$ with Dressed kinematics" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2011_I925932/d02-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="4.000000e+00" value="4.000000e+00" errorMinus="4.000000e+00"/> + <measurement errorPlus="1.439678e-03" value="5.386700e-02" errorMinus="1.439678e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.500000e+00" value="1.550000e+01" errorMinus="7.500000e+00"/> + <measurement errorPlus="6.549695e-04" value="2.538800e-02" errorMinus="6.549695e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.500000e+00" value="3.050000e+01" errorMinus="7.500000e+00"/> + <measurement errorPlus="3.380087e-04" value="7.061900e-03" errorMinus="3.380087e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="8.500000e+00" value="4.650000e+01" errorMinus="8.500000e+00"/> + <measurement errorPlus="1.064708e-04" value="2.554000e-03" errorMinus="1.064708e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="6.500000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="4.472192e-05" value="1.035300e-03" errorMinus="4.472192e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="8.500000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="2.244548e-05" value="4.373000e-04" errorMinus="2.244548e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.250000e+01" value="1.075000e+02" errorMinus="1.250000e+01"/> + <measurement errorPlus="1.252288e-05" value="1.901300e-04" errorMinus="1.252288e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.250000e+01" value="1.325000e+02" errorMinus="1.250000e+01"/> + <measurement errorPlus="6.730242e-06" value="7.291000e-05" errorMinus="6.730242e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.500000e+01" value="1.600000e+02" errorMinus="1.500000e+01"/> + <measurement errorPlus="4.548983e-06" value="3.955700e-05" errorMinus="4.548983e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.750000e+01" value="1.925000e+02" errorMinus="1.750000e+01"/> + <measurement errorPlus="3.260211e-06" value="1.863300e-05" errorMinus="3.260211e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.500000e+01" value="2.550000e+02" errorMinus="4.500000e+01"/> + <measurement errorPlus="1.007665e-06" value="4.658200e-06" errorMinus="1.007665e-06"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y02" dimension="2" + path="/REF/ATLAS_2011_I925932" title="$W\to \mu\nu_{\mu}$ $p_\perp$ with Bare kinematics"> + <dimenstion dim="0" title="$p_\perp(W)$ [GeV]"/> + <dimenstion dim="1" title="$1/\sigma_{fid}\mathrm{d}\sigma_{fid}/\mathrm{d}p_\perp$ [GeV$^{-1}$]"/> + <annotation> + <item key="Title" value="$W\to \mu\nu_{\mu}$ $p_\perp$ with Bare kinematics" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2011_I925932/d02-x01-y02" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="4.000000e+00" value="4.000000e+00" errorMinus="4.000000e+00"/> + <measurement errorPlus="1.429632e-03" value="5.350700e-02" errorMinus="1.429632e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.500000e+00" value="1.550000e+01" errorMinus="7.500000e+00"/> + <measurement errorPlus="6.590554e-04" value="2.554300e-02" errorMinus="6.590554e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.500000e+00" value="3.050000e+01" errorMinus="7.500000e+00"/> + <measurement errorPlus="3.403427e-04" value="7.109900e-03" errorMinus="3.403427e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="8.500000e+00" value="4.650000e+01" errorMinus="8.500000e+00"/> + <measurement errorPlus="1.063429e-04" value="2.549500e-03" errorMinus="1.063429e-04"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="6.500000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="4.449831e-05" value="1.029000e-03" errorMinus="4.449831e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e+01" value="8.500000e+01" errorMinus="1.000000e+01"/> + <measurement errorPlus="2.225947e-05" value="4.336500e-04" errorMinus="2.225947e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.250000e+01" value="1.075000e+02" errorMinus="1.250000e+01"/> + <measurement errorPlus="1.238083e-05" value="1.880700e-04" errorMinus="1.238083e-05"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.250000e+01" value="1.325000e+02" errorMinus="1.250000e+01"/> + <measurement errorPlus="6.655996e-06" value="7.210700e-05" errorMinus="6.655996e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.500000e+01" value="1.600000e+02" errorMinus="1.500000e+01"/> + <measurement errorPlus="4.491681e-06" value="3.906100e-05" errorMinus="4.491681e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.750000e+01" value="1.925000e+02" errorMinus="1.750000e+01"/> + <measurement errorPlus="3.212788e-06" value="1.836300e-05" errorMinus="3.212788e-06"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.500000e+01" value="2.550000e+02" errorMinus="4.500000e+01"/> + <measurement errorPlus="9.844233e-07" value="4.550700e-06" errorMinus="9.844233e-07"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Wed Nov 9 13:46:55 2011 (r3475) +++ trunk/data/refdata/Makefile.am Wed Nov 9 16:24:45 2011 (r3476) @@ -27,7 +27,9 @@ ATLAS_2011_S9126244.aida \ ATLAS_2011_S9128077.aida \ ATLAS_2011_S9131140.aida \ + ATLAS_2011_S9212183.aida \ ATLAS_2011_I919017.aida \ + ATLAS_2011_I925932.aida \ BABAR_2006_S6511112.aida \ BELLE_2006_S6265367.aida \ CLEO_1997_S3486664.aida \ Added: trunk/src/Analyses/ATLAS_2011_I925932.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/ATLAS_2011_I925932.cc Wed Nov 9 16:24:45 2011 (r3476) @@ -0,0 +1,135 @@ +// -*- C++ -*- +// ATLAS W pT analysis + +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/WFinder.hh" + +namespace Rivet { + + + class ATLAS_2011_I925932 : public Analysis { + public: + + /// Constructor + ATLAS_2011_I925932() : Analysis("ATLAS_2011_I925932") { } + + + /// @name Analysis methods + //@{ + + void init() { + // Set up projections + WFinder wfinder_dressed_el(-2.4, 2.4, 20.0*GeV, ELECTRON, 0.0*GeV, 1000.0*GeV, 25.0*GeV, 0.2,true,false); + addProjection(wfinder_dressed_el, "WFinder_dressed_el"); + WFinder wfinder_bare_el(-2.4, 2.4, 20.0*GeV, ELECTRON, 0.0*GeV, 1000.0*GeV, 25.0*GeV, 0.0,true,false); + addProjection(wfinder_bare_el, "WFinder_bare_el"); + WFinder wfinder_dressed_mu(-2.4, 2.4, 20.0*GeV, MUON, 0.0*GeV, 1000.0*GeV, 25.0*GeV, 0.2,true,false); + addProjection(wfinder_dressed_mu, "WFinder_dressed_mu"); + WFinder wfinder_bare_mu(-2.4, 2.4, 20.0*GeV, MUON, 0.0*GeV, 1000.0*GeV, 25.0*GeV, 0.0,true,false); + addProjection(wfinder_bare_mu, "WFinder_bare_mu"); + + // Counters + _sumW_e_dressed = 0; + _sumW_e_bare = 0; + _sumW_mu_dressed = 0; + _sumW_mu_bare = 0; + + // Book histograms + _hist_wpt_dressed_el = bookHistogram1D(1, 1, 1); + _hist_wpt_bare_el = bookHistogram1D(1, 1, 2); + _hist_wpt_dressed_mu = bookHistogram1D(2, 1, 1); + _hist_wpt_bare_mu = bookHistogram1D(2, 1, 2); + } + + + inline double calcMT(const FourMomentum& a, const FourMomentum& b) { + return sqrt(2.0 * a.pT() * b.pT() * (1.0 - cos(a.phi() - b.phi())) ); + } + + + /// Do the analysis + void analyze(const Event& event) { + + const WFinder& wfinder_dressed_el = applyProjection<WFinder>(event, "WFinder_dressed_el"); + const WFinder& wfinder_bare_el = applyProjection<WFinder>(event, "WFinder_bare_el"); + const WFinder& wfinder_dressed_mu = applyProjection<WFinder>(event, "WFinder_dressed_mu"); + const WFinder& wfinder_bare_mu = applyProjection<WFinder>(event, "WFinder_bare_mu"); + + if (wfinder_dressed_el.particles().empty() && wfinder_bare_el.particles().empty() && + wfinder_dressed_mu.particles().empty() && wfinder_bare_mu.particles().empty()) { + MSG_DEBUG("No W bosons found"); + vetoEvent; + } + + // "Dressed" electron + if (!wfinder_dressed_el.particles().empty()) { + const FourMomentum el = wfinder_dressed_el.constituentLeptons()[0].momentum(); + const FourMomentum nu = wfinder_dressed_el.constituentNeutrinos()[0].momentum(); + if (calcMT(el, nu) > 40.0*GeV && nu.pT() > 25.0*GeV) { + _sumW_e_dressed += event.weight(); + _hist_wpt_dressed_el->fill(wfinder_dressed_el.bosons()[0].momentum().pT()/GeV, event.weight()); + } + } + + // "Bare" electron + if (!wfinder_bare_el.particles().empty()) { + const FourMomentum el = wfinder_bare_el.constituentLeptons()[0].momentum(); + const FourMomentum nu = wfinder_bare_el.constituentNeutrinos()[0].momentum(); + if (calcMT(el, nu) > 40.0*GeV && nu.pT() > 25.0*GeV) { + _sumW_e_bare += event.weight(); + _hist_wpt_bare_el->fill(wfinder_bare_el.bosons()[0].momentum().pT()/GeV, event.weight()); + } + } + + // "Dressed" muon + if (!wfinder_dressed_mu.particles().empty()) { + const FourMomentum mu = wfinder_dressed_mu.constituentLeptons()[0].momentum(); + const FourMomentum nu = wfinder_dressed_mu.constituentNeutrinos()[0].momentum(); + if (calcMT(mu, nu) > 40.0*GeV && nu.pT() > 25.0*GeV) { + _sumW_mu_dressed += event.weight(); + _hist_wpt_dressed_mu->fill(wfinder_dressed_mu.bosons()[0].momentum().pT()/GeV, event.weight()); + } + } + + // "Bare" muon + if (!wfinder_bare_mu.particles().empty()) { + const FourMomentum mu = wfinder_bare_mu.constituentLeptons()[0].momentum(); + const FourMomentum nu = wfinder_bare_mu.constituentNeutrinos()[0].momentum(); + if (calcMT(mu, nu) > 40.0*GeV && nu.pT() > 25.0*GeV) { + _sumW_mu_bare += event.weight(); + _hist_wpt_bare_mu->fill(wfinder_bare_mu.bosons()[0].momentum().pT()/GeV, event.weight()); + } + } + + } + + + // Normalize histos + void finalize() { + scale(_hist_wpt_dressed_el, 1/_sumW_e_dressed); + scale(_hist_wpt_bare_el, 1/_sumW_e_bare); + scale(_hist_wpt_dressed_mu, 1/_sumW_mu_dressed); + scale(_hist_wpt_bare_mu, 1/_sumW_mu_bare); + } + + //@} + + + private: + + double _sumW_e_dressed, _sumW_e_bare, _sumW_mu_dressed, _sumW_mu_bare; + + AIDA::IHistogram1D* _hist_wpt_dressed_el; + AIDA::IHistogram1D* _hist_wpt_bare_el; + AIDA::IHistogram1D* _hist_wpt_dressed_mu; + AIDA::IHistogram1D* _hist_wpt_bare_mu; + + }; + + + // Hook for the plugin system + DECLARE_RIVET_PLUGIN(ATLAS_2011_I925932); + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Wed Nov 9 13:46:55 2011 (r3475) +++ trunk/src/Analyses/Makefile.am Wed Nov 9 16:24:45 2011 (r3476) @@ -60,8 +60,9 @@ ATLAS_2011_S9126244.cc \ ATLAS_2011_S9128077.cc \ ATLAS_2011_S9131140.cc \ + ATLAS_2011_S9212183.cc \ ATLAS_2011_I919017.cc \ - ATLAS_2011_S9212183.cc + ATLAS_2011_I925932.cc if ENABLE_UNVALIDATED RivetATLASAnalyses_la_SOURCES += \ ATLAS_2010_CONF_2010_049.cc \
More information about the Rivet-svn mailing list |