| 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 // 26 //
>> 27 // $Id: G4MultiFunctionalDetector.cc,v 1.4 2006/06/29 18:05:47 gunter Exp $
>> 28 // GEANT4 tag $Name: geant4-08-03-patch-02 $
27 // 29 //
28 // G4MultiFunctionalDetector 30 // G4MultiFunctionalDetector
29 // <<
30 // 2010-07-23 T.Aso Call PS if the step length <<
31 // <<
32 // <<
33 #include "G4MultiFunctionalDetector.hh" 31 #include "G4MultiFunctionalDetector.hh"
34 <<
35 #include "G4SDManager.hh" 32 #include "G4SDManager.hh"
36 #include "G4VPrimitiveScorer.hh" 33 #include "G4VPrimitiveScorer.hh"
37 34
38 G4MultiFunctionalDetector::G4MultiFunctionalDe << 35 G4MultiFunctionalDetector::G4MultiFunctionalDetector(G4String name)
39 : G4VSensitiveDetector(name) << 36 :G4VSensitiveDetector(name)
40 { << 37 {;}
41 ; <<
42 } <<
43 38
44 G4MultiFunctionalDetector::~G4MultiFunctionalD 39 G4MultiFunctionalDetector::~G4MultiFunctionalDetector()
45 { << 40 {;}
46 for (auto pr : primitives) <<
47 delete pr; <<
48 primitives.clear(); <<
49 } <<
50 41
51 G4bool G4MultiFunctionalDetector::ProcessHits( << 42 G4bool G4MultiFunctionalDetector::ProcessHits(G4Step* aStep,G4TouchableHistory* aTH)
52 { << 43 {
53 if (aStep->GetStepLength() > 0. || aStep->Ge << 44 G4int nPrim = primitives.size();
54 for (auto pr : primitives) { << 45 for(G4int iPrim=0;iPrim<nPrim;iPrim++)
55 pr->HitPrimitive(aStep, aTH); << 46 { primitives[iPrim]->HitPrimitive(aStep,aTH); }
56 } << 47 return true;
57 } <<
58 return true; <<
59 } 48 }
60 49
61 G4bool G4MultiFunctionalDetector::RegisterPrim 50 G4bool G4MultiFunctionalDetector::RegisterPrimitive(G4VPrimitiveScorer* aPS)
62 { 51 {
63 for (auto pr : primitives) { << 52 G4int nPrim = primitives.size();
64 if (pr == aPS) { << 53 for(G4int iPrim=0;iPrim<nPrim;iPrim++)
65 G4ExceptionDescription ED; << 54 {
66 ED << "Primitive <" << aPS->GetName() << << 55 if(primitives[iPrim]==aPS)
67 << ">." << G4endl << "Method Register << 56 {
68 G4Exception("G4MultiFunctionalDetector:: << 57 G4cerr << "Primitive <" << aPS->GetName() << "> is already defined in <" << SensitiveDetectorName
69 return false; << 58 << ">." << G4endl << "Method RegisterPrimitive() is ignored." << G4endl;
70 } << 59 return false;
71 } << 60 }
72 primitives.push_back(aPS); << 61 }
73 aPS->SetMultiFunctionalDetector(this); << 62 primitives.push_back(aPS);
74 collectionName.insert(aPS->GetName()); << 63 aPS->SetMultiFunctionalDetector(this);
75 if (G4SDManager::GetSDMpointer()->FindSensit << 64 collectionName.insert(aPS->GetName());
76 { << 65 if(G4SDManager::GetSDMpointer()->FindSensitiveDetector(SensitiveDetectorName,false))
77 // This G4MultiFunctionalDetector has alre << 66 {
78 // G4SDManager. Make sure this new primiti << 67 // This G4MultiFunctionalDetector has already been registered to G4SDManager.
79 G4SDManager::GetSDMpointer()->AddNewCollec << 68 // Make sure this new primitive is registered as well.
80 } << 69 G4SDManager::GetSDMpointer()->AddNewCollection(SensitiveDetectorName,aPS->GetName());
81 return true; << 70 }
>> 71 return true;
82 } 72 }
83 73
84 G4bool G4MultiFunctionalDetector::RemovePrimit 74 G4bool G4MultiFunctionalDetector::RemovePrimitive(G4VPrimitiveScorer* aPS)
85 { 75 {
86 auto pr = std::find(primitives.begin(), prim << 76 std::vector<G4VPrimitiveScorer*>::iterator iterPS;
87 if (pr != primitives.end()) { << 77 std::vector<G4String>::iterator iterName = collectionName.begin();
88 primitives.erase(pr); << 78 for(iterPS=primitives.begin();iterPS!=primitives.end();iterPS++)
89 aPS->SetMultiFunctionalDetector(nullptr); << 79 {
90 return true; << 80 if(*iterPS==aPS)
91 } << 81 {
92 G4cerr << "Primitive <" << aPS->GetName() << << 82 primitives.erase(iterPS);
93 << ">." << G4endl << "Method RemovePr << 83 aPS->SetMultiFunctionalDetector(0);
94 return false; << 84 return true;
95 } << 85 }
>> 86 iterName++;
>> 87 }
>> 88 G4cerr << "Primitive <" << aPS->GetName() << "> is not defined in <" << SensitiveDetectorName
>> 89 << ">." << G4endl << "Method RemovePrimitive() is ignored." << G4endl;
>> 90 return false;
>> 91 }
96 92
97 void G4MultiFunctionalDetector::Initialize(G4H 93 void G4MultiFunctionalDetector::Initialize(G4HCofThisEvent* HC)
98 { 94 {
99 for (auto pr : primitives) { << 95 G4int nPrim = primitives.size();
100 pr->Initialize(HC); << 96 for(G4int iPrim=0;iPrim<nPrim;iPrim++)
101 } << 97 { primitives[iPrim]->Initialize(HC); }
102 } 98 }
103 99
104 void G4MultiFunctionalDetector::EndOfEvent(G4H 100 void G4MultiFunctionalDetector::EndOfEvent(G4HCofThisEvent* HC)
105 { 101 {
106 for (auto pr : primitives) { << 102 G4int nPrim = primitives.size();
107 pr->EndOfEvent(HC); << 103 for(G4int iPrim=0;iPrim<nPrim;iPrim++)
108 } << 104 { primitives[iPrim]->EndOfEvent(HC); }
109 } 105 }
110 106
111 void G4MultiFunctionalDetector::clear() 107 void G4MultiFunctionalDetector::clear()
112 { 108 {
113 for (auto pr : primitives) { << 109 G4int nPrim = primitives.size();
114 pr->clear(); << 110 for(G4int iPrim=0;iPrim<nPrim;iPrim++)
115 } << 111 { primitives[iPrim]->clear(); }
116 } 112 }
117 113
118 void G4MultiFunctionalDetector::DrawAll() 114 void G4MultiFunctionalDetector::DrawAll()
119 { 115 {
120 for (auto pr : primitives) { << 116 G4int nPrim = primitives.size();
121 pr->DrawAll(); << 117 for(G4int iPrim=0;iPrim<nPrim;iPrim++)
122 } << 118 { primitives[iPrim]->DrawAll(); }
123 } 119 }
124 120
125 void G4MultiFunctionalDetector::PrintAll() 121 void G4MultiFunctionalDetector::PrintAll()
126 { 122 {
127 for (auto pr : primitives) { << 123 G4int nPrim = primitives.size();
128 pr->PrintAll(); << 124 for(G4int iPrim=0;iPrim<nPrim;iPrim++)
129 } << 125 { primitives[iPrim]->PrintAll(); }
130 } 126 }
>> 127
131 128