|
[Rivet-svn] r2476 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Jun 10 08:40:21 BST 2010
Author: holsch Date: Thu Jun 10 08:40:21 2010 New Revision: 2476 Log: Add missing UA5 multiplicity analysis at 546 GeV and mark the correlation analysis VALIDATED. Added: trunk/data/anainfo/UA5_1987_S1640666.info trunk/data/plotinfo/UA5_1987_S1640666.plot trunk/data/refdata/UA5_1987_S1640666.aida trunk/src/Analyses/UA5_1987_S1640666.cc Modified: trunk/data/anainfo/Makefile.am trunk/data/anainfo/UA5_1988_S1867512.info trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/src/Analyses/Makefile.am Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Wed Jun 9 22:21:03 2010 (r2475) +++ trunk/data/anainfo/Makefile.am Thu Jun 10 08:40:21 2010 (r2476) @@ -92,5 +92,6 @@ UA1_1990_S2044935.info \ UA5_1982_S875503.info \ UA5_1986_S1583476.info \ + UA5_1987_S1640666.info \ UA5_1988_S1867512.info \ UA5_1989_S1926373.info Added: trunk/data/anainfo/UA5_1987_S1640666.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/UA5_1987_S1640666.info Thu Jun 10 08:40:21 2010 (r2476) @@ -0,0 +1,35 @@ +Name: UA5_1987_S1640666 +Year: 1987 +Summary: UA5 charged multiplicity measurements at 546 GeV +Experiment: UA5 +Collider: CERN SPS +SpiresID: 1640666 +Status: VALIDATED +Authors: + - Holger Schulz <holger.schulz at physik.hu-berlin.de> +References: + - Phys.Rept.154:247-383,1987 +RunInfo: + QCD and diffractive events at 546 GeV +NumEvents: 1000000 +Beams: [p-, p+] +Energies: [546] +PtCuts: +Description: + 'Charged particle multiplicity measurement.' +BibKey: Alner:1987wb +BibTeX: '@Article{Alner:1987wb, + author = "Alner, G. J. and others", + collaboration = "UA5", + title = "{UA5: A general study of proton-antiproton physics at + $\sqrt{s}$ = 546-GeV}", + journal = "Phys. Rept.", + volume = "154", + year = "1987", + pages = "247-383", + doi = "10.1016/0370-1573(87)90130-X", + SLACcitation = "%%CITATION = PRPLC,154,247;%%" +}' +ToDo: + - Include more distributions + Modified: trunk/data/anainfo/UA5_1988_S1867512.info ============================================================================== --- trunk/data/anainfo/UA5_1988_S1867512.info Wed Jun 9 22:21:03 2010 (r2475) +++ trunk/data/anainfo/UA5_1988_S1867512.info Thu Jun 10 08:40:21 2010 (r2476) @@ -4,7 +4,7 @@ Experiment: UA5 Collider: CERN SPS SpiresID: 1867512 -Status: UNVALIDATED +Status: VALIDATED Authors: - Holger Schulz <holger.schulz at physik.hu-berlin.de> References: @@ -12,12 +12,14 @@ RunInfo: ppbar events. Non-single diffractive events need to be switched on. The trigger implementation is the same as in UA5_1989_S1926373. + Important: Only the correlation strengths with symmetric eta bins + should be used for tuning. NumEvents: 100000 Beams: [p-, p+] Energies: [200, 546, 900] PtCuts: [0] Description: - Data on two-particle pseudorapidity and multiplicity correlations + 'Data on two-particle pseudorapidity and multiplicity correlations of charged particles for non single-diffractive $\bar{p}p$ collisions at c.m. energies of 200, 546 and 900 GeV. Pseudorapidity correlations interpreted in terms of a cluster model, which has been motivated by this and other @@ -29,7 +31,7 @@ depends strongly on position and size of the pseudorapidity gap separating the forward and backward interval. All our correlation studies can be understood in terms of a cluster model in which clusters contain on average about two - charged particles, i.e. are of similar magnitude to earlier estimates from the ISR. + charged particles, i.e. are of similar magnitude to earlier estimates from the ISR.' BibKey: Ansorge:1988fg BibTeX: '@Article{Ansorge:1988fg, author = "Ansorge, R. E. and others", @@ -43,3 +45,6 @@ doi = "10.1007/BF01579906", SLACcitation = "%%CITATION = ZEPYA,C37,191;%%" }' +ToDo: + - Check sanity of implementation of b with assymetric eta bins + - Get rid of vectors, use temp. histos instead Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Wed Jun 9 22:21:03 2010 (r2475) +++ trunk/data/plotinfo/Makefile.am Thu Jun 10 08:40:21 2010 (r2476) @@ -90,5 +90,6 @@ UA1_1990_S2044935.plot \ UA5_1982_S875503.plot \ UA5_1986_S1583476.plot \ + UA5_1987_S1640666.plot \ UA5_1988_S1867512.plot \ UA5_1989_S1926373.plot Added: trunk/data/plotinfo/UA5_1987_S1640666.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/UA5_1987_S1640666.plot Thu Jun 10 08:40:21 2010 (r2476) @@ -0,0 +1,11 @@ +# BEGIN PLOT /UA5_1987_S1640666/d01-x01-y01 +Title=Mean charged multiplicity at $\sqrt{s} = 546~\text{GeV}$, $|\eta| < 5.0$ +XLabel=$\text{GeV}$ +YLabel=$N_\text{ch}$ +# END PLOT + +# BEGIN PLOT /UA5_1987_S1640666/d03-x01-y01 +Title=Charged multiplicity at $\sqrt{s} = 546~\text{GeV}$, $|\eta| < 5.0$ +XLabel=$N_\text{ch}$ +YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{N_\text{ch}}$ +# END PLOT Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Wed Jun 9 22:21:03 2010 (r2475) +++ trunk/data/refdata/Makefile.am Thu Jun 10 08:40:21 2010 (r2476) @@ -35,6 +35,7 @@ UA5_1986_S1583476.aida \ UA5_1989_S1926373.aida \ UA5_1988_S1867512.aida \ + UA5_1987_S1640666.aida \ CDF_1988_S1865951.aida \ CDF_1990_S2089246.aida \ CDF_1991_S2313472.aida \ Added: trunk/data/refdata/UA5_1987_S1640666.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/UA5_1987_S1640666.aida Thu Jun 10 08:40:21 2010 (r2476) @@ -0,0 +1,207 @@ +<?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/UA5_1987_S1640666" title="MULT(P=3)" > + <dataPoint> + <measurement value="546.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="29.4" errorPlus="0.9486832980505138" errorMinus="0.9486832980505138"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/UA5_1987_S1640666" title="(1/N)*D(N)/DETARAP" > + <dataPoint> + <measurement value="546.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="3.05" errorPlus="0.03" errorMinus="0.03"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y01" dimension="2" path="/REF/UA5_1987_S1640666" title="PROB" > + <dataPoint> + <measurement value="2.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0027" errorPlus="8.0E-4" errorMinus="8.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="4.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0077" errorPlus="0.0013" errorMinus="0.0013"/> + </dataPoint> + <dataPoint> + <measurement value="6.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0122" errorPlus="0.0014" errorMinus="0.0014"/> + </dataPoint> + <dataPoint> + <measurement value="8.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0195" errorPlus="0.0017" errorMinus="0.0017"/> + </dataPoint> + <dataPoint> + <measurement value="10.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0324" errorPlus="0.0021" errorMinus="0.0021"/> + </dataPoint> + <dataPoint> + <measurement value="12.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0453" errorPlus="0.0025" errorMinus="0.0025"/> + </dataPoint> + <dataPoint> + <measurement value="14.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0531" errorPlus="0.0026" errorMinus="0.0026"/> + </dataPoint> + <dataPoint> + <measurement value="16.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0557" errorPlus="0.0027" errorMinus="0.0027"/> + </dataPoint> + <dataPoint> + <measurement value="18.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0566" errorPlus="0.0027" errorMinus="0.0027"/> + </dataPoint> + <dataPoint> + <measurement value="20.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0589" errorPlus="0.0027" errorMinus="0.0027"/> + </dataPoint> + <dataPoint> + <measurement value="22.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0608" errorPlus="0.0028" errorMinus="0.0028"/> + </dataPoint> + <dataPoint> + <measurement value="24.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0588" errorPlus="0.0027" errorMinus="0.0027"/> + </dataPoint> + <dataPoint> + <measurement value="26.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0551" errorPlus="0.0026" errorMinus="0.0026"/> + </dataPoint> + <dataPoint> + <measurement value="28.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0506" errorPlus="0.0025" errorMinus="0.0025"/> + </dataPoint> + <dataPoint> + <measurement value="30.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0451" errorPlus="0.0024" errorMinus="0.0024"/> + </dataPoint> + <dataPoint> + <measurement value="32.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0408" errorPlus="0.0023" errorMinus="0.0023"/> + </dataPoint> + <dataPoint> + <measurement value="34.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0373" errorPlus="0.0022" errorMinus="0.0022"/> + </dataPoint> + <dataPoint> + <measurement value="36.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0344" errorPlus="0.0021" errorMinus="0.0021"/> + </dataPoint> + <dataPoint> + <measurement value="38.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0334" errorPlus="0.0020" errorMinus="0.0020"/> + </dataPoint> + <dataPoint> + <measurement value="40.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0321" errorPlus="0.0020" errorMinus="0.0020"/> + </dataPoint> + <dataPoint> + <measurement value="42.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.029" errorPlus="0.0019" errorMinus="0.0019"/> + </dataPoint> + <dataPoint> + <measurement value="44.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0254" errorPlus="0.0018" errorMinus="0.0018"/> + </dataPoint> + <dataPoint> + <measurement value="46.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0208" errorPlus="0.0016" errorMinus="0.0016"/> + </dataPoint> + <dataPoint> + <measurement value="48.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0171" errorPlus="0.0015" errorMinus="0.0015"/> + </dataPoint> + <dataPoint> + <measurement value="50.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0156" errorPlus="0.0014" errorMinus="0.0014"/> + </dataPoint> + <dataPoint> + <measurement value="52.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.014" errorPlus="0.0013" errorMinus="0.0013"/> + </dataPoint> + <dataPoint> + <measurement value="54.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0124" errorPlus="0.0012" errorMinus="0.0012"/> + </dataPoint> + <dataPoint> + <measurement value="56.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0109" errorPlus="0.0012" errorMinus="0.0012"/> + </dataPoint> + <dataPoint> + <measurement value="58.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0095" errorPlus="0.0011" errorMinus="0.0011"/> + </dataPoint> + <dataPoint> + <measurement value="60.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0083" errorPlus="0.0010" errorMinus="0.0010"/> + </dataPoint> + <dataPoint> + <measurement value="62.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0072" errorPlus="9.0E-4" errorMinus="9.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="64.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0063" errorPlus="9.0E-4" errorMinus="9.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="66.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0055" errorPlus="8.0E-4" errorMinus="8.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="68.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0048" errorPlus="8.0E-4" errorMinus="8.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="70.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0040" errorPlus="7.0E-4" errorMinus="7.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="72.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0032" errorPlus="6.0E-4" errorMinus="6.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="74.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0024" errorPlus="6.0E-4" errorMinus="6.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="76.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0019" errorPlus="5.0E-4" errorMinus="5.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="78.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.0015" errorPlus="4.0E-4" errorMinus="4.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="80.0" errorPlus="2.0" errorMinus="1.0"/> + <measurement value="0.0013" errorPlus="4.0E-4" errorMinus="4.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="83.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="0.00111" errorPlus="2.6E-4" errorMinus="2.6E-4"/> + </dataPoint> + <dataPoint> + <measurement value="87.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="8.4E-4" errorPlus="2.3E-4" errorMinus="2.3E-4"/> + </dataPoint> + <dataPoint> + <measurement value="91.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="5.9E-4" errorPlus="1.9E-4" errorMinus="1.9E-4"/> + </dataPoint> + <dataPoint> + <measurement value="95.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="3.4E-4" errorPlus="1.5E-4" errorMinus="1.5E-4"/> + </dataPoint> + <dataPoint> + <measurement value="99.0" errorPlus="1.0" errorMinus="1.0"/> + <measurement value="1.5E-4" errorPlus="1.0E-4" errorMinus="1.0E-4"/> + </dataPoint> + <dataPoint> + <measurement value="106.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="5.0E-5" errorPlus="4.0E-5" errorMinus="4.0E-5"/> + </dataPoint> + <dataPoint> + <measurement value="116.0" errorPlus="4.0" errorMinus="4.0"/> + <measurement value="5.0E-5" errorPlus="4.0E-5" errorMinus="4.0E-5"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Wed Jun 9 22:21:03 2010 (r2475) +++ trunk/src/Analyses/Makefile.am Thu Jun 10 08:40:21 2010 (r2476) @@ -33,7 +33,6 @@ SFM_1984_S1178091.cc \ STAR_2008_S7993412.cc \ STAR_2009_UE_HELEN.cc \ - UA5_1988_S1867512.cc \ ZEUS_2001_S4815815.cc endif @@ -121,6 +120,8 @@ UA1_1990_S2044935.cc \ UA5_1982_S875503.cc \ UA5_1986_S1583476.cc \ + UA5_1987_S1640666.cc \ + UA5_1988_S1867512.cc \ UA5_1989_S1926373.cc Added: trunk/src/Analyses/UA5_1987_S1640666.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/UA5_1987_S1640666.cc Thu Jun 10 08:40:21 2010 (r2476) @@ -0,0 +1,76 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/TriggerUA5.hh" +#include "Rivet/Projections/ChargedFinalState.hh" + +namespace Rivet { + + + class UA5_1987_S1640666 : public Analysis { + public: + + /// Constructor + UA5_1987_S1640666() + : Analysis("UA5_1987_S1640666") + { + setBeams(PROTON, ANTIPROTON); + _sumWPassed = 0; + + } + + /// Book histograms and initialise projections before the run + void init() { + addProjection(TriggerUA5(), "Trigger"); + addProjection(ChargedFinalState(-5.0, 5.0), "CFS"); + + _hist_mean_nch = bookHistogram1D(1, 1, 1); + _hist_nch = bookHistogram1D(3, 1, 1); + + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + // Trigger + const TriggerUA5& trigger = applyProjection<TriggerUA5>(event, "Trigger"); + if (!trigger.nsdDecision()) vetoEvent; + + const double weight = event.weight(); + _sumWPassed += weight; + + // Count final state particles in several eta regions + const int Nch = applyProjection<ChargedFinalState>(event, "CFS").size(); + + // Fill histograms + _hist_nch->fill(Nch, weight); + _hist_mean_nch->fill(_hist_mean_nch->binMean(0), Nch); + + } + + + /// Normalise histograms etc., after the run + void finalize() { + + scale(_hist_nch, _sumWPassed); + scale(_hist_mean_nch, 1.0/_sumWPassed); + + } + + private: + + double _sumWPassed; + + AIDA::IHistogram1D *_hist_mean_nch; + AIDA::IHistogram1D *_hist_nch; + + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<UA5_1987_S1640666> plugin_UA5_1987_S1640666; + + +}
More information about the Rivet-svn mailing list |