L1TBMTFAlgoSelector.cc

Go to the documentation of this file.

/*
/P.Katsoulis
/G.Karathanasis     
*/

#include "L1TBMTFAlgoSelector.h"

//The Constructor defines what to be consumed and produced
dqmBmtfAlgoSelector::L1TBMTFAlgoSelector::L1TBMTFAlgoSelector(const edm::ParameterSet & ps)
{
  bmtfKalmanToken = consumes<l1t::RegionalMuonCandBxCollection>(ps.getParameter<edm::InputTag>("bmtfKalman"));
  bmtfLegacyToken = consumes<l1t::RegionalMuonCandBxCollection>(ps.getParameter<edm::InputTag>("bmtfLegacy"));
  fedToken = consumes<FEDRawDataCollection>(ps.getParameter<edm::InputTag>("feds"));

  produces<l1t::RegionalMuonCandBxCollection>("BMTF");
  produces<l1t::RegionalMuonCandBxCollection>("BMTF2");
}


dqmBmtfAlgoSelector::L1TBMTFAlgoSelector::~L1TBMTFAlgoSelector() {
}

void dqmBmtfAlgoSelector::L1TBMTFAlgoSelector::produce(edm::Event & eve, const edm::EventSetup & eveSetup)
{

  //Declare the Handlers
  edm::Handle<FEDRawDataCollection> feds;
  eve.getByToken(fedToken, feds);
  edm::Handle<l1t::RegionalMuonCandBxCollection> bmtfKalman;
  eve.getByToken(bmtfKalmanToken, bmtfKalman);
  edm::Handle<l1t::RegionalMuonCandBxCollection> bmtfLegacy;
  eve.getByToken(bmtfLegacyToken, bmtfLegacy);


  //----> Make RegMuonCand Copies
  auto *bmtfKalman_copy = new l1t::RegionalMuonCandBxCollection(*bmtfKalman);
  auto *bmtfLegacy_copy = new l1t::RegionalMuonCandBxCollection(*bmtfLegacy);
  edm::LogInfo("L1TDQM") << "copy RegionalMuonCandBxCollections created";
  edm::LogInfo("L1TDQM") << "bmtfKalman_copy address: " << bmtfKalman_copy;
  edm::LogInfo("L1TDQM") << "bmtfLegacy_copy address: " << bmtfLegacy_copy;


  //-----> Get the fw-ver
  int nonEmptyFed=0;
  if (feds->FEDData(1376).size() > 0)
    nonEmptyFed = 1376;
  else if (feds->FEDData(1377).size() > 0)
    nonEmptyFed = 1377;
  else {
    edm::LogError("L1TDQM") << "[L1TBMTFAlgoSelector] Both BMTF feds (1376, 1377) seem empty."
                << "Using Default Configuration."
                << "Will be wrong if this is Kalman triggering and ZS is enabled";
    //Define the default configuration
    std::unique_ptr<l1t::RegionalMuonCandBxCollection> bmtfTriggering, bmtfSecondary;
    bmtfTriggering.reset(bmtfLegacy_copy);
    bmtfSecondary.reset(bmtfKalman_copy);
    //Print in log the addresses
    edm::LogInfo("L1TDQM") << "Triggering and Secondary pointers filled:";
    edm::LogInfo("L1TDQM") << "bmtfTriggering address: " << bmtfTriggering.get();
    edm::LogInfo("L1TDQM") << "bmtfSecondary address: " << bmtfSecondary.get();
    //Produce the products
    eve.put(std::move(bmtfTriggering),"BMTF");
    eve.put(std::move(bmtfSecondary),"BMTF2");
    return;
  }

  const FEDRawData& l1tRcd = feds->FEDData(nonEmptyFed);;
  edm::LogInfo("L1TDQM") << "L1T Rcd taken from the FEDData.";
  edm::LogInfo("L1TDQM") << "l1tRcd.size=" << l1tRcd.size() << "   for fed:" << nonEmptyFed;

  const unsigned char *data = l1tRcd.data();
  FEDHeader header(data);
  edm::LogInfo("L1TDQM") << "header and data extracted from the Rcd.";

  amc13::Packet packet;
  if (!packet.parse( (const uint64_t*) data, (const uint64_t*) (data + 8),
                     (l1tRcd.size()) / 8,
                     header.lvl1ID(), header.bxID(), false, false ) ) {

    edm::LogError("L1TDQM") << "Could not extract AMC13 Packet.";
    return;
  }

  edm::LogInfo("L1TDQM") << "AMC13-packet-payload size = " << packet.payload().size();
  unsigned algo_ver;
  if (!packet.payload().empty()) {
    auto payload64 = ( packet.payload().at(0) ).data();
    const uint32_t *start = (const uint32_t*) payload64.get();
    const uint32_t *end = start + (packet.payload().at(0).size() * 2);

    l1t::MP7Payload payload(start, end, false);
    algo_ver = payload.getAlgorithmFWVersion();

    edm::LogInfo("L1TDQM") << "algo-ver = " << algo_ver << std::endl;
  }
  else {
    edm::LogError("L1TDQM") << "amc13 payload is empty, cannot extract AMC13 Packet...";
    return;
  }


  //----->Make the Decision which Algo Triggers
  std::unique_ptr<l1t::RegionalMuonCandBxCollection> bmtfTriggering, bmtfSecondary;
  if ( algo_ver >= 2499805536) {//95000160(hex)
    // kalman triggers
    bmtfTriggering.reset(bmtfKalman_copy);
    bmtfSecondary.reset(bmtfLegacy_copy);
  }
  else {
    // legacy triggers
    bmtfTriggering.reset(bmtfLegacy_copy);
    bmtfSecondary.reset(bmtfKalman_copy);
  }

  edm::LogInfo("L1TDQM") << "Triggering and Secondary pointers filled:";
  edm::LogInfo("L1TDQM") << "bmtfTriggering address: " << bmtfTriggering.get();
  edm::LogInfo("L1TDQM") << "bmtfSecondary address: " << bmtfSecondary.get();
  
  eve.put(std::move(bmtfTriggering),"BMTF");
  eve.put(std::move(bmtfSecondary),"BMTF2");

  return;
}

using namespace dqmBmtfAlgoSelector;
DEFINE_FWK_MODULE(L1TBMTFAlgoSelector);