]>
Commit | Line | Data |
---|---|---|
9419eeef | 1 | // main01.cc is a part of the PYTHIA event generator. |
2 | // Copyright (C) 2010 Torbjorn Sjostrand. | |
3 | // PYTHIA is licenced under the GNU GPL version 2, see COPYING for details. | |
4 | // Please respect the MCnet Guidelines, see GUIDELINES for details. | |
5 | ||
6 | // This is a simple test program. It fits on one slide in a talk. | |
7 | // It studies the charged multiplicity distribution at the LHC. | |
8 | ||
9 | #include "Pythia.h" | |
10 | using namespace Pythia8; | |
11 | int main() { | |
12 | // Generator. Process selection. LHC initialization. Histogram. | |
13 | Pythia pythia; | |
14 | pythia.readString("HardQCD:all = on"); | |
15 | pythia.readString("PhaseSpace:pTHatMin = 20."); | |
16 | pythia.init( 2212, 2212, 14000.); | |
17 | Hist mult("charged multiplicity", 100, -0.5, 799.5); | |
18 | // Begin event loop. Generate event. Skip if error. List first one. | |
19 | for (int iEvent = 0; iEvent < 100; ++iEvent) { | |
20 | if (!pythia.next()) continue; | |
21 | if (iEvent < 1) {pythia.info.list(); pythia.event.list();} | |
22 | // Find number of all final charged particles and fill histogram. | |
23 | int nCharged = 0; | |
24 | for (int i = 0; i < pythia.event.size(); ++i) | |
25 | if (pythia.event[i].isFinal() && pythia.event[i].isCharged()) | |
26 | ++nCharged; | |
27 | mult.fill( nCharged ); | |
28 | // End of event loop. Statistics. Histogram. Done. | |
29 | } | |
30 | pythia.statistics(); | |
31 | cout << mult; | |
32 | return 0; | |
33 | } |