EmulateCPPF.cc

Go to the documentation of this file.

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

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()