| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 /// \file electromagnetic/TestEm11/include/doiPETRun.hh
27 /// \brief Definition of the doiPETRun class
28 //
29 // $Id: doiPETRun.hh 71375 2013-06-14 07:39:33Z maire $
30 //
31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33
34 #ifndef doiPETRun_h
35 #define doiPETRun_h 1
36
37 #include "G4Run.hh"
38 #include "G4VProcess.hh"
39 #include "globals.hh"
40
41 #include <iostream>
42 #include <fstream>
43 #include <sstream>
44 #include <map>
45 #include <set>
46 class doiPETDetectorConstruction;
47 class G4ParticleDefinition;
48 class doiPETRunAction;
49
50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
51 class InteractionInformation;
52 class doiPETAnalysis;
53 class doiPETRun : public G4Run
54 {
55 public:
56 doiPETRun(/*doiPETDetectorConstruction**/);
57 ~doiPETRun();
58
59 public:
60 //
61 void GetIntractionInfomation(InteractionInformation*);
62 void FindInteractingCrystal();
63 void Clear();
64 void OpenRun(G4String);
65 void GetEventIDRun(G4int);
66 void CalulateAcquisitionTime();
67 void SetAnnihilationTime(G4double);
68 void SetEventID(G4int);
69 virtual void Merge(const G4Run*);
70 G4ThreeVector CenterOfMassInteractionPos(const std::vector<G4int>&, const std::vector<G4double>&, G4double, const std::vector<G4ThreeVector>&);
71 //void EndOfRun();
72
73 private:
74
75 std::multimap< G4int, InteractionInformation* > mapBlockInteraction;
76 std::set<G4int> setBlockInteraction;
77 std::ofstream ofs;
78
79 std::vector<G4int> crystalID_vec;
80 std::vector<G4ThreeVector> posInter_vec;
81 std::vector<G4double>edepInCry_vec;
82
83 //center of mass interaction of position
84 G4ThreeVector posInterInCrystal;
85
86 G4int eventID;
87 doiPETRunAction* fRunAction;
88 doiPETAnalysis* fAnalysis;
89 G4double totalEdep;
90 G4int blockID, crystalID;
91
92 G4double totalTime;
93 G4double prev_totalTime;
94 G4double timeInterval;
95
96 G4ThreeVector interactionPos;
97 //
98 G4double interactionTime;
99 G4double timeOfAnnihil;
100
101 G4int numberofInteractions;
102 G4double edep;
103 G4double edepMax;
104 G4double edep_AfterCrystalBlurring;
105
106 };
107
108
109 //The following is to get interaction information
110 class InteractionInformation
111 {
112
113 public:
114 InteractionInformation(){;};
115 ~InteractionInformation(){;};
116
117 //set energy deposition
118 void SetEdep(G4double e) {edep=e;};
119
120 //set block number
121 void SetBlockNo(G4int n) {blockNo=n;};
122
123 //set crystal ID (this is continuous numbering of crystals)
124 void SetCrystalNo(G4int n) {crystalNo=n;};
125
126 //set global time
127 void SetGlobalTime(G4double t){globalTime=t;};
128
129 //set interaction position with respect to the crystal axis (local position)
130 void SetInteractionPositionInCrystal(G4ThreeVector pos){crystalPosition = pos;};
131
132
133 G4double GetEdep() {return edep;};
134 G4int GetBlockNo() {return blockNo;};
135 G4int GetCrystalNo() {return crystalNo;};
136 G4double GetGlobalTime() {return globalTime;};
137
138 //Interaction position in the detector
139 G4ThreeVector GetInteractionPositionInCrystal(){return crystalPosition;};
140
141
142 private:
143 G4double edep;
144 G4double globalTime;
145 G4int blockNo;
146 G4int crystalNo;
147 G4ThreeVector crystalPosition;
148
149 };
150
151 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
152
153 #endif
154
155