#include <MagGeoBuilderFromDDD.h>

Classes
class	bLayer

class	bRod

class	bSector

class	bSlab

class	eLayer

class	eSector

struct	ExtractAbsZ

struct	ExtractPhi

struct	ExtractPhiMax

struct	ExtractR

struct	ExtractRN

struct	ExtractZ

struct	LessDPhi

struct	LessZ

class	volumeHandle

Public Member Functions
std::vector< MagBLayer * >	barrelLayers () const
	Get barrel layers. More...

std::vector< MagESector * >	endcapSectors () const
	Get endcap layers. More...

	MagGeoBuilderFromDDD (std::string version_, bool debug=false, bool overrideMasterSector=false)

float	maxR () const

float	maxZ () const

void	setScaling (std::vector< int > keys, std::vector< double > values)

virtual	~MagGeoBuilderFromDDD ()
	Destructor. More...

Private Types
typedef std::vector < volumeHandle * >	handles

typedef ConstReferenceCountingPointer < Surface >	RCPS

typedef std::unary_function < const volumeHandle *, double >	uFcn

Private Member Functions
std::vector< MagVolume6Faces * >	barrelVolumes () const

virtual void	build (const DDCompactView &cpv)

void	buildInterpolator (const volumeHandle vol, std::map< std::string, MagProviderInterpol > &interpolators)

void	buildMagVolumes (const handles &volumes, std::map< std::string, MagProviderInterpol * > &interpolators)

std::vector< MagVolume6Faces * >	endcapVolumes () const

void	summary (handles &volumes)

void	testInside (handles &volumes)

Private Attributes
handles	bVolumes

handles	eVolumes

std::vector< MagBLayer * >	mBLayers

std::vector< MagESector * >	mESectors

bool	overrideMasterSector

std::map< int, double >	theScalingFactors

std::string	version

Static Private Attributes
static bool	debug

Friends
class	MagGeometry

class	magneticfield::AutoMagneticFieldESProducer

class	magneticfield::VolumeBasedMagneticFieldESProducer

class	TestMagVolume

Detailed Description

Parse the XML magnetic geometry, build individual volumes and match their shared surfaces. Build MagVolume6Faces and organise them in a hierarchical structure. Build MagGeometry out of it.

Date:: 2009/01/21 12:05:14

Revision:: 1.11

Author: N. Amapane - INFN Torino

Definition at line 32 of file MagGeoBuilderFromDDD.h.

Member Typedef Documentation

typedef std::vector<volumeHandle*> MagGeoBuilderFromDDD::handles

private

Definition at line 78 of file MagGeoBuilderFromDDD.h.

typedef ConstReferenceCountingPointer<Surface> MagGeoBuilderFromDDD::RCPS

private

Definition at line 60 of file MagGeoBuilderFromDDD.h.

typedef std::unary_function<const volumeHandle*, double> MagGeoBuilderFromDDD::uFcn

private

Definition at line 106 of file MagGeoBuilderFromDDD.h.

Constructor & Destructor Documentation

MagGeoBuilderFromDDD::MagGeoBuilderFromDDD	(	std::string	version_,
		bool	debug = `false`,
		bool	overrideMasterSector = `false`
	)

Constructor. overrideMasterSector is a hack to allow switching between phi-symmetric maps and maps with sector-specific tables. It won't be necessary anymore once the geometry is decoupled from the specification of tables, ie when tables will come from the DB.

Definition at line 60 of file MagGeoBuilderFromDDD.cc.

References gather_cfg::cout, and debug.

                                                                                                    :
   version (version_),
   overrideMasterSector(overrideMasterSector_)
 {  
   debug = debug_;
   if (debug) cout << "Constructing a MagGeoBuilderFromDDD" <<endl;
 }

MagGeoBuilderFromDDD::~MagGeoBuilderFromDDD ( )

virtual

Destructor.

