MTDNavigationPrinter.cc

Go to the documentation of this file.

#include "RecoMTD/Navigation/interface/MTDNavigationPrinter.h"

#include "TrackingTools/DetLayers/interface/BarrelDetLayer.h" 
#include "TrackingTools/DetLayers/interface/ForwardDetLayer.h"
#include "DataFormats/GeometrySurface/interface/BoundCylinder.h"
#include "DataFormats/GeometrySurface/interface/BoundDisk.h"
#include "RecoMTD/DetLayers/interface/MTDDetLayerGeometry.h"
#include "RecoTracker/TkDetLayers/interface/GeometricSearchTracker.h"
#include "RecoMTD/Navigation/interface/MTDNavigationSchool.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <iostream>
#include <iomanip>
using namespace std;

// #define VI_DEBUG

#ifdef VI_DEBUG
#define PRINT(x) std::cout << x << ' '
#else
#define PRINT(x) edm::LogInfo(x)
#endif

MTDNavigationPrinter::MTDNavigationPrinter(const MTDDetLayerGeometry * mtdLayout,  MTDNavigationSchool const & sh,   bool enableBTL, bool enableETL) :
  school(&sh) {

  PRINT("MTDNavigationPrinter")<< "MTDNavigationPrinter::MTDNavigationPrinter" << std::endl;
  PRINT("MTDNavigationPrinter")<<"================================" << std::endl;
  PRINT("MTDNavigationPrinter")<< "BARREL:" << std::endl;
  vector<const DetLayer*> barrel;
  if ( enableBTL ) barrel = mtdLayout->allBarrelLayers();
  else barrel = mtdLayout->allBarrelLayers();

  PRINT("MTDNavigationPrinter")<<"There are "<<barrel.size()<<" Barrel DetLayers";
  for (auto i: barrel ) printLayer(i);
  PRINT("MTDNavigationPrinter")<<"================================" << std::endl;
  PRINT("MTDNavigationPrinter")  << "BACKWARD:" << std::endl;

  vector<const DetLayer*> backward;

  if ( enableETL ) backward = mtdLayout->allBackwardLayers();  
  else backward = mtdLayout->allBackwardLayers();

  PRINT("MTDNavigationPrinter")<<"There are "<<backward.size()<<" Backward DetLayers";
  for (auto i : backward ) printLayer(i);
  PRINT("MTDNavigationPrinter") << "==============================" << std::endl;
  PRINT("MTDNavigationPrinter") << "FORWARD:" << std::endl;
  vector<const DetLayer*> forward;

  if ( enableETL ) forward = mtdLayout->allForwardLayers();  
  else forward = mtdLayout->allForwardLayers();

  PRINT("MTDNavigationPrinter")<<"There are "<<forward.size()<<" Forward DetLayers" << std::endl;
  for (auto i : forward ) printLayer(i);

}

MTDNavigationPrinter::MTDNavigationPrinter(const MTDDetLayerGeometry * mtdLayout,  MTDNavigationSchool const & sh, const GeometricSearchTracker * tracker)  :
  school(&sh){

  PRINT("MTDNavigationPrinter")<< "MTDNavigationPrinter::MTDNavigationPrinter" << std::endl ;
//  vector<BarrelDetLayer*>::const_iterator tkiter;
//  vector<ForwardDetLayer*>::const_iterator tkfiter;
  PRINT("MTDNavigationPrinter")<<"================================" << std::endl;
  PRINT("MTDNavigationPrinter")<< "BARREL:" << std::endl;
  const vector<const BarrelDetLayer*>& tkbarrel = tracker->barrelLayers();
  PRINT("MTDNavigationPrinter")<<"There are "<<tkbarrel.size()<<" Tk Barrel DetLayers" << std::endl;
//  for ( tkiter = tkbarrel.begin(); tkiter != tkbarrel.end(); tkiter++ ) printLayer(*tkiter);
  vector<const DetLayer*> barrel = mtdLayout->allBarrelLayers();
  PRINT("MTDNavigationPrinter")<<"There are "<<barrel.size()<<" Mu Barrel DetLayers";
  for ( auto i : barrel ) printLayer(i);
  PRINT("MTDNavigationPrinter")<<"================================" << std::endl;
  PRINT("MTDNavigationPrinter")  << "BACKWARD:" << std::endl;
  const vector<const ForwardDetLayer*>& tkbackward = tracker->negForwardLayers();
  PRINT("MTDNavigationPrinter")<<"There are "<<tkbackward.size()<<" Tk Backward DetLayers" << std::endl;
  vector<const DetLayer*> backward = mtdLayout->allBackwardLayers();
  PRINT("MTDNavigationPrinter")<<"There are "<<backward.size()<<" Mu Backward DetLayers << std::endl";
  for (auto i : backward ) printLayer(i);
  PRINT("MTDNavigationPrinter") << "==============================" << std::endl;
  PRINT("MTDNavigationPrinter") << "FORWARD:" << std::endl;
  const vector<const ForwardDetLayer*>& tkforward =  tracker->posForwardLayers();
  PRINT("MTDNavigationPrinter")<<"There are "<<tkforward.size()<<" Tk Forward DetLayers" << std::endl;
//  for ( tkfiter = tkforward.begin(); tkfiter != tkforward.end(); tkfiter++ ) printLayer(*tkfiter);

  vector<const DetLayer*> forward = mtdLayout->allForwardLayers();
  PRINT("MTDNavigationPrinter")<<"There are "<<forward.size()<<" Mu Forward DetLayers";
  for ( auto i : forward ) printLayer(i);

}

