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 (edm::ESHandle< 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
edm::ESHandle< CaloGeometry >	caloGeometry

Detailed Description

Definition at line 14 of file TrackDetMatchInfo.h.

Member Enumeration Documentation

enum TrackDetMatchInfo::EnergyType

Enumerator
EcalRecHits
HcalRecHits
HORecHits
TowerTotal
TowerEcal
TowerHcal
TowerHO

Definition at line 16 of file TrackDetMatchInfo.h.

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

Constructor & Destructor Documentation

TrackDetMatchInfo::TrackDetMatchInfo

(

)

Definition at line 314 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(0)
     , ecalTrueEnergy(-999)
     , hcalTrueEnergy(-999)
{
}

Member Function Documentation

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 108 of file TrackDetMatchInfo.cc.

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

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

{
   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;
}

double TrackDetMatchInfo::crossedEnergy

(

EnergyType

type

)

energy in detector elements crossed by the track by types

Definition at line 47 of file TrackDetMatchInfo.cc.

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

Referenced by BetaCalculatorECAL::addInfoToCandidate(), ecalCrossedEnergy(), ecalTowerEnergy(), MuonIdProducer::fillMuonId(), 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;
}

std::string TrackDetMatchInfo::dumpGeometry

(

const DetId &

id

)

Definition at line 17 of file TrackDetMatchInfo.cc.

References EZArrayFL< T >::begin(), caloGeometry, EZArrayFL< T >::end(), Exception, CaloCellGeometry::getCorners(), CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), edm::ESHandleBase::isValid(), point, and hiPixelPairStep_cff::points.

Referenced by setCaloGeometry().

{
   if ( ! caloGeometry.isValid() || 
    ! 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();
}

double TrackDetMatchInfo::ecalConeEnergy

(

)

Definition at line 581 of file TrackDetMatchInfo.cc.

References coneEnergy(), and EcalRecHits.

{
   return coneEnergy (999, EcalRecHits );
}

double TrackDetMatchInfo::ecalCrossedEnergy

(

)

Definition at line 596 of file TrackDetMatchInfo.cc.

References crossedEnergy(), and EcalRecHits.

Referenced by ecalEnergy().

{
   return crossedEnergy( EcalRecHits );
}

double TrackDetMatchInfo::ecalEnergy ( )

inline

Definition at line 103 of file TrackDetMatchInfo.h.

References ecalCrossedEnergy().

{ return ecalCrossedEnergy(); }

double TrackDetMatchInfo::ecalTowerConeEnergy ( )

inline

Definition at line 97 of file TrackDetMatchInfo.h.

References coneEnergy(), and TowerEcal.

{ return coneEnergy(999,TowerEcal); }

double TrackDetMatchInfo::ecalTowerEnergy ( )

inline

Definition at line 96 of file TrackDetMatchInfo.h.

References crossedEnergy(), and TowerEcal.

{ return crossedEnergy(TowerEcal); }

DetId TrackDetMatchInfo::findMaxDeposition

(

EnergyType

type

)

Find detector elements with highest energy deposition.

Definition at line 332 of file TrackDetMatchInfo.cc.

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

Referenced by MuonIdProducer::fillMuonId(), findMaxDeposition(), and InterestingTrackEcalDetIdProducer::produce().

{
   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:
    {
       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();
        }
         }
    }
    default:
      throw cms::Exception("FatalError") << "Maximal energy deposition: unkown or not implemented energy type requested, type:" << type;
   }
   return id;
}

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

Definition at line 546 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;
}

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

Definition at line 399 of file TrackDetMatchInfo.cc.

References funct::abs(), DetId::Calo, 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,0)<<"\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,0)<<"\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,0)<<"\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,0)<<"\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;
}

GlobalPoint TrackDetMatchInfo::getPosition

(

const DetId &

id

)

Definition at line 34 of file TrackDetMatchInfo.cc.

References caloGeometry, Exception, CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), and edm::ESHandleBase::isValid().

Referenced by BetaCalculatorECAL::addInfoToCandidate(), insideCone(), and setCaloGeometry().

