#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 39 of file ConversionTrackEcalImpactPoint.h.

Constructor & Destructor Documentation

ConversionTrackEcalImpactPoint::ConversionTrackEcalImpactPoint ( const MagneticField * field )

Definition at line 53 of file ConversionTrackEcalImpactPoint.cc.

References alongMomentum, dir_, forwardPropagator_, and theMF_.

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

ConversionTrackEcalImpactPoint::~ConversionTrackEcalImpactPoint ( )

Definition at line 62 of file ConversionTrackEcalImpactPoint.cc.

References forwardPropagator_.

                                                                 {
 
     delete    forwardPropagator_ ; 
     
 }

Member Function Documentation

static const BoundCylinder& ConversionTrackEcalImpactPoint::barrel ( )

inlinestaticprivate

Definition at line 72 of file ConversionTrackEcalImpactPoint.h.

References theBarrel_.

Referenced by find().

72 { return *theBarrel_;}

ConversionTrackEcalImpactPoint::theBarrel_

static const ReferenceCountingPointer< BoundCylinder > theBarrel_

Definition: ConversionTrackEcalImpactPoint.h:68

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

Definition at line 68 of file ConversionTrackEcalImpactPoint.cc.

References barrel(), particleFlow_cfi::dEta, particleFlow_cfi::dPhi, PV3DBase< T, PVType, FrameType >::eta(), forwardPropagator_, TrajectoryStateOnSurface::globalPosition(), mps_fire::i, TrajectoryStateOnSurface::isValid(), matchingBC(), matchingBC_, negativeEtaEndcap(), positiveEtaEndcap(), Propagator::propagate(), mps_fire::result, mathSSE::sqrt(), and stateAtECAL_.

Referenced by ConvertedPhotonProducer::buildCollections().

                                                                                                                                                                       {
 
   std::vector<math::XYZPointF> result;
   // 
   matchingBC_.clear();   
 
   std::vector<reco::CaloClusterPtr> matchingBC(2);
 
   // 
 
   int iTrk=0;
   for (    std::vector<reco::TransientTrack>::const_iterator iTk=tracks.begin(); iTk!=tracks.end(); ++iTk) {
 
     math::XYZPointF ecalImpactPosition(0.,0.,0.);
     const TrajectoryStateOnSurface myTSOS=(*iTk).innermostMeasurementState();
     if ( !( myTSOS.isValid() ) ) continue; 
 
     stateAtECAL_= forwardPropagator_->propagate( myTSOS, barrel() );
     
     if (!stateAtECAL_.isValid() || ( stateAtECAL_.isValid() && fabs(stateAtECAL_.globalPosition().eta() ) >1.479 )  ) {
     
       if ( (*iTk).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 (unsigned i = 0; i < bcHandle->size(); ++i ) {
     float dEta= bcHandle->ptrAt(i)->position().eta() - ecalImpactPosition.eta()  ;
     float dPhi= bcHandle->ptrAt(i)->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;
 
 }