OMTFinputMaker Class Reference

#include <OMTFinputMaker.h>

Public Member Functions
const 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)
void	processCSC (const CSCCorrelatedLCTDigiCollection *cscDigis, unsigned int iProcessor, l1t::tftype type)
void	processDT (const L1MuDTChambPhContainer *dtPhDigis, const L1MuDTChambThContainer *dtThDigis, unsigned int iProcessor, l1t::tftype type)
void	processRPC (const RPCDigiCollection *rpcDigis, unsigned int iProcessor, l1t::tftype type)

Private Attributes
std::unique_ptr< AngleConverter >	katownik
OMTFinput *	myInput
	Output object. More...

Detailed Description

Definition at line 20 of file OMTFinputMaker.h.

Constructor & Destructor Documentation

OMTFinputMaker::OMTFinputMaker

(

)

Definition at line 21 of file OMTFinputMaker.cc.

References katownik, and myInput.

                              { 

  myInput = new OMTFinput();

  katownik.reset(new AngleConverter());
}

OMTFinputMaker::~OMTFinputMaker

(

)

Definition at line 35 of file OMTFinputMaker.cc.

References myInput.

                               { 

  if(myInput) delete myInput;

}

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 42 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, 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::barrelMin[iProcessor];
  unsigned int aMax = OMTFConfiguration::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::endcap10DegMin[iProcessor];
      aMax = OMTFConfiguration::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::endcap10DegMin[iProcessor];
    aMax = OMTFConfiguration::endcap10DegMax[iProcessor];

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

  return false;
}

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

References OMTFinput::clear(), myInput, processCSC(), processDT(), and processRPC().

Referenced by L1TMuonOverlapTrackProducer::produce().

                                            {  
  myInput->clear(); 

  processDT(dtPhDigis, dtThDigis, iProcessor, type);
  processCSC(cscDigis, iProcessor, type);
  processRPC(rpcDigis, iProcessor, type);

  return myInput;

}

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

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 alsgo 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 145 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, 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::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 = 0;
      //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;

      if(type==l1t::tftype::bmtf)iRoll = 1;
    }
    if(rpc.region()!=0){
      aSector = (rpc.sector()-1)*6+rpc.subsector();
      aMin = OMTFConfiguration::endcap10DegMin[iProcessor];
      if(iProcessor==5 && aSector<5) aMin = -3;
    }    
    break;
  }
  case MuonSubdetId::DT: {
    DTChamberId dt(rawId);
    aSector = dt.sector();
    if(iProcessor==5 && aSector<3) aMin = 0;
    break;
  }
  case MuonSubdetId::CSC: {   
    CSCDetId csc(rawId);    
    aSector = csc.chamber();    
    aMin = OMTFConfiguration::endcap10DegMin[iProcessor];       
    if(iProcessor==5 && aSector<5) aMin = -3;
    if( (type==l1t::tftype::emtf_pos || type==l1t::tftype::emtf_neg) &&
    csc.station()>1 && csc.ring()==1){
      aMin = OMTFConfiguration::endcap20DegMin[iProcessor];
      if(iProcessor==5 && aSector<3) aMin = -1;
    }
    break;
  }
  }

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

  return iInput;
}

void OMTFinputMaker::initialize

(

const edm::EventSetup &

es

)

Definition at line 29 of file OMTFinputMaker.cc.

References katownik.

Referenced by L1TMuonOverlapTrackProducer::produce().

                                                      { 

  katownik->checkAndUpdateGeometry(es);
}

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

References funct::abs(), acceptDigi(), OMTFinput::addLayerHit(), end, spr::find(), getInputNumber(), OMTFConfiguration::getLayerNumber(), OMTFConfiguration::hwToLogicLayer, katownik, and myInput.

Referenced by buildInputForProcessor().

                          {

  if(!cscDigis) return;

  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::getLayerNumber(rawid);
      if(OMTFConfiguration::hwToLogicLayer.find(hwNumber)==OMTFConfiguration::hwToLogicLayer.end()) continue;

      unsigned int iLayer = OMTFConfiguration::hwToLogicLayer[hwNumber];      
      int iPhi = katownik->getGlobalPhi(rawid, *digi);
      int iEta = katownik->getGlobalEta(rawid, *digi);
      if(abs(iEta)>1.26/2.61*240) continue;
      unsigned int iInput= getInputNumber(rawid, iProcessor, type);      
      myInput->addLayerHit(iLayer,iInput,iPhi,iEta);     
    }
  }      
}

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

References acceptDigi(), OMTFinput::addLayerHit(), spr::find(), L1MuDTChambPhContainer::getContainer(), getInputNumber(), OMTFConfiguration::getLayerNumber(), OMTFConfiguration::hwToLogicLayer, katownik, and myInput.

