|
[Rivet-svn] r3940 - branches/2012-06-aidarivet/data/anainfo branches/2012-06-aidarivet/data/plotinfo branches/2012-06-aidarivet/data/refdata branches/2012-06-aidarivet/src/Analyses trunk/data/anainfo trunk/data/plotinfo trunk/data/refdata trunk/src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgWed Sep 26 20:46:45 BST 2012
Author: fsiegert Date: Wed Sep 26 20:46:45 2012 New Revision: 3940 Log: Add ATLAS_2011_I954993 to trunk and 2012-06-aidarivet. (Thanks, Roman.) Added: branches/2012-06-aidarivet/data/anainfo/ATLAS_2011_I954993.info branches/2012-06-aidarivet/data/plotinfo/ATLAS_2011_I954993.plot branches/2012-06-aidarivet/data/refdata/ATLAS_2011_I954993.aida branches/2012-06-aidarivet/src/Analyses/ATLAS_2011_I954993.cc trunk/data/anainfo/ATLAS_2011_I954993.info trunk/data/plotinfo/ATLAS_2011_I954993.plot trunk/data/refdata/ATLAS_2011_I954993.aida trunk/src/Analyses/ATLAS_2011_I954993.cc Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am branches/2012-06-aidarivet/data/plotinfo/Makefile.am branches/2012-06-aidarivet/data/refdata/Makefile.am branches/2012-06-aidarivet/src/Analyses/Makefile.am trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/src/Analyses/Makefile.am Added: branches/2012-06-aidarivet/data/anainfo/ATLAS_2011_I954993.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/anainfo/ATLAS_2011_I954993.info Wed Sep 26 20:46:45 2012 (r3940) @@ -0,0 +1,48 @@ +Name: ATLAS_2011_I954993 +Year: 2011 +Summary: WZ fiducial cross-section at 7 TeV in ATLAS +Experiment: ATLAS +Collider: LHC 7TeV +InspireID: 954993 +Status: VALIDATED +Authors: + - Lynn Marx <Lynn.Marx at hep.manchester.ac.uk> + - Roman Lysak <lysak at fzu.cz> +References: + - Phys.Lett. B709 (2012) 341-357 + - arXiv:1111.5570 +RunInfo: + pp WZ events at 7 TeV with direct e, mu W/Z boson decays (no taus from W/Z) +NumEvents: 1000000 +NeedCrossSection: yes +Beams: [p+, p+] +Energies: [7000] +PtCuts: [0] +Description: + 'This is a measurement of WZ production in 1.02 $fb^{-1}$ of + pp collision data at $\sqrt{s} = $7 TeV collected by the ATLAS experiment + in 2011. Doubly leptonic decay events are selected with electrons, + muons and missing transverse momentum in the final state. + The measurement of the combined fiducial cross section for the WZ + bosons decaying directly into electrons and muons is performed.' +BibKey: Aad:2011cx +BibTeX: '@article{Aad:2011cx, + author = "Aad, Georges and others", + title = "{Measurement of the WZ production cross section and + limits on anomalous triple gauge couplings in + proton-proton collisions at sqrt(s) = 7 TeV with + the ATLAS detector}", + collaboration = "ATLAS Collaboration", + journal = "Phys.Lett.", + volume = "B709", + pages = "341-357", + doi = "10.1016/j.physletb.2012.02.053", + year = "2012", + eprint = "1111.5570", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "CERN-PH-EP-2011-184", + SLACcitation = "%%CITATION = ARXIV:1111.5570;%%", +}' + + Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/anainfo/Makefile.am Wed Sep 26 20:01:04 2012 (r3939) +++ branches/2012-06-aidarivet/data/anainfo/Makefile.am Wed Sep 26 20:46:45 2012 (r3940) @@ -41,6 +41,7 @@ ATLAS_2011_I926145.info \ ATLAS_2011_I944826.info \ ATLAS_2011_I945498.info \ + ATLAS_2011_I954993.info \ ATLAS_2011_S9225137.info \ ATLAS_2011_S9212353.info \ ATLAS_2012_I1083318.info \ Added: branches/2012-06-aidarivet/data/plotinfo/ATLAS_2011_I954993.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/plotinfo/ATLAS_2011_I954993.plot Wed Sep 26 20:46:45 2012 (r3940) @@ -0,0 +1,5 @@ +# BEGIN PLOT /ATLAS_2011_I954993/d01-x01-y01 +Title=Total fiducial cross-section $WZ\to\ell\nu\ell\ell$ +XLabel=$\sqrt{s}$~[GeV] +YLabel=$\sigma_{WZ\to\ell\nu\ell\ell}$~[fb] +# END PLOT Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/plotinfo/Makefile.am Wed Sep 26 20:01:04 2012 (r3939) +++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am Wed Sep 26 20:46:45 2012 (r3940) @@ -34,6 +34,7 @@ ATLAS_2011_I926145.plot \ ATLAS_2011_I944826.plot \ ATLAS_2011_I945498.plot \ + ATLAS_2011_I954993.plot \ ATLAS_2011_S9225137.plot \ ATLAS_2011_S9212183.plot \ ATLAS_2011_CONF_2011_090.plot \ Added: branches/2012-06-aidarivet/data/refdata/ATLAS_2011_I954993.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/refdata/ATLAS_2011_I954993.aida Wed Sep 26 20:46:45 2012 (r3940) @@ -0,0 +1,16 @@ +<?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_I954993" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2011_I954993/d01-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement value="7.000000e+03" errorPlus="5.000000e-01" errorMinus="5.000000e-01"/> + <measurement value="1.020000e+02" errorPlus="1.700000e+01" errorMinus="1.570000e+01"/> + </dataPoint> + </dataPointSet> +</aida> Modified: branches/2012-06-aidarivet/data/refdata/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/refdata/Makefile.am Wed Sep 26 20:01:04 2012 (r3939) +++ branches/2012-06-aidarivet/data/refdata/Makefile.am Wed Sep 26 20:46:45 2012 (r3940) @@ -40,6 +40,7 @@ ATLAS_2011_I926145.aida \ ATLAS_2011_I944826.aida \ ATLAS_2011_I945498.aida \ + ATLAS_2011_I954993.aida \ ATLAS_2011_S9225137.aida \ ATLAS_2011_S9212183.aida \ ATLAS_2012_I1082936.aida \ Added: branches/2012-06-aidarivet/src/Analyses/ATLAS_2011_I954993.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/src/Analyses/ATLAS_2011_I954993.cc Wed Sep 26 20:46:45 2012 (r3940) @@ -0,0 +1,139 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/ZFinder.hh" +#include "Rivet/Projections/WFinder.hh" +#include "Rivet/Projections/VetoedFinalState.hh" + +namespace Rivet { + + + /// @brief MC validation analysis for WZ events + class ATLAS_2011_I954993 : public Analysis { + public: + + /// Default constructor + ATLAS_2011_I954993() + : Analysis("ATLAS_2011_I954993") + { + setNeedsCrossSection(true); + } + + + /// @name Analysis methods + //@{ + + /// Book histograms + void init() { + + //// ZFinder: etaMin,etaMax,pTmin,pid,m2_min,m2_max,dRmax,clusterPhotons,excludePhotonsFromRFS + ZFinder zfinder_e( -2.5, 2.5, 15.0*GeV, ELECTRON, 81.1876*GeV, 101.1876*GeV, 0.1, true, true); + addProjection(zfinder_e, "ZFinder_e"); + ZFinder zfinder_mu(-2.5, 2.5, 15.0*GeV, MUON, 81.1876*GeV, 101.1876*GeV, 0.1, true, true); + addProjection(zfinder_mu, "ZFinder_mu"); + + //// WFinder: etaRanges,pTmin,pid,m2_min,m2_max,missingET,dRmax + VetoedFinalState weinput; + weinput.addVetoOnThisFinalState(zfinder_e); + WFinder wfinder_e(weinput, -2.5, 2.5, 15.0*GeV, ELECTRON, 0.0*GeV, 1000.0*GeV, 25.0*GeV, 0.1); + addProjection(wfinder_e, "WFinder_e"); + + VetoedFinalState wminput; + wminput.addVetoOnThisFinalState(zfinder_mu); + WFinder wfinder_mu(wminput,-2.5, 2.5, 15.0*GeV, MUON, 0.0*GeV, 1000.0*GeV, 25.0*GeV, 0.1); + addProjection(wfinder_mu, "WFinder_mu"); + + //// Histograms + _h_fiducial = bookHistogram1D(1,1,1); + + } + + /// Do the analysis + void analyze(const Event & e) { + + const double weight = e.weight(); + + const ZFinder& zfinder_e = applyProjection<ZFinder>(e, "ZFinder_e"); + const ZFinder& zfinder_mu = applyProjection<ZFinder>(e, "ZFinder_mu"); + const WFinder& wfinder_e = applyProjection<WFinder>(e, "WFinder_e"); + const WFinder& wfinder_mu = applyProjection<WFinder>(e, "WFinder_mu"); + + + // Looking for a Z + if (zfinder_e.bosons().size()!= 1 && zfinder_mu.bosons().size() != 1) { + MSG_DEBUG("No Z boson found, vetoing event"); + vetoEvent; + } + + // Looking for a W + if (wfinder_e.bosons().size()!= 1 && wfinder_mu.bosons().size() != 1) { + MSG_DEBUG("No W boson found, vetoing event"); + vetoEvent; + } + + // If we find a W... + FourMomentum wmom_e(0.0,0.0,0.0,0.0), We(0.0,0.0,0.0,0.0), Wenu(0.0,0.0,0.0,0.0); + FourMomentum wmom_mu(0.0,0.0,0.0,0.0), Wmu(0.0,0.0,0.0,0.0), Wmunu(0.0,0.0,0.0,0.0); + if(wfinder_e.bosons().size()== 1){ + wmom_e = wfinder_e.bosons().front().momentum(); + We = wfinder_e.constituentLeptons()[0].momentum(); + Wenu = wfinder_e.constituentNeutrinos()[0].momentum(); + } + if(wfinder_mu.bosons().size()== 1){ + wmom_mu = wfinder_mu.bosons().front().momentum(); + Wmu = wfinder_mu.constituentLeptons()[0].momentum(); + Wmunu = wfinder_mu.constituentNeutrinos()[0].momentum(); + } + + // Applying remaining fiducial phase space requirements + double mT = 0; + if(wfinder_e.bosons().size() == 1){ + mT = sqrt(2*We.pT()*Wenu.Et()*(1.0-cos(We.phi()-Wenu.phi()))); + if (Wenu.pT()/GeV < 25.0 || We.pT()/GeV < 20.0 || mT/GeV < 20.0) { + MSG_DEBUG(" Wnu.pT()/GeV:" << Wenu.pT()/GeV<<" Wl.pT()/GeV:" << We.pT()/GeV<<" mT/GeV:" << mT/GeV); + vetoEvent; + } + } + else if(wfinder_mu.bosons().size() == 1){ + mT = sqrt(2*Wmu.pT()*Wmunu.Et()*(1.0-cos(Wmu.phi()-Wmunu.phi()))); + if (Wmunu.pT()/GeV < 25.0 || Wmu.pT()/GeV < 20.0 || mT/GeV < 20.0) { + MSG_DEBUG(" Wnu.pT()/GeV:" << Wmunu.pT()/GeV<<" Wl.pT()/GeV:" << Wmu.pT()/GeV<<" mT/GeV:" << mT/GeV); + vetoEvent; + } + } + else{ + MSG_DEBUG("No W boson found, can't make a transverse mass, vetoing event"); + vetoEvent; + } + + _h_fiducial->fill(7000.0, weight); + + } + + + /// Finalize + void finalize() { + + scale(_h_fiducial, crossSection()/femtobarn/sumOfWeights()); + + } + + //@} + + + private: + + /// @name Histograms + //@{ + AIDA::IHistogram1D* _h_fiducial; + + //@} + + }; + + + //// The hook for the plugin system + DECLARE_RIVET_PLUGIN(ATLAS_2011_I954993); + +} Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/src/Analyses/Makefile.am Wed Sep 26 20:01:04 2012 (r3939) +++ branches/2012-06-aidarivet/src/Analyses/Makefile.am Wed Sep 26 20:46:45 2012 (r3940) @@ -68,6 +68,7 @@ ATLAS_2011_I926145.cc \ ATLAS_2011_I944826.cc \ ATLAS_2011_I945498.cc \ + ATLAS_2011_I954993.cc \ ATLAS_2012_I943401.cc \ ATLAS_2012_I1082936.cc \ ATLAS_2012_I1084540.cc \ Added: trunk/data/anainfo/ATLAS_2011_I954993.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/ATLAS_2011_I954993.info Wed Sep 26 20:46:45 2012 (r3940) @@ -0,0 +1,48 @@ +Name: ATLAS_2011_I954993 +Year: 2011 +Summary: WZ fiducial cross-section at 7 TeV in ATLAS +Experiment: ATLAS +Collider: LHC 7TeV +InspireID: 954993 +Status: VALIDATED +Authors: + - Lynn Marx <Lynn.Marx at hep.manchester.ac.uk> + - Roman Lysak <lysak at fzu.cz> +References: + - Phys.Lett. B709 (2012) 341-357 + - arXiv:1111.5570 +RunInfo: + pp WZ events at 7 TeV with direct e, mu W/Z boson decays (no taus from W/Z) +NumEvents: 1000000 +NeedCrossSection: yes +Beams: [p+, p+] +Energies: [7000] +PtCuts: [0] +Description: + 'This is a measurement of WZ production in 1.02 $fb^{-1}$ of + pp collision data at $\sqrt{s} = $7 TeV collected by the ATLAS experiment + in 2011. Doubly leptonic decay events are selected with electrons, + muons and missing transverse momentum in the final state. + The measurement of the combined fiducial cross section for the WZ + bosons decaying directly into electrons and muons is performed.' +BibKey: Aad:2011cx +BibTeX: '@article{Aad:2011cx, + author = "Aad, Georges and others", + title = "{Measurement of the WZ production cross section and + limits on anomalous triple gauge couplings in + proton-proton collisions at sqrt(s) = 7 TeV with + the ATLAS detector}", + collaboration = "ATLAS Collaboration", + journal = "Phys.Lett.", + volume = "B709", + pages = "341-357", + doi = "10.1016/j.physletb.2012.02.053", + year = "2012", + eprint = "1111.5570", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "CERN-PH-EP-2011-184", + SLACcitation = "%%CITATION = ARXIV:1111.5570;%%", +}' + + Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Wed Sep 26 20:01:04 2012 (r3939) +++ trunk/data/anainfo/Makefile.am Wed Sep 26 20:46:45 2012 (r3940) @@ -41,6 +41,7 @@ ATLAS_2011_I926145.info \ ATLAS_2011_I944826.info \ ATLAS_2011_I945498.info \ + ATLAS_2011_I954993.info \ ATLAS_2011_S9225137.info \ ATLAS_2011_S9212353.info \ ATLAS_2012_I1083318.info \ Added: trunk/data/plotinfo/ATLAS_2011_I954993.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/ATLAS_2011_I954993.plot Wed Sep 26 20:46:45 2012 (r3940) @@ -0,0 +1,5 @@ +# BEGIN PLOT /ATLAS_2011_I954993/d01-x01-y01 +Title=Total fiducial cross-section $WZ\to\ell\nu\ell\ell$ +XLabel=$\sqrt{s}$~[GeV] +YLabel=$\sigma_{WZ\to\ell\nu\ell\ell}$~[fb] +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Wed Sep 26 20:01:04 2012 (r3939) +++ trunk/data/plotinfo/Makefile.am Wed Sep 26 20:46:45 2012 (r3940) @@ -34,6 +34,7 @@ ATLAS_2011_I926145.plot \ ATLAS_2011_I944826.plot \ ATLAS_2011_I945498.plot \ + ATLAS_2011_I954993.plot \ ATLAS_2011_S9225137.plot \ ATLAS_2011_S9212183.plot \ ATLAS_2011_CONF_2011_090.plot \ Added: trunk/data/refdata/ATLAS_2011_I954993.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/ATLAS_2011_I954993.aida Wed Sep 26 20:46:45 2012 (r3940) @@ -0,0 +1,16 @@ +<?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_I954993" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2011_I954993/d01-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement value="7.000000e+03" errorPlus="5.000000e-01" errorMinus="5.000000e-01"/> + <measurement value="1.020000e+02" errorPlus="1.700000e+01" errorMinus="1.570000e+01"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Wed Sep 26 20:01:04 2012 (r3939) +++ trunk/data/refdata/Makefile.am Wed Sep 26 20:46:45 2012 (r3940) @@ -40,6 +40,7 @@ ATLAS_2011_I926145.aida \ ATLAS_2011_I944826.aida \ ATLAS_2011_I945498.aida \ + ATLAS_2011_I954993.aida \ ATLAS_2011_S9225137.aida \ ATLAS_2011_S9212183.aida \ ATLAS_2012_I1082936.aida \ Added: trunk/src/Analyses/ATLAS_2011_I954993.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/ATLAS_2011_I954993.cc Wed Sep 26 20:46:45 2012 (r3940) @@ -0,0 +1,139 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetYODA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/ZFinder.hh" +#include "Rivet/Projections/WFinder.hh" +#include "Rivet/Projections/VetoedFinalState.hh" + +namespace Rivet { + + + /// @brief MC validation analysis for WZ events + class ATLAS_2011_I954993 : public Analysis { + public: + + /// Default constructor + ATLAS_2011_I954993() + : Analysis("ATLAS_2011_I954993") + { + setNeedsCrossSection(true); + } + + + /// @name Analysis methods + //@{ + + /// Book histograms + void init() { + + //// ZFinder: etaMin,etaMax,pTmin,pid,m2_min,m2_max,dRmax,clusterPhotons,excludePhotonsFromRFS + ZFinder zfinder_e( -2.5, 2.5, 15.0*GeV, ELECTRON, 81.1876*GeV, 101.1876*GeV, 0.1, true, true); + addProjection(zfinder_e, "ZFinder_e"); + ZFinder zfinder_mu(-2.5, 2.5, 15.0*GeV, MUON, 81.1876*GeV, 101.1876*GeV, 0.1, true, true); + addProjection(zfinder_mu, "ZFinder_mu"); + + //// WFinder: etaRanges,pTmin,pid,m2_min,m2_max,missingET,dRmax + VetoedFinalState weinput; + weinput.addVetoOnThisFinalState(zfinder_e); + WFinder wfinder_e(weinput, -2.5, 2.5, 15.0*GeV, ELECTRON, 0.0*GeV, 1000.0*GeV, 25.0*GeV, 0.1); + addProjection(wfinder_e, "WFinder_e"); + + VetoedFinalState wminput; + wminput.addVetoOnThisFinalState(zfinder_mu); + WFinder wfinder_mu(wminput,-2.5, 2.5, 15.0*GeV, MUON, 0.0*GeV, 1000.0*GeV, 25.0*GeV, 0.1); + addProjection(wfinder_mu, "WFinder_mu"); + + //// Histograms + _h_fiducial = bookHisto1D(1,1,1); + + } + + /// Do the analysis + void analyze(const Event & e) { + + const double weight = e.weight(); + + const ZFinder& zfinder_e = applyProjection<ZFinder>(e, "ZFinder_e"); + const ZFinder& zfinder_mu = applyProjection<ZFinder>(e, "ZFinder_mu"); + const WFinder& wfinder_e = applyProjection<WFinder>(e, "WFinder_e"); + const WFinder& wfinder_mu = applyProjection<WFinder>(e, "WFinder_mu"); + + + // Looking for a Z + if (zfinder_e.bosons().size()!= 1 && zfinder_mu.bosons().size() != 1) { + MSG_DEBUG("No Z boson found, vetoing event"); + vetoEvent; + } + + // Looking for a W + if (wfinder_e.bosons().size()!= 1 && wfinder_mu.bosons().size() != 1) { + MSG_DEBUG("No W boson found, vetoing event"); + vetoEvent; + } + + // If we find a W... + FourMomentum wmom_e(0.0,0.0,0.0,0.0), We(0.0,0.0,0.0,0.0), Wenu(0.0,0.0,0.0,0.0); + FourMomentum wmom_mu(0.0,0.0,0.0,0.0), Wmu(0.0,0.0,0.0,0.0), Wmunu(0.0,0.0,0.0,0.0); + if(wfinder_e.bosons().size()== 1){ + wmom_e = wfinder_e.bosons().front().momentum(); + We = wfinder_e.constituentLeptons()[0].momentum(); + Wenu = wfinder_e.constituentNeutrinos()[0].momentum(); + } + if(wfinder_mu.bosons().size()== 1){ + wmom_mu = wfinder_mu.bosons().front().momentum(); + Wmu = wfinder_mu.constituentLeptons()[0].momentum(); + Wmunu = wfinder_mu.constituentNeutrinos()[0].momentum(); + } + + // Applying remaining fiducial phase space requirements + double mT = 0; + if(wfinder_e.bosons().size() == 1){ + mT = sqrt(2*We.pT()*Wenu.Et()*(1.0-cos(We.phi()-Wenu.phi()))); + if (Wenu.pT()/GeV < 25.0 || We.pT()/GeV < 20.0 || mT/GeV < 20.0) { + MSG_DEBUG(" Wnu.pT()/GeV:" << Wenu.pT()/GeV<<" Wl.pT()/GeV:" << We.pT()/GeV<<" mT/GeV:" << mT/GeV); + vetoEvent; + } + } + else if(wfinder_mu.bosons().size() == 1){ + mT = sqrt(2*Wmu.pT()*Wmunu.Et()*(1.0-cos(Wmu.phi()-Wmunu.phi()))); + if (Wmunu.pT()/GeV < 25.0 || Wmu.pT()/GeV < 20.0 || mT/GeV < 20.0) { + MSG_DEBUG(" Wnu.pT()/GeV:" << Wmunu.pT()/GeV<<" Wl.pT()/GeV:" << Wmu.pT()/GeV<<" mT/GeV:" << mT/GeV); + vetoEvent; + } + } + else{ + MSG_DEBUG("No W boson found, can't make a transverse mass, vetoing event"); + vetoEvent; + } + + _h_fiducial->fill(7000.0, weight); + + } + + + /// Finalize + void finalize() { + + scale(_h_fiducial, crossSection()/femtobarn/sumOfWeights()); + + } + + //@} + + + private: + + /// @name Histograms + //@{ + Histo1DPtr _h_fiducial; + + //@} + + }; + + + //// The hook for the plugin system + DECLARE_RIVET_PLUGIN(ATLAS_2011_I954993); + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Wed Sep 26 20:01:04 2012 (r3939) +++ trunk/src/Analyses/Makefile.am Wed Sep 26 20:46:45 2012 (r3940) @@ -68,6 +68,7 @@ ATLAS_2011_I926145.cc \ ATLAS_2011_I944826.cc \ ATLAS_2011_I945498.cc \ + ATLAS_2011_I954993.cc \ ATLAS_2012_I943401.cc \ ATLAS_2012_I1082936.cc \ ATLAS_2012_I1084540.cc \
More information about the Rivet-svn mailing list |