TrackDetMatchInfo Class Reference

#include <TrackDetMatchInfo.h>

Public Types
enum	EnergyType {   EcalRecHits, HcalRecHits, HORecHits, TowerTotal,   TowerEcal, TowerHcal, TowerHO }

Public Member Functions
double	coneEnergy (double dR, EnergyType)
double	crossedEnergy (EnergyType)
	energy in detector elements crossed by the track by types More...
std::string	dumpGeometry (const DetId &)
double	ecalConeEnergy ()
double	ecalCrossedEnergy ()
double	ecalEnergy ()
double	ecalTowerConeEnergy ()
double	ecalTowerEnergy ()
DetId	findMaxDeposition (EnergyType)
	Find detector elements with highest energy deposition. More...
DetId	findMaxDeposition (EnergyType, int gridSize)
DetId	findMaxDeposition (const DetId &, EnergyType, int gridSize)
GlobalPoint	getPosition (const DetId &)
double	hcalConeEnergy ()
double	hcalCrossedEnergy ()
double	hcalEnergy ()
double	hcalTowerConeEnergy ()
double	hcalTowerEnergy ()
double	hoConeEnergy ()
double	hoCrossedEnergy ()
double	hoEnergy ()
double	hoTowerConeEnergy ()
double	hoTowerEnergy ()
int	numberOfSegments () const
int	numberOfSegmentsInDetector (int detector) const
int	numberOfSegmentsInStation (int station) const
int	numberOfSegmentsInStation (int station, int detector) const
double	nXnEnergy (const DetId &, EnergyType, int gridSize=1)
	get energy of the NxN shape (N = 2*gridSize + 1) around given detector element More...
double	nXnEnergy (EnergyType, int gridSize=1)
	get energy of the NxN shape (N = 2*gridSize + 1) around track projection More...
void	setCaloGeometry (const CaloGeometry *geometry)
	TrackDetMatchInfo ()

Public Attributes
std::vector< TAMuonChamberMatch >	chambers
std::vector< DetId >	crossedEcalIds
std::vector< const EcalRecHit * >	crossedEcalRecHits
	hits in detector elements crossed by a track More...
std::vector< DetId >	crossedHcalIds
std::vector< const HBHERecHit * >	crossedHcalRecHits
std::vector< DetId >	crossedHOIds
std::vector< const HORecHit * >	crossedHORecHits
std::vector< DetId >	crossedPreshowerIds
std::vector< DetId >	crossedTowerIds
std::vector< const CaloTower * >	crossedTowers
std::vector< const EcalRecHit * >	ecalRecHits
	hits in the cone More...
double	ecalTrueEnergy
std::vector< const HBHERecHit * >	hcalRecHits
double	hcalTrueEnergy
double	hcalTrueEnergyCorrected
std::vector< const HORecHit * >	hoRecHits
bool	isGoodCalo
bool	isGoodEcal
bool	isGoodHcal
bool	isGoodHO
bool	isGoodMuon
const SimTrack *	simTrack
	MC truth info. More...
SimTrackRef	simTrackRef_
FreeTrajectoryState	stateAtIP
	track info More...
std::vector< const CaloTower * >	towers
reco::TrackRef	trackRef_
math::XYZPoint	trkGlobPosAtEcal
	Track position at different parts of the calorimeter. More...
math::XYZPoint	trkGlobPosAtHcal
math::XYZPoint	trkGlobPosAtHO
GlobalVector	trkMomAtEcal
GlobalVector	trkMomAtHcal
GlobalVector	trkMomAtHO

Private Member Functions
bool	insideCone (const DetId &, const double)

Private Attributes
const CaloGeometry *	caloGeometry

Detailed Description

Definition at line 13 of file TrackDetMatchInfo.h.

Member Enumeration Documentation

EnergyType

enum TrackDetMatchInfo::EnergyType

Enumerator
EcalRecHits
HcalRecHits
HORecHits
TowerTotal
TowerEcal
TowerHcal
TowerHO

Definition at line 15 of file TrackDetMatchInfo.h.

{ EcalRecHits, HcalRecHits, HORecHits, TowerTotal, TowerEcal, TowerHcal, TowerHO };

Constructor & Destructor Documentation

TrackDetMatchInfo()

TrackDetMatchInfo::TrackDetMatchInfo

(

)

Definition at line 289 of file TrackDetMatchInfo.cc.

    : trkGlobPosAtEcal(0, 0, 0),
      trkGlobPosAtHcal(0, 0, 0),
      trkGlobPosAtHO(0, 0, 0),
      trkMomAtEcal(0, 0, 0),
      trkMomAtHcal(0, 0, 0),
      trkMomAtHO(0, 0, 0),
      isGoodEcal(false),
      isGoodHcal(false),
      isGoodCalo(false),
      isGoodHO(false),
      isGoodMuon(false),
      simTrack(nullptr),
      ecalTrueEnergy(-999),
      hcalTrueEnergy(-999) {}

Member Function Documentation

coneEnergy()

double TrackDetMatchInfo::coneEnergy	(	double	dR,
		EnergyType	type
	)

cone energy around the track direction at the origin (0,0,0) ( not well defined for tracks originating away from IP)

Definition at line 91 of file TrackDetMatchInfo.cc.

References crossedTowers, HGC3DClusterGenMatchSelector_cfi::dR, EcalRecHits, ecalRecHits, HCALHighEnergyHPDFilter_cfi::energy, Exception, HcalRecHits, hcalRecHits, HORecHits, hoRecHits, insideCone(), TowerEcal, TowerHcal, TowerHO, and TowerTotal.

