TauNNId.cc

Go to the documentation of this file.

#include "L1Trigger/Phase2L1ParticleFlow/interface/TauNNId.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"
#include "DataFormats/Math/interface/deltaPhi.h"
#include <cmath>

static constexpr unsigned int n_particles_max = 10;

TauNNId::TauNNId(const std::string &iInput,
                 const tensorflow::Session *session,
                 const std::string &iWeightFile,
                 int iNParticles)
    : session_(session) {
  NNvectorVar_.clear();
  edm::FileInPath fp(iWeightFile);
  fNParticles_ = iNParticles;

  fPt_ = std::make_unique<float[]>(fNParticles_);
  fEta_ = std::make_unique<float[]>(fNParticles_);
  fPhi_ = std::make_unique<float[]>(fNParticles_);
  fId_ = std::make_unique<float[]>(fNParticles_);
  fInput_ = iInput;
}

void TauNNId::setNNVectorVar() {
  NNvectorVar_.clear();
  for (int i0 = 0; i0 < fNParticles_; i0++) {
    NNvectorVar_.push_back(fPt_.get()[i0]);   //pT
    NNvectorVar_.push_back(fEta_.get()[i0]);  //dEta from jet axis
    NNvectorVar_.push_back(fPhi_.get()[i0]);  //dPhi from jet axis
    if (fPt_.get()[i0] == 0) {
      for (int i1 = 0; i1 < 5; i1++)
        NNvectorVar_.push_back(0);
      continue;
    }
    NNvectorVar_.push_back(fId_.get()[i0] == l1t::PFCandidate::Photon);         // Photon
    NNvectorVar_.push_back(fId_.get()[i0] == l1t::PFCandidate::Electron);       // Electron
    NNvectorVar_.push_back(fId_.get()[i0] == l1t::PFCandidate::Muon);           // Muon
    NNvectorVar_.push_back(fId_.get()[i0] == l1t::PFCandidate::NeutralHadron);  // Neutral Had
    NNvectorVar_.push_back(fId_.get()[i0] == l1t::PFCandidate::ChargedHadron);  // Charged Had
  }
}
float TauNNId::EvaluateNN() {
  tensorflow::Tensor input(tensorflow::DT_FLOAT,
                           {1, (unsigned int)NNvectorVar_.size()});  //was {1,35} but get size mismatch, CHECK
  for (unsigned int i = 0; i < NNvectorVar_.size(); i++) {
    input.matrix<float>()(0, i) = float(NNvectorVar_[i]);
  }
  std::vector<tensorflow::Tensor> outputs;
  tensorflow::run(session_, {{fInput_, input}}, {"dense_4/Sigmoid:0"}, &outputs);
  return outputs[0].matrix<float>()(0, 0);
}  //end EvaluateNN

float TauNNId::compute(const l1t::PFCandidate &iSeed, l1t::PFCandidateCollection &iParts) {
  for (int i0 = 0; i0 < fNParticles_; i0++) {
    fPt_.get()[i0] = 0;
    fEta_.get()[i0] = 0;
    fPhi_.get()[i0] = 0;
    fId_.get()[i0] = 0;
  }
  std::sort(iParts.begin(), iParts.end(), [](l1t::PFCandidate i, l1t::PFCandidate j) { return (i.pt() > j.pt()); });
  for (unsigned int i0 = 0; i0 < iParts.size(); i0++) {
    if (i0 > n_particles_max || i0 >= (unsigned int)fNParticles_)
      break;
    fPt_.get()[i0] = iParts[i0].pt();
    fEta_.get()[i0] = iSeed.eta() - iParts[i0].eta();
    fPhi_.get()[i0] = deltaPhi(iSeed.phi(), iParts[i0].phi());
    fId_.get()[i0] = iParts[i0].id();
  }
  setNNVectorVar();
  return EvaluateNN();
}