#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	(	handles::const_iterator	begin,
		handles::const_iterator	end,
		bool	debugFlag = `false`
	)

Constructor from list of volumes.

Definition at line 22 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(), L1TObjectsTimingClient_cff::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?.
 }

magneticfield::eSector::~eSector ( )

default

Destructor.

Member Function Documentation

MagESector * eSector::buildMagESector ( ) const

Construct the MagESector upon request.

Definition at line 75 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;
 }