#include <ConversionTrackProducer.h>

Inheritance diagram for ConversionTrackProducer:

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

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

virtual	~ConversionTrackProducer ()

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
	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
	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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Types
typedef edm::AssociationMap < edm::OneToOne< std::vector < Trajectory > , reco::GsfTrackCollection, unsigned short > >	TrajGsfTrackAssociationCollection

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotInputTag

edm::ParameterSet	conf_

IdealHelixParameters	ConvTrackPreSelector

bool	filterOnConvTrackHyp

edm::EDGetTokenT< edm::View < reco::Track > >	genericTracks

edm::EDGetTokenT < TrajGsfTrackAssociationCollection >	gsfTrajectories

edm::EDGetTokenT < TrajTrackAssociationCollection >	kfTrajectories

double	minConvRadius

std::auto_ptr < reco::ConversionTrackCollection >	outputTrks

bool	setArbitratedEcalSeeded

bool	setArbitratedMerged

bool	setArbitratedMergedEcalGeneral

bool	setTrackerOnly

std::string	trackProducer

bool	useTrajectory

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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

Definition at line 43 of file ConversionTrackProducer.h.

Member Typedef Documentation

typedef edm::AssociationMap<edm::OneToOne<std::vector<Trajectory>, reco::GsfTrackCollection,unsigned short> > ConversionTrackProducer::TrajGsfTrackAssociationCollection

private

Definition at line 48 of file ConversionTrackProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 36 of file ConversionTrackProducer.cc.

References genericTracks, gsfTrajectories, kfTrajectories, and trackProducer.

                                                                             : 
   conf_(conf),
   trackProducer ( conf.getParameter<std::string>("TrackProducer") ),
   useTrajectory ( conf.getParameter<bool>("useTrajectory") ),
   setTrackerOnly ( conf.getParameter<bool>("setTrackerOnly") ),
   setArbitratedEcalSeeded ( conf.getParameter<bool>("setArbitratedEcalSeeded") ),    
   setArbitratedMerged ( conf.getParameter<bool>("setArbitratedMerged") ),
   setArbitratedMergedEcalGeneral ( conf.getParameter<bool>("setArbitratedMergedEcalGeneral") ),
   beamSpotInputTag (  consumes<reco::BeamSpot>(conf.getParameter<edm::InputTag>("beamSpotInputTag")) ),
   filterOnConvTrackHyp( conf.getParameter<bool>("filterOnConvTrackHyp") ),
   minConvRadius( conf.getParameter<double>("minConvRadius") )
 {
   edm::InputTag thetp(trackProducer);
   genericTracks =
     consumes<edm::View<reco::Track> >(thetp);
   kfTrajectories = 
     consumes<TrajTrackAssociationCollection>(thetp);
   gsfTrajectories = 
     consumes<TrajGsfTrackAssociationCollection>(thetp);
   produces<reco::ConversionTrackCollection>();
   
 }

ConversionTrackProducer::~ConversionTrackProducer ( )

virtual

Definition at line 61 of file ConversionTrackProducer.cc.

61 { }

Member Function Documentation

void ConversionTrackProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

virtual

Implements edm::stream::EDProducerBase.

