#include <EgammaHLTPixelMatchElectronAlgo.h>

Public Member Functions
	EgammaHLTPixelMatchElectronAlgo (const edm::ParameterSet &conf)

void	run (edm::Event &, reco::ElectronCollection &)

void	setupES (const edm::EventSetup &setup)

	~EgammaHLTPixelMatchElectronAlgo ()

Private Member Functions
void	process (edm::Handle< reco::TrackCollection > tracksH, reco::ElectronCollection &outEle, Global3DPoint &bs)

Private Attributes
edm::InputTag	BSProducer_

edm::ESHandle < GeometricSearchTracker >	theGeomSearchTracker

edm::ESHandle< MagneticField >	theMagField

edm::InputTag	trackProducer_

Detailed Description

Class to reconstruct electron tracks from electron pixel seeds keep track of information about the initiating supercluster

Author: Monica Vazquez Acosta (CERN)

Description: top algorithm producing TrackCandidate and Electron objects from supercluster driven pixel seeded Ckf tracking for HLT

Definition at line 76 of file EgammaHLTPixelMatchElectronAlgo.h.

Constructor & Destructor Documentation

EgammaHLTPixelMatchElectronAlgo::EgammaHLTPixelMatchElectronAlgo ( const edm::ParameterSet & conf )

Definition at line 52 of file EgammaHLTPixelMatchElectronAlgo.cc.

                                                                                             :
   trackProducer_( conf.getParameter<edm::InputTag>("TrackProducer") ),
   BSProducer_(    conf.getParameter<edm::InputTag>("BSProducer") )
 {}

EgammaHLTPixelMatchElectronAlgo::~EgammaHLTPixelMatchElectronAlgo ( )

Definition at line 57 of file EgammaHLTPixelMatchElectronAlgo.cc.

                                                                   {
 
   // delete theInitialStateEstimator;
   //delete theNavigationSchool;
   //delete theTrajectoryCleaner; 
     
 }

Member Function Documentation

void EgammaHLTPixelMatchElectronAlgo::process	(	edm::Handle< reco::TrackCollection >	tracksH,
		reco::ElectronCollection &	outEle,
		Global3DPoint &	bs
	)

private

Definition at line 90 of file EgammaHLTPixelMatchElectronAlgo.cc.

References reco::TrackBase::charge(), i, TrajectoryStateTransform::innerFreeState(), PV3DBase< T, PVType, FrameType >::mag(), TrajectoryStateClosestToPoint::momentum(), edm::Handle< T >::product(), edm::ESHandle< class >::product(), reco::Electron::setSuperCluster(), matplotRender::t, theMagField, testEve_cfg::tracks, reco::TrackBase::vertex(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by run().

                                                                                                                                  {
   const TrackCollection *tracks=tracksH.product();
   for (unsigned int i=0;i<tracks->size();++i) {
     const Track & t=(*tracks)[i];
 
     const TrackRef trackRef = edm::Ref<TrackCollection>(tracksH,i);
     edm::RefToBase<TrajectorySeed> seed = trackRef->extra()->seedRef();
     ElectronSeedRef elseed=seed.castTo<ElectronSeedRef>();
 
     edm::RefToBase<CaloCluster> caloCluster = elseed->caloCluster() ;
     SuperClusterRef scRef = caloCluster.castTo<SuperClusterRef>() ;
 
     //const SuperClusterRef & scRef=elseed->superCluster();
     
         // Get the momentum at vertex (not at the innermost layer)
     TSCPBuilderNoMaterial tscpBuilder;
     TrajectoryStateTransform tsTransform;
     FreeTrajectoryState fts = tsTransform.innerFreeState(t,theMagField.product());
     TrajectoryStateClosestToPoint tscp = tscpBuilder(fts, bs );
     
     float scale = scRef->energy()/tscp.momentum().mag();
   
     const math::XYZTLorentzVector momentum(tscp.momentum().x()*scale,
                        tscp.momentum().y()*scale,
                        tscp.momentum().z()*scale,
                        scRef->energy());
 
     
     Electron ele(t.charge(),momentum, t.vertex() );
     ele.setSuperCluster(scRef);
     edm::Ref<TrackCollection> myRef(tracksH,i);
     ele.setTrack(myRef);
     outEle.push_back(ele);
 
   }  // loop over tracks
 }

void EgammaHLTPixelMatchElectronAlgo::run	(	edm::Event &	e,
		reco::ElectronCollection &	outEle
	)

Definition at line 71 of file EgammaHLTPixelMatchElectronAlgo.cc.

References BSProducer_, edm::Event::getByLabel(), process(), and trackProducer_.

Referenced by EgammaHLTPixelMatchElectronProducers::produce().

                                                                                 {
 
   // get the input 
   edm::Handle<TrackCollection> tracksH;
   //  e.getByLabel(trackLabel_,trackInstanceName_,tracksH);
  e.getByLabel(trackProducer_,tracksH);
 
   //Get the Beam Spot position
   edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
   // iEvent.getByType(recoBeamSpotHandle);
   e.getByLabel(BSProducer_,recoBeamSpotHandle);
   // gets its position
   const BeamSpot::Point& BSPosition = recoBeamSpotHandle->position(); 
   Global3DPoint bs(BSPosition.x(),BSPosition.y(),0);
   process(tracksH,outEle,bs);
 
   return;
 }

void EgammaHLTPixelMatchElectronAlgo::setupES ( const edm::EventSetup & setup )

Definition at line 65 of file EgammaHLTPixelMatchElectronAlgo.cc.

References edm::EventSetup::get(), theGeomSearchTracker, and theMagField.

Referenced by EgammaHLTPixelMatchElectronProducers::produce().

                                                                      {
   //services
   es.get<TrackerRecoGeometryRecord>().get(theGeomSearchTracker);
   es.get<IdealMagneticFieldRecord>().get(theMagField);
 }