|
[Rivet-svn] r2654 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgMon Aug 16 11:57:11 BST 2010
Author: fsiegert Date: Mon Aug 16 11:57:11 2010 New Revision: 2654 Log: Add analysis D0_2000_S4480767 for W pT in Tevatron Run 1. Holger: I think this makes D0_1998_S3711838 obsolete (not detector corrected, histo normalisation not implemented yet in finalize, UNVALIDATED, ...); shall we remove it? Added: trunk/data/anainfo/D0_2000_S4480767.info trunk/data/plotinfo/D0_2000_S4480767.plot trunk/data/refdata/D0_2000_S4480767.aida trunk/src/Analyses/D0_2000_S4480767.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/D0_2000_S4480767.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/D0_2000_S4480767.info Mon Aug 16 11:57:11 2010 (r2654) @@ -0,0 +1,40 @@ +Name: D0_2000_S4480767 +Year: 2000 +Summary: Transverse momentum of the W boson +Experiment: D0 +Collider: Tevatron Run 1 +SpiresID: 4480767 +Status: VALIDATED +Authors: + - Frank Siegert <frank.siegert at cern.ch> +References: + - Phys.Lett. B513 (2001) 292-300 + - doi:10.1016/S0370-2693(01)00628-1 + - arXiv:hep-ex/0010026 +RunInfo: + Production of W+ and W- decaying into the electron channel. +NumEvents: 1000000 +Beams: [p-, p+] +Energies: [1800] +PtCuts: [0] +Description: + Measurement of the differential cross section for W boson production as a + function of its transverse momentum. The data were collected by the D0 + experiment at the Fermilab Tevatron Collider during 1994-1995 and correspond + to an integrated luminosity of 85 pb$^{-1}$. +BibKey: Abbott:2000xv +BibTeX: '@Article{Abbott:2000xv, + author = "Abbott, B. and others", + collaboration = "D0", + title = "{Differential cross section for $W$ boson production as a + function of transverse momentum in $p\bar{p}$ collisions at + $\sqrt{s} = 1.8$ TeV}", + journal = "Phys. Lett.", + volume = "B513", + year = "2001", + pages = "292-300", + eprint = "hep-ex/0010026", + archivePrefix = "arXiv", + doi = "10.1016/S0370-2693(01)00628-1", + SLACcitation = "%%CITATION = HEP-EX/0010026;%%" +}' Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Thu Aug 12 23:09:30 2010 (r2653) +++ trunk/data/anainfo/Makefile.am Mon Aug 16 11:57:11 2010 (r2654) @@ -42,6 +42,7 @@ D0_1996_S3214044.info \ D0_1996_S3324664.info \ D0_1998_S3711838.info \ + D0_2000_S4480767.info \ D0_2001_S4674421.info \ D0_2004_S5992206.info \ D0_2006_S6438750.info \ Added: trunk/data/plotinfo/D0_2000_S4480767.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/D0_2000_S4480767.plot Mon Aug 16 11:57:11 2010 (r2654) @@ -0,0 +1,5 @@ +# BEGIN PLOT /D0_2000_S4480767/d01-x01-y01 +Title=W boson pT +XLabel=$p_{\perp}^W$ [GeV] +YLabel=$\text{d}\sigma/\text{d}p_\perp^W$ [pb/GeV] +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Thu Aug 12 23:09:30 2010 (r2653) +++ trunk/data/plotinfo/Makefile.am Mon Aug 16 11:57:11 2010 (r2654) @@ -41,6 +41,7 @@ D0_1996_S3214044.plot \ D0_1996_S3324664.plot \ D0_1998_S3711838.plot \ + D0_2000_S4480767.plot \ D0_2001_S4674421.plot \ D0_2004_S5992206.plot \ D0_2006_S6438750.plot \ Added: trunk/data/refdata/D0_2000_S4480767.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/D0_2000_S4480767.aida Mon Aug 16 11:57:11 2010 (r2654) @@ -0,0 +1,95 @@ +<?xml version="1.0" encoding="utf-8" ?> +<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd"> +<aida version="3.3"> + <implementation version="1.0" package="HepData"/> + <dataPointSet name="d01-x01-y01" dimension="2" path="/REF/D0_2000_S4480767" title="D(SIG)/DPT IN PB/GEV" > + <dataPoint> + <measurement value="0.92" errorPlus="1.08" errorMinus="0.92"/> + <measurement value="109.37" errorPlus="12.551001505792277" errorMinus="12.551001505792277"/> + </dataPoint> + <dataPoint> + <measurement value="3.4" errorPlus="0.6000000000000001" errorMinus="1.4"/> + <measurement value="205.91" errorPlus="25.469784545645457" errorMinus="25.469784545645457"/> + </dataPoint> + <dataPoint> + <measurement value="4.97" errorPlus="1.0300000000000002" errorMinus="0.9699999999999998"/> + <measurement value="171.28" errorPlus="13.134356441881728" errorMinus="13.134356441881728"/> + </dataPoint> + <dataPoint> + <measurement value="6.98" errorPlus="1.0199999999999996" errorMinus="0.9800000000000004"/> + <measurement value="133.62" errorPlus="12.346033100490216" errorMinus="12.346033100490216"/> + </dataPoint> + <dataPoint> + <measurement value="8.98" errorPlus="1.0199999999999996" errorMinus="0.9800000000000004"/> + <measurement value="103.3" errorPlus="9.397267849752927" errorMinus="9.397267849752927"/> + </dataPoint> + <dataPoint> + <measurement value="10.98" errorPlus="1.0199999999999996" errorMinus="0.9800000000000004"/> + <measurement value="77.58" errorPlus="8.649596452459502" errorMinus="8.649596452459502"/> + </dataPoint> + <dataPoint> + <measurement value="12.98" errorPlus="1.0199999999999996" errorMinus="0.9800000000000004"/> + <measurement value="63.66" errorPlus="5.9684441089449765" errorMinus="5.9684441089449765"/> + </dataPoint> + <dataPoint> + <measurement value="14.98" errorPlus="1.0199999999999996" errorMinus="0.9800000000000004"/> + <measurement value="47.88" errorPlus="5.324666107691637" errorMinus="5.324666107691637"/> + </dataPoint> + <dataPoint> + <measurement value="16.98" errorPlus="1.0199999999999996" errorMinus="0.9800000000000004"/> + <measurement value="37.72" errorPlus="3.861274103505215" errorMinus="3.861274103505215"/> + </dataPoint> + <dataPoint> + <measurement value="18.98" errorPlus="1.0199999999999996" errorMinus="0.9800000000000004"/> + <measurement value="30.65" errorPlus="3.043488452417719" errorMinus="3.043488452417719"/> + </dataPoint> + <dataPoint> + <measurement value="22.4" errorPlus="2.6000000000000014" errorMinus="2.3999999999999986"/> + <measurement value="22.02" errorPlus="1.9193041589075974" errorMinus="1.9193041589075974"/> + </dataPoint> + <dataPoint> + <measurement value="27.41" errorPlus="2.59" errorMinus="2.41"/> + <measurement value="13.94" errorPlus="1.483748795989402" errorMinus="1.483748795989402"/> + </dataPoint> + <dataPoint> + <measurement value="32.42" errorPlus="2.5799999999999983" errorMinus="2.4200000000000017"/> + <measurement value="9.47" errorPlus="1.1535259955458308" errorMinus="1.1535259955458308"/> + </dataPoint> + <dataPoint> + <measurement value="37.42" errorPlus="2.5799999999999983" errorMinus="2.4200000000000017"/> + <measurement value="6.84" errorPlus="0.8705612681483137" errorMinus="0.8705612681483137"/> + </dataPoint> + <dataPoint> + <measurement value="44.7" errorPlus="5.299999999999997" errorMinus="4.700000000000003"/> + <measurement value="3.95" errorPlus="0.505871960084763" errorMinus="0.505871960084763"/> + </dataPoint> + <dataPoint> + <measurement value="54.72" errorPlus="5.280000000000001" errorMinus="4.719999999999999"/> + <measurement value="1.81" errorPlus="0.34182236556433815" errorMinus="0.34182236556433815"/> + </dataPoint> + <dataPoint> + <measurement value="64.77" errorPlus="5.230000000000004" errorMinus="4.769999999999996"/> + <measurement value="1.15" errorPlus="0.3303942493446277" errorMinus="0.3303942493446277"/> + </dataPoint> + <dataPoint> + <measurement value="74.79" errorPlus="5.209999999999994" errorMinus="4.790000000000006"/> + <measurement value="0.75" errorPlus="0.2785480209945854" errorMinus="0.2785480209945854"/> + </dataPoint> + <dataPoint> + <measurement value="89.21" errorPlus="10.790000000000006" errorMinus="9.209999999999994"/> + <measurement value="0.313" errorPlus="0.10932368446041324" errorMinus="0.10932368446041324"/> + </dataPoint> + <dataPoint> + <measurement value="109.27" errorPlus="10.730000000000004" errorMinus="9.269999999999996"/> + <measurement value="0.084" errorPlus="0.034331624138685894" errorMinus="0.034331624138685894"/> + </dataPoint> + <dataPoint> + <measurement value="137.4" errorPlus="22.599999999999994" errorMinus="17.400000000000006"/> + <measurement value="0.044" errorPlus="0.01854044487060653" errorMinus="0.01854044487060653"/> + </dataPoint> + <dataPoint> + <measurement value="177.64" errorPlus="22.360000000000014" errorMinus="17.639999999999986"/> + <measurement value="0.0077" errorPlus="0.00703738484381805" errorMinus="0.00703738484381805"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Thu Aug 12 23:09:30 2010 (r2653) +++ trunk/data/refdata/Makefile.am Mon Aug 16 11:57:11 2010 (r2654) @@ -75,6 +75,7 @@ D0_1996_S3214044.aida \ D0_1996_S3324664.aida \ D0_1998_S3711838.aida \ + D0_2000_S4480767.aida \ D0_2001_S4674421.aida \ D0_2004_S5992206.aida \ D0_2006_S6438750.aida \ Added: trunk/src/Analyses/D0_2000_S4480767.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/D0_2000_S4480767.cc Mon Aug 16 11:57:11 2010 (r2654) @@ -0,0 +1,80 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/WFinder.hh" + +namespace Rivet { + + + class D0_2000_S4480767 : public Analysis { + public: + + /// @name Constructors etc. + //@{ + + /// Constructor + D0_2000_S4480767() + : Analysis("D0_2000_S4480767") + { + setNeedsCrossSection(true); + } + + //@} + + + public: + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + WFinder wf(-5, 5, 0.0*GeV, ELECTRON, 0.0*GeV, 200.0*GeV, 0.0*GeV, 0.2); + addProjection(wf, "WFinder"); + + _h_W_pT = bookHistogram1D(1, 1, 1); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const double weight = event.weight(); + + const WFinder& wf = applyProjection<WFinder>(event, "WFinder"); + if (wf.size() == 0) vetoEvent; + + _h_W_pT->fill(wf.particles()[0].momentum().pT()/GeV, weight); + } + + + /// Normalise histograms etc., after the run + void finalize() { + + scale(_h_W_pT, crossSection()/sumOfWeights()); + + } + + //@} + + + private: + + /// @name Histograms + //@{ + + AIDA::IHistogram1D *_h_W_pT; + //@} + + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<D0_2000_S4480767> plugin_D0_2000_S4480767; + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Thu Aug 12 23:09:30 2010 (r2653) +++ trunk/src/Analyses/Makefile.am Mon Aug 16 11:57:11 2010 (r2654) @@ -78,6 +78,7 @@ lib_LTLIBRARIES += RivetD0Analyses.la RivetD0Analyses_la_SOURCES = \ + D0_2000_S4480767.cc \ D0_2001_S4674421.cc \ D0_2004_S5992206.cc \ D0_2006_S6438750.cc \
More information about the Rivet-svn mailing list |