MagGeoBuilderFromDDD::bLayer Class Reference

#include <bLayer.h>

Public Member Functions
	bLayer (handles::const_iterator begin, handles::const_iterator end)
	Constructor from list of volumes. More...
MagBLayer *	buildMagBLayer () const
	Construct the MagBLayer upon request. More...
double	minR () const
	Return sector at i (handling periodicity) More...
const float	RN () const
	Distance from center along normal of sectors. More...
const handles &	volumes () const
	Return the list of all volumes. More...
	~bLayer ()
	Destructor. More...

Private Member Functions
int	bin (int i) const

Private Attributes
MagBLayer *	mlayer
std::vector< bSector >	sectors
int	size
handles	theVolumes

Detailed Description

A layer of barrel volumes. Holds a list of volumes and 12 sectors. It is assumed that the geometry is 12-fold periodic in phi!

Author

N. Amapane - INFN Torino

Definition at line 17 of file bLayer.h.

Constructor & Destructor Documentation

MagGeoBuilderFromDDD::bLayer::bLayer	(	handles::const_iterator	begin,
		handles::const_iterator	end
	)

Constructor from list of volumes.

Definition at line 19 of file bLayer.cc.

References bin(), gather_cfg::cout, ddbox, ddtrap, MagGeoBuilderFromDDD::debug, i, PFRecoTauDiscriminationByIsolation_cfi::offset, SurfaceOrientation::outer, precomputed_value_sort(), MagGeoBuilderFromDDD::volumeHandle::printUniqueNames(), svgfig::rotate(), GloballyPositioned< T >::rotation(), sectors, size, theVolumes, RecoTauCleanerPlugins::tolerance, TkRotation< T >::zx(), TkRotation< T >::zy(), and TkRotation< T >::zz().

                                               :
  size(end-begin),
  theVolumes(begin,end),
  mlayer(0) 
{
  // Sort in phi
  precomputed_value_sort(theVolumes.begin(), theVolumes.end(), ExtractPhi());
  
  if (MagGeoBuilderFromDDD::debug) {
    std::cout << " elements: " << theVolumes.size() << " unique volumes: ";
    volumeHandle::printUniqueNames(theVolumes.begin(), theVolumes.end());
  }
  

  // Find sectors in phi
  handles::iterator secBegin = theVolumes.begin();
  handles::iterator secEnd;
  int binOffset = 0;

  const Surface & refSurf = (*secBegin)->surface(outer);

  int newbegin=0;
  int newend=0;  

  // A sector is made of several volumes in R, and, for planar layers
  // (box and traps) also in phi, so it might cross the -phi boundary. 
  // For those, we have to look for the end of first sector and rotate the 
  // vector of volumes.
  // ASSUMPTION: all volumes in a layer must be of compatible type.


  if (size==1) { // Only one volume; this is the case for barrel
    // cylinders.
    // FIXME sectors.push_back(bSector(theVolumes.begin(),theVolumes.end());
    if (MagGeoBuilderFromDDD::debug) std::cout <<"      Sector is just one volume." << std::endl;

  } else if (size==12 || // In this case, each volume is a sector.
         (((*secBegin)->shape()!=ddtrap) && (*secBegin)->shape()!=ddbox)) {
    secEnd = secBegin+size/12;

  }  else { // there are more than one volume per sector.
    float tolerance = 0.025; // 250 micron
    do {
      if (MagGeoBuilderFromDDD::debug) std::cout << (*secBegin)->name 
                 << " " << (*secBegin)->copyno << std::endl;
      ++secBegin;
    } while ((secBegin != theVolumes.end()) &&
         (*secBegin)->sameSurface(refSurf,outer, tolerance)); // This works only if outer surface is a plane, otherwise sameSurface returns always true!
    
    secEnd = secBegin;
    secBegin = theVolumes.begin()+bin((secEnd-theVolumes.begin())-size/12);;
    newend   = secEnd-theVolumes.begin();
    newbegin = secBegin-theVolumes.begin();
    
    // Rotate the begin of the first sector to the vector beginning.
    rotate(theVolumes.begin(),secBegin,theVolumes.end());
    secBegin = theVolumes.begin();
    secEnd   = secBegin+size/12;

    // Test it is correct...
    if (!((*secBegin)->sameSurface((*(secEnd-1))->surface(outer),
                   outer, tolerance))) {
      std::cout << "*** ERROR: Big mess while looking for sectors "
        << (*secBegin)->name << " " << (*secBegin)->copyno << " "
        << (*(secEnd-1))->name << " " << (*(secEnd-1))->copyno
        << std::endl;
    }
  }

  if (MagGeoBuilderFromDDD::debug) {
    std::cout << "      First sector: volumes " << secEnd-theVolumes.begin()
     << " from " << newbegin
     << " (phi = " << (*secBegin)->center().phi() << ") "
     << " to " << newend
     << " (phi = " << (*secEnd)->center().phi() << ") "
     << " # " << (*secBegin)->copyno << " ";
    std::cout << GlobalVector( refSurf.rotation().zx(), refSurf.rotation().zy(),
              refSurf.rotation().zz()) << std::endl;
  }

  if (size!=1) { // Build the 12 sectors
    int offset = size/12;
    sectors.resize(12);
    for (int i = 0; i<12; ++i) {
      int isec = (i+binOffset)%12;
      sectors[isec>=0?isec:isec+12] = bSector(theVolumes.begin()+((i)*offset),
                          theVolumes.begin()+((i+1)*offset));
    }
  }

  if (MagGeoBuilderFromDDD::debug) std::cout << "-----------------------" << std::endl;

}