void MTDNavigationPrinter::printLayer(const DetLayer* layer) const {
  vector<const DetLayer*> nextLayers = school->nextLayers(*layer,insideOut);
  vector<const DetLayer*> compatibleLayers = school->compatibleLayers(*layer,insideOut);
  if (const BarrelDetLayer* bdl = dynamic_cast<const BarrelDetLayer*>(layer)) {
    PRINT("MTDNavigationPrinter") 
         << layer->location() << " " << layer->subDetector() << " layer at R: "
         << setiosflags(ios::showpoint | ios::fixed)
         << setw(8) << setprecision(2)
         << bdl->specificSurface().radius() << "  length: "
         << setw(6) << setprecision(2)
         << layer->surface().bounds().length() << std::endl;
          
  }
  else if (const ForwardDetLayer* fdl = dynamic_cast<const ForwardDetLayer*>(layer)) {
    PRINT("MTDNavigationPrinter") << endl
         << layer->location() << " " << layer->subDetector() << "layer at z: "
         << setiosflags(ios::showpoint | ios::fixed)
         << setw(8) << setprecision(2)
         << layer->surface().position().z() << "  inner r: "
         << setw(6) << setprecision(2)
         << fdl->specificSurface().innerRadius() << "  outer r: "
         << setw(6) << setprecision(2)
         << fdl->specificSurface().outerRadius() << std::endl;
  }
  PRINT("MTDNavigationPrinter") << " has " << nextLayers.size() << " next layers in the direction inside-out: " << std::endl;
  printLayers(nextLayers);

  nextLayers.clear();
  nextLayers = school->nextLayers(*layer,outsideIn);

   PRINT("MTDNavigationPrinter") << " has " << nextLayers.size() << " next layers in the direction outside-in: " << std::endl;
  printLayers(nextLayers);

  PRINT("MTDNavigationPrinter") << " has " << compatibleLayers.size() << " compatible layers in the direction inside-out:: " << std::endl;
  printLayers(compatibleLayers);
  compatibleLayers.clear();
  compatibleLayers = school->compatibleLayers(*layer,outsideIn);
  
  PRINT("MTDNavigationPrinter") << " has " << compatibleLayers.size() << " compatible layers in the direction outside-in: " << std::endl;
  printLayers(compatibleLayers);

}

void MTDNavigationPrinter::printLayers(const vector<const DetLayer*>& nextLayers) const {

  for ( vector<const DetLayer*>::const_iterator inext = nextLayers.begin();
      inext != nextLayers.end(); inext++ ) {

     PRINT("MTDNavigationPrinter") << " --> " << std::endl; 
     if ( (*inext)->location() == GeomDetEnumerators::barrel ) {
      const BarrelDetLayer* l = dynamic_cast<const BarrelDetLayer*>(&(**inext));
      PRINT("MTDNavigationPrinter") << (*inext)->location() << " "
           << (*inext)->subDetector()
           << " layer at R: "
           << setiosflags(ios::showpoint | ios::fixed)
           << setw(8) << setprecision(2)
           << l->specificSurface().radius() << "   " << std::endl;
    }
    else {
      const ForwardDetLayer* l = dynamic_cast<const ForwardDetLayer*>(&(**inext));
       PRINT("MTDNavigationPrinter") << (*inext)->location() << " "
           << (*inext)->subDetector()
           << " layer at z: "
           << setiosflags(ios::showpoint | ios::fixed)
           << setw(8) << setprecision(2)
           << l->surface().position().z() << "   " << std::endl;
    }
    PRINT("MTDNavigationPrinter") << setiosflags(ios::showpoint | ios::fixed)
         << setprecision(1)
         << setw(6) << (*inext)->surface().bounds().length() << ", "
         << setw(6) << (*inext)->surface().bounds().width() << ", "
         << setw(4) <<(*inext)->surface().bounds().thickness() << " : " 
         << (*inext)->surface().position() << std::endl;
  }

}