1 #ifndef RecoTrackerDeDx_ProductDeDxDiscriminator_h
2 #define RecoTrackerDeDx_ProductDeDxDiscriminator_h
22 std::vector<float> vect_probs;
23 for(
size_t i = 0;
i< Hits.size();
i ++){
24 float path = Hits[
i].pathLength() * 10.0;
31 if(prob>=0)vect_probs.push_back(prob);
34 size_t size = vect_probs.size();
35 if(size<=0)
return std::make_pair( -1 , -1);
38 if(vect_probs[
i]<=0.0001){TotalProb *=
pow(0.0001 , 1.0/size);}
39 else {TotalProb *=
pow(vect_probs[
i], 1.0/size);}
41 return std::make_pair( TotalProb , -1);
T getParameter(std::string const &) const
virtual std::pair< float, float > dedx(const reco::DeDxHitCollection &Hits)
std::string ProbabilityMode
std::vector< DeDxHit > DeDxHitCollection
ProductDeDxDiscriminator(const edm::ParameterSet &iConfig)
virtual void beginRun(edm::Run const &run, const edm::EventSetup &iSetup)
tuple size
Write out results.
Power< A, B >::type pow(const A &a, const B &b)