|
[Rivet-svn] rivet: New analyses: 1 ALICE, 2 ATLAS, 5 CMS, 2 LHCFRivet Mercurial rivet at projects.hepforge.orgWed 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 |