RPCRecHitProducer.cc

Go to the documentation of this file.

 
#include "RPCRecHitProducer.h"
 
#include "FWCore/Framework/interface/ESHandle.h"
 
#include "Geometry/RPCGeometry/interface/RPCRoll.h"
#include "Geometry/RPCGeometry/interface/RPCGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "DataFormats/MuonDetId/interface/RPCDetId.h"
#include "DataFormats/RPCRecHit/interface/RPCRecHit.h"
 
#include "RecoLocalMuon/RPCRecHit/interface/RPCRecHitAlgoFactory.h"
#include "DataFormats/RPCRecHit/interface/RPCRecHitCollection.h"
 
#include "CondFormats/DataRecord/interface/RPCMaskedStripsRcd.h"
#include "CondFormats/DataRecord/interface/RPCDeadStripsRcd.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include <string>
#include <fstream>
 
using namespace edm;
using namespace std;
 
RPCRecHitProducer::RPCRecHitProducer(const ParameterSet& config)
    : theRPCDigiLabel(consumes<RPCDigiCollection>(config.getParameter<InputTag>("rpcDigiLabel"))),
      // Get the concrete reconstruction algo from the factory
      theAlgo{RPCRecHitAlgoFactory::get()->create(config.getParameter<string>("recAlgo"),
                                                  config.getParameter<ParameterSet>("recAlgoConfig"))},
      maskSource_(MaskSource::EventSetup),
      deadSource_(MaskSource::EventSetup) {
  // Set verbose output
  produces<RPCRecHitCollection>();
 
  // Get masked- and dead-strip information
  theRPCMaskedStripsObj = std::make_unique<RPCMaskedStrips>();
  theRPCDeadStripsObj = std::make_unique<RPCDeadStrips>();
 
  const string maskSource = config.getParameter<std::string>("maskSource");
  if (maskSource == "File") {
    maskSource_ = MaskSource::File;
    edm::FileInPath fp = config.getParameter<edm::FileInPath>("maskvecfile");
    std::ifstream inputFile(fp.fullPath().c_str(), std::ios::in);
    if (!inputFile) {
      std::cerr << "Masked Strips File cannot not be opened" << std::endl;
      exit(1);
    }
    while (inputFile.good()) {
      RPCMaskedStrips::MaskItem Item;
      inputFile >> Item.rawId >> Item.strip;
      if (inputFile.good())
        MaskVec.push_back(Item);
    }
    inputFile.close();
  }
 
  const string deadSource = config.getParameter<std::string>("deadSource");
  if (deadSource == "File") {
    deadSource_ = MaskSource::File;
    edm::FileInPath fp = config.getParameter<edm::FileInPath>("deadvecfile");
    std::ifstream inputFile(fp.fullPath().c_str(), std::ios::in);
    if (!inputFile) {
      std::cerr << "Dead Strips File cannot not be opened" << std::endl;
      exit(1);
    }
    while (inputFile.good()) {
      RPCDeadStrips::DeadItem Item;
      inputFile >> Item.rawId >> Item.strip;
      if (inputFile.good())
        DeadVec.push_back(Item);
    }
    inputFile.close();
  }
}
 
void RPCRecHitProducer::beginRun(const edm::Run& r, const edm::EventSetup& setup) {
  // Getting the masked-strip information
  if (maskSource_ == MaskSource::EventSetup) {
    edm::ESHandle<RPCMaskedStrips> readoutMaskedStrips;
    setup.get<RPCMaskedStripsRcd>().get(readoutMaskedStrips);
    const RPCMaskedStrips* tmp_obj = readoutMaskedStrips.product();
    theRPCMaskedStripsObj->MaskVec = tmp_obj->MaskVec;
    delete tmp_obj;
  } else if (maskSource_ == MaskSource::File) {
    std::vector<RPCMaskedStrips::MaskItem>::iterator posVec;
    for (posVec = MaskVec.begin(); posVec != MaskVec.end(); ++posVec) {
      RPCMaskedStrips::MaskItem Item;
      Item.rawId = (*posVec).rawId;
      Item.strip = (*posVec).strip;
      theRPCMaskedStripsObj->MaskVec.push_back(Item);
    }
  }
 
  // Getting the dead-strip information
  if (deadSource_ == MaskSource::EventSetup) {
    edm::ESHandle<RPCDeadStrips> readoutDeadStrips;
    setup.get<RPCDeadStripsRcd>().get(readoutDeadStrips);
    const RPCDeadStrips* tmp_obj = readoutDeadStrips.product();
    theRPCDeadStripsObj->DeadVec = tmp_obj->DeadVec;
    delete tmp_obj;
  } else if (deadSource_ == MaskSource::File) {
    std::vector<RPCDeadStrips::DeadItem>::iterator posVec;
    for (posVec = DeadVec.begin(); posVec != DeadVec.end(); ++posVec) {
      RPCDeadStrips::DeadItem Item;
      Item.rawId = (*posVec).rawId;
      Item.strip = (*posVec).strip;
      theRPCDeadStripsObj->DeadVec.push_back(Item);
    }
  }
}
 
void RPCRecHitProducer::produce(Event& event, const EventSetup& setup) {
  // Get the RPC Geometry
  ESHandle<RPCGeometry> rpcGeom;
  setup.get<MuonGeometryRecord>().get(rpcGeom);
 
  // Get the digis from the event
  Handle<RPCDigiCollection> digis;
  event.getByToken(theRPCDigiLabel, digis);
 
  // Pass the EventSetup to the algo
  theAlgo->setES(setup);
 
  // Create the pointer to the collection which will store the rechits
  auto recHitCollection = std::make_unique<RPCRecHitCollection>();
 
  // Iterate through all digi collections ordered by LayerId
 
  for (auto rpcdgIt = digis->begin(); rpcdgIt != digis->end(); ++rpcdgIt) {
    // The layerId
    const RPCDetId& rpcId = (*rpcdgIt).first;
 
    // Get the GeomDet from the setup
    const RPCRoll* roll = rpcGeom->roll(rpcId);
    if (roll == nullptr) {
      edm::LogError("BadDigiInput") << "Failed to find RPCRoll for ID " << rpcId;
      continue;
    }
 
    // Get the iterators over the digis associated with this LayerId
    const RPCDigiCollection::Range& range = (*rpcdgIt).second;
 
    // Getting the roll mask, that includes dead strips, for the given RPCDet
    RollMask mask;
    const int rawId = rpcId.rawId();
    for (const auto& tomask : theRPCMaskedStripsObj->MaskVec) {
      if (tomask.rawId == rawId) {
        const int bit = tomask.strip;
        mask.set(bit - 1);
      }
    }
 
    for (const auto& tomask : theRPCDeadStripsObj->DeadVec) {
      if (tomask.rawId == rawId) {
        const int bit = tomask.strip;
        mask.set(bit - 1);
      }
    }
 
    // Call the reconstruction algorithm
    OwnVector<RPCRecHit> recHits = theAlgo->reconstruct(*roll, rpcId, range, mask);
 
    if (!recHits.empty())  //FIXME: is it really needed?
      recHitCollection->put(rpcId, recHits.begin(), recHits.end());
  }
 
  event.put(std::move(recHitCollection));
}