RecoTracker/ConversionSeedGenerators/plugins/PrintRecoObjects.cc - master - githubmirror/cmssw

RecoTracker/ConversionSeedGenerators/plugins/PrintRecoObjects.cc
76 строк3 KB
Jonas Rembser
improve TrajectorySeed interface to recHits
27 май 2020, 19:03
27 май 2020, 19:03453bba8
#include "PrintRecoObjects.h"
#include "DataFormats/SiStripDetId/interface/SiStripDetId.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"

void PrintRecoObjects::print(std::stringstream& ss, const SiStripCluster& clus) {
  ss << "\n\t"
     << " firstStrip " << clus.firstStrip() << " size " << clus.amplitudes().size();
  float charge = 0;
  size_t idx = 0;
  for (auto adc : clus.amplitudes()) {
    charge += adc;
    ss << "\n\t\t strip " << ++idx << " adc " << (size_t)adc;
  }
  ss << "\n\t\t charge " << charge << " barycenter " << clus.barycenter() << std::endl;
}

void PrintRecoObjects::print(std::stringstream& ss, const TrajectorySeed& tjS) {
  ss << "\n\t nHits " << tjS.nHits() << "\n\t PTrajectoryStateOnDet: detid " << tjS.startingState().detId()
     << " tsos momentum" << tjS.startingState().parameters().momentum() << " pt "
     << sqrt(tjS.startingState().parameters().momentum().mag2() -
             tjS.startingState().parameters().momentum().y() * tjS.startingState().parameters().momentum().y())
     << " charge " << tjS.startingState().parameters().charge() << "\n\t error ";
  for (size_t ie = 0; ie < 15; ++ie)
    ss << "\t " << tjS.startingState().error(ie);
  for (auto const& recHit : tjS.recHits()) {
    ss << "\n\t TrackingRecHit on detid " << recHit.geographicalId().rawId() << " \t localPos "
       << recHit.localPosition();
  }
}

void PrintRecoObjects::print(std::stringstream& ss, const uint32_t& detid, const TrackerTopology* tTopo) const {
  ss << getString(detid, tTopo);
}

std::string PrintRecoObjects::getString(uint32_t detid, const TrackerTopology* tTopo) const {
  std::string append = " ";
  char cindex[128];
  SiStripDetId a(detid);
  if (a.subdetId() == 3) {
    append += "_TIB_L";
    sprintf(cindex, "%d", tTopo->tibLayer(detid));
  } else if (a.subdetId() == 4) {
    if (tTopo->tidSide(detid) == 1) {
      append += "_M_D";
    } else {
      append += "_M_P";
    }
    sprintf(cindex, "%d", tTopo->tidWheel(detid));
  } else if (a.subdetId() == 5) {
    append += "_TOB_L";
    sprintf(cindex, "%d", tTopo->tobLayer(detid));
  } else if (a.subdetId() == 6) {
    if (tTopo->tecSide(detid) == 1) {
      append += "_TEC_M";
    } else {
      append += "_TEC_P";
    }
    sprintf(cindex, "%d", tTopo->tecWheel(detid));
  }

  append += std::string(cindex);
  return append;
}

void PrintRecoObjects::print(std::stringstream& ss, const reco::Track* track, const math::XYZPoint& vx) {
  ss << "[PrintObject] "
     << "\n\tcharge \t" << track->charge() << "\talgo \t" << track->algo() << "\n\treferencePoint \t"
     << track->referencePoint() << " \t r: " << track->referencePoint().rho() << "\n\tinnerpos \t"
     << track->innerPosition() << " \t r: " << track->innerPosition().rho();

  ss << "\n\tmomentum \t" << track->momentum() << " pt " << track->momentum().rho() << "\n\tinnerMom \t"
     << track->innerMomentum() << " pt " << track->innerMomentum().rho();

  ss << "\n\tinnerok \t" << track->innerOk() << "\n\tinnerdetid \t" << track->innerDetId() << "\n\tdxy    \t"
     << track->dxy() << "\n\tdxy(vx) \t" << track->dxy(vx) << "\t where vx \t" << vx << std::endl;
}