10 #include "TLorentzVector.h"
21 chargedPionMass = 0.139;
22 booleanOutput = iConfig.
getParameter<
bool>(
"BooleanOutput");
44 double dE = DeltaE(tau);
45 if(booleanOutput)
return ( dE > deltaEmin && dE < deltaEmax ? 1. : 0. );
52 for(
size_t i = 0;
i < signalTracks.
size(); ++
i){
54 p4.SetXYZM(signalTracks[
i]->px(),
55 signalTracks[
i]->py(),
56 signalTracks[
i]->pz(),
61 double hadrTauP = tau->momentum().r() * (1.0 - tau->emFraction());
62 if (tau->emFraction() >= 1.0) {
65 return tracksE / hadrTauP - 1.0;
T getParameter(std::string const &) const
PFRecoTauDiscriminationByDeltaE(const ParameterSet &iConfig)
double DeltaE(const PFTauRef &)
DEFINE_FWK_MODULE(CosmicTrackingParticleSelector)
double discriminate(const reco::PFTauRef &)
size_type size() const
Size of the RefVector.
void beginEvent(const edm::Event &, const edm::EventSetup &)
~PFRecoTauDiscriminationByDeltaE()