CmsTrackerDetIdBuilder.cc

Go to the documentation of this file.

#include "Geometry/TrackerNumberingBuilder/plugins/CmsTrackerDetIdBuilder.h"
#include "Geometry/TrackerNumberingBuilder/interface/GeometricDet.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "DataFormats/DetId/interface/DetId.h"

#include <iostream>
#include <fstream>
#include <sstream>
#include <string>
#include <bitset>

CmsTrackerDetIdBuilder::CmsTrackerDetIdBuilder( unsigned int layerNumberPXB )
  : m_layerNumberPXB( layerNumberPXB )
{}

GeometricDet*
CmsTrackerDetIdBuilder::buildId( GeometricDet* tracker )
{
  DetId t( DetId::Tracker, 0 );
  tracker->setGeographicalID( t );
  iterate( tracker, 0, tracker->geographicalID().rawId() );

  return tracker;
}

void
CmsTrackerDetIdBuilder::iterate( GeometricDet *in, int level, unsigned int ID )
{
  std::bitset<32> binary_ID(ID);

  // SubDetector (useful to know fron now on, valid only after level 0, where SubDetector is assigned)
  uint32_t mask = (7<<25);
  uint32_t iSubDet = ID & mask;
  iSubDet = iSubDet >> 25;
  //
  
  switch( level )
  {
    // level 0
  case 0:
    {  
      for( uint32_t i = 0; i<(in)->components().size(); i++ )
      {
    GeometricDet* component = in->component(i);
    uint32_t jSubDet = component->geographicalID().rawId();
    uint32_t temp = ID;
    temp |= (jSubDet<<25);
    component->setGeographicalID(temp); 
    
    switch( jSubDet )
    { 
      // PXF
    case 2:
      {
        // SubDetector Side start bit is 23 [3 unused and Side length is 2 bit]
        if(component->translation().z()<0. )
        {
          temp |= (1<<23); // PXF-
        }
        else
        {
          temp |= (2<<23); // PXF+
        }
        break;
      }
      // TID
    case 4:
      {
        temp|= (0<<15); // SubDetector Side start bit is 13 [10 unused and Side length is 2 bit]
        if((component)->components()[0]->translation().z()<0. )
        {
          temp |= (1<<13); // TIDB = TID-
        }
        else
        {
          temp |= (2<<13); // TIDF = TID+
        }
        break;
      }
      // TEC
    case 6:
      {
        temp |= (0<<20); // SubDetector Side start bit is 18 [5 unused and Side length is 2 bit]
        if(component->translation().z()<0. )
        {
          temp |= (1<<18); // TEC-
        }
        else
        {
          temp |= (2<<18); // TEC+
        }
        break;
      }
      // PXB, TIB, TOB (barrel)
    default:
      {
        // do nothing
      }
      // SubDetector switch ends
    }
    
    component->setGeographicalID(DetId(temp));  
    
    // next level
    iterate(component,level+1,component->geographicalID().rawId());
      }
      break;
    }
  // level 1
  case 1:
    {
      for( uint32_t i = 0; i < (in)->components().size(); i++ )
      {
    auto component = in->component(i);
    uint32_t temp = ID;
      
    switch( iSubDet )
    {
    // PXB
    case 1:
      {
        temp |= (component->geographicalID().rawId() << m_layerNumberPXB ); // Layer Number start bit is 16 [5 unused]
        break;
      }
    // PXF
    case 2:
      {
        temp |= (component->geographicalID().rawId() << 16 ); // Disk Number start bit is 16
        break;
      }
    // TIB
    case 3:
      {
        temp |= (component->geographicalID().rawId() << 14); // Layer Number start bit is 14 [8 unused]
        break;
      }
    
    // TID
    case 4:
      {
        temp |= (component->geographicalID().rawId() << 11); // Disk (Wheel) Number start bit is 11
        break;
      }
    
    // TOB
    case 5:
      {
        temp |= (component->geographicalID().rawId() << 14); // Layer Number start bit is 14 [8 unused]
        break;
      }
    
    // TEC
    case 6:
      {
        temp |= (component->geographicalID().rawId() << 14); // Wheel Number start bit is 14
        break;
      }
      
    // the rest
    default:
      {
        // do nothing
      }
    
    // SubDetector switch ends
    }
      
    component->setGeographicalID( temp );
      
    // next level
    iterate(component,level+1,component->geographicalID().rawId());      
      }
    
      break; 
    }    
    // level 2
  case 2: {
    
    for (uint32_t i=0;i<(in)->components().size();i++) {
      auto component = in->component(i);

      switch (iSubDet) {

    // PXB
      case 1:
    {
      uint32_t temp = ID;
      // Ladder Starting bit = 2 (last unused) + 6 (Module Number) = 8
      temp |= (component->geographicalID().rawId()<<8);
      component->setGeographicalID(temp);
      break;
    }
    
    // PXF
      case 2:
    {
      uint32_t temp = ID;
      // Blade Starting bit = 1 (last unused) + 5 (Module Number) + 2 (Plaquette part) = 8
      temp |= (component->geographicalID().rawId()<<8);
      component->setGeographicalID(temp);
      break;
    }
    
    // TIB
      case 3:
    {
      uint32_t temp = ID;
      // Side+Part+String Starting bit = 2 (Module Type) + 2 (Module Number) = 4
      temp |= (component->geographicalID().rawId()<<4);
      component->setGeographicalID(temp);
      break;
    }
    
    // TID
      case 4:
    {
      uint32_t temp = ID;
      // Ring+Part Starting bit = 2 (Module Type) + 5 (Module Number) + 2 (Disk Part)= 9
      temp |= (component->geographicalID().rawId()<<9);
      component->setGeographicalID(DetId(temp));
      break;
    }
    
    // TOB
      case 5:
    {
      uint32_t temp = ID;
      // Side+Rod Starting bit = 2 (Module Type) + 3 (Module Number) = 5
      temp |= (component->geographicalID().rawId()<<5);
      component->setGeographicalID(temp);
      break;
    }
    
    // TEC
      case 6:
    {
      uint32_t temp = ID;
      // Petal+Part Starting bit = 2 (Module Type) + 3 (Module Number) + 3 (Ring Number) = 8
      temp |= (component->geographicalID().rawId()<<8);
      component->setGeographicalID(temp);
      break;
    }
    
    // the rest
      default:
    {
      // do nothing
    }
    
    // SubDetector switch ends
      }

      // next level
      iterate(component,level+1,component->geographicalID().rawId());
    }
    
    break;    
  }
    
    // level 3
  case 3:
    {
      for (uint32_t i=0;i<(in)->components().size();i++) {
    auto component = in->component(i);

    switch (iSubDet) {
      
      // TEC
    case 6:
      {
        // Ring Starting bit = 2 (Module Type) + 3 (Module Number)
        uint32_t temp = ID;
        temp |= (component->geographicalID().rawId()<<5);
        component->setGeographicalID(temp);
        break;
      }
      
      // the others but TEC
    default:
      {
        uint32_t temp = ID;
        temp |= (component->geographicalID().rawId()<<2); // Starting bit = 2 (Module Type)
        component->setGeographicalID(temp);
      }
      
      // SubDetector switch ends
    }
    
    // next level
    iterate(component,level+1,component->geographicalID().rawId());
    
      }
      
      break;
    }
    
    // level 4
  case 4:
    {

      for (uint32_t i=0;i<(in)->components().size();i++) {
    
    switch (iSubDet) {
      
      // TEC
    case 6:
      {
        // Module Number bit = 2 (Module Type)
        uint32_t temp = ID;
        auto component = in->component(i);
        temp |= (component->geographicalID().rawId()<<2);
        component->setGeographicalID(temp);
        // next level
        iterate(component,level+1,component->geographicalID().rawId());
        break;
      }
      
      // the others stop here!
    default:
      {
        for (uint32_t i=0;i<(in)->components().size();i++) {
          uint32_t temp = ID;
          auto component = in->component(i);
          temp |= (component->geographicalID().rawId());
          component->setGeographicalID(temp);
        }
      }
      
      // SubDetector switch ends
    }

      }
      
      break;
    }
    
    // level 5
  case 5:
    {
      // TEC Module Type (only TEC arrives here)
      for (uint32_t i=0;i<(in)->components().size();i++) {
    uint32_t temp = ID;
    auto component = in->component(i);
    temp |= (component->geographicalID().rawId());
    component->setGeographicalID(temp);
      }
      break;
    }
    
    // throw exception
  default:
    {
      cms::Exception("LogicError") <<" CmsTrackerDetIdBuilder invalid level "<< level;
    }
    
    // level switch ends
  }
  
  return;

}