{
   // this part might be slow
   if ( ! caloGeometry.isValid() || 
    ! 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();
}

double TrackDetMatchInfo::hcalConeEnergy

(

)

Definition at line 586 of file TrackDetMatchInfo.cc.

References coneEnergy(), and HcalRecHits.

{
   return coneEnergy (999, HcalRecHits );
}

double TrackDetMatchInfo::hcalCrossedEnergy

(

)

Definition at line 601 of file TrackDetMatchInfo.cc.

References crossedEnergy(), and HcalRecHits.

Referenced by hcalEnergy().

{
   return crossedEnergy( HcalRecHits );
}

double TrackDetMatchInfo::hcalEnergy ( )

inline

Definition at line 104 of file TrackDetMatchInfo.h.

References hcalCrossedEnergy().

{ return hcalCrossedEnergy(); }

double TrackDetMatchInfo::hcalTowerConeEnergy ( )

inline

Definition at line 99 of file TrackDetMatchInfo.h.

References coneEnergy(), and TowerHcal.

{ return coneEnergy(999, TowerHcal); }

double TrackDetMatchInfo::hcalTowerEnergy ( )

inline

Definition at line 98 of file TrackDetMatchInfo.h.

References crossedEnergy(), and TowerHcal.

{ return crossedEnergy(TowerHcal); }

double TrackDetMatchInfo::hoConeEnergy

(

)

Definition at line 591 of file TrackDetMatchInfo.cc.

References coneEnergy(), and HcalRecHits.

{
   return coneEnergy (999, HcalRecHits );
}

double TrackDetMatchInfo::hoCrossedEnergy

(

)

Definition at line 606 of file TrackDetMatchInfo.cc.

References crossedEnergy(), and HORecHits.

Referenced by hoEnergy().

{
   return crossedEnergy( HORecHits );
}

double TrackDetMatchInfo::hoEnergy ( )

inline

Definition at line 105 of file TrackDetMatchInfo.h.

References gamEcalExtractorBlocks_cff::detector, hoCrossedEnergy(), numberOfSegments(), numberOfSegmentsInDetector(), numberOfSegmentsInStation(), and relativeConstraints::station.

Referenced by Jet.Jet::hoEnergyFraction().

{ return hoCrossedEnergy(); }

double TrackDetMatchInfo::hoTowerConeEnergy ( )

inline

Definition at line 101 of file TrackDetMatchInfo.h.

References coneEnergy().

{ return coneEnergy(999, TowerHO); }

double TrackDetMatchInfo::hoTowerEnergy ( )

inline

Definition at line 100 of file TrackDetMatchInfo.h.

References crossedEnergy().

{ return crossedEnergy(TowerHO); }

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

private

Definition at line 99 of file TrackDetMatchInfo.cc.

References HiRegitMuonDetachedTripletStep_cff::DeltaR, PFRecoTauDiscriminationAgainstElectronDeadECAL_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(), and setCaloGeometry().

                                                                   {
   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;
}

int TrackDetMatchInfo::numberOfSegments

(

)

const

Definition at line 612 of file TrackDetMatchInfo.cc.

References relativeConstraints::chamber, and chambers.

Referenced by hoEnergy().

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

int TrackDetMatchInfo::numberOfSegmentsInDetector

(

int

detector

)

const

Definition at line 633 of file TrackDetMatchInfo.cc.

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

Referenced by hoEnergy().

                                                                    {
   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;
}

int TrackDetMatchInfo::numberOfSegmentsInStation

(

int

station

)

const

Definition at line 619 of file TrackDetMatchInfo.cc.

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

Referenced by hoEnergy().

                                                                  {
   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;
}

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

Definition at line 626 of file TrackDetMatchInfo.cc.

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

                                                                                {
   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;
}

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 163 of file TrackDetMatchInfo.cc.

References funct::abs(), DetId::Calo, 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().

Referenced by BetaCalculatorECAL::addInfoToCandidate(), EopTreeWriter::analyze(), MuonIdProducer::fillMuonId(), and 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,0)<<"\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,0)<<"\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,0)<<"\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,0)<<"\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;
}

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

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

Definition at line 285 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;
}

void TrackDetMatchInfo::setCaloGeometry ( edm::ESHandle< CaloGeometry >  geometry )

inline

Definition at line 112 of file TrackDetMatchInfo.h.

