#include <LowPtGsfElectronSeedValueMapsProducer.h>

Inheritance diagram for LowPtGsfElectronSeedValueMapsProducer:

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

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

	~LowPtGsfElectronSeedValueMapsProducer () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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

Private Attributes
const edm::EDGetTokenT< reco::GsfTrackCollection >	gsfTracks_

const std::vector< std::string >	names_

const edm::EDGetTokenT< edm::ValueMap< reco::PreIdRef > >	preIdsValueMap_

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

Detailed Description

Definition at line 13 of file LowPtGsfElectronSeedValueMapsProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 15 of file LowPtGsfElectronSeedValueMapsProducer.cc.

References dataset::name, and names_.

                                                                                                           :
   gsfTracks_(consumes<reco::GsfTrackCollection>(conf.getParameter<edm::InputTag>("gsfTracks"))),
   preIdsValueMap_(consumes< edm::ValueMap<reco::PreIdRef> >(conf.getParameter<edm::InputTag>("preIdsValueMap"))),
   names_(conf.getParameter< std::vector<std::string> >("ModelNames"))
 {
   for ( const auto& name : names_ ) { produces< edm::ValueMap<float> >(name); }
 }

LowPtGsfElectronSeedValueMapsProducer::~LowPtGsfElectronSeedValueMapsProducer ( )

override

Definition at line 25 of file LowPtGsfElectronSeedValueMapsProducer.cc.

25 {}

Member Function Documentation

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

static

Definition at line 76 of file LowPtGsfElectronSeedValueMapsProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and DEFINE_FWK_MODULE.

 {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("gsfTracks",edm::InputTag("lowPtGsfEleGsfTracks"));
   desc.add<edm::InputTag>("preIdsValueMap",edm::InputTag("lowPtGsfElectronSeeds"));
   desc.add< std::vector<std::string> >("ModelNames",std::vector<std::string>());
   descriptions.add("defaultLowPtGsfElectronSeedValueMaps",desc);
 }

void LowPtGsfElectronSeedValueMapsProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	setup
	)

override

Definition at line 29 of file LowPtGsfElectronSeedValueMapsProducer.cc.

References begin, end, edm::helper::Filler< Map >::fill(), objects.autophobj::filler, ntupleEnum::gsf, gedGsfElectronCores_cfi::gsfTracks, gsfTracks_, edm::helper::Filler< Map >::insert(), edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), eostools::move(), names_, convertSQLitetoXML_cfg::output, preIdsValueMap_, and SurveyInfoScenario_cff::seed.

                                                                                                  {
 
   // Retrieve GsfTracks from Event
   edm::Handle<reco::GsfTrackCollection> gsfTracks;
   event.getByToken(gsfTracks_,gsfTracks);
   if ( !gsfTracks.isValid() ) { edm::LogError("Problem with gsfTracks handle"); }
 
   // Retrieve PreIds from Event
   edm::Handle< edm::ValueMap<reco::PreIdRef> > preIdsValueMap;
   event.getByToken(preIdsValueMap_,preIdsValueMap);
   if ( !preIdsValueMap.isValid() ) { edm::LogError("Problem with preIdsValueMap handle"); }
   
   // Iterate through GsfTracks, extract BDT output, and store result in ValueMap for each model
   std::vector< std::vector<float> > output;
   for ( unsigned int iname = 0; iname < names_.size(); ++iname ) { 
     output.push_back( std::vector<float>(gsfTracks->size(),-999.) );
   }
   for ( unsigned int igsf = 0; igsf < gsfTracks->size(); igsf++ ) {
     reco::GsfTrackRef gsf(gsfTracks,igsf);
     if ( gsf.isNonnull() && 
      gsf->extra().isNonnull() && 
      gsf->extra()->seedRef().isNonnull() ) {
       reco::ElectronSeedRef seed = gsf->extra()->seedRef().castTo<reco::ElectronSeedRef>();
       if ( seed.isNonnull() && seed->ctfTrack().isNonnull() ) {
     const reco::PreIdRef preid = (*preIdsValueMap)[seed->ctfTrack()];
     if ( preid.isNonnull() ) {
       for ( unsigned int iname = 0; iname < names_.size(); ++iname ) {
         output[iname][igsf] = preid->mva(iname);
       }
     }
       }
     }
   }
   
   // Create and put ValueMap in Event
   for ( unsigned int iname = 0; iname < names_.size(); ++iname ) {
     auto ptr = std::make_unique< edm::ValueMap<float> >( edm::ValueMap<float>() );
     edm::ValueMap<float>::Filler filler(*ptr);
     filler.insert(gsfTracks, output[iname].begin(), output[iname].end());
     filler.fill();
     event.put(std::move(ptr),names_[iname]);
   }
 
 }