#include <ConversionTrackEcalImpactPoint.h>

Public Member Functions
	ConversionTrackEcalImpactPoint (const MagneticField *field)

std::vector< math::XYZPointF >	find (const std::vector< reco::TransientTrack > &tracks, const edm::Handle< edm::View< reco::CaloCluster > > &bcHandle)

std::vector< reco::CaloClusterPtr >	matchingBC () const

void	setMagneticField (const MagneticField *mf)

	~ConversionTrackEcalImpactPoint ()

Static Private Member Functions
static const BoundCylinder &	barrel ()

static const BoundDisk &	negativeEtaEndcap ()

static const BoundDisk &	positiveEtaEndcap ()

Private Attributes
PropagationDirection	dir_

PropagatorWithMaterial *	forwardPropagator_

std::vector< reco::CaloClusterPtr >	matchingBC_

TrajectoryStateOnSurface	stateAtECAL_

const MagneticField *	theMF_

Static Private Attributes
static const ReferenceCountingPointer< BoundCylinder >	theBarrel_ = initBarrel()

static const ReferenceCountingPointer< BoundDisk >	theNegativeEtaEndcap_ = initNegative()

static const ReferenceCountingPointer< BoundDisk >	thePositiveEtaEndcap_ = initPositive()

Detailed Description

Author: N. Marinelli - Univ. of Notre Dame

Version

Definition at line 37 of file ConversionTrackEcalImpactPoint.h.

Constructor & Destructor Documentation

◆ ConversionTrackEcalImpactPoint()

ConversionTrackEcalImpactPoint::ConversionTrackEcalImpactPoint ( const MagneticField * field )

Definition at line 47 of file ConversionTrackEcalImpactPoint.cc.

References alongMomentum, dir_, forwardPropagator_, and theMF_.

     :
 
       theMF_(field) {
   forwardPropagator_ = new PropagatorWithMaterial(dir_ = alongMomentum, 0.000511, theMF_);
 }

◆ ~ConversionTrackEcalImpactPoint()

ConversionTrackEcalImpactPoint::~ConversionTrackEcalImpactPoint ( )

Definition at line 54 of file ConversionTrackEcalImpactPoint.cc.

References forwardPropagator_.

54 { delete forwardPropagator_; }

ConversionTrackEcalImpactPoint::forwardPropagator_

PropagatorWithMaterial * forwardPropagator_

Definition: ConversionTrackEcalImpactPoint.h:54

Member Function Documentation

◆ barrel()

static const BoundCylinder& ConversionTrackEcalImpactPoint::barrel ( )

inlinestaticprivate

Definition at line 62 of file ConversionTrackEcalImpactPoint.h.

References theBarrel_.

Referenced by find().

62 { return *theBarrel_; }

ConversionTrackEcalImpactPoint::theBarrel_

static const ReferenceCountingPointer< BoundCylinder > theBarrel_

Definition: ConversionTrackEcalImpactPoint.h:58

◆ find()

std::vector< math::XYZPointF > ConversionTrackEcalImpactPoint::find	(	const std::vector< reco::TransientTrack > &	tracks,
		const edm::Handle< edm::View< reco::CaloCluster > > &	bcHandle
	)

Definition at line 56 of file ConversionTrackEcalImpactPoint.cc.

References barrel(), HLT_2022v15_cff::dEta, HLT_2022v15_cff::dPhi, PV3DBase< T, PVType, FrameType >::eta(), forwardPropagator_, TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), matchingBC(), matchingBC_, negativeEtaEndcap(), positiveEtaEndcap(), PropagatorWithMaterial::propagate(), mps_fire::result, mathSSE::sqrt(), stateAtECAL_, HLT_2022v15_cff::track, and pwdgSkimBPark_cfi::tracks.

Referenced by ConvertedPhotonProducer::buildCollections().

                                                                                                        {
   std::vector<math::XYZPointF> result;
   //
   matchingBC_.clear();
 
   std::vector<reco::CaloClusterPtr> matchingBC(2);
 
   //
 
   int iTrk = 0;
   for (auto const& track : tracks) {
     math::XYZPointF ecalImpactPosition(0., 0., 0.);
     const TrajectoryStateOnSurface myTSOS = track.innermostMeasurementState();
     if (!(myTSOS.isValid()))
       continue;
 
     stateAtECAL_ = forwardPropagator_->propagate(myTSOS, barrel());
 
     if (!stateAtECAL_.isValid() || (stateAtECAL_.isValid() && fabs(stateAtECAL_.globalPosition().eta()) > 1.479)) {
       if (track.innermostMeasurementState().globalPosition().eta() > 0.) {
         stateAtECAL_ = forwardPropagator_->propagate(myTSOS, positiveEtaEndcap());
 
       } else {
         stateAtECAL_ = forwardPropagator_->propagate(myTSOS, negativeEtaEndcap());
       }
     }
 
     if (stateAtECAL_.isValid())
       ecalImpactPosition = stateAtECAL_.globalPosition();
 
     result.push_back(ecalImpactPosition);
 
     if (stateAtECAL_.isValid()) {
       int goodBC = 0;
       float bcDistanceToTrack = 9999;
       reco::CaloClusterPtr matchedBCItr;
       int ibc = 0;
       goodBC = 0;
 
       for (auto const& bc : *bcHandle) {
         float dEta = bc.position().eta() - ecalImpactPosition.eta();
         float dPhi = bc.position().phi() - ecalImpactPosition.phi();
         if (sqrt(dEta * dEta + dPhi * dPhi) < bcDistanceToTrack) {
           goodBC = ibc;
           bcDistanceToTrack = sqrt(dEta * dEta + dPhi * dPhi);
         }
         ibc++;
       }
 
       matchingBC[iTrk] = bcHandle->ptrAt(goodBC);
     }
 
     iTrk++;
   }
 
   matchingBC_ = matchingBC;
 
   return result;
 }