MagBLayer.cc

Go to the documentation of this file.

// #include "Utilities/Configuration/interface/Architecture.h"

/*
 *  See header file for a description of this class.
 *
 *  \author N. Amapane - INFN Torino
 */

#include "MagneticField/Layers/interface/MagBLayer.h"
#include "MagneticField/Layers/interface/MagBSector.h"
#include "MagneticField/VolumeGeometry/interface/MagVolume.h"
#include "Utilities/BinningTools/interface/PeriodicBinFinderInPhi.h"
#include "DataFormats/GeometryVector/interface/Phi.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <iostream>

using namespace std;

MagBLayer::MagBLayer(vector<MagBSector*>& sectors, double rMin)
    : theSectors(sectors), theSingleVolume(nullptr), theRMin(rMin), theBinFinder(nullptr) {
  //   LogTrace("MagGeometry") << "Building MagBLayer with " << theSectors.size()
  //          << " sectors, minR " << theRMin << endl;
  //FIXME: PeriodicBinFinderInPhi gets *center* of first bin
  theBinFinder = new PeriodicBinFinderInPhi<float>(theSectors.front()->minPhi() + Geom::pi() / 12., 12);
}

MagBLayer::MagBLayer(MagVolume* aVolume, double rMin) : theSingleVolume(nullptr), theRMin(rMin), theBinFinder(nullptr) {
  //   LogTrace("MagGeometry") << "Building MagBLayer with " << 0
  //          << " sectors, minR " << theRMin << endl;
}

MagBLayer::~MagBLayer() {
  delete theBinFinder;

  delete theSingleVolume;

  for (vector<MagBSector*>::const_iterator isec = theSectors.begin(); isec != theSectors.end(); ++isec) {
    delete (*isec);
  }
}

const MagVolume* MagBLayer::findVolume(const GlobalPoint& gp, double tolerance) const {
  const MagVolume* result = nullptr;

  //In case the layer is composed of a single volume...
  if (theSingleVolume) {
    //    LogTrace("MagGeometry") << "   Trying the unique volume " << endl;
    if (theSingleVolume->inside(gp, tolerance)) {
      result = theSingleVolume;
      //       LogTrace("MagGeometry") << "***In unique bsector"
      //          << (result==0? " failed " : " OK ") <<endl;
    }
    return result;
  }

  // Normal cases - query the sectors.

  Geom::Phi<float> phi = gp.phi();

  // FIXME assume sectors are sorted in phi!
  int bin = theBinFinder->binIndex(phi);
  LogTrace("MagGeometry") << "   Trying sector at phi " << theSectors[bin]->minPhi() << " " << phi << endl;
  result = theSectors[bin]->findVolume(gp, tolerance);
  LogTrace("MagGeometry") << "***In guessed bsector" << (result == nullptr ? " failed " : " OK ") << endl;

  if (result == nullptr) {  // If fails, can be in previous bin.
    LogTrace("MagGeometry") << "   Trying sector at phi " << theSectors[theBinFinder->binIndex(bin - 1)]->minPhi()
                            << " " << phi << endl;

    result = theSectors[theBinFinder->binIndex(bin - 1)]->findVolume(gp, tolerance);
    LogTrace("MagGeometry") << "***In previous bsector" << (result == nullptr ? " failed " : " OK ") << endl;
  }
  return result;
}