HiEgammaAlgos/interface/HiPhotonType.h

0001 #ifndef HiPhotonType_h
0002 #define HiPhotonType_h
0003
0004 #include "FWCore/Framework/interface/Event.h"
0005 #include "FWCore/Framework/interface/EventSetup.h"
0006 #include "FWCore/Framework/interface/ESHandle.h"
0007
0008 #include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
0009 #include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
0010 #include "DataFormats/Candidate/interface/CandidateFwd.h"
0011
0012 #include "Geometry/CaloGeometry/interface/CaloGeometry.h"
0013 #include "Geometry/Records/interface/IdealGeometryRecord.h"
0014
0015 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
0016 #include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h"
0017
0018 #include "FWCore/Framework/interface/ESHandle.h"
0019 #include "DataFormats/Common/interface/Handle.h"
0020
0021 #include <vector>
0022
0023 class HiGammaJetSignalDef {
0024 public:
0025   HiGammaJetSignalDef();
0026   HiGammaJetSignalDef(const reco::GenParticleCollection *sigPartic);
0027   bool IsIsolated(const reco::GenParticle &pp);
0028   bool IsIsolatedPP(const reco::GenParticle &pp);
0029   bool IsIsolatedJP(const reco::GenParticle &pp);
0030
0031   //  bool IsSignal(const reco::Candidate &pp, double dPhi, bool isIso);
0032   //  int getIndex(const reco::Candidate &pp);
0033   double getDeltaR(const reco::Candidate &track1, const reco::Candidate &track2);
0034   double getDeltaPhi(const reco::Candidate &track1, const reco::Candidate &track2);
0035   double PI;
0036
0037 private:
0038   const reco::GenParticleCollection *fSigParticles;
0039 };
0040
0041 class HiPhotonType {
0042 public:
0043   HiPhotonType(edm::Handle<reco::GenParticleCollection> inputHandle);
0044   bool IsPrompt(const reco::GenParticle &pp);
0045   bool IsIsolated(const reco::GenParticle &pp);
0046   //  bool IsIsolatedPP(const reco::GenParticle &pp);
0047   //  bool IsIsolatedJP(const reco::GenParticle &pp);
0048   double PI;
0049
0050 private:
0051   HiGammaJetSignalDef mcisocut;
0052 };
0053
0054 #endif