Definition at line 64 of file ConversionTrackProducer.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, beamSpotInputTag, edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::GsfTrackCollection, unsigned short > >::const_iterator, edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, ConvTrackPreSelector, filterOnConvTrackHyp, genericTracks, edm::EventSetup::get(), edm::Event::getByToken(), gsfTrajectories, IdealHelixParameters::isTangentPointDistanceLessThan(), kfTrajectories, minConvRadius, outputTrks, edm::ESHandle< class >::product(), edm::Event::put(), setArbitratedEcalSeeded, setArbitratedMerged, setArbitratedMergedEcalGeneral, reco::ConversionTrack::setIsArbitratedEcalSeeded(), reco::ConversionTrack::setIsArbitratedMerged(), reco::ConversionTrack::setIsArbitratedMergedEcalGeneral(), reco::ConversionTrack::setIsTrackerOnly(), IdealHelixParameters::setMagnField(), setTrackerOnly, reco::ConversionTrack::setTrajRef(), and useTrajectory.

   {
     //get input collection (through edm::View)
     edm::Handle<edm::View<reco::Track> > hTrks;
     e.getByToken(genericTracks, hTrks);
 
     //get association maps between trajectories and tracks and build temporary maps
     edm::Handle< TrajTrackAssociationCollection > hTTAss;
     edm::Handle< TrajGsfTrackAssociationCollection > hTTAssGsf;    
                                           
     std::map<reco::TrackRef,edm::Ref<std::vector<Trajectory> > > tracktrajmap;
     std::map<reco::GsfTrackRef,edm::Ref<std::vector<Trajectory> > > gsftracktrajmap;
                                           
     if (useTrajectory) {
       if (hTrks->size()>0) {
         if (dynamic_cast<const reco::GsfTrack*>(&hTrks->at(0))) {
           //fill map for gsf tracks
           e.getByToken(gsfTrajectories, hTTAssGsf);     
           for ( TrajGsfTrackAssociationCollection::const_iterator iPair = 
           hTTAssGsf->begin();
         iPair != hTTAssGsf->end(); ++iPair) {
         
             gsftracktrajmap[iPair->val] = iPair->key;
 
           }
                 
         }
         else {
           //fill map for standard tracks
           e.getByToken(kfTrajectories, hTTAss);
           for ( TrajTrackAssociationCollection::const_iterator iPair = hTTAss->begin();
             iPair != hTTAss->end();
             ++iPair) {
         
             tracktrajmap[iPair->val] = iPair->key;
 
           }
         }
       }
     }
 
     // Step B: create empty output collection
     outputTrks = std::auto_ptr<reco::ConversionTrackCollection>(new reco::ConversionTrackCollection);    
 
     //--------------------------------------------------
     //Added by D. Giordano
     // 2011/08/05
     // Reduction of the track sample based on geometric hypothesis for conversion tracks
  
     edm::Handle<reco::BeamSpot> beamSpotHandle;
     e.getByToken(beamSpotInputTag,beamSpotHandle);
    
     edm::ESHandle<MagneticField> magFieldHandle;
     es.get<IdealMagneticFieldRecord>().get( magFieldHandle );
 
 
     if(filterOnConvTrackHyp && !beamSpotHandle.isValid()) {
       edm::LogError("Invalid Collection") 
     << "invalid collection for the BeamSpot";
       throw;
     }
 
     ConvTrackPreSelector.setMagnField(magFieldHandle.product());
 
     //----------------------------------------------------------
    
  
     // Simple conversion of tracks to conversion tracks, setting appropriate flags from configuration
     for (edm::RefToBaseVector<reco::Track>::const_iterator it = hTrks->refVector().begin(); it != hTrks->refVector().end(); ++it) {
  
       //--------------------------------------------------
       //Added by D. Giordano
       // 2011/08/05
       // Reduction of the track sample based on geometric hypothesis for conversion tracks
       
       math::XYZVector beamSpot=  math::XYZVector(beamSpotHandle->position());
 
       if( filterOnConvTrackHyp && ConvTrackPreSelector.isTangentPointDistanceLessThan( minConvRadius, it->get(), beamSpot )  )
     continue;
       //--------------------------------------------------
 
       reco::ConversionTrack convTrack(*it);
       convTrack.setIsTrackerOnly(setTrackerOnly);
       convTrack.setIsArbitratedEcalSeeded(setArbitratedEcalSeeded);
       convTrack.setIsArbitratedMerged(setArbitratedMerged);
       convTrack.setIsArbitratedMergedEcalGeneral(setArbitratedMergedEcalGeneral);
             
       //fill trajectory association if configured, using correct map depending on track type
       if (useTrajectory) {
         if (gsftracktrajmap.size()) {
           convTrack.setTrajRef(gsftracktrajmap.find(it->castTo<reco::GsfTrackRef>())->second);
         }
         else {
           convTrack.setTrajRef(tracktrajmap.find(it->castTo<reco::TrackRef>())->second);
         }
       }
       
       outputTrks->push_back(convTrack);
     }
     
     e.put(outputTrks);
     return;
 
   }//end produce

Member Data Documentation

edm::EDGetTokenT<reco::BeamSpot> ConversionTrackProducer::beamSpotInputTag

private

Definition at line 77 of file ConversionTrackProducer.h.

Referenced by produce().

edm::ParameterSet ConversionTrackProducer::conf_

private

Definition at line 60 of file ConversionTrackProducer.h.

IdealHelixParameters ConversionTrackProducer::ConvTrackPreSelector

private

Definition at line 80 of file ConversionTrackProducer.h.

Referenced by produce().

bool ConversionTrackProducer::filterOnConvTrackHyp

private

Definition at line 78 of file ConversionTrackProducer.h.

Referenced by produce().

edm::EDGetTokenT<edm::View<reco::Track> > ConversionTrackProducer::genericTracks

private

Definition at line 63 of file ConversionTrackProducer.h.

Referenced by ConversionTrackProducer(), and produce().

edm::EDGetTokenT<TrajGsfTrackAssociationCollection> ConversionTrackProducer::gsfTrajectories

private

Definition at line 65 of file ConversionTrackProducer.h.

Referenced by ConversionTrackProducer(), and produce().

edm::EDGetTokenT<TrajTrackAssociationCollection> ConversionTrackProducer::kfTrajectories

private

Definition at line 64 of file ConversionTrackProducer.h.

Referenced by ConversionTrackProducer(), and produce().

double ConversionTrackProducer::minConvRadius

private

Definition at line 79 of file ConversionTrackProducer.h.

Referenced by produce().

std::auto_ptr<reco::ConversionTrackCollection> ConversionTrackProducer::outputTrks

private

Definition at line 83 of file ConversionTrackProducer.h.

Referenced by produce().

bool ConversionTrackProducer::setArbitratedEcalSeeded

private

Definition at line 68 of file ConversionTrackProducer.h.

Referenced by produce().

bool ConversionTrackProducer::setArbitratedMerged

private

Definition at line 69 of file ConversionTrackProducer.h.

Referenced by produce().

bool ConversionTrackProducer::setArbitratedMergedEcalGeneral

private

Definition at line 70 of file ConversionTrackProducer.h.

Referenced by produce().

bool ConversionTrackProducer::setTrackerOnly

private

Definition at line 67 of file ConversionTrackProducer.h.

Referenced by produce().

std::string ConversionTrackProducer::trackProducer

private

Definition at line 62 of file ConversionTrackProducer.h.

Referenced by ConversionTrackProducer().

bool ConversionTrackProducer::useTrajectory

private

Definition at line 66 of file ConversionTrackProducer.h.

Referenced by produce().

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation