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