#include <BTagLikeDeDxDiscriminator.h>

Inheritance diagram for BTagLikeDeDxDiscriminator:

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

	BTagLikeDeDxDiscriminator (const edm::ParameterSet &iConfig)

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

Public Member Functions inherited from BaseDeDxEstimator
virtual	~BaseDeDxEstimator ()

Private Attributes
float	meVperADCStrip

TH3F *	Prob_ChargePath

std::string	ProbabilityMode

std::string	Reccord

Detailed Description

Definition at line 7 of file BTagLikeDeDxDiscriminator.h.

Constructor & Destructor Documentation

BTagLikeDeDxDiscriminator::BTagLikeDeDxDiscriminator ( const edm::ParameterSet & iConfig )

inline

Definition at line 10 of file BTagLikeDeDxDiscriminator.h.

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

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

Member Function Documentation

virtual void BTagLikeDeDxDiscriminator::beginRun	(	edm::Run const &	run,
		const edm::EventSetup &	iSetup
	)

inlinevirtual

Reimplemented from BaseDeDxEstimator.

Definition at line 17 of file BTagLikeDeDxDiscriminator.h.

References DeDxTools::buildDiscrimMap(), Prob_ChargePath, ProbabilityMode, and Reccord.

                                                                       {
     DeDxTools::buildDiscrimMap(run, iSetup, Reccord,  ProbabilityMode, Prob_ChargePath);
  }

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

inlinevirtual

Implements BaseDeDxEstimator.

Definition at line 21 of file BTagLikeDeDxDiscriminator.h.

References reco_application_tbsim_DetSim-Digi_cfg::BinX, reco_application_tbsim_DetSim-Digi_cfg::BinY, RecoTauCleanerPlugins::charge, create_public_lumi_plots::exp, f, i, prof2calltree::l, cmsBatch::log, bookConverter::max, meVperADCStrip, min(), fed_dqm_sourceclient-live_cfg::path, Prob_ChargePath, findQualityFiles::size, and python.multivaluedict::sort().

                                                                      {
     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 SumJet = 0.;
     for(size_t i=0;i<size;i++){if(vect_probs[i]<=0.0001)vect_probs[i] = 0.0001;   SumJet+= log(vect_probs[i]);  }
 
     float Loginvlog=log(-SumJet);
     float Prob =1;
     float lfact=1;
 
     for(size_t l=1; l!=size; l++){
          lfact*=l;
          Prob+=exp(l*Loginvlog-log(lfact));
     }
 
     float LogProb=log(Prob);
     float ProbJet=std::min((float)exp(std::max(LogProb+SumJet,-30.0f)),1.0f);
     float TotalProb = -log10(ProbJet)/4;
     TotalProb = 1-TotalProb;
 
     return std::make_pair( TotalProb , -1);
  }