[Rivet-svn] r2933 - in trunk: data/anainfo data/plotinfo data/refdata src/Analyses

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Tue Feb 15 17:51:37 GMT 2011


Author: hoeth
Date: Tue Feb 15 17:51:37 2011
New Revision: 2933

Log:
Add CMS event shape analysis CMS_2011_S8957746.
Data is read off the plots, you know the deal ...

Added:
   trunk/data/anainfo/CMS_2011_S8957746.info
   trunk/data/plotinfo/CMS_2011_S8957746.plot
   trunk/data/refdata/CMS_2011_S8957746.aida
   trunk/src/Analyses/CMS_2011_S8957746.cc
Modified:
   trunk/data/plotinfo/Makefile.am
   trunk/data/refdata/Makefile.am
   trunk/src/Analyses/Makefile.am

Added: trunk/data/anainfo/CMS_2011_S8957746.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/anainfo/CMS_2011_S8957746.info	Tue Feb 15 17:51:37 2011	(r2933)
@@ -0,0 +1,24 @@
+Name: CMS_2011_S8957746
+Year: 2011
+Summary: Event shapes
+Experiment: CMS
+Collider: LHC
+SpiresID: 8957746
+Status: VALIDATED
+Authors:
+ - Hendrik Hoeth <hendrik.hoeth at cern.ch>
+References:
+ - arXiv:1102.0068
+RunInfo:
+  pp QCD interactions at 7000 GeV. Particles with c*tau>10mm are stable.
+NumEvents: 1000000000
+Beams: [p+, p+]
+Energies: [7000]
+PtCuts: [30.0]
+Description:
+  Central transverse Thrust and Minor measured by ATLAS at 7 TeV. Data is
+  read off the plots.
+BibKey: TODO
+BibTeX: 'TODO'
+ToDo:
+  Add BibKey and BibTeX as soon as it is fixed on Spires.

Added: trunk/data/plotinfo/CMS_2011_S8957746.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/plotinfo/CMS_2011_S8957746.plot	Tue Feb 15 17:51:37 2011	(r2933)
@@ -0,0 +1,38 @@
+# BEGIN PLOT /CMS_2011_S8957746/d
+LegendYPos=0.4
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /CMS_2011_S8957746/d01
+XLabel=$\ln(1-T_\mathcal{C})$
+YLabel=$1/N \,\, \text{d}N/\text{d}\,\ln(1-T_\mathcal{C})$
+# END PLOT
+
+# BEGIN PLOT /CMS_2011_S8957746/d02
+XLabel=$\ln(T_{m,\mathcal{C}})$
+YLabel=$1/N \,\, \text{d}N/\text{d}\,\ln(T_{m,\mathcal{C}})$
+# END PLOT
+
+# BEGIN PLOT /CMS_2011_S8957746/d01-x01-y01
+Title=Central Transv. Thrust, $90\,\text{GeV} < p_\perp^\text{jet 1} < 125\,\text{GeV}$, $\sqrt{s}=7\,\text{TeV}$
+# END PLOT
+
+# BEGIN PLOT /CMS_2011_S8957746/d01-x02-y01
+Title=Central Transv. Thrust, $125\,\text{GeV} < p_\perp^\text{jet 1} < 200\,\text{GeV}$, $\sqrt{s}=7\,\text{TeV}$
+# END PLOT
+
+# BEGIN PLOT /CMS_2011_S8957746/d01-x03-y01
+Title=Central Transv. Thrust, $p_\perp^\text{jet 1} > 200\,\text{GeV}$, $\sqrt{s}=7\,\text{TeV}$
+# END PLOT
+
+# BEGIN PLOT /CMS_2011_S8957746/d02-x01-y01
+Title=Central Transv. Minor, $90\,\text{GeV} < p_\perp^\text{jet 1} < 125\,\text{GeV}$, $\sqrt{s}=7\,\text{TeV}$
+# END PLOT
+
+# BEGIN PLOT /CMS_2011_S8957746/d02-x02-y01
+Title=Central Transv. Minor, $125\,\text{GeV} < p_\perp^\text{jet 1} < 200\,\text{GeV}$, $\sqrt{s}=7\,\text{TeV}$
+# END PLOT
+
+# BEGIN PLOT /CMS_2011_S8957746/d02-x03-y01
+Title=Central Transv. Minor, $p_\perp^\text{jet 1} > 200\,\text{GeV}$, $\sqrt{s}=7\,\text{TeV}$
+# END PLOT

Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am	Tue Feb 15 17:44:00 2011	(r2932)
+++ trunk/data/plotinfo/Makefile.am	Tue Feb 15 17:51:37 2011	(r2933)
@@ -48,6 +48,7 @@
   CDF_2009_NOTE_9936.plot \
   CDF_2009_S8383952.plot \
   CDF_2009_S8436959.plot \
