#include <ProductDeDxDiscriminator.h>

Inheritance diagram for ProductDeDxDiscriminator:

Public Member Functions
void	beginRun (edm::Run const &, const edm::EventSetup &iSetup) override

std::pair< float, float >	dedx (const reco::DeDxHitCollection &Hits) override

	ProductDeDxDiscriminator (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iCollector)

Public Member Functions inherited from BaseDeDxEstimator
virtual	~BaseDeDxEstimator ()

Private Attributes
float	meVperADCStrip

TH3F *	Prob_ChargePath

std::string	ProbabilityMode

DeDxTools::ESGetTokenH3DDVariant	token_

Detailed Description

Definition at line 8 of file ProductDeDxDiscriminator.h.

Constructor & Destructor Documentation

ProductDeDxDiscriminator::ProductDeDxDiscriminator	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &	iCollector
	)

inline

Definition at line 10 of file ProductDeDxDiscriminator.h.

References edm::ParameterSet::getParameter(), meVperADCStrip, Prob_ChargePath, ProbabilityMode, and AlCaHLTBitMon_QueryRunRegistry::string.

       : token_(DeDxTools::esConsumes(iConfig.getParameter<std::string>("Reccord"), iCollector)) {
     meVperADCStrip =
         iConfig.getParameter<double>("MeVperADCStrip");  //currently needed until the map on the database are redone
     ProbabilityMode = iConfig.getParameter<std::string>("ProbabilityMode");
     Prob_ChargePath = nullptr;
   }

Member Function Documentation

void ProductDeDxDiscriminator::beginRun	(	edm::Run const &	,
		const edm::EventSetup &	iSetup
	)

inlineoverridevirtual

Reimplemented from BaseDeDxEstimator.

Definition at line 18 of file ProductDeDxDiscriminator.h.

References DeDxTools::buildDiscrimMap(), DeDxTools::getHistogramD3D(), Prob_ChargePath, ProbabilityMode, and token_.

                                                                      {
     auto const& histD3D = DeDxTools::getHistogramD3D(iSetup, token_);
     DeDxTools::buildDiscrimMap(histD3D, ProbabilityMode, Prob_ChargePath);
   }

std::pair<float, float> ProductDeDxDiscriminator::dedx ( const reco::DeDxHitCollection & Hits )

inlineoverridevirtual

Implements BaseDeDxEstimator.

Definition at line 23 of file ProductDeDxDiscriminator.h.

References reco_application_tbsim_DetSim-Digi_cfg::BinX, reco_application_tbsim_DetSim-Digi_cfg::BinY, RecoTauCleanerPlugins::charge, mps_fire::i, meVperADCStrip, fed_dqm_sourceclient-live_cfg::path, funct::pow(), Prob_ChargePath, and findQualityFiles::size.

                                                                          {
     std::vector<float> vect_probs;
     for (size_t i = 0; i < Hits.size(); i++) {
       float path = Hits[i].pathLength() * 10.0;  //x10 in order to be compatible with the map content
       float charge =
           Hits[i].charge() /
           (10.0 *
            meVperADCStrip);  // 10/meVperADCStrip in order to be compatible with the map content in ADC/mm instead of MeV/cm
 
       int BinX = Prob_ChargePath->GetXaxis()->FindBin(Hits[i].momentum());
       int BinY = Prob_ChargePath->GetYaxis()->FindBin(path);
       int BinZ = Prob_ChargePath->GetZaxis()->FindBin(charge);
       float prob = Prob_ChargePath->GetBinContent(BinX, BinY, BinZ);
       if (prob >= 0)
         vect_probs.push_back(prob);
     }
 
     size_t size = vect_probs.size();
     if (size <= 0)
       return std::make_pair(-1, -1);
     float TotalProb = 1;
     for (size_t i = 0; i < size; i++) {
       if (vect_probs[i] <= 0.0001) {
         TotalProb *= pow(0.0001, 1.0 / size);
       } else {
         TotalProb *= pow(vect_probs[i], 1.0 / size);
       }
     }
     return std::make_pair(TotalProb, -1);
   }