EgammaPhotonTkIsolationProducer.cc

Go to the documentation of this file.

//*****************************************************************************
// File:      EgammaPhotonTkIsolationProducer.cc
// ----------------------------------------------------------------------------
// OrigAuth:  Matthias Mozer
// Institute: IIHE-VUB
//=============================================================================
//*****************************************************************************


#include "RecoEgamma/EgammaIsolationAlgos/plugins/EgammaPhotonTkIsolationProducer.h"

// Framework
#include "DataFormats/Common/interface/Handle.h"

#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Candidate/interface/CandAssociation.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/RecoCandidate/interface/RecoCandidate.h"
#include "DataFormats/EgammaReco/interface/SuperCluster.h"
#include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"

#include "RecoEgamma/EgammaIsolationAlgos/interface/PhotonTkIsolation.h"


EgammaPhotonTkIsolationProducer::EgammaPhotonTkIsolationProducer(const edm::ParameterSet& config) :

    // use configuration file to setup input/output collection names
    photonProducer_       (config.getParameter<edm::InputTag>("photonProducer")),
  
    trackProducer_        (config.getParameter<edm::InputTag>("trackProducer")),
    beamspotProducer_     (config.getParameter<edm::InputTag>("BeamspotProducer")),

    ptMin_                (config.getParameter<double>("ptMin")),
    intRadiusBarrel_      (config.getParameter<double>("intRadiusBarrel")),
    intRadiusEndcap_      (config.getParameter<double>("intRadiusEndcap")),
    stripBarrel_          (config.getParameter<double>("stripBarrel")),
    stripEndcap_          (config.getParameter<double>("stripEndcap")),
    extRadius_            (config.getParameter<double>("extRadius")),
    maxVtxDist_           (config.getParameter<double>("maxVtxDist")),
    drb_                  (config.getParameter<double>("maxVtxDistXY"))
{

  //register your products
  produces < edm::ValueMap<double> >();
}


EgammaPhotonTkIsolationProducer::~EgammaPhotonTkIsolationProducer(){}


//
// member functions
//

// ------------ method called to produce the data  ------------
void
EgammaPhotonTkIsolationProducer::produce(edm::StreamID sid, edm::Event& iEvent, const edm::EventSetup& iSetup) const
{
  
  // Get the  filtered objects
  edm::Handle< edm::View<reco::Candidate> > photonHandle;
  iEvent.getByLabel(photonProducer_,photonHandle);
  
  //get the tracks
  edm::Handle<reco::TrackCollection> tracks;
  iEvent.getByLabel(trackProducer_,tracks);
  const reco::TrackCollection* trackCollection = tracks.product();

  edm::Handle<reco::BeamSpot> beamSpotH;
  iEvent.getByLabel(beamspotProducer_,beamSpotH);
  reco::TrackBase::Point beamspot = beamSpotH->position();

  //prepare product
  auto isoMap = std::make_unique<edm::ValueMap<double>>();
  edm::ValueMap<double>::Filler filler(*isoMap);
  std::vector<double> retV(photonHandle->size(),0);

  PhotonTkIsolation myTkIsolation (extRadius_,intRadiusBarrel_,intRadiusEndcap_,stripBarrel_,stripEndcap_,ptMin_,maxVtxDist_,drb_,trackCollection,beamspot) ;

  for(unsigned int i = 0 ; i < photonHandle->size(); ++i ){
    double isoValue = myTkIsolation.getIso(&(photonHandle->at(i))).second;
    retV[i] = isoValue;
  }
   
  //fill and insert valuemap
  filler.insert(photonHandle,retV.begin(),retV.end());
  filler.fill();
  iEvent.put(std::move(isoMap));

}