Rivet_i.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef RIVET_I_H
#define RIVET_I_H
 
#include "AthenaBaseComps/AthAlgorithm.h"
#include "CxxUtils/checker_macros.h"
#include "xAODEventInfo/EventInfo.h"
#include "GaudiKernel/ToolHandle.h"
#include "GenInterfaces/IxAODtoHepMCTool.h"
 
#include "Rivet/AnalysisHandler.hh"
 
#include "AtlasHepMC/GenEvent.h"
 
#include <vector>
#include <string>
#include <memory>
 
class ISvcLocator;
class IxAODtoHepMCTool;
//class ITHistSvc;
 
 
class Rivet_i : public AthAlgorithm {
public:
 
  Rivet_i(const std::string &name, ISvcLocator *pSvcLocator);
 
  virtual StatusCode initialize ATLAS_NOT_THREAD_SAFE () override;
 
  virtual StatusCode execute() override;
 
  virtual StatusCode finalize() override;
 
 
private:
 
//  StatusCode regHist(const AIDA::IDataPointSet& dps, const std::string& path);
 
//  StatusCode regGraph(const AIDA::IDataPointSet& dps, const std::string& path);
 
  // Check and potentially modify events for correct units, beam particles, ...
  std::unique_ptr<HepMC::GenEvent> checkEvent(const HepMC::GenEvent& event, const EventContext& ctx);
 
  // Utility method to convert units of the event
  void MeV2GeV(HepMC::GenEvent& event);
 
  //ServiceHandle<ITHistSvc> m_histSvc;
 
  std::string m_stream;
 
  std::string m_file;
 
  //specify a pre-existing yoda file to initialize from
  std::string m_preload;
 
  std::string m_anapath;
 
  bool m_ignorebeams;
 
  bool m_needsConversion;
 
  bool m_patchBeams;
 
  ToolHandle<IxAODtoHepMCTool> m_xAODtoHepMCTool{this, "HepMCTool", "xAODtoHepMCTool"};
 
  bool m_doRootHistos;
 
  std::string m_runname;
 
  std::string m_genEventKey;
 
  Rivet::AnalysisHandler* m_analysisHandler;
 
  std::vector<std::string> m_analysisNames;
 
  double m_crossSection;
 
  double m_crossSection_uncert;
 
  bool m_init;
 
  bool m_skipweights;
 
  std::string m_matchWeights;
 
  std::string m_unmatchWeights;
 
  std::string m_nominalWeightName;
 
  double m_weightcap;
 
  SG::ReadHandleKey<xAOD::EventInfo> m_evtInfoKey{this
      , "EventInfo"
      , "EventInfo"
      , "ReadHandleKey for xAOD::EventInfo" };
};
 
#endif