PileupReweightingAlg.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/

 
 
#ifndef ASG_ANALYSIS_ALGORITHMS__PILEUP_REWEIGHTING_ALG_H
#define ASG_ANALYSIS_ALGORITHMS__PILEUP_REWEIGHTING_ALG_H
 
#include <xAODEventInfo/EventInfo.h>
#include <AnaAlgorithm/AnaAlgorithm.h>
#include <AsgAnalysisInterfaces/IPileupReweightingTool.h>
#include <SelectionHelpers/OutOfValidityHelper.h>
#include <SystematicsHandles/SysReadHandle.h>
#include <SystematicsHandles/SysListHandle.h>
#include <SystematicsHandles/SysWriteDecorHandle.h>
#include <AsgTools/PropertyWrapper.h>
#include <AsgDataHandles/ReadHandleKey.h>
#include <AsgDataHandles/ReadHandle.h>
#include <AsgDataHandles/WriteDecorHandleKey.h>
 
namespace CP
{
 
  class PileupReweightingAlg final : public EL::AnaAlgorithm
  {
  public:
    using EL::AnaAlgorithm::AnaAlgorithm;
    StatusCode initialize () override;
    StatusCode execute (const EventContext& ctx) override;
 
 

  private:
    ToolHandle<IPileupReweightingTool> m_pileupReweightingTool {this, "pileupReweightingTool", "PileupReweightingTool", "the pileup reweighting tool we apply"};

  private:
    SysListHandle m_systematicsList {this};

  private:
    CP::SysWriteDecorHandle<float> m_weightDecorator{
        this, "pileupWeightDecoration", "PileupWeight_%SYS%", "the decoration for the pileup weight"};

  private:
    CP::SysReadHandle<xAOD::EventInfo> m_eventInfoHandle{
        this, "eventInfo", "EventInfo", "the input EventInfo object"};

  private:
    SG::ReadHandleKey<xAOD::EventInfo> m_baseEventInfoName {this, "baseEventInfo", "EventInfo",       "The name of the original event info. The non-systematic dependent decorations will be applied to this "
      "object so it should be at least a base of the shallow copies read in by the 'eventInfo' handle. "
      "The default (and strongly recommended behaviour) is to leave all of these pointed at the central 'EventInfo' object!"};

  private:
    SG::WriteDecorHandleKey<xAOD::EventInfo> m_correctedScaledAverageMuDecorator
    { this, "correctedScaledAverageMuDecoration", m_baseEventInfoName, "",
        "the decoration for the corrected and scaled average interactions per crossing" };

  private:
    SG::WriteDecorHandleKey<xAOD::EventInfo> m_correctedActualMuDecorator
    { this, "correctedActualMuDecoration", m_baseEventInfoName, "",
        "the decoration for the corrected actual interactions per crossing" };

  private:
    SG::WriteDecorHandleKey<xAOD::EventInfo> m_correctedScaledActualMuDecorator
    { this, "correctedScaledActualMuDecoration", m_baseEventInfoName, "",
        "the decoration for the corrected and scaled actual interactions per crossing" };

  private:
    OutOfValidityHelper m_outOfValidity {this};
 
    SG::WriteDecorHandleKey<xAOD::EventInfo> m_decRRNKey
    { this, "RandomRunNumberKey", m_baseEventInfoName, "RandomRunNumber",
        "Name for the RandomRunNumber decoration" };
 
    SG::WriteDecorHandleKey<xAOD::EventInfo> m_decRLBNKey
    { this, "RandomLumiBlockNumberKey", m_baseEventInfoName, "RandomLumiBlockNumber",
        "Name for the RandomLumiBlockNumber decoration" }; // unsigned int
 
    SG::WriteDecorHandleKey<xAOD::EventInfo> m_decHashKey
    { this, "PRWHashKey", m_baseEventInfoName, "PRWHash",
        "Name for the PRWHash decoration" };  // uint64_t
  };
}
 
#endif