TrackingTools/TrackAssociator/interface/TrackDetMatchInfo.h - master - githubmirror/cmssw

TrackingTools/TrackAssociator/interface/TrackDetMatchInfo.h
119 строк4 KB
Matti Kortelainen
Migrate TrackDetectorAssociator to EventSetup consumes
15 окт 2020, 04:05
15 окт 2020, 04:05968b21d
#ifndef TrackAssociator_TrackDetMatchInfo_h
#define TrackAssociator_TrackDetMatchInfo_h

#include "DataFormats/CaloTowers/interface/CaloTower.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
#include "DataFormats/HcalRecHit/interface/HBHERecHit.h"
#include "DataFormats/HcalRecHit/interface/HORecHit.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "SimDataFormats/Track/interface/SimTrackContainer.h"
#include "TrackingTools/TrackAssociator/interface/TAMuonChamberMatch.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
class TrackDetMatchInfo {
public:
  enum EnergyType { EcalRecHits, HcalRecHits, HORecHits, TowerTotal, TowerEcal, TowerHcal, TowerHO };

  TrackDetMatchInfo();

  /// energy in detector elements crossed by the track by types
  double crossedEnergy(EnergyType);

  /// cone energy around the track direction at the origin (0,0,0)
  /// ( not well defined for tracks originating away from IP)
  double coneEnergy(double dR, EnergyType);

  /// get energy of the NxN shape (N = 2*gridSize + 1) around given detector element
  double nXnEnergy(const DetId&, EnergyType, int gridSize = 1);

  /// get energy of the NxN shape (N = 2*gridSize + 1) around track projection
  double nXnEnergy(EnergyType, int gridSize = 1);

  /// Find detector elements with highest energy deposition
  DetId findMaxDeposition(EnergyType);
  DetId findMaxDeposition(EnergyType, int gridSize);
  DetId findMaxDeposition(const DetId&, EnergyType, int gridSize);

  /// Track position at different parts of the calorimeter
  math::XYZPoint trkGlobPosAtEcal;
  math::XYZPoint trkGlobPosAtHcal;
  math::XYZPoint trkGlobPosAtHO;

  GlobalVector trkMomAtEcal;
  GlobalVector trkMomAtHcal;
  GlobalVector trkMomAtHO;

  bool isGoodEcal;
  bool isGoodHcal;
  bool isGoodCalo;
  bool isGoodHO;
  bool isGoodMuon;

  /// hits in the cone
  std::vector<const EcalRecHit*> ecalRecHits;
  std::vector<const HBHERecHit*> hcalRecHits;
  std::vector<const HORecHit*> hoRecHits;
  std::vector<const CaloTower*> towers;

  /// hits in detector elements crossed by a track
  std::vector<const EcalRecHit*> crossedEcalRecHits;
  std::vector<const HBHERecHit*> crossedHcalRecHits;
  std::vector<const HORecHit*> crossedHORecHits;
  std::vector<const CaloTower*> crossedTowers;

  /// detector elements crossed by a track
  /// (regardless of whether energy was deposited or not)
  std::vector<DetId> crossedEcalIds;
  std::vector<DetId> crossedHcalIds;
  std::vector<DetId> crossedHOIds;
  std::vector<DetId> crossedTowerIds;
  std::vector<DetId> crossedPreshowerIds;

  std::vector<TAMuonChamberMatch> chambers;

  /// track info
  FreeTrajectoryState stateAtIP;

  /// MC truth info
  const SimTrack* simTrack;
  double ecalTrueEnergy;
  double hcalTrueEnergy;
  double hcalTrueEnergyCorrected;

  /// Obsolete methods and data members for backward compatibility.
  /// Will be removed in future releases.
  reco::TrackRef trackRef_;
  SimTrackRef simTrackRef_;

  double ecalCrossedEnergy();
  double ecalConeEnergy();
  double hcalCrossedEnergy();
  double hcalConeEnergy();
  double hoCrossedEnergy();
  double hoConeEnergy();

  double ecalTowerEnergy() { return crossedEnergy(TowerEcal); }
  double ecalTowerConeEnergy() { return coneEnergy(999, TowerEcal); }
  double hcalTowerEnergy() { return crossedEnergy(TowerHcal); }
  double hcalTowerConeEnergy() { return coneEnergy(999, TowerHcal); }
  double hoTowerEnergy() { return crossedEnergy(TowerHO); }
  double hoTowerConeEnergy() { return coneEnergy(999, TowerHO); }

  double ecalEnergy() { return ecalCrossedEnergy(); }
  double hcalEnergy() { return hcalCrossedEnergy(); }
  double hoEnergy() { return hoCrossedEnergy(); }

  int numberOfSegments() const;
  int numberOfSegmentsInStation(int station) const;
  int numberOfSegmentsInStation(int station, int detector) const;
  int numberOfSegmentsInDetector(int detector) const;

  void setCaloGeometry(const CaloGeometry* geometry) { caloGeometry = geometry; }
  GlobalPoint getPosition(const DetId&);
  std::string dumpGeometry(const DetId&);

private:
  bool insideCone(const DetId&, const double);
  const CaloGeometry* caloGeometry;
};
#endif