/data/refman/pasoursint/CMSSW_5_3_10_patch1/src/DetectorDescription/Core/src/Specific.cc

Go to the documentation of this file.

#include "DetectorDescription/Core/src/Specific.h"
#include "DetectorDescription/Core/interface/DDCompactView.h"

// Message logger.
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <assert.h>

namespace DDI {
   
   Specific::Specific(const std::vector<std::string>& selections,
                      const DDsvalues_type & specs,
                      bool doRegex)
   : specifics_(specs), 
   partSelections_(0), 
   valid_(false),
   doRegex_(doRegex)                    
   {
      std::vector<std::string>::const_iterator it = selections.begin();
      for(; it != selections.end(); ++it) {
         createPartSelections(*it);
      }
   }
   
   Specific::Specific(const std::vector<DDPartSelection> & selections,
                      const DDsvalues_type & specs)
   : specifics_(specs), partSelections_(selections), valid_(false), doRegex_(false)        
   { }             
   
   void Specific::createPartSelections(const std::string & selString)
   {
      
      std::vector<DDPartSelRegExpLevel> regv;
      std::vector<DDPartSelection> temp;
      DDTokenize2(selString,regv);
      
      if (!regv.size()) throw cms::Exception("DDException") << "Could not evaluate the selection-std::string ->" << selString << "<-";
      std::vector<DDPartSelRegExpLevel>::const_iterator it = regv.begin();
      std::pair<bool,std::string> res;
      for (; it != regv.end(); ++it) {
         std::vector<DDLogicalPart> lpv;
         res = DDIsValid(it->ns_,it->nm_,lpv,doRegex_);
         if (!res.first) {
            std::string msg("Could not process q-name of a DDLogicalPart, reason:\n"+res.second);
            msg+="\nSpecPar selection is:\n" + selString + "\n";
            //throw cms::Exception("DDException") << "Could not process q-name of a DDLogicalPart, reason:\n" << res.second;
            edm::LogError("Specific") << msg;
            break; //EXIT for loop
         }
         //edm::LogInfo ("Specific") << "call addSelectionLevel" << std::endl;
         addSelectionLevel(lpv,it->copyno_,it->selectionType_,temp);
      }
      if ( res.first ) { // i.e. it wasn't "thrown" out of the loop
         std::vector<DDPartSelection>::const_iterator iit = temp.begin();
         partSelections_.reserve(temp.size() + partSelections_.size());
         for (; iit != temp.end(); ++iit) {
            partSelections_.push_back(*iit);
            //edm::LogInfo ("Specific") << *iit << std::endl;
         } 
      }
   }
   
   
   
   void Specific::addSelectionLevel(std::vector<DDLogicalPart> & lpv, int copyno, ddselection_type st, 
                                    std::vector<DDPartSelection> & selv)
   {
      //static int count =0;
      //++count;
      //edm::LogInfo ("Specific") << "count=" << count << " " << flush;
      if (!selv.size()) { // create one, no entry yet!
         selv.push_back(DDPartSelection());
      }
      typedef std::vector<DDLogicalPart>::size_type lpv_sizetype;
      typedef std::vector<DDPartSelection>::size_type ps_sizetype;
      ps_sizetype ps_sz = selv.size();
      lpv_sizetype lpv_sz = lpv.size();
      //edm::LogInfo ("Specific") << "lpv_sz=" << lpv_sz << std::endl;
      lpv_sizetype lpv_i = 0;
      std::vector<DDPartSelection> result;
      for (; lpv_i < lpv_sz; ++lpv_i) {
         std::vector<DDPartSelection>::const_iterator ps_it = selv.begin();
         for (; ps_it != selv.end(); ++ps_it) {
            result.push_back(*ps_it);
         }
      }
      //edm::LogInfo ("Specific") << "result-size=" << result.size() << std::endl;
      //ps_sizetype ps_sz = result.size();
      ps_sizetype ps_i = 0;
      for(lpv_i=0; lpv_i < lpv_sz; ++lpv_i) {
         for(ps_i = ps_sz*lpv_i; ps_i < ps_sz*(lpv_i+1); ++ps_i) {
            result[ps_i].push_back(DDPartSelectionLevel(lpv[lpv_i],copyno,st));
         }
      }    
      selv = result;
   }                   
   
