DTHitQualityUtils.h

Go to the documentation of this file.

#ifndef Validation_DTRecHits_DTHitQualityUtils_h
#define Validation_DTRecHits_DTHitQualityUtils_h

#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/GeometryVector/interface/LocalVector.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"

#include <atomic>
#include <map>

class PSimHit;
class DTGeometry;

namespace DTHitQualityUtils {

  std::map<DTWireId, edm::PSimHitContainer> mapSimHitsPerWire(const edm::PSimHitContainer &simhits);
  std::map<DTWireId, const PSimHit *> mapMuSimHitsPerWire(
      const std::map<DTWireId, edm::PSimHitContainer> &simHitWireMap);
  const PSimHit *findMuSimHit(const edm::PSimHitContainer &hits);
  std::pair<const PSimHit *, const PSimHit *> findMuSimSegment(
      const std::map<DTWireId, const PSimHit *> &mapWireAndMuSimHit);
  std::pair<LocalVector, LocalPoint> findMuSimSegmentDirAndPos(
      const std::pair<const PSimHit *, const PSimHit *> &inAndOutSimHit, const DetId detId, const DTGeometry *muonGeom);
  std::pair<double, double> findSegmentAlphaAndBeta(const LocalVector &direction);

  // Set the verbosity level
  extern std::atomic<bool> debug;

  // Find angle error
  double sigmaAngle(double Angle, double sigma2TanAngle);

}  // namespace DTHitQualityUtils
#endif