L1TMicroGMTLUTDumper.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    L1TMicroGMTLUTDumper
// Class:      L1TMicroGMTLUTDumper
//
//
// Original Author:  Joschka Philip Lingemann,40 3-B01,+41227671598,
//         Created:  Thu Oct  3 10:12:30 CEST 2013
// $Id$
//
//


// system include files
#include <memory>
#include <fstream>

// user include files
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"

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

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "L1Trigger/L1TMuon/interface/MicroGMTRankPtQualLUT.h"
#include "L1Trigger/L1TMuon/interface/MicroGMTAbsoluteIsolationCheckLUT.h"
#include "L1Trigger/L1TMuon/interface/MicroGMTRelativeIsolationCheckLUT.h"
#include "L1Trigger/L1TMuon/interface/MicroGMTCaloIndexSelectionLUT.h"
#include "L1Trigger/L1TMuon/interface/MicroGMTExtrapolationLUT.h"
#include "L1Trigger/L1TMuon/interface/MicroGMTMatchQualLUT.h"
#include "L1Trigger/L1TMuon/interface/MicroGMTLUTFactories.h"

#include "CondFormats/L1TObjects/interface/L1TMuonGlobalParams.h"
#include "CondFormats/DataRecord/interface/L1TMuonGlobalParamsRcd.h"
#include "L1Trigger/L1TMuon/interface/L1TMuonGlobalParamsHelper.h"

#include <iostream>
//
// class declaration
//
using namespace l1t;

class L1TMicroGMTLUTDumper : public edm::EDAnalyzer {
   public:
      explicit L1TMicroGMTLUTDumper(const edm::ParameterSet&);
      ~L1TMicroGMTLUTDumper();
      virtual void analyze(const edm::Event&, const edm::EventSetup&);

   private:
      virtual void beginRun(edm::Run const&, edm::EventSetup const&);

      void dumpLut(MicroGMTLUT*, const std::string&);

      // ----------member data ---------------------------
      std::unique_ptr<L1TMuonGlobalParamsHelper> microGMTParamsHelper;
      std::string m_foldername;

      std::shared_ptr<MicroGMTRankPtQualLUT> m_rankLUT;

      std::shared_ptr<MicroGMTAbsoluteIsolationCheckLUT> m_absIsoCheckMemLUT;
      std::shared_ptr<MicroGMTRelativeIsolationCheckLUT> m_relIsoCheckMemLUT;

      std::shared_ptr<MicroGMTCaloIndexSelectionLUT> m_idxSelMemPhiLUT;
      std::shared_ptr<MicroGMTCaloIndexSelectionLUT> m_idxSelMemEtaLUT;

      std::shared_ptr<MicroGMTExtrapolationLUT> m_bPhiExtrapolationLUT;
      std::shared_ptr<MicroGMTExtrapolationLUT> m_oPhiExtrapolationLUT;
      std::shared_ptr<MicroGMTExtrapolationLUT> m_fPhiExtrapolationLUT;
      std::shared_ptr<MicroGMTExtrapolationLUT> m_bEtaExtrapolationLUT;
      std::shared_ptr<MicroGMTExtrapolationLUT> m_oEtaExtrapolationLUT;
      std::shared_ptr<MicroGMTExtrapolationLUT> m_fEtaExtrapolationLUT;

      std::shared_ptr<MicroGMTMatchQualLUT> m_boPosMatchQualLUT;
      std::shared_ptr<MicroGMTMatchQualLUT> m_boNegMatchQualLUT;
      std::shared_ptr<MicroGMTMatchQualLUT> m_foPosMatchQualLUT;
      std::shared_ptr<MicroGMTMatchQualLUT> m_foNegMatchQualLUT;
      //std::shared_ptr<MicroGMTMatchQualLUT> m_brlSingleMatchQualLUT;
      std::shared_ptr<MicroGMTMatchQualLUT> m_ovlPosSingleMatchQualLUT;
      std::shared_ptr<MicroGMTMatchQualLUT> m_ovlNegSingleMatchQualLUT;
      std::shared_ptr<MicroGMTMatchQualLUT> m_fwdPosSingleMatchQualLUT;
      std::shared_ptr<MicroGMTMatchQualLUT> m_fwdNegSingleMatchQualLUT;
};

//
// constants, enums and typedefs
//


//
// static data member definitions
//

//
// constructors and destructor
//
L1TMicroGMTLUTDumper::L1TMicroGMTLUTDumper(const edm::ParameterSet& iConfig)
{
  //now do what ever other initialization is needed
  m_foldername = iConfig.getParameter<std::string> ("out_directory");

  microGMTParamsHelper = std::unique_ptr<L1TMuonGlobalParamsHelper>(new L1TMuonGlobalParamsHelper());
}


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


//
// member functions
//
void
L1TMicroGMTLUTDumper::dumpLut(MicroGMTLUT* lut, const std::string& oName) {
  std::ofstream fStream(m_foldername+oName);
  lut->save(fStream);
  fStream.close();
}



