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");
   produceMuonParticles_ = iConfig.getUntrackedParameter<bool>("produceMuonParticles");  
   centralBxOnly_ = iConfig.getUntrackedParameter<bool>("centralBxOnly");

   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 ~~~~~~~~~~~~~~~~~~~~
   // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

   std::unique_ptr< MuonBxCollection > imdMuonsLegacy( new MuonBxCollection() );

   LogDebug("L1TMuonLegacyConverter")
        << "\nWarning: L1MuGMTReadoutCollection with " << muonSource_InputTag
        << "\nrequested in configuration, but not found in the event."
        << std::endl;
    
   if( produceMuonParticles_ )
   {
      ESHandle< L1MuTriggerScales > muScales ;
      iSetup.get< L1MuTriggerScalesRcd >().get( muScales ) ;

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

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

      vector< L1MuGMTExtendedCand > simMuCands ;


      if( !simMuCollection.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
          simMuCands = simMuCollection->getRecord().getGMTCands() ;
        }
      else
        {
          // Get GMT candidates from all bunch crossings
          vector< L1MuGMTReadoutRecord > records = simMuCollection->getRecords();
          vector< L1MuGMTReadoutRecord >::const_iterator rItr = records.begin();
          vector< L1MuGMTReadoutRecord >::const_iterator rEnd = records.end();

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

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

      vector< L1MuGMTExtendedCand >::const_iterator muItr = simMuCands.begin() ;
      vector< L1MuGMTExtendedCand >::const_iterator muEnd = simMuCands.end() ;
      for( int i = 0 ; muItr != muEnd ; ++muItr, ++i )
        {
  
          if( !muItr->empty() ){
            // keep x and y components non-zero and protect against roundoff.
            double pt = muPtScale->getPtScale()->getLowEdge( muItr->ptIndex() ) + 1.e-6 ;
            std::cout << "pt from muPtScale = " << pt << std::endl; 
            double eta = muScales->getGMTEtaScale()->getCenter( muItr->etaIndex() ) ;
            double phi = muScales->getPhiScale()->getLowEdge( muItr->phiIndex() ) ;

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

            Muon outMu{p4, (int)0, (int)0, (int)0, (int)muItr->quality(), (int)muItr->charge(), (int)muItr->charge_valid() , (int)muItr->isol(), (int)0 , (int)0, true, (int)0, (int)0, (int)0 , (int)muItr->rank() };
            imdMuonsLegacy->push_back( muItr->bx(), outMu ) ;
        }
        }
    }
   }
   
   iEvent.put( std::move(imdMuonsLegacy), "imdMuonsLegacy" );

} // closing produce

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