#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.

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

References alongMomentum, dir_, forwardPropagator_, and theMF_.

◆ ~ConversionTrackEcalImpactPoint()

ConversionTrackEcalImpactPoint::~ConversionTrackEcalImpactPoint ( )

Definition at line 54 of file ConversionTrackEcalImpactPoint.cc.

54 { delete forwardPropagator_; }

References forwardPropagator_.

Member Function Documentation

◆ barrel()

static const BoundCylinder& ConversionTrackEcalImpactPoint::barrel ( )

inlinestaticprivate

Definition at line 62 of file ConversionTrackEcalImpactPoint.h.

62 { return *theBarrel_; }

References theBarrel_.

Referenced by find().

◆ 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.

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