ConversionTrackEcalImpactPoint Class Reference

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

{ return *theBarrel_;}

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(), TrajectoryStateOnSurface::isValid(), matchingBC(), matchingBC_, negativeEtaEndcap(), positiveEtaEndcap(), Propagator::propagate(), mps_fire::result, mathSSE::sqrt(), stateAtECAL_, and HiIsolationCommonParameters_cff::track.

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;

}

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

inline

Definition at line 51 of file ConversionTrackEcalImpactPoint.h.

References matchingBC_.

Referenced by ConvertedPhotonProducer::buildCollections(), and find().

{return matchingBC_;}

static const BoundDisk& ConversionTrackEcalImpactPoint::negativeEtaEndcap ( )

inlinestaticprivate

Definition at line 73 of file ConversionTrackEcalImpactPoint.h.

References theNegativeEtaEndcap_.

Referenced by find().

{ return *theNegativeEtaEndcap_;}

static const BoundDisk& ConversionTrackEcalImpactPoint::positiveEtaEndcap ( )

inlinestaticprivate

Definition at line 74 of file ConversionTrackEcalImpactPoint.h.

References thePositiveEtaEndcap_.

Referenced by find().

{ return *thePositiveEtaEndcap_;}

void ConversionTrackEcalImpactPoint::setMagneticField ( const MagneticField *  mf )

inline

Definition at line 52 of file ConversionTrackEcalImpactPoint.h.

References theMF_.

{ theMF_=mf;}

Member Data Documentation

PropagationDirection ConversionTrackEcalImpactPoint::dir_

private

Definition at line 65 of file ConversionTrackEcalImpactPoint.h.

Referenced by ConversionTrackEcalImpactPoint().

PropagatorWithMaterial* ConversionTrackEcalImpactPoint::forwardPropagator_

private

Definition at line 64 of file ConversionTrackEcalImpactPoint.h.

Referenced by ConversionTrackEcalImpactPoint(), find(), and ~ConversionTrackEcalImpactPoint().

std::vector<reco::CaloClusterPtr> ConversionTrackEcalImpactPoint::matchingBC_

private

Definition at line 66 of file ConversionTrackEcalImpactPoint.h.

Referenced by find(), and matchingBC().

TrajectoryStateOnSurface ConversionTrackEcalImpactPoint::stateAtECAL_

private

Definition at line 62 of file ConversionTrackEcalImpactPoint.h.

Referenced by find().

const ReferenceCountingPointer< BoundCylinder > ConversionTrackEcalImpactPoint::theBarrel_ = initBarrel()

staticprivate

Definition at line 68 of file ConversionTrackEcalImpactPoint.h.

Referenced by barrel(), and initPositive().

const MagneticField* ConversionTrackEcalImpactPoint::theMF_

private

Definition at line 60 of file ConversionTrackEcalImpactPoint.h.

Referenced by ConversionTrackEcalImpactPoint(), and setMagneticField().

const ReferenceCountingPointer< BoundDisk > ConversionTrackEcalImpactPoint::theNegativeEtaEndcap_ = initNegative()

staticprivate

Definition at line 69 of file ConversionTrackEcalImpactPoint.h.

Referenced by initPositive(), and negativeEtaEndcap().

const ReferenceCountingPointer< BoundDisk > ConversionTrackEcalImpactPoint::thePositiveEtaEndcap_ = initPositive()

staticprivate

Definition at line 70 of file ConversionTrackEcalImpactPoint.h.

Referenced by initPositive(), and positiveEtaEndcap().