|
[Rivet-svn] r1962 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgMon Oct 26 17:04:07 GMT 2009
Author: fsiegert Date: Mon Oct 26 17:04:07 2009 New Revision: 1962 Log: Add CDF_2009_S8436959. Compiles, but everything else remains to be tested. Added: trunk/data/anainfo/CDF_2009_S8436959.info trunk/data/plotinfo/CDF_2009_S8436959.plot trunk/data/refdata/CDF_2009_S8436959.aida trunk/src/Analyses/CDF_2009_S8436959.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_2009_S8436959.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/CDF_2009_S8436959.info Mon Oct 26 17:04:07 2009 (r1962) @@ -0,0 +1,22 @@ +Name: CDF_2009_S8436959 +Year: 2009 +Summary: Measurement of the Inclusive Isolated Prompt Photon Cross Section +Experiment: CDF +Collider: Tevatron Run 2 +SpiresID: 8436959 +Status: UNVALIDATED +Authors: + - Frank Siegert <frank.siegert at durham.ac.uk> +References: + - arXiv:0910.3623 [hep-ex] +RunInfo: + $\gamma$ + jet processes in ppbar collisions at $\sqrt{s} = 1960$ GeV. + Minimum pT cut on the photon in the analysis is 30 GeV. +NumEvents: 1000000 +PtCuts: [30] +Description: + A measurement of the cross section for the inclusive production of isolated + photons. The measurement covers the pseudorapidity region $|\eta^\gamma|<1.0$ + and the transverse energy range $E_T^\gamma>30$ GeV and is based on 2.5/fb of + integrated luminosity. The cross section is measured differential in + $E_\perp(\gamma)$. Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Mon Oct 26 16:15:34 2009 (r1961) +++ trunk/data/anainfo/Makefile.am Mon Oct 26 17:04:07 2009 (r1962) @@ -26,6 +26,7 @@ CDF_2008_S8095620.info \ CDF_2009_S8233977.info \ CDF_2009_S8383952.info \ + CDF_2009_S8436959.info \ D0_1996_S3214044.info \ D0_1996_S3324664.info \ D0_2001_S4674421.info \ Added: trunk/data/plotinfo/CDF_2009_S8436959.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/CDF_2009_S8436959.plot Mon Oct 26 17:04:07 2009 (r1962) @@ -0,0 +1,5 @@ +# BEGIN PLOT /CDF_2009_S8436959/d01-x01-y01 +Title=Transverse energy of isolated prompt photon +XLabel=$E_\perp^\gamma$ [GeV] +YLabel=$\mathrm{d}\sigma/\mathrm{d}\E_\perp^\gamma \mathrm{d}\eta^\gamma$ [pb/GeV] +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Mon Oct 26 16:15:34 2009 (r1961) +++ trunk/data/plotinfo/Makefile.am Mon Oct 26 17:04:07 2009 (r1962) @@ -26,6 +26,7 @@ CDF_2008_S8095620.plot \ CDF_2009_S8233977.plot \ CDF_2009_S8383952.plot \ + CDF_2009_S8436959.plot \ D0_1996_S3214044.plot \ D0_1996_S3324664.plot \ D0_2001_S4674421.plot \ Added: trunk/data/refdata/CDF_2009_S8436959.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/CDF_2009_S8436959.aida Mon Oct 26 17:04:07 2009 (r1962) @@ -0,0 +1,71 @@ +<?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/CDF_2009_S8436959" title="D(SIG)/DET(P=3)/DETARAP(P=3) IN PB/GEV" > + <dataPoint> + <measurement value="32.0" errorPlus="2.0" errorMinus="2.0"/> + <measurement value="123.0" errorPlus="17.86301276380891" errorMinus="17.86301276380891"/> + </dataPoint> + <dataPoint> + <measurement value="36.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="62.1" errorPlus="6.093189775478851" errorMinus="6.093189775478851"/> + </dataPoint> + <dataPoint> + <measurement value="41.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="31.0" errorPlus="2.611669198041743" errorMinus="2.611669198041743"/> + </dataPoint> + <dataPoint> + <measurement value="47.0" errorPlus="3.0" errorMinus="3.0"/> + <measurement value="17.2" errorPlus="1.4074822343461388" errorMinus="1.4074822343461388"/> + </dataPoint> + <dataPoint> + <measurement value="55.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="7.93" errorPlus="0.6709067404639784" errorMinus="0.6709067404639784"/> + </dataPoint> + <dataPoint> + <measurement value="65.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="3.54" errorPlus="0.3050259169316601" errorMinus="0.3050259169316601"/> + </dataPoint> + <dataPoint> + <measurement value="75.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.76" errorPlus="0.16294546817877445" errorMinus="0.16294546817877445"/> + </dataPoint> + <dataPoint> + <measurement value="85.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.908" errorPlus="0.07308542826035845" errorMinus="0.07308542826035845"/> + </dataPoint> + <dataPoint> + <measurement value="100.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.441" errorPlus="0.038691429141348606" errorMinus="0.038691429141348606"/> + </dataPoint> + <dataPoint> + <measurement value="120.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.168" errorPlus="0.014920705613341481" errorMinus="0.014920705613341481"/> + </dataPoint> + <dataPoint> + <measurement value="140.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.0725" errorPlus="0.0060166435825965295" errorMinus="0.0060166435825965295"/> + </dataPoint> + <dataPoint> + <measurement value="160.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.0341" errorPlus="0.0035025847598594956" errorMinus="0.0035025847598594956"/> + </dataPoint> + <dataPoint> + <measurement value="185.0" errorPlus="15.0" errorMinus="15.0"/> + <measurement value="0.0146" errorPlus="0.0013875078234013673" errorMinus="0.0013875078234013673"/> + </dataPoint> + <dataPoint> + <measurement value="215.0" errorPlus="15.0" errorMinus="15.0"/> + <measurement value="0.00566" errorPlus="6.461826379592692E-4" errorMinus="6.461826379592692E-4"/> + </dataPoint> + <dataPoint> + <measurement value="265.0" errorPlus="35.0" errorMinus="35.0"/> + <measurement value="0.00138" errorPlus="1.679389043670346E-4" errorMinus="1.679389043670346E-4"/> + </dataPoint> + <dataPoint> + <measurement value="350.0" errorPlus="50.0" errorMinus="50.0"/> + <measurement value="1.49E-4" errorPlus="2.8976562183944456E-5" errorMinus="2.8976562183944456E-5"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Mon Oct 26 16:15:34 2009 (r1961) +++ trunk/data/refdata/Makefile.am Mon Oct 26 17:04:07 2009 (r1962) @@ -54,6 +54,7 @@ CDF_2008_LEADINGJETS.aida \ CDF_2009_S8233977.aida \ CDF_2009_S8383952.aida \ + CDF_2009_S8436959.aida \ D0_1996_S3214044.aida \ D0_1996_S3324664.aida \ D0_2001_S4674421.aida \ Added: trunk/src/Analyses/CDF_2009_S8436959.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/CDF_2009_S8436959.cc Mon Oct 26 17:04:07 2009 (r1962) @@ -0,0 +1,95 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/IdentifiedFinalState.hh" + +namespace Rivet { + + + class CDF_2009_S8436959 : public Analysis { + public: + + /// @name Constructors etc. + //@{ + + /// Constructor + CDF_2009_S8436959() + : Analysis("CDF_2009_S8436959") + { + setBeams(PROTON, ANTIPROTON); + setNeedsCrossSection(true); + } + + //@} + + + public: + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + FinalState fs; + addProjection(fs, "FS"); + + IdentifiedFinalState ifs(-1.0, 1.0, 30.0*GeV); + ifs.acceptId(PHOTON); + addProjection(ifs, "IFS"); + + _h_Et_photon = bookHistogram1D(1, 1, 1); + + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const double weight = event.weight(); + + ParticleVector photons; + ParticleVector fs = applyProjection<IdentifiedFinalState>(event, "FS").particles(); + foreach (const Particle& photon, applyProjection<IdentifiedFinalState>(event, "IFS").particles()) { + FourMomentum mom_in_cone; + foreach (const Particle& p, fs) { + if (deltaR(p.momentum(), photon.momentum()) < 0.4) { + mom_in_cone += p.momentum(); + } + } + if (mom_in_cone.Et()-photon.momentum().Et() < 2.0*GeV) { + photons.push_back(photon); + } + } + if (photons.size() != 1) { + vetoEvent; + } + + _h_Et_photon->fill(photons[0].momentum().Et(), weight); + } + + + /// Normalise histograms etc., after the run + void finalize() { + scale(_h_Et_photon, crossSection()/sumOfWeights()/2.0); + } + + //@} + + + private: + + /// @name Histograms + //@{ + AIDA::IHistogram1D *_h_Et_photon; + //@} + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<CDF_2009_S8436959> plugin_CDF_2009_S8436959; + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Mon Oct 26 16:15:34 2009 (r1961) +++ trunk/src/Analyses/Makefile.am Mon Oct 26 17:04:07 2009 (r1962) @@ -34,7 +34,8 @@ CDF_2008_S8093652.cc \ CDF_2008_S8095620.cc \ CDF_2009_S8233977.cc \ - CDF_2009_S8383952.cc + CDF_2009_S8383952.cc \ + CDF_2009_S8436959.cc lib_LTLIBRARIES += RivetD0Analyses.la
More information about the Rivet-svn mailing list |