IsolationProducerForTracks.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/ValueMap.h"

class IsolationProducerForTracks : public edm::EDProducer {
public:
  IsolationProducerForTracks(const edm::ParameterSet & );
private:
  void produce(edm::Event& event, const edm::EventSetup& setup);

  edm::InputTag tracks_;
  edm::InputTag highPtTracks_;
  edm::InputTag isoDeps_;
  double trackPtMin_;
  double coneSize_;

};


#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Common/interface/Association.h"
#include "DataFormats/Common/interface/Ref.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "DataFormats/Candidate/interface/OverlapChecker.h"
#include "DataFormats/RecoCandidate/interface/IsoDeposit.h"
#include "DataFormats/RecoCandidate/interface/IsoDepositFwd.h"
#include <iostream>
#include <iterator>
#include <vector>

using namespace edm;
using namespace reco;

typedef edm::ValueMap<float> TkIsoMap;

IsolationProducerForTracks::IsolationProducerForTracks(const ParameterSet & pset) : 
  tracks_( pset.getParameter<InputTag>( "tracks" ) ),
  highPtTracks_( pset.getParameter<InputTag>( "highPtTracks" ) ),
  isoDeps_( pset.getParameter<InputTag>( "isoDeps" ) ),
  trackPtMin_( pset.getParameter<double>( "trackPtMin" ) ),
  coneSize_( pset.getParameter<double>( "coneSize" ) )
 {
  produces<TkIsoMap>();
 }

void IsolationProducerForTracks::produce(Event & event, const EventSetup & setup) {
  std::auto_ptr<TkIsoMap> caloIsolations(new TkIsoMap);
  TkIsoMap::Filler filler(*caloIsolations);
  {
    Handle<CandidateView> tracks;
    event.getByLabel(tracks_, tracks);
    
    Handle<CandidateView> highPtTracks;
    event.getByLabel(highPtTracks_, highPtTracks);
    
    Handle<IsoDepositMap> isoDeps;
    event.getByLabel(isoDeps_, isoDeps);
    
    int nTracks = tracks->size();
    int nHighPtTracks = highPtTracks->size();
    std::vector<double> iso(nTracks);   

    OverlapChecker overlap;
    
    for(int i = 0; i < nTracks; ++i ) {
      const Candidate & tkCand = (*tracks)[ i ];
      double caloIso = - 1.0;
      if( tkCand.pt() > trackPtMin_) {
    for(int j = 0; j < nHighPtTracks; ++j ) {
      const Candidate & highPtTkCand = (*highPtTracks)[ j ];
      if(overlap(tkCand, highPtTkCand) ) {
        CandidateBaseRef tkRef = highPtTracks->refAt(j);
        const IsoDeposit &isoDep = (*isoDeps)[tkRef];
        caloIso = isoDep.depositWithin(coneSize_);
        break;
      }
    }
      } 
      iso[i] = caloIso;
    }
    filler.insert(tracks, iso.begin(), iso.end());
  }

  // really fill the association map
  filler.fill();
  event.put(caloIsolations);
}  
  
#include "FWCore/Framework/interface/MakerMacros.h"

DEFINE_FWK_MODULE(IsolationProducerForTracks);