L1Trigger/RPCTrigger/plugins/RPCTrigger.cc - master - githubmirror/cmssw

L1Trigger/RPCTrigger/plugins/RPCTrigger.cc
253 строки10 KB
Ben Jones
added esConsumes to modules in L1Trigger/RPCTrigger
30 июн 2020, 18:36
30 июн 2020, 18:36b79db81
#include "L1Trigger/RPCTrigger/interface/RPCTrigger.h"
//#define ML_DEBUG

#include "L1Trigger/RPCTrigger/interface/MuonsGrabber.h"

RPCTrigger::RPCTrigger(const edm::ParameterSet& iConfig)
    : m_trigConfig(),
      m_pacTrigger(),
      // 0 - no debug
      // 1 - dump to xml
      m_triggerDebug{iConfig.getUntrackedParameter<int>("RPCTriggerDebug", 0) == 1 ? 1 : 0},
      m_cacheID{0},
      m_label{iConfig.getParameter<std::string>("label")},
      m_rpcDigiToken{consumes<RPCDigiCollection>(m_label)},

      m_brlCandPutToken{produces<std::vector<L1MuRegionalCand> >("RPCb")},
      m_fwdCandPutToken{produces<std::vector<L1MuRegionalCand> >("RPCf")},

      m_brlLinksPutToken{produces<std::vector<RPCDigiL1Link> >("RPCb")},
      m_fwdLinksPutToken{produces<std::vector<RPCDigiL1Link> >("RPCf")},
      m_configToken(esConsumes<L1RPCConfig, L1RPCConfigRcd>()),
      m_coneBuilderToken(esConsumes<L1RPCConeBuilder, L1RPCConeBuilderRcd>()),
      m_coneDefinitionToken(esConsumes<L1RPCConeDefinition, L1RPCConeDefinitionRcd>()),
      m_hwConfigToken(esConsumes<L1RPCHwConfig, L1RPCHwConfigRcd>()),
      m_BxOrConfigToken(esConsumes<L1RPCBxOrConfig, L1RPCBxOrConfigRcd>()),
      m_hsbConfigToken(esConsumes<L1RPCHsbConfig, L1RPCHsbConfigRcd>()) {
  //MuonsGrabber is a singleton
  usesResource("MuonsGrabber");
}

