|
[Rivet-svn] r4022 - in branches/2012-06-aidarivet: . data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Nov 15 15:17:03 GMT 2012
Author: hoeth Date: Thu Nov 15 15:17:03 2012 New Revision: 4022 Log: Added LHCF_2012_I1115479 analysis Added: branches/2012-06-aidarivet/data/anainfo/LHCF_2012_I1115479.info branches/2012-06-aidarivet/data/plotinfo/LHCF_2012_I1115479.plot branches/2012-06-aidarivet/data/refdata/LHCF_2012_I1115479.aida branches/2012-06-aidarivet/src/Analyses/LHCF_2012_I1115479.cc Modified: branches/2012-06-aidarivet/ChangeLog branches/2012-06-aidarivet/data/anainfo/Makefile.am branches/2012-06-aidarivet/data/plotinfo/Makefile.am branches/2012-06-aidarivet/data/refdata/Makefile.am branches/2012-06-aidarivet/src/Analyses/Makefile.am Modified: branches/2012-06-aidarivet/ChangeLog ============================================================================== --- branches/2012-06-aidarivet/ChangeLog Thu Nov 15 13:39:16 2012 (r4021) +++ branches/2012-06-aidarivet/ChangeLog Thu Nov 15 15:17:03 2012 (r4022) @@ -1,3 +1,7 @@ +2012-11-15 Hendrik Hoeth <hendrik.hoeth at cern.ch> + + * Added LHCF_2012_I1115479 analysis + 2012-11-14 Hendrik Hoeth <hendrik.hoeth at cern.ch> * Added CMS_2012_I1102908 analysis Added: branches/2012-06-aidarivet/data/anainfo/LHCF_2012_I1115479.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/anainfo/LHCF_2012_I1115479.info Thu Nov 15 15:17:03 2012 (r4022) @@ -0,0 +1,35 @@ +Name: LHCF_2012_I1115479 +Year: 2012 +Summary: Measurement of forward neutral pion transverse momentum spectra for $\sqrt{s}$ = 7 TeV proton-proton collisions at LHC +Experiment: LHCF +Collider: LHC +SpireID: 9602062 +InspireID: 1115479 +Status: VALIDATED +Authors: + - Sercan Sen <ssen at cern.ch> +References: + - arXiv:1205.4578 [hep-ex] +RunInfo: + Inelastic events (ND+SD+DD) at $\sqrt{s}$ = 7 TeV. +NumEvents: 500000 +Beams: [p+, p+] +Energies: [7000] +Description: + The inclusive production rate of neutral pions has been measured by + LHCf experiment during $\sqrt{s}=7$ TeV pp collision operation in early + 2010. In order to ensure good event reconstruction efficiency, the + range of the $\pi^0$ rapidity and $p_\perp$ are limited to $8.9 < y < 11.0$ + and $p_\perp < 0.6$ GeV, respectively. +BibKey: Adriani:2012ap +BibTeX: '@Article{Adriani:2012ap, + author = "Adriani, O. and others", + title = "{Measurement of forward neutral pion transverse momentum + spectra for $\sqrt{s}$ = 7TeV proton-proton collisions at + LHC}", + collaboration = "LHCf Collaboration", + year = "2012", + eprint = "1205.4578", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1205.4578;%%",}' Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/anainfo/Makefile.am Thu Nov 15 13:39:16 2012 (r4021) +++ branches/2012-06-aidarivet/data/anainfo/Makefile.am Thu Nov 15 15:17:03 2012 (r4022) @@ -163,6 +163,7 @@ LHCB_2010_S8758301.info \ LHCB_2011_I917009.info \ LHCB_2011_I919315.info \ + LHCF_2012_I1115479.info \ MC_DIJET.info \ MC_DIPHOTON.info \ MC_GENERIC.info \ Added: branches/2012-06-aidarivet/data/plotinfo/LHCF_2012_I1115479.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/plotinfo/LHCF_2012_I1115479.plot Thu Nov 15 15:17:03 2012 (r4022) @@ -0,0 +1,28 @@ +# BEGIN PLOT /LHCF_2012_I1115479/d0 +XLabel=$p_\perp$ [GeV] +YLabel=$1/\sigma_\text{inel}$ $E \mathrm{d}^3\sigma/\mathrm{d}p^3$ [GeV$^{-2}$] +# END PLOT + +# BEGIN PLOT /LHCF_2012_I1115479/d01-x01-y01 +Title=$\pi^0$ $p_\perp$ spectrum, $8.9 < y < 9.0$ +# END PLOT + +# BEGIN PLOT /LHCF_2012_I1115479/d02-x01-y01 +Title=$\pi^0$ $p_\perp$ spectrum, $9.0 < y < 9.2$ +# END PLOT + +# BEGIN PLOT /LHCF_2012_I1115479/d03-x01-y01 +Title=$\pi^0$ $p_\perp$ spectrum, $9.2 < y < 9.4$ +# END PLOT + +# BEGIN PLOT /LHCF_2012_I1115479/d04-x01-y01 +Title=$\pi^0$ $p_\perp$ spectrum, $9.4 < y < 9.6$ +# END PLOT + +# BEGIN PLOT /LHCF_2012_I1115479/d05-x01-y01 +Title=$\pi^0$ $p_\perp$ spectrum, $9.6 < y < 10.0$ +# END PLOT + +# BEGIN PLOT /LHCF_2012_I1115479/d06-x01-y01 +Title=$\pi^0$ $p_\perp$ spectrum, $10.0 < y < 11.0$ +# END PLOT Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/plotinfo/Makefile.am Thu Nov 15 13:39:16 2012 (r4021) +++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am Thu Nov 15 15:17:03 2012 (r4022) @@ -160,6 +160,7 @@ LHCB_2010_S8758301.plot \ LHCB_2011_I917009.plot \ LHCB_2011_I919315.plot \ + LHCF_2012_I1115479.plot \ MC_GENERIC.plot \ MC_IDENTIFIED.plot \ MC_LEADJETUE.plot \ Added: branches/2012-06-aidarivet/data/refdata/LHCF_2012_I1115479.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/refdata/LHCF_2012_I1115479.aida Thu Nov 15 15:17:03 2012 (r4022) @@ -0,0 +1,213 @@ +<?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="LWH"/> + <dataPointSet name="d01-x01-y01" dimension="2" path="/REF/LHCF_2012_I1115479" title="1/SIG E**D3(SIG)/DP**3" > + <dataPoint> + <measurement value="0.125" errorPlus="0.024999999999999994" errorMinus="0.024999999999999994"/> + <measurement value="0.271" errorPlus="0.158" errorMinus="0.141"/> + </dataPoint> + <dataPoint> + <measurement value="0.175" errorPlus="0.025000000000000022" errorMinus="0.024999999999999994"/> + <measurement value="0.195" errorPlus="0.0995" errorMinus="0.085"/> + </dataPoint> + <dataPoint> + <measurement value="0.225" errorPlus="0.024999999999999994" errorMinus="0.024999999999999994"/> + <measurement value="0.125" errorPlus="0.0566" errorMinus="0.0498"/> + </dataPoint> + <dataPoint> + <measurement value="0.275" errorPlus="0.024999999999999967" errorMinus="0.025000000000000022"/> + <measurement value="0.0715" errorPlus="0.029" errorMinus="0.0254"/> + </dataPoint> + <dataPoint> + <measurement value="0.32499999999999996" errorPlus="0.025000000000000022" errorMinus="0.024999999999999967"/> + <measurement value="0.0434" errorPlus="0.00422" errorMinus="0.00421"/> + </dataPoint> + <dataPoint> + <measurement value="0.375" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.0236" errorPlus="0.00366" errorMinus="0.00365"/> + </dataPoint> + <dataPoint> + <measurement value="0.42500000000000004" errorPlus="0.024999999999999967" errorMinus="0.025000000000000022"/> + <measurement value="0.015" errorPlus="0.00326" errorMinus="0.00325"/> + </dataPoint> + <dataPoint> + <measurement value="0.475" errorPlus="0.025000000000000022" errorMinus="0.024999999999999967"/> + <measurement value="0.00673" errorPlus="0.00216" errorMinus="0.00215"/> + </dataPoint> + <dataPoint> + <measurement value="0.55" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/> + <measurement value="0.00342" errorPlus="0.00127" errorMinus="0.00127"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/LHCF_2012_I1115479" title="1/SIG E**D3(SIG)/DP**3" > + <dataPoint> + <measurement value="0.125" errorPlus="0.024999999999999994" errorMinus="0.024999999999999994"/> + <measurement value="0.23" errorPlus="0.0812" errorMinus="0.0806"/> + </dataPoint> + <dataPoint> + <measurement value="0.175" errorPlus="0.025000000000000022" errorMinus="0.024999999999999994"/> + <measurement value="0.151" errorPlus="0.0419" errorMinus="0.0399"/> + </dataPoint> + <dataPoint> + <measurement value="0.225" errorPlus="0.024999999999999994" errorMinus="0.024999999999999994"/> + <measurement value="0.0892" errorPlus="0.00314" errorMinus="0.00314"/> + </dataPoint> + <dataPoint> + <measurement value="0.275" errorPlus="0.024999999999999967" errorMinus="0.025000000000000022"/> + <measurement value="0.0543" errorPlus="0.00335" errorMinus="0.00335"/> + </dataPoint> + <dataPoint> + <measurement value="0.32499999999999996" errorPlus="0.025000000000000022" errorMinus="0.024999999999999967"/> + <measurement value="0.0321" errorPlus="0.00332" errorMinus="0.00331"/> + </dataPoint> + <dataPoint> + <measurement value="0.375" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.018" errorPlus="0.00273" errorMinus="0.00272"/> + </dataPoint> + <dataPoint> + <measurement value="0.42500000000000004" errorPlus="0.024999999999999967" errorMinus="0.025000000000000022"/> + <measurement value="0.00891" errorPlus="0.00184" errorMinus="0.00183"/> + </dataPoint> + <dataPoint> + <measurement value="0.475" errorPlus="0.025000000000000022" errorMinus="0.024999999999999967"/> + <measurement value="0.00359" errorPlus="0.00101" errorMinus="0.00101"/> + </dataPoint> + <dataPoint> + <measurement value="0.55" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/> + <measurement value="9.72E-4" errorPlus="3.65E-4" errorMinus="3.66E-4"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y01" dimension="2" path="/REF/LHCF_2012_I1115479" title="1/SIG E**D3(SIG)/DP**3" > + <dataPoint> + <measurement value="0.025" errorPlus="0.025" errorMinus="0.025"/> + <measurement value="0.331" errorPlus="0.184" errorMinus="0.158"/> + </dataPoint> + <dataPoint> + <measurement value="0.07500000000000001" errorPlus="0.024999999999999994" errorMinus="0.02500000000000001"/> + <measurement value="0.231" errorPlus="0.107" errorMinus="0.0905"/> + </dataPoint> + <dataPoint> + <measurement value="0.125" errorPlus="0.024999999999999994" errorMinus="0.024999999999999994"/> + <measurement value="0.166" errorPlus="0.0196" errorMinus="0.0196"/> + </dataPoint> + <dataPoint> + <measurement value="0.175" errorPlus="0.025000000000000022" errorMinus="0.024999999999999994"/> + <measurement value="0.106" errorPlus="0.00376" errorMinus="0.00376"/> + </dataPoint> + <dataPoint> + <measurement value="0.225" errorPlus="0.024999999999999994" errorMinus="0.024999999999999994"/> + <measurement value="0.0571" errorPlus="0.00269" errorMinus="0.00269"/> + </dataPoint> + <dataPoint> + <measurement value="0.275" errorPlus="0.024999999999999967" errorMinus="0.025000000000000022"/> + <measurement value="0.0358" errorPlus="0.00298" errorMinus="0.00297"/> + </dataPoint> + <dataPoint> + <measurement value="0.32499999999999996" errorPlus="0.025000000000000022" errorMinus="0.024999999999999967"/> + <measurement value="0.0177" errorPlus="0.00229" errorMinus="0.00229"/> + </dataPoint> + <dataPoint> + <measurement value="0.375" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.00807" errorPlus="0.00149" errorMinus="0.00149"/> + </dataPoint> + <dataPoint> + <measurement value="0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="0.00135" errorPlus="3.72E-4" errorMinus="3.71E-4"/> + </dataPoint> + <dataPoint> + <measurement value="0.55" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/> + <measurement value="1.47E-4" errorPlus="9.18E-5" errorMinus="9.17E-5"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y01" dimension="2" path="/REF/LHCF_2012_I1115479" title="1/SIG E**D3(SIG)/DP**3" > + <dataPoint> + <measurement value="0.025" errorPlus="0.025" errorMinus="0.025"/> + <measurement value="0.203" errorPlus="0.0309" errorMinus="0.0309"/> + </dataPoint> + <dataPoint> + <measurement value="0.07500000000000001" errorPlus="0.024999999999999994" errorMinus="0.02500000000000001"/> + <measurement value="0.173" errorPlus="0.0164" errorMinus="0.0164"/> + </dataPoint> + <dataPoint> + <measurement value="0.125" errorPlus="0.024999999999999994" errorMinus="0.024999999999999994"/> + <measurement value="0.107" errorPlus="0.0046" errorMinus="0.00461"/> + </dataPoint> + <dataPoint> + <measurement value="0.175" errorPlus="0.025000000000000022" errorMinus="0.024999999999999994"/> + <measurement value="0.063" errorPlus="0.00245" errorMinus="0.00245"/> + </dataPoint> + <dataPoint> + <measurement value="0.225" errorPlus="0.024999999999999994" errorMinus="0.024999999999999994"/> + <measurement value="0.032" errorPlus="0.00264" errorMinus="0.00263"/> + </dataPoint> + <dataPoint> + <measurement value="0.275" errorPlus="0.024999999999999967" errorMinus="0.025000000000000022"/> + <measurement value="0.0145" errorPlus="0.00215" errorMinus="0.00214"/> + </dataPoint> + <dataPoint> + <measurement value="0.32499999999999996" errorPlus="0.025000000000000022" errorMinus="0.024999999999999967"/> + <measurement value="0.00364" errorPlus="9.29E-4" errorMinus="9.28E-4"/> + </dataPoint> + <dataPoint> + <measurement value="0.375" errorPlus="0.025000000000000022" errorMinus="0.025000000000000022"/> + <measurement value="0.00154" errorPlus="6.21E-4" errorMinus="6.2E-4"/> + </dataPoint> + <dataPoint> + <measurement value="0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="5.43E-5" errorPlus="6.41E-5" errorMinus="6.41E-5"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d05-x01-y01" dimension="2" path="/REF/LHCF_2012_I1115479" title="1/SIG E**D3(SIG)/DP**3" > + <dataPoint> + <measurement value="0.025" errorPlus="0.025" errorMinus="0.025"/> + <measurement value="0.12" errorPlus="0.00968" errorMinus="0.00966"/> + </dataPoint> + <dataPoint> + <measurement value="0.07500000000000001" errorPlus="0.024999999999999994" errorMinus="0.02500000000000001"/> + <measurement value="0.0828" errorPlus="0.0029" errorMinus="0.00289"/> + </dataPoint> + <dataPoint> + <measurement value="0.125" errorPlus="0.024999999999999994" errorMinus="0.024999999999999994"/> + <measurement value="0.0449" errorPlus="0.00188" errorMinus="0.00188"/> + </dataPoint> + <dataPoint> + <measurement value="0.175" errorPlus="0.025000000000000022" errorMinus="0.024999999999999994"/> + <measurement value="0.021" errorPlus="0.00128" errorMinus="0.00128"/> + </dataPoint> + <dataPoint> + <measurement value="0.225" errorPlus="0.024999999999999994" errorMinus="0.024999999999999994"/> + <measurement value="0.00743" errorPlus="9.76E-4" errorMinus="9.73E-4"/> + </dataPoint> + <dataPoint> + <measurement value="0.275" errorPlus="0.024999999999999967" errorMinus="0.025000000000000022"/> + <measurement value="0.00184" errorPlus="5.16E-4" errorMinus="5.15E-4"/> + </dataPoint> + <dataPoint> + <measurement value="0.35" errorPlus="0.050000000000000044" errorMinus="0.04999999999999999"/> + <measurement value="2.17E-4" errorPlus="1.33E-4" errorMinus="1.33E-4"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d06-x01-y01" dimension="2" path="/REF/LHCF_2012_I1115479" title="1/SIG E**D3(SIG)/DP**3" > + <dataPoint> + <measurement value="0.025" errorPlus="0.025" errorMinus="0.025"/> + <measurement value="0.0128" errorPlus="0.00142" errorMinus="0.00142"/> + </dataPoint> + <dataPoint> + <measurement value="0.07500000000000001" errorPlus="0.024999999999999994" errorMinus="0.02500000000000001"/> + <measurement value="0.00755" errorPlus="8.85E-4" errorMinus="8.88E-4"/> + </dataPoint> + <dataPoint> + <measurement value="0.125" errorPlus="0.024999999999999994" errorMinus="0.024999999999999994"/> + <measurement value="0.00237" errorPlus="3.76E-4" errorMinus="3.77E-4"/> + </dataPoint> + <dataPoint> + <measurement value="0.175" errorPlus="0.025000000000000022" errorMinus="0.024999999999999994"/> + <measurement value="1.91E-4" errorPlus="6.98E-5" errorMinus="6.99E-5"/> + </dataPoint> + <dataPoint> + <measurement value="0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/> + <measurement value="8.37E-6" errorPlus="1.05E-5" errorMinus="1.05E-5"/> + </dataPoint> + </dataPointSet> +</aida> Modified: branches/2012-06-aidarivet/data/refdata/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/refdata/Makefile.am Thu Nov 15 13:39:16 2012 (r4021) +++ branches/2012-06-aidarivet/data/refdata/Makefile.am Thu Nov 15 15:17:03 2012 (r4022) @@ -2,6 +2,7 @@ LHCB_2010_S8758301.aida \ LHCB_2011_I917009.aida \ LHCB_2011_I919315.aida \ + LHCF_2012_I1115479.aida \ JADE_OPAL_2000_S4300807.aida \ JADE_1998_S3612880.aida \ TASSO_1990_S2148048.aida \ Added: branches/2012-06-aidarivet/src/Analyses/LHCF_2012_I1115479.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/src/Analyses/LHCF_2012_I1115479.cc Thu Nov 15 15:17:03 2012 (r4022) @@ -0,0 +1,66 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Projections/UnstableFinalState.hh" +#include "Rivet/Tools/BinnedHistogram.hh" + +namespace Rivet { + + + class LHCF_2012_I1115479 : public Analysis { + public: + + LHCF_2012_I1115479() + : Analysis("LHCF_2012_I1115479") + { } + + + public: + + void init() { + addProjection(UnstableFinalState(),"UFS"); + + _binnedHistos_y_pT.addHistogram( 8.9, 9.0, bookHistogram1D(1, 1, 1)); + _binnedHistos_y_pT.addHistogram( 9.0, 9.2, bookHistogram1D(2, 1, 1)); + _binnedHistos_y_pT.addHistogram( 9.2, 9.4, bookHistogram1D(3, 1, 1)); + _binnedHistos_y_pT.addHistogram( 9.4, 9.6, bookHistogram1D(4, 1, 1)); + _binnedHistos_y_pT.addHistogram( 9.6, 10.0, bookHistogram1D(5, 1, 1)); + _binnedHistos_y_pT.addHistogram(10.0, 11.0, bookHistogram1D(6, 1, 1)); + } + + + void analyze(const Event& event) { + const UnstableFinalState& ufs = applyProjection<UnstableFinalState>(event, "UFS"); + + const double weight = event.weight(); + const double dphi = TWOPI; + + foreach (const Particle& p, ufs.particles()) { + if (p.pdgId() == 111) { + double pT = p.momentum().pT(); + double y = p.momentum().rapidity(); + + if (pT > 0.6*GeV) continue; + + const double scaled_weight = weight/(dphi*pT/GeV); + _binnedHistos_y_pT.fill(y, pT/GeV, scaled_weight); + } + } + } + + + void finalize() { + _binnedHistos_y_pT.scale( 1./sumOfWeights() , this); + } + + private: + + BinnedHistogram<double> _binnedHistos_y_pT; + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(LHCF_2012_I1115479); + +} Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/src/Analyses/Makefile.am Thu Nov 15 13:39:16 2012 (r4021) +++ branches/2012-06-aidarivet/src/Analyses/Makefile.am Thu Nov 15 15:17:03 2012 (r4022) @@ -140,6 +140,11 @@ endif +lib_LTLIBRARIES += RivetLHCfAnalyses.la +RivetLHCfAnalyses_la_SOURCES = \ + LHCF_2012_I1115479.cc + + lib_LTLIBRARIES += RivetTOTEMAnalyses.la RivetTOTEMAnalyses_la_SOURCES = \ TOTEM_2012_I1115294.cc
More information about the Rivet-svn mailing list |