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