MuonSeedOrcaPatternRecognition.h

Go to the documentation of this file.

#ifndef MuonSeedGenerator_MuonSeedOrcaPatternRecognition_h
#define MuonSeedGenerator_MuonSeedOrcaPatternRecognition_h

#include "RecoMuon/MuonSeedGenerator/src/MuonSeedVPatternRecognition.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Utilities/interface/ESGetToken.h"
#include "RecoMuon/MeasurementDet/interface/MuonDetLayerMeasurements.h"
#include "RecoMuon/DetLayers/interface/MuonDetLayerGeometry.h"
#include "RecoMuon/Records/interface/MuonRecoGeometryRecord.h"

class MuonSeedOrcaPatternRecognition : public MuonSeedVPatternRecognition {
public:
  explicit MuonSeedOrcaPatternRecognition(const edm::ParameterSet& pset, edm::ConsumesCollector& iC);

  void produce(const edm::Event& event,
               const edm::EventSetup& eSetup,
               std::vector<MuonRecHitContainer>& result) override;

private:
  // aalocates a zeroed array of a given size
  bool* zero(unsigned listSize);

  void endcapPatterns(const MuonRecHitContainer& me11,
                      const MuonRecHitContainer& me12,
                      const MuonRecHitContainer& me2,
                      const MuonRecHitContainer& me3,
                      const MuonRecHitContainer& me4,
                      const MuonRecHitContainer& me0,
                      const MuonRecHitContainer& mb1,
                      const MuonRecHitContainer& mb2,
                      const MuonRecHitContainer& mb3,
                      bool* MB1,
                      bool* MB2,
                      bool* MB3,
                      std::vector<MuonRecHitContainer>& result);

  void complete(MuonRecHitContainer& seedSegments, const MuonRecHitContainer& recHits, bool* used = nullptr) const;

  MuonRecHitPointer bestMatch(const ConstMuonRecHitPointer& first, MuonRecHitContainer& good_rhit) const;
  // some score to measure how well the two hits match
  double discriminator(const ConstMuonRecHitPointer& first, MuonRecHitPointer& other) const;
  // see if it's OK to add
  bool check(const MuonRecHitContainer& segments);
  bool isCrack(const ConstMuonRecHitPointer& segment) const;
  void rememberCrackSegments(const MuonRecHitContainer& segments, MuonRecHitContainer& crackSegments) const;

  void dumpLayer(const char* name, const MuonRecHitContainer& segments) const;

  MuonRecHitContainer filterSegments(const MuonRecHitContainer& segments, double dThetaCut) const;
  void filterOverlappingChambers(MuonRecHitContainer& segments) const;
  bool isME1A(const ConstMuonRecHitPointer& segment) const;
  int countHits(const MuonRecHitPointer& segment) const;
  // can mark other ME1A as used if one is.
  void markAsUsed(int nr, const MuonRecHitContainer& recHits, bool* used) const;
  std::vector<double> theCrackEtas;
  double theCrackWindow;
  MuonDetLayerMeasurements* muonMeasurements;
  double theDeltaPhiWindow;
  double theDeltaEtaWindow;
  double theDeltaCrackWindow;
  edm::ESGetToken<MuonDetLayerGeometry, MuonRecoGeometryRecord> muonLayersToken;
};

#endif