AsgCutBookkeeperAlg.h

Go to the documentation of this file.

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

 
 
#ifndef ASG_ANALYSIS_ALGORITHMS__ASG_CUT_BOOKKEEPER_ALG_H
#define ASG_ANALYSIS_ALGORITHMS__ASG_CUT_BOOKKEEPER_ALG_H
 
#include <AnaAlgorithm/AnaAlgorithm.h>
#include <AsgServices/ServiceHandle.h>
#include <PATInterfaces/SystematicSet.h>
#include <PMGAnalysisInterfaces/IPMGTruthWeightTool.h>
#include <SystematicsHandles/ISystematicsSvc.h>
#include <xAODCutFlow/CutBookkeeper.h>
#include <AsgTools/PropertyWrapper.h>
 
namespace CP
{
 
  class AsgCutBookkeeperAlg final : public EL::AnaAlgorithm
  {
  public:
    using EL::AnaAlgorithm::AnaAlgorithm;
 
  public:
    StatusCode initialize () override;
 
  public:
    StatusCode fileExecute () override;
 
  public:
    StatusCode finalize () override;
 
  private:
    void processCutBookkeeper (const xAOD::CutBookkeeper *cbk,
                               size_t index);

  private:
    Gaudi::Property<bool> m_enableSystematics {this, "enableSystematics", false, "enable systematics"};

  private:
    ToolHandle<PMGTools::IPMGTruthWeightTool> m_truthWeightTool {this, "truthWeightTool", "PMGTools::PMGTruthWeightTool", "the truth weight tool"};

  private:
    ServiceHandle<ISystematicsSvc> m_systematics {this, "systematics", "SystematicsSvc", "systematics service"};

  private:
    Gaudi::Property<uint32_t> m_runNumber {this, "runNumber", 0, "the run number we are processing"};

  private:
    uint32_t m_mcChannelNumber {};

  private:
    struct WeightsGroup
    {
      float nEventsProcessed {};
      float sumOfWeights {};
      float sumOfWeightsSquared {};
    };

  private:
    std::vector<std::string> m_allowed_streams{"StreamRDO", "StreamAOD", "StreamEVGEN", "StreamEVNT"};

  private:
    std::unordered_map<size_t, WeightsGroup> m_weights;

  private:
    Gaudi::Property<std::string> m_histPattern {this, "histPattern", "CutBookkeeper_%DSID%_%RUN%_%SYS%", "the pattern for histogram names"};
  };
}
 
#endif