/data/refman/pasoursint/CMSSW_4_1_8_patch13/src/DetectorDescription/Core/src/DDComparator.cc

Go to the documentation of this file.

#include "DetectorDescription/Core/interface/DDComparator.h"
//#include "DetectorDescription/Core/interface/DDPartSelection.h"
//#include "DetectorDescription/Base/interface/DDException.h"

#include<map>
#include<iostream>

namespace {

  struct Counter {
    int old;
    int diff;
    int t;
    int f;
    std::map<int,int> res;
    std::map<int,int> siz;
    Counter() :old(0),diff(0),t(0),f(0){}
    ~Counter() {
      for (std::map<int,int>::const_iterator p=res.begin();p!=res.end(); ++p)
        std::cout << (*p).first <<","<<(*p).second <<" "; 
      if (!res.empty()) std::cout << std::endl;
       for (std::map<int,int>::const_iterator p=siz.begin();p!=siz.end(); ++p)
        std::cout << (*p).first <<","<<(*p).second <<" "; 
       if (!siz.empty()) std::cout << std::endl;
    }
    void add(bool r, int id, int im) {
      if(r) {
        ++t;
        return;
      }
      else ++f;
      /*
      std::map<int,int>::iterator p;
      p= res.find(im-id);
      if (p==res.end()) res[im-id]=1;
      else ++(*p).second;
      p= siz.find(id);
      if (p==siz.end()) siz[id]=1;
      else ++(*p).second;
      */
    }

  };

  inline Counter & counter() {
    static Counter local;
    return local;
  }

}

// reason for the ctor: reference initialization at construction.
// FIXME: DDCompareEqual: use pointers instead of references, initialize to 0 and 
// FIXME: do the check in operator() instead of in the ctor

bool DDCompareEqual::operator() (const DDGeoHistory &, const DDPartSelection &) 
{
  return (*this)();
}

bool DDCompareEqual::operator() () 
{

  // don't compare, if history or partsel is empty! (see ctor) 
  bool result(absResult_);
  
  /*
     sIndex_  =  running index in the part-selection-std::vector
     sMax_    =  max. value + 1 of sIndex_
     hIndex_  =  runninig index in the geo-history-std::vector
     hMax_    =  max. value + 1 of hIndex_
     sLp_     =  current LogicalPart (the redir-ptr!) in the part-selection-std::vector
     hLp_     =  current LogicalPart (the redir-ptr!) in the geo-history-std::vector
     sCopyno_ =  current copy-no in the part-selection-std::vector
  */
  //DCOUT('U', "DDCompareEqual: comparing");
    
  while(result && sIndex_ < sMax_) {
    sLp_ = partsel_[sIndex_].lp_;
    sCopyno_ = partsel_[sIndex_].copyno_;
    ddselection_type stype = partsel_[sIndex_].selectionType_; 
    switch (stype) {
   
     case ddanylogp:
        result=nextAnylogp(); 
        break;
   
     case ddanyposp:
        result=nextAnyposp();
        break;
         
     case ddchildlogp:
        result=nextChildlogp();
        break;
        
     case ddchildposp:
        result=nextChildposp();
        break;
        
     case ddanychild:
        ++sIndex_;
        ++hIndex_;
        result=true;
        break;
     
     // ddanynode IS NOT SUPPORTED IN PROTOTYPE SW !!!!
     case ddanynode:
        result=false;
        break;
                
     default:
      result=false;
      //throw DDException("DDCompareEqual: undefined state!");
    }
    ++sIndex_;
  }
  counter().add(result,sIndex_, sMax_);
  return result;
}


bool DDCompareEqual::nextAnylogp()
{
  /* does not help, most of the time is spent when FALSE....
  static size_t hIndexOld=0;
  // hope same position that previous
  {
    size_t oldH = hIndexOld; // thread safe??? 
    if (oldH<hMax_ && sLp_== hist_[oldH].logicalPart()) {
      hIndex_ = oldH+1;
      ++counter.old;
      return true;
    }
  }
  */
  register size_t hi = hIndex_;
  while (hi < hMax_) {
    if (sLp_==hist_[hi].logicalPart()) {
      // hIndexOld=hIndex_;
      // ++counter.diff;
      hIndex_ = hi+1;
      return true;
    }
    ++hi;  
  }
  hIndex_ = hi;
  return false;
}


bool DDCompareEqual::nextAnyposp()
{
  bool result(false);
  while (hIndex_ < hMax_) {
    if (sLp_ == hist_[hIndex_].logicalPart() && 
        sCopyno_ == hist_[hIndex_].copyno() ) 
     { result=true;
       ++hIndex_; 
       break; 
     }
    ++hIndex_;
  }    
  return result;
}


bool DDCompareEqual::nextChildlogp()
{
  bool result(false);
  if (hIndex_ < hMax_) {
    if (sLp_ == hist_[hIndex_].logicalPart()) {
      ++hIndex_;
      result=true;
    }
  }
  return result;
}


bool DDCompareEqual::nextChildposp()
{
  bool result(false);
  if (hIndex_ < hMax_) {
    if (sLp_ == hist_[hIndex_].logicalPart() &&
        sCopyno_ == hist_[hIndex_].copyno() ) {
      ++hIndex_;
      result=true;
    }
  }
  return result;
}