/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/DetectorDescription/RegressionTest/src/DDExpandedViewDump.cc

Go to the documentation of this file.

namespace std { } using namespace std;
#include <sstream>
#include <cstdlib>
#include <set>
#include <algorithm>
#include <string>
#include "DetectorDescription/RegressionTest/interface/DDExpandedViewDump.h"

void DDExpandedViewDump(ostream & os, DDExpandedView & ex, size_t skip, size_t sto)
{
  typedef set<string>::iterator s_iter;
  set<string> result;
  bool go(true);
  int  count(0);
  bool dotrans(true);
  if (getenv("DDNOTRANS")) dotrans=false;
  ++skip;
  while(go) {
    if ((count % 5000)==0) cout << count << ' ' << flush;
    if (sto > 0) if ((count % sto)==0) break;
    ++count;
    if((count % skip) ==0){
      stringstream s;
      s.setf(ios_base::fixed,ios_base::floatfield);
      s.precision(2);
      s << ex.logicalPart().name() << ' '
        << ex.copyno() << ' ' 
        << ex.geoHistory() << " r="
        << ex.geoHistory().back().posdata()->rot_.name() << "\n"; 
      DDRotationMatrix rm = ex.rotation();
      {
        double v[9]; rm.GetComponents(v,v+9);
        s << "R=(";
        s << v[0] << ' ' << v[1] << ' ' << v[2] << endl;
        s << v[3] << ' ' << v[4] << ' ' << v[5] << endl;
        s << v[6] << ' ' << v[7] << ' ' << v[7] << endl;
      }
      rm = ex.geoHistory().back().posdata()->rotation();
      {
        double v[9]; rm.GetComponents(v,v+9);
        s << "r=(";
        s << v[0] << ' ' << v[1] << ' ' << v[2] << endl;
        s << v[3] << ' ' << v[4] << ' ' << v[5] << endl;
        s << v[6] << ' ' << v[7] << ' ' << v[7] << endl;
      }
      if (dotrans) {
        s << "T=("
          << ex.translation().x() << ','
          << ex.translation().y() << ','
          << ex.translation().z() << ") "
          ;
      }      
      pair<s_iter,bool> ins = result.insert(s.str());
      if(!ins.second) {
        cout << "DDExpandedViewDump:ERROR: duplicated=" << s.str() << endl;
      } 
    }
    go = ex.next();
  }
  s_iter it(result.begin()), ed(result.end());
  for(;it != ed;++it) {
    os << *it << endl;
  }
}