|
|
[Rivet-svn] r2802 - in trunk: data/anainfo data/plotinfo data/refdata src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Thu Dec 2 15:07:31 GMT 2010
Author: holsch
Date: Thu Dec 2 15:07:30 2010
New Revision: 2802
Log:
Add LHCB analysis for differential cross-section measurement of prompt K0-short production, there are some issues with the normalisation, hence the status UNVALIDATED. Now only CMS is missing ;)
Added:
trunk/data/anainfo/LHCB_2010_S8758301.info
trunk/data/plotinfo/LHCB_2010_S8758301.plot
trunk/data/refdata/LHCB_2010_S8758301.aida
trunk/src/Analyses/LHCB_2010_S8758301.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/LHCB_2010_S8758301.info
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/data/anainfo/LHCB_2010_S8758301.info Thu Dec 2 15:07:30 2010 (r2802)
@@ -0,0 +1,42 @@
+Name: LHCB_2010_S8758301
+Year: 2010
+Summary: Differential cross section measurement of $K^0_S$ production in three
+ rapidity windows at $\sqrt{s}=\unit{0.9}{\TeV}$
+Experiment: LHCB
+Collider: LHCB
+SpiresID: 8758301
+Status: UNVALIDATED
+Authors:
+ - Holger Schulz holger.schulz at physik.hu-berlin.de
+References:
+ - Phys.Lett.B693:69-80,2010
+ - arXiv:1008.3105[hep-ex]
+RunInfo:
+ QCD events. See paper for MC discussion.
+NumEvents: 1000000
+Beams: [p+, p+]
+Energies: [900]
+PtCuts:
+Description:
+ 'Differential cross-section measurement of prompt $K^0_S$ production in pp collisions
+ at $\sqrt{s}=\unit{0.9}{\TeV}$ in the rapidity windows $2.5<y<3.0$, $3.0<y<3.5$ and $3.5<y<4.0$.
+ Unresolved issues at the moment are the normalisations. Not clear to me how to do this from the
+ paper so I normalise to reference data area.'
+BibKey: Aaij:2010nx
+BibTeX: '@Article{Aaij:2010nx,
+ author = "Aaij, R and others",
+ collaboration = "LHCb",
+ title = "{Prompt K_short production in pp collisions at sqrt(s)=0.9
+ TeV}",
+ journal = "Phys. Lett.",
+ volume = "B693",
+ year = "2010",
+ pages = "69-80",
+ eprint = "1008.3105",
+ archivePrefix = "arXiv",
+ primaryClass = "Unknown",
+ doi = "10.1016/j.physletb.2010.08.055",
+ SLACcitation = "%%CITATION = 1008.3105;%%"
+}'
+ToDo
+ - Do not normalise to data!
Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am Thu Dec 2 14:46:40 2010 (r2801)
+++ trunk/data/anainfo/Makefile.am Thu Dec 2 15:07:30 2010 (r2802)
@@ -78,6 +78,7 @@
H1_1995_S3167097.info \
H1_2000_S4129130.info \
JADE_OPAL_2000_S4300807.info \
+ LHCB_2010_S8758301.info \
MC_GENERIC.info \
MC_LEADINGJETS.info \
MC_DIJET.info \
Added: trunk/data/plotinfo/LHCB_2010_S8758301.plot
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/data/plotinfo/LHCB_2010_S8758301.plot Thu Dec 2 15:07:30 2010 (r2802)
@@ -0,0 +1,17 @@
+# BEGIN PLOT /LHCB_2010_S8758301/d01-x01-y01
+Title=Diff. cross-section for $K_S^0$ in production ($2.5<y<3.0$)
+XLabel=$p_\perp$ [GeV]
+YLabel=$\text{d}^2\sigma(pp\rightarrow\text{K}_S^0 X) / (\text{d}y\text{d}p_\perp)$
+# END PLOT
+
+# BEGIN PLOT /LHCB_2010_S8758301/d01-x01-y02
+Title=Diff. cross-section for $K_S^0$ in production ($3.0<y<3.5$)
+XLabel=$p_\perp$ [GeV]
+YLabel=$\text{d}^2\sigma(pp\rightarrow\text{K}_S^0 X) / (\text{d}y\text{d}p_\perp)$
+# END PLOT
+
+# BEGIN PLOT /LHCB_2010_S8758301/d01-x01-y03
+Title=Diff. cross-section for $K_S^0$ in production ($3.5<y<4.0$)
+XLabel=$p_\perp$ [GeV]
+YLabel=$\text{d}^2\sigma(pp\rightarrow\text{K}_S^0 X) / (\text{d}y\text{d}p_\perp)$
+# END PLOT
Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am Thu Dec 2 14:46:40 2010 (r2801)
+++ trunk/data/plotinfo/Makefile.am Thu Dec 2 15:07:30 2010 (r2802)
@@ -76,6 +76,7 @@
H1_1995_S3167097.plot \
H1_2000_S4129130.plot \
JADE_OPAL_2000_S4300807.plot \
+ LHCB_2010_S8758301.plot \
MC_GENERIC.plot \
MC_LEADINGJETS.plot \
MC_DIJET.plot \
Added: trunk/data/refdata/LHCB_2010_S8758301.aida
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/data/refdata/LHCB_2010_S8758301.aida Thu Dec 2 15:07:30 2010 (r2802)
@@ -0,0 +1,107 @@
+<?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/LHCB_2010_S8758301" title="SIG IN MUB" >
+ <dataPoint>
+ <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/>
+ <measurement value="294.0" errorPlus="126.26955294131677" errorMinus="126.26955294131677"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.30000000000000004" errorPlus="0.09999999999999998" errorMinus="0.10000000000000003"/>
+ <measurement value="649.0" errorPlus="263.9583300447251" errorMinus="263.9583300447251"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.5" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/>
+ <measurement value="618.0" errorPlus="133.16906547693424" errorMinus="133.16906547693424"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.7" errorPlus="0.10000000000000009" errorMinus="0.09999999999999998"/>
+ <measurement value="401.0" errorPlus="73.78346698278686" errorMinus="73.78346698278686"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.9" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/>
+ <measurement value="232.0" errorPlus="42.24926034855522" errorMinus="42.24926034855522"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="1.1" errorPlus="0.09999999999999987" errorMinus="0.10000000000000009"/>
+ <measurement value="115.0" errorPlus="21.470910553583888" errorMinus="21.470910553583888"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="1.2999999999999998" errorPlus="0.10000000000000009" errorMinus="0.09999999999999987"/>
+ <measurement value="85.0" errorPlus="17.46424919657298" errorMinus="17.46424919657298"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="1.5" errorPlus="0.10000000000000009" errorMinus="0.10000000000000009"/>
+ <measurement value="46.0" errorPlus="10.099504938362077" errorMinus="10.099504938362077"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet name="d01-x01-y02" dimension="2" path="/REF/LHCB_2010_S8758301" title="SIG IN MUB" >
+ <dataPoint>
+ <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/>
+ <measurement value="316.0" errorPlus="94.70480452437458" errorMinus="94.70480452437458"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.30000000000000004" errorPlus="0.09999999999999998" errorMinus="0.10000000000000003"/>
+ <measurement value="562.0" errorPlus="111.57508682497182" errorMinus="111.57508682497182"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.5" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/>
+ <measurement value="534.0" errorPlus="91.86947262284681" errorMinus="91.86947262284681"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.7" errorPlus="0.10000000000000009" errorMinus="0.09999999999999998"/>
+ <measurement value="371.0" errorPlus="63.26926584053272" errorMinus="63.26926584053272"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.9" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/>
+ <measurement value="183.0" errorPlus="32.7566787083184" errorMinus="32.7566787083184"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="1.1" errorPlus="0.09999999999999987" errorMinus="0.10000000000000009"/>
+ <measurement value="136.0" errorPlus="24.779023386727733" errorMinus="24.779023386727733"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="1.2999999999999998" errorPlus="0.10000000000000009" errorMinus="0.09999999999999987"/>
+ <measurement value="70.0" errorPlus="13.74772708486752" errorMinus="13.74772708486752"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="1.5" errorPlus="0.10000000000000009" errorMinus="0.10000000000000009"/>
+ <measurement value="49.0" errorPlus="10.44030650891055" errorMinus="10.44030650891055"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet name="d01-x01-y03" dimension="2" path="/REF/LHCB_2010_S8758301" title="SIG IN MUB" >
+ <dataPoint>
+ <measurement value="0.1" errorPlus="0.1" errorMinus="0.1"/>
+ <measurement value="196.0" errorPlus="66.9776081985614" errorMinus="66.9776081985614"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.30000000000000004" errorPlus="0.09999999999999998" errorMinus="0.10000000000000003"/>
+ <measurement value="571.0" errorPlus="124.03628501370073" errorMinus="124.03628501370073"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.5" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/>
+ <measurement value="477.0" errorPlus="82.46817568007674" errorMinus="82.46817568007674"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.7" errorPlus="0.10000000000000009" errorMinus="0.09999999999999998"/>
+ <measurement value="323.0" errorPlus="55.51576352712804" errorMinus="55.51576352712804"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="0.9" errorPlus="0.09999999999999998" errorMinus="0.09999999999999998"/>
+ <measurement value="201.0" errorPlus="37.16180835212409" errorMinus="37.16180835212409"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="1.1" errorPlus="0.09999999999999987" errorMinus="0.10000000000000009"/>
+ <measurement value="108.0" errorPlus="21.400934559032695" errorMinus="21.400934559032695"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="1.2999999999999998" errorPlus="0.10000000000000009" errorMinus="0.09999999999999987"/>
+ <measurement value="35.0" errorPlus="11.224972160321824" errorMinus="11.224972160321824"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement value="1.5" errorPlus="0.10000000000000009" errorMinus="0.10000000000000009"/>
+ <measurement value="23.0" errorPlus="7.483314773547883" errorMinus="7.483314773547883"/>
+ </dataPoint>
+ </dataPointSet>
+</aida>
Modified: trunk/data/refdata/Makefile.am
==============================================================================
--- trunk/data/refdata/Makefile.am Thu Dec 2 14:46:40 2010 (r2801)
+++ trunk/data/refdata/Makefile.am Thu Dec 2 15:07:30 2010 (r2802)
@@ -1,4 +1,5 @@
dist_pkgdata_DATA = \
+ LHCB_2010_S8758301.aida \
JADE_OPAL_2000_S4300807.aida \
ALEPH_1991_S2435284.aida \
ALEPH_1996_S3486095.aida \
Added: trunk/src/Analyses/LHCB_2010_S8758301.cc
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/src/Analyses/LHCB_2010_S8758301.cc Thu Dec 2 15:07:30 2010 (r2802)
@@ -0,0 +1,128 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/UnstableFinalState.hh"
+#include "Rivet/Tools/ParticleIdUtils.hh"
+/// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder...
+
+namespace Rivet {
+
+
+ class LHCB_2010_S8758301 : public Analysis {
+ public:
+
+ /// @name Constructors etc.
+ //@{
+
+ /// Constructor
+ LHCB_2010_S8758301()
+ : Analysis("LHCB_2010_S8758301"),
+ sumKs0_30(0.0), sumKs0_35(0.0), sumKs0_40(0.0)
+
+ {
+ /// @todo Set whether your finalize method needs the generator cross section
+ setNeedsCrossSection(false);
+ setBeams(PROTON, PROTON);
+ }
+
+ //@}
+
+
+ public:
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+
+ addProjection(UnstableFinalState(), "UFS");
+
+ _h_K0s_pt_y_30 = bookHistogram1D(1,1,1);
+ _h_K0s_pt_y_35 = bookHistogram1D(1,1,2);
+ _h_K0s_pt_y_40 = bookHistogram1D(1,1,3);
+
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const double weight = event.weight();
+ const UnstableFinalState& ufs = applyProjection<UnstableFinalState>(event, "UFS");
+
+ foreach (const Particle& p, ufs.particles()) {
+ const PdgId id = abs(p.pdgId());
+
+ if (id == 310) {
+ double y = p.momentum().rapidity();
+ double pT = p.momentum().perp();
+ if (y > 2.5 && y < 3.0) {
+ //_h_K0s_pt_y_30->fill(pT, weight/(0.5*0.2));
+ _h_K0s_pt_y_30->fill(pT, weight);
+ sumKs0_30 += weight;
+ }
+ else if (y > 3.0 && y < 3.5) {
+ //_h_K0s_pt_y_35->fill(pT, weight/(0.5*0.2));
+ _h_K0s_pt_y_35->fill(pT, weight);
+ sumKs0_35 += weight;
+ }
+ else if (y > 3.5 && y < 4.0) {
+ //_h_K0s_pt_y_40->fill(pT, weight/(0.5*0.2));
+ _h_K0s_pt_y_40->fill(pT, weight);
+ sumKs0_40 += weight;
+ }
+ }
+ }
+
+
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+
+
+ // here we normalise to reference data...
+ scale(_h_K0s_pt_y_30, 4.880000e+02/sumKs0_30);
+ scale(_h_K0s_pt_y_35, 4.442000e+02/sumKs0_35);
+ scale(_h_K0s_pt_y_40, 3.868000e+02/sumKs0_40);
+
+
+
+ }
+
+ //@}
+
+
+ private:
+
+ // Data members like post-cuts event weight counters go here
+
+
+ private:
+
+ /// @name Histograms
+ //@{
+
+
+ AIDA::IHistogram1D *_h_K0s_pt_y_30;
+ AIDA::IHistogram1D *_h_K0s_pt_y_35;
+ AIDA::IHistogram1D *_h_K0s_pt_y_40;
+
+ double sumKs0_30;
+ double sumKs0_35;
+ double sumKs0_40;
+
+ //@}
+
+
+ };
+
+
+
+ // This global object acts as a hook for the plugin system
+ AnalysisBuilder<LHCB_2010_S8758301> plugin_LHCB_2010_S8758301;
+
+
+}
Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am Thu Dec 2 14:46:40 2010 (r2801)
+++ trunk/src/Analyses/Makefile.am Thu Dec 2 15:07:30 2010 (r2802)
@@ -37,6 +37,7 @@
STAR_2008_S7869363.cc \
STAR_2008_S7993412.cc \
STAR_2009_UE_HELEN.cc \
+ LHCB_2010_S8758301.cc \
ZEUS_2001_S4815815.cc
endif
More information about the Rivet-svn
mailing list
|