CaloMETProducer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    METProducers
// Class:      CaloMETProducer
//
//
 
//____________________________________________________________________________||
#include "RecoMET/METProducers/interface/CaloMETProducer.h"
 
#include "FWCore/Framework/interface/ConsumesCollector.h"
 
#include "FWCore/Utilities/interface/InputTag.h"
 
#include "DataFormats/Common/interface/Handle.h"
 
#include "DataFormats/METReco/interface/METFwd.h"
#include "DataFormats/METReco/interface/CaloMETFwd.h"
#include "DataFormats/METReco/interface/CaloMET.h"
#include "DataFormats/METReco/interface/CommonMETData.h"
 
#include "RecoMET/METAlgorithms/interface/METAlgo.h"
#include "RecoMET/METAlgorithms/interface/SignAlgoResolutions.h"
#include "RecoMET/METAlgorithms/interface/CaloSpecificAlgo.h"
#include "RecoMET/METAlgorithms/interface/SignCaloSpecificAlgo.h"
 
#include <cstring>
 
//____________________________________________________________________________||
namespace cms {
 
  //____________________________________________________________________________||
  CaloMETProducer::CaloMETProducer(const edm::ParameterSet& iConfig)
      : inputToken_(consumes<edm::View<reco::Candidate> >(iConfig.getParameter<edm::InputTag>("src"))),
        calculateSignificance_(iConfig.getParameter<bool>("calculateSignificance")),
        resolutions_(nullptr),
        globalThreshold_(iConfig.getParameter<double>("globalThreshold")) {
    noHF_ = iConfig.getParameter<bool>("noHF");
 
    std::string alias = iConfig.exists("alias") ? iConfig.getParameter<std::string>("alias") : "";
 
    produces<reco::CaloMETCollection>().setBranchAlias(alias);
 
    if (calculateSignificance_)
      resolutions_ = new metsig::SignAlgoResolutions(iConfig);
  }
 
  //____________________________________________________________________________||
  CaloMETProducer::~CaloMETProducer() {
    if (resolutions_)
      delete resolutions_;
  }
 
  //____________________________________________________________________________||
  void CaloMETProducer::produce(edm::Event& event, const edm::EventSetup& setup) {
    edm::Handle<edm::View<reco::Candidate> > input;
    event.getByToken(inputToken_, input);
 
    METAlgo algo;
    CommonMETData commonMETdata = algo.run(*input.product(), globalThreshold_);
 
    CaloSpecificAlgo calospecalgo;
    reco::CaloMET calomet = calospecalgo.addInfo(input, commonMETdata, noHF_, globalThreshold_);
 
    /*
    if( calculateSignificance_ )
      {
    SignCaloSpecificAlgo signcalospecalgo;
    signcalospecalgo.calculateBaseCaloMET(input, commonMETdata, *resolutions_, noHF_, globalThreshold_);
    calomet.SetMetSignificance(signcalospecalgo.getSignificance());
    calomet.setSignificanceMatrix(signcalospecalgo.getSignificanceMatrix());
      }
*/
    auto calometcoll = std::make_unique<reco::CaloMETCollection>();
    calometcoll->push_back(calomet);
    event.put(std::move(calometcoll));
  }
 
  //____________________________________________________________________________||
  DEFINE_FWK_MODULE(CaloMETProducer);
}  // namespace cms
 
//____________________________________________________________________________||