OMTFinputMaker Class Reference

#include <OMTFinputMaker.h>

Inheritance diagram for OMTFinputMaker:

Public Member Functions
OMTFinput	buildInputForProcessor (const L1MuDTChambPhContainer *dtPhDigis, const L1MuDTChambThContainer *dtThDigis, const CSCCorrelatedLCTDigiCollection *cscDigis, const RPCDigiCollection *rpcDigis, unsigned int iProcessor, l1t::tftype type=l1t::tftype::omtf_pos, int bx=0)
	Method translating trigger digis into input matrix with global phi coordinates. More...
int	getFlag () const
void	initialize (const edm::EventSetup &, const OMTFConfiguration *)
void	initialize (const edm::ParameterSet &edmCfg, const edm::EventSetup &es, const MuonGeometryTokens &muonGeometryTokens) override
	OMTFinputMaker (edm::ConsumesCollector &, bool getDuringEvent=true)
	OMTFinputMaker (const edm::ParameterSet &edmParameterSet, MuStubsInputTokens &muStubsInputTokens, const OMTFConfiguration *config, std::unique_ptr< OmtfAngleConverter > angleConverter)
void	setFlag (int aFlag)
	~OMTFinputMaker ()
	~OMTFinputMaker () override
Public Member Functions inherited from MuonStubMakerBase
void	buildInputForProcessor (MuonStubPtrs2D &muonStubsInLayers, unsigned int iProcessor, l1t::tftype procTyp, int bxFrom=0, int bxTo=0)
	Method translating trigger digis into input matrix with global phi coordinates, fills the muonStubsInLayers. More...
void	loadAndFilterDigis (const edm::Event &event)
	MuonStubMakerBase (const ProcConfigurationBase *procConf)
virtual	~MuonStubMakerBase ()

Static Public Member Functions
static bool	acceptDtDigi (const OMTFConfiguration *config, const DTChamberId &dTChamberId, unsigned int iProcessor, l1t::tftype procType)
static void	addStub (const OMTFConfiguration *config, MuonStubPtrs2D &muonStubsInLayers, unsigned int iLayer, unsigned int iInput, MuonStub &stub)
static unsigned int	getInputNumber (const OMTFConfiguration *config, unsigned int rawId, unsigned int iProcessor, l1t::tftype type)
static int	getProcessorPhiZero (const OMTFConfiguration *config, unsigned int iProcessor)

Protected Attributes
std::unique_ptr< OmtfAngleConverter >	angleConverter
const OMTFConfiguration *	config = nullptr
Protected Attributes inherited from MuonStubMakerBase
const ProcConfigurationBase *	config = nullptr
std::vector< std::unique_ptr< DigiToStubsConverterBase > >	digiToStubsConverters
RpcClusterization	rpcClusterization

Private Member Functions
bool	acceptDigi (uint32_t rawId, unsigned int iProcessor, l1t::tftype type)
unsigned int	getInputNumber (unsigned int rawId, unsigned int iProcessor, l1t::tftype type)
OMTFinput	processCSC (const CSCCorrelatedLCTDigiCollection *cscDigis, unsigned int iProcessor, l1t::tftype type, int bx)
OMTFinput	processDT (const L1MuDTChambPhContainer *dtPhDigis, const L1MuDTChambThContainer *dtThDigis, unsigned int iProcessor, l1t::tftype type, int bx)
OMTFinput	processRPC (const RPCDigiCollection *rpcDigis, unsigned int iProcessor, l1t::tftype type, int bx)

Private Attributes
int	flag
std::unique_ptr< AngleConverter >	myAngleConverter
const OMTFConfiguration *	myOmtfConfig

Detailed Description

Definition at line 19 of file OMTFinputMaker.h.

Constructor & Destructor Documentation

OMTFinputMaker() [1/2]

OMTFinputMaker::OMTFinputMaker	(	edm::ConsumesCollector &	iC,
		bool	getDuringEvent = true
	)

Definition at line 20 of file OMTFinputMaker.cc.

    : myAngleConverter(std::make_unique<AngleConverter>(iC, getDuringEvent)) {}

~OMTFinputMaker() [1/2]

OMTFinputMaker::~OMTFinputMaker

(

)

Definition at line 31 of file OMTFinputMaker.cc.

{}

OMTFinputMaker() [2/2]

OMTFinputMaker::OMTFinputMaker	(	const edm::ParameterSet &	edmParameterSet,
		MuStubsInputTokens &	muStubsInputTokens,
		const OMTFConfiguration *	config,
		std::unique_ptr< OmtfAngleConverter >	angleConverter
	)

Definition at line 311 of file OMTFinputMaker.cc.