Referenced by ecalConeEnergy(), ecalTowerConeEnergy(), hcalConeEnergy(), hcalTowerConeEnergy(), hoConeEnergy(), and hoTowerConeEnergy().

                                                               {
  double energy(0);
  switch (type) {
    case EcalRecHits: {
      for (std::vector<const EcalRecHit*>::const_iterator hit = ecalRecHits.begin(); hit != ecalRecHits.end(); hit++)
        if (insideCone((*hit)->detid(), dR))
          energy += (*hit)->energy();
    } break;
    case HcalRecHits: {
      for (std::vector<const HBHERecHit*>::const_iterator hit = hcalRecHits.begin(); hit != hcalRecHits.end(); hit++)
        if (insideCone((*hit)->detid(), dR))
          energy += (*hit)->energy();
    } break;
    case HORecHits: {
      for (std::vector<const HORecHit*>::const_iterator hit = hoRecHits.begin(); hit != hoRecHits.end(); hit++)
        if (insideCone((*hit)->detid(), dR))
          energy += (*hit)->energy();
    } break;
    case TowerTotal: {
      for (std::vector<const CaloTower*>::const_iterator hit = crossedTowers.begin(); hit != crossedTowers.end(); hit++)
        if (insideCone((*hit)->id(), dR))
          energy += (*hit)->energy();
    } break;
    case TowerEcal: {
      for (std::vector<const CaloTower*>::const_iterator hit = crossedTowers.begin(); hit != crossedTowers.end(); hit++)
        if (insideCone((*hit)->id(), dR))
          energy += (*hit)->emEnergy();
    } break;
    case TowerHcal: {
      for (std::vector<const CaloTower*>::const_iterator hit = crossedTowers.begin(); hit != crossedTowers.end(); hit++)
        if (insideCone((*hit)->id(), dR))
          energy += (*hit)->hadEnergy();
    } break;
    case TowerHO: {
      for (std::vector<const CaloTower*>::const_iterator hit = crossedTowers.begin(); hit != crossedTowers.end(); hit++)
        if (insideCone((*hit)->id(), dR))
          energy += (*hit)->outerEnergy();
    } break;
    default:
      throw cms::Exception("FatalError") << "Unknown calo energy type: " << type;
  }
  return energy;
}

crossedEnergy()

double TrackDetMatchInfo::crossedEnergy

(

EnergyType

type

)

energy in detector elements crossed by the track by types

Definition at line 39 of file TrackDetMatchInfo.cc.

References crossedEcalRecHits, crossedHcalRecHits, crossedHORecHits, crossedTowers, EcalRecHits, HCALHighEnergyHPDFilter_cfi::energy, Exception, HcalRecHits, HORecHits, TowerEcal, TowerHcal, TowerHO, and TowerTotal.

Referenced by ecalCrossedEnergy(), ecalTowerEnergy(), hcalCrossedEnergy(), hcalTowerEnergy(), hoCrossedEnergy(), and hoTowerEnergy().

                                                       {
  double energy(0);
  switch (type) {
    case EcalRecHits: {
      for (std::vector<const EcalRecHit*>::const_iterator hit = crossedEcalRecHits.begin();
           hit != crossedEcalRecHits.end();
           hit++)
        energy += (*hit)->energy();
    } break;
    case HcalRecHits: {
      for (std::vector<const HBHERecHit*>::const_iterator hit = crossedHcalRecHits.begin();
           hit != crossedHcalRecHits.end();
           hit++)
        energy += (*hit)->energy();
    } break;
    case HORecHits: {
      for (std::vector<const HORecHit*>::const_iterator hit = crossedHORecHits.begin(); hit != crossedHORecHits.end();
           hit++)
        energy += (*hit)->energy();
    } break;
    case TowerTotal: {
      for (std::vector<const CaloTower*>::const_iterator hit = crossedTowers.begin(); hit != crossedTowers.end(); hit++)
        energy += (*hit)->energy();
    } break;
    case TowerEcal: {
      for (std::vector<const CaloTower*>::const_iterator hit = crossedTowers.begin(); hit != crossedTowers.end(); hit++)
        energy += (*hit)->emEnergy();
    } break;
    case TowerHcal: {
      for (std::vector<const CaloTower*>::const_iterator hit = crossedTowers.begin(); hit != crossedTowers.end(); hit++)
        energy += (*hit)->hadEnergy();
    } break;
    case TowerHO: {
      for (std::vector<const CaloTower*>::const_iterator hit = crossedTowers.begin(); hit != crossedTowers.end(); hit++)
        energy += (*hit)->outerEnergy();
    } break;
    default:
      throw cms::Exception("FatalError") << "Unknown calo energy type: " << type;
  }
  return energy;
}

dumpGeometry()

std::string TrackDetMatchInfo::dumpGeometry

(

const DetId &

id

)

Definition at line 16 of file TrackDetMatchInfo.cc.

References caloGeometry, Exception, CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), point, and hiPixelPairStep_cff::points.

                                                         {
  if (!caloGeometry || !caloGeometry->getSubdetectorGeometry(id) ||
      !caloGeometry->getSubdetectorGeometry(id)->getGeometry(id)) {
    throw cms::Exception("FatalError") << "Failed to access geometry for DetId: " << id.rawId();
  }
  std::ostringstream oss;

  const CaloCellGeometry::CornersVec& points = caloGeometry->getSubdetectorGeometry(id)->getGeometry(id)->getCorners();
  for (CaloCellGeometry::CornersVec::const_iterator point = points.begin(); point != points.end(); ++point)
    oss << "(" << point->z() << ", " << point->perp() << ", " << point->eta() << ", " << point->phi() << "), \t";
  return oss.str();
}

ecalConeEnergy()

double TrackDetMatchInfo::ecalConeEnergy

(

)

Definition at line 547 of file TrackDetMatchInfo.cc.

References coneEnergy(), and EcalRecHits.

{ return coneEnergy(999, EcalRecHits); }

ecalCrossedEnergy()

double TrackDetMatchInfo::ecalCrossedEnergy

(

)

Definition at line 553 of file TrackDetMatchInfo.cc.

References crossedEnergy(), and EcalRecHits.

Referenced by ecalEnergy().

{ return crossedEnergy(EcalRecHits); }

ecalEnergy()

