Back to home page

EIC code displayed by LXR

 
 

    


Warning, file /geant4/examples/advanced/xray_telescope/include/XrayTelStepCut.hh was not indexed or was modified since last indexation (in which case cross-reference links may be missing, inaccurate or erroneous).

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 //
0027 // **********************************************************************
0028 // *                                                                    *
0029 // *                    GEANT 4 xray_telescope advanced example         *
0030 // *                                                                    *
0031 // * MODULE:            XrayTelStepCut.hh                               *
0032 // * -------                                                            *
0033 // *                                                                    *
0034 // * Version:           0.4                                             *
0035 // * Date:              06/11/00                                        *
0036 // * Author:            R Nartallo                                      *
0037 // * Organisation:      ESA/ESTEC, Noordwijk, THe Netherlands           *
0038 // *                                                                    *
0039 // **********************************************************************
0040 //
0041 // CHANGE HISTORY
0042 // --------------
0043 //
0044 // 06.11.2000 R.Nartallo
0045 // - First implementation of X-ray Telescope advanced example.
0046 // - Based on Chandra and XMM models
0047 //
0048 //
0049 // **********************************************************************
0050 
0051 #ifndef XrayTelStepCut_h
0052 #define XrayTelStepCut_h 1
0053 
0054 #include "G4ios.hh"
0055 #include "G4VDiscreteProcess.hh"
0056 #include "G4Step.hh"
0057 #include "globals.hh"
0058 
0059 class XrayTelStepCut : public G4VDiscreteProcess
0060 {
0061 public:     
0062 
0063   XrayTelStepCut(const G4String& processName ="UserStepCut" );
0064   XrayTelStepCut(XrayTelStepCut &);
0065 
0066   ~XrayTelStepCut();
0067 
0068   G4double PostStepGetPhysicalInteractionLength(
0069                         const G4Track& track,
0070                         G4double   previousStepSize,
0071                         G4ForceCondition* condition
0072                         );
0073 
0074   G4VParticleChange* PostStepDoIt(
0075                   const G4Track& ,
0076                   const G4Step& 
0077                   );
0078 
0079   void SetMaxStep(G4double);
0080 
0081 protected:
0082 
0083   // it is not needed here !
0084   G4double GetMeanFreePath(const G4Track& aTrack,
0085                G4double   previousStepSize,
0086                G4ForceCondition* condition
0087                );
0088 
0089                 
0090 private:
0091   
0092   // hide assignment operator as private 
0093   XrayTelStepCut & operator=(const XrayTelStepCut &right);
0094 
0095 private:
0096 
0097   G4double MaxChargedStep ;
0098 };
0099 
0100 // inlined function members implementation
0101 
0102 #include "G4Step.hh"
0103 #include "G4Track.hh"
0104 #include "G4UserLimits.hh"
0105 #include "G4VParticleChange.hh"
0106 #include "G4EnergyLossTables.hh"
0107 
0108 inline G4double XrayTelStepCut::PostStepGetPhysicalInteractionLength(
0109                                      const G4Track& aTrack,
0110                                      G4double,
0111                                      G4ForceCondition* condition
0112                                      )
0113 {
0114   // condition is set to "Not Forced"
0115   *condition = NotForced;
0116 
0117   G4double ProposedStep = DBL_MAX;
0118 
0119   if((MaxChargedStep > 0.) &&
0120      (aTrack.GetVolume() != NULL) &&
0121      (aTrack.GetVolume()->GetName() == "Absorber") &&
0122      (aTrack.GetDynamicParticle()->GetDefinition()->GetPDGCharge() != 0.))
0123     ProposedStep = MaxChargedStep ;
0124 
0125   return ProposedStep;
0126 }
0127 
0128 inline G4VParticleChange* XrayTelStepCut::PostStepDoIt(
0129                                const G4Track& aTrack,
0130                                const G4Step&
0131                                )
0132 {
0133   // do nothing
0134   aParticleChange.Initialize(aTrack);
0135   return &aParticleChange;
0136 }
0137 
0138 inline G4double XrayTelStepCut::GetMeanFreePath(const G4Track&,
0139                         G4double,
0140                         G4ForceCondition*
0141                         )
0142 {
0143   return 0.;
0144 }
0145 
0146 #endif
0147