CharmTagger.cc

Go to the documentation of this file.

#include "RecoBTag/CTagging/interface/CharmTagger.h"

#include "DataFormats/BTauReco/interface/CandSoftLeptonTagInfo.h"
#include "DataFormats/BTauReco/interface/CandIPTagInfo.h"
#include "DataFormats/BTauReco/interface/SecondaryVertexTagInfo.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "CondFormats/DataRecord/interface/GBRWrapperRcd.h"

#include <iostream>
#include <vector>
#include <algorithm>
#include <map>
#include <iostream>

CharmTagger::CharmTagger(const edm::ParameterSet & configuration):
    sl_computer_(configuration.getParameter<edm::ParameterSet>("slComputerCfg")),
  mva_name_( configuration.getParameter<std::string >("mvaName") ),
  use_condDB_(configuration.getParameter<bool>("useCondDB")),
  gbrForest_label_(configuration.getParameter<std::string>("gbrForestLabel")),
  weight_file_(configuration.getParameter<edm::FileInPath>("weightFile")),
  use_GBRForest_(configuration.getParameter<bool>("useGBRForest")),
  use_adaBoost_(configuration.getParameter<bool>("useAdaBoost")),
  defaultValueNoTracks_(configuration.getParameter<bool>("defaultValueNoTracks"))
{
    vpset vars_definition = configuration.getParameter<vpset>("variables");
    for(auto &var : vars_definition) {
        MVAVar mva_var;
        mva_var.name = var.getParameter<std::string>("name");
        mva_var.id = reco::getTaggingVariableName(
            var.getParameter<std::string>("taggingVarName")
            );
        mva_var.has_index = var.existsAs<int>("idx") ;
        mva_var.index = mva_var.has_index ? var.getParameter<int>("idx") : 0;
        mva_var.default_value = var.getParameter<double>("default");

        variables_.push_back(mva_var);
    }

    uses(0, "pfImpactParameterTagInfos");
    uses(1, "pfInclusiveSecondaryVertexFinderCvsLTagInfos");
    uses(2, "softPFMuonsTagInfos");
    uses(3, "softPFElectronsTagInfos");
}

void CharmTagger::initialize(const JetTagComputerRecord & record)
{
    mvaID_.reset(new TMVAEvaluator());

    std::vector<std::string> variable_names;
    variable_names.reserve(variables_.size());

    for(auto &var : variables_) {
        variable_names.push_back(var.name);
    }
    std::vector<std::string> spectators;

  if(use_condDB_) {
        const GBRWrapperRcd & gbrWrapperRecord = record.getRecord<GBRWrapperRcd>();

        edm::ESHandle<GBRForest> gbrForestHandle;
        gbrWrapperRecord.get(gbrForest_label_.c_str(), gbrForestHandle);

        mvaID_->initializeGBRForest(
            gbrForestHandle.product(), variable_names, 
            spectators, use_adaBoost_
            );
  }
  else {
    mvaID_->initialize(
            "Color:Silent:Error", mva_name_,
            weight_file_.fullPath(), variable_names, 
            spectators, use_GBRForest_, use_adaBoost_
            );
  }
}

CharmTagger::~CharmTagger()
{
}

float CharmTagger::discriminator(const TagInfoHelper & tagInfo) const {
  // default value, used if there are no leptons associated to this jet
  const reco::CandIPTagInfo & ip_info = tagInfo.get<reco::CandIPTagInfo>(0);
    const reco::CandSecondaryVertexTagInfo & sv_info = tagInfo.get<reco::CandSecondaryVertexTagInfo>(1);
    const reco::CandSoftLeptonTagInfo& softmu_info = tagInfo.get<reco::CandSoftLeptonTagInfo>(2);
    const reco::CandSoftLeptonTagInfo& softel_info = tagInfo.get<reco::CandSoftLeptonTagInfo>(3);
    reco::TaggingVariableList vars = sl_computer_(ip_info, sv_info, softmu_info, softel_info);

    // Loop over input variables
    std::map<std::string, float> inputs;
    for(auto &mva_var : variables_){
        //vectorial tagging variable
        if(mva_var.has_index){
            std::vector<float> vals = vars.getList(mva_var.id, false);
            inputs[mva_var.name] = (vals.size() > mva_var.index) ? vals[mva_var.index] : mva_var.default_value;
        }
        //single value tagging var
        else {
            inputs[mva_var.name] = vars.get(mva_var.id, mva_var.default_value);
        }
    }

    //get the MVA output
    bool notracks = (vars.get(reco::btau::jetNTracks,0) == 0);
    float tag = mvaID_->evaluate(inputs);
    if (defaultValueNoTracks_ && notracks){tag = -2;} // if no tracks available, put value at -2 (only for Phase I)
    return tag;
}