double TrackDetMatchInfo::ecalEnergy ( )

inline

Definition at line 102 of file TrackDetMatchInfo.h.

References ecalCrossedEnergy().

{ return ecalCrossedEnergy(); }

ecalTowerConeEnergy()

double TrackDetMatchInfo::ecalTowerConeEnergy ( )

inline

Definition at line 96 of file TrackDetMatchInfo.h.

References coneEnergy(), and TowerEcal.

{ return coneEnergy(999, TowerEcal); }

ecalTowerEnergy()

double TrackDetMatchInfo::ecalTowerEnergy ( )

inline

Definition at line 95 of file TrackDetMatchInfo.h.

References crossedEnergy(), and TowerEcal.

{ return crossedEnergy(TowerEcal); }

findMaxDeposition() [1/3]

DetId TrackDetMatchInfo::findMaxDeposition

(

EnergyType

type

)

Find detector elements with highest energy deposition.

Definition at line 305 of file TrackDetMatchInfo.cc.

References EcalRecHits, ecalRecHits, HCALHighEnergyHPDFilter_cfi::energy, Exception, HcalRecHits, hcalRecHits, HORecHits, hoRecHits, l1ctLayer2EG_cff::id, particleFlowClusterECALTimeSelected_cfi::maxEnergy, TowerEcal, TowerHcal, TowerHO, towers, and TowerTotal.

Referenced by findMaxDeposition().

                                                          {
  DetId id;
  float maxEnergy = -9999;
  switch (type) {
    case EcalRecHits: {
      for (std::vector<const EcalRecHit*>::const_iterator hit = ecalRecHits.begin(); hit != ecalRecHits.end(); hit++)
        if ((*hit)->energy() > maxEnergy) {
          maxEnergy = (*hit)->energy();
          id = (*hit)->detid();
        }
    } break;
    case HcalRecHits: {
      for (std::vector<const HBHERecHit*>::const_iterator hit = hcalRecHits.begin(); hit != hcalRecHits.end(); hit++)
        if ((*hit)->energy() > maxEnergy) {
          maxEnergy = (*hit)->energy();
          id = (*hit)->detid();
        }
    } break;
    case HORecHits: {
      for (std::vector<const HORecHit*>::const_iterator hit = hoRecHits.begin(); hit != hoRecHits.end(); hit++)
        if ((*hit)->energy() > maxEnergy) {
          maxEnergy = (*hit)->energy();
          id = (*hit)->detid();
        }
    } break;
    case TowerTotal:
    case TowerEcal:
    case TowerHcal:
    case TowerHO: {
      for (std::vector<const CaloTower*>::const_iterator hit = towers.begin(); hit != towers.end(); hit++) {
        double energy = 0;
        switch (type) {
          case TowerTotal:
            energy = (*hit)->energy();
            break;
          case TowerEcal:
            energy = (*hit)->emEnergy();
            break;
          case TowerHcal:
            energy = (*hit)->hadEnergy();
            break;
          case TowerHO:
            energy = (*hit)->energy();
            break;
          default:
            throw cms::Exception("FatalError") << "Unknown calo tower energy type: " << type;
        }
        if (energy > maxEnergy) {
          maxEnergy = energy;
          id = (*hit)->id();
        }
      }
    } break;
    default:
      throw cms::Exception("FatalError")
          << "Maximal energy deposition: unkown or not implemented energy type requested, type:" << type;
  }
  return id;
}

findMaxDeposition() [2/3]

DetId TrackDetMatchInfo::findMaxDeposition	(	EnergyType	type,
		int	gridSize
	)

Definition at line 511 of file TrackDetMatchInfo.cc.

References crossedEcalIds, crossedHcalIds, crossedHOIds, crossedTowerIds, EcalRecHits, Exception, findMaxDeposition(), HcalRecHits, HORecHits, TowerEcal, TowerHcal, TowerHO, and TowerTotal.

                                                                        {
  DetId id_max;
  switch (type) {
    case TowerTotal:
    case TowerHcal:
    case TowerEcal:
    case TowerHO:
      if (crossedTowerIds.empty())
        return id_max;
      return findMaxDeposition(crossedTowerIds.front(), type, gridSize);
      break;
    case EcalRecHits:
      if (crossedEcalIds.empty())
        return id_max;
      return findMaxDeposition(crossedEcalIds.front(), type, gridSize);
      break;
    case HcalRecHits:
      if (crossedHcalIds.empty())
        return id_max;
      return findMaxDeposition(crossedHcalIds.front(), type, gridSize);
      break;
    case HORecHits:
      if (crossedHOIds.empty())
        return id_max;
      return findMaxDeposition(crossedHOIds.front(), type, gridSize);
      break;
    default:
      throw cms::Exception("FatalError") << "Unkown or not implemented energy type requested, type:" << type;
  }
  return id_max;
}

findMaxDeposition() [3/3]

DetId TrackDetMatchInfo::findMaxDeposition	(	const DetId &	id,
		EnergyType	type,
		int	gridSize
	)

Definition at line 365 of file TrackDetMatchInfo.cc.

