#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.
std::vector< MagESector * >	endcapSectors () const
	Get endcap layers.
	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.
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/05/23 22:36:43

Revision:: 1.12

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 ) [private, virtual]

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, Exception, f, ClusterizingHistogram::fill(), first, DDExpandedView::firstChild(), DDExpandedView::geoHistory(), i, prof2calltree::last, DDExpandedView::logicalPart(), MagGeoBuilderFromDDD::volumeHandle::magFile, MagGeoBuilderFromDDD::volumeHandle::masterSector, mBLayers, mESectors, DDName::name(), DDBase< N, C >::name(), mergeVDriftHistosByStation::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 buildInterpolator(), and 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 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;
}