An example of how to quickly make some plots at a certain point in the cutflow. More...

#include <JetResponsePlots.h>

Inheritance diagram for top::JetResponsePlots:

Collaboration diagram for top::JetResponsePlots:

Public Member Functions
	JetResponsePlots (const std::string &name, TFile outputFile, const std::string &params, std::shared_ptr< top::TopConfig > config, EL::Worker wk=nullptr)
	Setup some example plots. More...

bool	apply (const top::Event &event) const override
	Fill the histograms. More...

void	FillHistograms (const double w_event, const top::Event &event) const
	Helper function to fill the histograms. More...

std::string	name () const override
	Return the name for the cutflow table. More...

virtual bool	applyParticleLevel (const top::ParticleLevelEvent &) const
	This does stuff based on the information in a particle level event. More...

Private Attributes
std::shared_ptr< PlotManager >	m_hists

std::size_t	m_nominalHashValue

float	m_deltaR

int	m_bins

float	m_min

float	m_max

std::vector< double >	m_ptBinning

std::shared_ptr< top::TopConfig >	m_config

ToolHandle< PMGTools::IPMGTruthWeightTool >	m_PMGTruthWeights

Detailed Description

An example of how to quickly make some plots at a certain point in the cutflow.

Definition at line 31 of file JetResponsePlots.h.

Constructor & Destructor Documentation

◆ JetResponsePlots()

top::JetResponsePlots::JetResponsePlots	(	const std::string &	name,
		TFile *	outputFile,
		const std::string &	params,
		std::shared_ptr< top::TopConfig >	config,
		EL::Worker *	wk = `nullptr`
	)

Setup some example plots.

Add a bunch of histograms.

Parameters

name	The name of the directory to store histograms in, in the output file. e.g. you might have ee, mumu and emu.
outputFile	The output file. Needs setting up at the very start so that we can attach the files.
params	The arguments, e.g. for the binning of the plots.
config	Instance of TopConfig
wk	Only used by EventLoop, ok as nullptr as default.

Definition at line 24 of file JetResponsePlots.cxx.

                                                    :
     m_nominalHashValue(0),
     m_deltaR(0.3),
     m_bins(30),
     m_min(-3),
     m_max(3),
     m_config(config),
     m_PMGTruthWeights(nullptr) {
  
     CP::SystematicSet nominal;
     m_nominalHashValue = nominal.hash();
  
     //retrieve PMGTruthWeights
     static const std::string truthWeightToolName = "PMGTruthWeightTool";
     if (asg::ToolStore::contains<PMGTools::IPMGTruthWeightTool>(truthWeightToolName)) m_PMGTruthWeights =
         asg::ToolStore::get<PMGTools::IPMGTruthWeightTool>(truthWeightToolName);
     else {
       std::unique_ptr<PMGTools::PMGTruthWeightTool> pmgtruthtool = std::make_unique<PMGTools::PMGTruthWeightTool>(truthWeightToolName);
       top::check(pmgtruthtool->initialize(), "Failed to initialize " + truthWeightToolName);
       m_PMGTruthWeights = pmgtruthtool.release();
     }
     //decoding arguments
     std::istringstream stream(params);
     std::string s;
     while (stream >> s) {
       if (s.substr(0,7) == "deltaR:") {
         m_deltaR = std::stof(s.substr(7, s.size() - 7));
         m_config->setJetResponseMatchingDeltaR(m_deltaR);
       } else if (s.substr(0,5) == "bins:") {
         m_bins = std::stoi(s.substr(5, s.size() - 5));
       } else if (s.substr(0,4) == "min:") {
         m_min = std::stof(s.substr(4, s.size() - 4));
       } else if (s.substr(0,4) == "max:") {
         m_max = std::stof(s.substr(4, s.size() - 4));
       } else if (s.substr(0,10) == "ptBinning:") {
         const std::string tmp = s.substr(10, s.size() - 4);
         // split by the delimiter
         std::vector<std::string> split;
         tokenize(tmp, split, ",");
         for (const std::string& istring : split) {
           m_ptBinning.emplace_back(std::stof(istring));
         }
       } else {
         throw std::runtime_error{"ERROR: Can't understand argument " + s + "For JETRESPONSEPLOTS"};
       }
     }
  
     if (m_ptBinning.size() < 2) {
       throw std::invalid_argument{"pT binning (x axis) size is smaller than 2"};
     }
  
     m_hists = std::make_shared<PlotManager>(name + "/JetResponsePlots", outputFile, wk);
     m_hists->addHist("JetResponse",      ";p_{T}^{truth} [GeV];p_{T}^{reco}/p_{T}^{truth};Events", m_ptBinning.size() - 1, m_ptBinning.data(), m_bins, m_min, m_max);
     m_hists->addHist("JetResponse_bb",   ";p_{T}^{truth} [GeV];p_{T}^{reco}/p_{T}^{truth};Events", m_ptBinning.size() - 1, m_ptBinning.data(), m_bins, m_min, m_max);
     m_hists->addHist("JetResponse_b",    ";p_{T}^{truth} [GeV];p_{T}^{reco}/p_{T}^{truth};Events", m_ptBinning.size() - 1, m_ptBinning.data(), m_bins, m_min, m_max);
     m_hists->addHist("JetResponse_c",    ";p_{T}^{truth} [GeV];p_{T}^{reco}/p_{T}^{truth};Events", m_ptBinning.size() - 1, m_ptBinning.data(), m_bins, m_min, m_max);
     m_hists->addHist("JetResponse_tau",  ";p_{T}^{truth} [GeV];p_{T}^{reco}/p_{T}^{truth};Events", m_ptBinning.size() - 1, m_ptBinning.data(), m_bins, m_min, m_max);
     m_hists->addHist("JetResponse_q",    ";p_{T}^{truth} [GeV];p_{T}^{reco}/p_{T}^{truth};Events", m_ptBinning.size() - 1, m_ptBinning.data(), m_bins, m_min, m_max);
     m_hists->addHist("JetResponse_g",    ";p_{T}^{truth} [GeV];p_{T}^{reco}/p_{T}^{truth};Events", m_ptBinning.size() - 1, m_ptBinning.data(), m_bins, m_min, m_max);
     m_hists->addHist("JetResponse_other",";p_{T}^{truth} [GeV];p_{T}^{reco}/p_{T}^{truth};Events", m_ptBinning.size() - 1, m_ptBinning.data(), m_bins, m_min, m_max);
   }

