MagGeoBuilderFromDDD::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)
	Constructor from list of volumes. More...
	~eSector ()
	Destructor. More...

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

Detailed Description

A sector of volumes in the endcap. One sector is composed of several layers (eLayer)

Author

N. Amapane - INFN Torino

Definition at line 17 of file eSector.h.

Constructor & Destructor Documentation

MagGeoBuilderFromDDD::eSector::eSector	(	handles::const_iterator	begin,
		handles::const_iterator	end
	)

Constructor from list of volumes.

Definition at line 21 of file eSector.cc.

References ClusterizingHistogram::clusterize(), gather_cfg::cout, MagGeoBuilderFromDDD::debug, f, ClusterizingHistogram::fill(), first, i, prof2calltree::last, layers, precomputed_value_sort(), MagGeoBuilderFromDDD::volumeHandle::printUniqueNames(), dtDQMClient_cfg::resolution, theVolumes, SiStripMonitorClusterAlca_cfi::zmax, and SiStripMonitorClusterAlca_cfi::zmin.

                                                  :
  theVolumes(begin,end),
  msector(0)
{
  //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 (MagGeoBuilderFromDDD::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 (MagGeoBuilderFromDDD::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 (MagGeoBuilderFromDDD::debug) {
      cout << "     Layer at: " << zClust[i]
       << " elements: " << separ-layStart << " unique volumes: ";
      volumeHandle::printUniqueNames(layStart, separ);
    }
    
    layers.push_back(eLayer(layStart, separ));
    layStart = separ;
  }
  {
    if (MagGeoBuilderFromDDD::debug) {
      cout << "     Layer at: " << zClust.back() <<" elements: " << last-separ
       << " unique volumes: ";
      volumeHandle::printUniqueNames(separ,last);
    }
    layers.push_back(eLayer(separ, last));
  }

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

MagGeoBuilderFromDDD::eSector::~eSector

(

)

Destructor.

Definition at line 79 of file eSector.cc.

{}

Member Function Documentation

MagESector * MagGeoBuilderFromDDD::eSector::buildMagESector

(

)

const

Construct the MagESector upon request.

Definition at line 82 of file eSector.cc.

References LayerTriplets::layers().

                                                              {
  if (msector==0) {
    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

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

private

Definition at line 32 of file eSector.h.

Referenced by eSector().

MagESector* MagGeoBuilderFromDDD::eSector::msector

mutableprivate

Definition at line 34 of file eSector.h.

handles MagGeoBuilderFromDDD::eSector::theVolumes

private

Definition at line 33 of file eSector.h.

Referenced by eSector().