7 |
|
#include <iostream> |
8 |
|
|
9 |
|
// user include files |
10 |
– |
#include "FWCore/Framework/interface/Frameworkfwd.h" |
11 |
– |
#include "FWCore/Framework/interface/EDAnalyzer.h" |
10 |
|
#include "FWCore/Framework/interface/Event.h" |
13 |
– |
//#include "FWCore/Framework/interface/Handle.h" |
14 |
– |
#include "FWCore/Framework/interface/MakerMacros.h" |
11 |
|
#include "FWCore/ParameterSet/interface/ParameterSet.h" |
12 |
|
#include "FWCore/Framework/interface/ESHandle.h" |
13 |
|
|
18 |
– |
#include "SimDataFormats/HepMCProduct/interface/HepMCProduct.h" |
19 |
– |
|
14 |
|
#include "DataFormats/JetReco/interface/GenJet.h" |
15 |
< |
#include "DataFormats/JetReco/interface/GenJetCollection.h" |
15 |
> |
#include "DataFormats/METReco/interface/GenMET.h" |
16 |
|
#include "DataFormats/HepMCCandidate/interface/GenParticle.h" |
17 |
|
#include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h" |
18 |
|
|
19 |
+ |
#include "SimDataFormats/HepMCProduct/interface/HepMCProduct.h" |
20 |
+ |
|
21 |
|
#include "SimDataFormats/Track/interface/SimTrack.h" |
22 |
|
#include "SimDataFormats/Track/interface/SimTrackContainer.h" |
23 |
|
#include "SimDataFormats/Vertex/interface/SimVertex.h" |
32 |
|
#include "../interface/TRootSignalEvent.h" |
33 |
|
#include "../interface/TRootTopTop.h" |
34 |
|
#include "../interface/TRootParticle.h" |
35 |
+ |
#include "../interface/TRootMCParticle.h" |
36 |
|
#include "../interface/TRootMCPhoton.h" |
37 |
|
#include "../interface/TRootJet.h" |
38 |
|
|
39 |
|
#include "TClonesArray.h" |
40 |
|
|
41 |
|
|
42 |
+ |
class MCAnalyzer |
43 |
+ |
{ |
44 |
+ |
|
45 |
+ |
public: |
46 |
+ |
MCAnalyzer(const edm::ParameterSet& config, const edm::ParameterSet& producersNames); |
47 |
+ |
~MCAnalyzer(); |
48 |
+ |
void setVerbosity(int verbosity) {verbosity_ = verbosity; }; |
49 |
+ |
void drawMCTree(const edm::Event& iEvent, const edm::EventSetup& iSetup, const edm::ParameterSet& config, const edm::ParameterSet& producersNames); |
50 |
+ |
bool pdfInfo(const edm::Event& iEvent, TRootEvent* rootEvent); |
51 |
+ |
bool processMCParticle(const edm::Event& iEvent, TClonesArray* rootMCParticles); |
52 |
+ |
bool processGenJets(const edm::Event& iEvent, TClonesArray* rootGenJets); |
53 |
+ |
bool processGenMETs(const edm::Event& iEvent, TClonesArray* rootGenMETs); |
54 |
+ |
bool processConvertedPhoton(const edm::Event& iEvent, TClonesArray* rootMCPhotons); |
55 |
+ |
bool processMuMuGammaEvent(const edm::Event& iEvent, TRootSignalEvent* rootSignalEvent); |
56 |
+ |
bool processTopTopEvent(const edm::Event& iEvent, TClonesArray* rootMCTopTop); |
57 |
+ |
|
58 |
+ |
private: |
59 |
+ |
|
60 |
+ |
int verbosity_; |
61 |
+ |
|
62 |
+ |
bool doPhotonMC_; |
63 |
+ |
double photonMC_etaMax_; |
64 |
+ |
double photonMC_ptMin_; |
65 |
+ |
bool doElectronMC_; |
66 |
+ |
double electronMC_etaMax_; |
67 |
+ |
double electronMC_ptMin_; |
68 |
+ |
bool doMuonMC_; |
69 |
+ |
double muonMC_etaMax_; |
70 |
+ |
double muonMC_ptMin_; |
71 |
+ |
bool doJetMC_; |
72 |
+ |
double jetMC_etaMax_; |
73 |
+ |
double jetMC_ptMin_; |
74 |
+ |
bool doMETMC_; |
75 |
+ |
bool doUnstablePartsMC_; |
76 |
+ |
|
77 |
+ |
std::string signalGenerator_; |
78 |
+ |
edm::InputTag genParticlesProducer_; |
79 |
+ |
edm::InputTag genJetsProducer_; |
80 |
+ |
edm::InputTag genMETsProducer_; |
81 |
+ |
|
82 |
+ |
bool allowMissingCollection_; |
83 |
|
|
46 |
– |
class MCAnalyzer{ |
47 |
– |
|
48 |
– |
public: |
49 |
– |
MCAnalyzer(); |
50 |
– |
MCAnalyzer(const edm::ParameterSet& config, const edm::ParameterSet& producersNames); |
51 |
– |
~MCAnalyzer(); |
52 |
– |
void SetVerbosity(int verbosity) {verbosity_ = verbosity; }; |
53 |
– |
void DrawMCTree(const edm::Event& iEvent, const edm::EventSetup& iSetup, const edm::ParameterSet& config, const edm::ParameterSet& producersNames); |
54 |
– |
void PDFInfo(const edm::Event& iEvent, TRootEvent* rootEvent); |
55 |
– |
void ProcessMCParticle(const edm::Event& iEvent, TClonesArray* rootMCParticles); |
56 |
– |
//void ProcessMCJets(const edm::Event& iEvent, TClonesArray* rootMCJets); |
57 |
– |
void ProcessConvertedPhoton(const edm::Event& iEvent, TClonesArray* rootMCPhotons); |
58 |
– |
void ProcessMuMuGammaEvent(const edm::Event& iEvent, TRootSignalEvent* rootSignalEvent); |
59 |
– |
void ProcessTopTopEvent(const edm::Event& iEvent, TClonesArray* rootMCTopTop); |
60 |
– |
|
61 |
– |
private: |
62 |
– |
|
63 |
– |
int verbosity_; |
64 |
– |
|
65 |
– |
bool doPhotonMC_; |
66 |
– |
double photonMC_etaMax_; |
67 |
– |
double photonMC_ptMin_; |
68 |
– |
bool doElectronMC_; |
69 |
– |
double electronMC_etaMax_; |
70 |
– |
double electronMC_ptMin_; |
71 |
– |
bool doMuonMC_; |
72 |
– |
double muonMC_etaMax_; |
73 |
– |
double muonMC_ptMin_; |
74 |
– |
bool doUnstablePartsMC_; |
75 |
– |
|
76 |
– |
std::string signalGenerator_; |
77 |
– |
edm::InputTag genParticlesProducer_; |
78 |
– |
edm::InputTag genJetsProducer_; |
84 |
|
}; |
85 |
|
|
86 |
|
#endif |