Member Function Documentation

◆ apply()

bool top::JetResponsePlots::apply ( const top::Event & event ) const

overridevirtual

Fill the histograms.

Returns: True because it doesn't select any events.

Implements top::EventSelectorBase.

Definition at line 90 of file JetResponsePlots.cxx.

                                                           {
     //only MC
     if (!top::isSimulation(event)) return true;
  
     // only nominal
     if (event.m_hashValue != m_nominalHashValue) return true;
  
     // do loose or tight events only if requested
     if (event.m_isLoose && !m_config->doLooseEvents()) return true;
  
     if (!event.m_isLoose && !m_config->doTightEvents()) return true;
  
     const double nominalWeight = event.m_info->auxdata<float>("AnalysisTop_eventWeight");
  
     FillHistograms(nominalWeight, event);
       
     return true;
   }

◆ applyParticleLevel()

virtual bool top::EventSelectorBase::applyParticleLevel ( const top::ParticleLevelEvent & ) const

inlinevirtualinherited

This does stuff based on the information in a particle level event.

The idea is that you implement this to return either true or false, based on the information held within the top::ParticleLevelEvent. If this function returns true, then the event is kept, otherwise it is removed. The function has a default implementation (which returns true) because it is expected that many EventSelector objects do not operate on ParticleLevelEvent objects.

Parameters

top::ParticleLevelEvent	the current particle level event.
true	if the event should be kept (i.e. it passed the selector criteria), false otherwise.

Reimplemented in top::JetNGhostSelector, top::PrintEventSelector, top::PseudoTopRecoRun, top::NElectronNMuonTightSelector, top::NElectronNMuonSelector, top::NFwdElectronSelector, top::HTSelector, top::OSLeptonTightSelector, top::MLLSelector, top::MWTSelector, top::NElectronTightSelector, top::NFwdElectronTightSelector, top::NMuonTightSelector, top::OSLeptonSelector, top::METMWTSelector, top::METSelector, top::MLLWindow, top::NElectronSelector, top::NJetSelector, top::NMuonSelector, top::NPhotonSelector, top::NSoftMuonSelector, top::NTauSelector, top::SSLeptonTightSelector, top::SSLeptonSelector, top::ParticleLevelSelector, top::RecoLevelSelector, top::NVarRCJetSelector, top::NLargeJetSelector, and top::NRCJetSelector.

Definition at line 73 of file EventSelectorBase.h.

73 {return true;}

◆ FillHistograms()

