DTCollector.cc

Go to the documentation of this file.

#include <memory>

#include "DataFormats/L1DTTrackFinder/interface/L1MuDTChambPhContainer.h"
#include "DataFormats/L1DTTrackFinder/interface/L1MuDTChambThContainer.h"
#include "DataFormats/MuonDetId/interface/DTChamberId.h"
#include "L1Trigger/L1TMuon/interface/deprecate/DTCollector.h"

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

using namespace L1TMuon;

DTCollector::DTCollector(const edm::ParameterSet& ps)
    : SubsystemCollector(ps), bx_min(ps.getParameter<int>("BX_min")), bx_max(ps.getParameter<int>("BX_max")) {
  if (ps.getParameter<bool>("runBunchCrossingCleaner")) {
    const edm::ParameterSet& bxccfg = ps.getParameterSet("bxCleanerCfg");
    _bxc = std::make_unique<DTBunchCrossingCleaner>(bxccfg);
  } else {
    _bxc.reset(nullptr);
  }
}

void DTCollector::extractPrimitives(const edm::Event& ev,
                                    const edm::EventSetup& es,
                                    TriggerPrimitiveCollection& out) const {
  TriggerPrimitiveCollection cleaned, temp, chamb_list;
  edm::Handle<L1MuDTChambPhContainer> phiDigis;
  edm::Handle<L1MuDTChambThContainer> thetaDigis;
  ev.getByLabel(_src, phiDigis);
  ev.getByLabel(_src, thetaDigis);
  for (int wheel = -2; wheel <= 2; ++wheel) {
    for (int station = 1; station <= 4; ++station) {
      for (int sector = 0; sector <= 11; ++sector) {
        chamb_list.clear();
        for (int bx = bx_min; bx <= bx_max; ++bx) {
          std::unique_ptr<const L1MuDTChambPhDigi> phi_segm_1(phiDigis->chPhiSegm1(wheel, station, sector, bx));
          std::unique_ptr<const L1MuDTChambPhDigi> phi_segm_2(phiDigis->chPhiSegm2(wheel, station, sector, bx));
          std::unique_ptr<const L1MuDTChambThDigi> theta_segm(thetaDigis->chThetaSegm(wheel, station, sector, bx));

          int bti_group_1 = -1, bti_group_2 = -1;

          if (theta_segm) {
            bti_group_1 = findBTIGroupForThetaDigi(*theta_segm, 1);
            bti_group_2 = findBTIGroupForThetaDigi(*theta_segm, 2);
          }

          if (phi_segm_1 && bti_group_1 != -1) {
            chamb_list.push_back(processDigis(*phi_segm_1, *theta_segm, bti_group_1));
          } else if (phi_segm_1 && bti_group_1 == -1) {
            chamb_list.push_back(processDigis(*phi_segm_1, 1));
          } else if (!phi_segm_1 && bti_group_1 != -1) {
            chamb_list.push_back(processDigis(*theta_segm, bti_group_1));
          }

          if (phi_segm_2 && bti_group_2 != -1) {
            chamb_list.push_back(processDigis(*phi_segm_2, *theta_segm, bti_group_2));
          } else if (phi_segm_2 && bti_group_2 == -1) {
            chamb_list.push_back(processDigis(*phi_segm_2, 2));
          } else if (!phi_segm_2 && bti_group_2 != -1) {
            chamb_list.push_back(processDigis(*phi_segm_2, bti_group_2));
          }

          phi_segm_1.release();
          phi_segm_2.release();
          theta_segm.release();
        }
        if (_bxc) {
          temp = _bxc->clean(chamb_list);
          cleaned.insert(cleaned.end(), temp.begin(), temp.end());
        } else {
          cleaned.insert(cleaned.end(), chamb_list.begin(), chamb_list.end());
        }
      }
    }
  }
  out.insert(out.end(), cleaned.begin(), cleaned.end());
}

TriggerPrimitive DTCollector::processDigis(const L1MuDTChambPhDigi& digi, const int& segment_number) const {
  DTChamberId detid(digi.whNum(), digi.stNum(), digi.scNum() + 1);
  return TriggerPrimitive(detid, digi, segment_number);
}

TriggerPrimitive DTCollector::processDigis(const L1MuDTChambThDigi& digi_th, const int bti_group) const {
  DTChamberId detid(digi_th.whNum(), digi_th.stNum(), digi_th.scNum() + 1);
  return TriggerPrimitive(detid, digi_th, bti_group);
}

TriggerPrimitive DTCollector::processDigis(const L1MuDTChambPhDigi& digi_phi,
                                           const L1MuDTChambThDigi& digi_theta,
                                           const int bti_group) const {
  DTChamberId detid(digi_phi.whNum(), digi_phi.stNum(), digi_phi.scNum() + 1);
  return TriggerPrimitive(detid, digi_phi, digi_theta, bti_group);
}

int DTCollector::findBTIGroupForThetaDigi(const L1MuDTChambThDigi& digi, const int pos) const {
  //if( digi.stNum() == 4 ) return -1; // there is no theta layer there
  int result = -1;
  for (int i = 0; i < 7; ++i) {
    if (digi.position(i) == pos)
      result = i;
  }
  return result;
}

#include "L1Trigger/L1TMuon/interface/deprecate/SubsystemCollectorFactory.h"
DEFINE_EDM_PLUGIN(SubsystemCollectorFactory, DTCollector, "DTCollector");