FastSimModel.h

Go to the documentation of this file.

// Copyright 2019-2026 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
#ifndef O2_FASTSIM_MODEL_H_
#define O2_FASTSIM_MODEL_H_
 
 
#include "G4VFastSimulationModel.hh"
 
#include <vector>
 
class G4FastStep;
class G4FastTrack;
class G4ParticleDefinition;
 
namespace o2::fastsim
{
 
struct FastSimInput {
  int pdg = 0;
  double position[3] = {};   
  double direction[3] = {};  
  double kineticEnergy = 0.; 
  double mass = 0.;          
  double time = 0.;          
  double exitDistance = 0.;  
};
 
struct FastSimOutput {
  int pdg = 0;
  double position[3] = {}; 
  double momentum[3] = {}; 
  double time = 0.;        
};
 
constexpr double kSurfaceEpsilonCm = 1e-5;
 
class FastSimModel : public G4VFastSimulationModel
{
 public:
  FastSimModel(const G4String& name, const G4String& envelopeVolume, double minEnergyGeV);
 
  G4bool IsApplicable(const G4ParticleDefinition& particle) override;
  G4bool ModelTrigger(const G4FastTrack& fastTrack) override;
 
  void DoIt(const G4FastTrack& fastTrack, G4FastStep& fastStep) override;
 
 protected:
  virtual std::vector<FastSimOutput> sample(const FastSimInput& input) const = 0;
 
 private:
  int envelopeDepth(const G4Track* track) const;
 
  G4String mEnvelope;     
  double mMinEnergy = 0.; 
  mutable bool mWarned = false;
};
 
} // namespace o2::fastsim
 
#endif // O2_FASTSIM_MODEL_H_