References funct::abs(), DetId::Calo, HLT_2023v12_cff::dEta, HLT_2023v12_cff::dPhi, DetId::Ecal, EcalBarrel, EcalEndcap, EcalRecHits, ecalRecHits, Exception, DetId::Hcal, HcalBarrel, HcalEndcap, HcalOuter, HcalRecHits, hcalRecHits, HORecHits, hoRecHits, CaloTowerDetId::ieta(), EBDetId::ieta(), HcalDetId::ieta(), DetIdInfo::info(), CaloTowerDetId::iphi(), EBDetId::iphi(), HcalDetId::iphi(), EEDetId::ix(), EEDetId::iy(), TowerEcal, TowerHcal, TowerHO, towers, TowerTotal, and EEDetId::zside().

                                                                                         {
  double energy_max(0);
  DetId id_max;
  if (id.rawId() == 0)
    return id_max;
  switch (type) {
    case TowerTotal:
    case TowerHcal:
    case TowerEcal:
    case TowerHO: {
      if (id.det() != DetId::Calo) {
        throw cms::Exception("FatalError") << "Wrong DetId. Expected CaloTower, but found:\n"
                                           << DetIdInfo::info(id, nullptr) << "\n";
      }
      CaloTowerDetId centerId(id);
      for (std::vector<const CaloTower*>::const_iterator hit = towers.begin(); hit != towers.end(); hit++) {
        CaloTowerDetId neighborId((*hit)->id());
        int dEta = abs((centerId.ieta() < 0 ? centerId.ieta() + 1 : centerId.ieta()) -
                       (neighborId.ieta() < 0 ? neighborId.ieta() + 1 : neighborId.ieta()));
        int dPhi = abs(centerId.iphi() - neighborId.iphi());
        if (abs(72 - dPhi) < dPhi)
          dPhi = 72 - dPhi;
        if (dEta <= gridSize && dPhi <= gridSize) {
          switch (type) {
            case TowerTotal:
              if (energy_max < (*hit)->energy()) {
                energy_max = (*hit)->energy();
                id_max = (*hit)->id();
              }
              break;
            case TowerEcal:
              if (energy_max < (*hit)->emEnergy()) {
                energy_max = (*hit)->emEnergy();
                id_max = (*hit)->id();
              }
              break;
            case TowerHcal:
              if (energy_max < (*hit)->hadEnergy()) {
                energy_max = (*hit)->hadEnergy();
                id_max = (*hit)->id();
              }
              break;
            case TowerHO:
              if (energy_max < (*hit)->outerEnergy()) {
                energy_max = (*hit)->outerEnergy();
                id_max = (*hit)->id();
              }
              break;
            default:
              throw cms::Exception("FatalError") << "Unknown calo tower energy type: " << type;
          }
        }
      }
    } break;
    case EcalRecHits: {
      if (id.det() != DetId::Ecal || (id.subdetId() != EcalBarrel && id.subdetId() != EcalEndcap)) {
        throw cms::Exception("FatalError") << "Wrong DetId. Expected EcalBarrel or EcalEndcap, but found:\n"
                                           << DetIdInfo::info(id, nullptr) << "\n";
      }
      // Since the ECAL granularity is small and the gap between EE and EB is significant,
      // energy is computed only within the system that contains the central element
      if (id.subdetId() == EcalBarrel) {
        EBDetId centerId(id);
        for (std::vector<const EcalRecHit*>::const_iterator hit = ecalRecHits.begin(); hit != ecalRecHits.end();
             hit++) {
          if ((*hit)->id().subdetId() != EcalBarrel)
            continue;
          EBDetId neighborId((*hit)->id());
          int dEta = abs((centerId.ieta() < 0 ? centerId.ieta() + 1 : centerId.ieta()) -
                         (neighborId.ieta() < 0 ? neighborId.ieta() + 1 : neighborId.ieta()));
          int dPhi = abs(centerId.iphi() - neighborId.iphi());
          if (abs(360 - dPhi) < dPhi)
            dPhi = 360 - dPhi;
          if (dEta <= gridSize && dPhi <= gridSize) {
            if (energy_max < (*hit)->energy()) {
              energy_max = (*hit)->energy();
              id_max = (*hit)->id();
            }
          }
        }
      } else {
        // Endcap
        EEDetId centerId(id);
        for (std::vector<const EcalRecHit*>::const_iterator hit = ecalRecHits.begin(); hit != ecalRecHits.end();
             hit++) {
          if ((*hit)->id().subdetId() != EcalEndcap)
            continue;
          EEDetId neighborId((*hit)->id());
          if (centerId.zside() == neighborId.zside() && abs(centerId.ix() - neighborId.ix()) <= gridSize &&
              abs(centerId.iy() - neighborId.iy()) <= gridSize) {
            if (energy_max < (*hit)->energy()) {
              energy_max = (*hit)->energy();
              id_max = (*hit)->id();
            }
          }
        }
      }
    } break;
    case HcalRecHits: {
      if (id.det() != DetId::Hcal || (id.subdetId() != HcalBarrel && id.subdetId() != HcalEndcap)) {
        throw cms::Exception("FatalError") << "Wrong DetId. Expected HE or HB, but found:\n"
                                           << DetIdInfo::info(id, nullptr) << "\n";
      }
      HcalDetId centerId(id);
      for (std::vector<const HBHERecHit*>::const_iterator hit = hcalRecHits.begin(); hit != hcalRecHits.end(); hit++) {
        HcalDetId neighborId((*hit)->id());
        int dEta = abs((centerId.ieta() < 0 ? centerId.ieta() + 1 : centerId.ieta()) -
                       (neighborId.ieta() < 0 ? neighborId.ieta() + 1 : neighborId.ieta()));
        int dPhi = abs(centerId.iphi() - neighborId.iphi());
        if (abs(72 - dPhi) < dPhi)
          dPhi = 72 - dPhi;
        if (dEta <= gridSize && dPhi <= gridSize) {
          if (energy_max < (*hit)->energy()) {
            energy_max = (*hit)->energy();
            id_max = (*hit)->id();
          }
        }
      }
    } break;
    case HORecHits: {
      if (id.det() != DetId::Hcal || (id.subdetId() != HcalOuter)) {
        throw cms::Exception("FatalError") << "Wrong DetId. Expected HO, but found:\n"
                                           << DetIdInfo::info(id, nullptr) << "\n";
      }
      HcalDetId centerId(id);
      for (std::vector<const HORecHit*>::const_iterator hit = hoRecHits.begin(); hit != hoRecHits.end(); hit++) {
        HcalDetId neighborId((*hit)->id());
        int dEta = abs((centerId.ieta() < 0 ? centerId.ieta() + 1 : centerId.ieta()) -
                       (neighborId.ieta() < 0 ? neighborId.ieta() + 1 : neighborId.ieta()));
        int dPhi = abs(centerId.iphi() - neighborId.iphi());
        if (abs(72 - dPhi) < dPhi)
          dPhi = 72 - dPhi;
        if (dEta <= gridSize && dPhi <= gridSize) {
          if (energy_max < (*hit)->energy()) {
            energy_max = (*hit)->energy();
            id_max = (*hit)->id();
          }
        }
      }
    } break;
    default:
      throw cms::Exception("FatalError") << "Unkown or not implemented energy type requested, type:" << type;
  }
  return id_max;
}

