ProductDeDxDiscriminator Class Reference

#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::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>("Record"), 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

beginRun()

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);
  }

dedx()

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

inlineoverridevirtual

Implements BaseDeDxEstimator.

Definition at line 23 of file ProductDeDxDiscriminator.h.

References VtxSmearedBeamProfile_cfi::BinX, VtxSmearedBeamProfile_cfi::BinY, ALCARECOTkAlJpsiMuMu_cff::charge, mps_fire::i, meVperADCStrip, castor_dqm_sourceclient_file_cfg::path, conifer::pow(), TtFullHadEvtBuilder_cfi::prob, 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);
  }

Member Data Documentation

meVperADCStrip

float ProductDeDxDiscriminator::meVperADCStrip

private

Definition at line 55 of file ProductDeDxDiscriminator.h.

Referenced by dedx(), and ProductDeDxDiscriminator().

Prob_ChargePath

TH3F* ProductDeDxDiscriminator::Prob_ChargePath

private

Definition at line 58 of file ProductDeDxDiscriminator.h.

Referenced by beginRun(), dedx(), and ProductDeDxDiscriminator().

ProbabilityMode

std::string ProductDeDxDiscriminator::ProbabilityMode

private

Definition at line 57 of file ProductDeDxDiscriminator.h.

Referenced by beginRun(), and ProductDeDxDiscriminator().

token_

deDxTools::ESGetTokenH3DDVariant ProductDeDxDiscriminator::token_

private

Definition at line 56 of file ProductDeDxDiscriminator.h.

Referenced by beginRun().