#include <SmirnovDeDxDiscriminator.h>

Inheritance diagram for SmirnovDeDxDiscriminator:

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

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

	SmirnovDeDxDiscriminator (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 SmirnovDeDxDiscriminator.h.

Constructor & Destructor Documentation

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

inline

Definition at line 10 of file SmirnovDeDxDiscriminator.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 SmirnovDeDxDiscriminator::beginRun	(	edm::Run const &	run,
		const edm::EventSetup &	iSetup
	)

inlineoverridevirtual

Reimplemented from BaseDeDxEstimator.

Definition at line 18 of file SmirnovDeDxDiscriminator.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> SmirnovDeDxDiscriminator::dedx ( const reco::DeDxHitCollection & Hits )

inlineoverridevirtual

Implements BaseDeDxEstimator.

Definition at line 23 of file SmirnovDeDxDiscriminator.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);
     std::sort(vect_probs.begin(), vect_probs.end(), std::less<float>());
     float TotalProb = 1.0 / (12 * size);
     for (size_t i = 1; i <= size; i++) {
       TotalProb += pow(vect_probs[i - 1] - ((2.0 * i - 1.0) / (2.0 * size)), 2);
     }
     TotalProb *= (3.0 / size);
     return std::make_pair(TotalProb, -1);
   }