References MuonStubMakerBase::digiToStubsConverters, edm::ParameterSet::getParameter(), MuStubsInputTokens::inputTokenCSC, MuStubsInputTokens::inputTokenDtPh, MuStubsInputTokens::inputTokenDtTh, MuStubsInputTokens::inputTokenRPC, and MuonStubMakerBase::rpcClusterization.

    : MuonStubMakerBase(config), config(config), angleConverter(std::move(angleConverter)) {
  if (!edmParameterSet.getParameter<bool>("dropDTPrimitives"))
    digiToStubsConverters.emplace_back(std::make_unique<DtDigiToStubsConverterOmtf>(
        config, this->angleConverter.get(), muStubsInputTokens.inputTokenDtPh, muStubsInputTokens.inputTokenDtTh));

  if (!edmParameterSet.getParameter<bool>("dropCSCPrimitives"))
    digiToStubsConverters.emplace_back(std::make_unique<CscDigiToStubsConverterOmtf>(
        config, this->angleConverter.get(), muStubsInputTokens.inputTokenCSC));

  if (!edmParameterSet.getParameter<bool>("dropRPCPrimitives"))
    digiToStubsConverters.emplace_back(std::make_unique<RpcDigiToStubsConverterOmtf>(
        config, this->angleConverter.get(), &rpcClusterization, muStubsInputTokens.inputTokenRPC));
}

~OMTFinputMaker() [2/2]

OMTFinputMaker::~OMTFinputMaker ( )

override

Member Function Documentation

acceptDigi()

bool OMTFinputMaker::acceptDigi ( uint32_t  rawId, unsigned int  iProcessor, l1t::tftype  type  )

private

Check if digis are within a give processor input. Simply checks sectors range.

Clean up digis. Remove unconnected detectors

Select RPC chambers connected to OMTF

RPC RE1/2 temporarily not used (aId.region()==1 && aId.station()==1 && aId.ring()<2) ||

Definition at line 34 of file OMTFinputMaker.cc.

References l1t::bmtf, MuonSubdetId::CSC, omtf::DataWord64::csc, DetId::det(), MuonSubdetId::DT, dt, l1t::emtf_neg, l1t::emtf_pos, OMTFConfiguration::getBarrelMax(), OMTFConfiguration::getBarrelMin(), OMTFConfiguration::getEndcap10DegMax(), OMTFConfiguration::getEndcap10DegMin(), OMTFConfiguration::getEndcap20DegMax(), OMTFConfiguration::getEndcap20DegMin(), RPCDetId::layer(), DetId::Muon, myOmtfConfig, l1t::omtf_neg, l1t::omtf_pos, RPCDetId::region(), RPCDetId::ring(), RPCDetId::roll(), MuonSubdetId::RPC, RPCDetId::sector(), RPCDetId::station(), DetId::subdetId(), and RPCDetId::subsector().

Referenced by processCSC(), processDT(), and processRPC().

                                                                                       {
  unsigned int aMin = myOmtfConfig->getBarrelMin()[iProcessor];
  unsigned int aMax = myOmtfConfig->getBarrelMax()[iProcessor];
  unsigned int aSector = 99;

  DetId detId(rawId);
  if (detId.det() != DetId::Muon)
    edm::LogError("Critical OMTFinputMaker") << "PROBLEM: hit in unknown Det, detID: " << detId.det() << std::endl;
  switch (detId.subdetId()) {
    case MuonSubdetId::RPC: {
      RPCDetId aId(rawId);

      if (type == l1t::tftype::omtf_pos &&
          (aId.region() < 0 || (aId.region() == 0 && aId.ring() != 2) || (aId.region() == 0 && aId.station() == 4) ||
           (aId.region() == 0 && aId.station() == 2 && aId.layer() == 2 && aId.roll() == 1) ||
           (aId.region() == 0 && aId.station() == 3 && aId.roll() == 1) || (aId.region() == 1 && aId.station() == 4) ||
           (aId.region() == 1 && aId.station() > 0 && aId.ring() < 3)))
        return false;

      if (type == l1t::tftype::omtf_neg &&
          (aId.region() > 0 || (aId.region() == 0 && aId.ring() != -2) || (aId.region() == 0 && aId.station() == 4) ||
           (aId.region() == 0 && aId.station() == 2 && aId.layer() == 2 && aId.roll() == 1) ||
           (aId.region() == 0 && aId.station() == 3 && aId.roll() == 1) || (aId.region() == -1 && aId.station() == 4) ||
           //RPC RE1/2 temporarily not used (aId.region()==1 && aId.station()==1 && aId.ring()<2) ||
           (aId.region() == -1 && aId.station() > 0 && aId.ring() < 3)))
        return false;

      if (type == l1t::tftype::bmtf && aId.region() != 0)
        return false;

      if (type == l1t::tftype::emtf_pos && (aId.region() <= 0 || (aId.station() == 1 && aId.ring() == 3)))
        return false;
      if (type == l1t::tftype::emtf_neg && (aId.region() >= 0 || (aId.station() == 1 && aId.ring() == 3)))
        return false;
      if (aId.region() == 0)
        aSector = aId.sector();
      if (aId.region() != 0) {
        aSector = (aId.sector() - 1) * 6 + aId.subsector();
        aMin = myOmtfConfig->getEndcap10DegMin()[iProcessor];
        aMax = myOmtfConfig->getEndcap10DegMax()[iProcessor];
      }

      break;
    }
    case MuonSubdetId::DT: {
      DTChamberId dt(rawId);

      if (type == l1t::tftype::omtf_pos && dt.wheel() != 2)
        return false;
      if (type == l1t::tftype::omtf_neg && dt.wheel() != -2)
        return false;
      if (type == l1t::tftype::emtf_pos || type == l1t::tftype::emtf_neg)
        return false;

      aSector = dt.sector();
      break;
    }
    case MuonSubdetId::CSC: {
      CSCDetId csc(rawId);
      if (type == l1t::tftype::omtf_pos && (csc.endcap() == 2 || csc.ring() == 1 || csc.station() == 4))
        return false;
      if (type == l1t::tftype::omtf_neg && (csc.endcap() == 1 || csc.ring() == 1 || csc.station() == 4))
        return false;

      if (type == l1t::tftype::emtf_pos && (csc.endcap() == 2 || (csc.station() == 1 && csc.ring() == 3)))
        return false;
      if (type == l1t::tftype::emtf_neg && (csc.endcap() == 1 || (csc.station() == 1 && csc.ring() == 3)))
        return false;

      aSector = csc.chamber();
      aMin = myOmtfConfig->getEndcap10DegMin()[iProcessor];
      aMax = myOmtfConfig->getEndcap10DegMax()[iProcessor];

      if ((type == l1t::tftype::emtf_pos || type == l1t::tftype::emtf_neg) && csc.station() > 1 && csc.ring() == 1) {
        aMin = myOmtfConfig->getEndcap20DegMin()[iProcessor];
        aMax = myOmtfConfig->getEndcap20DegMax()[iProcessor];
      }
      break;
    }
  }

  if (aMax > aMin && aSector >= aMin && aSector <= aMax)
    return true;
  if (aMax < aMin && (aSector >= aMin || aSector <= aMax))
    return true;

  return false;
}

