/scratch0/jsnuveri/BDSIM/BDSIMgit/bdsim/src/BDSTerminatorUserLimits.cc

/* BDSIM
   Author: L. Nevay,
   Last modified 15/04/2014
   Copyright (c) 2014.  ALL RIGHTS RESERVED.
*/

/*

This class is to give dynamic user limits to a certain volume.
In the case of a ring, we want to cut ALL particles at the end 
of the ring after a certain number of turns.
*/

#include "G4ios.hh"
#include "BDSTerminatorUserLimits.hh"
#include "G4UserLimits.hh"
#include "BDSGlobalConstants.hh"
#include "BDSExecOptions.hh"
#include "G4Track.hh"
#include "BDSDebug.hh"

// basic inheritance - just use everything normally from G4UserLimits but 
// replace one function in inherited class
// default values are defined in G4UserLimits so all particles continue
BDSTerminatorUserLimits::BDSTerminatorUserLimits(G4double ustepMax,
                                                 G4double utrakMax,
                                                 G4double utimeMax,
                                                 G4double uekinMin,
                                                 G4double urangMin):
  G4UserLimits(ustepMax,
                 utrakMax,
                 utimeMax,
                 uekinMin,
                 urangMin), keeprunningEK(0.0), stoprunningEK(DBL_MAX)
{
  verbose = BDSExecOptions::Instance()->GetVerbose();
}

BDSTerminatorUserLimits::BDSTerminatorUserLimits(const G4String& type,
                                                 G4double ustepMax,
                                                 G4double utrakMax,
                                                 G4double utimeMax,
                                                 G4double uekinMin,
                                                 G4double urangMin):
  G4UserLimits(type,
               ustepMax,
               utrakMax,
               utimeMax,
               uekinMin,
               urangMin), keeprunningEK(0.0), stoprunningEK(DBL_MAX)
{
  verbose = BDSExecOptions::Instance()->GetVerbose();
}

inline G4double BDSTerminatorUserLimits::GetUserMinEkine(const G4Track& trk)
{
  // does the number of turns passed == number of turns to take
#ifdef BDSDEBUG
  G4cout << __METHOD_NAME__ << " turns taken : " << BDSGlobalConstants::Instance()->GetTurnsTaken() << G4endl;
#endif
  if ((BDSGlobalConstants::Instance()->GetTurnsTaken() == BDSGlobalConstants::Instance()->GetTurnsToTake())
      && trk.GetTrackLength()/CLHEP::m > 1.0*CLHEP::m)
    {
      //only stop if it's travelled at least one metre - protects against starting distribution
      //starting inside terminator
      if (verbose){
        G4cout << "Requested number of turns completed - stopping all particles" << G4endl;
        G4cout << "Track length: " << trk.GetTrackLength()/CLHEP::m << " m" << G4endl;
      }
      return stoprunningEK;
    } // yes - stop: return DBL_MAX eV so no particles will be tracked
  else
    {
      return keeprunningEK;
    } // no  - keep going: return 0 eV so all are tracked
}

---