getPosition()

GlobalPoint TrackDetMatchInfo::getPosition

(

const DetId &

id

)

Definition at line 29 of file TrackDetMatchInfo.cc.

References caloGeometry, Exception, CaloSubdetectorGeometry::getGeometry(), and CaloGeometry::getSubdetectorGeometry().

Referenced by insideCone().

                                                          {
  // this part might be slow
  if (!caloGeometry || !caloGeometry->getSubdetectorGeometry(id) ||
      !caloGeometry->getSubdetectorGeometry(id)->getGeometry(id)) {
    throw cms::Exception("FatalError") << "Failed to access geometry for DetId: " << id.rawId();
    return GlobalPoint(0, 0, 0);
  }
  return caloGeometry->getSubdetectorGeometry(id)->getGeometry(id)->getPosition();
}

hcalConeEnergy()

double TrackDetMatchInfo::hcalConeEnergy

(

)

Definition at line 549 of file TrackDetMatchInfo.cc.

References coneEnergy(), and HcalRecHits.

{ return coneEnergy(999, HcalRecHits); }

hcalCrossedEnergy()

double TrackDetMatchInfo::hcalCrossedEnergy

(

)

Definition at line 555 of file TrackDetMatchInfo.cc.

References crossedEnergy(), and HcalRecHits.

Referenced by hcalEnergy().

{ return crossedEnergy(HcalRecHits); }

hcalEnergy()

double TrackDetMatchInfo::hcalEnergy ( )

inline

Definition at line 103 of file TrackDetMatchInfo.h.

References hcalCrossedEnergy().

{ return hcalCrossedEnergy(); }

hcalTowerConeEnergy()

double TrackDetMatchInfo::hcalTowerConeEnergy ( )

inline

Definition at line 98 of file TrackDetMatchInfo.h.

References coneEnergy(), and TowerHcal.

{ return coneEnergy(999, TowerHcal); }

hcalTowerEnergy()

double TrackDetMatchInfo::hcalTowerEnergy ( )

inline

Definition at line 97 of file TrackDetMatchInfo.h.

References crossedEnergy(), and TowerHcal.

{ return crossedEnergy(TowerHcal); }

hoConeEnergy()

double TrackDetMatchInfo::hoConeEnergy

(

)

Definition at line 551 of file TrackDetMatchInfo.cc.

References coneEnergy(), and HcalRecHits.

{ return coneEnergy(999, HcalRecHits); }

hoCrossedEnergy()

double TrackDetMatchInfo::hoCrossedEnergy

(

)

Definition at line 557 of file TrackDetMatchInfo.cc.

References crossedEnergy(), and HORecHits.

Referenced by hoEnergy().

{ return crossedEnergy(HORecHits); }

hoEnergy()

double TrackDetMatchInfo::hoEnergy ( )

inline

Definition at line 104 of file TrackDetMatchInfo.h.

References hoCrossedEnergy().

Referenced by Jet.Jet::hoEnergyFraction().

{ return hoCrossedEnergy(); }

hoTowerConeEnergy()

double TrackDetMatchInfo::hoTowerConeEnergy ( )

inline

Definition at line 100 of file TrackDetMatchInfo.h.

References coneEnergy(), and TowerHO.

{ return coneEnergy(999, TowerHO); }

hoTowerEnergy()

double TrackDetMatchInfo::hoTowerEnergy ( )

inline

Definition at line 99 of file TrackDetMatchInfo.h.

References crossedEnergy(), and TowerHO.

{ return crossedEnergy(TowerHO); }

insideCone()

bool TrackDetMatchInfo::insideCone ( const DetId &  id, const double  dR  )

private

Definition at line 81 of file TrackDetMatchInfo.cc.

References electronAnalyzer_cfi::DeltaR, HGC3DClusterGenMatchSelector_cfi::dR, getPosition(), PV3DBase< T, PVType, FrameType >::mag(), FreeTrajectoryState::momentum(), stateAtIP, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by coneEnergy().

                                                                   {
  GlobalPoint idPosition = getPosition(id);
  if (idPosition.mag() < 0.01)
    return false;

  math::XYZVector idPositionRoot(idPosition.x(), idPosition.y(), idPosition.z());
  math::XYZVector trackP3(stateAtIP.momentum().x(), stateAtIP.momentum().y(), stateAtIP.momentum().z());
  return ROOT::Math::VectorUtil::DeltaR(trackP3, idPositionRoot) < dR;
}

numberOfSegments()

int TrackDetMatchInfo::numberOfSegments

(

)

const

Definition at line 559 of file TrackDetMatchInfo.cc.

References relativeConstraints::chamber, and chambers.

                                              {
  int numSegments = 0;
  for (std::vector<TAMuonChamberMatch>::const_iterator chamber = chambers.begin(); chamber != chambers.end(); chamber++)
    numSegments += chamber->segments.size();
  return numSegments;
}

numberOfSegmentsInDetector()

int TrackDetMatchInfo::numberOfSegmentsInDetector

(

int

detector

)

const

Definition at line 582 of file TrackDetMatchInfo.cc.

References relativeConstraints::chamber, chambers, and hgcalTestNeighbor_cfi::detector.

                                                                    {
  int numSegments = 0;
  for (std::vector<TAMuonChamberMatch>::const_iterator chamber = chambers.begin(); chamber != chambers.end(); chamber++)
    if (chamber->detector() == detector)
      numSegments += chamber->segments.size();
  return numSegments;
}

