#include <SmirnovDeDxDiscriminator.h>

Inheritance diagram for SmirnovDeDxDiscriminator:

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

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

	SmirnovDeDxDiscriminator (const edm::ParameterSet &iConfig)

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

Constructor & Destructor Documentation

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

inline

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

inlinevirtual

Reimplemented from BaseDeDxEstimator.

Definition at line 17 of file SmirnovDeDxDiscriminator.h.

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

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

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

inlinevirtual

Implements BaseDeDxEstimator.

Definition at line 21 of file SmirnovDeDxDiscriminator.h.

References reco_application_tbsim_DetSim-Digi_cfg::BinX, reco_application_tbsim_DetSim-Digi_cfg::BinY, i, meVperADCStrip, cmsHarvester::path, funct::pow(), mix_2012_Summer_50ns_PoissonOOTPU_cfi::prob, 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 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);
  }