MuonPatternRecoDumper.cc

Go to the documentation of this file.

 
// This Class Header
#include "RecoMuon/TrackingTools/interface/MuonPatternRecoDumper.h"
 
// Collaborating Class Header
#include "DataFormats/GeometrySurface/interface/BoundCylinder.h"
#include "DataFormats/GeometrySurface/interface/BoundDisk.h"
#include "TrackingTools/DetLayers/interface/DetLayer.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
 
#include "DataFormats/MuonDetId/interface/MuonSubdetId.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "DataFormats/MuonDetId/interface/RPCDetId.h"
#include "DataFormats/MuonDetId/interface/GEMDetId.h"
#include "DataFormats/MuonDetId/interface/ME0DetId.h"
 
#include <sstream>
 
using namespace std;
 
// Constructor
MuonPatternRecoDumper::MuonPatternRecoDumper() {}
 
// Destructor
MuonPatternRecoDumper::~MuonPatternRecoDumper() {}
 
// Operations
 
string MuonPatternRecoDumper::dumpLayer(const DetLayer* layer) const {
  stringstream output;
 
  const BoundSurface* sur = nullptr;
  const BoundCylinder* bc = nullptr;
  const BoundDisk* bd = nullptr;
 
  sur = &(layer->surface());
  if ((bc = dynamic_cast<const BoundCylinder*>(sur))) {
    output << "  Cylinder of radius: " << bc->radius() << endl;
  } else if ((bd = dynamic_cast<const BoundDisk*>(sur))) {
    output << "  Disk at: " << bd->position().z() << endl;
  }
  return output.str();
}
 
string MuonPatternRecoDumper::dumpFTS(const FreeTrajectoryState& fts) const {
  stringstream output;
 
  output << " pos: " << fts.position() << " radius: " << fts.position().perp() << endl
         << " charge*pt: " << fts.momentum().perp() * fts.parameters().charge() << " eta: " << fts.momentum().eta()
         << " phi: " << fts.momentum().phi() << endl;
 
  return output.str();
}
 
string MuonPatternRecoDumper::dumpTSOS(const TrajectoryStateOnSurface& tsos) const {
  stringstream output;
 
  output << tsos << endl;
  output << "dir: " << tsos.globalDirection() << endl;
  output << dumpFTS(*tsos.freeTrajectoryState());
 
  return output.str();
}
 
string MuonPatternRecoDumper::dumpMuonId(const DetId& id) const {
  stringstream output;
 
  if (id.subdetId() == MuonSubdetId::DT) {
    DTWireId wireId(id.rawId());
 
    output << "(DT): " << wireId << endl;
  } else if (id.subdetId() == MuonSubdetId::CSC) {
    CSCDetId chamberId(id.rawId());
    output << "(CSC): " << chamberId << endl;
  } else if (id.subdetId() == MuonSubdetId::GEM) {
    GEMDetId chamberId(id.rawId());
    output << "(GEM): " << chamberId << endl;
  } else if (id.subdetId() == MuonSubdetId::ME0) {
    ME0DetId chamberId(id.rawId());
    output << "(ME0): " << chamberId << endl;
  } else if (id.subdetId() == MuonSubdetId::RPC) {
    RPCDetId chamberId(id.rawId());
    output << "(RPC): " << chamberId << endl;
  } else
    output << "The DetLayer is not a valid Muon DetLayer. ";
 
  return output.str();
}