[Rivet] [Rivet-svn] r2212 - in trunk: data/anainfo data/plotinfo data/refdata src/Analyses

Andy Buckley andy.buckley at ed.ac.uk
Sun Jan 17 20:07:32 GMT 2010


Looks nice!

Hendrik, can you please add your analysis validation scripts to a
suitable area of the SVN repo, and then we can all start providing
similar scripts for other analyses. Including this one and any future
ones... best to do it while the memory is fresh!

Andy


On 17/01/10 15:34, blackhole at projects.hepforge.org wrote:
> 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 \
> _______________________________________________
> Rivet-svn mailing list
> Rivet-svn at projects.hepforge.org
> http://www.hepforge.org/lists/listinfo/rivet-svn


-- 
Dr Andy Buckley
SUPA Advanced Research Fellow
Particle Physics Experiment Group, University of Edinburgh


More information about the Rivet mailing list