acceptDtDigi()

bool OMTFinputMaker::acceptDtDigi ( const OMTFConfiguration *  config, const DTChamberId &  dTChamberId, unsigned int  iProcessor, l1t::tftype  procType  )

static

Definition at line 285 of file OMTFinputMaker.cc.

References l1t::emtf_neg, l1t::emtf_pos, l1t::omtf_neg, l1t::omtf_pos, DTChamberId::sector(), and DTChamberId::wheel().

Referenced by DtDigiToStubsConverterOmtf::acceptDigi().

                                                      {
  unsigned int aMin = config->getBarrelMin()[iProcessor];
  unsigned int aMax = config->getBarrelMax()[iProcessor];

  if (procType == l1t::tftype::omtf_pos && dTChamberId.wheel() != 2)
    return false;
  if (procType == l1t::tftype::omtf_neg && dTChamberId.wheel() != -2)
    return false;
  if (procType == l1t::tftype::emtf_pos || procType == l1t::tftype::emtf_neg)
    return false;

  unsigned int aSector = dTChamberId.sector();

  if (aMax > aMin && aSector >= aMin && aSector <= aMax)
    return true;
  if (aMax < aMin && (aSector >= aMin || aSector <= aMax))
    return true;

  return false;
}

addStub()

void OMTFinputMaker::addStub ( const OMTFConfiguration *  config, MuonStubPtrs2D &  muonStubsInLayers, unsigned int  iLayer, unsigned int  iInput, MuonStub &  stub  )

static

Definition at line 440 of file OMTFinputMaker.cc.

References MuonStub::etaHw, LogTrace, MuonStub::phiBHw, and MuonStub::phiHw.

Referenced by CscDigiToStubsConverterOmtf::addCSCstubs(), DtDigiToStubsConverterOmtf::addDTphiDigi(), and RpcDigiToStubsConverterOmtf::addRPCstub().

                                             {
  //LogTrace("l1tOmtfEventPrint") << __FUNCTION__ << ":" << __LINE__ << " iInput " << iInput << " " << stub << endl;
  //there is a small rate of duplicated digis in the real data in the DT and CSC, the reason for this duplicates is not understood
  //in case of RPC the duplicated digis appear in data only for the layer 17 (RE3), where the rolls 2 and 3 has the same eta = 115 assigned, and the muon can hit both rolls
  //the reason of the  duplicates cannot be the fact that the same data (links) goes to neighboring boards, because this is  removed in the OMTF unpacker
  //the duplicates cannot be dropped, because if they appear in data, are also duplicated on the input of the algorithm, and both can be used as the reference hits
  if (muonStubsInLayers[iLayer][iInput] && muonStubsInLayers[iLayer][iInput]->phiHw == stub.phiHw &&
      muonStubsInLayers[iLayer][iInput]->phiBHw == stub.phiBHw &&
      muonStubsInLayers[iLayer][iInput]->etaHw == stub.etaHw) {
    LogTrace("OMTFReconstruction") << "addStub: the stub with exactly the same phi, phiB and eta was already added:\n"
                                   << "incomnig stub " << stub << "\n"
                                   << "existing stub " << *(muonStubsInLayers[iLayer][iInput]) << std::endl;
    //return;
  }

  if (muonStubsInLayers[iLayer][iInput] && muonStubsInLayers[iLayer][iInput]->phiHw != (int)config->nPhiBins())
    ++iInput;

  if (muonStubsInLayers[iLayer][iInput] && muonStubsInLayers[iLayer][iInput]->phiHw != (int)config->nPhiBins())
    return;
  //in this implementation only two first stubs are added for a given iInput

  muonStubsInLayers.at(iLayer).at(iInput) = std::make_shared<MuonStub>(stub);
}

