1 | #ifndef DelphesLHEFReader_h
|
---|
2 | #define DelphesLHEFReader_h
|
---|
3 |
|
---|
4 | /** \class DelphesLHEFReader
|
---|
5 | *
|
---|
6 | * Reads LHEF file
|
---|
7 | *
|
---|
8 | *
|
---|
9 | * $Date$
|
---|
10 | * $Revision$
|
---|
11 | *
|
---|
12 | *
|
---|
13 | * \author P. Demin - UCL, Louvain-la-Neuve
|
---|
14 | *
|
---|
15 | */
|
---|
16 |
|
---|
17 | #include <stdio.h>
|
---|
18 |
|
---|
19 | class TObjArray;
|
---|
20 | class TStopwatch;
|
---|
21 | class TDatabasePDG;
|
---|
22 | class ExRootTreeBranch;
|
---|
23 | class DelphesFactory;
|
---|
24 |
|
---|
25 | class DelphesLHEFReader
|
---|
26 | {
|
---|
27 | public:
|
---|
28 |
|
---|
29 | DelphesLHEFReader();
|
---|
30 | ~DelphesLHEFReader();
|
---|
31 |
|
---|
32 | void SetInputFile(FILE *inputFile);
|
---|
33 |
|
---|
34 | void Clear();
|
---|
35 | bool EventReady();
|
---|
36 |
|
---|
37 | bool ReadBlock(DelphesFactory *factory,
|
---|
38 | TObjArray *allParticleOutputArray,
|
---|
39 | TObjArray *stableParticleOutputArray,
|
---|
40 | TObjArray *partonOutputArray);
|
---|
41 |
|
---|
42 | void AnalyzeEvent(ExRootTreeBranch *branch, long long eventNumber,
|
---|
43 | TStopwatch *readStopWatch, TStopwatch *procStopWatch);
|
---|
44 |
|
---|
45 | private:
|
---|
46 |
|
---|
47 | void AnalyzeParticle(DelphesFactory *factory,
|
---|
48 | TObjArray *allParticleOutputArray,
|
---|
49 | TObjArray *stableParticleOutputArray,
|
---|
50 | TObjArray *partonOutputArray);
|
---|
51 |
|
---|
52 | FILE *fInputFile;
|
---|
53 |
|
---|
54 | char *fBuffer;
|
---|
55 |
|
---|
56 | TDatabasePDG *fPDG;
|
---|
57 |
|
---|
58 | int fEventCounter;
|
---|
59 |
|
---|
60 | int fParticleCounter, fProcessID;
|
---|
61 | double fWeight, fScalePDF, fAlphaQCD, fAlphaQED;
|
---|
62 |
|
---|
63 | int fPID, fStatus, fM1, fM2, fC1, fC2;
|
---|
64 | double fPx, fPy, fPz, fE, fMass;
|
---|
65 | };
|
---|
66 |
|
---|
67 | #endif // DelphesLHEFReader_h
|
---|
68 |
|
---|
69 |
|
---|