| 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 /// \file eventgenerator/exgps/src/TrackingAct
28 /// \brief Implementation of the TrackingActio
29 //
30 //....oooOO0OOooo........oooOO0OOooo........oo
31 //....oooOO0OOooo........oooOO0OOooo........oo
32
33 #include "TrackingAction.hh"
34
35 #include "HistoManager.hh"
36
37 #include "G4PhysicalConstants.hh"
38 #include "G4SystemOfUnits.hh"
39 #include "G4Track.hh"
40
41 //....oooOO0OOooo........oooOO0OOooo........oo
42
43 TrackingAction::TrackingAction() : G4UserTrack
44
45 //....oooOO0OOooo........oooOO0OOooo........oo
46
47 void TrackingAction::PreUserTrackingAction(con
48 {
49 G4int pid = track->GetDynamicParticle()->Get
50 G4double ekin = track->GetKineticEnergy();
51 G4ThreeVector vertex = track->GetPosition();
52 G4ThreeVector direction = track->GetMomentum
53 G4double weight = track->GetWeight();
54
55 G4double x = vertex.x(), y = vertex.y(), z =
56 G4double theta = direction.theta(), phi = di
57 if (phi < 0.) phi += twopi;
58 G4double cost = std::cos(theta);
59
60 G4AnalysisManager* analysisManager = G4Analy
61
62 G4double r = vertex.mag();
63 G4double dr = analysisManager->GetH1Width(2)
64 G4double dv = 2 * twopi * r * r * dr;
65
66 // fill histograms
67 analysisManager->FillH1(1, ekin);
68 if (dv > 0.) analysisManager->FillH1(2, r, 1
69 analysisManager->FillH1(3, cost);
70 analysisManager->FillH1(4, phi);
71 analysisManager->FillH2(1, x, y);
72 analysisManager->FillH2(2, y, z);
73 analysisManager->FillH2(3, z, x);
74 analysisManager->FillH2(4, phi, cost);
75 analysisManager->FillH2(5, phi, theta);
76
77 // fill ntuple
78 analysisManager->FillNtupleIColumn(0, pid);
79 analysisManager->FillNtupleDColumn(1, ekin);
80 analysisManager->FillNtupleDColumn(2, x);
81 analysisManager->FillNtupleDColumn(3, y);
82 analysisManager->FillNtupleDColumn(4, z);
83 analysisManager->FillNtupleDColumn(5, theta)
84 analysisManager->FillNtupleDColumn(6, phi);
85 analysisManager->FillNtupleDColumn(7, weight
86 analysisManager->AddNtupleRow();
87 }
88
89 //....oooOO0OOooo........oooOO0OOooo........oo
90