buildInputForProcessor()

OMTFinput OMTFinputMaker::buildInputForProcessor	(	const L1MuDTChambPhContainer *	dtPhDigis,
		const L1MuDTChambThContainer *	dtThDigis,
		const CSCCorrelatedLCTDigiCollection *	cscDigis,
		const RPCDigiCollection *	rpcDigis,
		unsigned int	iProcessor,
		l1t::tftype	type = l1t::tftype::omtf_pos,
		int	bx = 0
	)

Method translating trigger digis into input matrix with global phi coordinates.

Definition at line 383 of file OMTFinputMaker.cc.

References simKBmtfDigis_cfi::bx, myOmtfConfig, processCSC(), processDT(), processRPC(), and mps_fire::result.

Referenced by OMTFPatternMaker::analyze(), OMTFReconstruction::getProcessorCandidates(), and OMTFProcessor< GoldenPatternType >::run().

                                                         {
  OMTFinput result(myOmtfConfig);
  result += processDT(dtPhDigis, dtThDigis, iProcessor, type, bx);
  result += processCSC(cscDigis, iProcessor, type, bx);
  result += processRPC(rpcDigis, iProcessor, type, bx);
  return result;
}

getFlag()

int OMTFinputMaker::getFlag ( ) const

inline

Definition at line 37 of file OMTFinputMaker.h.

References flag.

Referenced by OMTFReconstruction::getProcessorCandidates().

{ return flag; }

getInputNumber() [1/2]

unsigned int OMTFinputMaker::getInputNumber ( unsigned int  rawId, unsigned int  iProcessor, l1t::tftype  type  )

private

Give input number for givedn processor, using the chamber sector number. Result is modulo allowed number of hits per chamber

on the 0-2pi border we need to add 1 30 deg sector to get the correct index








Set roll number by hand to keep common input number shift formula for all stations







Only one roll from station 3 is connected.

At the moment do not use RPC chambers splitting into rolls for bmtf part

on the 0-2pi border we need to add 4 10 deg sectors to get the correct index




on the 0-2pi border we need to add 1 30 deg sector to get the correct index



on the 0-2pi border we need to add 4 10deg sectors to get the correct index


Endcap region covers algo 10 deg sectors on the 0-2pi border we need to add 2 20deg sectors to get the correct index

Assume 2 hits per chamber

Chambers divided into two rolls have rolls number 1 and 3

Definition at line 128 of file OMTFinputMaker.cc.

References l1t::bmtf, MuonSubdetId::CSC, omtf::DataWord64::csc, DetId::det(), MuonSubdetId::DT, dt, l1t::emtf_neg, l1t::emtf_pos, OMTFConfiguration::getBarrelMin(), OMTFConfiguration::getEndcap10DegMin(), OMTFConfiguration::getEndcap20DegMin(), DetId::Muon, myOmtfConfig, MuonSubdetId::RPC, omtf::DataWord64::rpc, and DetId::subdetId().

Referenced by RpcDigiToStubsConverterOmtf::acceptDigi(), CscDigiToStubsConverterOmtf::addCSCstubs(), DtDigiToStubsConverterOmtf::addDTphiDigi(), RpcDigiToStubsConverterOmtf::addRPCstub(), processCSC(), processDT(), and processRPC().

                                                                                                       {
  unsigned int iInput = 99;
  unsigned int aSector = 99;
  int aMin = myOmtfConfig->getBarrelMin()[iProcessor];
  int iRoll = 1;
  int nInputsPerSector = 2;

  DetId detId(rawId);
  if (detId.det() != DetId::Muon)
    edm::LogError("Critical OMTFinputMaker") << "PROBLEM: hit in unknown Det, detID: " << detId.det() << std::endl;
  switch (detId.subdetId()) {
    case MuonSubdetId::RPC: {
      RPCDetId rpc(rawId);
      if (rpc.region() == 0) {
        nInputsPerSector = 4;
        aSector = rpc.sector();
        if (iProcessor == 5 && aSector < 3)
          aMin = -1;
        //Use division into rolls
        iRoll = rpc.roll();
        if (rpc.station() == 2 && rpc.layer() == 2 && rpc.roll() == 2)
          iRoll = 1;
        if (rpc.station() == 3) {
          iRoll = 1;
          nInputsPerSector = 2;
        }
        if (type == l1t::tftype::bmtf)
          iRoll = 1;
      }
      if (rpc.region() != 0) {
        aSector = (rpc.sector() - 1) * 6 + rpc.subsector();
        aMin = myOmtfConfig->getEndcap10DegMin()[iProcessor];
        if (iProcessor == 5 && aSector < 5)
          aMin = -4;
      }
      break;
    }
    case MuonSubdetId::DT: {
      DTChamberId dt(rawId);
      aSector = dt.sector();
      if (iProcessor == 5 && aSector < 3)
        aMin = -1;
      break;
    }
    case MuonSubdetId::CSC: {
      CSCDetId csc(rawId);
      aSector = csc.chamber();
      aMin = myOmtfConfig->getEndcap10DegMin()[iProcessor];
      if (iProcessor == 5 && aSector < 5)
        aMin = -4;
      if ((type == l1t::tftype::emtf_pos || type == l1t::tftype::emtf_neg) && csc.station() > 1 && csc.ring() == 1) {
        aMin = myOmtfConfig->getEndcap20DegMin()[iProcessor];
        if (iProcessor == 5 && aSector < 3)
          aMin = -2;
      }
      break;
    }
  }

  iInput = (aSector - aMin) * nInputsPerSector;
  iInput += iRoll - 1;

  return iInput;
}

getInputNumber() [2/2]

unsigned int OMTFinputMaker::getInputNumber ( const OMTFConfiguration *  config, unsigned int  rawId, unsigned int  iProcessor, l1t::tftype  type  )

static

Give input number for given processor, using the chamber sector number. Result is modulo allowed number of hits per chamber

on the 0-2pi border we need to add 1 30 deg sector to get the correct index








Set roll number by hand to keep common input number shift formula for all stations







Only one roll from station 3 is connected.

At the moment do not use RPC chambers splitting into rolls for bmtf part

on the 0-2pi border we need to add 4 10 deg sectors to get the correct index




on the 0-2pi border we need to add 1 30 deg sector to get the correct index



on the 0-2pi border we need to add 4 10deg sectors to get the correct index


Endcap region covers algo 10 deg sectors on the 0-2pi border we need to add 2 20deg sectors to get the correct index

Assume 2 hits per chamber

Chambers divided into two rolls have rolls number 1 and 3

Definition at line 341 of file OMTFinputMaker.cc.

References l1t::bmtf, MuonSubdetId::CSC, omtf::DataWord64::csc, DetId::det(), MuonSubdetId::DT, dt, l1t::emtf_neg, l1t::emtf_pos, DetId::Muon, MuonSubdetId::RPC, omtf::DataWord64::rpc, and DetId::subdetId().

                                                            {
  unsigned int iInput = 99;
  unsigned int aSector = 99;
  int aMin = config->getBarrelMin()[iProcessor];
  int iRoll = 1;
  int nInputsPerSector = 2;

  DetId detId(rawId);
  if (detId.det() != DetId::Muon)
    edm::LogError("Critical OMTFinputMaker") << "PROBLEM: hit in unknown Det, detID: " << detId.det() << std::endl;
  switch (detId.subdetId()) {
    case MuonSubdetId::RPC: {
      RPCDetId rpc(rawId);
      if (rpc.region() == 0) {
        nInputsPerSector = 4;
        aSector = rpc.sector();
        if (iProcessor == 5 && aSector < 3)
          aMin = -1;
        //Use division into rolls
        iRoll = rpc.roll();
        if (rpc.station() == 2 && rpc.layer() == 2 && rpc.roll() == 2)
          iRoll = 1;
        if (rpc.station() == 3) {
          iRoll = 1;
          nInputsPerSector = 2;
        }
        if (type == l1t::tftype::bmtf)
          iRoll = 1;
      }
      if (rpc.region() != 0) {
        aSector = (rpc.sector() - 1) * 6 + rpc.subsector();
        aMin = config->getEndcap10DegMin()[iProcessor];
        if (iProcessor == 5 && aSector < 5)
          aMin = -4;
      }
      break;
    }
    case MuonSubdetId::DT: {
      DTChamberId dt(rawId);
      aSector = dt.sector();
      if (iProcessor == 5 && aSector < 3)
        aMin = -1;
      break;
    }
    case MuonSubdetId::CSC: {
      CSCDetId csc(rawId);
      aSector = csc.chamber();
      aMin = config->getEndcap10DegMin()[iProcessor];
      if (iProcessor == 5 && aSector < 5)
        aMin = -4;
      if ((type == l1t::tftype::emtf_pos || type == l1t::tftype::emtf_neg) && csc.station() > 1 && csc.ring() == 1) {
        aMin = config->getEndcap20DegMin()[iProcessor];
        if (iProcessor == 5 && aSector < 3)
          aMin = -2;
      }
      break;
    }
  }

  iInput = (aSector - aMin) * nInputsPerSector;
  iInput += iRoll - 1;

  return iInput;
}

getProcessorPhiZero()

int OMTFinputMaker::getProcessorPhiZero ( const OMTFConfiguration *  config, unsigned int  iProcessor  )

static

iProcessor - from 0 to 5 returns the global phi in hardware scale (myOmtfConfig->nPhiBins() ) at which the scale starts for give processor

Definition at line 429 of file OMTFinputMaker.cc.

References ecaldqm::binning::nPhiBins, and l1ctLayer1_cff::phiZero.

Referenced by CscDigiToStubsConverterOmtf::addCSCstubs(), DtDigiToStubsConverterOmtf::addDTphiDigi(), RpcDigiToStubsConverterOmtf::addRPCstub(), StubsSimHitsMatcher::match(), and EventCapture::observeEventEnd().

                                                                                                {
  unsigned int nPhiBins = config->nPhiBins();

  int phiZero = nPhiBins / 6. * (iProcessor) + nPhiBins / 24;
  // "0" is 15degree moved cyclically to each processor, note [0,2pi]

  return config->foldPhi(phiZero);
}

initialize() [1/2]

void OMTFinputMaker::initialize	(	const edm::EventSetup &	es,
		const OMTFConfiguration *	omtfConfig
	)

Definition at line 24 of file OMTFinputMaker.cc.

References myAngleConverter, myOmtfConfig, and OMTFConfiguration::nPhiBins().

Referenced by OMTFPatternMaker::analyze(), and OMTFReconstruction::beginRun().

                                                                                            {
  myAngleConverter->checkAndUpdateGeometry(es, omtfConfig->nPhiBins());

  myOmtfConfig = omtfConfig;
}

initialize() [2/2]

void OMTFinputMaker::initialize ( const edm::ParameterSet &  edmCfg, const edm::EventSetup &  es, const MuonGeometryTokens &  muonGeometryTokens  )

overridevirtual

Reimplemented from MuonStubMakerBase.

Definition at line 329 of file OMTFinputMaker.cc.

References angleConverter, and MuonStubMakerBase::initialize().

                                                                              {
  MuonStubMakerBase::initialize(edmCfg, es, muonGeometryTokens);
  angleConverter->checkAndUpdateGeometry(es, config, muonGeometryTokens);
}

processCSC()

OMTFinput OMTFinputMaker::processCSC ( const CSCCorrelatedLCTDigiCollection *  cscDigis, unsigned int  iProcessor, l1t::tftype  type, int  bx  )

private

Take the CSC digis, select chambers connected to given processor, convers logal angles to global scale. For CSC do NOT take the bending angle.

Check it the data fits into given processor input range

Check if LCT trigger primitive has the right BX.

Accept CSC digis only up to eta=1.26. The nominal OMTF range is up to 1.24, but cutting at 1.24 kill efficnency at the edge. 1.26 is one eta bin above nominal.

Definition at line 263 of file OMTFinputMaker.cc.

References acceptDigi(), relativeConstraints::chamber, OMTFConfiguration::getHwToLogicLayer(), getInputNumber(), OMTFConfiguration::getLayerNumber(), l1tTowerCalibrationProducer_cfi::iEta, CSCConstants::LCT_CENTRAL_BX, myAngleConverter, myOmtfConfig, and mps_fire::result.

Referenced by buildInputForProcessor().

                                                {
  OMTFinput result(myOmtfConfig);
  if (!cscDigis)
    return result;

  auto chamber = cscDigis->begin();
  auto chend = cscDigis->end();
  for (; chamber != chend; ++chamber) {
    unsigned int rawid = (*chamber).first;
    if (!acceptDigi(rawid, iProcessor, type))
      continue;
    auto digi = (*chamber).second.first;
    auto dend = (*chamber).second.second;
    for (; digi != dend; ++digi) {
      if (digi->getBX() - CSCConstants::LCT_CENTRAL_BX != bxTrg)
        continue;

      unsigned int hwNumber = myOmtfConfig->getLayerNumber(rawid);
      if (myOmtfConfig->getHwToLogicLayer().find(hwNumber) == myOmtfConfig->getHwToLogicLayer().end())
        continue;

      unsigned int iLayer = myOmtfConfig->getHwToLogicLayer().at(hwNumber);
      int iPhi = myAngleConverter->getProcessorPhi(iProcessor, type, CSCDetId(rawid), *digi);
      int iEta = myAngleConverter->getGlobalEta(rawid, *digi);
      //if(abs(iEta)>1.26/2.61*240) continue;
      //if (abs(iEta) > 115) continue;
      unsigned int iInput = getInputNumber(rawid, iProcessor, type);
      //    std::cout <<" ADDING CSC hit, proc: "<<iProcessor<<" iPhi : " << iPhi <<" iEta: "<< iEta << std::endl;
      bool allowOverwrite = false;
      result.addLayerHit(iLayer, iInput, iPhi, iEta, allowOverwrite);
    }
  }
  return result;
}