MagGeoBuilderFromDDD::bLayer::~bLayer

(

)

Destructor.

Definition at line 114 of file bLayer.cc.

{}

Member Function Documentation

int MagGeoBuilderFromDDD::bLayer::bin ( int  i ) const

private

Definition at line 116 of file bLayer.cc.

References size.

Referenced by bLayer().

                                               {
  i = i%size;
  return (i>=0?i:i+size);
}

MagBLayer * MagGeoBuilderFromDDD::bLayer::buildMagBLayer

(

)

const

Construct the MagBLayer upon request.

Definition at line 138 of file bLayer.cc.

References gather_cfg::cout, MagGeoBuilderFromDDD::debug, i, minR(), mlayer, sectors, size, and theVolumes.

Referenced by volumes().

                                                             {
  if (mlayer==0) {

    // If we have only one volume, do not build any MagBSector.
    if (sectors.size()==0) {
      if (MagGeoBuilderFromDDD::debug && size!=0) {
    std::cout << "ERROR: bLayer::buildMagBLayer, 0 sectors but "
         << size << " volumes" << std::endl;
      }
      // Technically we might have only one bSector built and we would
      // not need a separate MagBLayer constructor...
      mlayer = new MagBLayer(theVolumes.front()->magVolume, minR());
    }

    // If we have several sectors, create the MagBSector
    std::vector<MagBSector*> mSectors;
    for (unsigned int i=0; i<sectors.size(); ++i) {
      mSectors.push_back(sectors[i].buildMagBSector());
    }
    mlayer = new MagBLayer(mSectors, minR());
  }
  return mlayer;
}

double MagGeoBuilderFromDDD::bLayer::minR

(

)

const

Return sector at i (handling periodicity)

Min R (conservative guess).

Definition at line 128 of file bLayer.cc.

References theVolumes.

Referenced by buildMagBLayer(), and volumes().

                                              {
  // ASSUMPTION: a layer is only 1 volume thick (by construction). 
  return theVolumes.front()->minR();
}

const float MagGeoBuilderFromDDD::bLayer::RN ( ) const

inline

Distance from center along normal of sectors.

Definition at line 26 of file bLayer.h.

References theVolumes.

                         {
    return theVolumes.front()->RN();
  }

const handles& MagGeoBuilderFromDDD::bLayer::volumes ( ) const

inline

Return the list of all volumes.

Definition at line 31 of file bLayer.h.

References buildMagBLayer(), minR(), and theVolumes.

{return theVolumes;}

Member Data Documentation

MagBLayer* MagGeoBuilderFromDDD::bLayer::mlayer

mutableprivate

Definition at line 55 of file bLayer.h.

Referenced by buildMagBLayer().

std::vector<bSector> MagGeoBuilderFromDDD::bLayer::sectors

private

Definition at line 52 of file bLayer.h.

Referenced by bLayer(), and buildMagBLayer().

int MagGeoBuilderFromDDD::bLayer::size

private

Definition at line 47 of file bLayer.h.

Referenced by ntuple._Collection::__iter__(), ntuple._Collection::__len__(), bin(), bLayer(), buildMagBLayer(), and cuy.FindIssue::startElement().

handles MagGeoBuilderFromDDD::bLayer::theVolumes

private

Definition at line 53 of file bLayer.h.

Referenced by bLayer(), buildMagBLayer(), minR(), RN(), and volumes().