OMTFinputMaker Class Reference

#include <OMTFinputMaker.h>

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)
	Method translating trigger digis into input matrix with global phi coordinates. More...
void	initialize (const edm::EventSetup &es)
	OMTFinputMaker ()
	~OMTFinputMaker ()

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)
OMTFinput	processDT (const L1MuDTChambPhContainer *dtPhDigis, const L1MuDTChambThContainer *dtThDigis, unsigned int iProcessor, l1t::tftype type)
OMTFinput	processRPC (const RPCDigiCollection *rpcDigis, unsigned int iProcessor, l1t::tftype type)

Private Attributes
AngleConverter	myAngleConverter

Detailed Description

Definition at line 20 of file OMTFinputMaker.h.

Constructor & Destructor Documentation

OMTFinputMaker::OMTFinputMaker

(

)

Definition at line 21 of file OMTFinputMaker.cc.

{}

OMTFinputMaker::~OMTFinputMaker

(

)

Definition at line 30 of file OMTFinputMaker.cc.

{ }

Member Function Documentation

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 33 of file OMTFinputMaker.cc.

References OMTFConfiguration::barrelMax, OMTFConfiguration::barrelMin, l1t::bmtf, CSCDetId::chamber(), MuonSubdetId::CSC, DetId::det(), MuonSubdetId::DT, dt, l1t::emtf_neg, l1t::emtf_pos, CSCDetId::endcap(), OMTFConfiguration::endcap10DegMax, OMTFConfiguration::endcap10DegMin, OMTFConfiguration::endcap20DegMax, OMTFConfiguration::endcap20DegMin, OMTFConfiguration::instance(), RPCDetId::layer(), DetId::Muon, l1t::omtf_neg, l1t::omtf_pos, RPCDetId::region(), RPCDetId::ring(), CSCDetId::ring(), RPCDetId::roll(), MuonSubdetId::RPC, DTChamberId::sector(), RPCDetId::sector(), CSCDetId::station(), RPCDetId::station(), DetId::subdetId(), RPCDetId::subsector(), and DTChamberId::wheel().

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

                                {

  unsigned int aMin = OMTFConfiguration::instance()->barrelMin[iProcessor];
  unsigned int aMax = OMTFConfiguration::instance()->barrelMax[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 = OMTFConfiguration::instance()->endcap10DegMin[iProcessor];
      aMax = OMTFConfiguration::instance()->endcap10DegMax[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 = OMTFConfiguration::instance()->endcap10DegMin[iProcessor];
    aMax = OMTFConfiguration::instance()->endcap10DegMax[iProcessor];

    if( (type==l1t::tftype::emtf_pos || type==l1t::tftype::emtf_neg) &&
    csc.station()>1 && csc.ring()==1){
      aMin = OMTFConfiguration::instance()->endcap20DegMin[iProcessor];
      aMax = OMTFConfiguration::instance()->endcap20DegMax[iProcessor];
    }
    break;
  }    
  }
  
  if(aMax>aMin && aSector>=aMin && aSector<=aMax) return true;
  if(aMax<aMin && (aSector>=aMin || aSector<=aMax)) return true;

  return false;
}

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
	)

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

Definition at line 360 of file OMTFinputMaker.cc.

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

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

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

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 136 of file OMTFinputMaker.cc.

References OMTFConfiguration::barrelMin, l1t::bmtf, CSCDetId::chamber(), MuonSubdetId::CSC, DetId::det(), MuonSubdetId::DT, dt, l1t::emtf_neg, l1t::emtf_pos, OMTFConfiguration::endcap10DegMin, OMTFConfiguration::endcap20DegMin, OMTFConfiguration::instance(), RPCDetId::layer(), DetId::Muon, RPCDetId::region(), CSCDetId::ring(), RPCDetId::roll(), MuonSubdetId::RPC, DTChamberId::sector(), RPCDetId::sector(), CSCDetId::station(), RPCDetId::station(), DetId::subdetId(), and RPCDetId::subsector().

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

                                       {

  unsigned int iInput = 99;
  unsigned int aSector = 99;
  int aMin = OMTFConfiguration::instance()->barrelMin[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 = OMTFConfiguration::instance()->endcap10DegMin[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 = OMTFConfiguration::instance()->endcap10DegMin[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 = OMTFConfiguration::instance()->endcap20DegMin[iProcessor];
      if(iProcessor==5 && aSector<3) aMin = -2;
    }
    break;
  }
  }

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

void OMTFinputMaker::initialize

(

const edm::EventSetup &

es

)

Definition at line 24 of file OMTFinputMaker.cc.

References AngleConverter::checkAndUpdateGeometry(), and myAngleConverter.

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

{ 
  myAngleConverter.checkAndUpdateGeometry(es);
}

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

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 Trigger primitive quality. CSC central BX is 6 for some reason.

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 257 of file OMTFinputMaker.cc.

References funct::abs(), acceptDigi(), OMTFinput::addLayerHit(), CSCDetId, AngleConverter::getGlobalEta(), getInputNumber(), OMTFConfiguration::getLayerNumber(), AngleConverter::getProcessorPhi(), OMTFConfiguration::hwToLogicLayer, OMTFConfiguration::instance(), myAngleConverter, and mps_fire::result.

Referenced by buildInputForProcessor().

                          {

  OMTFinput result;
  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 (abs(digi->getBX()- 6)>0) continue;
      
      unsigned int hwNumber = OMTFConfiguration::instance()->getLayerNumber(rawid);
      if(OMTFConfiguration::instance()->hwToLogicLayer.find(hwNumber)==OMTFConfiguration::instance()->hwToLogicLayer.end()) continue;


      //unsigned int iLayer = OMTFConfiguration::instance()->hwToLogicLayer[hwNumber];      
      auto iter = OMTFConfiguration::instance()->hwToLogicLayer.find(hwNumber);
      unsigned int iLayer = iter->second;

      int iPhi = myAngleConverter.getProcessorPhi(iProcessor, type, CSCDetId(rawid), *digi);
      int iEta = myAngleConverter.getGlobalEta(rawid, *digi);
      if(abs(iEta)>1.26/2.61*240) continue;
      unsigned int iInput= getInputNumber(rawid, iProcessor, type);      
      result.addLayerHit(iLayer,iInput,iPhi,iEta);     
    }
  }      
  return result;
}

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

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 215 of file OMTFinputMaker.cc.

