Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 // ------------------------------------------- 28 // GEANT 4 - Underground Dark Matter Detecto 29 // 30 // For information related to this code c 31 // e-mail: alexander.howard@cern.ch 32 // ------------------------------------------- 33 // Comments 34 // 35 // Underground Advanced 36 // by A. Howard and H. Araujo 37 // (27th November 2001) 38 // 39 // MaxTimeCuts program 40 // ------------------------------------------- 41 42 #include "DMXMaxTimeCuts.hh" 43 44 #include "G4PhysicalConstants.hh" 45 #include "G4Step.hh" 46 #include "G4UserLimits.hh" 47 #include "G4VParticleChange.hh" 48 49 DMXMaxTimeCuts::DMXMaxTimeCuts(const G4String& 50 : DMXSpecialCuts(aName) 51 { 52 if (verboseLevel>1) { 53 G4cout << GetProcessName() << " is create 54 } 55 SetProcessType(fUserDefined); 56 } 57 58 DMXMaxTimeCuts::~DMXMaxTimeCuts() 59 {} 60 61 DMXMaxTimeCuts::DMXMaxTimeCuts(DMXMaxTimeCuts& 62 : DMXSpecialCuts() 63 {} 64 65 66 G4double DMXMaxTimeCuts::PostStepGetPhysicalIn 67 const G4Track& aT 68 G4double , 69 G4ForceCondition* condition 70 ) 71 { 72 // condition is set to "Not Forced" 73 *condition = NotForced; 74 75 G4double proposedStep = DBL_MAX; 76 // get the pointer to UserLimits 77 G4UserLimits* pUserLimits = aTrack.GetVolum 78 const G4DynamicParticle* aParticle = aTrack 79 80 // can apply cuts for specific particles - 81 // G4ParticleDefinition* aParticleDef = a 82 83 // G4cout << " Time: " << pUserLimits->Ge 84 85 if (pUserLimits) { 86 G4double temp = DBL_MAX; 87 //max time limit 88 G4double dTime= (pUserLimits->GetUserMaxT 89 if (dTime < 0. ) { 90 proposedStep = 0.; 91 } else { 92 G4double beta = (aParticle->GetTotalMom 93 temp = beta*c_light*dTime; 94 if (proposedStep > temp) proposedStep = 95 } 96 97 } 98 return proposedStep; 99 } 100