|
[Rivet-svn] r2400 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgTue Apr 13 12:42:12 BST 2010
Author: fsiegert Date: Tue Apr 13 12:42:10 2010 New Revision: 2400 Log: Add first unvalidated version of CDF_1993_S2742446. Added: trunk/data/anainfo/CDF_1993_S2742446.info trunk/data/plotinfo/CDF_1993_S2742446.plot trunk/data/refdata/CDF_1993_S2742446.aida trunk/src/Analyses/CDF_1993_S2742446.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_1993_S2742446.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/CDF_1993_S2742446.info Tue Apr 13 12:42:10 2010 (r2400) @@ -0,0 +1,40 @@ +Name: CDF_1993_S2742446 +Year: 1993 +Summary: Angular distribution of prompt photon +Experiment: CDF +Collider: Tevatron Run 1 +SpiresID: 2742446 +Status: UNVALIDATED +Authors: + - Frank Siegert <frank.siegert at durham.ac.uk> +References: + - Phys.Rev.Lett.71:679-683,1993 + - doi:10.1103/PhysRevLett.71.679 +RunInfo: + All prompt photon production processes in $p \bar{p}$ at 1800 GeV. + Hadronisation should be switched off, because non-prompt photon production + has been corrected for. +NumEvents: 10000000 +Beams: [p-, p+] +Energies: [1800] +PtCuts: [15] +Description: + Data taken with the Collider Detector at Fermilab (CDF) during the 1988-1989 + run of the Tevatron are used to measure the distribution of the center-of-mass + (rest frame of the initial state partons) angle between isolated prompt + photons and the beam direction. +BibKey: Abe:1993cv +BibTeX: '@Article{Abe:1993cv, + author = "Abe, F. and others", + collaboration = "CDF", + title = "{The Center-of-mass angular distribution of prompt photons + produced in $p\bar{p}$ collisions at $\sqrt{s} = 1.8$ + TeV}", + journal = "Phys. Rev. Lett.", + volume = "71", + year = "1993", + pages = "679-683", + doi = "10.1103/PhysRevLett.71.679", + SLACcitation = "%%CITATION = PRLTA,71,679;%%" +}' +ToDo: Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Tue Apr 13 10:23:00 2010 (r2399) +++ trunk/data/anainfo/Makefile.am Tue Apr 13 12:42:10 2010 (r2400) @@ -8,6 +8,7 @@ CDF_1988_S1865951.info \ CDF_1990_S2089246.info \ CDF_1991_S2313472.info \ + CDF_1993_S2742446.info \ CDF_1994_S2952106.info \ CDF_1996_S3108457.info \ CDF_1996_S3349578.info \ Added: trunk/data/plotinfo/CDF_1993_S2742446.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/CDF_1993_S2742446.plot Tue Apr 13 12:42:10 2010 (r2400) @@ -0,0 +1,6 @@ +# BEGIN PLOT /CDF_1993_S2742446/d01-x01-y01 +Title=Angular distribution of prompt photons +XLabel=$|\cos\theta^*|=|\tanh((\eta_j-\eta_\gamma)/2.0)|$ +YLabel=$\mathrm{d}N/\mathrm{d}\cos\theta^*$ +LogY=0 +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Tue Apr 13 10:23:00 2010 (r2399) +++ trunk/data/plotinfo/Makefile.am Tue Apr 13 12:42:10 2010 (r2400) @@ -7,6 +7,7 @@ CDF_1988_S1865951.plot \ CDF_1990_S2089246.plot \ CDF_1991_S2313472.plot \ + CDF_1993_S2742446.plot \ CDF_1994_S2952106.plot \ CDF_1996_S3108457.plot \ CDF_1996_S3349578.plot \ Added: trunk/data/refdata/CDF_1993_S2742446.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/CDF_1993_S2742446.aida Tue Apr 13 12:42:10 2010 (r2400) @@ -0,0 +1,39 @@ +<?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_1993_S2742446" title="D(N)/DCOS(THETA(P=3,RF=CM))" > + <dataPoint> + <measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/> + <measurement value="1.086" errorPlus="0.2659661632614194" errorMinus="0.2659661632614194"/> + </dataPoint> + <dataPoint> + <measurement value="0.15000000000000002" errorPlus="0.04999999999999999" errorMinus="0.05000000000000002"/> + <measurement value="0.964" errorPlus="0.250569351677335" errorMinus="0.250569351677335"/> + </dataPoint> + <dataPoint> + <measurement value="0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="0.949" errorPlus="0.23992707225321613" errorMinus="0.23992707225321613"/> + </dataPoint> + <dataPoint> + <measurement value="0.35" errorPlus="0.050000000000000044" errorMinus="0.04999999999999999"/> + <measurement value="0.915" errorPlus="0.23787601812709075" errorMinus="0.23787601812709075"/> + </dataPoint> + <dataPoint> + <measurement value="0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="0.947" errorPlus="0.3097256850827842" errorMinus="0.3097256850827842"/> + </dataPoint> + <dataPoint> + <measurement value="0.55" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/> + <measurement value="1.838" errorPlus="0.6919609815589315" errorMinus="0.6919609815589315"/> + </dataPoint> + <dataPoint> + <measurement value="0.6499999999999999" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/> + <measurement value="2.782" errorPlus="0.9687543548289216" errorMinus="0.9687543548289216"/> + </dataPoint> + <dataPoint> + <measurement value="0.75" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/> + <measurement value="4.79" errorPlus="1.6587914275158284" errorMinus="1.6587914275158284"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Tue Apr 13 10:23:00 2010 (r2399) +++ trunk/data/refdata/Makefile.am Tue Apr 13 12:42:10 2010 (r2400) @@ -37,6 +37,7 @@ CDF_1988_S1865951.aida \ CDF_1990_S2089246.aida \ CDF_1991_S2313472.aida \ + CDF_1993_S2742446.aida \ CDF_1994_S2952106.aida \ CDF_1996_S3108457.aida \ CDF_1996_S3349578.aida \ Added: trunk/src/Analyses/CDF_1993_S2742446.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/CDF_1993_S2742446.cc Tue Apr 13 12:42:10 2010 (r2400) @@ -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/FastJets.hh" +#include "Rivet/Projections/LeadingParticlesFinalState.hh" +#include "Rivet/Projections/VetoedFinalState.hh" +#include "Rivet/Tools/ParticleIdUtils.hh" + +namespace Rivet { + + + class CDF_1993_S2742446 : public Analysis { + public: + + CDF_1993_S2742446() + : Analysis("CDF_1993_S2742446") + { + setBeams(PROTON, ANTIPROTON); + setNeedsCrossSection(false); + } + + + public: + + void init() { + + // The photon selection has been corrected to pTmin=22 GeV (vs. 23 in the trigger) + LeadingParticlesFinalState photonfs(FinalState(-0.9, 0.9, 22.0*GeV)); + photonfs.addParticleId(PHOTON); + addProjection(photonfs, "LeadingPhoton"); + + // FS excluding the leading photon + VetoedFinalState vfs(FinalState(-4.2, 4.2)); + vfs.addVetoOnThisFinalState(photonfs); + addProjection(vfs, "VFS"); + + // Jets + addProjection(FastJets(vfs, FastJets::CDFJETCLU, 0.7), "Jets"); + + _h_costheta = bookHistogram1D(1, 1, 1); + + } + + + void analyze(const Event& event) { + + const double weight = event.weight(); + + ParticleVector photons = applyProjection<LeadingParticlesFinalState>(event, "LeadingPhoton").particles(); + if (photons.size()!=1 || photons[0].momentum().pT()>45.0*GeV) { + vetoEvent; + } + FourMomentum leadingPhoton = photons[0].momentum(); + double eta_P = leadingPhoton.eta(); + double phi_P = leadingPhoton.phi(); + + // photon isolation: less than 2 GeV EM E_T + double Etsum=0.0; + foreach (const Particle& p, applyProjection<VetoedFinalState>(event, "VFS").particles()) { + if (PID::threeCharge(p.pdgId())!=0 && + deltaR(eta_P, phi_P, p.momentum().eta(), p.momentum().phi()) < 0.7) { + Etsum += p.momentum().Et(); + } + } + if (Etsum > 2.0*GeV) { + vetoEvent; + } + + // sum all jets in the opposite hemisphere in phi from the photon + FourMomentum jetsum; + foreach (const Jet& jet, applyProjection<FastJets>(event, "Jets").jets(10.0*GeV)) { + if (jet.momentum().phi()-phi_P > M_PI) { + jetsum+=jet.momentum(); + } + } + + double costheta = fabs(tanh((jetsum.eta()-eta_P)/2.0)); + + _h_costheta->fill(costheta, weight); + + } + + + void finalize() { + + normalize(_h_costheta, 1.4271); // fixed norm ok + + } + + + private: + + AIDA::IHistogram1D *_h_costheta; + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<CDF_1993_S2742446> plugin_CDF_1993_S2742446; + + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Tue Apr 13 10:23:00 2010 (r2399) +++ trunk/src/Analyses/Makefile.am Tue Apr 13 12:42:10 2010 (r2400) @@ -12,6 +12,7 @@ lib_LTLIBRARIES += RivetUnvalidated.la RivetUnvalidated_la_SOURCES = \ CDF_1991_S2313472.cc \ + CDF_1993_S2742446.cc \ CDF_1994_S2952106.cc \ CDF_1996_S3108457.cc \ CDF_1996_S3349578.cc \
More information about the Rivet-svn mailing list |