Definition at line 68 of file MagGeoBuilderFromDDD.cc.

References bVolumes, eVolumes, and i.

                                            {
   for (handles::const_iterator i=bVolumes.begin();
        i!=bVolumes.end(); ++i){
     delete (*i);
   }
   
   for (handles::const_iterator i=eVolumes.begin();
        i!=eVolumes.end(); ++i){
     delete (*i);
   }
 }

Member Function Documentation

vector< MagBLayer * > MagGeoBuilderFromDDD::barrelLayers ( ) const

Get barrel layers.

Definition at line 640 of file MagGeoBuilderFromDDD.cc.

References mBLayers.

Referenced by magneticfield::VolumeBasedMagneticFieldESProducer::produce().

                                                            {
   return mBLayers;
 }

vector< MagVolume6Faces * > MagGeoBuilderFromDDD::barrelVolumes ( ) const

private

Definition at line 648 of file MagGeoBuilderFromDDD.cc.

References bVolumes, i, and v.

Referenced by magneticfield::VolumeBasedMagneticFieldESProducer::produce().

                                                                   {
   vector<MagVolume6Faces*> v;
   v.reserve(bVolumes.size());
   for (handles::const_iterator i=bVolumes.begin();
        i!=bVolumes.end(); ++i){
     v.push_back((*i)->magVolume);
   }
   return v;
 }

void MagGeoBuilderFromDDD::build ( const DDCompactView & cpv )

privatevirtual

Definition at line 126 of file MagGeoBuilderFromDDD.cc.