References acceptDigi(), OMTFinput::addLayerHit(), L1MuDTChambPhContainer::getContainer(), AngleConverter::getGlobalEta(), getInputNumber(), OMTFConfiguration::getLayerNumber(), AngleConverter::getProcessorPhi(), OMTFConfiguration::hwToLogicLayer, OMTFConfiguration::instance(), myAngleConverter, and mps_fire::result.

Referenced by buildInputForProcessor().

{

  OMTFinput result;
  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;

    unsigned int hwNumber = OMTFConfiguration::instance()->getLayerNumber(detid.rawId());
    if(OMTFConfiguration::instance()->hwToLogicLayer.find(hwNumber)==OMTFConfiguration::instance()->hwToLogicLayer.end()) continue;
    
    //unsigned int iLayer = OMTFConfiguration::instance()->hwToLogicLayer[hwNumber];
    auto iter = OMTFConfiguration::instance()->hwToLogicLayer.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);    
    result.addLayerHit(iLayer,iInput,iPhi,iEta);
    result.addLayerHit(iLayer+1,iInput,digiIt.phiB(),iEta);    
  }

  return result;
  
}

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

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==0 //FIXME: find a better place/way to filtering digi against quality/BX etc.

Definition at line 304 of file OMTFinputMaker.cc.

References acceptDigi(), OMTFinput::addLayerHit(), RPCDigi::bx(), HLT_FULL_cff::clusters, AngleConverter::getGlobalEta(), getInputNumber(), OMTFConfiguration::getLayerNumber(), AngleConverter::getProcessorPhi(), OMTFConfiguration::hwToLogicLayer, OMTFConfiguration::instance(), myAngleConverter, DetId::rawId(), mps_fire::result, reco::return(), and rpcPrimitiveCmp().

Referenced by buildInputForProcessor().

                               {

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

  const RPCDigiCollection & rpcDigiCollection = *rpcDigis;
  for (auto rollDigis : rpcDigiCollection) {
    RPCDetId roll = rollDigis.first;    
    unsigned int rawid = roll.rawId();
    if(!acceptDigi(rawid, iProcessor, type)) continue;    
    std::vector<RPCDigi> digisCopy;
    std::copy_if(rollDigis.second.first, rollDigis.second.second, std::back_inserter(digisCopy), [](const RPCDigi & aDigi){return (aDigi.bx()==0);});
    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 = (iPhiHalfStrip1+iPhiHalfStrip2)/2;
      int iEta =  myAngleConverter.getGlobalEta(rawid, cluster.first);      
      unsigned int hwNumber = OMTFConfiguration::instance()->getLayerNumber(rawid);
      //unsigned int iLayer = OMTFConfiguration::instance()->hwToLogicLayer[hwNumber];
      auto iter = OMTFConfiguration::instance()->hwToLogicLayer.find(hwNumber);
      unsigned int iLayer = iter->second;

      unsigned int iInput= getInputNumber(rawid, iProcessor, type);
      result.addLayerHit(iLayer,iInput,iPhi,iEta);

      str<<" RPC halfDigi "
           <<" begin: "<<cluster.first<<" end: "<<cluster.second
           <<" iPhi: "<<iPhi
           <<" iEta: "<<iEta
           <<" hwNumber: "<<hwNumber
           <<" iInput: "<<iInput
           <<" iLayer: "<<iLayer
           <<std::endl;      
    }
  }

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

Member Data Documentation

AngleConverter OMTFinputMaker::myAngleConverter

private

Definition at line 76 of file OMTFinputMaker.h.

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