magneticfield::bLayer Class Reference

#include <bLayer.h>

Public Member Functions
	bLayer (handles::const_iterator begin, handles::const_iterator end, bool debugFlag=false)
	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 ()=default
	Destructor. More...

Private Member Functions
int	bin (int i) const

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

Detailed Description

Definition at line 17 of file bLayer.h.

Constructor & Destructor Documentation

bLayer()

bLayer::bLayer	(	handles::const_iterator	begin,
		handles::const_iterator	end,
		bool	debugFlag = false
	)

Constructor from list of volumes.

Definition at line 20 of file bLayer.cc.

References bin(), gather_cfg::cout, ddbox, ddtrap, debug, mps_fire::i, hltrates_dqm_sourceclient-live_cfg::offset, SurfaceOrientation::outer, precomputed_value_sort(), magneticfield::printUniqueNames(), svgfig::rotate(), GloballyPositioned< T >::rotation(), sectors, size, theVolumes, tolerance, TkRotation< T >::zx(), TkRotation< T >::zy(), and TkRotation< T >::zz().

    : size(end - begin), theVolumes(begin, end), mlayer(nullptr), debug(debugFlag) {
  // Sort in phi
  precomputed_value_sort(theVolumes.begin(), theVolumes.end(), ExtractPhi());

  if (debug) {
    std::cout << " elements: " << theVolumes.size() << " unique volumes: ";
    printUniqueNames(theVolumes.begin(), theVolumes.end());
  }

  // Find sectors in phi
  handles::iterator secBegin = theVolumes.begin();
  handles::iterator secEnd = secBegin;
  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 (debug)
      std::cout << "      Sector is just one volume." << std::endl;

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

  } else {                    // there are more than one volume per sector.
    float tolerance = 0.025;  // 250 micron
    do {
      if (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 (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), debug);
    }
  }

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

~bLayer()

magneticfield::bLayer::~bLayer ( )

default

Destructor.

Member Function Documentation

bin()

int bLayer::bin ( int  i ) const

private

Definition at line 109 of file bLayer.cc.

References mps_fire::i, and size.

Referenced by bLayer().

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

buildMagBLayer()

MagBLayer * bLayer::buildMagBLayer

(

)

const

Construct the MagBLayer upon request.

Definition at line 130 of file bLayer.cc.

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

                                        {
  if (mlayer == nullptr) {
    // If we have only one volume, do not build any MagBSector.
    if (sectors.empty()) {
      if (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;
}

minR()

double bLayer::minR

(

)

const

Return sector at i (handling periodicity)

Min R (conservative guess).

Definition at line 120 of file bLayer.cc.

References theVolumes.

Referenced by buildMagBLayer().

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

RN()

const float magneticfield::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(); }

volumes()

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

inline

Return the list of all volumes.

Definition at line 29 of file bLayer.h.

References theVolumes.

{ return theVolumes; }

Member Data Documentation

debug

const bool magneticfield::bLayer::debug

private

Definition at line 54 of file bLayer.h.

Referenced by bLayer(), buildMagBLayer(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), and runTauIdMVA.TauIDEmbedder::runTauID().

mlayer

MagBLayer* magneticfield::bLayer::mlayer

mutableprivate

Definition at line 53 of file bLayer.h.

Referenced by buildMagBLayer().

sectors

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

private

Definition at line 50 of file bLayer.h.

Referenced by bLayer(), and buildMagBLayer().

size

int magneticfield::bLayer::size

private

Definition at line 45 of file bLayer.h.

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

theVolumes

handles magneticfield::bLayer::theVolumes

private

Definition at line 51 of file bLayer.h.

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