LDMXCascadeInterface.h

Go to the documentation of this file.

 
#ifndef SIMCORE_BERTINI_LDMXCASCADEINTERFACE_H
#define SIMCORE_BERTINI_LDMXCASCADEINTERFACE_H
 
// Include the hack header first to expose private Geant4 members
#include "Framework/Logger.h"
#include "SimCore/Bertini/G4BertiniHack.h"
#include "SimCore/Event/CascadeHistory.h"
 
class G4HadProjectile;
class G4Nucleus;
class G4HadFinalState;
 
namespace simcore {
namespace bertini {
 
class LDMXCascadeInterface : public G4CascadeInterface {
 public:
  LDMXCascadeInterface(const G4String& name = "LDMXBertiniCascade");
  virtual ~LDMXCascadeInterface();
 
  G4HadFinalState* ApplyYourself(const G4HadProjectile& projectile,
                                 G4Nucleus& targetNucleus) override;
 
  void setEnergyThreshold(double threshold) { energy_threshold_ = threshold; }
  double getEnergyThreshold() const { return energy_threshold_; }
 
  const ldmx::CascadeHistory* getLastCascadeHistory() const {
    return last_history_.empty() ? nullptr : &last_history_;
  }
 
  ldmx::CascadeHistory extractHistory() { return std::move(last_history_); }
 
  bool hasHistory() const { return !last_history_.empty(); }
 
  void setIncidentTrackId(int trackId) { incident_track_id_ = trackId; }
 
 private:
  void ensureCascadeHistoryExists();
 
  void captureHistory();
 
  void captureDeexcitationProducts(G4HadFinalState* finalState);
 
  double energy_threshold_{5000.0};
 
  int incident_track_id_{-1};
 
  ldmx::CascadeHistory last_history_;
 
  enableLogging("LDMXCascadeInterface")
};
 
}  // namespace bertini
}  // namespace simcore
 
#endif  // SIMCORE_BERTINI_LDMXCASCADEINTERFACE_H