|
[Rivet-svn] r2212 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgSun Jan 17 15:34:03 GMT 2010
Author: holsch Date: Sun Jan 17 15:34:02 2010 New Revision: 2212 Log: Adding CDF run 1 measurement of W-pT in electronic W decays. Plots (scaled to area of refhisto) can be found here: http://users.hepforge.org/~holsch/Rivet/Validation/CDF_1991_S2313472-check. I used the fpythia-D0_2008_S7837160.params from AGILe for W production and electronic decays. Added: trunk/data/anainfo/CDF_1991_S2313472.info trunk/data/plotinfo/CDF_1991_S2313472.plot trunk/data/refdata/CDF_1991_S2313472.aida trunk/src/Analyses/CDF_1991_S2313472.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/CDF_1991_S2313472.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/CDF_1991_S2313472.info Sun Jan 17 15:34:02 2010 (r2212) @@ -0,0 +1,22 @@ +Name: CDF_1991_S2313472 +Year: 1991 +Summary: W-boson pT measurement in $p\bar{p}$ collisions at $\sqrt{s}=1.8~\TeV$ +Experiment: CDF +Collider: Tevatron +SpiresID: 2313472 +Status: UNVALIDATED +Authors: + - Holger Schulz <hschulz at physik.hu-berlin.de> +References: + - Phys.Rev.Lett.66:2951-2955,1991 +RunInfo: + QCD events with W+- production and electronic decays +NumEvents: The measurement had 2500 events. Maybe 100k MC events are enough. +PtCuts: +Description: + This is a CDF analysis from run 1, where the distribution of the transverse + momentum of W candidates that decay electronically, is measured. The + electron is required to be within $\left|\eta\right| < 1.1$, to have + a transverse energy of $E_\perp > 20~\GeV$ and a $p_\perp > 12~\GeV$. + The neutrino is required to produce a missing energy of + $E_{\perp, \text{ miss}}>20~\GeV$. Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Sat Jan 16 21:57:06 2010 (r2211) +++ trunk/data/anainfo/Makefile.am Sun Jan 17 15:34:02 2010 (r2212) @@ -4,6 +4,7 @@ ALEPH_1996_S3196992.info \ CDF_1988_S1865951.info \ CDF_1990_S2089246.info \ + CDF_1991_S2313472.info \ CDF_1994_S2952106.info \ CDF_1996_S3108457.info \ CDF_1996_S3349578.info \ Added: trunk/data/plotinfo/CDF_1991_S2313472.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/CDF_1991_S2313472.plot Sun Jan 17 15:34:02 2010 (r2212) @@ -0,0 +1,5 @@ +# BEGIN PLOT /CDF_1991_S2313472/d01-x01-y01 +Title=$p_\perp$ of W boson in $\PW \to e\nu$ decays +XLabel=$p_\perp(\PW)$ / GeV +YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{p_\perp(\PW)}$ +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Sat Jan 16 21:57:06 2010 (r2211) +++ trunk/data/plotinfo/Makefile.am Sun Jan 17 15:34:02 2010 (r2212) @@ -4,6 +4,7 @@ ALEPH_1996_S3196992.plot \ CDF_1988_S1865951.plot \ CDF_1990_S2089246.plot \ + CDF_1991_S2313472.plot \ CDF_1994_S2952106.plot \ CDF_1996_S3108457.plot \ CDF_1996_S3349578.plot \ Added: trunk/data/refdata/CDF_1991_S2313472.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/CDF_1991_S2313472.aida Sun Jan 17 15:34:02 2010 (r2212) @@ -0,0 +1,95 @@ +<?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.0" package="HepData"/> + <dataPointSet name="d01-x01-y01" dimension="2" + path="/REF/CDF_1991_S2313472" title="D(SIG)/D(PT) IN PB/GEV"> + <dimension dim="0" title="None" /> + <dimension dim="1" title="None" /> + <annotation> + <item key="Title" value="D(SIG)/DPT" sticky="true"/> + <item key="AidaPath" value="/REF/CDF_1991_S2313472/d01-x01-y01" sticky="true"/> + <item key="FullPath" value="/REF/CDF_1991_S2313472/d01-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement value="1.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="694.0" errorPlus="457.19361325372864" errorMinus="457.19361325372864"/> + </dataPoint> + <dataPoint> + <measurement value="3.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="1562.0" errorPlus="684.6407817242557" errorMinus="684.6407817242557"/> + </dataPoint> + <dataPoint> + <measurement value="5.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="1419.0" errorPlus="391.1265779770022" errorMinus="391.1265779770022"/> + </dataPoint> + <dataPoint> + <measurement value="7.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="1084.0" errorPlus="229.31419493786248" errorMinus="229.31419493786248"/> + </dataPoint> + <dataPoint> + <measurement value="9.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="963.0" errorPlus="205.27298896834918" errorMinus="205.27298896834918"/> + </dataPoint> + <dataPoint> + <measurement value="11.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="762.0" errorPlus="181.23189564753773" errorMinus="181.23189564753773"/> + </dataPoint> + <dataPoint> + <measurement value="13.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="684.0" errorPlus="171.74690681348528" errorMinus="171.74690681348528"/> + </dataPoint> + <dataPoint> + <measurement value="15.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="521.0" errorPlus="130.034610777285" errorMinus="130.034610777285"/> + </dataPoint> + <dataPoint> + <measurement value="17.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="451.0" errorPlus="112.53888216967503" errorMinus="112.53888216967503"/> + </dataPoint> + <dataPoint> + <measurement value="19.0" errorPlus="1.75" errorMinus="1.0"/> + <measurement value="388.0" errorPlus="94.847245611035" errorMinus="94.847245611035"/> + </dataPoint> + <dataPoint> + <measurement value="22.5" errorPlus="2.5" errorMinus="1.75"/> + <measurement value="291.0" errorPlus="68.62215385719105" errorMinus="68.62215385719105"/> + </dataPoint> + <dataPoint> + <measurement value="27.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="154.0" errorPlus="36.6742416417845" errorMinus="36.6742416417845"/> + </dataPoint> + <dataPoint> + <measurement value="32.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="115.0" errorPlus="27.784887978899608" errorMinus="27.784887978899608"/> + </dataPoint> + <dataPoint> + <measurement value="37.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="61.1" errorPlus="16.18177987738061" errorMinus="16.18177987738061"/> + </dataPoint> + <dataPoint> + <measurement value="42.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="51.5" errorPlus="14.035668847618199" errorMinus="14.035668847618199"/> + </dataPoint> + <dataPoint> + <measurement value="47.5" errorPlus="3.6000000000000014" errorMinus="2.5"/> + <measurement value="40.4" errorPlus="11.597413504743201" errorMinus="11.597413504743201"/> + </dataPoint> + <dataPoint> + <measurement value="54.7" errorPlus="7.100000000000001" errorMinus="3.6000000000000014"/> + <measurement value="19.6" errorPlus="5.946427498927402" errorMinus="5.946427498927402"/> + </dataPoint> + <dataPoint> + <measurement value="68.9" errorPlus="15.349999999999994" errorMinus="7.100000000000001"/> + <measurement value="7.3" errorPlus="2.701851217221259" errorMinus="2.701851217221259"/> + </dataPoint> + <dataPoint> + <measurement value="99.6" errorPlus="25.799999999999997" errorMinus="15.349999999999994"/> + <measurement value="1.18" errorPlus="0.6860029154456998" errorMinus="0.6860029154456998"/> + </dataPoint> + <dataPoint> + <measurement value="151.2" errorPlus="25.799999999999997" errorMinus="25.799999999999997"/> + <measurement value="0.44" errorPlus="0.2830194339616981" errorMinus="0.2830194339616981"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Sat Jan 16 21:57:06 2010 (r2211) +++ trunk/data/refdata/Makefile.am Sun Jan 17 15:34:02 2010 (r2212) @@ -33,6 +33,7 @@ UA5_1988_S1867512.aida \ CDF_1988_S1865951.aida \ CDF_1990_S2089246.aida \ + CDF_1991_S2313472.aida \ CDF_1994_S2952106.aida \ CDF_1996_S3108457.aida \ CDF_1996_S3349578.aida \ Added: trunk/src/Analyses/CDF_1991_S2313472.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/CDF_1991_S2313472.cc Sun Jan 17 15:34:02 2010 (r2212) @@ -0,0 +1,105 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/ChargedFinalState.hh" +#include "Rivet/Projections/WFinder.hh" +#include "Rivet/Tools/ParticleIdUtils.hh" +/// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder... + +namespace Rivet { + + + class CDF_1991_S2313472 : public Analysis { + public: + + /// @name Constructors etc. + //@{ + + /// Constructor + CDF_1991_S2313472() + : Analysis("CDF_1991_S2313472") + { + /// @todo Set approriate for your analysis + setBeams(PROTON, ANTIPROTON); + + /// @todo Set whether your finalize method needs the generator cross section + setNeedsCrossSection(false); + } + + //@} + + + public: + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + WFinder wfe(-5, 5, 0.0*GeV, ELECTRON, 60.0*GeV, 100.0*GeV, 0.2); + addProjection(wfe, "WFe"); + + // Book histogram + _hist_wpt = 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, "WFe"); + if (wf.size() == 0) { + getLog() << Log::DEBUG << "No W candidates found: vetoing" << endl; + vetoEvent; + } + + // Require the electron to have ET > 12 GeV, pT > 6 GeV and |eta| < 1.1 + FourMomentum p_e; + int chg_e = 0; + + foreach (const Particle& l, wf.constituentsFinalState().particles()) { + const FourMomentum pl = l.momentum(); + if (abs(l.pdgId()) == ELECTRON) { + chg_e = PID::threeCharge(l.pdgId()); + p_e = pl; + const double eta_e = fabs(p_e.pseudorapidity()); + if ( (pl.Et()/GeV < 12.0) || (pl.pT()/GeV < 6.0) || (eta_e > 1.1) ) { + getLog() << Log::DEBUG << l.pdgId() << " ET,pT,eta:" << pl.Et()/GeV << "," << pl.pT()/GeV << "," << eta_e << " fails electron cut" << endl; + vetoEvent; + } + } + } + assert(chg_e != 0); + + FourMomentum pW = wf.particles()[0].momentum(); + getLog() << Log::DEBUG << "Dilepton mass = " << pW.mass()/GeV << " GeV" << endl; + getLog() << Log::DEBUG << "Dilepton pT = " << pW.pT()/GeV << " GeV" << endl; + _hist_wpt->fill(pW.pT()/GeV, weight); + } + + + /// Normalise histograms etc., after the run + void finalize() { + + + } + + //@} + + + private: + + /// @name Histograms + AIDA::IHistogram1D *_hist_wpt; + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<CDF_1991_S2313472> plugin_CDF_1991_S2313472; + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Sat Jan 16 21:57:06 2010 (r2211) +++ trunk/src/Analyses/Makefile.am Sun Jan 17 15:34:02 2010 (r2212) @@ -11,6 +11,7 @@ if ENABLE_UNVALIDATED lib_LTLIBRARIES += RivetUnvalidated.la RivetUnvalidated_la_SOURCES = \ + CDF_1991_S2313472.cc \ CDF_1996_S3108457.cc \ CDF_1996_S3349578.cc \ CDF_1996_S3418421.cc \
More information about the Rivet-svn mailing list |