References caloGeometry, dumpGeometry(), geometry, getPosition(), insideCone(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by TrackDetectorAssociator::associate().

{ caloGeometry = geometry; }

Member Data Documentation

edm::ESHandle<CaloGeometry> TrackDetMatchInfo::caloGeometry

private

Definition at line 117 of file TrackDetMatchInfo.h.

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

std::vector<TAMuonChamberMatch> TrackDetMatchInfo::chambers

Definition at line 73 of file TrackDetMatchInfo.h.

Referenced by TrackDetectorAssociator::fillMuon(), MuonIdProducer::fillMuonId(), numberOfSegments(), numberOfSegmentsInDetector(), and numberOfSegmentsInStation().

std::vector<DetId> TrackDetMatchInfo::crossedEcalIds

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

Definition at line 67 of file TrackDetMatchInfo.h.

Referenced by EcalCosmicsHists::analyze(), muonisolation::CaloExtractorByAssociator::deposits(), TrackDetectorAssociator::fillEcal(), MuonIdProducer::fillMuonId(), findMaxDeposition(), nXnEnergy(), pat::PATIsolatedTrackProducer::produce(), and HighPtTrackEcalDetIdProducer::produce().

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

hits in detector elements crossed by a track

Definition at line 60 of file TrackDetMatchInfo.h.

Referenced by crossedEnergy(), CaloCleaner< T >::fill_correction_map(), and TrackDetectorAssociator::fillEcal().

std::vector<DetId> TrackDetMatchInfo::crossedHcalIds

Definition at line 68 of file TrackDetMatchInfo.h.

Referenced by muonisolation::CaloExtractorByAssociator::deposits(), TrackDetectorAssociator::fillHcal(), MuonIdProducer::fillMuonId(), findMaxDeposition(), nXnEnergy(), pat::PATIsolatedTrackProducer::produce(), and ReduceHcalRecHitCollectionProducer::produce().

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

Definition at line 61 of file TrackDetMatchInfo.h.

Referenced by BetaCalculatorECAL::addInfoToCandidate(), crossedEnergy(), CaloCleaner< T >::fill_correction_map(), and TrackDetectorAssociator::fillHcal().

std::vector<DetId> TrackDetMatchInfo::crossedHOIds

Definition at line 69 of file TrackDetMatchInfo.h.

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

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

Definition at line 62 of file TrackDetMatchInfo.h.

Referenced by crossedEnergy(), CaloCleaner< T >::fill_correction_map(), and TrackDetectorAssociator::fillHO().

std::vector<DetId> TrackDetMatchInfo::crossedPreshowerIds

Definition at line 71 of file TrackDetMatchInfo.h.

Referenced by CaloCleaner< T >::fill_correction_map(), and TrackDetectorAssociator::fillPreshower().

std::vector<DetId> TrackDetMatchInfo::crossedTowerIds

Definition at line 70 of file TrackDetMatchInfo.h.

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

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

Definition at line 63 of file TrackDetMatchInfo.h.

Referenced by coneEnergy(), crossedEnergy(), muonisolation::JetExtractor::deposit(), TrackDetectorAssociator::fillCaloTowers(), and MuonMETAlgo::GetMuDepDeltas().

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

hits in the cone

Definition at line 54 of file TrackDetMatchInfo.h.

Referenced by coneEnergy(), muonisolation::CaloExtractorByAssociator::deposits(), TrackDetectorAssociator::fillEcal(), MuonIdProducer::fillMuonId(), findMaxDeposition(), nXnEnergy(), and ProduceIsolationMap::produce().

double TrackDetMatchInfo::ecalTrueEnergy

Definition at line 80 of file TrackDetMatchInfo.h.

Referenced by TrackDetectorAssociator::fillCaloTruth().

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

Definition at line 55 of file TrackDetMatchInfo.h.

Referenced by coneEnergy(), muonisolation::CaloExtractorByAssociator::deposits(), TrackDetectorAssociator::fillHcal(), MuonIdProducer::fillMuonId(), findMaxDeposition(), nXnEnergy(), ReduceHcalRecHitCollectionProducer::produce(), and ProduceIsolationMap::produce().

double TrackDetMatchInfo::hcalTrueEnergy

Definition at line 81 of file TrackDetMatchInfo.h.

Referenced by TrackDetectorAssociator::fillCaloTruth().

double TrackDetMatchInfo::hcalTrueEnergyCorrected

Definition at line 82 of file TrackDetMatchInfo.h.

Referenced by TrackDetectorAssociator::fillCaloTruth().

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

Definition at line 56 of file TrackDetMatchInfo.h.

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

bool TrackDetMatchInfo::isGoodCalo

Definition at line 49 of file TrackDetMatchInfo.h.

Referenced by TrackDetectorAssociator::fillCaloTowers().

bool TrackDetMatchInfo::isGoodEcal

Definition at line 47 of file TrackDetMatchInfo.h.

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

bool TrackDetMatchInfo::isGoodHcal

Definition at line 48 of file TrackDetMatchInfo.h.

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

bool TrackDetMatchInfo::isGoodHO

Definition at line 50 of file TrackDetMatchInfo.h.

Referenced by TrackDetectorAssociator::fillHO().

bool TrackDetMatchInfo::isGoodMuon

Definition at line 51 of file TrackDetMatchInfo.h.

const SimTrack* TrackDetMatchInfo::simTrack

MC truth info.

Definition at line 79 of file TrackDetMatchInfo.h.

Referenced by TrackDetectorAssociator::fillCaloTruth().

SimTrackRef TrackDetMatchInfo::simTrackRef_

Definition at line 87 of file TrackDetMatchInfo.h.

FreeTrajectoryState TrackDetMatchInfo::stateAtIP

track info

Definition at line 76 of file TrackDetMatchInfo.h.

Referenced by TrackDetectorAssociator::associate(), TrackDetectorAssociator::fillCaloTruth(), and insideCone().

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

Definition at line 57 of file TrackDetMatchInfo.h.

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

reco::TrackRef TrackDetMatchInfo::trackRef_

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

Definition at line 86 of file TrackDetMatchInfo.h.

math::XYZPoint TrackDetMatchInfo::trkGlobPosAtEcal

Track position at different parts of the calorimeter.

Definition at line 39 of file TrackDetMatchInfo.h.

Referenced by EopTreeWriter::analyze(), ValidIsoTrkCalib::analyze(), IsolatedTracksCone::analyze(), TrackDetectorAssociator::associate(), spr::chargeIsolation(), spr::chargeIsolationEcal(), muonisolation::CaloExtractorByAssociator::deposits(), MuonIdProducer::fillMuonId(), MuonMETAlgo::GetMuDepDeltas(), and pat::PATIsolatedTrackProducer::produce().

math::XYZPoint TrackDetMatchInfo::trkGlobPosAtHcal

Definition at line 40 of file TrackDetMatchInfo.h.

Referenced by ValidIsoTrkCalib::analyze(), IsolatedTracksCone::analyze(), TrackDetectorAssociator::associate(), spr::chargeIsolationHcal(), spr::coneChargeIsolation(), muonisolation::JetExtractor::deposit(), muonisolation::CaloExtractorByAssociator::deposits(), TrackDetectorAssociator::fillCaloTruth(), MuonIdProducer::fillMuonId(), and MuonMETAlgo::GetMuDepDeltas().

math::XYZPoint TrackDetMatchInfo::trkGlobPosAtHO

Definition at line 41 of file TrackDetMatchInfo.h.

Referenced by TrackDetectorAssociator::associate(), muonisolation::CaloExtractorByAssociator::deposits(), and MuonMETAlgo::GetMuDepDeltas().

GlobalVector TrackDetMatchInfo::trkMomAtEcal

Definition at line 43 of file TrackDetMatchInfo.h.

Referenced by IsolatedTracksCone::analyze(), and TrackDetectorAssociator::associate().

GlobalVector TrackDetMatchInfo::trkMomAtHcal

Definition at line 44 of file TrackDetMatchInfo.h.

Referenced by IsolatedTracksCone::analyze(), and TrackDetectorAssociator::associate().

GlobalVector TrackDetMatchInfo::trkMomAtHO

Definition at line 45 of file TrackDetMatchInfo.h.

Referenced by TrackDetectorAssociator::associate().