magneticfield::InterpolatorBuilder Class Reference

#include <InterpolatorBuilder.h>

Public Member Functions
std::unique_ptr< MagProviderInterpol >	build (volumeHandle const *)
	InterpolatorBuilder (std::string iTableSet, bool useMergeFileIfAvailable=true)
	InterpolatorBuilder (const InterpolatorBuilder &)=delete
const InterpolatorBuilder &	operator= (const InterpolatorBuilder &)=delete

Private Attributes
std::unordered_map< std::string, unsigned int >	offsets_
std::optional< interpolation::binary_ifstream >	stream_
std::string	tableSet_

Detailed Description

Definition at line 34 of file InterpolatorBuilder.h.

Constructor & Destructor Documentation

InterpolatorBuilder() [1/2]

InterpolatorBuilder::InterpolatorBuilder	(	std::string	iTableSet,
		bool	useMergeFileIfAvailable = true
	)

Definition at line 40 of file InterpolatorBuilder.cc.

References edm::FileInPath::findFile(), eostools::move(), HLT_IsoTrack_cff::offset, offsets_, stream_, AlCaHLTBitMon_QueryRunRegistry::string, and tableSet_.

      : tableSet_(std::move(iTableSet)) {
    if (not useMergeFileIfAvailable)
      return;
    auto indexFileName = edm::FileInPath::findFile("MagneticField/Interpolation/data/" + tableSet_ + "/merged.index");
    if (not indexFileName.empty()) {
      auto binaryFileName = edm::FileInPath::findFile("MagneticField/Interpolation/data/" + tableSet_ + "/merged.bin");
      if (not binaryFileName.empty()) {
        std::ifstream indexFile(indexFileName);
        while (indexFile) {
          std::string magFile;
          unsigned int offset;
          indexFile >> magFile >> offset;
          offsets_.emplace(std::move(magFile), offset);
        }
        stream_ = interpolation::binary_ifstream(binaryFileName);
      }
    }
  }

InterpolatorBuilder() [2/2]

magneticfield::InterpolatorBuilder::InterpolatorBuilder ( const InterpolatorBuilder &  )

delete

Member Function Documentation

build()

std::unique_ptr< MagProviderInterpol > InterpolatorBuilder::build

(

volumeHandle const *

vol

)

Definition at line 63 of file InterpolatorBuilder.cc.

References MFGridFactory::build(), Exception, spr::find(), edm::FileInPath::findFile(), contentValuesFiles::fullPath, magneticfield::BaseVolumeHandle::magFile, magneticfield::BaseVolumeHandle::masterSector, offsets_, magneticfield::BaseVolumeHandle::placement(), GloballyPositioned< T >::position(), hcal_runs::rf, makeMuonMisalignmentScenario::rot, GloballyPositioned< T >::rotation(), stream_, and tableSet_.

Referenced by magneticfield::MagGeoBuilder::buildInterpolator().

                                                                                       {
    if (tableSet_ == "fake" || vol->magFile == "fake") {
      return std::make_unique<magneticfield::FakeInterpolator>();
    }

    // If the table is in "local" coordinates, must create a reference
    // frame that is appropriately rotated along the CMS Z axis.

    GloballyPositioned<float> rf = *(vol->placement());

    if (vol->masterSector != 1) {
      typedef Basic3DVector<float> Vector;

      // Phi of the master sector
      double masterSectorPhi = (vol->masterSector - 1) * 1._pi / 6.;

      GloballyPositioned<float>::RotationType rot(Vector(0, 0, 1), -masterSectorPhi);
      Vector vpos(vol->placement()->position());

      rf = GloballyPositioned<float>(GloballyPositioned<float>::PositionType(rot.multiplyInverse(vpos)),
                                     vol->placement()->rotation() * rot);
    }

    if (not stream_) {
      auto fullPath = edm::FileInPath::findFile("MagneticField/Interpolation/data/" + tableSet_ + "/" + vol->magFile);
      if (fullPath.empty()) {
        //cause the exception to happen
        edm::FileInPath mydata("MagneticField/Interpolation/data/" + tableSet_ + "/" + vol->magFile);
        return {};
      }

      magneticfield::interpolation::binary_ifstream strm(fullPath);
      return std::unique_ptr<MagProviderInterpol>(MFGridFactory::build(strm, rf));
    }

    auto find = offsets_.find(vol->magFile);
    if (find == offsets_.end()) {
      throw cms::Exception("MissingMagFileEntry") << vol->magFile << " was not an entry in the index file";
    }
    stream_->seekg(find->second);
    if (stream_->fail()) {
      throw cms::Exception("SeekMagFileEntry") << " failed seekg within merged binary file";
    }
    return std::unique_ptr<MagProviderInterpol>(MFGridFactory::build(*stream_, rf));
  }

operator=()

const InterpolatorBuilder& magneticfield::InterpolatorBuilder::operator= ( const InterpolatorBuilder &  )

delete

Member Data Documentation

offsets_

std::unordered_map<std::string, unsigned int> magneticfield::InterpolatorBuilder::offsets_

private

Definition at line 51 of file InterpolatorBuilder.h.

Referenced by build(), and InterpolatorBuilder().

stream_

std::optional<interpolation::binary_ifstream> magneticfield::InterpolatorBuilder::stream_

private

Definition at line 52 of file InterpolatorBuilder.h.

Referenced by build(), and InterpolatorBuilder().

tableSet_

std::string magneticfield::InterpolatorBuilder::tableSet_

private

Definition at line 50 of file InterpolatorBuilder.h.

Referenced by build(), and InterpolatorBuilder().