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 // Hadrontherapy advanced example for Geant4 27 // See more at: https://twiki.cern.ch/twiki/bi 28 29 #include "G4SteppingManager.hh" 30 #include "G4TrackVector.hh" 31 #include "HadrontherapySteppingAction.hh" 32 #include "G4ios.hh" 33 #include "G4SteppingManager.hh" 34 #include "G4Track.hh" 35 #include "G4Step.hh" 36 #include "G4StepPoint.hh" 37 #include "G4TrackStatus.hh" 38 #include "G4TrackVector.hh" 39 #include "G4ParticleDefinition.hh" 40 #include "G4ParticleTypes.hh" 41 #include "G4UserEventAction.hh" 42 #include "G4TransportationManager.hh" 43 #include "G4VSensitiveDetector.hh" 44 #include "HadrontherapyRunAction.hh" 45 #include "HadrontherapyMatrix.hh" 46 #include "G4SystemOfUnits.hh" 47 48 ////////////////////////////////////////////// 49 HadrontherapySteppingAction::HadrontherapyStep 50 { 51 runAction = run; 52 } 53 54 ////////////////////////////////////////////// 55 HadrontherapySteppingAction::~HadrontherapySte 56 { 57 } 58 59 ////////////////////////////////////////////// 60 void HadrontherapySteppingAction::UserStepping 61 { 62 63 // The followings are calls to usefuls inf 64 // Please, comment out them if want to use 65 66 // G4Track* theTrack = aStep->GetTrack(); 67 68 G4StepPoint* PreStep = aStep->GetPreStepPo 69 G4StepPoint* PostStep = aStep->GetPostStep 70 71 G4TouchableHandle touchPreStep = PreStep-> 72 G4TouchableHandle touchPostStep = PostStep 73 74 //G4double PreStepX =PreStep->GetPosition( 75 //G4double PreStepY =PreStep->GetPosition( 76 //G4double PreStepZ =PreStep->GetPosition( 77 78 //G4double PostStepX =PostStep->GetPositio 79 //G4double PostStepY =PostStep->GetPositio 80 //G4double PostStepZ =PostStep->GetPositi 81 82 //To get the current volume: 83 G4VPhysicalVolume* volumePre = touchPreSte 84 //G4VPhysicalVolume* volumePost =touchPost 85 86 //To get its name: 87 G4String namePre = volumePre->GetName(); 88 89 90 // positions in the global coordinate syst 91 //G4ThreeVector posPreStep = PreStep->Get 92 //G4ThreeVector posPostStep = PostStep->Ge 93 94 //G4int eventNum = G4RunManager::GetRunMan 95 96 //G4double parentID =aStep->GetTrack()->Ge 97 //G4double trackID =aStep->GetTrack()->Get 98 99 G4double eKin = aStep -> GetPreStepPoint() 100 101 G4double PosX = aStep->GetTrack()->GetPosi 102 G4double PosY = aStep->GetTrack()->GetPosi 103 G4double PosZ = aStep->GetTrack()->GetPosi 104 105 G4String volume= aStep->GetTrack()->GetVo 106 G4Track* theTrack = aStep->GetTrack(); 107 108 //G4String material= aStep -> GetTrack() - 109 //G4cout << "material " << material << G 110 //G4String volume= aStep->GetTrack()->Get 111 //G4String pvname= pv-> GetName(); 112 113 G4String particleName = aStep->GetTrack()- 114 115 G4double momentumX = aStep->GetTrack()->G 116 G4double momentumY = aStep->GetTrack()->G 117 G4double momentumZ = aStep->GetTrack()->G 118 119 120 G4ParticleDefinition *particleDef = theTra 121 G4int pdg = particleDef ->GetPDGEncoding() 122 123 if(namePre == "VirtualLayer") 124 { 125 std::ofstream WriteDataIn("Virtual_Lay 126 WriteDataIn 127 128 << eKin <<" " // 1 129 << PosX <<" " // 2 130 << PosY <<" " // 3 131 << PosZ <<" " // 4 132 << momentumX <<" " // 5 133 << momentumY <<" " // 6 134 << momentumZ <<" " // 7 135 << pdg 136 //<< theTrack << '\t' << " 137 138 << G4endl; 139 140 theTrack -> SetTrackStatus(fKillTrackA 141 142 143 } 144 145 146 147 148 } 149 150 151 152