DTHitQualityUtils.h

Go to the documentation of this file.

#ifndef DTHitQualityUtils_H
#define 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>

class PSimHit;
class DTGeometry;

class DTHitQualityUtils {
public:

  DTHitQualityUtils();


  virtual ~DTHitQualityUtils();

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

  // Set the verbosity level
  static bool debug; 

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

protected:

private:

};
#endif