#include <hugues/SeedToTrackProducer/plugins/SeedToTrackProducer.cc>

Inheritance diagram for SeedToTrackProducer:

Public Member Functions
	SeedToTrackProducer (const edm::ParameterSet &)

	~SeedToTrackProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
virtual void	beginJob () override

virtual void	endJob () override

virtual void	produce (edm::Event &, const edm::EventSetup &) override

virtual TrajectoryStateOnSurface	seedTransientState (const TrajectorySeed &)

Private Attributes
edm::EDGetTokenT< edm::View < TrajectorySeed > >	L2seedsTagS_

edm::EDGetTokenT < TrajectorySeedCollection >	L2seedsTagT_

edm::ESHandle< MagneticField >	theMGField

edm::ESHandle < GlobalTrackingGeometry >	theTrackingGeometry

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Description:

Definition at line 53 of file SeedToTrackProducer.h.

Constructor & Destructor Documentation

SeedToTrackProducer::SeedToTrackProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 47 of file SeedToTrackProducer.cc.

References edm::ParameterSet::getParameter(), L2seedsTagS_, and L2seedsTagT_.

 {
     
   L2seedsTagT_ = consumes<TrajectorySeedCollection>(iConfig.getParameter<edm::InputTag>("L2seedsCollection"));
   L2seedsTagS_ = consumes<edm::View<TrajectorySeed> >(iConfig.getParameter<edm::InputTag>("L2seedsCollection"));
 
     produces<reco::TrackCollection>();
     produces<reco::TrackExtraCollection>();
     produces<TrackingRecHitCollection>();
     
 }

SeedToTrackProducer::~SeedToTrackProducer ( )

Definition at line 64 of file SeedToTrackProducer.cc.

 {
  
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

void SeedToTrackProducer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 182 of file SeedToTrackProducer.cc.

183 {

184 }

void SeedToTrackProducer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 188 of file SeedToTrackProducer.cc.

188 {

189 }

void SeedToTrackProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 225 of file SeedToTrackProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                 {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

void SeedToTrackProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 79 of file SeedToTrackProducer.cc.

References reco::TrackBase::appendHitPattern(), clone(), funct::cos(), TrajectoryStateOnSurface::curvilinearError(), PV3DBase< T, PVType, FrameType >::eta(), reco::TrackBase::eta(), edm::EventSetup::get(), edm::Event::getByToken(), edm::Event::getRefBeforePut(), TrajectoryStateOnSurface::globalMomentum(), susybsm::HSCParticleType::globalMuon, i, customizeTrackingMonitorSeedNumber::idx, L2seedsTagS_, L2seedsTagT_, CurvilinearTrajectoryError::matrix(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), reco::TrackBase::phi(), reco::TrackBase::phiError(), reco::TrackBase::pt(), reco::TrackBase::ptError(), edm::Event::put(), HLT_25ns14e33_v3_cff::recHits, fileCollector::seed, seedTransientState(), TrackCollections2monitor_cff::selectedTracks, reco::Track::setExtra(), funct::sin(), theMGField, and theTrackingGeometry.

 {
     using namespace edm;
     using namespace std;
 
     std::auto_ptr<reco::TrackCollection> selectedTracks(new reco::TrackCollection);
     std::auto_ptr<reco::TrackExtraCollection> selectedTrackExtras( new reco::TrackExtraCollection() );
     std::auto_ptr<TrackingRecHitCollection> selectedTrackHits( new TrackingRecHitCollection() );
     
     reco::TrackRefProd rTracks = iEvent.getRefBeforePut<reco::TrackCollection>();
     reco::TrackExtraRefProd rTrackExtras = iEvent.getRefBeforePut<reco::TrackExtraCollection>();
     TrackingRecHitRefProd rHits = iEvent.getRefBeforePut<TrackingRecHitCollection>();
     
     edm::Ref<reco::TrackExtraCollection>::key_type hidx = 0;
     edm::Ref<reco::TrackExtraCollection>::key_type idx = 0;
     
     // magnetic fied and detector geometry
     iSetup.get<IdealMagneticFieldRecord>().get(theMGField);
     iSetup.get<GlobalTrackingGeometryRecord>().get(theTrackingGeometry);
 
     edm::ESHandle<TrackerTopology> httopo;
     iSetup.get<TrackerTopologyRcd>().get(httopo);
     const TrackerTopology& ttopo = *httopo;
 
     // now read the L2 seeds collection :
     edm::Handle<TrajectorySeedCollection> L2seedsCollection;
     iEvent.getByToken(L2seedsTagT_,L2seedsCollection);
     const std::vector<TrajectorySeed>* L2seeds = 0;
     if (L2seedsCollection.isValid()) L2seeds = L2seedsCollection.product();
     else  edm::LogError("SeedToTrackProducer") << "L2 seeds collection not found !! " << endl;
     
     edm::Handle<edm::View<TrajectorySeed> > seedHandle;
     iEvent.getByToken(L2seedsTagS_, seedHandle);
     
     
     int countRH = 0;
 
     // now  loop on the seeds :
     for (unsigned int i = 0; i < L2seeds->size() ; i++){
         
         
         //get the kinematic extrapolation from the seed
         TrajectoryStateOnSurface theTrajectory = seedTransientState(L2seeds->at(i));
         float seedEta =  theTrajectory.globalMomentum().eta();
         float seedPhi =  theTrajectory.globalMomentum().phi();
         float seedPt =  theTrajectory.globalMomentum().perp();
         CovarianceMatrix matrixSeedErr = theTrajectory.curvilinearError().matrix();
         edm::LogVerbatim("SeedToTrackProducer") <<  "seedPt=" << seedPt << " seedEta=" << seedEta << " seedPhi=" << seedPhi << endl;
         /*AlgebraicSymMatrix66 errors = theTrajectory.cartesianError().matrix();
         double partialPterror = errors(3,3)*pow(theTrajectory.globalMomentum().x(),2) + errors(4,4)*pow(theTrajectory.globalMomentum().y(),2);
     edm::LogVerbatim("SeedToTrackProducer") <<  "seedPtError=" << sqrt(partialPterror)/theTrajectory.globalMomentum().perp() << "seedPhiError=" << theTrajectory.curvilinearError().matrix()(2,2) << endl;*/ 
         //fill the track in a way that its pt, phi and eta will be the same as the seed
         math::XYZPoint initPoint(0,0,0);
         math::XYZVector initMom(seedPt*cos(seedPhi),seedPt*sin(seedPhi),seedPt*sinh(seedEta));
         reco::Track theTrack(1, 1, //dummy Chi2 and ndof
                              initPoint, initMom,
                              1, matrixSeedErr,
                              reco::TrackBase::TrackAlgorithm::globalMuon, reco::TrackBase::TrackQuality::tight);
 
         //fill the extra track with dummy information
         math::XYZPoint dummyFinalPoint(1,1,1);
         math::XYZVector dummyFinalMom(0,0,10);
         edm::RefToBase<TrajectorySeed> seed(seedHandle, i);
         CovarianceMatrix matrixExtra = ROOT::Math::SMatrixIdentity();
         reco::TrackExtra theTrackExtra(dummyFinalPoint, dummyFinalMom, true, initPoint, initMom, true,
                                   matrixSeedErr, 1,
                                   matrixExtra, 2,
                                   (L2seeds->at(i)).direction(), seed);
         theTrack.setExtra( reco::TrackExtraRef( rTrackExtras, idx++ ) );
         edm::LogVerbatim("SeedToTrackProducer") << "trackPt=" << theTrack.pt() << " trackEta=" << theTrack.eta() << " trackPhi=" << theTrack.phi() << endl;
         edm::LogVerbatim("SeedToTrackProducer") << "trackPtError=" << theTrack.ptError() << "trackPhiError=" << theTrack.phiError() << endl;
  
         //fill the seed segments in the track
         unsigned int nHitsAdded = 0;
         for(TrajectorySeed::recHitContainer::const_iterator itRecHits=(L2seeds->at(i)).recHits().first; itRecHits!=(L2seeds->at(i)).recHits().second; ++itRecHits, ++countRH) {
             TrackingRecHit* hit = (itRecHits)->clone();
             theTrack.appendHitPattern(*hit, ttopo);
             selectedTrackHits->push_back(hit);
             nHitsAdded++;
         }
         theTrackExtra.setHits( rHits, hidx, nHitsAdded );
         hidx += nHitsAdded;
         selectedTracks->push_back(theTrack);
         selectedTrackExtras->push_back(theTrackExtra);
 
     }
     iEvent.put(selectedTracks);
     iEvent.put(selectedTrackExtras);
     iEvent.put(selectedTrackHits);
     
 }

TrajectoryStateOnSurface SeedToTrackProducer::seedTransientState ( const TrajectorySeed & tmpSeed )

privatevirtual

Definition at line 171 of file SeedToTrackProducer.cc.

References PTrajectoryStateOnDet::detId(), TrajectorySeed::startingState(), theMGField, theTrackingGeometry, and trajectoryStateTransform::transientState().

Referenced by produce().

                                                                                              {
     
     PTrajectoryStateOnDet tmpTSOD = tmpSeed.startingState();
     DetId tmpDetId(tmpTSOD.detId());
     const GeomDet* tmpGeomDet = theTrackingGeometry->idToDet(tmpDetId);
     TrajectoryStateOnSurface tmpTSOS = trajectoryStateTransform::transientState(tmpTSOD, &(tmpGeomDet->surface()), &(*theMGField));
     return tmpTSOS;
 }