EmulateCPPF.cc

Go to the documentation of this file.

#include "CondFormats/RPCObjects/interface/RPCDeadStrips.h"
#include "CondFormats/RPCObjects/interface/RPCMaskedStrips.h"
#include "L1Trigger/L1TMuonCPPF/interface/EmulateCPPF.h"
#include <fstream>
#include <string>

EmulateCPPF::EmulateCPPF(const edm::ParameterSet &iConfig,
                         edm::ConsumesCollector &&iConsumes)
    :  // rpcDigi_processors_(),
      recHit_processors_(),
      // rpcDigiToken_(
      // iConsumes.consumes<RPCTag::digi_collection>(iConfig.getParameter<edm::InputTag>("recHitLabel"))
      // ),
      recHitToken_(iConsumes.consumes<RPCRecHitCollection>(iConfig.getParameter<edm::InputTag>("recHitLabel"))),
      cppfSource_(CppfSource::EventSetup),
      MaxClusterSize_(0) {
  MaxClusterSize_ = iConfig.getParameter<int>("MaxClusterSize");

  const std::string cppfSource = iConfig.getParameter<std::string>("cppfSource");
  // Look up table
  if (cppfSource == "File") {
    cppfSource_ = CppfSource::File;
    edm::FileInPath fp = iConfig.getParameter<edm::FileInPath>("cppfvecfile");
    std::ifstream inputFile(fp.fullPath().c_str(), std::ios::in);
    if (!inputFile) {
      throw cms::Exception("No LUT") << "Error: CPPF look up table file cannot not be opened";
      exit(1);
    }
    while (inputFile.good()) {
      RecHitProcessor::CppfItem Item;
      inputFile >> Item.rawId >> Item.strip >> Item.lb >> Item.halfchannel >> Item.int_phi >> Item.int_theta;
      if (inputFile.good())
        CppfVec_1.push_back(Item);
    }
    inputFile.close();
  }

  // RPC Geometry
  else if (cppfSource == "Geo") {
    cppfSource_ = CppfSource::EventSetup;
  }
  // Error for wrong input
  else {
    throw cms::Exception("Invalid option") << "Error: Specify in python/emulatorCppfDigis_cfi 'File' for look up  "
                                              "table or 'Geo' for RPC Geometry";
    exit(1);
  }
}

EmulateCPPF::~EmulateCPPF() {}

void EmulateCPPF::process(const edm::Event &iEvent,
                          const edm::EventSetup &iSetup,
                          // l1t::CPPFDigiCollection& cppf_rpcDigi,
                          l1t::CPPFDigiCollection &cppf_recHit) {
  if (cppfSource_ == CppfSource::File) {
    // Using the look up table to fill the information
    cppf_recHit.clear();
    for (auto &recHit_processor : recHit_processors_) {
      recHit_processor.processLook(iEvent, iSetup, recHitToken_, CppfVec_1, cppf_recHit, MaxClusterSize_);
      //  recHit_processors_.at(recHit_processor).processLook( iEvent, iSetup,
      //  recHitToken_, CppfVec_1, cppf_recHit );
    }
  } else if (cppfSource_ == CppfSource::EventSetup) {
    // Clear output collections
    // cppf_rpcDigi.clear();
    cppf_recHit.clear();

    // // Get the RPCDigis from the event
    // edm::Handle<RPCTag::digi_collection> rpcDigis;
    // iEvent.getByToken(rpcDigiToken_, rpcDigis);

    // _________________________________________________________________________________
    // Run the CPPF clusterization+coordinate conversion algo on RPCDigis and
    // RecHits

    // For now, treat CPPF as single board
    // In the future, may want to treat the 4 CPPF boards in each endcap as
    // separate entities

    // for (unsigned int iBoard = 0; iBoard < rpcDigi_processors_.size();
    // iBoard++) { rpcDigi_processors_.at(iBoard).process( iSetup, rpcDigis,
    // cppf_rpcDigi );
    // }
    for (auto &recHit_processor : recHit_processors_) {
      recHit_processor.process(iEvent, iSetup, recHitToken_, cppf_recHit);
      // recHit_processors_.at(recHit_processor).process( iEvent, iSetup,
      // recHitToken_, cppf_recHit );
    }
  }
}  // End void EmulateCPPF::process()