|
[Rivet-svn] r3061 - in trunk: . data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgTue Apr 19 18:52:58 BST 2011
Author: hoeth Date: Tue Apr 19 18:52:58 2011 New Revision: 3061 Log: Added ATLAS_2011_S9002537 W asymmetry analysis Added: trunk/data/anainfo/ATLAS_2011_S9002537.info trunk/data/plotinfo/ATLAS_2011_S9002537.plot trunk/data/refdata/ATLAS_2011_S9002537.aida trunk/src/Analyses/ATLAS_2011_S9002537.cc Modified: trunk/ChangeLog trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/src/Analyses/Makefile.am Modified: trunk/ChangeLog ============================================================================== --- trunk/ChangeLog Thu Apr 14 17:56:34 2011 (r3060) +++ trunk/ChangeLog Tue Apr 19 18:52:58 2011 (r3061) @@ -1,3 +1,7 @@ +2011-04-19 Hendrik Hoeth <hendrik.hoeth at cern.ch> + + * Added ATLAS_2011_S9002537 W asymmetry analysis + 2011-04-14 Hendrik Hoeth <hendrik.hoeth at cern.ch> * deltaR, deltaPhi, deltaEta now available in all combinations of Added: trunk/data/anainfo/ATLAS_2011_S9002537.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/ATLAS_2011_S9002537.info Tue Apr 19 18:52:58 2011 (r3061) @@ -0,0 +1,21 @@ +Name: ATLAS_2011_S9002537 +Year: 2011 +Summary: Muon charge asymmetry in W events at 7 TeV in ATLAS +Experiment: ATLAS +Collider: LHC +SpiresID: 9002537 +Status: VALIDATED +Authors: + - Frank Krauss <frank.krauss at durham.ac.uk> + - Hendrik Hoeth <hendrik.hoeth at cern.ch> +References: + - arXiv:1103.2929 +RunInfo: + W -> munu events at 7 TeV +NumEvents: 10000000 +Description: + Measurement of the muon charge asymmetry from W bosons produced in + proton-proton collisions at a centre-of-mass energy of 7 TeV with ATLAS. + The asymmetry is measured in the $W \rightarrow \mu$ decay mode as a function + of the muon pseudorapidity using a data sample corresponding to a total + integrated luminosity of 31pb$^-1$. Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Thu Apr 14 17:56:34 2011 (r3060) +++ trunk/data/anainfo/Makefile.am Tue Apr 19 18:52:58 2011 (r3061) @@ -17,6 +17,7 @@ ATLAS_2011_S8983313.info \ ATLAS_2010_CONF_2010_049.info \ ATLAS_2010_S8994773.info \ + ATLAS_2011_S9002537.info \ BELLE_2006_S6265367.info \ CDF_1988_S1865951.info \ CDF_1990_S2089246.info \ Added: trunk/data/plotinfo/ATLAS_2011_S9002537.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/ATLAS_2011_S9002537.plot Tue Apr 19 18:52:58 2011 (r3061) @@ -0,0 +1,7 @@ +# BEGIN PLOT /ATLAS_2011_S9002537/d01-x01-y01 +LogY=0 +LegendXPos=0.1 +Title=Muon charge asymmetry in $W$ decays +XLabel=$|\eta_\mu|$ +YLabel=$\mathcal{A}_\mu$ +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Thu Apr 14 17:56:34 2011 (r3060) +++ trunk/data/plotinfo/Makefile.am Tue Apr 19 18:52:58 2011 (r3061) @@ -14,8 +14,9 @@ ATLAS_2010_S8919674.plot \ ATLAS_2011_S8924791.plot \ ATLAS_2011_S8971293.plot \ - ATLAS_2010_CONF_2010_049.plot \ ATLAS_2010_S8994773.plot \ + ATLAS_2011_S9002537.plot \ + ATLAS_2010_CONF_2010_049.plot \ CDF_1988_S1865951.plot \ CDF_1990_S2089246.plot \ CDF_1991_S2313472.plot \ Added: trunk/data/refdata/ATLAS_2011_S9002537.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/ATLAS_2011_S9002537.aida Tue Apr 19 18:52:58 2011 (r3061) @@ -0,0 +1,73 @@ +<?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_2011_S9002537" title="Muon charge asymmetry in $W$ decays"> + <annotation> + <item key="Title" value="Muon charge asymmetry in $W$ decays" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2011_S9002537/d01-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2011_S9002537.aida/REF/ATLAS_2011_S9002537/d01-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.050000e-01" value="1.050000e-01" errorMinus="1.050000e-01"/> + <measurement errorPlus="2.024000e-02" value="1.470000e-01" errorMinus="2.024000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.050000e-01" value="3.150000e-01" errorMinus="1.050000e-01"/> + <measurement errorPlus="1.562000e-02" value="1.500000e-01" errorMinus="1.562000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.050000e-01" value="5.250000e-01" errorMinus="1.050000e-01"/> + <measurement errorPlus="1.562000e-02" value="1.580000e-01" errorMinus="1.562000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.050000e-01" value="7.350000e-01" errorMinus="1.050000e-01"/> + <measurement errorPlus="1.802000e-02" value="1.840000e-01" errorMinus="1.802000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.050000e-01" value="9.450000e-01" errorMinus="1.050000e-01"/> + <measurement errorPlus="2.024000e-02" value="1.860000e-01" errorMinus="2.024000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.600000e-01" value="1.210000e+00" errorMinus="1.600000e-01"/> + <measurement errorPlus="1.360000e-02" value="2.400000e-01" errorMinus="1.360000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.500000e-02" value="1.445000e+00" errorMinus="7.500000e-02"/> + <measurement errorPlus="1.486000e-02" value="2.500000e-01" errorMinus="1.486000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.100000e-01" value="1.630000e+00" errorMinus="1.100000e-01"/> + <measurement errorPlus="1.345000e-02" value="2.690000e-01" errorMinus="1.345000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.050000e-01" value="1.845000e+00" errorMinus="1.050000e-01"/> + <measurement errorPlus="1.345000e-02" value="2.730000e-01" errorMinus="1.345000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.150000e-01" value="2.065000e+00" errorMinus="1.150000e-01"/> + <measurement errorPlus="1.499000e-02" value="2.760000e-01" errorMinus="1.499000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.100000e-01" value="2.290000e+00" errorMinus="1.100000e-01"/> + <measurement errorPlus="1.562000e-02" value="2.730000e-01" errorMinus="1.562000e-02"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" + path="/REF/ATLAS_2011_S9002537" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2011_S9002537/d02-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2011_S9002537.aida/REF/ATLAS_2011_S9002537/d02-x01-y01" sticky="true"/> + </annotation> + </dataPointSet> + <dataPointSet name="d02-x02-y01" dimension="2" + path="/REF/ATLAS_2011_S9002537" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2011_S9002537/d02-x02-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2011_S9002537.aida/REF/ATLAS_2011_S9002537/d02-x02-y01" sticky="true"/> + </annotation> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Thu Apr 14 17:56:34 2011 (r3060) +++ trunk/data/refdata/Makefile.am Tue Apr 19 18:52:58 2011 (r3061) @@ -18,8 +18,9 @@ ATLAS_2010_S8919674.aida \ ATLAS_2011_S8924791.aida \ ATLAS_2011_S8971293.aida \ + ATLAS_2010_S8994773.aida \ + ATLAS_2011_S9002537.aida \ ATLAS_2010_CONF_2010_049.aida \ - ATLAS_2010_S8994773.aida \ CMS_2011_S8957746.aida \ DELPHI_1994_S3021912.aida \ DELPHI_1995_S3137023.aida \ Added: trunk/src/Analyses/ATLAS_2011_S9002537.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/ATLAS_2011_S9002537.cc Tue Apr 19 18:52:58 2011 (r3061) @@ -0,0 +1,109 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/IdentifiedFinalState.hh" +#include "Rivet/Projections/ChargedFinalState.hh" +#include "Rivet/Projections/MissingMomentum.hh" +#include "Rivet/Projections/FinalState.hh" + +namespace Rivet { + + class ATLAS_2011_S9002537 : public Analysis { + + public: + ATLAS_2011_S9002537(): Analysis("ATLAS_2011_S9002537") + { + } + + public: + void init() { + IdentifiedFinalState Muons(-2.4,2.4,20.*GeV); + Muons.acceptIdPair(MUON); + addProjection(Muons,"muons"); + + ChargedFinalState CFS(-2.8,2.8,0.*GeV); + addProjection(CFS,"tracks"); + + MissingMomentum missmom(FinalState(-5.,5.,0.*GeV)); + addProjection(missmom,"MissingMomentum"); + + _h_plus = bookHistogram1D("_h_plus", binEdges(1,1,1)); + _h_minus = bookHistogram1D("_h_minus", binEdges(1,1,1)); + _h_asym = bookDataPointSet(1,1,1); + } + + void analyze(const Event& event) { + const IdentifiedFinalState& muons = + applyProjection<IdentifiedFinalState>(event, "muons"); + + const ChargedFinalState& tracks = + applyProjection<ChargedFinalState>(event, "tracks"); + + if (muons.size()<1) vetoEvent; + ParticleVector selected_muons; + foreach (Particle muon, muons.particles()) { + FourMomentum testmom = muon.momentum(); + double ptmu(testmom.pT()), ptsum(-ptmu), ratio(0.); + foreach (Particle track,tracks.particles()) { + FourMomentum trackmom = track.momentum(); + if (deltaR(testmom,trackmom)<0.4) { + ptsum += trackmom.pT(); + ratio = ptsum/ptmu; + if (ratio>0.2) + break; + } + } + if (ratio<0.2) + selected_muons.push_back(muon); + } + if (selected_muons.size()!=1) vetoEvent; + + const FourMomentum muonmom = selected_muons.begin()->momentum(); + const MissingMomentum& missmom = applyProjection<MissingMomentum>(event, "MissingMomentum"); + FourMomentum missvec = -missmom.visibleMomentum(); + if (fabs(missvec.Et())<25) vetoEvent; + + double MTW = sqrt(2.*missvec.pT()*muonmom.pT()*(1.-cos(deltaPhi(missvec.phi(),muonmom.phi())))); + if (MTW<40.*GeV) vetoEvent; + + if (selected_muons.begin()->pdgId()>0) + _h_minus->fill(muonmom.eta(),event.weight()); + else + _h_plus->fill(muonmom.eta(),event.weight()); + } + + + /// Normalise histograms etc., after the run + void finalize() { + int Nbins = _h_plus->axis().bins(); + std::vector<double> asym, asym_err; + for (int i=0; i<Nbins; i++) { + double num(0.), denom(0.), err(0.); + num = _h_plus->binHeight(i) - _h_minus->binHeight(i); + denom = _h_plus->binHeight(i) + _h_minus->binHeight(i); + err = _h_plus->binError(i) + _h_minus->binError(i); + + if (num==0 || denom==0) { + asym.push_back(0.); + asym_err.push_back(0.); + } + else { + asym.push_back(num/denom); + asym_err.push_back(num/denom*((err/num)+(err/denom))); + } + } + _h_asym->setCoordinate(1, asym, asym_err); + + histogramFactory().destroy(_h_plus); + histogramFactory().destroy(_h_minus); + } + + private: + AIDA::IHistogram1D *_h_plus, *_h_minus; + AIDA::IDataPointSet *_h_asym; + + }; + + AnalysisBuilder<ATLAS_2011_S9002537> plugin_ATLAS_2011_S9002537; +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Thu Apr 14 17:56:34 2011 (r3060) +++ trunk/src/Analyses/Makefile.am Tue Apr 19 18:52:58 2011 (r3061) @@ -52,7 +52,8 @@ ATLAS_2011_S8924791.cc \ ATLAS_2011_S8971293.cc \ ATLAS_2010_S8994773.cc \ - ATLAS_2011_S8983313.cc + ATLAS_2011_S8983313.cc \ + ATLAS_2011_S9002537.cc if ENABLE_UNVALIDATED RivetATLASAnalyses_la_SOURCES += \ ATLAS_2010_CONF_2010_049.cc
More information about the Rivet-svn mailing list |