rivet is hosted by Hepforge, IPPP Durham

## Rivet analyses reference

### ATLAS_2015_I1386475

Centrality dependent of pseudorapidity distributions in p-Pb collisions $\sqrt{s_{_\text{NN}}}=5.02$ TeV.
Experiment: ATLAS (LHC)
Inspire ID: 1386475
Status: UNVALIDATED
Authors:
• Christian Bierlich
References:
• Eur.Phys.J. C76 (2016) no.4, 199
• DOI:0.1140/epjc/s10052-016-4002-3
• arXiv: 1508.00848
Beams: p+ 1000822080
Beam energies: (4000.0, 326560.0) GeV
Run details:
• p-Pb minimum bias events with $c\tau > 10 cm$ stable.

The centrality dependence of the mean charged-particle multiplicity as a function of pseudorapidity in proton-lead collisions. The collision energy is asymmetric, 4000 GeV for the proton and 1570/N for the lead. The primary particle definition is $c\tau > 10 cm$ stable, and should be enforced at generator level.

Source code: ATLAS_2015_I1386475.cc
  1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Tools/AtlasCommon.hh" namespace Rivet { /// @brief Add a short analysis description here class ATLAS_2015_I1386475 : public Analysis { public: /// Constructor DEFAULT_RIVET_ANALYSIS_CTOR(ATLAS_2015_I1386475); /// Book histograms and initialise projections before the run void init() { // The centrality projection. declareCentrality(ATLAS::SumET_PB_Centrality(), "ATLAS_pPb_Calib", "SumETPb", "CENT"); // The trigger projection. declare(ATLAS::MinBiasTrigger(), "Trigger"); // The particles to be analysed. declare(ChargedFinalState(Cuts::eta > -2.7 && Cuts::eta < 2.7 && Cuts::pT > 0.1*GeV), "CFS"); // The centrality bins' upper edges. centralityBins = {90., 60., 40., 30., 20., 10., 5., 1.}; for (int i = 0; i < 8; ++i) { histEta[centralityBins[i]] = bookHisto1D(2, 1, i + 1); sow[centralityBins[i]] = 0.; } } /// Perform the per-event analysis void analyze(const Event& event) { const double weight = event.weight(); // Apply event triggers. if ( !apply(event, "Trigger")() ) vetoEvent; // We must have direct acces to the centrality projection. const CentralityProjection& cent = apply(event,"CENT"); double c = cent(); // Find the correct centrality histogram auto hItr = histEta.upper_bound(c); if (hItr == histEta.end()) return; // Find the correct sow. auto sItr = sow.upper_bound(c); if (sItr == sow.end()) return; sItr->second += weight; for ( const auto &p : apply(event,"CFS").particles() ) hItr->second->fill(p.eta(), weight); } /// Finalize void finalize() { // Scale by the inverse sum of event weights in each centrality // bin. for (int i = 0; i < 8; ++i) histEta[centralityBins[i]]->scaleW(1./sow[centralityBins[i]]); } private: /// The histograms binned in centrality. vector centralityBins; map histEta; map sow; }; // The hook for the plugin system DECLARE_RIVET_PLUGIN(ATLAS_2015_I1386475); }