CmsTrackerDiskBuilder Class Reference

#include <CmsTrackerDiskBuilder.h>

Inheritance diagram for CmsTrackerDiskBuilder:

List of all members.

Public Member Functions
	CmsTrackerDiskBuilder (unsigned int totalBlade)
Private Member Functions
virtual void	buildComponent (DDFilteredView &, GeometricDet *, std::string)
void	PhiPosNegSplit_innerOuter (std::vector< GeometricDet const * >::iterator begin, std::vector< GeometricDet const * >::iterator end)
virtual void	sortNS (DDFilteredView &, GeometricDet *)
Private Attributes
unsigned int	m_totalBlade

---

Detailed Description

Class which contructs PixelForward/Disk.

Definition at line 11 of file CmsTrackerDiskBuilder.h.

---

Constructor & Destructor Documentation

CmsTrackerDiskBuilder::CmsTrackerDiskBuilder

(

unsigned int

totalBlade

)

Definition at line 14 of file CmsTrackerDiskBuilder.cc.

  : m_totalBlade( totalBlade )
{}

---

Member Function Documentation

void CmsTrackerDiskBuilder::buildComponent	(	DDFilteredView &	fv,
		GeometricDet *	g,
		std::string	s
	)		[private, virtual]

Implements CmsTrackerLevelBuilder.

Definition at line 71 of file CmsTrackerDiskBuilder.cc.

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

{
  CmsTrackerPanelBuilder theCmsTrackerPanelBuilder;
  GeometricDet * subdet = new GeometricDet( &fv, theCmsTrackerStringToEnum.type( ExtractStringFromDDD::getString( s, &fv )));

  switch( theCmsTrackerStringToEnum.type( ExtractStringFromDDD::getString( s, &fv )))
  {
  case GeometricDet::panel:
    theCmsTrackerPanelBuilder.build( fv, subdet, s );
    break;
  default:
    edm::LogError( "CmsTrackerDiskBuilder" ) << " ERROR - I was expecting a Panel, I got a " << ExtractStringFromDDD::getString( s, &fv );   
  }  
  g->addComponent( subdet );
}

void CmsTrackerDiskBuilder::PhiPosNegSplit_innerOuter	(	std::vector< GeometricDet const * >::iterator	begin,
		std::vector< GeometricDet const * >::iterator	end
	)		[private]

Definition at line 25 of file CmsTrackerDiskBuilder.cc.

References begin, filterCSVwithJSON::copy, end, max(), min, PhiSort(), and python::multivaluedict::sort().

Referenced by sortNS().

{
  // first sort in phi, lowest first (-pi to +pi)
  std::sort( begin, end, PhiSort );

  // now put positive phi (in order) ahead of negative phi as in std geometry
  std::vector<const GeometricDet*> theCompsPosNeg;
  theCompsPosNeg.empty();
  theCompsPosNeg.clear();
  // also find the average radius (used to split inner and outer disk panels)
  double theRmin = (**begin).rho();
  double theRmax = theRmin;
  for(vector<const GeometricDet*>::const_iterator it=begin;
      it!=end;it++){
    if((**it).phi() >= 0) theCompsPosNeg.push_back(*it);
    theRmin = std::min( theRmin, (**it).rho());
    theRmax = std::max( theRmax, (**it).rho());
  }
  for(vector<const GeometricDet*>::const_iterator it=begin;
      it!=end;it++){
    if((**it).phi() < 0) theCompsPosNeg.push_back(*it);
  }

  // now put inner disk panels first
  double radius_split = 0.5 * (theRmin + theRmax);
  std::vector<const GeometricDet*> theCompsInnerOuter;
  theCompsInnerOuter.empty();
  theCompsInnerOuter.clear();
  //unsigned int num_inner = 0;
  for(vector<const GeometricDet*>::const_iterator it=theCompsPosNeg.begin();
      it!=theCompsPosNeg.end();it++){
    if((**it).rho() <= radius_split) {
      theCompsInnerOuter.push_back(*it);
      //num_inner++;
    }
  }
  for(vector<const GeometricDet*>::const_iterator it=theCompsPosNeg.begin();
      it!=theCompsPosNeg.end();it++){
    if((**it).rho() > radius_split) theCompsInnerOuter.push_back(*it);
  }
  //std::cout << "num of inner = " << num_inner << " with radius less than " << radius_split << std::endl;
  std::copy(theCompsInnerOuter.begin(), theCompsInnerOuter.end(), begin);
}

