RPCConeBuilder.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    RPCConeBuilder
// Class:      RPCConeBuilder
// 
//
// Original Author:  Tomasz Maciej Frueboes
//         Created:  Fri Feb 22 13:57:06 CET 2008
//
//


// system include files

// user include files

#include "L1Trigger/RPCTrigger/interface/RPCConeBuilder.h"
#include "L1Trigger/RPCTrigger/interface/RPCStripsRing.h"

//#include "L1TriggerConfig/RPCConeBuilder/interface/RPCConeBuilder.h"
//#include "L1TriggerConfig/RPCConeBuilder/interface/RPCStripsRing.h"

#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "FWCore/Framework/interface/ModuleFactory.h"


#include "CondFormats/DataRecord/interface/L1RPCConeDefinitionRcd.h"

#include <sstream>
#include <vector>

#include "DataFormats/MuonDetId/interface/RPCDetId.h"


RPCConeBuilder::RPCConeBuilder(const edm::ParameterSet& iConfig) :
      m_towerBeg(iConfig.getParameter<int>("towerBeg")),
      m_towerEnd(iConfig.getParameter<int>("towerEnd")),
      m_runOnceBuildCones(false)
      //m_rollBeg(iConfig.getParameter<int>("rollBeg")),
      //m_rollEnd(iConfig.getParameter<int>("rollEnd")),
      //m_hwPlaneBeg(iConfig.getParameter<int>("hwPlaneBeg")),
      //m_hwPlaneEnd(iConfig.getParameter<int>("hwPlaneEnd"))
{
   
  setWhatProduced(this, (dependsOn (&RPCConeBuilder::geometryCallback) &  
                                   (&RPCConeBuilder::coneDefCallback)     )
                        );

  /* TT
   for (int i = m_towerBeg; i <= m_towerEnd; ++i){
      
      std::stringstream name;
      name << "lpSizeTower" << i;
      
      L1RPCConeBuilder::TLogPlaneSize newSizes = 
            iConfig.getParameter<std::vector<int> >(name.str().c_str());
      
      m_LPSizesInTowers.push_back(newSizes);
      
   }
  
  */
   
   //  hw planes numbered from 0 to 5
   // rolls from 0 to 17 (etaPartition)
   //
   //  rollConnLP_[roll]_[hwPlane-1]
   //  rollConnLP_5_3 = cms.vint32(6, 0, 0),
   //     ----- roll 5, hwPlane 4 (3+1) is logplane 6 (OK)
   //
   //  rollConnT_[roll]_[hwPlane-1]
   //  rollConnT_5_3 = cms.vint32(4, -1, -1),
   //     ----- roll 5, hwPlane 4 (3+1) contirubtes to tower 4 (OK)
  
  /*
   for (int roll = m_rollBeg; roll <= m_rollEnd; ++roll){
      L1RPCConeDefinition::THWplaneToTower newHwPlToTower;
      L1RPCConeDefinition::THWplaneToLP newHWplaneToLP;
      for (int hwpl = m_hwPlaneBeg; hwpl <= m_hwPlaneEnd; ++hwpl){
         std::stringstream name;
         name << "rollConnLP_" << roll << "_" << hwpl;
         
         L1RPCConeDefinition::TTowerList newListLP = 
               iConfig.getParameter<std::vector<int> >(name.str().c_str());
         newHWplaneToLP.push_back(newListLP);
         
         
         std::stringstream name1;
         name1 << "rollConnT_" << roll << "_" << hwpl;
         
         L1RPCConeDefinition::TLPList newListT = 
               iConfig.getParameter<std::vector<int> >(name1.str().c_str());
         newHwPlToTower.push_back(newListT);
      }
      m_RingsToTowers.push_back(newHwPlToTower);
      m_RingsToLP.push_back(newHWplaneToLP);
   }
  */
}




// 
// member functions
//