processDT()

OMTFinput OMTFinputMaker::processDT ( const L1MuDTChambPhContainer *  dtPhDigis, const L1MuDTChambThContainer *  dtThDigis, unsigned int  iProcessor, l1t::tftype  type, int  bx  )

private

Take the DT digis, select chambers connected to given processor, convers logal angles to global scale. For DT take also the bending angle.

Check it the data fits into given processor input range

Check Trigger primitive quality Ts2Tag() == 0 - take only first track from DT Trigger Server BxCnt() == 0 - ?? code()>=3 - take only double layer hits, HH, HL and LL

Definition at line 211 of file OMTFinputMaker.cc.

References acceptDigi(), OMTFConfiguration::fwVersion(), L1MuDTChambPhContainer::getContainer(), OMTFConfiguration::getHwToLogicLayer(), getInputNumber(), OMTFConfiguration::getLayerNumber(), l1tTowerCalibrationProducer_cfi::iEta, myAngleConverter, myOmtfConfig, and mps_fire::result.

Referenced by buildInputForProcessor().

                                               {
  OMTFinput result(myOmtfConfig);
  if (!dtPhDigis)
    return result;

  for (const auto &digiIt : *dtPhDigis->getContainer()) {
    DTChamberId detid(digiIt.whNum(), digiIt.stNum(), digiIt.scNum() + 1);

    if (!acceptDigi(detid.rawId(), iProcessor, type))
      continue;

    // FIXME (MK): at least Ts2Tag selection is not correct! Check it
    //    if (digiIt.bxNum()!= 0 || digiIt.BxCnt()!= 0 || digiIt.Ts2Tag()!= 0 || digiIt.code()<4) continue;

    if (digiIt.bxNum() != bxTrg)
      continue;

    if (myOmtfConfig->fwVersion() <= 4) {
      if (digiIt.code() != 4 && digiIt.code() != 5 && digiIt.code() != 6)
        continue;
    } else {
      if (digiIt.code() != 2 && digiIt.code() != 3 && digiIt.code() != 4 && digiIt.code() != 5 && digiIt.code() != 6)
        continue;
    }

    unsigned int hwNumber = myOmtfConfig->getLayerNumber(detid.rawId());
    if (myOmtfConfig->getHwToLogicLayer().find(hwNumber) == myOmtfConfig->getHwToLogicLayer().end())
      continue;

    auto iter = myOmtfConfig->getHwToLogicLayer().find(hwNumber);
    unsigned int iLayer = iter->second;
    int iPhi = myAngleConverter->getProcessorPhi(iProcessor, type, digiIt);
    int iEta = myAngleConverter->getGlobalEta(detid.rawId(), digiIt, dtThDigis);
    unsigned int iInput = getInputNumber(detid.rawId(), iProcessor, type);
    bool allowOverwrite = false;
    result.addLayerHit(iLayer, iInput, iPhi, iEta, allowOverwrite);
    result.addLayerHit(iLayer + 1, iInput, digiIt.phiB(), iEta, allowOverwrite);
  }

  return result;
}

processRPC()

OMTFinput OMTFinputMaker::processRPC ( const RPCDigiCollection *  rpcDigis, unsigned int  iProcessor, l1t::tftype  type, int  bx  )

private

Decluster nearby hits in single chamber, by taking average cluster position, expressed in half RPC strip: pos = (cluster_begin + cluster_end)

Find clusters of consecutive fired strips. Have to copy the digis in chamber to sort them (not optimal). NOTE: when copying I select only digis with bx== //FIXME: find a better place/way to filtering digi against quality/BX etc.

Definition at line 310 of file OMTFinputMaker.cc.

References funct::abs(), acceptDigi(), bsc_activity_cfg::clusters, OMTFConfiguration::getHwToLogicLayer(), getInputNumber(), OMTFConfiguration::getLayerNumber(), l1tTowerCalibrationProducer_cfi::iEta, myAngleConverter, myOmtfConfig, heavyionUCCDQM_cfi::nClusters, DetId::rawId(), mps_fire::result, rpcPrimitiveCmp(), jetUpdater_cfi::sort, and str.

Referenced by buildInputForProcessor().

                                                {
  OMTFinput result(myOmtfConfig);
  if (!rpcDigis)
    return result;
  std::stringstream str;

  //  std::cout <<" RPC HITS, processor : " << iProcessor << std::endl;

  const RPCDigiCollection &rpcDigiCollection = *rpcDigis;
  for (auto rollDigis : rpcDigiCollection) {
    RPCDetId roll = rollDigis.first;
    unsigned int rawid = roll.rawId();
    int nClusters = 0;
    if (!acceptDigi(rawid, iProcessor, type))
      continue;
    //  for (auto tdigi = rollDigis.second.first; tdigi != rollDigis.second.second; tdigi++) { std::cout << "RPC DIGIS: " << roll.rawId()<< " "<<roll<<" digi: " << tdigi->strip() <<" bx: " << tdigi->bx() << std::endl; }
    std::vector<RPCDigi> digisCopy;
    //  std::copy_if(rollDigis.second.first, rollDigis.second.second, std::back_inserter(digisCopy), [](const RPCDigi & aDigi){return (aDigi.bx()==0);} );
    for (auto pDigi = rollDigis.second.first; pDigi != rollDigis.second.second; pDigi++) {
      if (pDigi->bx() == bxTrg)
        digisCopy.push_back(*pDigi);
    }
    std::sort(digisCopy.begin(), digisCopy.end(), rpcPrimitiveCmp);
    typedef std::pair<unsigned int, unsigned int> Cluster;
    std::vector<Cluster> clusters;
    for (auto &digi : digisCopy) {
      if (clusters.empty())
        clusters.push_back(Cluster(digi.strip(), digi.strip()));
      else if (digi.strip() - clusters.back().second == 1)
        clusters.back().second = digi.strip();
      else if (digi.strip() - clusters.back().second > 1)
        clusters.push_back(Cluster(digi.strip(), digi.strip()));
    }

    for (auto &cluster : clusters) {
      //      int iPhiHalfStrip1 = myAngleConverter->getProcessorPhi(iProcessor, type, roll, cluster.first);
      //      int iPhiHalfStrip2 = myAngleConverter->getProcessorPhi(iProcessor, type, roll, cluster.second);
      int iPhi = myAngleConverter->getProcessorPhi(iProcessor, type, roll, cluster.first, cluster.second);
      int cSize = abs(int(cluster.first) - int(cluster.second)) + 1;
      //      std::cout << " HStrip_1: " << iPhiHalfStrip1 <<" HStrip_2: "<<iPhiHalfStrip2<<" iPhi: " << iPhi << " cluster: ["<< cluster.first << ", "<<  cluster.second <<"]"<< std::endl;
      if (cSize > 3)
        continue;
      int iEta = myAngleConverter->getGlobalEta(rawid, cluster.first);
      unsigned int hwNumber = myOmtfConfig->getLayerNumber(rawid);
      unsigned int iLayer = myOmtfConfig->getHwToLogicLayer().at(hwNumber);
      unsigned int iInput = getInputNumber(rawid, iProcessor, type);
      //      std::cout <<"ADDING HIT: iLayer = " << iLayer << " iInput: " << iInput << " iPhi: " << iPhi << std::endl;
      if (iLayer == 17 && (iInput == 0 || iInput == 1))
        continue;  // FIXME (MK) there is no RPC link for that input, because it is taken by DAQ link
      bool outres = result.addLayerHit(iLayer, iInput, iPhi, iEta);
      //      if (cSize>2) flag |= 2;
      //      if (!outres) flag |= 1;
      nClusters++;

      str << " RPC halfDigi "
          << " begin: " << cluster.first << " end: " << cluster.second << " iPhi: " << iPhi << " iEta: " << iEta
          << " hwNumber: " << hwNumber << " iInput: " << iInput << " iLayer: " << iLayer << " out: " << outres
          << std::endl;
    }
    //    if (nClusters > 2) flag=1;
  }

  edm::LogInfo("OMTFInputMaker") << str.str();
  return result;
}

setFlag()

void OMTFinputMaker::setFlag ( int  aFlag )

inline

Definition at line 36 of file OMTFinputMaker.h.

References flag.

Referenced by OMTFReconstruction::getProcessorCandidates().

{ flag = aFlag; }

Member Data Documentation

angleConverter

std::unique_ptr<OmtfAngleConverter> OMTFinputMaker::angleConverter

protected

Definition at line 138 of file OMTFinputMaker.h.

Referenced by initialize().

config

const OMTFConfiguration* OMTFinputMaker::config = nullptr

protected

Definition at line 137 of file OMTFinputMaker.h.

flag

int OMTFinputMaker::flag

private

Definition at line 75 of file OMTFinputMaker.h.

Referenced by getFlag(), and setFlag().

myAngleConverter

std::unique_ptr<AngleConverter> OMTFinputMaker::myAngleConverter

private

Definition at line 71 of file OMTFinputMaker.h.

Referenced by initialize(), processCSC(), processDT(), and processRPC().

myOmtfConfig

const OMTFConfiguration* OMTFinputMaker::myOmtfConfig

private

Definition at line 73 of file OMTFinputMaker.h.

Referenced by acceptDigi(), buildInputForProcessor(), getInputNumber(), initialize(), processCSC(), processDT(), and processRPC().