Referenced by buildInputForProcessor().

                          {

  if(!dtPhDigis) return;

  for (const auto digiIt: *dtPhDigis->getContainer()) {

    DTChamberId detid(digiIt.whNum(),digiIt.stNum(),digiIt.scNum()+1);

    if(!acceptDigi(detid.rawId(), iProcessor, type)) continue;
    if (digiIt.bxNum()!= 0 || digiIt.BxCnt()!= 0 || digiIt.Ts2Tag()!= 0 || digiIt.code()<4) continue;

    unsigned int hwNumber = OMTFConfiguration::getLayerNumber(detid.rawId());
    if(OMTFConfiguration::hwToLogicLayer.find(hwNumber)==OMTFConfiguration::hwToLogicLayer.end()) continue;
    
    unsigned int iLayer = OMTFConfiguration::hwToLogicLayer[hwNumber];   
    int iPhi =  katownik->getGlobalPhi(detid.rawId(), digiIt);
    int iEta =  katownik->getGlobalEta(detid.rawId(), digiIt, dtThDigis);
    unsigned int iInput= getInputNumber(detid.rawId(), iProcessor, type);

    myInput->addLayerHit(iLayer,iInput,iPhi,iEta);
    myInput->addLayerHit(iLayer+1,iInput,digiIt.phiB(),iEta);
  }
  
}

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

private

Take the CSC digis, select chambers connected to given processor, convers logal angles to global scale. Decluster nearby hits in single chamber, by taking average clister position, expressed in half RPC strip: pos = cluster_begin + cluster_end)/2

Check it the data fits into given processor input range

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

If phi1 is close to Pi, and phi2 close to -Pi the results phi is 0 instead -pi

Definition at line 299 of file OMTFinputMaker.cc.

References acceptDigi(), OMTFinput::addLayerHit(), RPCDigi::bx(), getInputNumber(), OMTFConfiguration::getLayerNumber(), OMTFConfiguration::hwToLogicLayer, katownik, M_PI, myInput, OMTFConfiguration::nPhiBins, phi, query::result, reco::return(), and rpcPrimitiveCmp().

Referenced by buildInputForProcessor().

                               {

  if(!rpcDigis) return;

  std::ostringstream myStr;

  typedef std::pair<RPCDigi *, RPCDigi *> halfDigi;
  
  auto chamber = rpcDigis->begin();
  auto chend  = rpcDigis->end();
  for( ; chamber != chend; ++chamber ) {
    unsigned int rawid = (*chamber).first;
    
    if(!acceptDigi(rawid, iProcessor, type)) continue;

    std::vector<RPCDigi> digisCopy;
    std::copy_if((*chamber).second.first, (*chamber).second.second, std::back_inserter(digisCopy), [](const RPCDigi & aDigi){return (aDigi.bx()==0);});
    std::sort(digisCopy.begin(),digisCopy.end(),rpcPrimitiveCmp);
    std::vector<halfDigi> result;
    for(auto &stripIt: digisCopy) {
      if(result.empty()) result.push_back(halfDigi(&stripIt,&stripIt));
      else if (stripIt.strip() - result.back().second->strip() == 1) result.back().second = &stripIt;
      else if (stripIt.strip() - result.back().second->strip() > 1) result.push_back(halfDigi(&stripIt,&stripIt));
    }
      for(auto halfDigiIt:result){
    /* This code should be used when LUT for RPC angle converiosn will be implemented.
    int strip1 = halfDigiIt.first->strip();
    int strip2 = halfDigiIt.second->strip();
    int clusterHalfStrip = strip1 + strip2;
    int iPhi = katownik->getGlobalPhi(rawid,clusterHalfStrip);
    */
    float phi1 =  katownik->getGlobalPhi(rawid,*halfDigiIt.first);
    float phi2 =  katownik->getGlobalPhi(rawid,*halfDigiIt.second);
    float phi = (phi1+phi2)/2.0;
    if(phi1*phi2<0 && fabs(phi1)>M_PI/2.0) phi = (M_PI-phi)*(1 - 2*std::signbit(phi));
    int iPhi =  phi/(2.0*M_PI)*OMTFConfiguration::nPhiBins;
    int iEta =  katownik->getGlobalEta(rawid,*halfDigiIt.first);
    unsigned int hwNumber = OMTFConfiguration::getLayerNumber(rawid);
    unsigned int iLayer = OMTFConfiguration::hwToLogicLayer[hwNumber];
    unsigned int iInput= getInputNumber(rawid, iProcessor, type);
    
    myInput->addLayerHit(iLayer,iInput,iPhi,iEta);
    
    myStr<<" RPC halfDigi "
         <<" begin: "<<halfDigiIt.first->strip()<<" end: "<<halfDigiIt.second->strip()
         <<" iPhi: "<<iPhi
         <<" iEta: "<<iEta
         <<" hwNumber: "<<hwNumber
         <<" iInput: "<<iInput
         <<" iLayer: "<<iLayer
         <<std::endl;
      }    
  }
  edm::LogInfo("OMTFInputMaker")<<myStr.str();
}

Member Data Documentation

std::unique_ptr<AngleConverter> OMTFinputMaker::katownik

private

Definition at line 81 of file OMTFinputMaker.h.

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

OMTFinput* OMTFinputMaker::myInput

private

Output object.

Definition at line 79 of file OMTFinputMaker.h.

Referenced by buildInputForProcessor(), OMTFinputMaker(), processCSC(), processDT(), processRPC(), and ~OMTFinputMaker().