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 biasing/B02/src/B02RunAction.cc 27 /// \brief Implementation of the B02RunAction class 28 // 29 // 30 // 31 #include "B02RunAction.hh" 32 33 #include "B02Run.hh" 34 35 //-- In order to obtain detector information. 36 #include "B02DetectorConstruction.hh" 37 #include "B02ImportanceDetectorConstruction.hh" 38 39 #include "G4RunManager.hh" 40 #include "G4THitsMap.hh" 41 #include "G4UnitsTable.hh" 42 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 44 // 45 // B02RunAction 46 // 47 // 48 // 49 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 50 51 // Constructor 52 B02RunAction::B02RunAction() 53 : G4UserRunAction(), 54 // fFieldName(15), 55 fFieldValue(14) 56 { 57 // - Prepare data member for B02Run. 58 // vector represents a list of MultiFunctionalDetector names. 59 fSDName.push_back(G4String("ConcreteSD")); 60 } 61 62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 63 64 // Destructor. 65 B02RunAction::~B02RunAction() 66 { 67 fSDName.clear(); 68 } 69 70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 71 72 G4Run* B02RunAction::GenerateRun() 73 { 74 // Generate new RUN object, which is specially 75 // dedicated for MultiFunctionalDetector scheme. 76 // Detail description can be found in B02Run.hh/cc. 77 return new B02Run(fSDName); 78 } 79 80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 81 82 void B02RunAction::BeginOfRunAction(const G4Run* aRun) 83 { 84 G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl; 85 } 86 87 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 88 89 void B02RunAction::EndOfRunAction(const G4Run* aRun) 90 { 91 G4cout << " ###### EndOfRunAction " << G4endl; 92 //- B02Run object. 93 B02Run* b02Run = (B02Run*)aRun; 94 //--- Dump all socred quantities involved in B02Run. 95 // re02Run->DumpAllScorer(); 96 //--- 97 G4RunManager* mgr = G4RunManager::GetRunManager(); 98 // 99 100 for (G4int i = 0; i < (G4int)fSDName.size(); i++) { 101 const G4VUserDetectorConstruction* vdet = mgr->GetUserDetectorConstruction(); 102 B02ImportanceDetectorConstruction* bdet = (B02ImportanceDetectorConstruction*)vdet; 103 // 104 105 //--------------------------------------------- 106 // Dump accumulated quantities for this RUN. 107 // (Display only central region of x-y plane) 108 // 0 ConcreteSD/Collisions 109 // 1 ConcreteSD/CollWeight 110 // 2 ConcreteSD/Population 111 // 3 ConcreteSD/TrackEnter 112 // 4 ConcreteSD/SL 113 // 5 ConcreteSD/SLW 114 // 6 ConcreteSD/SLWE 115 // 7 ConcreteSD/SLW_V 116 // 8 ConcreteSD/SLWE_V 117 //--------------------------------------------- 118 G4THitsMap<G4double>* Collisions = b02Run->GetHitsMap(fSDName[i] + "/Collisions"); 119 G4THitsMap<G4double>* CollWeight = b02Run->GetHitsMap(fSDName[i] + "/CollWeight"); 120 G4THitsMap<G4double>* Population = b02Run->GetHitsMap(fSDName[i] + "/Population"); 121 G4THitsMap<G4double>* TrackEnter = b02Run->GetHitsMap(fSDName[i] + "/TrackEnter"); 122 G4THitsMap<G4double>* SL = b02Run->GetHitsMap(fSDName[i] + "/SL"); 123 G4THitsMap<G4double>* SLW = b02Run->GetHitsMap(fSDName[i] + "/SLW"); 124 G4THitsMap<G4double>* SLWE = b02Run->GetHitsMap(fSDName[i] + "/SLWE"); 125 G4THitsMap<G4double>* SLW_V = b02Run->GetHitsMap(fSDName[i] + "/SLW_V"); 126 G4THitsMap<G4double>* SLWE_V = b02Run->GetHitsMap(fSDName[i] + "/SLWE_V"); 127 128 if (IsMaster()) { 129 G4cout << "\n--------------------End of Global Run-----------------------" << G4endl; 130 G4cout << " Number of event processed : " << aRun->GetNumberOfEvent() << G4endl; 131 } 132 else { 133 G4cout << "\n--------------------End of Local Run------------------------" << G4endl; 134 G4cout << " Number of event processed : " << aRun->GetNumberOfEvent() << G4endl; 135 } 136 137 G4cout << "=============================================================" << G4endl; 138 G4cout << "=============================================================" << G4endl; 139 140 std::ostream* myout = &G4cout; 141 PrintHeader(myout); 142 143 for (G4int iz = 0; iz < 20; iz++) { 144 G4double* SumCollisions = (*Collisions)[iz]; 145 G4double* SumCollWeight = (*CollWeight)[iz]; 146 G4double* Populations = (*Population)[iz]; 147 G4double* TrackEnters = (*TrackEnter)[iz]; 148 G4double* SLs = (*SL)[iz]; 149 G4double* SLWs = (*SLW)[iz]; 150 G4double* SLWEs = (*SLWE)[iz]; 151 G4double* SLW_Vs = (*SLW_V)[iz]; 152 G4double* SLWE_Vs = (*SLWE_V)[iz]; 153 if (!SumCollisions) SumCollisions = new G4double(0.0); 154 if (!SumCollWeight) SumCollWeight = new G4double(0.0); 155 if (!Populations) Populations = new G4double(0.0); 156 if (!TrackEnters) TrackEnters = new G4double(0.0); 157 if (!SLs) SLs = new G4double(0.0); 158 if (!SLWs) SLWs = new G4double(0.0); 159 if (!SLWEs) SLWEs = new G4double(0.0); 160 if (!SLW_Vs) SLW_Vs = new G4double(0.0); 161 if (!SLWE_Vs) SLWE_Vs = new G4double(0.0); 162 G4double NumWeightedEnergy = 0.0; 163 G4double FluxWeightedEnergy = 0.0; 164 G4double AverageTrackWeight = 0.0; 165 if (*SLW_Vs != 0.) NumWeightedEnergy = (*SLWE_Vs) / (*SLW_Vs); 166 if (*SLWs != 0.) FluxWeightedEnergy = (*SLWEs) / (*SLWs); 167 if (*SLs != 0.) AverageTrackWeight = (*SLWs) / (*SLs); 168 G4String cname = bdet->GetCellName(iz); 169 G4cout << std::setw(fFieldValue) << cname << " |" << std::setw(fFieldValue) << (*TrackEnters) 170 << " |" << std::setw(fFieldValue) << (*Populations) << " |" << std::setw(fFieldValue) 171 << (*SumCollisions) << " |" << std::setw(fFieldValue) << (*SumCollWeight) << " |" 172 << std::setw(fFieldValue) << NumWeightedEnergy << " |" << std::setw(fFieldValue) 173 << FluxWeightedEnergy << " |" << std::setw(fFieldValue) << AverageTrackWeight << " |" 174 << std::setw(fFieldValue) << (*SLs) << " |" << std::setw(fFieldValue) << (*SLWs) 175 << " |" << std::setw(fFieldValue) << (*SLW_Vs) << " |" << std::setw(fFieldValue) 176 << (*SLWEs) << " |" << std::setw(fFieldValue) << (*SLWE_Vs) << " |" << G4endl; 177 } 178 G4cout << "=============================================" << G4endl; 179 } 180 } 181 182 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 183 184 void B02RunAction::PrintHeader(std::ostream* out) 185 { 186 std::vector<G4String> vecScoreName; 187 vecScoreName.push_back("Tr.Entering"); 188 vecScoreName.push_back("Population"); 189 vecScoreName.push_back("Collisions"); 190 vecScoreName.push_back("Coll*WGT"); 191 vecScoreName.push_back("NumWGTedE"); 192 vecScoreName.push_back("FluxWGTedE"); 193 vecScoreName.push_back("Av.Tr.WGT"); 194 vecScoreName.push_back("SL"); 195 vecScoreName.push_back("SLW"); 196 vecScoreName.push_back("SLW_v"); 197 vecScoreName.push_back("SLWE"); 198 vecScoreName.push_back("SLWE_v"); 199 200 // head line 201 // std::string vname; 202 // vname = FillString("Volume", ' ', fFieldName+1); 203 //*out << vname << '|'; 204 *out << std::setw(fFieldValue) << "Volume" 205 << " |"; 206 for (std::vector<G4String>::iterator it = vecScoreName.begin(); it != vecScoreName.end(); it++) { 207 // vname = FillString((*it), 208 // ' ', 209 // fFieldValue+1, 210 // false); 211 // *out << vname << '|'; 212 *out << std::setw(fFieldValue) << (*it) << " |"; 213 } 214 *out << G4endl; 215 } 216 217 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 218 219 std::string B02RunAction::FillString(const std::string& name, char c, G4int n, G4bool back) 220 { 221 std::string fname(""); 222 G4int k = n - name.size(); 223 if (k > 0) { 224 if (back) { 225 fname = name; 226 fname += std::string(k, c); 227 } 228 else { 229 fname = std::string(k, c); 230 fname += name; 231 } 232 } 233 else { 234 fname = name; 235 } 236 return fname; 237 } 238 239 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 240