DQMTauProducer.cc

Go to the documentation of this file.

#include "RecoTauTag/HLTProducers/interface/DQMTauProducer.h"
#include "DataFormats/TrackReco/interface/Track.h"
//
// class decleration
//


DQMTauProducer::DQMTauProducer(const edm::ParameterSet& iConfig)
{
  trackIsolatedJets_ = consumes<reco::IsolatedTauTagInfoCollection>(iConfig.getParameter<edm::InputTag>("TrackIsoJets"));
  matchingCone_ = iConfig.getParameter<double>("MatchingCone");
  signalCone_ = iConfig.getParameter<double>("SignalCone");
  ptMin_ = iConfig.getParameter<double>("MinPtTracks");
  
  isolationCone_ = iConfig.getParameter<double>("IsolationCone");
  produces<reco::HLTTauCollection>();
}

DQMTauProducer::~DQMTauProducer(){ }

void DQMTauProducer::produce(edm::Event& iEvent, const edm::EventSetup& iES)
{

  using namespace reco;
  using namespace edm;
  using namespace std;
  

  HLTTauCollection * jetCollection = new HLTTauCollection;
 

  edm::Handle<IsolatedTauTagInfoCollection> tauL25Jets;
  iEvent.getByToken(trackIsolatedJets_, tauL25Jets );
  

  IsolatedTauTagInfoCollection tau = *(tauL25Jets.product());

  
  float eta_, phi_, pt_;
  float ptLeadTk=0.;
  int trackIsolation=1000.;  
  int nTracks=1000.;

 for(unsigned int i=0; i<tau.size();i++)
   {
     JetTracksAssociationRef jetTracks = tau[i].jtaRef();
     math::XYZVector jetDir(jetTracks->first->px(),jetTracks->first->py(),jetTracks->first->pz());   
     eta_ = jetDir.eta();
     phi_ = jetDir.phi();
     pt_ = jetTracks->first->pt();
     
     const TrackRef leadTk = tau[i].leadingSignalTrack(jetDir,matchingCone_, 1.);
     if(!leadTk) 
       {}else{
     trackIsolation = (int)tau[i].discriminator(jetDir,matchingCone_, signalCone_, isolationCone_,1.,1.,0);
     ptLeadTk = (*leadTk).pt();
     nTracks = (tau[i].tracksInCone((*leadTk).momentum(), isolationCone_,   ptMin_ )).size() - (tau[i].tracksInCone((*leadTk).momentum(), signalCone_, ptMin_ )).size();
       }
     HLTTau pippo(eta_,phi_,pt_,-1,trackIsolation,ptLeadTk,trackIsolation,ptLeadTk);
     pippo.setNL25TrackIsolation(nTracks);
     pippo.setNL3TrackIsolation(nTracks);
     jetCollection->push_back(pippo);
     
   }
 
   
  std::auto_ptr<reco::HLTTauCollection> selectedTaus(jetCollection);
  
  iEvent.put(selectedTaus);
  



}