[Rivet-svn] rivet: New analyses: 1 ALICE, 2 ATLAS, 5 CMS, 2 LHCF

Rivet Mercurial rivet at projects.hepforge.org
Wed May 23 17:15:03 BST 2018


details:   https://rivet.hepforge.org/hg/rivet/rev/1e0935647aa3
branches:  release-2-6-x
changeset: 6304:1e0935647aa3
user:      Jon Butterworth <j.butterworth at cern.ch>
date:      Wed May 23 17:03:15 2018 +0100
description:
New analyses: 1 ALICE, 2 ATLAS, 5 CMS, 2 LHCF

diffs (truncated from 67129 to 50 lines):

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/analyses/pluginALICE/ALICE_2017_I1620477.cc	Wed May 23 17:03:15 2018 +0100
@@ -0,0 +1,90 @@
+//-*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/UnstableFinalState.hh"
+#include "Rivet/Tools/ParticleUtils.hh"
+
+namespace Rivet {
+
+
+  class ALICE_2017_I1620477 : public Analysis {
+  public:
+
+    /// Constructor
+    ALICE_2017_I1620477()
+      : Analysis("ALICE_2017_I1620477"),
+        _rapmax(0.8)
+    {    }
+
+
+    void init() {
+
+      const UnstableFinalState ufs(Cuts::absrap < _rapmax);
+      addProjection(ufs, "UFS");
+
+      _h_pi0 = bookHisto1D(1,1,1);
+      _h_eta = bookHisto1D(2,1,1);
+      _h_etaToPion = bookScatter2D(8,1,1);
+
+      // temporary plots with the binning of _h_etaToPion
+      // to construct the eta/pi0 ratio in the end
+      _temp_h_pion = bookHisto1D("TMP/h_pion",refData(8,1,1));
+      _temp_h_eta = bookHisto1D("TMP/h_eta",refData(8,1,1));
+    }
+
+
+    void analyze(const Event& event) {
+
+      const double weight = event.weight();
+      const UnstableFinalState& ufs = applyProjection<UnstableFinalState>(event, "UFS");
+
+      for(auto p: ufs.particles()) {
+
+        if (p.pid() == 111) {
+          // neutral pion; ALICE corrects for pi0 feed-down
+          if ( !(p.hasAncestor(310)  || p.hasAncestor(130)   || // K0_s, K0_l
+                 p.hasAncestor(321)  || p.hasAncestor(-321)  || // K+,K-
+                 p.hasAncestor(3122) || p.hasAncestor(-3122) || // Lambda, Anti-Lambda
+                 p.hasAncestor(3212) || p.hasAncestor(-3212) || // Sigma0


More information about the Rivet-svn mailing list