15 |
|
} |
16 |
|
|
17 |
|
#include "fastjet/PseudoJet.hh"; |
18 |
+ |
|
19 |
|
#include "TParticle.h"; |
20 |
+ |
#include "TLorentzVector.h"; |
21 |
|
|
22 |
|
class UltraFastSim { |
23 |
|
|
29 |
|
|
30 |
|
bool run(Pythia8::Event &event); |
31 |
|
|
32 |
< |
const std::vector<TParticle>& electronList() {return selectedElectrons;} |
33 |
< |
const std::vector<TParticle>& muonList() {return selectedMuons;} |
34 |
< |
const std::vector<TParticle>& genTauList() {return selectedTaus;} |
35 |
< |
const std::vector<TParticle>& bQuarkList() {return selectedBQuarks;} |
32 |
> |
const std::vector<TParticle>& genTauList() {return genTaus;} |
33 |
> |
const std::vector<TParticle>& visTauList() {return visTaus;} |
34 |
> |
const std::vector<TParticle>& bQuarkList() {return bQuarks;} |
35 |
> |
const std::vector<TParticle>& cQuarkList() {return cQuarks;} |
36 |
> |
const std::vector<TParticle>& photonList() {return photons;} |
37 |
> |
const std::vector<TParticle>& electronList() {return electrons;} |
38 |
> |
const std::vector<TParticle>& muonList() {return muons;} |
39 |
|
const std::vector<TParticle>& tauList() {return taus;} |
40 |
< |
const std::vector<fastjet::PseudoJet>& particleList() {return selectedParticles;} |
40 |
> |
const std::vector<TParticle>& chargedHadronList() {return chargedHadrons;} |
41 |
> |
const std::vector<TParticle>& neutralHadronList() {return neutralHadrons;} |
42 |
|
const std::vector<fastjet::PseudoJet>& jetList() {return jets;} |
37 |
– |
const std::vector<fastjet::PseudoJet>& sortedJetList() {return sortedJets;} |
43 |
|
const std::vector<fastjet::PseudoJet>& bJetList() {return bJets;} |
44 |
+ |
const std::vector<fastjet::PseudoJet>& bJetListStdGeom() {return bJetsStdGeom;} |
45 |
+ |
|
46 |
+ |
const TLorentzVector& getMET() {return MET;} |
47 |
+ |
const TLorentzVector& getMHT() {return MHT;} |
48 |
+ |
|
49 |
+ |
const double getET() {return ET;} |
50 |
+ |
const double getHT() {return HT;} |
51 |
|
|
52 |
|
private: |
53 |
|
|
58 |
|
void makeJets(); |
59 |
|
void makeBJets(); |
60 |
|
void makeTaus(); |
61 |
+ |
void makeETSums(); |
62 |
|
|
63 |
< |
void tkSmear(Pythia8::Particle particle); |
64 |
< |
void emSmear(Pythia8::Particle particle); |
65 |
< |
void hdSmear(Pythia8::Particle particle); |
63 |
> |
void setCommon(Pythia8::Particle& particle, TParticle& smearedParticle); |
64 |
> |
void tkSmear(Pythia8::Particle& particle, TParticle& smearedParticle); |
65 |
> |
void emSmear(Pythia8::Particle& particle, TParticle& smearedParticle); |
66 |
> |
void hdSmear(Pythia8::Particle& particle, TParticle& smearedParticle); |
67 |
|
|
68 |
|
Pythia8::Rndm* rndmPtr; |
69 |
|
|
70 |
|
fastjet::JetDefinition* jetDefPtr; |
71 |
|
fastjet::ClusterSequence* cs; |
72 |
|
|
73 |
< |
float trackerResolution; |
74 |
< |
float ecalResolution; |
75 |
< |
float ecalConstantTerm; |
76 |
< |
float hcalResolution; |
77 |
< |
float hcalConstantTerm; |
78 |
< |
|
79 |
< |
std::vector<TParticle> selectedElectrons; |
66 |
< |
std::vector<TParticle> selectedMuons; |
67 |
< |
std::vector<TParticle> selectedTaus; |
68 |
< |
std::vector<TParticle> selectedBQuarks; |
73 |
> |
std::vector<TParticle> genTaus; |
74 |
> |
std::vector<TParticle> visTaus; |
75 |
> |
std::vector<TParticle> bQuarks; |
76 |
> |
std::vector<TParticle> cQuarks; |
77 |
> |
std::vector<TParticle> photons; |
78 |
> |
std::vector<TParticle> electrons; |
79 |
> |
std::vector<TParticle> muons; |
80 |
|
std::vector<TParticle> taus; |
81 |
< |
std::vector<fastjet::PseudoJet> selectedParticles; |
81 |
> |
std::vector<TParticle> chargedHadrons; |
82 |
> |
std::vector<TParticle> neutralHadrons; |
83 |
|
std::vector<fastjet::PseudoJet> jets; |
72 |
– |
std::vector<fastjet::PseudoJet> sortedJets; |
84 |
|
std::vector<fastjet::PseudoJet> bJets; |
85 |
+ |
std::vector<fastjet::PseudoJet> bJetsStdGeom; |
86 |
+ |
|
87 |
+ |
TLorentzVector MET; |
88 |
+ |
TLorentzVector MHT; |
89 |
+ |
|
90 |
+ |
double ET; |
91 |
+ |
double HT; |
92 |
|
|
93 |
|
}; |
94 |
|
|