ECALRegFEDSelector.cc

Go to the documentation of this file.

 
#include "Calibration/HcalIsolatedTrackReco/interface/ECALRegFEDSelector.h"
#include "EventFilter/EcalRawToDigi/interface/EcalRegionCabling.h"
#include "DataFormats/HcalIsolatedTrack/interface/IsolatedPixelTrackCandidateFwd.h"
#include "DataFormats/HcalIsolatedTrack/interface/IsolatedPixelTrackCandidate.h"
#include "DataFormats/Math/interface/RectangularEtaPhiRegion.h"
 
ECALRegFEDSelector::ECALRegFEDSelector(const edm::ParameterSet& iConfig) {
  tok_seed_ = consumes<trigger::TriggerFilterObjectWithRefs>(iConfig.getParameter<edm::InputTag>("regSeedLabel"));
  delta_ = iConfig.getParameter<double>("delta");
 
  tok_raw_ = consumes<FEDRawDataCollection>(iConfig.getParameter<edm::InputTag>("rawInputLabel"));
 
  ec_mapping = std::make_unique<EcalElectronicsMapping>();
 
  produces<FEDRawDataCollection>();
  produces<EcalListOfFEDS>();
 
  for (int p = 0; p < 1200; p++) {
    fedSaved[p] = false;
  }
}
 
ECALRegFEDSelector::~ECALRegFEDSelector() {}
 
void ECALRegFEDSelector::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
  for (int p = 0; p < 1200; p++) {
    fedSaved[p] = false;
  }
 
  auto producedData = std::make_unique<FEDRawDataCollection>();
 
  auto fedList = std::make_unique<EcalListOfFEDS>();
 
  edm::Handle<trigger::TriggerFilterObjectWithRefs> trigSeedTrks;
  iEvent.getByToken(tok_seed_, trigSeedTrks);
 
  std::vector<edm::Ref<reco::IsolatedPixelTrackCandidateCollection> > isoPixTrackRefs;
  trigSeedTrks->getObjects(trigger::TriggerTrack, isoPixTrackRefs);
 
  edm::Handle<FEDRawDataCollection> rawIn;
  iEvent.getByToken(tok_raw_, rawIn);
 
  //  std::vector<int> EC_FED_IDs;
 
  for (uint32_t p = 0; p < isoPixTrackRefs.size(); p++) {
    double etaObj_ = isoPixTrackRefs[p]->track()->eta();
    double phiObj_ = isoPixTrackRefs[p]->track()->phi();
 
    RectangularEtaPhiRegion ecEtaPhi(etaObj_ - delta_, etaObj_ + delta_, phiObj_ - delta_, phiObj_ + delta_);
 
    const std::vector<int> EC_FED_IDs = ec_mapping->GetListofFEDs(ecEtaPhi);
 
    const FEDRawDataCollection* rdc = rawIn.product();
 
    for (int j = 0; j < FEDNumbering::MAXFEDID; j++) {
      bool rightFED = false;
      for (uint32_t k = 0; k < EC_FED_IDs.size(); k++) {
        if (j == EcalRegionCabling::fedIndex(EC_FED_IDs[k])) {
          if (!fedSaved[j]) {
            fedList->AddFED(j);
            rightFED = true;
            fedSaved[j] = true;
          }
        }
      }
      if (j >= FEDNumbering::MINPreShowerFEDID && j <= FEDNumbering::MAXPreShowerFEDID) {
        fedSaved[j] = true;
        rightFED = true;
      }
      if (!rightFED)
        continue;
      const FEDRawData& fedData = rdc->FEDData(j);
      size_t size = fedData.size();
 
      if (size > 0) {
        // this fed has data -- lets copy it
        FEDRawData& fedDataProd = producedData->FEDData(j);
        if (fedDataProd.size() != 0) {
          //          std::cout << " More than one FEDRawDataCollection with data in FED ";
          //          std::cout << j << " Skipping the 2nd\n";
          continue;
        }
        fedDataProd.resize(size);
        unsigned char* dataProd = fedDataProd.data();
        const unsigned char* data = fedData.data();
        for (unsigned int k = 0; k < size; ++k) {
          dataProd[k] = data[k];
        }
      }
    }
  }
 
  iEvent.put(std::move(producedData));
  iEvent.put(std::move(fedList));
}
 
void ECALRegFEDSelector::beginJob() {}
 
void ECALRegFEDSelector::endJob() {}