CmsTrackerLayerBuilder Class Reference

#include <CmsTrackerLayerBuilder.h>

Inheritance diagram for CmsTrackerLayerBuilder:

Private Member Functions
virtual void	buildComponent (DDFilteredView &, GeometricDet *, std::string)
virtual void	sortNS (DDFilteredView &, GeometricDet *)

Additional Inherited Members
Public Member Functions inherited from CmsTrackerLevelBuilder
virtual void	build (DDFilteredView &, GeometricDet *, std::string)
virtual	~CmsTrackerLevelBuilder ()
Protected Attributes inherited from CmsTrackerLevelBuilder
CmsTrackerStringToEnum	theCmsTrackerStringToEnum

Detailed Description

Class which contructs TIB/TOB layers

Definition at line 10 of file CmsTrackerLayerBuilder.h.

Member Function Documentation

void CmsTrackerLayerBuilder::buildComponent ( DDFilteredView &  fv, GeometricDet *  g, std::string  s  )

privatevirtual

Implements CmsTrackerLevelBuilder.

Definition at line 15 of file CmsTrackerLayerBuilder.cc.

References GeometricDet::addComponent(), CmsTrackerLevelBuilder::build(), ExtractStringFromDDD::getString(), GeometricDet::ladder, GeometricDet::rod, GeometricDet::strng, CmsTrackerLevelBuilder::theCmsTrackerStringToEnum, and CmsTrackerStringToEnum::type().

                                                                                           {

  CmsTrackerStringBuilder theCmsTrackerStringBuilder ;
  CmsTrackerRodBuilder theCmsTrackerRodBuilder;
  CmsTrackerLadderBuilder theCmsTrackerLadderBuilder;

  GeometricDet * subdet = new GeometricDet(&fv,theCmsTrackerStringToEnum.type(ExtractStringFromDDD::getString(s,&fv)));
  switch (theCmsTrackerStringToEnum.type(ExtractStringFromDDD::getString(s,&fv))){
  case GeometricDet::strng:
    theCmsTrackerStringBuilder.build(fv,subdet,s);      
    break;
  case GeometricDet::rod:
    theCmsTrackerRodBuilder.build(fv,subdet,s);      
    break;
  case GeometricDet::ladder:
    theCmsTrackerLadderBuilder.build(fv,subdet,s);      
    break;
  default:
    edm::LogError("CmsTrackerLayerBuilder")<<" ERROR - I was expecting a String, Rod or Ladder, I got a "<<ExtractStringFromDDD::getString(s,&fv);

  }  
  g->addComponent(subdet);

}

void CmsTrackerLayerBuilder::sortNS ( DDFilteredView &  fv, GeometricDet *  det  )

privatevirtual

Reimplemented from CmsTrackerLevelBuilder.

Definition at line 40 of file CmsTrackerLayerBuilder.cc.

