Fork me on GitHub

source: git/classes/DelphesHepMC3Reader.h@ c6f9311

Last change on this file since c6f9311 was b4786d3, checked in by Pavel Demin <pavel.demin@…>, 3 years ago

add DelphesHepMC3

  • Property mode set to 100644
File size: 2.5 KB
Line 
1/*
2 * Delphes: a framework for fast simulation of a generic collider experiment
3 * Copyright (C) 2012-2021 Universite catholique de Louvain (UCL), Belgium
4 *
5 * This program is free software: you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation, either version 3 of the License, or
8 * (at your option) any later version.
9 *
10 * This program is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU General Public License for more details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with this program. If not, see <http://www.gnu.org/licenses/>.
17 */
18
19#ifndef DelphesHepMC3Reader_h
20#define DelphesHepMC3Reader_h
21
22/** \class DelphesHepMC3Reader
23 *
24 * Reads HepMC file
25 *
26 * \author P. Demin - UCL, Louvain-la-Neuve
27 *
28 */
29
30#include <map>
31#include <vector>
32
33#include <stdio.h>
34
35class TObjArray;
36class TStopwatch;
37class TDatabasePDG;
38class ExRootTreeBranch;
39class DelphesFactory;
40
41class DelphesHepMC3Reader
42{
43public:
44 DelphesHepMC3Reader();
45 ~DelphesHepMC3Reader();
46
47 void SetInputFile(FILE *inputFile);
48
49 void Clear();
50 bool EventReady();
51
52 bool ReadBlock(DelphesFactory *factory,
53 TObjArray *allParticleOutputArray,
54 TObjArray *stableParticleOutputArray,
55 TObjArray *partonOutputArray);
56
57 void AnalyzeEvent(ExRootTreeBranch *branch, long long eventNumber,
58 TStopwatch *readStopWatch, TStopwatch *procStopWatch);
59
60 void AnalyzeWeight(ExRootTreeBranch *branch);
61
62private:
63 void AnalyzeParticle(DelphesFactory *factory,
64 TObjArray *allParticleOutputArray,
65 TObjArray *stableParticleOutputArray,
66 TObjArray *partonOutputArray);
67
68 void FinalizeParticles(TObjArray *allParticleOutputArray);
69
70 FILE *fInputFile;
71
72 char *fBuffer;
73
74 TDatabasePDG *fPDG;
75
76 int fEventNumber, fMPI, fProcessID, fSignalCode, fVertexCounter, fParticleCounter;
77 double fScale, fAlphaQCD, fAlphaQED;
78
79 double fMomentumCoefficient, fPositionCoefficient;
80
81 std::vector<double> fWeight;
82
83 double fCrossSection, fCrossSectionError;
84
85 int fID1, fID2;
86 double fX1, fX2, fScalePDF, fPDF1, fPDF2;
87
88 int fVertexCode, fVertexStatus;
89 double fX, fY, fZ, fT;
90
91 int fParticleCode, fPID, fParticleStatus, fOutVertexCode;
92 double fPx, fPy, fPz, fE, fMass;
93
94 int fM1, fM2;
95
96 std::map<int, int> fVertexMap;
97 std::map<int, std::pair<int, int> > fDaughterMap;
98};
99
100#endif // DelphesHepMC3Reader_h
Note: See TracBrowser for help on using the repository browser.