void RPCTrigger::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
  //static int ev=1;
  //std::cout << "----------------------------------- "  << ev++ << std::endl;
  if (m_triggerDebug == 1)
    MuonsGrabber::Instance().startNewEvent(iEvent.id().event(), iEvent.bunchCrossing());

  if (m_cacheID != iSetup.get<L1RPCConfigRcd>().cacheIdentifier()) {
    //std::cout << " New pats: " << iSetup.get<L1RPCConfigRcd>().cacheIdentifier() << std::endl ;
    m_cacheID = iSetup.get<L1RPCConfigRcd>().cacheIdentifier();

    edm::ESHandle<L1RPCConfig> conf = iSetup.getHandle(m_configToken);
    const L1RPCConfig* rpcconf = conf.product();

    m_pacManager.init(rpcconf);
    m_trigConfig = std::make_unique<RPCBasicTrigConfig>(&m_pacManager);
    m_trigConfig->setDebugLevel(m_triggerDebug);

    m_pacTrigger = std::make_unique<RPCPacTrigger>(m_trigConfig.get());

    if (m_triggerDebug == 1)
      MuonsGrabber::Instance().setRPCBasicTrigConfig(m_trigConfig.get());
  }

  edm::Handle<RPCDigiCollection> rpcDigis;
  iEvent.getByToken(m_rpcDigiToken, rpcDigis);

  std::vector<L1MuRegionalCand> candBarell;
  std::vector<L1MuRegionalCand> candForward;
  if (!rpcDigis.isValid()) {
    LogDebug("RPCTrigger")
        << "\nWarning: RPCDigiCollection with input tag " << m_label
        << "\nrequested in configuration, but not found in the event. Emulator will produce empty collection \n ";

    iEvent.emplace(m_brlCandPutToken, std::move(candBarell));
    iEvent.emplace(m_fwdCandPutToken, std::move(candForward));

    return;
  }

  if (rpcDigis->begin() == rpcDigis->end()) {
    LogDebug("RPCTrigger") << "\nWarning: RPCDigiCollection with input tag " << m_label
                           << "\n seems to be empty for this event. Emulator will run on empty collection ";
  }

  std::vector<RPCDigiL1Link> brlLinks;
  std::vector<RPCDigiL1Link> fwdLinks;

  edm::ESHandle<L1RPCConeBuilder> coneBuilder = iSetup.getHandle(m_coneBuilderToken);

  edm::ESHandle<L1RPCConeDefinition> l1RPCConeDefinition = iSetup.getHandle(m_coneDefinitionToken);

  edm::ESHandle<L1RPCHwConfig> hwConfig = iSetup.getHandle(m_hwConfigToken);

  edm::ESHandle<L1RPCBxOrConfig> bxOrConfig = iSetup.getHandle(m_BxOrConfigToken);

  edm::ESHandle<L1RPCHsbConfig> hsbConfig = iSetup.getHandle(m_hsbConfigToken);

  for (int iBx = -1; iBx < 2; ++iBx) {
    L1RpcLogConesVec ActiveCones =
        m_theLinksystemFromES.getConesFromES(rpcDigis, coneBuilder, l1RPCConeDefinition, bxOrConfig, hwConfig, iBx);

    L1RpcTBMuonsVec2 finalMuons = m_pacTrigger->runEvent(ActiveCones, hsbConfig);

    //int maxFiredPlanes = 0;

    /*
    for (unsigned int i=0;i<ActiveCones.size();i++){
        int fpCnt = ActiveCones[i].getFiredPlanesCnt();
        if (fpCnt > maxFiredPlanes)
          maxFiredPlanes = fpCnt;
    }
  
    // Fill out the products
    // finalMuons[0]=barell, finalMuons[1]=endcap
    LogDebug("RPCTrigger") << "---Filling candindates in new event--- " 
                          << maxFiredPlanes << std::endl;
<<<<<<< RPCTrigger.cc
    */

    //    std::vector<L1MuRegionalCand> RPCb = giveFinallCandindates(finalMuons[0],1, iBx);
    //    std::vector<L1MuRegionalCand> RPCf = giveFinallCandindates(finalMuons[1],3, iBx);
    std::vector<RPCDigiL1Link> dlBrl;
    std::vector<RPCDigiL1Link> dlFwd;
    std::vector<L1MuRegionalCand> RPCb = giveFinallCandindates(finalMuons[0], 1, iBx, rpcDigis, dlBrl);
    std::vector<L1MuRegionalCand> RPCf = giveFinallCandindates(finalMuons[1], 3, iBx, rpcDigis, dlFwd);

    brlLinks.insert(brlLinks.end(), dlBrl.begin(), dlBrl.end());
    fwdLinks.insert(fwdLinks.end(), dlFwd.begin(), dlFwd.end());

    candBarell.insert(candBarell.end(), RPCb.begin(), RPCb.end());
    candForward.insert(candForward.end(), RPCf.begin(), RPCf.end());

    if (m_triggerDebug == 1)
      MuonsGrabber::Instance().writeDataForRelativeBX(iBx);
  }

  iEvent.emplace(m_fwdLinksPutToken, std::move(fwdLinks));
  iEvent.emplace(m_brlLinksPutToken, std::move(brlLinks));
  iEvent.emplace(m_brlCandPutToken, std::move(candBarell));
  iEvent.emplace(m_fwdCandPutToken, std::move(candForward));
}
///////////////////////////////////////////////////////////////////////////////
/**
 *
 * \brief Returns vector of L1MuRegionalCand (input of L1GMT)
 * \note - type is defined in L1MuRegionalCand 1 - barell, 3 - forward
 * \todo - we use offset value of 5 deegres. It should be stored centrally.
 *
 */
