| 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 /// \file biasing/B02/src/B02Run.cc
27 /// \brief Implementation of the B02Run class
28 //
29 //
30 //
31
32 //....oooOO0OOooo........oooOO0OOooo........oo
33 //
34 // (Description)
35 // B02Run Class is for accumulating scored
36 // scored using G4MutiFunctionalDetector and
37 // Accumulation is done using G4THitsMap obje
38 //
39 // The constructor B02Run(const std::vector
40 // needs a vector filled with MultiFunctional
41 // was assigned at instantiation of MultiFunc
42 // Then B02Run constructor automatically scan
43 // in the MFD, and obtains collectionIDs of a
44 // to those primitive scorers. Futhermore, th
45 // for accumulating during a RUN are automati
46 // (*) Collection Name is same as primitive s
47 //
48 // The resultant information is kept inside
49 // data members.
50 // std::vector<G4String> fCollName;
51 // std::vector<G4int> fCollID;
52 // std::vector<G4THitsMap<G4double>*> fRunMap
53 //
54 // The resualtant HitsMap objects are obtain
55 // GetHitsMap(..).
56 //
57 //....oooOO0OOooo........oooOO0OOooo........oo
58
59 #include "B02Run.hh"
60
61 #include "G4MultiFunctionalDetector.hh"
62 #include "G4SDManager.hh"
63 #include "G4VPrimitiveScorer.hh"
64
65 //
66 //....oooOO0OOooo........oooOO0OOooo........oo
67 // Constructor.
68 // (The vector of MultiFunctionalDetector na
69 B02Run::B02Run(const std::vector<G4String> mfd
70 {
71 G4SDManager* SDman = G4SDManager::GetSDMpoin
72 //==========================================
73 // Initalize RunMaps for accumulation.
74 // Get CollectionIDs for HitCollections.
75 //==========================================
76 G4int Nmfd = mfdName.size();
77 for (G4int idet = 0; idet < Nmfd; idet++) {
78 G4String detName = mfdName[idet];
79 //--- Seek and Obtain MFD objects from SDm
80 G4MultiFunctionalDetector* mfd =
81 (G4MultiFunctionalDetector*)(SDman->Find
82 //
83 if (mfd) {
84 //--- Loop over the registered primitive
85 for (G4int icol = 0; icol < mfd->GetNumb
86 // Get Primitive Scorer object.
87 G4VPrimitiveScorer* scorer = mfd->GetP
88 // collection name and collectionID fo
89 // where type of HitsCollection is G4T
90 // The collection name is given by <MF
91 G4String collectionName = scorer->GetN
92 G4String fullCollectionName = detName
93 G4int collectionID = SDman->GetCollect
94
95 if (collectionID >= 0) {
96 G4cout << "++ " << fullCollectionNam
97 // Store obtained HitsCollection inf
98 // And, creates new G4THitsMap for a
99 fCollName.push_back(fullCollectionNa
100 fCollID.push_back(collectionID);
101 fRunMap.push_back(new G4THitsMap<G4d
102 }
103 else {
104 G4cout << "** collection " << fullCo
105 }
106 }
107 }
108 }
109 }
110
111 //....oooOO0OOooo........oooOO0OOooo........oo
112
113 // Destructor
114 // clear all data members.
115 B02Run::~B02Run()
116 {
117 //--- Clear HitsMap for RUN
118 G4int Nmap = fRunMap.size();
119 for (G4int i = 0; i < Nmap; i++) {
120 if (fRunMap[i]) fRunMap[i]->clear();
121 }
122 fCollName.clear();
123 fCollID.clear();
124 fRunMap.clear();
125 }
126
127 //....oooOO0OOooo........oooOO0OOooo........oo
128
129 // RecordEvent is called at end of event.
130 // For scoring purpose, the resultant quantit
131 // is accumulated during a Run.
132 void B02Run::RecordEvent(const G4Event* aEvent
133 {
134 numberOfEvent++; // This is an original lin
135
136 //=============================
137 // HitsCollection of This Event
138 //============================
139 G4HCofThisEvent* HCE = aEvent->GetHCofThisEv
140 if (!HCE) return;
141
142 //==========================================
143 // Sum up HitsMap of this Event into HitsMa
144 //==========================================
145 G4int Ncol = fCollID.size();
146 for (G4int i = 0; i < Ncol; i++) { // Loop
147 G4THitsMap<G4double>* EvtMap = 0;
148 if (fCollID[i] >= 0) { // Collection is a
149 EvtMap = (G4THitsMap<G4double>*)(HCE->Ge
150 }
151 else {
152 G4cout << " Error EvtMap Not Found " <<
153 }
154 if (EvtMap) {
155 //=== Sum up HitsMap of this event to Hi
156 *fRunMap[i] += *EvtMap;
157 //======================================
158 }
159 }
160 }
161
162 //....oooOO0OOooo........oooOO0OOooo........oo
163 //
164 // Access method for HitsMap of the RUN
165 //
166 //-----
167 // Access HitsMap.
168 // By MultiFunctionalDetector name and Colle
169 G4THitsMap<G4double>* B02Run::GetHitsMap(const
170 {
171 G4String fullName = detName + "/" + colName;
172 return GetHitsMap(fullName);
173 }
174
175 //-----
176 // Access HitsMap.
177 // By full description of collection name, th
178 // <MultiFunctional Detector Name>/<Primiti
179 G4THitsMap<G4double>* B02Run::GetHitsMap(const
180 {
181 G4int Ncol = fCollName.size();
182 for (G4int i = 0; i < Ncol; i++) {
183 if (fCollName[i] == fullName) {
184 return fRunMap[i];
185 }
186 }
187 return NULL;
188 }
189
190 //....oooOO0OOooo........oooOO0OOooo........oo
191
192 // - Dump All HitsMap of this RUN. (for debugi
193 // This method calls G4THisMap::PrintAll() f
194 void B02Run::DumpAllScorer()
195 {
196 // - Number of HitsMap in this RUN.
197 G4int n = GetNumberOfHitsMap();
198 // - GetHitsMap and dump values.
199 for (G4int i = 0; i < n; i++) {
200 G4THitsMap<G4double>* RunMap = GetHitsMap(
201 if (RunMap) {
202 G4cout << " PrimitiveScorer RUN " << Run
203 << G4endl;
204 G4cout << " Number of entries " << RunMa
205 std::map<G4int, G4double*>::iterator itr
206 for (; itr != RunMap->GetMap()->end(); i
207 G4cout << " copy no.: " << itr->first
208 }
209 }
210 }
211 }
212
213 //....oooOO0OOooo........oooOO0OOooo........oo
214
215 void B02Run::Merge(const G4Run* aRun)
216 {
217 const B02Run* localRun = static_cast<const B
218 //==========================================
219 // Merge HitsMap of working threads
220 //==========================================
221 G4int nCol = localRun->fCollID.size();
222 for (G4int i = 0; i < nCol; i++) { // Loop
223 if (localRun->fCollID[i] >= 0) {
224 *fRunMap[i] += *localRun->fRunMap[i];
225 }
226 }
227
228 G4Run::Merge(aRun);
229 }
230
231 //....oooOO0OOooo........oooOO0OOooo........oo
232