PileupJetIdAlgo.h

Go to the documentation of this file.

//--------------------------------------------------------------------------------------------------
//
// PileupJetIdAlgo
//
// Author: P. Musella, P. Harris
//--------------------------------------------------------------------------------------------------

#ifndef RecoJets_JetProducers_plugins_PileupJetIdAlgo_h
#define RecoJets_JetProducers_plugins_PileupJetIdAlgo_h

#include "DataFormats/JetReco/interface/Jet.h"
#include "DataFormats/JetReco/interface/PFJetCollection.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/PatCandidates/interface/Jet.h"
#include "DataFormats/PatCandidates/interface/PackedCandidate.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "TMVA/Tools.h"
#include "TMVA/Reader.h"
#include "TMVA/Tools.h"
#include "TMVA/Reader.h"

#include "DataFormats/JetReco/interface/PileupJetIdentifier.h"


// ----------------------------------------------------------------------------------------------------
class PileupJetIdAlgo {
public:
    enum version_t { USER=-1, PHILv0=0 };
    
    PileupJetIdAlgo(int version=PHILv0, const std::string & tmvaWeight="", const std::string & tmvaMethod="", 
            Float_t impactParTkThreshod_=1., const std::vector<std::string> & tmvaVariables= std::vector<std::string>(), bool runMvas=true);
    PileupJetIdAlgo(const edm::ParameterSet & ps, bool runMvas); 
    ~PileupJetIdAlgo(); 
    
    PileupJetIdentifier computeIdVariables(const reco::Jet * jet, 
                           float jec, const reco::Vertex *, const reco::VertexCollection &, double rho);

    void set(const PileupJetIdentifier &);
    PileupJetIdentifier computeMva();
    const std::string method() const { return tmvaMethod_; }
    
    std::string dumpVariables() const;

    typedef std::map<std::string,std::pair<float *,float> > variables_list_t;

    std::pair<int,int> getJetIdKey(float jetPt, float jetEta);
    int computeCutIDflag(float betaStarClassic,float dR2Mean,float nvtx, float jetPt, float jetEta);
    int computeIDflag   (float mva, float jetPt, float jetEta);
    int computeIDflag   (float mva,int ptId,int etaId);

    const variables_list_t & getVariables() const { return variables_; };
    
protected:

    void setup(); 
    void runMva(); 
    void bookReader();  
    void resetVariables();
    void initVariables();

    
    PileupJetIdentifier internalId_;
    variables_list_t variables_;

    std::unique_ptr<TMVA::Reader> reader_, reader_jteta_0_2_, reader_jteta_2_2p5_, reader_jteta_2p5_3_, reader_jteta_3_5_;
    std::string    tmvaWeights_, tmvaWeights_jteta_0_2_, tmvaWeights_jteta_2_2p5_, tmvaWeights_jteta_2p5_3_, tmvaWeights_jteta_3_5_, tmvaMethod_; 
    std::vector<std::string>  tmvaVariables_, tmvaVariables_jteta_0_3_, tmvaVariables_jteta_3_5_;
    std::vector<std::string>  tmvaSpectators_;
    std::map<std::string,std::string>  tmvaNames_;
    
    Int_t   version_;
    Float_t impactParTkThreshod_;
    bool    cutBased_; 
    bool    etaBinnedWeights_;
    bool runMvas_;
    Float_t mvacut_     [3][4][4]; //Keep the array fixed
    Float_t rmsCut_     [3][4][4]; //Keep the array fixed
    Float_t betaStarCut_[3][4][4]; //Keep the array fixed
};

#endif