References GeometricDet::addComponents(), GeometricDet::clearComponents(), AlCaHLTBitMon_QueryRunRegistry::comp, GeometricDet::component(), GeometricDet::components(), i, GeometricDet::ladder, GeometricDet::params(), GeometricDet::rod, GeometricDet::setGeographicalID(), GeometricDet::strng, groupFilesInBlocks::temp, TrackerStablePhiSort(), and GeometricDet::translation().

                                                                        {

  GeometricDet::ConstGeometricDetContainer& comp = det->components();

  // TIB
  // SubDetector Side: 2 bits [TIB-:1 TIB+:2]
  // Layer Part      : 2 bits [internal:1 external:0]
  // String Number   : 6 bits [1,...,56 (at most)]
  //
  if(comp.front()->type()== GeometricDet::strng){
    float layerRadius = (det->params()[2]+det->params()[1])/2.;

    GeometricDet::GeometricDetContainer neg;
    GeometricDet::GeometricDetContainer pos;
    GeometricDet::GeometricDetContainer extneg;
    GeometricDet::GeometricDetContainer intneg;
    GeometricDet::GeometricDetContainer extpos;
    GeometricDet::GeometricDetContainer intpos;
    neg.clear();
    pos.clear();
    extneg.clear();
    intneg.clear();
    extpos.clear();
    intpos.clear();

    for(size_t i = 0; i< comp.size(); ++i) {
      auto component = det->component(i);
      if(component->translation().z()<0.){
    neg.push_back(component);
      }else{
    pos.push_back(component);
      }
    }

    for(GeometricDet::GeometricDetContainer::iterator i=neg.begin();i!=neg.end();i++){
      double rPos = (*i)->translation().Rho();
      if(rPos > layerRadius ){ 
    extneg.push_back(*i);
      }else{
    intneg.push_back(*i);
      }
    }

    for(GeometricDet::GeometricDetContainer::iterator i=pos.begin();i!=pos.end();i++){
      double rPos = (*i)->translation().Rho();
      if(rPos > layerRadius ){ 
    extpos.push_back(*i);
      }else{
    intpos.push_back(*i);
      }
    }

    TrackerStablePhiSort(extneg.begin(), extneg.end(), ExtractPhi());
    TrackerStablePhiSort(extpos.begin(), extpos.end(), ExtractPhi());
    TrackerStablePhiSort(intneg.begin(), intneg.end(), ExtractPhi());
    TrackerStablePhiSort(intpos.begin(), intpos.end(), ExtractPhi());

    for(uint32_t i=0;i<intneg.size();i++){
      uint32_t temp=i+1;
      temp|=(1<<8); // 1 : SubDetector Side TIB-
      temp|=(1<<6); // 1 : Layer Part int
      intneg[i]->setGeographicalID(DetId(temp));
    }

    for(uint32_t i=0;i<extneg.size();i++){
      uint32_t temp=i+1;
      temp|=(1<<8); // 1 : SubDetector Side TIB-
      temp|=(2<<6); // 2 : Layer Part ext
      extneg[i]->setGeographicalID(DetId(temp));
    }

    for(uint32_t i=0;i<intpos.size();i++){
      uint32_t temp=i+1;
      temp|=(2<<8); // 2 : SubDetector Side TIB+
      temp|=(1<<6); // 1 : Layer Part int 
      intpos[i]->setGeographicalID(DetId(temp));
    }

    for(uint32_t i=0;i<extpos.size();i++){
      uint32_t temp=i+1;
      temp|=(2<<8); // 2 : SubDetector Side TIB+
      temp|=(2<<6); // 2 : Layer Part ext 
      extpos[i]->setGeographicalID(DetId(temp));
    }
    
    
    det->clearComponents();
    det->addComponents(intneg);
    det->addComponents(extneg);
    det->addComponents(intpos);
    det->addComponents(extpos);
    
  }else if(comp.front()->type()== GeometricDet::rod){
    GeometricDet::GeometricDetContainer neg;
    GeometricDet::GeometricDetContainer pos;
    neg.clear();
    pos.clear();
    
    for(size_t i=0; i<comp.size(); ++i) {
      auto component = det->component(i);
      if(component->translation().z()<0.){
    neg.push_back(component);
      }else{
    pos.push_back(component);
      }
    }

    TrackerStablePhiSort(neg.begin(), neg.end(), ExtractPhi());
    TrackerStablePhiSort(pos.begin(), pos.end(), ExtractPhi());
    
    for(uint32_t i=0; i<neg.size();i++){      
      uint32_t temp = i+1;
      temp|=(1<<7);
      neg[i]->setGeographicalID(DetId(temp));
    }
    
    for(uint32_t i=0; i<pos.size();i++){
      uint32_t temp = i+1;
      temp|=(2<<7);
      pos[i]->setGeographicalID(DetId(temp));
    }
    
    det->clearComponents();
    det->addComponents(neg);
    det->addComponents(pos);
    
  }else if(det->components().front()->type()== GeometricDet::ladder){

    TrackerStablePhiSort(comp.begin(), comp.end(), ExtractPhi());
    
    for(uint32_t i=0; i<comp.size();i++){
      det->component(i)->setGeographicalID(DetId(i+1));
    }    

  }else{
    edm::LogError("CmsTrackerLayerBuilder")<<"ERROR - wrong SubDet to sort..... "<<det->components().front()->type();
  }

}