|
|
Rivet analyses reference
UA5_1989_S1926373
UA5 charged multiplicity measurements
Experiment: UA5 (CERN SPS)
Inspire ID: 267179
Status: VALIDATED
Authors:
- Holger Schulz
- Christophe L. J. Vaillant
- Andy Buckley
References:
Beams: p- p+
Beam energies: (100.0, 100.0); (450.0, 450.0) GeV
Run details:
- Minimum bias events at $\sqrt{s} = 200$ and 900 GeV. Enable single and double diffractive events in addition to non-diffractive processes. Beam energy must be specified as analysis option "ENERGY" (200 or 900 GeV) when rivet-merge'ing samples.
Multiplicity distributions of charged particles produced in non-single-diffractive collisions between protons and antiprotons at centre-of-mass energies of 200 and 900 GeV. The data were recorded in the UA5 streamer chambers at the CERN collider, which was operated in a pulsed mode between the two energies. This analysis confirms the violation of KNO scaling in full phase space found by the UA5 group at an energy of 546 GeV, with similar measurements at 200 and 900 GeV. Beam energy must be specified as analysis option "ENERGY" (200 or 900 GeV) when rivet-merging samples.
Source code:
UA5_1989_S1926373.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
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106 | // -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/TriggerUA5.hh"
namespace Rivet {
/// UA5 min bias charged multiplicities in central \f$ \eta \f$ ranges
class UA5_1989_S1926373 : public Analysis {
public:
/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(UA5_1989_S1926373);
/// @name Analysis methods
/// @{
/// Book histograms and projections
void init() {
declare(TriggerUA5(), "Trigger");
declare(ChargedFinalState((Cuts::etaIn(-0.5, 0.5))), "CFS05");
declare(ChargedFinalState((Cuts::etaIn(-1.5, 1.5))), "CFS15");
declare(ChargedFinalState((Cuts::etaIn(-3.0, 3.0))), "CFS30");
declare(ChargedFinalState((Cuts::etaIn(-5.0, 5.0))), "CFS50");
// NB. _hist_nch and _hist_ncheta50 use the same data but different binning
if (isCompatibleWithSqrtS(200)) {
book(_hist_nch ,1, 1, 1);
book(_hist_nch_eta05 ,3, 1, 1);
book(_hist_nch_eta15 ,4, 1, 1);
book(_hist_nch_eta30 ,5, 1, 1);
book(_hist_nch_eta50 ,6, 1, 1);
book(_hist_mean_nch ,11, 1, 1);
} else if (isCompatibleWithSqrtS(900)) {
book(_hist_nch ,2, 1, 1);
book(_hist_nch_eta05 ,7, 1, 1);
book(_hist_nch_eta15 ,8, 1, 1);
book(_hist_nch_eta30 ,9, 1, 1);
book(_hist_nch_eta50 ,10, 1, 1);
book(_hist_mean_nch ,12, 1, 1);
}
book(_sumWPassed, "SumW");
/// @todo Moments of distributions
}
/// Do the analysis
void analyze(const Event& event) {
// Trigger
const TriggerUA5& trigger = apply<TriggerUA5>(event, "Trigger");
if (!trigger.nsdDecision()) vetoEvent;
_sumWPassed->fill();
// Count final state particles in several eta regions
const int numP05 = apply<ChargedFinalState>(event, "CFS05").size();
const int numP15 = apply<ChargedFinalState>(event, "CFS15").size();
const int numP30 = apply<ChargedFinalState>(event, "CFS30").size();
const int numP50 = apply<ChargedFinalState>(event, "CFS50").size();
// Fill histograms
_hist_nch->fill(numP50);
_hist_nch_eta05->fill(numP05);
_hist_nch_eta15->fill(numP15);
_hist_nch_eta30->fill(numP30);
_hist_nch_eta50->fill(numP50);
_hist_mean_nch->fill(_hist_mean_nch->bin(0).xMid(), numP50);
}
void finalize() {
scale(_hist_nch, 1.0 / *_sumWPassed);
scale(_hist_nch_eta05, 1.0 / *_sumWPassed);
scale(_hist_nch_eta15, 1.0 / *_sumWPassed);
scale(_hist_nch_eta30, 1.0 / *_sumWPassed);
scale(_hist_nch_eta50, 1.0 / *_sumWPassed);
scale(_hist_mean_nch, 1.0 / *_sumWPassed);
}
/// @}
private:
/// Weight counter
CounterPtr _sumWPassed;
/// @name Histograms
/// @{
Histo1DPtr _hist_nch;
Histo1DPtr _hist_nch_eta05;
Histo1DPtr _hist_nch_eta15;
Histo1DPtr _hist_nch_eta30;
Histo1DPtr _hist_nch_eta50;
Histo1DPtr _hist_mean_nch;
/// @}
};
RIVET_DECLARE_ALIASED_PLUGIN(UA5_1989_S1926373, UA5_1989_I267179);
}
|
|