magneticfield::eSector Class Reference

#include <eSector.h>

Public Member Functions
MagESector *	buildMagESector () const
	Construct the MagESector upon request. More...
	eSector (handles::const_iterator begin, handles::const_iterator end, bool debugFlag=false)
	Constructor from list of volumes. More...
	~eSector ()=default
	Destructor. More...

Private Attributes
const bool	debug
std::vector< eLayer >	layers
MagESector *	msector
handles	theVolumes

Detailed Description

Definition at line 16 of file eSector.h.

Constructor & Destructor Documentation

eSector()

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

Constructor from list of volumes.

Definition at line 21 of file eSector.cc.

References ClusterizingHistogram::clusterize(), gather_cfg::cout, debug, f, ClusterizingHistogram::fill(), dqmdumpme::first, mps_fire::i, dqmdumpme::last, layers, precomputed_value_sort(), magneticfield::printUniqueNames(), dtDQMClient_cfg::resolution, theVolumes, SiStripMonitorCluster_cfi::zmax, and SiStripMonitorCluster_cfi::zmin.

    : theVolumes(begin, end), msector(nullptr), debug(debugFlag) {
  //FIXME!!!
  //precomputed_value_sort(theVolumes.begin(), theVolumes.end(), ExtractAbsZ());
  precomputed_value_sort(theVolumes.begin(), theVolumes.end(), ExtractZ());

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

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

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

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

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

  handles::const_iterator layStart = 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) {
      cout << "     Layer at: " << zClust[i] << " elements: " << separ - layStart << " unique volumes: ";
      printUniqueNames(layStart, separ);
    }

    layers.push_back(eLayer(layStart, separ));
    layStart = separ;
  }
  {
    if (debug) {
      cout << "     Layer at: " << zClust.back() << " elements: " << last - separ << " unique volumes: ";
      printUniqueNames(separ, last);
    }
    layers.push_back(eLayer(separ, last));
  }

  // FIXME: Check that all layers have the same dz?.
}

~eSector()

magneticfield::eSector::~eSector ( )

default

Destructor.

Member Function Documentation

buildMagESector()

MagESector * eSector::buildMagESector

(

)

const

Construct the MagESector upon request.

Definition at line 74 of file eSector.cc.

References layers, msector, and theVolumes.

                                           {
  if (msector == nullptr) {
    vector<MagELayer*> mLayers;
    for (vector<eLayer>::const_iterator lay = layers.begin(); lay != layers.end(); ++lay) {
      mLayers.push_back((*lay).buildMagELayer());
    }
    msector = new MagESector(mLayers, theVolumes.front()->minPhi());  //FIXME
  }
  return msector;
}

Member Data Documentation

debug

const bool magneticfield::eSector::debug

private

Definition at line 34 of file eSector.h.

Referenced by eSector(), runTauIdMVA.TauIDEmbedder::load_againstElectronMVA6(), and runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017().

layers

std::vector<eLayer> magneticfield::eSector::layers

private

Definition at line 31 of file eSector.h.

Referenced by buildMagESector(), and eSector().

msector

MagESector* magneticfield::eSector::msector

mutableprivate

Definition at line 33 of file eSector.h.

Referenced by buildMagESector().

theVolumes

handles magneticfield::eSector::theVolumes

private

Definition at line 32 of file eSector.h.

Referenced by buildMagESector(), and eSector().