magneticfield::bRod Class Reference

#include <bRod.h>

Public Member Functions
	bRod (handles::const_iterator begin, handles::const_iterator end, bool debugVal=false)
	Constructor from list of volumes. More...
MagBRod *	buildMagBRod () const
	Construct the MagBRod upon request. More...
const float	RN () const
	Distance from center along sector normal. More...
	~bRod ()=default
	Destructor. More...

Private Attributes
bool	debug
MagBRod *	mrod
std::vector< bSlab >	slabs
handles	volumes

Detailed Description

Definition at line 17 of file bRod.h.

Constructor & Destructor Documentation

bRod::bRod	(	handles::const_iterator	begin,
		handles::const_iterator	end,
		bool	debugVal = false
	)

Constructor from list of volumes.

Definition at line 19 of file bRod.cc.

References ClusterizingHistogram::clusterize(), gather_cfg::cout, debug, validate-o2o-wbm::f, ClusterizingHistogram::fill(), first, mps_fire::i, dqmdumpme::last, phimax, phimin, precomputed_value_sort(), magneticfield::printUniqueNames(), dtDQMClient_cfg::resolution, slabs, and volumes.

    : volumes(begin, end), mrod(nullptr), debug(debugVal) {
  precomputed_value_sort(volumes.begin(), volumes.end(), ExtractZ());

  // Clusterize in Z
  const float resolution = 5.;  // cm
  float zmin = volumes.front()->center().z() - resolution;
  float zmax = volumes.back()->center().z() + resolution;
  ClusterizingHistogram hisZ(int((zmax - zmin) / resolution) + 1, zmin, zmax);

  if (debug)
    std::cout << "     Z slabs: " << zmin << " " << zmax << std::endl;

  handles::const_iterator first = volumes.begin();
  handles::const_iterator last = volumes.end();

  for (handles::const_iterator i = first; i != last; ++i) {
    hisZ.fill((*i)->center().z());
  }
  std::vector<float> zClust = hisZ.clusterize(resolution);

  if (debug)
    std::cout << "     Found " << zClust.size() << " clusters in Z, "
              << " slabs: " << std::endl;

  handles::const_iterator slabStart = first;
  handles::const_iterator separ = first;

  for (unsigned int i = 0; i < zClust.size() - 1; ++i) {
    float zSepar = (zClust[i] + zClust[i + 1]) / 2.f;
    while ((*separ)->center().z() < zSepar)
      ++separ;
    if (debug) {
      std::cout << "     Slab at: " << zClust[i] << " elements: " << separ - slabStart << " unique volumes: ";
      printUniqueNames(slabStart, separ);
    }

    slabs.push_back(bSlab(slabStart, separ, debug));
    slabStart = separ;
  }
  {
    if (debug) {
      std::cout << "     Slab at: " << zClust.back() << " elements: " << last - separ << " unique volumes: ";
      printUniqueNames(separ, last);
    }
    slabs.push_back(bSlab(separ, last, debug));
  }

  // Check that all slabs have the same dphi.
  std::vector<bSlab>::const_iterator i = slabs.begin();
  Geom::Phi<float> phimax = (*i).maxPhi();
  Geom::Phi<float> phimin = (*i).minPhi();
  for (++i; i != slabs.end(); ++i) {
    if (fabs(phimax - (*i).maxPhi()) > 0.001 || fabs(phimin - (*i).minPhi()) > 0.001) {
      if (debug) {
        std::cout << "*** WARNING: slabs in this rod have different dphi! minphi " << phimin;
        std::cout << " != " << (*i).minPhi() << " or maxphi " << phimax << " != " << (*i).maxPhi() << std::endl;
      }
    }
  }
}

magneticfield::bRod::~bRod ( )

default

Destructor.

Member Function Documentation

MagBRod * bRod::buildMagBRod

(

)

const

Construct the MagBRod upon request.

Definition at line 81 of file bRod.cc.

References mrod, and slabs.

                                  {
  if (mrod == nullptr) {
    std::vector<MagBSlab*> mSlabs;
    for (std::vector<bSlab>::const_iterator slab = slabs.begin(); slab != slabs.end(); ++slab) {
      mSlabs.push_back((*slab).buildMagBSlab());
    }
    mrod = new MagBRod(mSlabs, slabs.front().minPhi());  //FIXME
  }
  return mrod;
}

const float magneticfield::bRod::RN ( ) const

inline

Distance from center along sector normal.

Definition at line 26 of file bRod.h.

References volumes.

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

Member Data Documentation

bool magneticfield::bRod::debug

private

Definition at line 35 of file bRod.h.

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

MagBRod* magneticfield::bRod::mrod

mutableprivate

Definition at line 34 of file bRod.h.

Referenced by buildMagBRod().

std::vector<bSlab> magneticfield::bRod::slabs

private

Definition at line 32 of file bRod.h.

Referenced by bRod(), and buildMagBRod().

handles magneticfield::bRod::volumes

private

Definition at line 33 of file bRod.h.

Referenced by bRod(), and RN().