Merge pull request cms-sw#103 from MRD2F/CMSSW_10_2_X_tau_pog_DNNTauID

DNN tau IDs for CMSSW_10_2_X

RecoTauTag/RecoTau/BuildFile.xml

@@ -27,6 +27,7 @@
 <use   name="FastSimulation/BaseParticlePropagator"/>
 <use   name="FastSimulation/Particle"/>
 <use   name="roottmva"/>
+<use   name="PhysicsTools/TensorFlow"/>
 <export>
   <lib   name="1"/>
 </export>

RecoTauTag/RecoTau/interface/DeepTauBase.h

@@ -0,0 +1,80 @@
+#ifndef RecoTauTag_RecoTau_DeepTauBase_h
+#define RecoTauTag_RecoTau_DeepTauBase_h
+
+/*
+ * \class DeepTauBase
+ *
+ * Definition of the base class for tau identification using Deep NN.
+ *
+ * \author Konstantin Androsov, INFN Pisa
+ * \author Maria Rosaria Di Domenico, University of Siena & INFN Pisa
+ */
+
+#include <Math/VectorUtil.h>
+#include "FWCore/Framework/interface/stream/EDProducer.h"
+#include "FWCore/ParameterSet/interface/ParameterSet.h"
+#include "PhysicsTools/TensorFlow/interface/TensorFlow.h"
+#include "DataFormats/PatCandidates/interface/Electron.h"
+#include "DataFormats/PatCandidates/interface/Muon.h"
+#include "DataFormats/PatCandidates/interface/Tau.h"
+#include "DataFormats/PatCandidates/interface/PATTauDiscriminator.h"
+#include "CommonTools/Utils/interface/StringObjectFunction.h"
+#include "RecoTauTag/RecoTau/interface/PFRecoTauClusterVariables.h"
+#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
+#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
+
+namespace deep_tau {
+
+class DeepTauBase : public edm::stream::EDProducer<> {
+public:
+    using TauType = pat::Tau;
+    using TauDiscriminator = pat::PATTauDiscriminator;
+    using TauCollection = std::vector<TauType>;
+    using TauRef = edm::Ref<TauCollection>;
+    using TauRefProd = edm::RefProd<TauCollection>;
+    using ElectronCollection = pat::ElectronCollection;
+    using MuonCollection = pat::MuonCollection;
+    using LorentzVectorXYZ = ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double>>;
+    using GraphPtr = std::shared_ptr<tensorflow::GraphDef>;
+    using Cutter = StringObjectFunction<TauType>;
+    using CutterPtr = std::unique_ptr<Cutter>;
+    using WPMap = std::map<std::string, CutterPtr>;
+
+
+    struct Output {
+        using ResultMap = std::map<std::string, std::unique_ptr<TauDiscriminator>>;
+        std::vector<size_t> num, den;
+
+        Output(const std::vector<size_t>& _num, const std::vector<size_t>& _den) : num(_num), den(_den) {}
+
+        ResultMap get_value(const edm::Handle<TauCollection>& taus, const tensorflow::Tensor& pred,
+                            const WPMap& working_points) const;
+    };
+
+    using OutputCollection = std::map<std::string, Output>;
+
+
+    DeepTauBase(const edm::ParameterSet& cfg, const OutputCollection& outputs);
+    virtual ~DeepTauBase();
+
+    virtual void produce(edm::Event& event, const edm::EventSetup& es) override;
+
+private:
+    virtual tensorflow::Tensor GetPredictions(edm::Event& event, const edm::EventSetup& es,
+                                              edm::Handle<TauCollection> taus) = 0;
+    virtual void CreateOutputs(edm::Event& event, const tensorflow::Tensor& pred, edm::Handle<TauCollection> taus);
+
+protected:
+    edm::EDGetTokenT<TauCollection> taus_token;
+    std::string graphName;
+    GraphPtr graph;
+    tensorflow::Session* session;
+    std::map<std::string, WPMap> working_points;
+    OutputCollection outputs;
+};
+
+} // namespace deep_tau
+
+
+
+#endif

RecoTauTag/RecoTau/plugins/BuildFile.xml

@@ -36,5 +36,6 @@
   <use   name="MagneticField/Engine"/>
   <use   name="MagneticField/Records"/>
   <use   name="FastSimulation/BaseParticlePropagator"/>
+  <use   name="PhysicsTools/TensorFlow"/>
   <use   name="root"/>
 </library>