References buildInterpolator(), buildMagVolumes(), bVolumes, MagGeoBuilderFromDDD::volumeHandle::center(), ClusterizingHistogram::clusterize(), MagGeoBuilderFromDDD::volumeHandle::copyno, gather_cfg::cout, debug, eVolumes, edm::hlt::Exception, f, ClusterizingHistogram::fill(), first, DDExpandedView::firstChild(), DDExpandedView::geoHistory(), i, prof2calltree::last, DDExpandedView::logicalPart(), MagGeoBuilderFromDDD::volumeHandle::magFile, MagGeoBuilderFromDDD::volumeHandle::masterSector, mBLayers, mESectors, DDName::name(), mergeVDriftHistosByStation::name, DDBase< N, C >::name(), DDExpandedView::nextSibling(), evf::evtn::offset(), overrideMasterSector, PV3DBase< T, PVType, FrameType >::perp(), precomputed_value_sort(), dttmaxenums::R, dtT0WireCalibration_cfg::resolution, DDExpandedView::rotation(), mergeVDriftHistosByStation::sectors, findQualityFiles::size, summary(), testInside(), DDExpandedView::translation(), v, Gflash::Z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by magneticfield::VolumeBasedMagneticFieldESProducer::produce().

 {
 //    DDCompactView cpv;
   DDExpandedView fv(cpva);
 
   if (debug) cout << "**********************************************************" <<endl;
 
   // The actual field interpolators
   map<string, MagProviderInterpol*> bInterpolators;
   map<string, MagProviderInterpol*> eInterpolators;
   
   // Counter of different volumes
   int bVolCount = 0;
   int eVolCount = 0;
 
   if (fv.logicalPart().name().name()!="MAGF") {
      std::string topNodeName(fv.logicalPart().name().name());
 
      //see if one of the children is MAGF
      bool doSubDets = fv.firstChild();
      
      bool go=true;
      while(go&& doSubDets) {
     if (fv.logicalPart().name().name()=="MAGF")
        break;
     else
        go = fv.nextSibling();
      }
      if (!go) {
     throw cms::Exception("NoMAGFinDDD")<<" Neither he top node, nor any child node of the DDCompactView is \"MAGF\" but the top node is instead \""<<topNodeName<<"\"";
      }
   }
   // Loop over MAGF volumes and create volumeHandles. 
   if (debug) { cout << endl << "*** MAGF: " << fv.geoHistory() << endl
             << "translation: " << fv.translation() << endl
             << " rotation: " << fv.rotation() << endl;
   }
   
   bool doSubDets = fv.firstChild();
   while (doSubDets){
     
     string name = fv.logicalPart().name().name();
     if (debug) cout << endl << "Name: " << name << endl
                    << "      " << fv.geoHistory() <<endl;
 
     // FIXME: special handling of version 85l_030919. This version is no 
     //        longer supported and special handling for it may not work
     //        anymore - it will eventually be removed.
     // Build only the z-negative volumes, assuming symmetry
     if (name.substr(2,2)=="ZP") {
       doSubDets = fv.nextSibling();
       continue;
     }
     
     // FIXME: single-volyme cylinders - this is feature has been removed 
     //        and should be revisited.
     //    bool mergeCylinders=false;
 
     // If true, In the barrel, cylinders sectors will be skipped to build full 
     // cylinders out of sector copyno #1.
     bool expand = false;
 
 //     if (mergeCylinders) {
 //       if (name == "V_ZN_1"
 //    || name == "V_ZN_2") {
 //  if (debug && fv.logicalPart().solid().shape()!=ddtubs) {
 //    cout << "ERROR: MagGeoBuilderFromDDD::build: volume " << name
 //         << " should be a cylinder" << endl;
 //  }
 //  if(fv.copyno()==1) {
 //    expand = true;
 //  } else {
 //    //cout << "... to be skipped: "
 //    //     << name << " " << fv.copyno() << endl;
 //  }
 //       }
 //     }
 
     volumeHandle* v = new volumeHandle(fv, expand);
     //FIXME: overrideMasterSector is a hack to allow switching between 
     // phi-symmetric maps and maps with sector-specific tables. 
     // It won't be necessary anymore once the geometry is decoupled from 
     // the specification of tables, ie when tables will come from the DB.
     if (overrideMasterSector && v->masterSector!=1) {
       v->masterSector=1;
       std::string::size_type ibeg, iend;
       ibeg = (v->magFile).rfind('/');
       iend = (v->magFile).size();
       v->magFile = (v->magFile).substr(ibeg+1,iend-ibeg-1);
     }
 
     // Select volumes, build volume handles.
     float Z = v->center().z();
     float R = v->center().perp();
 
     // v 85l: Barrel is everything up to |Z| = 661.0, excluding 
     // volume #7, centered at 6477.5
     // v 1103l: same numbers work fine. #16 instead of #7, same coords;
     // see comment below for V6,7
     //ASSUMPTION: no misalignment is applied to mag volumes.
     //FIXME: implement barrel/endcap flags as DDD SpecPars.
     if ((fabs(Z)<647. || (R>350. && fabs(Z)<662.)) &&
     !(fabs(Z)>480 && R<172) // in 1103l we place V_6 and V_7 in the 
                             // endcaps to preserve nice layer structure
                             // in the barrel. This does not hurt in v85l
                             // where there is a single V1 
     ) { // Barrel
       if (debug) cout << " (Barrel)" <<endl;
       bVolumes.push_back(v);
 
 
       // Build the interpolator of the "master" volume (the one which is
       // not replicated in phi)
       // ASSUMPTION: copyno == sector. 
       //             Note this is not the case for obsolete grid_85l_030919 - 
       // In case in  future models this should no more be the case, can 
       // implement secotor numbers as SpecPars in the XML      
       if (v->copyno==v->masterSector) {
     buildInterpolator(v, bInterpolators);
     ++bVolCount;
       }
     } else {               // Endcaps
       if (debug) cout << " (Endcaps)" <<endl;
       eVolumes.push_back(v);
       if (v->copyno==v->masterSector) { 
     buildInterpolator(v, eInterpolators);
     ++eVolCount;
       }
     }
 
     doSubDets = fv.nextSibling(); // end of loop over MAGF
   }
     
   if (debug) {
     cout << "Number of volumes (barrel): " << bVolumes.size() <<endl
           << "Number of volumes (endcap): " << eVolumes.size() <<endl;
     cout << "**********************************************************" <<endl;
   }
 
   // Now all volumeHandles are there, and parameters for each of the planes
   // are calculated.
 
   //----------------------------------------------------------------------
   // Print summary information
 
   if (debug) {
     cout << "-----------------------" << endl;
     cout << "SUMMARY: Barrel " << endl;
     summary(bVolumes);
     
     cout << endl << "SUMMARY: Endcaps " << endl;
     summary(eVolumes);
     cout << "-----------------------" << endl;
   }
 
 
   //----------------------------------------------------------------------
   // Find barrel layers.
 
   vector<bLayer> layers; // the barrel layers
   precomputed_value_sort(bVolumes.begin(), bVolumes.end(), ExtractRN());
 
   // Find the layers (in R)
   const float resolution = 1.; // cm
   float rmin = bVolumes.front()->RN()-resolution;
   float rmax = bVolumes.back()->RN()+resolution;
   ClusterizingHistogram  hisR( int((rmax-rmin)/resolution) + 1, rmin, rmax);
 
   if (debug) cout << " R layers: " << rmin << " " << rmax << endl;
 
   handles::const_iterator first = bVolumes.begin();
   handles::const_iterator last = bVolumes.end();  
 
   for (handles::const_iterator i=first; i!=last; ++i){
     hisR.fill((*i)->RN());
   }
   vector<float> rClust = hisR.clusterize(resolution);
 
   handles::const_iterator ringStart = first;
   handles::const_iterator separ = first;
 
   for (unsigned int i=0; i<rClust.size() - 1; ++i) {
     if (debug) cout << " Layer at RN = " << rClust[i];
     float rSepar = (rClust[i] + rClust[i+1])/2.f;
     while ((*separ)->RN() < rSepar) ++separ;
 
     bLayer thislayer(ringStart, separ);
     layers.push_back(thislayer);
     ringStart = separ;
   }
   {
     if (debug) cout << " Layer at RN = " << rClust.back();
     bLayer thislayer(separ, last);
     layers.push_back(thislayer);
   }
 
   if (debug) cout << "Barrel: Found " << rClust.size() << " clusters in R, "
           << layers.size() << " layers " << endl << endl;
 
 
   //----------------------------------------------------------------------
   // Find endcap sectors
 
   vector<eSector> sectors; // the endcap sectors
   precomputed_value_sort(eVolumes.begin(), eVolumes.end(), ExtractPhi()); 
  
 
   // ASSUMPTION: There are 12 sectors and each sector is 30 deg wide.
   for (int i = 0; i<12; ++i) {
     int offset = eVolumes.size()/12;
     //    int isec = (i+binOffset)%12;
     if (debug) cout << " Sector at phi = "
             << (*(eVolumes.begin()+((i)*offset)))->center().phi()
             << endl;
     sectors.push_back(eSector(eVolumes.begin()+((i)*offset),
                   eVolumes.begin()+((i+1)*offset)));
   }
    
   if (debug) cout << "Endcap: Found " 
           << sectors.size() << " sectors " << endl;
 
 
   //----------------------------------------------------------------------  
   // Match surfaces.
 
 //  cout << "------------------" << endl << "Now associating planes..." << endl;
 
 //   // Loop on layers
 //   for (vector<bLayer>::const_iterator ilay = layers.begin();
 //        ilay!= layers.end(); ++ilay) {
 //     cout << "On Layer: " << ilay-layers.begin() << " RN: " << (*ilay).RN()
 //   <<endl;     
 
 //     // Loop on wheels
 //     for (vector<bWheel>::const_iterator iwheel = (*ilay).wheels.begin();
 //   iwheel != (*ilay).wheels.end(); ++iwheel) {
 //       cout << "  On Wheel: " << iwheel- (*ilay).wheels.begin()<< " Z: "
 //     << (*iwheel).minZ() << " " << (*iwheel).maxZ() << " " 
 //     << ((*iwheel).minZ()+(*iwheel).maxZ())/2. <<endl;
 
 //       // Loop on sectors.
 //       for (int isector = 0; isector<12; ++isector) {
 //  // FIXME: create new constructor...
 //  bSectorNavigator navy(layers,
 //                ilay-layers.begin(),
 //                iwheel-(*ilay).wheels.begin(),isector);
     
 //  const bSector & isect = (*iwheel).sector(isector);
     
 //  isect.matchPlanes(navy); //FIXME refcount
 //       }
 //     }
 //   }
 
 
   //----------------------------------------------------------------------
   // Build MagVolumes and the MagGeometry hierarchy.
 
   //--- Barrel
 
   // Build MagVolumes and associate interpolators to them
   buildMagVolumes(bVolumes, bInterpolators);
 
   // Build MagBLayers
   for (vector<bLayer>::const_iterator ilay = layers.begin();
        ilay!= layers.end(); ++ilay) {
     mBLayers.push_back((*ilay).buildMagBLayer());
   }
 
   if (debug) {  
     cout << "*** BARREL ********************************************" << endl
      << "Number of different volumes   = " << bVolCount << endl
      << "Number of interpolators built = " << bInterpolators.size() << endl
          << "Number of MagBLayers built    = " << mBLayers.size() << endl;
 
     testInside(bVolumes); // FIXME: all volumes should be checked in one go.
   }
   
   //--- Endcap
   // Build MagVolumes  and associate interpolators to them
   buildMagVolumes(eVolumes, eInterpolators);
 
   // Build the MagESectors
   for (vector<eSector>::const_iterator isec = sectors.begin();
        isec!= sectors.end(); ++isec) {
     mESectors.push_back((*isec).buildMagESector());
   }
 
   if (debug) {
     cout << "*** ENDCAP ********************************************" << endl
      << "Number of different volumes   = " << eVolCount << endl
      << "Number of interpolators built = " << eInterpolators.size() << endl
          << "Number of MagESector built    = " << mESectors.size() << endl;
 
     testInside(eVolumes); // FIXME: all volumes should be checked in one go.
   }
 }