+  CMS_2011_S8957746.plot \
   D0_1996_S3214044.plot \
   D0_1996_S3324664.plot \
   D0_1998_S3711838.plot \

Added: trunk/data/refdata/CMS_2011_S8957746.aida
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/refdata/CMS_2011_S8957746.aida	Tue Feb 15 17:51:37 2011	(r2933)
@@ -0,0 +1,365 @@
+<?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.1" package="FreeHEP"/>
+  <dataPointSet name="d01-x01-y01" dimension="2"
+    path="/REF/CMS_2011_S8957746" title="transverse Thrust, $90\,\text{GeV} < p_\perp(leading jet) \le 125\,\text{GeV}$">
+    <annotation>
+      <item key="Title" value="transverse Thrust, $90\,\text{GeV} < p_\perp(leading jet) \le 125\,\text{GeV}$" sticky="true"/>
+      <item key="AidaPath" value="/REF/CMS_2011_S8957746/d01-x01-y01" sticky="true"/>
+      <item key="FullPath" value="/CMS_2011_S8957746.aida/REF/CMS_2011_S8957746/d01-x01-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="1.300000e+00" value="-1.245000e+01" errorMinus="1.300000e+00"/>
+      <measurement errorPlus="2.360000e-03" value="1.277000e-02" errorMinus="2.360000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="9.750000e-01" value="-1.017500e+01" errorMinus="9.750000e-01"/>
+      <measurement errorPlus="5.140000e-03" value="2.828000e-02" errorMinus="5.140000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="6.500000e-01" value="-8.550000e+00" errorMinus="6.500000e-01"/>
+      <measurement errorPlus="5.140000e-03" value="5.817000e-02" errorMinus="5.140000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="6.500000e-01" value="-7.250000e+00" errorMinus="6.500000e-01"/>
+      <measurement errorPlus="4.120000e-03" value="9.763000e-02" errorMinus="4.120000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-6.275000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="5.580000e-03" value="1.361700e-01" errorMinus="5.580000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-5.625000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="4.980000e-03" value="1.573400e-01" errorMinus="4.980000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-4.975000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="4.450000e-03" value="1.588000e-01" errorMinus="4.450000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-4.325000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="4.360000e-03" value="1.613100e-01" errorMinus="4.360000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-3.675000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="4.430000e-03" value="1.492900e-01" errorMinus="4.430000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-3.025000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="4.280000e-03" value="1.302900e-01" errorMinus="4.280000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-2.375000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="4.120000e-03" value="1.125300e-01" errorMinus="4.120000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-1.725000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="3.350000e-03" value="7.636000e-02" errorMinus="3.350000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-1.075000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="7.700000e-04" value="1.000000e-02" errorMinus="7.700000e-04"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d02-x01-y01" dimension="2"
+    path="/REF/CMS_2011_S8957746" title="transverse minor, $90\,\text{GeV} < p_\perp(leading jet) \le 125\,\text{GeV}$">
+    <annotation>
+      <item key="Title" value="transverse minor, $90\,\text{GeV} < p_\perp(leading jet) \le 125\,\text{GeV}$" sticky="true"/>
+      <item key="AidaPath" value="/REF/CMS_2011_S8957746/d02-x01-y01" sticky="true"/>
+      <item key="FullPath" value="/CMS_2011_S8957746.aida/REF/CMS_2011_S8957746/d02-x01-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="6.000000e-01" value="-5.650000e+00" errorMinus="6.000000e-01"/>
+      <measurement errorPlus="4.570000e-03" value="3.390000e-02" errorMinus="4.570000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.500000e-01" value="-4.600000e+00" errorMinus="4.500000e-01"/>
+      <measurement errorPlus="9.330000e-03" value="6.530000e-02" errorMinus="9.330000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.000000e-01" value="-3.850000e+00" errorMinus="3.000000e-01"/>
+      <measurement errorPlus="1.011000e-02" value="1.271000e-01" errorMinus="1.011000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.000000e-01" value="-3.250000e+00" errorMinus="3.000000e-01"/>
+      <measurement errorPlus="8.460000e-03" value="2.021000e-01" errorMinus="8.460000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-2.800000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="1.215000e-02" value="2.818000e-01" errorMinus="1.215000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-2.500000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="1.108000e-02" value="3.304000e-01" errorMinus="1.108000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-2.200000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="1.050000e-02" value="3.440000e-01" errorMinus="1.050000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-1.900000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="1.001000e-02" value="3.546000e-01" errorMinus="1.001000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-1.600000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="1.079000e-02" value="3.542000e-01" errorMinus="1.079000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-1.300000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="1.274000e-02" value="3.307000e-01" errorMinus="1.274000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-1.000000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="1.215000e-02" value="2.287000e-01" errorMinus="1.215000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-7.000000e-01" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="7.190000e-03" value="1.034000e-01" errorMinus="7.190000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-4.000000e-01" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="1.550000e-03" value="1.530000e-02" errorMinus="1.550000e-03"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d01-x02-y01" dimension="2"
+    path="/REF/CMS_2011_S8957746" title="transverse Thrust, $125\,\text{GeV} < p_\perp(leading jet) \le 200\,\text{GeV}$">
+    <annotation>
+      <item key="Title" value="transverse Thrust, $125\,\text{GeV} < p_\perp(leading jet) \le 200\,\text{GeV}$" sticky="true"/>
+      <item key="AidaPath" value="/REF/CMS_2011_S8957746/d01-x02-y01" sticky="true"/>
+      <item key="FullPath" value="/CMS_2011_S8957746.aida/REF/CMS_2011_S8957746/d01-x02-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="1.300000e+00" value="-1.245000e+01" errorMinus="1.300000e+00"/>
+      <measurement errorPlus="9.900000e-04" value="1.436000e-02" errorMinus="9.900000e-04"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="9.750000e-01" value="-1.017500e+01" errorMinus="9.750000e-01"/>
+      <measurement errorPlus="2.060000e-03" value="3.036000e-02" errorMinus="2.060000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="6.500000e-01" value="-8.550000e+00" errorMinus="6.500000e-01"/>
+      <measurement errorPlus="2.700000e-03" value="6.200000e-02" errorMinus="2.700000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-7.575000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="4.190000e-03" value="9.065000e-02" errorMinus="4.190000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-6.925000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="3.830000e-03" value="1.080700e-01" errorMinus="3.830000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-6.275000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="3.410000e-03" value="1.240000e-01" errorMinus="3.410000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-5.625000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="3.190000e-03" value="1.337400e-01" errorMinus="3.190000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-4.975000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="3.120000e-03" value="1.358100e-01" errorMinus="3.120000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-4.325000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="3.480000e-03" value="1.464000e-01" errorMinus="3.480000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-3.675000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="4.050000e-03" value="1.560700e-01" errorMinus="4.050000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-3.025000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="4.190000e-03" value="1.474000e-01" errorMinus="4.190000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-2.375000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="4.400000e-03" value="1.375100e-01" errorMinus="4.400000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-1.725000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="2.630000e-03" value="7.729000e-02" errorMinus="2.630000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.250000e-01" value="-1.075000e+00" errorMinus="3.250000e-01"/>
+      <measurement errorPlus="4.200000e-04" value="8.240000e-03" errorMinus="4.200000e-04"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d02-x02-y01" dimension="2"
+    path="/REF/CMS_2011_S8957746" title="transverse minor, $125\,\text{GeV} < p_\perp(leading jet) \le 200\,\text{GeV}$">
+    <annotation>
+      <item key="Title" value="transverse minor, $125\,\text{GeV} < p_\perp(leading jet) \le 200\,\text{GeV}$" sticky="true"/>
+      <item key="AidaPath" value="/REF/CMS_2011_S8957746/d02-x02-y01" sticky="true"/>
+      <item key="FullPath" value="/CMS_2011_S8957746.aida/REF/CMS_2011_S8957746/d02-x02-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="6.000000e-01" value="-5.650000e+00" errorMinus="6.000000e-01"/>
+      <measurement errorPlus="1.910000e-03" value="3.710300e-02" errorMinus="1.910000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.500000e-01" value="-4.600000e+00" errorMinus="4.500000e-01"/>
+      <measurement errorPlus="3.830000e-03" value="7.037400e-02" errorMinus="3.830000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="3.000000e-01" value="-3.850000e+00" errorMinus="3.000000e-01"/>
+      <measurement errorPlus="5.220000e-03" value="1.358700e-01" errorMinus="5.220000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-3.400000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="8.180000e-03" value="1.924830e-01" errorMinus="8.180000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-3.100000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="7.830000e-03" value="2.269740e-01" errorMinus="7.830000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-2.800000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="7.490000e-03" value="2.656450e-01" errorMinus="7.490000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-2.500000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="6.960000e-03" value="2.940380e-01" errorMinus="6.960000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-2.200000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="7.830000e-03" value="3.280060e-01" errorMinus="7.830000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-1.900000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="8.530000e-03" value="3.729480e-01" errorMinus="8.530000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-1.600000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="9.920000e-03" value="3.823540e-01" errorMinus="9.920000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-1.300000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="1.167000e-02" value="3.283540e-01" errorMinus="1.167000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-1.000000e+00" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="1.045000e-02" value="2.107740e-01" errorMinus="1.045000e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-7.000000e-01" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="5.390000e-03" value="8.831600e-02" errorMinus="5.390000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.500000e-01" value="-4.000000e-01" errorMinus="1.500000e-01"/>
+      <measurement errorPlus="8.700000e-04" value="1.132200e-02" errorMinus="8.700000e-04"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d01-x03-y01" dimension="2"
+    path="/REF/CMS_2011_S8957746" title="transverse Thrust, $p_\perp(leading jet) > 200\,\text{GeV}$">
+    <annotation>
+      <item key="Title" value="transverse Thrust, $p_\perp(leading jet) > 200\,\text{GeV}$" sticky="true"/>
+      <item key="AidaPath" value="/REF/CMS_2011_S8957746/d01-x03-y01" sticky="true"/>
+      <item key="FullPath" value="/CMS_2011_S8957746.aida/REF/CMS_2011_S8957746/d01-x03-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="1.225000e+00" value="-1.252500e+01" errorMinus="1.225000e+00"/>
+      <measurement errorPlus="1.362000e-03" value="1.736800e-02" errorMinus="1.362000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.220000e+00" value="-1.008000e+01" errorMinus="1.220000e+00"/>
+      <measurement errorPlus="2.081000e-03" value="3.873300e-02" errorMinus="2.081000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.050000e-01" value="-8.455000e+00" errorMinus="4.050000e-01"/>
+      <measurement errorPlus="4.321000e-03" value="6.946200e-02" errorMinus="4.321000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.050000e-01" value="-7.645000e+00" errorMinus="4.050000e-01"/>
+      <measurement errorPlus="4.244000e-03" value="8.635000e-02" errorMinus="4.244000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.050000e-01" value="-6.835000e+00" errorMinus="4.050000e-01"/>
+      <measurement errorPlus="3.681000e-03" value="9.747100e-02" errorMinus="3.681000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.050000e-01" value="-6.025000e+00" errorMinus="4.050000e-01"/>
+      <measurement errorPlus="3.441000e-03" value="1.020330e-01" errorMinus="3.441000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.050000e-01" value="-5.215000e+00" errorMinus="4.050000e-01"/>
+      <measurement errorPlus="3.443000e-03" value="1.092360e-01" errorMinus="3.443000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.050000e-01" value="-4.405000e+00" errorMinus="4.050000e-01"/>
+      <measurement errorPlus="4.083000e-03" value="1.399660e-01" errorMinus="4.083000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.050000e-01" value="-3.595000e+00" errorMinus="4.050000e-01"/>
+      <measurement errorPlus="4.883000e-03" value="1.726160e-01" errorMinus="4.883000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.050000e-01" value="-2.785000e+00" errorMinus="4.050000e-01"/>
+      <measurement errorPlus="5.043000e-03" value="1.751780e-01" errorMinus="5.043000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.100000e-01" value="-1.970000e+00" errorMinus="4.100000e-01"/>
+      <measurement errorPlus="3.361000e-03" value="9.843100e-02" errorMinus="3.361000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="4.050000e-01" value="-1.155000e+00" errorMinus="4.050000e-01"/>
+      <measurement errorPlus="9.620000e-04" value="1.200500e-02" errorMinus="9.620000e-04"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d02-x03-y01" dimension="2"
+    path="/REF/CMS_2011_S8957746" title="transverse minor, $p_\perp(leading jet) > 200\,\text{GeV}$">
+    <annotation>
+      <item key="Title" value="transverse minor, $p_\perp(leading jet) > 200\,\text{GeV}$" sticky="true"/>
+      <item key="AidaPath" value="/REF/CMS_2011_S8957746/d02-x03-y01" sticky="true"/>
+      <item key="FullPath" value="/CMS_2011_S8957746.aida/REF/CMS_2011_S8957746/d02-x03-y01" sticky="true"/>
+    </annotation>
+    <dataPoint>
+      <measurement errorPlus="5.625000e-01" value="-5.687500e+00" errorMinus="5.625000e-01"/>
+      <measurement errorPlus="3.107000e-03" value="4.624000e-02" errorMinus="3.107000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="5.625000e-01" value="-4.562500e+00" errorMinus="5.625000e-01"/>
+      <measurement errorPlus="4.569000e-03" value="8.754000e-02" errorMinus="4.569000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.875000e-01" value="-3.812500e+00" errorMinus="1.875000e-01"/>
+      <measurement errorPlus="9.138000e-03" value="1.506000e-01" errorMinus="9.138000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.875000e-01" value="-3.437500e+00" errorMinus="1.875000e-01"/>
+      <measurement errorPlus="8.772000e-03" value="1.875100e-01" errorMinus="8.772000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.875000e-01" value="-3.062500e+00" errorMinus="1.875000e-01"/>
+      <measurement errorPlus="8.224000e-03" value="2.224200e-01" errorMinus="8.224000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.875000e-01" value="-2.687500e+00" errorMinus="1.875000e-01"/>
+      <measurement errorPlus="8.407000e-03" value="2.551400e-01" errorMinus="8.407000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.875000e-01" value="-2.312500e+00" errorMinus="1.875000e-01"/>
+      <measurement errorPlus="9.321000e-03" value="3.211200e-01" errorMinus="9.321000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.875000e-01" value="-1.937500e+00" errorMinus="1.875000e-01"/>
+      <measurement errorPlus="1.114800e-02" value="4.004400e-01" errorMinus="1.114800e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.875000e-01" value="-1.562500e+00" errorMinus="1.875000e-01"/>
+      <measurement errorPlus="1.151400e-02" value="3.649800e-01" errorMinus="1.151400e-02"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.875000e-01" value="-1.187500e+00" errorMinus="1.875000e-01"/>
+      <measurement errorPlus="9.321000e-03" value="2.487400e-01" errorMinus="9.321000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.875000e-01" value="-8.125000e-01" errorMinus="1.875000e-01"/>
+      <measurement errorPlus="4.934000e-03" value="1.025300e-01" errorMinus="4.934000e-03"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement errorPlus="1.875000e-01" value="-4.375000e-01" errorMinus="1.875000e-01"/>
+      <measurement errorPlus="1.096000e-03" value="1.133000e-02" errorMinus="1.096000e-03"/>
+    </dataPoint>
+  </dataPointSet>
+</aida>

