L1TMuonLegacyConverter.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    L1TMuonLegacyConverter
// Class:      L1TMuonLegacyConverter
// 
//
// Original Author:  Bortignon Pierluigi
//         Created:  Sun March 6 EDT 2016
//
//


// system include files
#include <memory>

// user include files
#include "L1Trigger/L1TCommon/plugins/L1TMuonLegacyConverter.h"

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/Common/interface/OrphanHandle.h"

#include "CondFormats/L1TObjects/interface/L1CaloGeometry.h"
#include "CondFormats/DataRecord/interface/L1CaloGeometryRecord.h"

#include "CondFormats/L1TObjects/interface/L1CaloEtScale.h"
#include "CondFormats/DataRecord/interface/L1EmEtScaleRcd.h"
#include "CondFormats/DataRecord/interface/L1JetEtScaleRcd.h"
#include "CondFormats/DataRecord/interface/L1HtMissScaleRcd.h"
#include "CondFormats/DataRecord/interface/L1HfRingEtScaleRcd.h"
#include "CondFormats/L1TObjects/interface/L1GctJetFinderParams.h"
#include "CondFormats/DataRecord/interface/L1GctJetFinderParamsRcd.h"

#include "CondFormats/L1TObjects/interface/L1MuTriggerScales.h"
#include "CondFormats/DataRecord/interface/L1MuTriggerScalesRcd.h"
#include "CondFormats/L1TObjects/interface/L1MuTriggerPtScale.h"
#include "CondFormats/DataRecord/interface/L1MuTriggerPtScaleRcd.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"


// #include "FWCore/Utilities/interface/EDMException.h"

//
// class decleration
//


//
// constants, enums and typedefs
//
using namespace l1t;


//
// static data member definitions
//

double const L1TMuonLegacyConverter::muonMassGeV_ = 0.105658369 ; // PDG06

//
// constructors and destructor
//
L1TMuonLegacyConverter::L1TMuonLegacyConverter(const edm::ParameterSet& iConfig)
{
   using namespace l1extra ;

   // moving inputTag here
   muonSource_InputTag = iConfig.getParameter<edm::InputTag>("muonSource");

   produces<MuonBxCollection>("imdMuonsLegacy");
   muonSource_InputToken = consumes<L1MuGMTReadoutCollection>(muonSource_InputTag);
}


L1TMuonLegacyConverter::~L1TMuonLegacyConverter()
{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}


//
// member functions
//

// ------------ method called to produce the data  ------------
void
L1TMuonLegacyConverter::produce( edm::Event& iEvent,
                   const edm::EventSetup& iSetup)
{
   using namespace edm ;
   using namespace l1extra ;
   using namespace std ;

   // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   // ~~~~~~~~~~~~~~~~~~~~ Muons ~~~~~~~~~~~~~~~~~~~~
   // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

   auto_ptr< MuonBxCollection > imdMuonsLegacy( new MuonBxCollection() );

   if( produceMuonParticles_ )
   {
      ESHandle< L1MuTriggerScales > muScales ;
      iSetup.get< L1MuTriggerScalesRcd >().get( muScales ) ;

      ESHandle< L1MuTriggerPtScale > muPtScale ;
      iSetup.get< L1MuTriggerPtScaleRcd >().get( muPtScale ) ;

      Handle< L1MuGMTReadoutCollection > hwMuCollection ;
      iEvent.getByToken( muonSource_InputToken, hwMuCollection ) ;

      vector< L1MuGMTExtendedCand > hwMuCands ;

      if( !hwMuCollection.isValid() )
    {
      LogDebug("L1TMuonLegacyConverter")
        << "\nWarning: L1MuGMTReadoutCollection with " << muonSource_InputTag
        << "\nrequested in configuration, but not found in the event."
        << std::endl;
    }
      else
    {
      if( centralBxOnly_ )
        {
          // Get GMT candidates from central bunch crossing only
          hwMuCands = hwMuCollection->getRecord().getGMTCands() ;
        }
      else
        {
          // Get GMT candidates from all bunch crossings
          vector< L1MuGMTReadoutRecord > records = hwMuCollection->getRecords();
          vector< L1MuGMTReadoutRecord >::const_iterator rItr = records.begin();
          vector< L1MuGMTReadoutRecord >::const_iterator rEnd = records.end();

          for( ; rItr != rEnd ; ++rItr )
        {
          vector< L1MuGMTExtendedCand > tmpCands = rItr->getGMTCands() ;

          hwMuCands.insert( hwMuCands.end(),
                    tmpCands.begin(),
                    tmpCands.end() ) ;
        }
        }

//       cout << "HW muons" << endl ;
      vector< L1MuGMTExtendedCand >::const_iterator muItr = hwMuCands.begin() ;
      vector< L1MuGMTExtendedCand >::const_iterator muEnd = hwMuCands.end() ;
      for( int i = 0 ; muItr != muEnd ; ++muItr, ++i )
        {
//   cout << "#" << i
//        << " name " << muItr->name()
//        << " empty " << muItr->empty()
//        << " pt " << muItr->ptIndex()
//        << " eta " << muItr->etaIndex()
//        << " phi " << muItr->phiIndex()
//        << " iso " << muItr->isol()
//        << " mip " << muItr->mip()
//        << " bx " << muItr->bx()
//        << endl ;

          if( !muItr->empty() ){
          // keep x and y components non-zero and protect against roundoff.
          double pt = muPtScale->getPtScale()->getLowEdge( muItr->ptIndex() ) + 1.e-6 ;
          double eta = muScales->getGMTEtaScale()->getCenter( muItr->etaIndex() ) ;
          double phi = muScales->getPhiScale()->getLowEdge( muItr->phiIndex() ) ;

          math::PtEtaPhiMLorentzVector p4( pt,
                           eta,
                           phi,
                           muonMassGeV_ ) ;

        // from L1TMuonProducer.cc - which is the stage2 muon producer:
        // Muon outMu{vec, mu->hwPt(), mu->hwEta(), mu->hwGlobalPhi(), outMuQual, mu->hwSign(), mu->hwSignValid(), iso, mu->tfMuonIndex(), 0, true, mu->hwIsoSum(), mu->hwDPhi(), mu->hwDEta(), mu->hwRank()};

          Muon outMu{p4, (int)muItr->ptIndex(), (int)muItr->etaIndex(), (int)muItr->phiIndex(), (int)muItr->quality(), muItr->charge(), muItr->isol(), (int)muItr->etaRegionIndex(), 0, true, 0, 0, 0 , (int)muItr->rank() };
          imdMuonsLegacy->push_back( muItr->bx(), outMu ) ;
        }
        }
    }
   }
   
   iEvent.put( imdMuonsLegacy, "imdMuonsLegacy" );

} // closing produce

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