// ------------ method called to produce the data  ------------
void
L1TMicroGMTLUTDumper::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
{
  using namespace edm;
  dumpLut(m_rankLUT.get(), std::string("/SortRank.txt"));
  dumpLut(m_absIsoCheckMemLUT.get(), std::string("/AbsIsoCheckMem.txt"));
  dumpLut(m_relIsoCheckMemLUT.get(), std::string("/RelIsoCheckMem.txt"));
  dumpLut(m_idxSelMemPhiLUT.get(), std::string("/IdxSelMemPhi.txt"));
  dumpLut(m_idxSelMemEtaLUT.get(), std::string("/IdxSelMemEta.txt"));
  dumpLut(m_bPhiExtrapolationLUT.get(), std::string("/BPhiExtrapolation.txt"));
  dumpLut(m_oPhiExtrapolationLUT.get(), std::string("/OPhiExtrapolation.txt"));
  dumpLut(m_fPhiExtrapolationLUT.get(), std::string("/EPhiExtrapolation.txt"));
  dumpLut(m_bEtaExtrapolationLUT.get(), std::string("/BEtaExtrapolation.txt"));
  dumpLut(m_oEtaExtrapolationLUT.get(), std::string("/OEtaExtrapolation.txt"));
  dumpLut(m_fEtaExtrapolationLUT.get(), std::string("/EEtaExtrapolation.txt"));
  dumpLut(m_boPosMatchQualLUT.get(), std::string("/BOPosMatchQual.txt"));
  dumpLut(m_boNegMatchQualLUT.get(), std::string("/BONegMatchQual.txt"));
  dumpLut(m_foPosMatchQualLUT.get(), std::string("/EOPosMatchQual.txt"));
  dumpLut(m_foNegMatchQualLUT.get(), std::string("/EONegMatchQual.txt"));
  //dumpLut(m_brlSingleMatchQualLUT.get(), std::string("/BmtfSingleMatchQual.txt"));
  dumpLut(m_ovlPosSingleMatchQualLUT.get(), std::string("/OmtfPosSingleMatchQual.txt"));
  dumpLut(m_ovlNegSingleMatchQualLUT.get(), std::string("/OmtfNegSingleMatchQual.txt"));
  dumpLut(m_fwdPosSingleMatchQualLUT.get(), std::string("/EmtfPosSingleMatchQual.txt"));
  dumpLut(m_fwdNegSingleMatchQualLUT.get(), std::string("/EmtfNegSingleMatchQual.txt"));

}

