Inheritance diagram for l1t::HGC3DClusterTMVASelector:

Classes
class	Var

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

	~HGC3DClusterTMVASelector () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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

Private Attributes
std::string	method_

StringCutObjectSelector< l1t::HGCalMulticluster >	preselection_

std::unique_ptr< TMVA::Reader >	reader_

edm::EDGetTokenT< l1t::HGCalMulticlusterBxCollection >	src_

std::vector< Var >	variables_

std::string	weightsFile_

StringObjectFunction< l1t::HGCalMulticluster >	wp_

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 16 of file HGC3DClusterTMVASelector.cc.

Constructor & Destructor Documentation

l1t::HGC3DClusterTMVASelector::HGC3DClusterTMVASelector ( const edm::ParameterSet & iConfig )

explicit

Definition at line 46 of file HGC3DClusterTMVASelector.cc.

References edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), reco::details::loadTMVAWeights(), method_, reader_, AlCaHLTBitMon_QueryRunRegistry::string, JetChargeProducer_cfi::var, variables_, and weightsFile_.

     : src_(consumes<l1t::HGCalMulticlusterBxCollection>(iConfig.getParameter<edm::InputTag>("src"))),
       preselection_(iConfig.getParameter<std::string>("preselection")),
       method_(iConfig.getParameter<std::string>("method")),
       weightsFile_(iConfig.getParameter<std::string>("weightsFile")),
       reader_(new TMVA::Reader()),
       wp_(iConfig.getParameter<std::string>("wp")) {
   // first create all the variables
   for (const auto &psvar : iConfig.getParameter<std::vector<edm::ParameterSet>>("variables")) {
     variables_.emplace_back(psvar.getParameter<std::string>("name"), psvar.getParameter<std::string>("value"));
   }
   // then declare them
   for (auto &var : variables_)
     var.declare(*reader_);
   // then read the weights
   if (weightsFile_[0] != '/' && weightsFile_[0] != '.') {
     weightsFile_ = edm::FileInPath(weightsFile_).fullPath();
   }
   reco::details::loadTMVAWeights(&*reader_, method_, weightsFile_);
   // finally, declare outputs
   produces<l1t::HGCalMulticlusterBxCollection>();
   produces<l1t::HGCalMulticlusterBxCollection>("fail");
 }

l1t::HGC3DClusterTMVASelector::~HGC3DClusterTMVASelector ( )

inlineoverride

Definition at line 19 of file HGC3DClusterTMVASelector.cc.

19 {}

Member Function Documentation

void l1t::HGC3DClusterTMVASelector::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &
	)

overrideprivate

Definition at line 70 of file HGC3DClusterTMVASelector.cc.

References EnergyCorrector::c, DEFINE_FWK_MODULE, cmsPerfPublish::fail(), edm::Event::getByToken(), method_, eostools::move(), MillePedeFileConverter_cfg::out, preselection_, edm::Event::put(), reader_, src_, JetChargeProducer_cfi::var, variables_, and wp_.

                                                                                  {
   std::unique_ptr<l1t::HGCalMulticlusterBxCollection> out = std::make_unique<l1t::HGCalMulticlusterBxCollection>();
   std::unique_ptr<l1t::HGCalMulticlusterBxCollection> fail = std::make_unique<l1t::HGCalMulticlusterBxCollection>();
 
   edm::Handle<l1t::HGCalMulticlusterBxCollection> multiclusters;
   iEvent.getByToken(src_, multiclusters);
 
   for (int bx = multiclusters->getFirstBX(); bx <= multiclusters->getLastBX(); ++bx) {
     for (auto it = multiclusters->begin(bx), ed = multiclusters->end(bx); it != ed; ++it) {
       const auto &c = *it;
       if (preselection_(c)) {
         for (auto &var : variables_)
           var.fill(c);
         float mvaOut = reader_->EvaluateMVA(method_);
         if (mvaOut > wp_(c)) {
           out->push_back(bx, c);
         } else {
           fail->push_back(bx, c);
         }
       }
     }
   }
 
   iEvent.put(std::move(out));
   iEvent.put(std::move(fail), "fail");
 }