// ------------ method called to produce the data  ------------
RPCConeBuilder::ReturnType
RPCConeBuilder::produce(const L1RPCConeBuilderRcd& iRecord)
//RPCConeBuilder::produce(const L1RPCConfigRcd& iRecord)
{

  
   //std::cout << " RPCConeBuilder::produce called " << std::endl;
   using namespace edm::es;
   auto pL1RPCConeBuilder = std::make_shared<L1RPCConeBuilder>();
   
   pL1RPCConeBuilder->setFirstTower(m_towerBeg);
   pL1RPCConeBuilder->setLastTower(m_towerEnd);
   
   /*
   pL1RPCConeBuilder->setLPSizeForTowers(m_LPSizesInTowers);
   pL1RPCConeBuilder->setRingsToLP(m_RingsToLP);
   */
   
   //iRecord.get(m_rpcGeometry);
   //iRecord.get(m_L1RPCConeDefinition);
   buildCones(m_rpcGeometry);

   // Compress all connections. Since members of this class are shared
   // pointers this call will compress all data
   m_ringsMap.begin()->second.compressConnections();
      
   pL1RPCConeBuilder->setConeConnectionMap(m_ringsMap.begin()->second.getConnectionsMap());
   
   pL1RPCConeBuilder->setCompressedConeConnectionMap(
           m_ringsMap.begin()->second.getCompressedConnectionsMap());
           
   m_ringsMap.clear(); // free mem
      
   return pL1RPCConeBuilder;
   
}

// ----------------------------------------------------------
void RPCConeBuilder::geometryCallback( const MuonGeometryRecord& record ){

  //std::cout << " Geometry callback called " << std::endl; 
  m_runOnceBuildCones = false; // allow re-running of buildCones
  record.get(m_rpcGeometry);
  
  
}

void RPCConeBuilder::coneDefCallback( const L1RPCConeDefinitionRcd& record ){

  //std::cout << " ConeDef callback called " << std::endl; 
  m_runOnceBuildCones = false; // allow re-running of buildCones

  //edm::ESHandle<RPCGeometry> rpcGeom;
  record.get(m_L1RPCConeDefinition);
  
  //std::cout << " ConeDef callback exit " << std::endl; 
  //std::cout.flush();
  //buildCones(rpcGeom);
  
}



void RPCConeBuilder::buildCones(const edm::ESHandle<RPCGeometry> & rpcGeom ){
  

  if (!m_runOnceBuildCones){
    m_runOnceBuildCones = true;
  } else {
    throw cms::Exception("RPCInternal") << "buildCones called twice \n";
  }
  
  //std::cout << "    ---> buildCones called " << std::endl; 
  
  // fetch geometricall data
  auto uncompressedCons = std::make_shared<L1RPCConeBuilder::TConMap>();
  
  
  int rolls = 0;
  for(TrackingGeometry::DetContainer::const_iterator it = rpcGeom->dets().begin();
      it != rpcGeom->dets().end();
      ++it)
  {
  
      if( dynamic_cast< RPCRoll const * >( *it ) == 0 ) continue;
      
      ++rolls;
      RPCRoll const* roll = dynamic_cast< RPCRoll const*>( *it );
      
      int ringId = RPCStripsRing::getRingId(roll);
      if ( m_ringsMap.find(ringId)  == m_ringsMap.end() ) {
        m_ringsMap[ringId]=RPCStripsRing(roll, uncompressedCons);
      } else {
         m_ringsMap[ringId].addRoll(roll);
      }
      //RPCStripsRing sr(roll);
      
  }

  //std::cout << " found: " << rolls << " dets" << std::endl;
  
  RPCStripsRing::TIdToRindMap::iterator it = m_ringsMap.begin();
  
  // filtermixed strips, fill gaps with virtual strips
  for (;it != m_ringsMap.end(); ++it){
     
    //int before = it->second.size();
    
    it->second.filterOverlapingChambers();    
    it->second.fillWithVirtualStrips();

    
    //std::cout << it->first << " " <<  it->second.isReferenceRing() << " " 
     //   << before << " -> " 
      //  << it->second.size() << std::endl;
    
    // In reference plane we should have 144*8 = 1152 strips
    //int plane = it->first/1000;
    int etaPart =  it->first%100; 
    if (it->second.isReferenceRing() && (it->second.size() != 1152)){

      if (std::abs(etaPart)>=14 || std::abs(etaPart)<=17 ) {
        //std::cout << "etaPart "  << etaPart << " size " << it->second.size() << std::endl;
      }
      else {
        throw cms::Exception("RPCInternal") << "Problem: refring " << it->first 
          << " has " << it->second.size() << " strips \n";
      }
    }
    
    
  }
    
  // Xcheck, if rings are symettrical 
  it = m_ringsMap.begin();
  for (;it != m_ringsMap.end(); ++it){
    int key = it->first;
    int sign = key/100 - (key/1000)*10;
    
    if (sign == 0) {
      key += 100;
    } else {
      key -= 100;
    }
    
    if  (key != 2000){// Hey 2100 has no counterring
      if (it->second.size() !=  m_ringsMap[key].size())  
      {
         throw cms::Exception("RPCInternal") << " Size differs for ring " << key << " +- 100 \n";
      }
    }
  
  
  }
      
  buildConnections();
}


