#include <bSector.h>

Public Member Functions
	bSector ()
	Default ctor is needed to have arrays.
	bSector (handles::const_iterator begin, handles::const_iterator end)
	Constructor from list of volumes.
MagBSector *	buildMagBSector () const
	Construct the MagBSector upon request.
const handles &	getVolumes () const
	Return all volumes in this sector.
const float	RN () const
	Distance from center along normal of sectors.
	~bSector ()
	Destructor.
Private Attributes
MagBSector *	msector
std::vector< bRod >	rods
handles	volumes

Detailed Description

A sector of volumes in a barrel layer (i.e. only 1 element in R) One sector is composed of 1 or more rods.

Date:: 2005/09/27 15:15:52

Revision:: 1.2

Author:: N. Amapane - INFN Torino

Definition at line 19 of file bSector.h.

Constructor & Destructor Documentation

MagGeoBuilderFromDDD::bSector::bSector ( )

Default ctor is needed to have arrays.

Definition at line 25 of file bSector.cc.

{}

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

Constructor from list of volumes.

FIXME: (float) resolution; ??

Definition at line 30 of file bSector.cc.

References ClusterizingHistogram::clusterize(), gather_cfg::cout, MagGeoBuilderFromDDD::debug, f, ClusterizingHistogram::fill(), first, i, prof2calltree::last, MagGeoBuilderFromDDD::maxZ(), cmsRun_displayProdMFGeom_cfg::minZ, jptDQMConfig_cff::phiMax, jptDQMConfig_cff::phiMin, precomputed_value_sort(), MagGeoBuilderFromDDD::volumeHandle::printUniqueNames(), dtDQMClient_cfg::resolution, rods, and volumes.

                                                                   :
  volumes(begin,end),
  msector(0)
{
  if (MagGeoBuilderFromDDD::debug) cout << "   Sector at Phi  " <<  volumes.front()->center().phi() << " " 
                  << volumes.back()->center().phi() <<  endl;

  if (volumes.size() == 1) {
    if (MagGeoBuilderFromDDD::debug) { 
      cout << "   Rod at: 0 elements: " << end-begin
           << " unique volumes: ";
      volumeHandle::printUniqueNames(begin,end);
    }
    rods.push_back(bRod(begin,end));
  } else {
    // Clusterize in phi. Use bin edge so that complete clusters can be 
    // easily found (not trivial using bin centers!)
    // Unfortunately this makes the result more sensitive to the 
    // "resolution" parameter...
    // To avoid +-pi boundary problems, take phi distance from min phi.
    // Caveat of implicit conversions of Geom::Phi!!!

    // Sort volumes in DELTA phi - i.e. phi(j)-phi(i) > 0 if j>1.
    precomputed_value_sort(volumes.begin(), volumes.end(),
                           ExtractPhiMax(), LessDPhi());

    const Geom::Phi<float> resolution(0.01); // rad
    Geom::Phi<float> phi0 = volumes.front()->maxPhi();
    float phiMin = -(float) resolution;
    float phiMax = volumes.back()->maxPhi() - phi0 + resolution; 

    ClusterizingHistogram hisPhi( int((phiMax-phiMin)/resolution) + 1,
                                  phiMin, phiMax);
    
    handles::const_iterator first = volumes.begin();
    handles::const_iterator last = volumes.end();  

    for (handles::const_iterator i=first; i!=last; ++i){
      hisPhi.fill((*i)->maxPhi()-phi0);
    }
    vector<float> phiClust = hisPhi.clusterize(resolution);

    if (MagGeoBuilderFromDDD::debug) cout << "     Found " << phiClust.size() << " clusters in Phi, "
                    << " rods: " << endl;

    handles::const_iterator rodStart = first;
    handles::const_iterator separ = first;
    
    float DZ = (*max_element(first,last,LessZ()))->maxZ() -
      (*min_element(first,last,LessZ()))->minZ();    

    float DZ1 = 0.;
    for (unsigned int i=0; i<phiClust.size(); ++i) {
      float phiSepar;
      if (i<phiClust.size()-1) {
        phiSepar = (phiClust[i] + phiClust[i+1])/2.f;
      } else {
        phiSepar = phiMax;
      }
      if (MagGeoBuilderFromDDD::debug) cout << "       cluster " << i
                      << " phisepar " << phiSepar <<endl;
      while (separ < last && (*separ)->maxPhi()-phi0 < phiSepar ) {
        DZ1 += ((*separ)->maxZ() - (*separ)->minZ());
        if (MagGeoBuilderFromDDD::debug) cout << "         " << (*separ)->name << " "
                        << (*separ)->maxPhi()-phi0  << " "
                        << (*separ)->maxZ() << " " << (*separ)->minZ() << " "
                        << DZ1 << endl;
        ++separ;
      }

      // FIXME: print warning for small discrepancies. Tolerance (below)
      // had to be increased since discrepancies sum to up to ~ 2 mm.
      if (fabs(DZ-DZ1) > 0.001 && fabs(DZ-DZ1) < 0.5) {
        if (MagGeoBuilderFromDDD::debug) cout << "*** WARNING: Z lenght mismatch by " << DZ-DZ1
                        << " " << DZ << " " << DZ1 << endl;

      }
      if (fabs(DZ-DZ1) > 0.25 ) { // FIXME hardcoded tolerance
        if (MagGeoBuilderFromDDD::debug) cout << "       Incomplete, use also next cluster: " 
                        << DZ << " " << DZ1 << " " << DZ-DZ1 << endl;
        DZ1 = 0.;
        continue;
      } else if (DZ1>DZ+0.05) { // Wrong: went past max lenght // FIXME hardcoded tolerance
        cout << " *** ERROR: bSector finding messed up." << endl;
        volumeHandle::printUniqueNames(rodStart, separ);
        DZ1 = 0.;
      } else {
        if (MagGeoBuilderFromDDD::debug) {
          cout << "       Rod at: " << phiClust[i] <<" elements: "
               << separ-rodStart << " unique volumes: ";
          volumeHandle::printUniqueNames(rodStart, separ);
        }
        
        rods.push_back(bRod(rodStart, separ));
        rodStart = separ;
        DZ1 = 0.;
      }
    }

    if (MagGeoBuilderFromDDD::debug) cout << "-----------------------" << endl;

  }
}