void top::JetResponsePlots::FillHistograms	(	const double	w_event,
		const top::Event &	event
	)		const

Helper function to fill the histograms.

Definition at line 109 of file JetResponsePlots.cxx.

                                                                      {
  
     for (const auto* const jetPtr : event.m_jets) {
       int jet_flavor = -99;
       bool status = jetPtr->getAttribute<int>("HadronConeExclExtendedTruthLabelID", jet_flavor);
       if (!status) continue;
  
       float matchedPt = -9999;
       status = jetPtr->getAttribute<float>("AnalysisTop_MatchedTruthJetPt", matchedPt);
       if (!status) continue;
       
       // protection against division by zero
       if (matchedPt < 1e-6) continue;
  
       const float response = jetPtr->pt() / matchedPt;
  
       static_cast<TH2D*>(m_hists->hist("JetResponse"))->Fill(matchedPt/1e3, response, w_event);
       if (jet_flavor == 55) {
         static_cast<TH2D*>(m_hists->hist("JetResponse_bb"))->Fill(matchedPt/1e3, response, w_event);
       } else if (jet_flavor == 5 || jet_flavor == 54) {
         static_cast<TH2D*>(m_hists->hist("JetResponse_b"))->Fill(matchedPt/1e3, response, w_event);
       } else if (jet_flavor == 4 || jet_flavor == 44) {
         static_cast<TH2D*>(m_hists->hist("JetResponse_c"))->Fill(matchedPt/1e3, response, w_event);
       } else if (jet_flavor == 15) {
         static_cast<TH2D*>(m_hists->hist("JetResponse_tau"))->Fill(matchedPt/1e3, response, w_event);
       } else {
         status = jetPtr->getAttribute<int>("PartonTruthLabelID", jet_flavor);
         if (!status) continue;
         if (jet_flavor >= 1 && jet_flavor <= 3) {
           static_cast<TH2D*>(m_hists->hist("JetResponse_q"))->Fill(matchedPt/1e3, response, w_event);
         } else if (jet_flavor == 21) {
           static_cast<TH2D*>(m_hists->hist("JetResponse_g"))->Fill(matchedPt/1e3, response, w_event);
         } else {
           static_cast<TH2D*>(m_hists->hist("JetResponse_other"))->Fill(matchedPt/1e3, response, w_event);
         }
       }
     }
   }

◆ name()

std::string top::JetResponsePlots::name ( ) const

overridevirtual

Return the name for the cutflow table.

Returns: The word JETRESPONSEPLOTS.

Implements top::EventSelectorBase.

Definition at line 149 of file JetResponsePlots.cxx.

                                          {
     return "JETRESPONSEPLOTS";
   }

Member Data Documentation

◆ m_bins

int top::JetResponsePlots::m_bins

private

Definition at line 82 of file JetResponsePlots.h.

◆ m_config

std::shared_ptr<top::TopConfig> top::JetResponsePlots::m_config

private

Definition at line 88 of file JetResponsePlots.h.

◆ m_deltaR

float top::JetResponsePlots::m_deltaR

private

Definition at line 81 of file JetResponsePlots.h.

◆ m_hists

std::shared_ptr<PlotManager> top::JetResponsePlots::m_hists

private

Definition at line 75 of file JetResponsePlots.h.

◆ m_max

float top::JetResponsePlots::m_max

private

Definition at line 84 of file JetResponsePlots.h.

◆ m_min

float top::JetResponsePlots::m_min

private

Definition at line 83 of file JetResponsePlots.h.

◆ m_nominalHashValue

std::size_t top::JetResponsePlots::m_nominalHashValue

private

Definition at line 78 of file JetResponsePlots.h.

◆ m_PMGTruthWeights

ToolHandle<PMGTools::IPMGTruthWeightTool> top::JetResponsePlots::m_PMGTruthWeights

private

Definition at line 91 of file JetResponsePlots.h.

◆ m_ptBinning

std::vector<double> top::JetResponsePlots::m_ptBinning

private

Definition at line 85 of file JetResponsePlots.h.

The documentation for this class was generated from the following files:

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ JetResponsePlots()

Member Function Documentation

◆ apply()

◆ applyParticleLevel()

◆ FillHistograms()

◆ name()

Member Data Documentation

◆ m_bins

◆ m_config

◆ m_deltaR

◆ m_hists

◆ m_max

◆ m_min

◆ m_nominalHashValue

◆ m_PMGTruthWeights

◆ m_ptBinning