void MagGeoBuilderFromDDD::buildInterpolator	(	const volumeHandle *	vol,
		std::map< std::string, MagProviderInterpol * > &	interpolators
	)

private

Definition at line 487 of file MagGeoBuilderFromDDD.cc.

Referenced by build().

                                                                                                                        {
 
 
   // FIXME: special handling of version 85l_030919. This version is no 
   //        longer supported and special handling for it may not work
   //        anymore - it will eventually be removed.
   // In version grid_85l_030919, interpolators should be built only 
   // for volumes on NEGATIVE z 
   // (Z symmetry in field tables)
   if (version=="grid_85l_030919" && vol->center().z()>0) return;
 
   // Phi of the master sector
   double masterSectorPhi = (vol->masterSector-1)*Geom::pi()/6.;
   //FIXME: special handling of version 85l_030919 (see FIXME above). 
   // In ver. grid_85l_030919, the master sector was sector 4 
   // (along Y axis).
   if (version=="grid_85l_030919") {
     masterSectorPhi=Geom::pi()/2.;
   } 
 
   if (debug) {
     cout << "Building interpolator from "
      << vol->name << " copyno " << vol->copyno
      << " at " << vol->center()
      << " phi: " << vol->center().phi()
      << " file: " << vol->magFile
      << endl;
 
     if ( fabs(vol->center().phi() - masterSectorPhi) > Geom::pi()/9.) {
       cout << "***WARNING wrong sector? " << endl;
     }
   }
 
   if (version == "fake") {
     interpolators[vol->magFile] = new magneticfield::FakeInterpolator();
     return;
   }
 
   string fullPath;
 
   try {
     edm::FileInPath mydata("MagneticField/Interpolation/data/"+version+"/"+vol->magFile);
     fullPath = mydata.fullPath();
   } catch (edm::Exception& exc) {
     cerr << "MagGeoBuilderFromDDD: exception in reading table; " << exc.what() << endl;
     if (!debug) throw;
     return;
   }
   
   
   try{
     if (vol->toExpand()){
       //FIXME: see discussion on mergeCylinders above.
 //       interpolators[vol->magFile] =
 //  MFGridFactory::build( fullPath, *(vol->placement()), vol->minPhi(), vol->maxPhi());
     } else {
       // If the table is in "local" coordinates, must create a reference 
       // frame that is appropriately rotated along the CMS Z axis.
 
       GloballyPositioned<float> rf = *(vol->placement());
 
       if (vol->masterSector != 1) {
     typedef Basic3DVector<float> Vector;
 
     GloballyPositioned<float>::RotationType rot(Vector(0,0,1), -masterSectorPhi);
     Vector vpos(vol->placement()->position());
     
 
     rf = GloballyPositioned<float>(GloballyPositioned<float>::PositionType(rot.multiplyInverse(vpos)), vol->placement()->rotation()*rot);
       }
 
       interpolators[vol->magFile] =
     MFGridFactory::build( fullPath, rf);
     }
   } catch (MagException& exc) {
     cout << exc.what() << endl;
     interpolators.erase(vol->magFile);
     if (!debug) throw; 
     return;
   }
 
 
     if (debug) {
     // Check that all grid points of the interpolator are inside the volume.
       const MagVolume6Faces tempVolume(vol->placement()->position(),
                  vol->placement()->rotation(),
                  vol->shape(),
                  vol->sides(), 
                  interpolators[vol->magFile]);
 
       const MFGrid3D* grid = dynamic_cast<const MFGrid3D*>(interpolators[vol->magFile]);
       if (grid!=0) {
     
     Dimensions sizes = grid->dimensions();
     cout << "Grid has 3 dimensions " 
          << " number of nodes is " << sizes.w << " " << sizes.h
          << " " << sizes.d << endl;
       
     const double tolerance = 0.03;
 
 
     int dumpCount = 0;
     for (int j=0; j < sizes.h; j++) {
       for (int k=0; k < sizes.d; k++) {
         for (int i=0; i < sizes.w; i++) {
           MFGrid::LocalPoint lp = grid->nodePosition( i, j, k);
           if (! tempVolume.inside(lp, tolerance)) {
         if (++dumpCount < 2) {
           MFGrid::GlobalPoint gp = tempVolume.toGlobal(lp);
           cout << "GRID ERROR: " << i << " " << j << " " << k
                << " local: " << lp
                << " global: " << gp
                << " R= " << gp.perp() << " phi=" << gp.phi() << endl;
         }
           }
         }
       }
     }
     
     cout << vol->name << " : Number of grid points outside the MagVolume: " << dumpCount << "/" << sizes.w*sizes.h*sizes.d << endl;
       }
     }
 }

void MagGeoBuilderFromDDD::buildMagVolumes	(	const handles &	volumes,
		std::map< std::string, MagProviderInterpol * > &	interpolators
	)

private

Definition at line 425 of file MagGeoBuilderFromDDD.cc.

References gather_cfg::cout, combine::key, mergeVDriftHistosByStation::name, MagVolume::ownsFieldProvider(), GloballyPositioned< T >::position(), GloballyPositioned< T >::rotation(), and theScalingFactors.

Referenced by build().

                                                                                                                      {
   // Build all MagVolumes setting the MagProviderInterpol
   for (handles::const_iterator vol=volumes.begin(); vol!=volumes.end();
        ++vol){
     const MagProviderInterpol* mp = 0;
     if (interpolators.find((*vol)->magFile)!=interpolators.end()) {
       mp = interpolators[(*vol)->magFile];
     } else {
       cout << "No interpolator found for file " << (*vol)->magFile
        << " vol: " << (*vol)->name << endl;
       cout << interpolators.size() <<endl;
     }
       
 
     // Get the volume number from the volume name.
     // ASSUMPTION: volume names ends with "_NUM" where NUM is the volume number
     int volNum;
     string name = (*vol)->name;
     name.erase(0,name.rfind('_')+1);
     stringstream str;    
     str << name;
     str >> volNum;
 
 
     // ASSUMPTION: copyno == sector. 
     //             Note this is not the case for obsolete grid_85l_030919 - 
     // In case in  future models this should no more be the case, can 
     // implement secotor numbers as SpecPars in the XML
     int key = volNum*100+(*vol)->copyno; 
     map<int, double>::const_iterator isf = theScalingFactors.find(key);
     if (isf == theScalingFactors.end()) {
       key = volNum*100;
       isf = theScalingFactors.find(key);
     }
     
     double sf = 1.;
     if (isf != theScalingFactors.end()) {
       sf = (*isf).second;
 
       edm::LogInfo("MagneticField|VolumeBasedMagneticFieldESProducer") << "Applying scaling factor " << sf << " to "<< (*vol)->name << "["<< (*vol)->copyno << "] (key:" << key << ")" << endl;
     }
 
     const GloballyPositioned<float> * gpos = (*vol)->placement();
     (*vol)->magVolume = new MagVolume6Faces(gpos->position(),
                         gpos->rotation(),
                         (*vol)->shape(),
                         (*vol)->sides(),
                         mp, sf);
 
     if ((*vol)->copyno==(*vol)->masterSector) {
       (*vol)->magVolume->ownsFieldProvider(true);
     }
 
     (*vol)->magVolume->setIsIron((*vol)->isIron());
 
     // The name and sector of the volume are saved for debug purposes only. They may be removed at some point...
     (*vol)->magVolume->name = (*vol)->name;
     (*vol)->magVolume->copyno = (*vol)->copyno;
   }
 }

vector< MagESector * > MagGeoBuilderFromDDD::endcapSectors ( ) const

Get endcap layers.

Definition at line 644 of file MagGeoBuilderFromDDD.cc.

References mESectors.

Referenced by magneticfield::VolumeBasedMagneticFieldESProducer::produce().

                                                              {
   return mESectors;
 }

vector< MagVolume6Faces * > MagGeoBuilderFromDDD::endcapVolumes ( ) const

private

Definition at line 658 of file MagGeoBuilderFromDDD.cc.

References eVolumes, i, and v.

Referenced by magneticfield::VolumeBasedMagneticFieldESProducer::produce().

                                                                   {
   vector<MagVolume6Faces*> v;
   v.reserve(eVolumes.size());
   for (handles::const_iterator i=eVolumes.begin();
        i!=eVolumes.end(); ++i){
     v.push_back((*i)->magVolume);
   }
   return v;
 }

float MagGeoBuilderFromDDD::maxR ( ) const

Definition at line 669 of file MagGeoBuilderFromDDD.cc.

References version.

Referenced by magneticfield::VolumeBasedMagneticFieldESProducer::produce().

                                       {
   //FIXME: should get it from the actual geometry - MAGF is an option, 
   //       but that is not changed together the geometry itself 
   //       (it lives in cmsMagneticField.xml in CMSCommonData)
   if (version=="grid_85l_030919") return 1000.;
   else return 900.;
 }

float MagGeoBuilderFromDDD::maxZ ( ) const

Definition at line 677 of file MagGeoBuilderFromDDD.cc.

Referenced by MagGeoBuilderFromDDD::bSector::bSector(), and magneticfield::VolumeBasedMagneticFieldESProducer::produce().

                                       {
   //FIXME: should get it from the actual geometry - see above
   return 1600.;
 }

void MagGeoBuilderFromDDD::setScaling	(	std::vector< int >	keys,
		std::vector< double >	values
	)

Set scaling factors for individual volumes. "keys" is a vector of 100*volume number + sector (sector 0 = all sectors) "values" are the corresponding scaling factors

Definition at line 683 of file MagGeoBuilderFromDDD.cc.

References edm::hlt::Exception, i, and theScalingFactors.

Referenced by magneticfield::VolumeBasedMagneticFieldESProducer::produce().

 {
   if (keys.size() != values.size()) {
     throw cms::Exception("InvalidParameter") << "Invalid field scaling parameters 'scalingVolumes' and 'scalingFactors' ";
   }
   for (unsigned int i=0; i<keys.size(); ++i) {
     theScalingFactors[keys[i]] = values[i];
   }
 }

void MagGeoBuilderFromDDD::summary ( handles & volumes )

private

Definition at line 81 of file MagGeoBuilderFromDDD.cc.

