![]() |
![]() |
00001 #ifndef BASICGENPARICLEVALIDATION_H 00002 #define BASICGENPARICLEVALIDATION_H 00003 00004 /*class BasicGenParticleValidation 00005 * 00006 * Class to fill Event Generator dqm monitor elements; works on HepMCProduct 00007 * 00008 * $Date: 2010/07/02 13:34:22 $ 00009 * $Revision: 1.2 $ 00010 * 00011 */ 00012 00013 // framework & common header files 00014 #include "FWCore/Framework/interface/EDAnalyzer.h" 00015 #include "FWCore/Framework/interface/Event.h" 00016 #include "FWCore/Framework/interface/EventSetup.h" 00017 #include "FWCore/Framework/interface/Run.h" 00018 00019 #include "DataFormats/Common/interface/Handle.h" 00020 #include "FWCore/Framework/interface/ESHandle.h" 00021 #include "FWCore/ParameterSet/interface/ParameterSet.h" 00022 #include "FWCore/Utilities/interface/InputTag.h" 00023 00024 #include "FWCore/MessageLogger/interface/MessageLogger.h" 00025 00026 //DQM services 00027 #include "DQMServices/Core/interface/DQMStore.h" 00028 #include "FWCore/ServiceRegistry/interface/Service.h" 00029 #include "DQMServices/Core/interface/MonitorElement.h" 00030 00031 #include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h" 00032 #include "DataFormats/HepMCCandidate/interface/GenParticle.h" 00033 #include "DataFormats/JetReco/interface/GenJetCollection.h" 00034 00035 #include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h" 00036 00037 00038 class BasicGenParticleValidation : public edm::EDAnalyzer 00039 { 00040 public: 00041 explicit BasicGenParticleValidation(const edm::ParameterSet&); 00042 virtual ~BasicGenParticleValidation(); 00043 virtual void beginJob(); 00044 virtual void endJob(); 00045 virtual void analyze(const edm::Event&, const edm::EventSetup&); 00046 virtual void beginRun(const edm::Run&, const edm::EventSetup&); 00047 virtual void endRun(const edm::Run&, const edm::EventSetup&); 00048 00049 bool matchParticles(const HepMC::GenParticle*&, const reco::GenParticle*&); 00050 00051 private: 00052 00053 edm::InputTag hepmcCollection_; 00054 edm::InputTag genparticleCollection_; 00055 edm::InputTag genjetCollection_; 00056 double matchPr_; 00057 00058 unsigned int verbosity_; 00059 00061 edm::ESHandle<HepPDT::ParticleDataTable> fPDGTable ; 00062 00064 DQMStore *dbe; 00065 00066 MonitorElement *nEvt; 00067 00068 // Basic reco::GenParticle test 00069 00070 MonitorElement *genPMultiplicity; 00071 MonitorElement *genMatched; 00072 MonitorElement *multipleMatching; 00073 MonitorElement *matchedResolution; 00074 00075 // Basci GenJets analysis 00076 00077 MonitorElement *genJetMult; 00078 MonitorElement *genJetEnergy; 00079 MonitorElement *genJetPt; 00080 MonitorElement *genJetEta; 00081 MonitorElement *genJetPhi; 00082 MonitorElement *genJetDeltaEtaMin; 00083 00084 MonitorElement *genJetPto1; 00085 MonitorElement *genJetPto10; 00086 MonitorElement *genJetPto100; 00087 MonitorElement *genJetCentral; 00088 00089 MonitorElement *genJetTotPt; 00090 00091 }; 00092 00093 #endif