1 |
#ifndef UltraFastSim_H
|
2 |
#define UltraFastSim_H
|
3 |
|
4 |
#include <vector>
|
5 |
|
6 |
namespace Pythia8 {
|
7 |
class Rndm;
|
8 |
class Event;
|
9 |
class Particle;
|
10 |
}
|
11 |
|
12 |
namespace fastjet {
|
13 |
class JetDefinition;
|
14 |
class ClusterSequence;
|
15 |
}
|
16 |
|
17 |
#include "fastjet/PseudoJet.hh";
|
18 |
|
19 |
class UltraFastSim {
|
20 |
|
21 |
public:
|
22 |
|
23 |
UltraFastSim(Pythia8::Rndm *r);
|
24 |
|
25 |
virtual ~UltraFastSim() {;}
|
26 |
|
27 |
bool run(Pythia8::Event &event);
|
28 |
|
29 |
private:
|
30 |
|
31 |
UltraFastSim();
|
32 |
|
33 |
void clear();
|
34 |
|
35 |
void makeJets();
|
36 |
void makeBJets();
|
37 |
void makeTaus();
|
38 |
|
39 |
void tkSmear(Pythia8::Particle particle);
|
40 |
void emSmear(Pythia8::Particle particle);
|
41 |
void hdSmear(Pythia8::Particle particle);
|
42 |
|
43 |
Pythia8::Rndm* rndmPtr;
|
44 |
|
45 |
fastjet::JetDefinition* jetDefPtr;
|
46 |
fastjet::ClusterSequence* cs;
|
47 |
|
48 |
float trackerResolution;
|
49 |
float ecalResolution;
|
50 |
float ecalConstantTerm;
|
51 |
float hcalResolution;
|
52 |
float hcalConstantTerm;
|
53 |
|
54 |
std::vector<fastjet::PseudoJet> selectedParticles;
|
55 |
std::vector<fastjet::PseudoJet> selectedElectrons;
|
56 |
std::vector<fastjet::PseudoJet> selectedMuons;
|
57 |
std::vector<fastjet::PseudoJet> selectedTaus;
|
58 |
std::vector<fastjet::PseudoJet> selectedBQuarks;
|
59 |
std::vector<fastjet::PseudoJet> jets;
|
60 |
std::vector<fastjet::PseudoJet> sortedJets;
|
61 |
std::vector<fastjet::PseudoJet> bJets;
|
62 |
std::vector<fastjet::PseudoJet> taus;
|
63 |
|
64 |
};
|
65 |
|
66 |
#endif
|