numberOfSegmentsInStation() [1/2]

int TrackDetMatchInfo::numberOfSegmentsInStation

(

int

station

)

const

Definition at line 566 of file TrackDetMatchInfo.cc.

References relativeConstraints::chamber, chambers, and relativeConstraints::station.

                                                                  {
  int numSegments = 0;
  for (std::vector<TAMuonChamberMatch>::const_iterator chamber = chambers.begin(); chamber != chambers.end(); chamber++)
    if (chamber->station() == station)
      numSegments += chamber->segments.size();
  return numSegments;
}

numberOfSegmentsInStation() [2/2]

int TrackDetMatchInfo::numberOfSegmentsInStation	(	int	station,
		int	detector
	)		const

Definition at line 574 of file TrackDetMatchInfo.cc.

References relativeConstraints::chamber, chambers, hgcalTestNeighbor_cfi::detector, and relativeConstraints::station.

                                                                                {
  int numSegments = 0;
  for (std::vector<TAMuonChamberMatch>::const_iterator chamber = chambers.begin(); chamber != chambers.end(); chamber++)
    if (chamber->station() == station && chamber->detector() == detector)
      numSegments += chamber->segments.size();
  return numSegments;
}

nXnEnergy() [1/2]

double TrackDetMatchInfo::nXnEnergy	(	const DetId &	id,
		EnergyType	type,
		int	gridSize = 1
	)

get energy of the NxN shape (N = 2*gridSize + 1) around given detector element

Definition at line 137 of file TrackDetMatchInfo.cc.

References funct::abs(), DetId::Calo, HLT_2023v12_cff::dEta, HLT_2023v12_cff::dPhi, DetId::Ecal, EcalBarrel, EcalEndcap, EcalRecHits, ecalRecHits, HCALHighEnergyHPDFilter_cfi::energy, Exception, DetId::Hcal, HcalBarrel, HcalEndcap, HcalOuter, HcalRecHits, hcalRecHits, HORecHits, hoRecHits, CaloTowerDetId::ieta(), EBDetId::ieta(), HcalDetId::ieta(), DetIdInfo::info(), CaloTowerDetId::iphi(), EBDetId::iphi(), HcalDetId::iphi(), EEDetId::ix(), EEDetId::iy(), TowerEcal, TowerHcal, TowerHO, towers, TowerTotal, and EEDetId::zside().

Referenced by nXnEnergy().

                                                                                  {
  double energy(0);
  if (id.rawId() == 0)
    return 0.;
  switch (type) {
    case TowerTotal:
    case TowerHcal:
    case TowerEcal:
    case TowerHO: {
      if (id.det() != DetId::Calo) {
        throw cms::Exception("FatalError") << "Wrong DetId. Expected CaloTower, but found:\n"
                                           << DetIdInfo::info(id, nullptr) << "\n";
      }
      CaloTowerDetId centerId(id);
      for (std::vector<const CaloTower*>::const_iterator hit = towers.begin(); hit != towers.end(); hit++) {
        CaloTowerDetId neighborId((*hit)->id());
        int dEta = abs((centerId.ieta() < 0 ? centerId.ieta() + 1 : centerId.ieta()) -
                       (neighborId.ieta() < 0 ? neighborId.ieta() + 1 : neighborId.ieta()));
        int dPhi = abs(centerId.iphi() - neighborId.iphi());
        if (abs(72 - dPhi) < dPhi)
          dPhi = 72 - dPhi;
        if (dEta <= gridSize && dPhi <= gridSize) {
          switch (type) {
            case TowerTotal:
              energy += (*hit)->energy();
              break;
            case TowerEcal:
              energy += (*hit)->emEnergy();
              break;
            case TowerHcal:
              energy += (*hit)->hadEnergy();
              break;
            case TowerHO:
              energy += (*hit)->outerEnergy();
              break;
            default:
              throw cms::Exception("FatalError") << "Unknown calo tower energy type: " << type;
          }
        }
      }
    } break;
    case EcalRecHits: {
      if (id.det() != DetId::Ecal || (id.subdetId() != EcalBarrel && id.subdetId() != EcalEndcap)) {
        throw cms::Exception("FatalError") << "Wrong DetId. Expected EcalBarrel or EcalEndcap, but found:\n"
                                           << DetIdInfo::info(id, nullptr) << "\n";
      }
      // Since the ECAL granularity is small and the gap between EE and EB is significant,
      // energy is computed only within the system that contains the central element
      if (id.subdetId() == EcalBarrel) {
        EBDetId centerId(id);
        for (std::vector<const EcalRecHit*>::const_iterator hit = ecalRecHits.begin(); hit != ecalRecHits.end();
             hit++) {
          if ((*hit)->id().subdetId() != EcalBarrel)
            continue;
          EBDetId neighborId((*hit)->id());
          int dEta = abs((centerId.ieta() < 0 ? centerId.ieta() + 1 : centerId.ieta()) -
                         (neighborId.ieta() < 0 ? neighborId.ieta() + 1 : neighborId.ieta()));
          int dPhi = abs(centerId.iphi() - neighborId.iphi());
          if (abs(360 - dPhi) < dPhi)
            dPhi = 360 - dPhi;
          if (dEta <= gridSize && dPhi <= gridSize) {
            energy += (*hit)->energy();
          }
        }
      } else {
        // Endcap
        EEDetId centerId(id);
        for (std::vector<const EcalRecHit*>::const_iterator hit = ecalRecHits.begin(); hit != ecalRecHits.end();
             hit++) {
          if ((*hit)->id().subdetId() != EcalEndcap)
            continue;
          EEDetId neighborId((*hit)->id());
          if (centerId.zside() == neighborId.zside() && abs(centerId.ix() - neighborId.ix()) <= gridSize &&
              abs(centerId.iy() - neighborId.iy()) <= gridSize) {
            energy += (*hit)->energy();
          }
        }
      }
    } break;
    case HcalRecHits: {
      if (id.det() != DetId::Hcal || (id.subdetId() != HcalBarrel && id.subdetId() != HcalEndcap)) {
        throw cms::Exception("FatalError") << "Wrong DetId. Expected HE or HB, but found:\n"
                                           << DetIdInfo::info(id, nullptr) << "\n";
      }
      HcalDetId centerId(id);
      for (std::vector<const HBHERecHit*>::const_iterator hit = hcalRecHits.begin(); hit != hcalRecHits.end(); hit++) {
        HcalDetId neighborId((*hit)->id());
        int dEta = abs((centerId.ieta() < 0 ? centerId.ieta() + 1 : centerId.ieta()) -
                       (neighborId.ieta() < 0 ? neighborId.ieta() + 1 : neighborId.ieta()));
        int dPhi = abs(centerId.iphi() - neighborId.iphi());
        if (abs(72 - dPhi) < dPhi)
          dPhi = 72 - dPhi;
        if (dEta <= gridSize && dPhi <= gridSize) {
          energy += (*hit)->energy();
        }
      }
    } break;
    case HORecHits: {
      if (id.det() != DetId::Hcal || (id.subdetId() != HcalOuter)) {
        throw cms::Exception("FatalError") << "Wrong DetId. Expected HO, but found:\n"
                                           << DetIdInfo::info(id, nullptr) << "\n";
      }
      HcalDetId centerId(id);
      for (std::vector<const HORecHit*>::const_iterator hit = hoRecHits.begin(); hit != hoRecHits.end(); hit++) {
        HcalDetId neighborId((*hit)->id());
        int dEta = abs((centerId.ieta() < 0 ? centerId.ieta() + 1 : centerId.ieta()) -
                       (neighborId.ieta() < 0 ? neighborId.ieta() + 1 : neighborId.ieta()));
        int dPhi = abs(centerId.iphi() - neighborId.iphi());
        if (abs(72 - dPhi) < dPhi)
          dPhi = 72 - dPhi;
        if (dEta <= gridSize && dPhi <= gridSize) {
          energy += (*hit)->energy();
        }
      }
    } break;
    default:
      throw cms::Exception("FatalError") << "Unkown or not implemented energy type requested, type:" << type;
  }
  return energy;
}