// ------------ method called when starting to processes a run  ------------
void
L1TMicroGMTLUTDumper::beginRun(edm::Run const& run, edm::EventSetup const& iSetup)
{
  const L1TMuonGlobalParamsRcd& microGMTParamsRcd = iSetup.get<L1TMuonGlobalParamsRcd>();
  edm::ESHandle<L1TMuonGlobalParams> microGMTParamsHandle;
  microGMTParamsRcd.get(microGMTParamsHandle);

  microGMTParamsHelper = std::unique_ptr<L1TMuonGlobalParamsHelper>(new L1TMuonGlobalParamsHelper(*microGMTParamsHandle.product()));
  if (!microGMTParamsHelper) {
    edm::LogError("L1TMicroGMTLUTDumper") << "Could not retrieve parameters from Event Setup" << std::endl;
  }

  int fwVersion = microGMTParamsHelper->fwVersion();
  m_rankLUT = MicroGMTRankPtQualLUTFactory::create(microGMTParamsHelper->sortRankLUTPath(), fwVersion, microGMTParamsHelper->sortRankLUTPtFactor(), microGMTParamsHelper->sortRankLUTQualFactor());

  m_absIsoCheckMemLUT = MicroGMTAbsoluteIsolationCheckLUTFactory::create(microGMTParamsHelper->absIsoCheckMemLUTPath(), fwVersion);
  m_relIsoCheckMemLUT = MicroGMTRelativeIsolationCheckLUTFactory::create(microGMTParamsHelper->relIsoCheckMemLUTPath(), fwVersion);
  m_idxSelMemPhiLUT = MicroGMTCaloIndexSelectionLUTFactory::create(microGMTParamsHelper->idxSelMemPhiLUTPath(), l1t::MicroGMTConfiguration::PHI, fwVersion);
  m_idxSelMemEtaLUT = MicroGMTCaloIndexSelectionLUTFactory::create(microGMTParamsHelper->idxSelMemEtaLUTPath(), l1t::MicroGMTConfiguration::ETA, fwVersion);

  m_bPhiExtrapolationLUT = MicroGMTExtrapolationLUTFactory::create(microGMTParamsHelper->bPhiExtrapolationLUTPath(), l1t::MicroGMTConfiguration::PHI_OUT, fwVersion);
  m_oPhiExtrapolationLUT = MicroGMTExtrapolationLUTFactory::create(microGMTParamsHelper->oPhiExtrapolationLUTPath(), l1t::MicroGMTConfiguration::PHI_OUT, fwVersion);
  m_fPhiExtrapolationLUT = MicroGMTExtrapolationLUTFactory::create(microGMTParamsHelper->fPhiExtrapolationLUTPath(), l1t::MicroGMTConfiguration::PHI_OUT, fwVersion);
  m_bEtaExtrapolationLUT = MicroGMTExtrapolationLUTFactory::create(microGMTParamsHelper->bEtaExtrapolationLUTPath(), l1t::MicroGMTConfiguration::ETA_OUT, fwVersion);
  m_oEtaExtrapolationLUT = MicroGMTExtrapolationLUTFactory::create(microGMTParamsHelper->oEtaExtrapolationLUTPath(), l1t::MicroGMTConfiguration::ETA_OUT, fwVersion);
  m_fEtaExtrapolationLUT = MicroGMTExtrapolationLUTFactory::create(microGMTParamsHelper->fEtaExtrapolationLUTPath(), l1t::MicroGMTConfiguration::ETA_OUT, fwVersion);

  m_boPosMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->bOPosMatchQualLUTPath(),
                                                                 microGMTParamsHelper->bOPosMatchQualLUTMaxDR(),
                                                                 microGMTParamsHelper->bOPosMatchQualLUTfEta(),
                                                                 microGMTParamsHelper->bOPosMatchQualLUTfEtaCoarse(),
                                                                 microGMTParamsHelper->bOPosMatchQualLUTfPhi(),
                                                                 cancel_t::omtf_bmtf_pos,
                                                                 fwVersion);
  m_boNegMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->bONegMatchQualLUTPath(),
                                                                 microGMTParamsHelper->bONegMatchQualLUTMaxDR(),
                                                                 microGMTParamsHelper->bONegMatchQualLUTfEta(),
                                                                 microGMTParamsHelper->bONegMatchQualLUTfEtaCoarse(),
                                                                 microGMTParamsHelper->bONegMatchQualLUTfPhi(),
                                                                 cancel_t::omtf_bmtf_neg,
                                                                 fwVersion);
  m_foPosMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->fOPosMatchQualLUTPath(),
                                                                 microGMTParamsHelper->fOPosMatchQualLUTMaxDR(),
                                                                 microGMTParamsHelper->fOPosMatchQualLUTfEta(),
                                                                 microGMTParamsHelper->fOPosMatchQualLUTfEtaCoarse(),
                                                                 microGMTParamsHelper->fOPosMatchQualLUTfPhi(),
                                                                 cancel_t::omtf_emtf_pos,
                                                                 fwVersion);
  m_foNegMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->fONegMatchQualLUTPath(),
                                                                 microGMTParamsHelper->fONegMatchQualLUTMaxDR(),
                                                                 microGMTParamsHelper->fONegMatchQualLUTfEta(),
                                                                 microGMTParamsHelper->fONegMatchQualLUTfEtaCoarse(),
                                                                 microGMTParamsHelper->fONegMatchQualLUTfPhi(),
                                                                 cancel_t::omtf_emtf_neg,
                                                                 fwVersion);
  m_ovlPosSingleMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->ovlPosSingleMatchQualLUTPath(),
                                                                        microGMTParamsHelper->ovlPosSingleMatchQualLUTMaxDR(),
                                                                        microGMTParamsHelper->ovlPosSingleMatchQualLUTfEta(),
                                                                        microGMTParamsHelper->ovlPosSingleMatchQualLUTfEtaCoarse(),
                                                                        microGMTParamsHelper->ovlPosSingleMatchQualLUTfPhi(),
                                                                        cancel_t::omtf_omtf_pos,
                                                                        fwVersion);
  m_ovlNegSingleMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->ovlNegSingleMatchQualLUTPath(),
                                                                        microGMTParamsHelper->ovlNegSingleMatchQualLUTMaxDR(),
                                                                        microGMTParamsHelper->ovlNegSingleMatchQualLUTfEta(),
                                                                        microGMTParamsHelper->ovlNegSingleMatchQualLUTfEtaCoarse(),
                                                                        microGMTParamsHelper->ovlNegSingleMatchQualLUTfPhi(),
                                                                        cancel_t::omtf_omtf_neg,
                                                                        fwVersion);
  m_fwdPosSingleMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->fwdPosSingleMatchQualLUTPath(),
                                                                        microGMTParamsHelper->fwdPosSingleMatchQualLUTMaxDR(),
                                                                        microGMTParamsHelper->fwdPosSingleMatchQualLUTfEta(),
                                                                        microGMTParamsHelper->fwdPosSingleMatchQualLUTfEta(),
                                                                        microGMTParamsHelper->fwdPosSingleMatchQualLUTfPhi(),
                                                                        cancel_t::emtf_emtf_pos,
                                                                        fwVersion);
  m_fwdNegSingleMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->fwdNegSingleMatchQualLUTPath(),
                                                                        microGMTParamsHelper->fwdNegSingleMatchQualLUTMaxDR(),
                                                                        microGMTParamsHelper->fwdNegSingleMatchQualLUTfEta(),
                                                                        microGMTParamsHelper->fwdNegSingleMatchQualLUTfEta(),
                                                                        microGMTParamsHelper->fwdNegSingleMatchQualLUTfPhi(),
                                                                        cancel_t::emtf_emtf_neg,
                                                                        fwVersion);
}

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