DTHitQualityUtils.h

Go to the documentation of this file.

#ifndef Validation_DTRecHits_DTHitQualityUtils_h
#define Validation_DTRecHits_DTHitQualityUtils_h

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

#include <map>
#include <atomic>

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);

}
#endif