void CmsTrackerDiskBuilder::sortNS	(	DDFilteredView &	fv,
		GeometricDet *	det
	)		[private, virtual]

Reimplemented from CmsTrackerLevelBuilder.

Definition at line 88 of file CmsTrackerDiskBuilder.cc.

References GeometricDet::addComponents(), GeometricDet::clearComponents(), AlCaHLTBitMon_QueryRunRegistry::comp, GeometricDet::components(), j, m_totalBlade, GeometricDet::panel, PhiPosNegSplit_innerOuter(), RingESSource_cfi::rings, groupFilesInBlocks::temp, TrackerStablePhiSort(), and z.

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

// NP** BIG switch between Phase 1 and Outer Tracker Pixels
if( fabs( comp[0]->translation().z() ) < 1000 ) {

  //std::cerr<<"PHASE1"<<std::endl;


  switch( det->components().front()->type())
  {
  case GeometricDet::panel:
    if( m_totalBlade == 24 )
      TrackerStablePhiSort( comp.begin(), comp.end(), ExtractPhi());
    else
      PhiPosNegSplit_innerOuter( comp.begin(), comp.end());
    break;
  default:
    edm::LogError( "CmsTrackerDiskBuilder" ) << "ERROR - wrong SubDet to sort..... " << det->components().front()->type();
  }

  GeometricDet::GeometricDetContainer zminpanels;  // Here z refers abs(z);
  GeometricDet::GeometricDetContainer zmaxpanels;  // So, zmin panel is always closer to ip.

  uint32_t totalblade = comp.size()/2;
  if( totalblade != m_totalBlade )
    edm::LogError( "CmsTrackerDiskBuilder" ) << "ERROR, The Total Number of Blade in one disk is " << totalblade << "; configured " << m_totalBlade;

  zminpanels.reserve( totalblade );
  zmaxpanels.reserve( totalblade );
  for( uint32_t j = 0; j < totalblade; j++ )
  {
    if( fabs( comp[2*j]->translation().z()) > fabs( comp[ 2*j +1 ]->translation().z()))
    {
      zmaxpanels.push_back( comp[2*j] );
      zminpanels.push_back( comp[2*j+1] );

    }
    else if( fabs( comp[2*j]->translation().z()) < fabs( comp[ 2*j +1 ]->translation().z()))
    {
      zmaxpanels.push_back( comp[2*j+1] );
      zminpanels.push_back( comp[2*j] );
    }
    else
    {
      edm::LogWarning( "CmsTrackerDiskBuilder" ) << "WARNING - The Z of  both panels are equal! ";
    }
  }

  for( uint32_t fn = 0; fn < zminpanels.size(); fn++ )
  {
    uint32_t blade = fn + 1;
    uint32_t panel = 1;
    uint32_t temp = ( blade << 2 ) | panel;
    zminpanels[fn]->setGeographicalID( temp );
  }
  
  for( uint32_t bn = 0; bn < zmaxpanels.size(); bn++)
  {
    uint32_t blade = bn + 1;
    uint32_t panel = 2;
    uint32_t temp = ( blade << 2) | panel;
    zmaxpanels[bn]->setGeographicalID( temp );
  }
  
  det->clearComponents();
  det->addComponents( zminpanels );
  det->addComponents( zmaxpanels );
}
// NP** BIG switch between Phase 1 and Outer Tracker Pixels
else {

 //std::cerr<< "PHASE2!"<<std::endl;

  switch(det->components().front()->type()){
    case GeometricDet::panel:
   //   std::cerr<<comp.size()<<" RINGS!!"<<std::endl;
   //PhiPosNegSplit_innerOuter(comp.begin(),comp.end());
   //TrackerStablePhiSort(comp.begin(),comp.end(), ExtractPhi());
      break;
    default:
      edm::LogError("CmsTrackerDiskBuilder")<<"ERROR - wrong SubDet to sort..... "<<det->components().front()->type();
  }

  GeometricDet::GeometricDetContainer rings;
  uint32_t  totalrings = comp.size();

  for ( uint32_t rn=0; rn<totalrings; rn++) {
    rings.push_back(comp[rn]);
    uint32_t blade = rn+1;
    uint32_t panel = 1;
    uint32_t temp = (blade<<2) | panel;
    rings[rn]->setGeographicalID(temp);

  }

  det->clearComponents();
  det->addComponents(rings);
}
// NP** End of the BIG switch


}

---

Member Data Documentation

unsigned int CmsTrackerDiskBuilder::m_totalBlade [private]

Definition at line 22 of file CmsTrackerDiskBuilder.h.

Referenced by sortNS().