/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/JetMETCorrections/Type1MET/src/Type1MET.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    Type1MET
// Class:      Type1MET
// 
//
// Original Author:  Oct 12 08:23
//         Created:  Wed Oct 12 12:16:04 CDT 2005
// $Id: Type1MET.cc,v 1.21 2010/09/15 16:01:21 lacroix Exp $
//
//


// user include files
#include "JetMETCorrections/Type1MET/interface/Type1MET.h"

#include "DataFormats/Common/interface/View.h"
#include "DataFormats/METReco/interface/MET.h"
#include "DataFormats/METReco/interface/METCollection.h"
#include "DataFormats/METReco/interface/PFMETCollection.h"
#include "DataFormats/METReco/interface/CaloMET.h"
#include "DataFormats/METReco/interface/CaloMETCollection.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "JetMETCorrections/Objects/interface/JetCorrector.h"
#include "DataFormats/JetReco/interface/PFJetCollection.h"
#include "DataFormats/PatCandidates/interface/Jet.h"

#include "DataFormats/MuonReco/interface/MuonMETCorrectionData.h"


using namespace reco;

namespace cms 
{
  // PRODUCER CONSTRUCTORS ------------------------------------------
  Type1MET::Type1MET( const edm::ParameterSet& iConfig ) : alg_() 
  {
    metType             = iConfig.getParameter<std::string>("metType");
    
    inputUncorMetLabel  = iConfig.getParameter<std::string>("inputUncorMetLabel");
    inputUncorJetsLabel = iConfig.getParameter<std::string>("inputUncorJetsLabel");
    correctorLabel   = iConfig.getParameter<std::string>("corrector");
    jetPTthreshold      = iConfig.getParameter<double>("jetPTthreshold");
    jetEMfracLimit      = iConfig.getParameter<double>("jetEMfracLimit");
    UscaleA      = iConfig.getParameter<double>("UscaleA");
    UscaleB      = iConfig.getParameter<double>("UscaleB");
    UscaleC      = iConfig.getParameter<double>("UscaleC");
    useTypeII      = iConfig.getParameter<bool>("useTypeII");
    hasMuonsCorr   = iConfig.getParameter<bool>("hasMuonsCorr");
    if( metType == "CaloMET" )
      produces<CaloMETCollection>();
    else
      produces<METCollection>();
  }
  Type1MET::Type1MET() : alg_() {}
  // PRODUCER DESTRUCTORS -------------------------------------------
  Type1MET::~Type1MET() {}

  // PRODUCER METHODS -----------------------------------------------
  void Type1MET::produce( edm::Event& iEvent, const edm::EventSetup& iSetup )
  {
    using namespace edm;

// Remove lable from the code
       Handle<View<reco::Muon> > inputMuons;
       iEvent.getByLabel( "muons", inputMuons );
// Remove lable from the code
       Handle<ValueMap<reco::MuonMETCorrectionData> > vm_muCorrData_h;
       iEvent.getByLabel( "muonMETValueMapProducer","muCorrData", vm_muCorrData_h);

    if( metType == "CaloMET")
      {
    Handle<CaloJetCollection> inputUncorJets;
    iEvent.getByLabel( inputUncorJetsLabel, inputUncorJets );
    const JetCorrector* corrector = JetCorrector::getJetCorrector (correctorLabel, iSetup);

        Handle<CaloMETCollection> inputUncorMet;                     //Define Inputs
        iEvent.getByLabel( inputUncorMetLabel,  inputUncorMet );     //Get Inputs
        std::auto_ptr<CaloMETCollection> output( new CaloMETCollection() );  //Create empty output
        alg_.run( *(inputUncorMet.product()), *corrector, *(inputUncorJets.product()), 
                  jetPTthreshold, jetEMfracLimit, UscaleA, UscaleB, UscaleC, useTypeII, hasMuonsCorr,
                  *(inputMuons.product()), *(vm_muCorrData_h.product()),
                  &*output );                                         //Invoke the algorithm
        iEvent.put( output );                                        //Put output into Event
      }
    else if (metType == "PFMET")
      {
      Handle<PFJetCollection> inputUncorJets;
      iEvent.getByLabel( inputUncorJetsLabel, inputUncorJets );
      const JetCorrector* corrector = JetCorrector::getJetCorrector (correctorLabel, iSetup);

        Handle<PFMETCollection> inputUncorMet;                     //Define Inputs
        iEvent.getByLabel( inputUncorMetLabel,  inputUncorMet );     //Get Inputs
        std::auto_ptr<METCollection> output( new METCollection() );  //Create empty output
        alg_.run( *(inputUncorMet.product()), *corrector, *(inputUncorJets.product()), 
                  jetPTthreshold, jetEMfracLimit, UscaleA, UscaleB, UscaleC, useTypeII, hasMuonsCorr,
                  *(inputMuons.product()), *(vm_muCorrData_h.product()),
                  &*output );                                         //Invoke the algorithm
        iEvent.put( output );                                        //Put output into Event
      }
    else
      {
        Handle<pat::JetCollection> inputUncorJets;
      iEvent.getByLabel( inputUncorJetsLabel, inputUncorJets );
      const JetCorrector* corrector = JetCorrector::getJetCorrector (correctorLabel, iSetup);

        Handle<PFMETCollection> inputUncorMet;                     //Define Inputs
        iEvent.getByLabel( inputUncorMetLabel,  inputUncorMet );     //Get Inputs
        std::auto_ptr<METCollection> output( new METCollection() );  //Create empty output
        alg_.run( *(inputUncorMet.product()), *corrector, *(inputUncorJets.product()), 
                  jetPTthreshold, jetEMfracLimit, UscaleA, UscaleB, UscaleC, useTypeII, hasMuonsCorr,
                  *(inputMuons.product()), *(vm_muCorrData_h.product()),
                  &*output );                                         //Invoke the algorithm
        iEvent.put( output );                                        //Put output into Event
      }
  }

  //  DEFINE_FWK_MODULE(Type1MET);  //define this as a plug-in

}//end namespace cms