Message: how to divide a run into two parts? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None how to divide a run into two parts? 

Forum: Event and Track Management
Date: 23 May, 2007
From: Xinjian Tan <Xinjian Tan>

Hello everybody
    in order to save the calculation time, i'd like to divide a run into 
two parts. I save the particles' information (positon, kinetic Energy, 
momentum direction, etc) in the first part when some condition are satisfied.
Reading the file in the second part, generate the particle according the particle
information saved in the file, and continue the simulation. So when I modified the
detector geometry, only the second part should be rerun. 

    but the simulated result by this method is always nearly 1.25 times of that by the normal 
simulation. I don't know why.

    realization detail:
1. save particle information:
    In G4UserSteppingAction::UserSteppingAction(aStep)

    position: aStep->GetPostStepPoint()->GetPosition();
    momentum Direction: aStep->GetTrack()->GetMomentumDirction();
    energy : aStep->GetTrack()->GetKineticEnergy();
    particle name: aStep->GetTrack->GetDefinition->GetParticleName();

2. read file
    dirive my class HEGammaDetRunManager from class G4RunManager, rewrite the funciton
DoEventLoop(). in this Function, i read the file and set the primary information according
the particles that saved in the file.

    Could you help me to identify what is the problem? Thanks^_^

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: how to divide a run into two parts?_ additional explaination   (Xinjian Tan - 23 May, 2007)
 Add Message Add Message
to: "how to divide a run into two parts?"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews

[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]