In This Package:

DsRpcModel.cc

Go to the documentation of this file.

#include "DsRpcModel.h"

//#include "G4OpticalPhoton.hh"
#include "G4MuonPlus.hh"
#include "G4MuonMinus.hh"
#include "G4FastTrack.hh"
#include "G4TransportationManager.hh"
#include "G4Navigator.hh"
#include "G4VSensitiveDetector.hh"

class G4Hooks {
public:
    G4Hooks()
        : fakeStep(new G4Step)
        , touchable(new G4TouchableHistory )
        , navi(new G4Navigator) 
        , setup(false) { }

    ~G4Hooks() {
        delete fakeStep;
        delete navi;
        
    }

    G4Step * fakeStep;
    G4TouchableHandle touchable;
    G4Navigator* navi;
    bool setup;
};

DsRpcModel::DsRpcModel(const G4String& name, G4Envelope* volume, G4bool unique)
    : G4VFastSimulationModel(name,volume,unique)
    , m_hooks(new G4Hooks)
{
    G4cout << "DsRpcModel " << name <<  " created" << G4endl;
}

DsRpcModel::DsRpcModel(const G4String& name)
    : G4VFastSimulationModel(name)
    , m_hooks(new G4Hooks)
{
    G4cout << "DsRpcModel " << name <<  " created" << G4endl;
}

DsRpcModel::~DsRpcModel ()
{
    delete m_hooks; m_hooks = 0;
}

    
G4bool DsRpcModel::IsApplicable(const G4ParticleDefinition& particle)
{
    G4cout << "DsRpcModel " << particle.GetParticleName() << G4endl;
    return &particle == G4MuonPlus::Definition()|| &particle == G4MuonMinus::Definition();
}

G4bool DsRpcModel::ModelTrigger(const G4FastTrack & /*track*/)
{
    // Veni, vidi, trigi
    return true;
}

void DsRpcModel::DoIt(const G4FastTrack& fastTrack, G4FastStep& fastStep)
{
    const G4Track* track = fastTrack.GetPrimaryTrack();
    double energy = track->GetKineticEnergy();

    fastStep.ProposeTrackStatus(fStopAndKill);
    fastStep.ProposePrimaryTrackPathLength(0.0);
    fastStep.ProposeTotalEnergyDeposited(energy);
}

---