void RPCConeBuilder::buildConnections(){



   RPCStripsRing::TIdToRindMap::iterator itRef = m_ringsMap.begin();
   for (;itRef != m_ringsMap.end(); ++itRef){ // iterate over reference rings
      
      
      RPCStripsRing::TOtherConnStructVec ringsToConnect; 
      
      if (!itRef->second.isReferenceRing()) continue; // iterate over reference rings
      
      RPCStripsRing::TIdToRindMap::iterator itOther = m_ringsMap.begin();
      for (;itOther != m_ringsMap.end(); ++itOther){ // iterate over nonreference rings
         
         if (itOther->second.isReferenceRing()) continue; // iterate over nonreference rings
         
         std::pair<int,int> pr = areConnected(itRef, itOther);
         if ( pr.first != -1 ) {
            RPCStripsRing::TOtherConnStruct newOtherConn;
            newOtherConn.m_it = itOther;
            newOtherConn.m_logplane = pr.first;
            newOtherConn.m_logplaneSize = pr.second;
            ringsToConnect.push_back(newOtherConn);
         }
         
         
      } // OtherRings iteration ends
      
      // 
      std::pair<int,int> prRef = areConnected(itRef, itRef);
      if (prRef.first == -1){
        throw cms::Exception("RPCConfig") << " Cannot determine logplane for reference ring "
            << itRef->first << "\n ";
      }

      /*&
      if (prRef.second != 8){
        // XXX        
        throw cms::Exception("RPCConfig") << " logplaneSize for reference ring "
            << itRef->first << " wrong "
            << " logplane: " << prRef.first
            << " etaPart: " << itRef->second.getEtaPartition()
            << " tower: " << itRef->second.getTowerForRefRing()
            << " hwPlane: " << itRef->second.getHwPlane()
            << " strips " << prRef.second << "\n";
      }*/
      
      itRef->second.createRefConnections(ringsToConnect, prRef.first, prRef.second);
      
   } // RefRings iteration ends

   
   // Fetch connection data, and save in one place
   /*
   RPCStripsRing::TIdToRindMap::iterator it = m_ringsMap.begin();
   for (;it != m_ringsMap.end(); ++it) { 
     
     L1RPCConeBuilder::TConMap nmap = it->second.getConnectionsMap();
     L1RPCConeBuilder::TConMap::iterator newMapIt = nmap.begin();
     for (; newMapIt != nmap.end(); ++ newMapIt) {
       uint32_t raw = newMapIt->first;
       TStrip2ConVec stripsVec = newMapIt->second;
       TStrip2ConVec::iterator stripIt = stripsVec.first
       //unsigned char strip = 
     
     
     }
     
     
  }*/
   
}