MagGeoBuilderFromDDD::bSector::~bSector ( )

Destructor.

Definition at line 27 of file bSector.cc.

{}

Member Function Documentation

MagBSector * MagGeoBuilderFromDDD::bSector::buildMagBSector ( ) const

Construct the MagBSector upon request.

Definition at line 136 of file bSector.cc.

                                                              {
  if (msector==0) {
    vector<MagBRod*> mRods;
    for (vector<bRod>::const_iterator rod = rods.begin();
         rod!=rods.end(); ++rod) {
      mRods.push_back((*rod).buildMagBRod());
    }
    msector = new MagBSector(mRods, volumes.front()->minPhi()); //FIXME
  }
  return msector;
}

const handles& MagGeoBuilderFromDDD::bSector::getVolumes ( ) const [inline]

Return all volumes in this sector.

Definition at line 36 of file bSector.h.

References volumes.

{return volumes;}

const float MagGeoBuilderFromDDD::bSector::RN ( ) const [inline]

Distance from center along normal of sectors.

Definition at line 31 of file bSector.h.

References volumes.

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

Member Data Documentation

MagBSector* MagGeoBuilderFromDDD::bSector::msector [mutable, private]

Definition at line 44 of file bSector.h.

std::vector<bRod> MagGeoBuilderFromDDD::bSector::rods [private]

Definition at line 42 of file bSector.h.

Referenced by bSector().

handles MagGeoBuilderFromDDD::bSector::volumes [private]

Definition at line 43 of file bSector.h.

Referenced by bSector(), getVolumes(), and RN().

MagGeoBuilderFromDDD::bSector Class Reference

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation