| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 /// \file medical/dna/range/src/TrackingAction 26 /// \file medical/dna/range/src/TrackingAction.cc
27 /// \brief Implementation of the TrackingActio 27 /// \brief Implementation of the TrackingAction class
28 // 28 //
29 // $Id: TrackingAction.cc 78723 2014-01-20 10: 29 // $Id: TrackingAction.cc 78723 2014-01-20 10:32:17Z gcosmo $
30 // 30 //
31 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 33
34 #include "TrackingAction.hh" 34 #include "TrackingAction.hh"
35 <<
36 #include "DetectorConstruction.hh" 35 #include "DetectorConstruction.hh"
37 #include "HistoManager.hh" <<
38 #include "PrimaryGeneratorAction.hh" <<
39 #include "Run.hh" 36 #include "Run.hh"
>> 37 #include "PrimaryGeneratorAction.hh"
>> 38 #include "HistoManager.hh"
40 39
41 #include "G4RunManager.hh" 40 #include "G4RunManager.hh"
42 41
43 //....oooOO0OOooo........oooOO0OOooo........oo 42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44 43
45 TrackingAction::TrackingAction(PrimaryGenerato 44 TrackingAction::TrackingAction(PrimaryGeneratorAction* /*prim*/)
46 : G4UserTrackingAction(), << 45 :G4UserTrackingAction(),
47 // fPrimary(prim), << 46 //fPrimary(prim),
48 fpDetector(0) << 47 fpDetector(0)
49 { 48 {
50 fpDetector = dynamic_cast<const DetectorCons << 49 fpDetector =
51 G4RunManager::GetRunManager()->GetUserDete << 50 dynamic_cast<const DetectorConstruction*>(G4RunManager::GetRunManager()
>> 51 ->GetUserDetectorConstruction());
52 } 52 }
53 53
54 //....oooOO0OOooo........oooOO0OOooo........oo 54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
55 55
56 void TrackingAction::PreUserTrackingAction(con 56 void TrackingAction::PreUserTrackingAction(const G4Track* track)
57 { 57 {
58 fpDetector = dynamic_cast<const DetectorCons << 58 fpDetector =
59 G4RunManager::GetRunManager()->GetUserDete << 59 dynamic_cast<const DetectorConstruction*>(G4RunManager::GetRunManager()
>> 60 ->GetUserDetectorConstruction());
60 61
61 G4int trackID = track->GetTrackID(); << 62 G4int trackID = track->GetTrackID();
62 if (trackID == 1) fTrackLength = 0; << 63 if (trackID == 1) fTrackLength = 0;
63 64
64 G4AnalysisManager* analysisManager = G4Analy 65 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
65 66
66 G4ThreeVector pos = track->GetPosition(); << 67 G4ThreeVector pos = track->GetPosition();
67 G4double RNP = fpDetector->GetNPRadius() / C << 68 G4double RNP = fpDetector->GetNPRadius()/CLHEP::nm;
68 G4double R = std::sqrt(pos.x() * pos.x() + p << 69 G4double R = std::sqrt(pos.x()*pos.x()+pos.y()*pos.y()+pos.z()*pos.z())/CLHEP::nm;
69 // G4double ene = track->GetKineticEnergy( << 70 //G4double ene = track->GetKineticEnergy();
70 G4double trackE = track->GetKineticEnergy() << 71 G4double trackE = track->GetKineticEnergy()/CLHEP::eV;
71 if (RNP > R && track->GetTrackID() != 1) { << 72 if(RNP>R && track->GetTrackID() != 1 ){
72 if (track->GetDefinition()->GetPDGCharge() << 73 if(track->GetDefinition()->GetPDGCharge() != 0){
73 analysisManager->FillH1(2, trackE); << 74 analysisManager->FillH1(2,trackE);
74 } << 75 }else {
75 else { << 76 analysisManager->FillH1(3,trackE);
76 analysisManager->FillH1(3, trackE); << 77 }
77 } << 78 }
78 } << 79 if(RNP<R && track->GetTrackID() != 1 ){
79 if (RNP < R && track->GetTrackID() != 1) { << 80 if(track->GetDefinition()->GetPDGCharge() != 0){
80 if (track->GetDefinition()->GetPDGCharge() << 81 analysisManager->FillH1(6,R);
81 analysisManager->FillH1(6, R); << 82 analysisManager->FillH2(1,R,trackE);
82 analysisManager->FillH2(1, R, trackE); << 83 }else {
83 } << 84 analysisManager->FillH1(7,R);
84 else { << 85 analysisManager->FillH2(2,R,trackE);
85 analysisManager->FillH1(7, R); << 86 }
86 analysisManager->FillH2(2, R, trackE); << 87
87 } << 88 }
88 } <<
89 } 89 }
90 90
91 //....oooOO0OOooo........oooOO0OOooo........oo 91 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
92 92
93 void TrackingAction::PostUserTrackingAction(co << 93 void TrackingAction::PostUserTrackingAction(const G4Track* /*track*/)
>> 94 {
>> 95
>> 96 }
94 97
95 //....oooOO0OOooo........oooOO0OOooo........oo 98 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
96 99