// first - logplane
// second - logplanesize
std::pair<int, int> RPCConeBuilder::areConnected(RPCStripsRing::TIdToRindMap::iterator ref,
                                  RPCStripsRing::TIdToRindMap::iterator other){

  int logplane = -1;
  
  //std::cout << "Checking " << ref->first << " and " << other->first << std::endl;
  
  // Do not connect  rolls lying on the oposite side of detector
  if ( ref->second.getEtaPartition()*other->second.getEtaPartition()<0  )
    return std::make_pair(-1,0);  
  
  
  /*std::cout << "Ref " << ref->second.getEtaPartition() << " " <<ref->second.getHwPlane() << std::endl;
  std::cout << "Other " << other->second.getEtaPartition() << " " <<other->second.getHwPlane() << std::endl;
  std::cout.flush();*/
  
  // refRing and otherRing areConnected, if they contribute to the same tower
  /*
  L1RPCConeDefinition::TTowerList refTowList 
      = m_L1RPCConeDefinition->getRingsToTowers().at(std::abs(ref->second.getEtaPartition()))
                         .at(ref->second.getHwPlane()-1);
      

  L1RPCConeDefinition::TTowerList otherTowList 
      = m_L1RPCConeDefinition->getRingsToTowers().at(std::abs(other->second.getEtaPartition()))
                         .at(other->second.getHwPlane()-1);
  */

  L1RPCConeDefinition::TRingToTowerVec::const_iterator itRef
      = m_L1RPCConeDefinition->getRingToTowerVec().begin();
  
  const L1RPCConeDefinition::TRingToTowerVec::const_iterator itEnd 
      = m_L1RPCConeDefinition->getRingToTowerVec().end();
    
  L1RPCConeDefinition::TRingToTowerVec::const_iterator itOther = itRef;
  
  int refTowerCnt = 0;
  int index = -1;
  int refTower = -1;
  
  for (;itRef != itEnd; ++itRef){
    if ( itRef->m_etaPart != std::abs(ref->second.getEtaPartition())
        || itRef->m_hwPlane != std::abs(ref->second.getHwPlane()-1) // -1?
       ) continue;
      
    ++refTowerCnt;
    refTower = itRef->m_tower;
    
    for (;itOther != itEnd; ++itOther){
      if ( itOther->m_etaPart != std::abs(other->second.getEtaPartition())
        || itOther->m_hwPlane != std::abs(other->second.getHwPlane()-1) // -1?
        ) continue;  
      
      if (itOther->m_tower == refTower) index = itOther->m_index;
      
    }
    
  }
  
  if(refTowerCnt>1){
    throw cms::Exception("RPCConeBuilder") << " Reference(?) ring "
        << ref->first << " "
        << "wants to be connected to " << refTowerCnt << " towers \n";
  
  }

  if(refTowerCnt==0){
    throw cms::Exception("RPCConeBuilder") << " Reference(?) ring "
        << ref->first << " "
        << " is not connected anywhere \n";
  
  }
  
  /*  
  if(index == -1){
    throw cms::Exception("RPCConeBuilder") << "Wrong Index -1 \n"
                                  }*/
  
  
  /*
  int refTower = -1;
  
  L1RPCConeDefinition::TTowerList::iterator rtlIt = refTowList.begin();
  for (; rtlIt != refTowList.end(); ++rtlIt){
  
     if ( *rtlIt >= 0 && refTower < 0){
        refTower = *rtlIt;
     }
     else if ( *rtlIt >= 0 && refTower >= 0) {
      throw cms::Exception("RPCConfig") << " Reference(?) ring "
            << ref->first << " "
            << "wants to be connected more than one tower: "
            << refTower << " "
            << *rtlIt << "\n";
     
     }
  
  }
  
  if (refTower < 0) {
     throw cms::Exception("RPCConfig") << " Reference(?) ring "
           << ref->first
           << " is not connected anywhere \n";
  }
  
  L1RPCConeDefinition::TTowerList::iterator otlIt = otherTowList.begin();
  
  int index = -1, i = 0;
  for (; otlIt != otherTowList.end(); ++otlIt){
     if (*otlIt == refTower) {
        index = i;
     }
     ++i;
  }
  */
  
  int lpSize = 0;
  if (index != -1){
    /*
    logplane = m_L1RPCConeDefinition->getRingsToLP().at(std::abs(other->second.getEtaPartition()))
           .at(other->second.getHwPlane()-1)
    .at(index);*/
    {
      L1RPCConeDefinition::TRingToLPVec::const_iterator it = m_L1RPCConeDefinition->getRingToLPVec().begin();
      L1RPCConeDefinition::TRingToLPVec::const_iterator itEnd = m_L1RPCConeDefinition->getRingToLPVec().end();
      for (;it!=itEnd;++it){
        
        if (it->m_etaPart != std::abs(other->second.getEtaPartition())
            || it->m_hwPlane != std::abs(other->second.getHwPlane()-1) 
            || it->m_index != index) continue;
        
        logplane = it->m_LP;  
        
      }
    }    
    //lpSize = m_L1RPCConeDefinition->getLPSizeForTowers().at(refTower).at(logplane-1);
    
    {
      L1RPCConeDefinition::TLPSizeVec::const_iterator it = m_L1RPCConeDefinition->getLPSizeVec().begin();
      L1RPCConeDefinition::TLPSizeVec::const_iterator itEnd = m_L1RPCConeDefinition->getLPSizeVec().end();
      for (;it!=itEnd;++it){
              
        //std::cout << it->m_LP  << " " << logplane << std::endl;
        if (it->m_tower != std::abs(refTower) || it->m_LP != logplane-1) continue;
        lpSize = it->m_size;
              
      }
  
              //FIXME
      if (lpSize==-1) {
                //throw cms::Exception("getLogStrip") << " lpSize==-1\n";
      }
    }
  }
  
  
  /*
  if (logplane != -1){
    
     std::cout << ref->first << " <-> " << other->first 
           << " logplane " << logplane
           << " lpsize " << lpSize 
           << std::endl;
  }//*/
  
  return std::make_pair(logplane,lpSize);

}