Modified: trunk/data/refdata/Makefile.am
==============================================================================
--- trunk/data/refdata/Makefile.am	Tue Feb 15 17:44:00 2011	(r2932)
+++ trunk/data/refdata/Makefile.am	Tue Feb 15 17:51:37 2011	(r2933)
@@ -16,6 +16,7 @@
   ATLAS_2011_S8971293.aida       \
   ATLAS_2010_CONF_2010_049.aida  \
   ATLAS_2010_CONF_2010_083.aida  \
+  CMS_2011_S8957746.aida         \
   DELPHI_1994_S3021912.aida      \
   DELPHI_1995_S3137023.aida      \
   DELPHI_1996_S3430090.aida      \

Added: trunk/src/Analyses/CMS_2011_S8957746.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/src/Analyses/CMS_2011_S8957746.cc	Tue Feb 15 17:51:37 2011	(r2933)
@@ -0,0 +1,96 @@
+// -*- C++ -*-
+
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/FastJets.hh"
+#include "Rivet/Projections/Thrust.hh"
+#include "Rivet/Tools/Logging.hh"
+
+namespace Rivet {
+
+  /// Rivet analysis class for CMS_2011_S8957746 dataset
+  class CMS_2011_S8957746 : public Analysis {
+  public:
+
+    /// Default constructor
+    CMS_2011_S8957746() : Analysis("CMS_2011_S8957746") {
+      setBeams(PROTON, PROTON);
+    }
+
+    /// Initialization, called once before running
+    void init() {
+      // Projections
+      const FastJets jets(FinalState(-5.0, 5.0, 0.0*GeV), FastJets::ANTIKT, 0.5);
+      addProjection(jets, "Jets");
+
+      // Book histograms
+      _hist_T_90  = bookHistogram1D(1, 1, 1);
+      _hist_m_90  = bookHistogram1D(2, 1, 1);
+      _hist_T_125 = bookHistogram1D(1, 2, 1);
+      _hist_m_125 = bookHistogram1D(2, 2, 1);
+      _hist_T_200 = bookHistogram1D(1, 3, 1);
+      _hist_m_200 = bookHistogram1D(2, 3, 1);
+    }
+
+    void analyze(const Event& event) {
+      const double weight = event.weight();
+      const Jets& jets = applyProjection<FastJets>(event, "Jets").jetsByPt(30.0*GeV);
+      if (jets.size() < 2 ||
+          fabs(jets[0].momentum().eta())>=1.3 ||
+          fabs(jets[1].momentum().eta())>=1.3 ||
+          jets[0].momentum().pT()<90) {
+        vetoEvent;
+      }
+      std::vector<Vector3> momenta;
+      foreach (Jet j, jets) {
+        if (fabs(j.momentum().eta())<1.3) {
+          Vector3 mom = j.momentum().vector3();
+          mom.setZ(0.0);
+          momenta.push_back(mom);
+        }
+      }
+      if (momenta.size()==2) {
+        // We need to use a ghost so that Thrust.calc() doesn't return 1.
+        momenta.push_back(Vector3(1e-10*MeV, 0., 0.));
+      }
+      Thrust thrust;
+      thrust.calc(momenta);
+
+      // The lowest bin also includes the underflow:
+      const double T=max(log(1-thrust.thrust()), -12.0);
+      const double M=max(log(thrust.thrustMajor()), -6.0);
+      if (jets[0].momentum().pT()>200) {
+        _hist_T_200->fill(T, weight);
+        _hist_m_200->fill(M, weight);
+      } else if (jets[0].momentum().pT()>125) {
+        _hist_T_125->fill(T, weight);
+        _hist_m_125->fill(M, weight);
+      } else if (jets[0].momentum().pT()>90) {
+        _hist_T_90->fill(T, weight);
+        _hist_m_90->fill(M, weight);
+      }
+    }
+
+    void finalize() {
+      normalize(_hist_T_90);
+      normalize(_hist_m_90);
+      normalize(_hist_T_125);
+      normalize(_hist_m_125);
+      normalize(_hist_T_200);
+      normalize(_hist_m_200);
+    }
+
+  private:
+
+    AIDA::IHistogram1D *_hist_T_90;
+    AIDA::IHistogram1D *_hist_m_90;
+    AIDA::IHistogram1D *_hist_T_125;
+    AIDA::IHistogram1D *_hist_m_125;
+    AIDA::IHistogram1D *_hist_T_200;
+    AIDA::IHistogram1D *_hist_m_200;
+  };
+
+  // This global object acts as a hook for the plugin system
+  AnalysisBuilder<CMS_2011_S8957746> plugin_CMS_2011_S8957746;
+}

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Tue Feb 15 17:44:00 2011	(r2932)
+++ trunk/src/Analyses/Makefile.am	Tue Feb 15 17:51:37 2011	(r2933)
@@ -59,6 +59,9 @@
     ATLAS_2011_S8971293.cc \
     ATLAS_2010_CONF_2010_083.cc
 
+lib_LTLIBRARIES += RivetCMSAnalyses.la
+RivetCMSAnalyses_la_SOURCES = \
+    CMS_2011_S8957746.cc
 
 lib_LTLIBRARIES += RivetCDFAnalyses.la
 RivetCDFAnalyses_la_SOURCES = \


More information about the Rivet-svn mailing list