   const std::vector<DDPartSelection> & Specific::selection() const 
   {
      return partSelections_; 
   }  
   
   
   void Specific::stream(std::ostream &  os) const
   {
      //  os << " no output available yet, sorry. ";
      os << " Size: " << specifics_.size() << std::endl;
      os << "\tSelections:" << std::endl;
      partsel_type::const_iterator pit(partSelections_.begin()), pend(partSelections_.end());
      for (;pit!=pend;++pit) {
         os << *pit << std::endl;
      }
      
      DDsvalues_type::const_iterator vit(specifics_.begin()), ved(specifics_.end());
      for (;vit!=ved;++vit) {
         const DDValue & v = vit->second;
         os << "\tParameter name= \"" << v.name() << "\" " << std::endl;
         os << "\t\t Value pairs: " << std::endl;
         size_t s=v.size();
         size_t i=0;
         if ( v.isEvaluated() ) {
            for (; i<s; ++i) {
               os << "\t\t\t\"" << v[i].first << "\"" << ", " << v[i].second << std::endl;
            }
         } else { // v is not evaluated
            const std::vector<std::string>& vs =  v.strings();
            for (; i<s; ++i) {
               os << "\t\t\t\"" << vs[i] << "\"" << ", not evaluated" << std::endl;
            }
         }
      } 
   }
   
   void Specific::updateLogicalPart(std::vector<std::pair<DDLogicalPart, std::pair<DDPartSelection*,DDsvalues_type*> > >& result) const
   {
      if (partSelections_.size()) {
         partsel_type::const_iterator it = partSelections_.begin();
         DDsvalues_type* sv = const_cast<DDsvalues_type*>(&specifics_);
         for (; it != partSelections_.end(); ++it) {
            DDLogicalPart logp = it->back().lp_; 
            /*if (!logp.isDefined().second) {
             throw DDException("Specific::updateLogicalPart(..): LogicalPart not defined, name=" + std::string(logp.ddname()));
             }*/
            DDPartSelection * ps = const_cast<DDPartSelection*>(&(*it));
            assert(ps); 
            assert(sv);
            std::pair<DDPartSelection*,DDsvalues_type*> pssv(ps,sv);
            result.push_back(std::make_pair(logp,pssv));
         }  
      }  
   }
   
   std::pair<bool,DDExpandedView> Specific::node() const
   {
      DDCompactView c;
      DDExpandedView e(c);
      
      if (partSelections_.size() != 1) {
         edm::LogError("Specific") << " >> more or less than one part-selector, currently NOT SUPPORTED! <<" << std::endl;
         return std::make_pair(false,e);
      }
      const DDPartSelection & ps = partSelections_[0];
      
      DDPartSelection::const_iterator it = ps.begin();
      DDPartSelection::const_iterator ed = ps.end();
      if ( (it != ed) && ( it->selectionType_ != ddanyposp) ) {
         edm::LogError("Specific") << " >> part-selector must start with //Name[no] ! << " << std::endl;
         return std::make_pair(false,e);
      }
      ++it;
      for (; it != ps.end(); ++it) {
         if ( it->selectionType_ != ddchildposp ) {
            edm::LogError("Specific") << " >> part-selector must be a concatenation of direct children\n"
            << "   including their copy-number only, CURRENT LIMITATION! <<" << std::endl;
            return std::make_pair(false,e);             
         }  
      } 
      
      it = ps.begin();
      bool result = true;
      for (; it != ed; ++it) {
         while(result) {
            if( (it->copyno_ == e.copyno()) && (it->lp_ == e.logicalPart())) {
               break;
            }
            else {
               result = e.nextSibling();
            }
         }
         if ((ed-it)>1) {
            result = e.firstChild();  
         }  
      }
      return std::make_pair(result,e);
   }
   
   
   Specific::~Specific()
   {
      //   DDsvalues_type::iterator it = specifics_.begin();
      //   for (; it != specifics_.end(); ++it) {
      //     it->second.clear();
      //   }
   }
   
}