Fork me on GitHub

source: git/classes/DelphesClasses.h@ 5a69ba1d

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 5a69ba1d was 5a69ba1d, checked in by Ulrike Schnoor <schnooru@…>, 7 years ago

Merge branch 'master' of github.com:delphes/delphes

  • Property mode set to 100644
File size: 19.1 KB
Line 
1/*
2 * Delphes: a framework for fast simulation of a generic collider experiment
3 * Copyright (C) 2012-2014 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 DelphesClasses_h
20#define DelphesClasses_h
21
22/**
23 *
24 * Definition of classes to be stored in the root tree.
25 * Function CompareXYZ sorts objects by the variable XYZ that MUST be
26 * present in the data members of the root tree class of the branch.
27 *
28 * \author P. Demin - UCL, Louvain-la-Neuve
29 *
30 */
31
32// Dependencies (#includes)
33
34#include "TRef.h"
35#include "TObject.h"
36#include "TRefArray.h"
37#include "TLorentzVector.h"
38
39#include "classes/SortableObject.h"
40
41class DelphesFactory;
42
43//---------------------------------------------------------------------------
44
45class Event: public TObject
46{
47public:
48
49 Long64_t Number; // event number
50
51 Float_t ReadTime; // read time
52 Float_t ProcTime; // processing time
53
54 ClassDef(Event, 1)
55};
56
57//---------------------------------------------------------------------------
58
59class LHCOEvent: public Event
60{
61public:
62
63 Int_t Trigger; // trigger word
64
65 ClassDef(LHCOEvent, 1)
66};
67
68//---------------------------------------------------------------------------
69
70class LHEFEvent: public Event
71{
72public:
73
74 Int_t ProcessID; // subprocess code for the event | hepup.IDPRUP
75
76 Float_t Weight; // weight for the event | hepup.XWGTUP
77 Float_t ScalePDF; // scale in GeV used in the calculation of the PDFs in the event | hepup.SCALUP
78 Float_t AlphaQED; // value of the QED coupling used in the event | hepup.AQEDUP
79 Float_t AlphaQCD; // value of the QCD coupling used in the event | hepup.AQCDUP
80
81 ClassDef(LHEFEvent, 2)
82};
83
84//---------------------------------------------------------------------------
85
86class LHEFWeight: public TObject
87{
88public:
89 Int_t ID; // weight ID
90 Float_t Weight; // weight value
91
92 ClassDef(LHEFWeight, 1)
93};
94
95//---------------------------------------------------------------------------
96
97class HepMCEvent: public Event
98{
99public:
100
101 Int_t ProcessID; // unique signal process id | signal_process_id()
102 Int_t MPI; // number of multi parton interactions | mpi ()
103
104 Float_t Weight; // weight for the event
105
106 Float_t Scale; // energy scale, see hep-ph/0109068 | event_scale()
107 Float_t AlphaQED; // QED coupling, see hep-ph/0109068 | alphaQED()
108 Float_t AlphaQCD; // QCD coupling, see hep-ph/0109068 | alphaQCD()
109
110 Int_t ID1; // flavour code of first parton | pdf_info()->id1()
111 Int_t ID2; // flavour code of second parton | pdf_info()->id2()
112
113 Float_t X1; // fraction of beam momentum carried by first parton ("beam side") | pdf_info()->x1()
114 Float_t X2; // fraction of beam momentum carried by second parton ("target side") | pdf_info()->x2()
115
116 Float_t ScalePDF; // Q-scale used in evaluation of PDF's (in GeV) | pdf_info()->scalePDF()
117
118 Float_t PDF1; // PDF (id1, x1, Q) | pdf_info()->pdf1()
119 Float_t PDF2; // PDF (id2, x2, Q) | pdf_info()->pdf2()
120
121 ClassDef(HepMCEvent, 2)
122};
123
124//---------------------------------------------------------------------------
125
126class GenParticle: public SortableObject
127{
128public:
129 Int_t PID; // particle HEP ID number | hepevt.idhep[number]
130
131 Int_t Status; // particle status | hepevt.isthep[number]
132 Int_t IsPU; // 0 or 1 for particles from pile-up interactions
133
134 Int_t M1; // particle 1st mother | hepevt.jmohep[number][0] - 1
135 Int_t M2; // particle 2nd mother | hepevt.jmohep[number][1] - 1
136
137 Int_t D1; // particle 1st daughter | hepevt.jdahep[number][0] - 1
138 Int_t D2; // particle last daughter | hepevt.jdahep[number][1] - 1
139
140 Int_t Charge; // particle charge
141
142 Float_t Mass; // particle mass
143
144 Float_t E; // particle energy | hepevt.phep[number][3]
145 Float_t Px; // particle momentum vector (x component) | hepevt.phep[number][0]
146 Float_t Py; // particle momentum vector (y component) | hepevt.phep[number][1]
147 Float_t Pz; // particle momentum vector (z component) | hepevt.phep[number][2]
148
149 Float_t P; // particle momentum
150 Float_t PT; // particle transverse momentum
151 Float_t Eta; // particle pseudorapidity
152 Float_t Phi; // particle azimuthal angle
153
154 Float_t Rapidity; // particle rapidity
155 Float_t CtgTheta; // particle cotangent of theta
156
157 Float_t D0; // particle transverse impact parameter
158 Float_t DZ; // particle longitudinal impact parameter
159
160 Float_t T; // particle vertex position (t component) | hepevt.vhep[number][3]
161 Float_t X; // particle vertex position (x component) | hepevt.vhep[number][0]
162 Float_t Y; // particle vertex position (y component) | hepevt.vhep[number][1]
163 Float_t Z; // particle vertex position (z component) | hepevt.vhep[number][2]
164
165 static CompBase *fgCompare; //!
166 const CompBase *GetCompare() const { return fgCompare; }
167
168 TLorentzVector P4() const;
169
170 ClassDef(GenParticle, 2)
171};
172
173//---------------------------------------------------------------------------
174
175class Vertex: public SortableObject
176{
177public:
178
179 Float_t T; // vertex position (t component)
180 Float_t X; // vertex position (x component)
181 Float_t Y; // vertex position (y component)
182 Float_t Z; // vertex position (z component)
183
184 Double_t ErrorT; // vertex position error (t component)
185 Double_t ErrorX; // vertex position error (x component)
186 Double_t ErrorY; // vertex position error (y component)
187 Double_t ErrorZ; // vertex position error (z component)
188
189 Int_t Index; // vertex index
190 Int_t NDF; // number of degrees of freedom
191
192 Double_t Sigma; // vertex position (z component) error
193 Double_t SumPT2; // sum pt^2 of tracks attached to the vertex
194 Double_t GenSumPT2; // sum pt^2 of gen tracks attached to the vertex
195
196 Double_t GenDeltaZ; // distance in z to closest generated vertex
197 Double_t BTVSumPT2; // sum pt^2 of tracks attached to the secondary vertex
198
199 TRefArray Constituents; // references to constituents
200
201 static CompBase *fgCompare; //!
202 const CompBase *GetCompare() const { return fgCompare; }
203
204 ClassDef(Vertex, 3)
205};
206
207//---------------------------------------------------------------------------
208
209class MissingET: public TObject
210{
211public:
212 Float_t MET; // mising transverse energy
213 Float_t Eta; // mising energy pseudorapidity
214 Float_t Phi; // mising energy azimuthal angle
215
216 TLorentzVector P4() const;
217
218 ClassDef(MissingET, 1)
219};
220
221//---------------------------------------------------------------------------
222
223class ScalarHT: public TObject
224{
225public:
226 Float_t HT; // scalar sum of transverse momenta
227
228 ClassDef(ScalarHT, 1)
229};
230
231//---------------------------------------------------------------------------
232
233class Rho: public TObject
234{
235public:
236 Float_t Rho; // rho energy density
237 Float_t Edges[2]; // pseudorapidity range edges
238
239 ClassDef(Rho, 1)
240};
241
242//---------------------------------------------------------------------------
243
244class Weight: public TObject
245{
246public:
247 Float_t Weight; // weight for the event
248
249 ClassDef(Weight, 1)
250};
251
252//---------------------------------------------------------------------------
253
254class Photon: public SortableObject
255{
256public:
257
258 Float_t PT; // photon transverse momentum
259 Float_t Eta; // photon pseudorapidity
260 Float_t Phi; // photon azimuthal angle
261
262 Float_t E; // photon energy
263
264 Float_t T; // particle arrival time of flight
265
266 Float_t EhadOverEem; // ratio of the hadronic versus electromagnetic energy deposited in the calorimeter
267
268 TRefArray Particles; // references to generated particles
269
270 Float_t IsolationVar; // isolation variable
271 Float_t IsolationVarRhoCorr; // isolation variable
272 Float_t SumPtCharged; // isolation variable
273 Float_t SumPtNeutral; // isolation variable
274 Float_t SumPtChargedPU; // isolation variable
275 Float_t SumPt; // isolation variable
276
277 static CompBase *fgCompare; //!
278 const CompBase *GetCompare() const { return fgCompare; }
279
280 TLorentzVector P4() const;
281
282 ClassDef(Photon, 3)
283};
284
285//---------------------------------------------------------------------------
286
287class Electron: public SortableObject
288{
289public:
290
291 Float_t PT; // electron transverse momentum
292 Float_t Eta; // electron pseudorapidity
293 Float_t Phi; // electron azimuthal angle
294
295 Float_t T; // particle arrival time of flight
296
297 Int_t Charge; // electron charge
298
299 Float_t EhadOverEem; // ratio of the hadronic versus electromagnetic energy deposited in the calorimeter
300
301 TRef Particle; // reference to generated particle
302
303 Float_t IsolationVar; // isolation variable
304 Float_t IsolationVarRhoCorr; // isolation variable
305 Float_t SumPtCharged; // isolation variable
306 Float_t SumPtNeutral; // isolation variable
307 Float_t SumPtChargedPU; // isolation variable
308 Float_t SumPt; // isolation variable
309
310 static CompBase *fgCompare; //!
311 const CompBase *GetCompare() const { return fgCompare; }
312
313 TLorentzVector P4() const;
314
315 ClassDef(Electron, 3)
316};
317
318//---------------------------------------------------------------------------
319
320class Muon: public SortableObject
321{
322public:
323
324 Float_t PT; // muon transverse momentum
325 Float_t Eta; // muon pseudorapidity
326 Float_t Phi; // muon azimuthal angle
327
328 Float_t T; // particle arrival time of flight
329
330 Int_t Charge; // muon charge
331
332 TRef Particle; // reference to generated particle
333
334 Float_t IsolationVar; // isolation variable
335 Float_t IsolationVarRhoCorr; // isolation variable
336 Float_t SumPtCharged; // isolation variable
337 Float_t SumPtNeutral; // isolation variable
338 Float_t SumPtChargedPU; // isolation variable
339 Float_t SumPt; // isolation variable
340
341 static CompBase *fgCompare; //!
342 const CompBase *GetCompare() const { return fgCompare; }
343
344 TLorentzVector P4() const;
345
346 ClassDef(Muon, 3)
347};
348
349//---------------------------------------------------------------------------
350
351class Jet: public SortableObject
352{
353public:
354
355 Float_t PT; // jet transverse momentum
356 Float_t Eta; // jet pseudorapidity
357 Float_t Phi; // jet azimuthal angle
358
359 Float_t T; //particle arrival time of flight
360
361 Float_t Mass; // jet invariant mass
362
363 Float_t DeltaEta; // jet radius in pseudorapidity
364 Float_t DeltaPhi; // jet radius in azimuthal angle
365
366 UInt_t Flavor; // jet flavor
367 UInt_t FlavorAlgo; // jet flavor
368 UInt_t FlavorPhys; // jet flavor
369
370 UInt_t BTag; // 0 or 1 for a jet that has been tagged as containing a heavy quark
371 UInt_t BTagAlgo; // 0 or 1 for a jet that has been tagged as containing a heavy quark
372 UInt_t BTagPhys; // 0 or 1 for a jet that has been tagged as containing a heavy quark
373
374 UInt_t TauTag; // 0 or 1 for a jet that has been tagged as a tau
375
376 Int_t Charge; // tau charge
377
378 Float_t EhadOverEem; // ratio of the hadronic versus electromagnetic energy deposited in the calorimeter
379
380 Int_t NCharged; // number of charged constituents
381 Int_t NNeutrals; // number of neutral constituents
382 Float_t Beta; // (sum pt of charged pile-up constituents)/(sum pt of charged constituents)
383 Float_t BetaStar; // (sum pt of charged constituents coming from hard interaction)/(sum pt of charged constituents)
384 Float_t MeanSqDeltaR; // average distance (squared) between constituent and jet weighted by pt (squared) of constituent
385 Float_t PTD; // average pt between constituent and jet weighted by pt of constituent
386 Float_t FracPt[5]; // (sum pt of constituents within a ring 0.1*i < DeltaR < 0.1*(i+1))/(sum pt of constituents)
387
388 Float_t Tau[5]; // N-subjettiness
389
390 TLorentzVector SoftDroppedJet;
391 TLorentzVector SoftDroppedSubJet1;
392 TLorentzVector SoftDroppedSubJet2;
393
394 TLorentzVector TrimmedP4[5]; // first entry (i = 0) is the total Trimmed Jet 4-momenta and from i = 1 to 4 are the trimmed subjets 4-momenta
395 TLorentzVector PrunedP4[5]; // first entry (i = 0) is the total Pruned Jet 4-momenta and from i = 1 to 4 are the pruned subjets 4-momenta
396 TLorentzVector SoftDroppedP4[5]; // first entry (i = 0) is the total SoftDropped Jet 4-momenta and from i = 1 to 4 are the pruned subjets 4-momenta
397
398
399 Int_t NSubJetsTrimmed; // number of subjets trimmed
400 Int_t NSubJetsPruned; // number of subjets pruned
401 Int_t NSubJetsSoftDropped; // number of subjets soft-dropped
402
403 Double_t ExclYmerge23;
404 Double_t ExclYmerge34;
405 Double_t ExclYmerge45;
406 Double_t ExclYmerge56;
407
408 TRefArray Constituents; // references to constituents
409 TRefArray Particles; // references to generated particles
410
411 static CompBase *fgCompare; //!
412 const CompBase *GetCompare() const { return fgCompare; }
413
414 TLorentzVector P4() const;
415 TLorentzVector Area;
416
417 ClassDef(Jet, 3)
418};
419
420//---------------------------------------------------------------------------
421
422class Track: public SortableObject
423{
424public:
425 Int_t PID; // HEP ID number
426
427 Int_t Charge; // track charge
428
429 Float_t P; // track momentum
430 Float_t PT; // track transverse momentum
431 Float_t Eta; // track pseudorapidity
432 Float_t Phi; // track azimuthal angle
433 Float_t CtgTheta; // track cotangent of theta
434
435 Float_t EtaOuter; // track pseudorapidity at the tracker edge
436 Float_t PhiOuter; // track azimuthal angle at the tracker edge
437
438 Float_t T; // track vertex position (t component)
439 Float_t X; // track vertex position (x component)
440 Float_t Y; // track vertex position (y component)
441 Float_t Z; // track vertex position (z component)
442
443 Float_t TOuter; // track position (t component) at the tracker edge
444 Float_t XOuter; // track position (x component) at the tracker edge
445 Float_t YOuter; // track position (y component) at the tracker edge
446 Float_t ZOuter; // track position (z component) at the tracker edge
447
448 Float_t Xd; // X coordinate of point of closest approach to vertex
449 Float_t Yd; // Y coordinate of point of closest approach to vertex
450 Float_t Zd; // Z coordinate of point of closest approach to vertex
451
452 Float_t L; // track path length
453 Float_t D0; // track transverse impact parameter
454 Float_t DZ; // track longitudinal impact parameter
455
456 Float_t ErrorP; // track momentum error
457 Float_t ErrorPT; // track transverse momentum error
458 Float_t ErrorPhi; // track azimuthal angle error
459 Float_t ErrorCtgTheta; // track cotangent of theta error
460
461 Float_t ErrorT; // time measurement error
462 Float_t ErrorD0; // track transverse impact parameter error
463 Float_t ErrorDZ; // track longitudinal impact parameter error
464
465 TRef Particle; // reference to generated particle
466
467 Int_t VertexIndex; // reference to vertex
468
469 static CompBase *fgCompare; //!
470 const CompBase *GetCompare() const { return fgCompare; }
471
472 TLorentzVector P4() const;
473
474 ClassDef(Track, 3)
475};
476
477//---------------------------------------------------------------------------
478
479class Tower: public SortableObject
480{
481public:
482 Float_t ET; // calorimeter tower transverse energy
483 Float_t Eta; // calorimeter tower pseudorapidity
484 Float_t Phi; // calorimeter tower azimuthal angle
485
486 Float_t E; // calorimeter tower energy
487
488 Float_t T; // ecal deposit time, averaged by sqrt(EM energy) over all particles, not smeared
489 Int_t NTimeHits; // number of hits contributing to time measurement
490
491 Float_t Eem; // calorimeter tower electromagnetic energy
492 Float_t Ehad; // calorimeter tower hadronic energy
493
494 Float_t Edges[4]; // calorimeter tower edges
495
496 TRefArray Particles; // references to generated particles
497
498 static CompBase *fgCompare; //!
499 const CompBase *GetCompare() const { return fgCompare; }
500
501 TLorentzVector P4() const;
502
503 ClassDef(Tower, 2)
504};
505
506//---------------------------------------------------------------------------
507
508class HectorHit: public SortableObject
509{
510public:
511 Float_t E; // reconstructed energy [GeV]
512
513 Float_t Tx; // angle of the momentum in the horizontal (x,z) plane [urad]
514 Float_t Ty; // angle of the momentum in the verical (y,z) plane [urad]
515
516 Float_t T; // time of flight to the detector [s]
517
518 Float_t X; // horizontal distance to the beam [um]
519 Float_t Y; // vertical distance to the beam [um]
520 Float_t S; // distance to the interaction point [m]
521
522 TRef Particle; // reference to generated particle
523
524 static CompBase *fgCompare; //!
525 const CompBase *GetCompare() const { return fgCompare; }
526
527 ClassDef(HectorHit, 1)
528};
529
530//---------------------------------------------------------------------------
531
532class Candidate: public SortableObject
533{
534 friend class DelphesFactory;
535
536public:
537 Candidate();
538
539 Int_t PID;
540
541 Int_t Status;
542 Int_t M1, M2, D1, D2;
543
544 Int_t Charge;
545
546 Float_t Mass;
547
548 Int_t IsPU;
549 Int_t IsRecoPU;
550
551 Int_t IsConstituent;
552
553 Int_t IsFromConversion;
554
555 UInt_t Flavor;
556 UInt_t FlavorAlgo;
557 UInt_t FlavorPhys;
558
559 UInt_t BTag;
560 UInt_t BTagAlgo;
561 UInt_t BTagPhys;
562
563 UInt_t TauTag;
564
565 Float_t Eem;
566 Float_t Ehad;
567
568 Float_t Edges[4];
569 Float_t DeltaEta;
570 Float_t DeltaPhi;
571
572 TLorentzVector Momentum, Position, InitialPosition, PositionError, Area;
573
574 Float_t L; // path length
575 Float_t ErrorT; // path length
576 Float_t D0;
577 Float_t ErrorD0;
578 Float_t DZ;
579 Float_t ErrorDZ;
580 Float_t P;
581 Float_t ErrorP;
582 Float_t PT;
583 Float_t ErrorPT;
584 Float_t CtgTheta;
585 Float_t ErrorCtgTheta;
586 Float_t Phi;
587 Float_t ErrorPhi;
588
589 Float_t Xd;
590 Float_t Yd;
591 Float_t Zd;
592
593 // tracking resolution
594
595 Float_t TrackResolution;
596
597 // PileUpJetID variables
598
599 Int_t NCharged;
600 Int_t NNeutrals;
601 Float_t Beta;
602 Float_t BetaStar;
603 Float_t MeanSqDeltaR;
604 Float_t PTD;
605 Float_t FracPt[5];
606
607 // Timing information
608
609 Int_t NTimeHits;
610 std::vector< std::pair< Float_t, Float_t > > ECalEnergyTimePairs;
611
612 // Isolation variables
613
614 Float_t IsolationVar;
615 Float_t IsolationVarRhoCorr;
616 Float_t SumPtCharged;
617 Float_t SumPtNeutral;
618 Float_t SumPtChargedPU;
619 Float_t SumPt;
620
621 // vertex variables
622
623 Int_t ClusterIndex;
624 Int_t ClusterNDF;
625 Double_t ClusterSigma;
626 Double_t SumPT2;
627 Double_t BTVSumPT2;
628 Double_t GenDeltaZ;
629 Double_t GenSumPT2;
630
631 // N-subjettiness variables
632
633 Float_t Tau[5];
634
635 // Other Substructure variables
636
637 TLorentzVector SoftDroppedJet;
638 TLorentzVector SoftDroppedSubJet1;
639 TLorentzVector SoftDroppedSubJet2;
640
641 TLorentzVector TrimmedP4[5]; // first entry (i = 0) is the total Trimmed Jet 4-momenta and from i = 1 to 4 are the trimmed subjets 4-momenta
642 TLorentzVector PrunedP4[5]; // first entry (i = 0) is the total Pruned Jet 4-momenta and from i = 1 to 4 are the pruned subjets 4-momenta
643 TLorentzVector SoftDroppedP4[5]; // first entry (i = 0) is the total SoftDropped Jet 4-momenta and from i = 1 to 4 are the pruned subjets 4-momenta
644
645 Int_t NSubJetsTrimmed; // number of subjets trimmed
646 Int_t NSubJetsPruned; // number of subjets pruned
647 Int_t NSubJetsSoftDropped; // number of subjets soft-dropped
648
649 // Exclusive clustering variables
650 Double_t ExclYmerge23;
651 Double_t ExclYmerge34;
652 Double_t ExclYmerge45;
653 Double_t ExclYmerge56;
654
655 static CompBase *fgCompare; //!
656 const CompBase *GetCompare() const { return fgCompare; }
657
658 void AddCandidate(Candidate *object);
659 TObjArray *GetCandidates();
660
661 Bool_t Overlaps(const Candidate *object) const;
662
663 virtual void Copy(TObject &object) const;
664 virtual TObject *Clone(const char *newname = "") const;
665 virtual void Clear(Option_t* option = "");
666
667private:
668 DelphesFactory *fFactory; //!
669 TObjArray *fArray; //!
670
671 void SetFactory(DelphesFactory *factory) { fFactory = factory; }
672
673 ClassDef(Candidate, 5)
674};
675
676#endif // DelphesClasses_h
677
678
Note: See TracBrowser for help on using the repository browser.