Phase2EndcapRingBuilder.cc

Go to the documentation of this file.

#include "Phase2EndcapRingBuilder.h"
#include "TrackingTools/DetLayers/interface/DetLayerException.h"

using namespace edm;
using namespace std;

Phase2EndcapRing* Phase2EndcapRingBuilder::build(const GeometricDet* aPhase2EndcapRing,
                             const TrackerGeometry* theGeomDetGeometry,
                             const bool useBrothers)
{
  vector<const GeometricDet*>  allGeometricDets = aPhase2EndcapRing->components();
  vector<const GeometricDet*>  compGeometricDets;
  LogDebug("TkDetLayers") << "Phase2EndcapRingBuilder with #Modules: " << allGeometricDets.size() << std::endl;

  vector<const GeomDet*> frontGeomDets;
  vector<const GeomDet*> backGeomDets;
  double meanZ = 0;

  if(!useBrothers){
  
    //---- to evaluate meanZ
    for(vector<const GeometricDet*>::const_iterator compGeometricDets=allGeometricDets.begin(); compGeometricDets!=allGeometricDets.end();compGeometricDets++){
        LogTrace("TkDetLayers") << " compGeometricDets->positionBounds().perp() " << (*compGeometricDets)->positionBounds().z() << std::endl;
        meanZ = meanZ + (*compGeometricDets)->positionBounds().z();
    }
    meanZ = meanZ/allGeometricDets.size();
    LogDebug("TkDetLayers") << " meanZ " << meanZ << std::endl;
  
    for(vector<const GeometricDet*>::const_iterator compGeometricDets=allGeometricDets.begin();
        compGeometricDets!=allGeometricDets.end();compGeometricDets++){
      const GeomDet* theGeomDet = theGeomDetGeometry->idToDet( (*compGeometricDets)->geographicalID() );
  
      if( fabs( (*compGeometricDets)->positionBounds().z() ) < fabs(meanZ) )
        frontGeomDets.push_back(theGeomDet);
  
      if( fabs( (*compGeometricDets)->positionBounds().z() ) > fabs(meanZ) )
        backGeomDets.push_back(theGeomDet);      

      if( fabs( (*compGeometricDets)->positionBounds().z() ) == fabs(meanZ) )
        throw DetLayerException("Not possible to assiciate this GeometricDet in front or back");
  
    }
  
    LogDebug("TkDetLayers") << "frontGeomDets.size(): " << frontGeomDets.size() ;
    LogDebug("TkDetLayers") << "backGeomDets.size(): " << backGeomDets.size() ;
  
    return new Phase2EndcapRing(frontGeomDets,backGeomDets);

  } else {

    vector<const GeomDet*> frontGeomDetBrothers;
    vector<const GeomDet*> backGeomDetBrothers;
    vector<const GeometricDet*>  compGeometricDets;
  
    //---- to evaluate meanZ
    double meanZ = 0;
    double meanZBrothers = 0;
    for(vector<const GeometricDet*>::const_iterator it=allGeometricDets.begin(); it!=allGeometricDets.end();it++){
      compGeometricDets = (*it)->components();
      if (compGeometricDets.size() != 2){
        throw DetLayerException("Phase2OTEndcapRing is considered as a stack but does not have two components");
      } else {
        LogTrace("TkDetLayers") << " compGeometricDets[0]->positionBounds().perp() " << compGeometricDets[0]->positionBounds().z() << std::endl;
        LogTrace("TkDetLayers") << " compGeometricDets[1]->positionBounds().perp() " << compGeometricDets[1]->positionBounds().z() << std::endl;
        meanZ = meanZ + compGeometricDets[0]->positionBounds().z();
        meanZBrothers = meanZBrothers + compGeometricDets[1]->positionBounds().z();
      }

    }
    meanZ = meanZ/allGeometricDets.size();
    meanZBrothers = meanZBrothers/allGeometricDets.size();
    LogDebug("TkDetLayers") << " meanZ Lower " << meanZ << std::endl;
    LogDebug("TkDetLayers") << " meanZ Upper " << meanZBrothers << std::endl;
  
    for(vector<const GeometricDet*>::const_iterator it=allGeometricDets.begin();
        it!=allGeometricDets.end();it++){
      compGeometricDets = (*it)->components(); 
      const GeomDet* theGeomDet = theGeomDetGeometry->idToDet( compGeometricDets[0]->geographicalID() );
  
      if( fabs( compGeometricDets[0]->positionBounds().z() ) < fabs(meanZ) )
        frontGeomDets.push_back(theGeomDet);
  
      if( fabs( compGeometricDets[0]->positionBounds().z() ) > fabs(meanZ) )
        backGeomDets.push_back(theGeomDet);      
  
      const GeomDet* theGeomDetBrother = theGeomDetGeometry->idToDet( compGeometricDets[1]->geographicalID() );
  
      if( fabs( compGeometricDets[1]->positionBounds().z() ) < fabs(meanZBrothers) )
        frontGeomDetBrothers.push_back(theGeomDetBrother);
  
      if( fabs( compGeometricDets[1]->positionBounds().z() ) > fabs(meanZBrothers) )
        backGeomDetBrothers.push_back(theGeomDetBrother);

      if( fabs( compGeometricDets[0]->positionBounds().z() ) == fabs(meanZ) || fabs( compGeometricDets[1]->positionBounds().z() ) == fabs(meanZBrothers) )
        throw DetLayerException("Not possible to assiciate components of this GeometricDet in front or back");
    }
  
    LogDebug("TkDetLayers") << "frontGeomDets.size(): " << frontGeomDets.size() ;
    LogDebug("TkDetLayers") << "backGeomDets.size(): " << backGeomDets.size() ;
    LogDebug("TkDetLayers") << "frontGeomDetsBro.size(): " << frontGeomDetBrothers.size() ;
    LogDebug("TkDetLayers") << "backGeomDetsBro.size(): " << backGeomDetBrothers.size() ;
  
    return new Phase2EndcapRing(frontGeomDets,backGeomDets,frontGeomDetBrothers,backGeomDetBrothers);
  }
  return new Phase2EndcapRing(frontGeomDets,backGeomDets);

}