References gather_cfg::cout, first, i, prof2calltree::last, and edm::second().

Referenced by build().

                                                    {  
   // The final countdown.
   int ivolumes  = volumes.size();  // number of volumes
   int isurfaces = ivolumes*6;       // number of individual surfaces
   int iassigned = 0;                // How many have been assigned
   int iunique   = 0;                // number of unique surfaces
   int iref_ass  = 0;
   int iref_nass = 0;
 
   set<const void *> ptrs;
 
   handles::const_iterator first = volumes.begin();
   handles::const_iterator last = volumes.end();
 
   for (handles::const_iterator i=first; i!=last; ++i){
     if (int((*i)->shape())>4) continue; // FIXME: implement test for missing shapes...
     for (int side = 0; side < 6; ++side) {
       int references =  (*i)->references(side);
       if ((*i)->isPlaneMatched(side)) {
     ++iassigned;
     bool firstOcc = (ptrs.insert(&((*i)->surface(side)))).second;
     if (firstOcc) iref_ass+=references;
     if (references<2){  
       cout << "*** Only 1 ref, vol: " << (*i)->name << " # "
            << (*i)->copyno << " side: " << side << endl;
     }   
       } else {
     iref_nass+=references;
     if (references>1){
       cout << "*** Ref_nass >1 " <<endl;
     }
       }
     }
   }
   iunique = ptrs.size();
 
   cout << "    volumes   " << ivolumes  << endl
        << "    surfaces  " << isurfaces << endl
        << "    assigned  " << iassigned << endl
        << "    unique    " << iunique << endl
        << "    iref_ass  " << iref_ass << endl
        << "    iref_nass " << iref_nass << endl;
 }

void MagGeoBuilderFromDDD::testInside ( handles & volumes )

private

Definition at line 613 of file MagGeoBuilderFromDDD.cc.

References gather_cfg::cout, and i.

Referenced by build().

                                                        {
   // test inside() for all volumes.
   cout << "--------------------------------------------------" << endl;
   cout << " inside(center) test" << endl;
   for (handles::const_iterator vol=volumes.begin(); vol!=volumes.end();
        ++vol){
     for (handles::const_iterator i=volumes.begin(); i!=volumes.end();
      ++i){
       if ((*i)==(*vol)) continue;
       //if ((*i)->magVolume == 0) continue;
       if ((*i)->magVolume->inside((*vol)->center())) {
     cout << "*** ERROR: center of " << (*vol)->name << " is inside " 
          << (*i)->name <<endl;
       }
     }
     
     if ((*vol)->magVolume->inside((*vol)->center())) {
       cout << (*vol)->name << " OK " << endl;
     } else {
       cout << "*** ERROR: center of volume is not inside it, "
        << (*vol)->name << endl;
     }
   }
   cout << "--------------------------------------------------" << endl;
 }