nXnEnergy() [2/2]

double TrackDetMatchInfo::nXnEnergy	(	EnergyType	type,
		int	gridSize = 1
	)

get energy of the NxN shape (N = 2*gridSize + 1) around track projection

Definition at line 258 of file TrackDetMatchInfo.cc.

References crossedEcalIds, crossedHcalIds, crossedHOIds, crossedTowerIds, EcalRecHits, Exception, HcalRecHits, HORecHits, nXnEnergy(), TowerEcal, TowerHcal, TowerHO, and TowerTotal.

                                                                 {
  switch (type) {
    case TowerTotal:
    case TowerHcal:
    case TowerEcal:
    case TowerHO:
      if (crossedTowerIds.empty())
        return 0;
      return nXnEnergy(crossedTowerIds.front(), type, gridSize);
      break;
    case EcalRecHits:
      if (crossedEcalIds.empty())
        return 0;
      return nXnEnergy(crossedEcalIds.front(), type, gridSize);
      break;
    case HcalRecHits:
      if (crossedHcalIds.empty())
        return 0;
      return nXnEnergy(crossedHcalIds.front(), type, gridSize);
      break;
    case HORecHits:
      if (crossedHOIds.empty())
        return 0;
      return nXnEnergy(crossedHOIds.front(), type, gridSize);
      break;
    default:
      throw cms::Exception("FatalError") << "Unkown or not implemented energy type requested, type:" << type;
  }
  return -999;
}

setCaloGeometry()

void TrackDetMatchInfo::setCaloGeometry ( const CaloGeometry *  geometry )

inline

Definition at line 111 of file TrackDetMatchInfo.h.

References caloGeometry, and relativeConstraints::geometry.

{ caloGeometry = geometry; }

Member Data Documentation

caloGeometry

const CaloGeometry* TrackDetMatchInfo::caloGeometry

private

Definition at line 117 of file TrackDetMatchInfo.h.

Referenced by dumpGeometry(), getPosition(), and setCaloGeometry().

chambers

std::vector<TAMuonChamberMatch> TrackDetMatchInfo::chambers

Definition at line 72 of file TrackDetMatchInfo.h.

Referenced by numberOfSegments(), numberOfSegmentsInDetector(), and numberOfSegmentsInStation().

crossedEcalIds

std::vector<DetId> TrackDetMatchInfo::crossedEcalIds

detector elements crossed by a track (regardless of whether energy was deposited or not)

Definition at line 66 of file TrackDetMatchInfo.h.

Referenced by muonisolation::CaloExtractorByAssociator::deposits(), findMaxDeposition(), nXnEnergy(), and pat::PATIsolatedTrackProducer::produce().

crossedEcalRecHits

std::vector<const EcalRecHit*> TrackDetMatchInfo::crossedEcalRecHits

hits in detector elements crossed by a track

Definition at line 59 of file TrackDetMatchInfo.h.

Referenced by crossedEnergy().

crossedHcalIds

std::vector<DetId> TrackDetMatchInfo::crossedHcalIds

Definition at line 67 of file TrackDetMatchInfo.h.

Referenced by muonisolation::CaloExtractorByAssociator::deposits(), findMaxDeposition(), nXnEnergy(), and pat::PATIsolatedTrackProducer::produce().

crossedHcalRecHits

std::vector<const HBHERecHit*> TrackDetMatchInfo::crossedHcalRecHits

Definition at line 60 of file TrackDetMatchInfo.h.

Referenced by crossedEnergy().

crossedHOIds

std::vector<DetId> TrackDetMatchInfo::crossedHOIds

Definition at line 68 of file TrackDetMatchInfo.h.

Referenced by muonisolation::CaloExtractorByAssociator::deposits(), findMaxDeposition(), and nXnEnergy().