///////////////////////////////////////////////////////////////////////////////
std::vector<L1MuRegionalCand> RPCTrigger::giveFinallCandindates(const L1RpcTBMuonsVec& finalMuons,
                                                                int type,
                                                                int bx,
                                                                edm::Handle<RPCDigiCollection> rpcDigis,
                                                                std::vector<RPCDigiL1Link>& retRPCDigiLink) {
  std::vector<L1MuRegionalCand> RPCCand;

  for (unsigned int iMu = 0; iMu < finalMuons.size(); iMu++) {
    if (finalMuons[iMu].getPtCode() == 0) {
      continue;
    }

    RPCDigiL1Link newDigiLink;

    //std::cout << "######################################## " << std::endl;
    //std::cout << finalMuons[iMu].getPhiAddr() << " " << finalMuons[iMu].getEtaAddr() << std::endl;
    RPCMuon::TDigiLinkVec digiIVec = finalMuons[iMu].getDigiIdxVec();
    // Here the iteration has to be the same as in
    short int digiIndex = 0;
    RPCDigiCollection::DigiRangeIterator detUnitIt;
    for (detUnitIt = rpcDigis->begin(); detUnitIt != rpcDigis->end(); ++detUnitIt) {
      const RPCDetId& id = (*detUnitIt).first;
      uint32_t rawId = id.rawId();
      const RPCDigiCollection::Range& range = (*detUnitIt).second;

      for (RPCDigiCollection::const_iterator digiIt = range.first; digiIt != range.second; ++digiIt) {
        ++digiIndex;

        RPCMuon::TDigiLinkVec::iterator it = digiIVec.begin();
        for (; it != digiIVec.end(); ++it) {
          if (digiIndex == it->m_digiIdx) {
            newDigiLink.setLink(it->m_layer + 1, rawId, digiIt->strip(), digiIt->bx());
            //std::cout << type << " " << iMu << " layer: " << it->m_layer <<  " index " << it->m_digiIdx << std::endl;
            //std::cout << "         " << id << " " << " |bx " << digiIt->bx() << " strip " << digiIt->strip() << std::endl;
          }
        }
      }
    }
    retRPCDigiLink.push_back(newDigiLink);

    L1MuRegionalCand l1Cand;

    l1Cand.setBx(bx);

    l1Cand.setQualityPacked(finalMuons[iMu].getQuality());
    l1Cand.setPtPacked(finalMuons[iMu].getPtCode());

    l1Cand.setType(type);

    int charge = finalMuons[iMu].getSign();

    if (charge == 0)  // negative
      l1Cand.setChargePacked(1);
    else
      l1Cand.setChargePacked(0);

    //RPCConst::l1RpcConeCrdnts cone = finalMuons[iMu].getConeCrdnts();

    /*
    int pac = cone.m_LogSector*12+cone.m_LogSegment;
    const float pi = 3.14159265;
    const float offset = 5*(2*pi/360); // redefinition! Defined also in RPCRingFromRolls::phiMapCompare
    float phi = 2*pi*pac/144-offset;
    if (phi<0)
      phi+=2*pi;
    
    l1Cand.setPhiValue(phi);
    */

    //Note: pac numbering begins at 5 deg and goes from 1 to 144.
    // we want phi values from 0 to 2.5 deg to be phiPacked=0
    // max phiPacked value is 143 (see CMS IN 2004-022)
    //int phiPacked = (finalMuons[iMu].getPhiAddr()+2)%144;
    int phiPacked = finalMuons[iMu].getPhiAddr();
    l1Cand.setPhiPacked(phiPacked);
    /*
    float eta = RPCConst::etaFromTowerNum(cone.m_Tower);
    l1Cand.setEtaValue(eta);
*/
    //Note: etaAddr is packed in special way: see CMS IN 2004-022
    signed short etaAddr = finalMuons[iMu].getEtaAddr();  //
    //    signed short etaAddr = finalMuons[iMu].getEtaAddr()-16; // -16..16
    //    bool etaNegative = false;
    //    if (etaAddr < 0){
    //      etaNegative = true;
    //      etaAddr = ~(-etaAddr)+1; // convert to negative :)
    //    }

    //    etaAddr &= 63; // 6 bits only

    l1Cand.setEtaPacked(etaAddr);
    l1Cand.setChargeValid(true);

    /*    
    std::cout<< std::endl << "RBMuon::" << finalMuons[iMu].getEtaAddr() << " " 
             << finalMuons[iMu].getPhiAddr() << std::endl ;
    std::cout<< "cand " <<  l1Cand.eta_packed() << " " 
             << l1Cand.phi_packed() << std::endl ;
    */

    RPCCand.push_back(l1Cand);

    LogDebug("RPCTrigger") << "Found muonf of pt " << finalMuons[iMu].getPtCode() << " bx " << l1Cand.bx()
                           << " L1Charge " << l1Cand.charge_packed() << " ql " << l1Cand.quality() << " fp "
                           << finalMuons[iMu].getFiredPlanes() << " b/f " << l1Cand.type_idx() << " phi "
                           << l1Cand.phi_packed() << " eta "
                           << l1Cand.eta_packed()
                           //<< " eta l1 " << l1Cand.etaValue() // will drop out soon
                           << " killed " << finalMuons[iMu].wasKilled();
  }

  return RPCCand;
}