crossedHORecHits

std::vector<const HORecHit*> TrackDetMatchInfo::crossedHORecHits

Definition at line 61 of file TrackDetMatchInfo.h.

Referenced by crossedEnergy().

crossedPreshowerIds

std::vector<DetId> TrackDetMatchInfo::crossedPreshowerIds

Definition at line 70 of file TrackDetMatchInfo.h.

crossedTowerIds

std::vector<DetId> TrackDetMatchInfo::crossedTowerIds

Definition at line 69 of file TrackDetMatchInfo.h.

Referenced by muonisolation::JetExtractor::deposit(), muonisolation::CaloExtractorByAssociator::deposits(), findMaxDeposition(), and nXnEnergy().

crossedTowers

std::vector<const CaloTower*> TrackDetMatchInfo::crossedTowers

Definition at line 62 of file TrackDetMatchInfo.h.

Referenced by coneEnergy(), crossedEnergy(), and muonisolation::JetExtractor::deposit().

ecalRecHits

std::vector<const EcalRecHit*> TrackDetMatchInfo::ecalRecHits

hits in the cone

Definition at line 53 of file TrackDetMatchInfo.h.

Referenced by coneEnergy(), muonisolation::CaloExtractorByAssociator::deposits(), findMaxDeposition(), and nXnEnergy().

ecalTrueEnergy

double TrackDetMatchInfo::ecalTrueEnergy

Definition at line 79 of file TrackDetMatchInfo.h.

hcalRecHits

std::vector<const HBHERecHit*> TrackDetMatchInfo::hcalRecHits

Definition at line 54 of file TrackDetMatchInfo.h.

Referenced by coneEnergy(), muonisolation::CaloExtractorByAssociator::deposits(), findMaxDeposition(), and nXnEnergy().

hcalTrueEnergy

double TrackDetMatchInfo::hcalTrueEnergy

Definition at line 80 of file TrackDetMatchInfo.h.

hcalTrueEnergyCorrected

double TrackDetMatchInfo::hcalTrueEnergyCorrected

Definition at line 81 of file TrackDetMatchInfo.h.

hoRecHits

std::vector<const HORecHit*> TrackDetMatchInfo::hoRecHits

Definition at line 55 of file TrackDetMatchInfo.h.

Referenced by coneEnergy(), muonisolation::CaloExtractorByAssociator::deposits(), findMaxDeposition(), and nXnEnergy().

isGoodCalo

bool TrackDetMatchInfo::isGoodCalo

Definition at line 48 of file TrackDetMatchInfo.h.

isGoodEcal

bool TrackDetMatchInfo::isGoodEcal

Definition at line 46 of file TrackDetMatchInfo.h.

Referenced by spr::chargeIsolation(), and spr::chargeIsolationEcal().

isGoodHcal

bool TrackDetMatchInfo::isGoodHcal

Definition at line 47 of file TrackDetMatchInfo.h.

Referenced by spr::chargeIsolationHcal(), and spr::coneChargeIsolation().

isGoodHO

bool TrackDetMatchInfo::isGoodHO

Definition at line 49 of file TrackDetMatchInfo.h.

isGoodMuon

bool TrackDetMatchInfo::isGoodMuon

Definition at line 50 of file TrackDetMatchInfo.h.

simTrack

const SimTrack* TrackDetMatchInfo::simTrack

MC truth info.

Definition at line 78 of file TrackDetMatchInfo.h.

simTrackRef_

SimTrackRef TrackDetMatchInfo::simTrackRef_

Definition at line 86 of file TrackDetMatchInfo.h.

stateAtIP

FreeTrajectoryState TrackDetMatchInfo::stateAtIP

track info

Definition at line 75 of file TrackDetMatchInfo.h.

Referenced by insideCone().

towers

std::vector<const CaloTower*> TrackDetMatchInfo::towers

Definition at line 56 of file TrackDetMatchInfo.h.

Referenced by muonisolation::CaloExtractorByAssociator::deposits(), findMaxDeposition(), and nXnEnergy().

trackRef_

reco::TrackRef TrackDetMatchInfo::trackRef_

Obsolete methods and data members for backward compatibility. Will be removed in future releases.

Definition at line 85 of file TrackDetMatchInfo.h.

trkGlobPosAtEcal

math::XYZPoint TrackDetMatchInfo::trkGlobPosAtEcal

Track position at different parts of the calorimeter.

Definition at line 38 of file TrackDetMatchInfo.h.

Referenced by EopTreeWriter::analyze(), IsolatedTracksCone::analyze(), spr::chargeIsolation(), spr::chargeIsolationEcal(), muonisolation::CaloExtractorByAssociator::deposits(), and pat::PATIsolatedTrackProducer::produce().

trkGlobPosAtHcal

math::XYZPoint TrackDetMatchInfo::trkGlobPosAtHcal

Definition at line 39 of file TrackDetMatchInfo.h.

Referenced by IsolatedTracksCone::analyze(), spr::chargeIsolationHcal(), spr::coneChargeIsolation(), muonisolation::JetExtractor::deposit(), and muonisolation::CaloExtractorByAssociator::deposits().

trkGlobPosAtHO

math::XYZPoint TrackDetMatchInfo::trkGlobPosAtHO

Definition at line 40 of file TrackDetMatchInfo.h.

Referenced by muonisolation::CaloExtractorByAssociator::deposits().

trkMomAtEcal

GlobalVector TrackDetMatchInfo::trkMomAtEcal

Definition at line 42 of file TrackDetMatchInfo.h.

Referenced by IsolatedTracksCone::analyze().

trkMomAtHcal

GlobalVector TrackDetMatchInfo::trkMomAtHcal

Definition at line 43 of file TrackDetMatchInfo.h.

Referenced by IsolatedTracksCone::analyze().

trkMomAtHO

GlobalVector TrackDetMatchInfo::trkMomAtHO

Definition at line 44 of file TrackDetMatchInfo.h.