Geant4 Cross Reference |
1 // ******************************************* 1 // ******************************************************************** 2 // * License and Disclaimer 2 // * License and Disclaimer * 3 // * 3 // * * 4 // * The Geant4 software is copyright of th 4 // * The Geant4 software is copyright of the Copyright Holders of * 5 // * the Geant4 Collaboration. It is provided 5 // * the Geant4 Collaboration. It is provided under the terms and * 6 // * conditions of the Geant4 Software License 6 // * conditions of the Geant4 Software License, included in the file * 7 // * LICENSE and available at http://cern.ch/ 7 // * LICENSE and available at http://cern.ch/geant4/license . These * 8 // * include a list of copyright holders. 8 // * include a list of copyright holders. * 9 // * 9 // * * 10 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 11 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 12 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 13 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 14 // * use. Please see the license in the file 14 // * use. Please see the license in the file LICENSE and URL above * 15 // * for the full disclaimer and the limitatio 15 // * for the full disclaimer and the limitation of liability. * 16 // * 16 // * * 17 // * This code implementation is the result 17 // * This code implementation is the result of the scientific and * 18 // * technical work of the GEANT4 collaboratio 18 // * technical work of the GEANT4 collaboration. * 19 // * By using, copying, modifying or distri 19 // * By using, copying, modifying or distributing the software (or * 20 // * any work based on the software) you ag 20 // * any work based on the software) you agree to acknowledge its * 21 // * use in resulting scientific publicati 21 // * use in resulting scientific publications, and indicate your * 22 // * acceptance of all terms of the Geant4 Sof 22 // * acceptance of all terms of the Geant4 Software license. * 23 // ******************************************* 23 // ******************************************************************** 24 // 24 // 25 #include <memory> 25 #include <memory> 26 #include "G4DNAEventSet.hh" 26 #include "G4DNAEventSet.hh" 27 #include "globals.hh" 27 #include "globals.hh" 28 #include "G4UnitsTable.hh" 28 #include "G4UnitsTable.hh" 29 #include "G4DNAMolecularReactionTable.hh" 29 #include "G4DNAMolecularReactionTable.hh" 30 30 31 Event::Event(const G4double& time, const Index 31 Event::Event(const G4double& time, const Index& index, ReactionData* pReactionData) 32 : fTimeStep(time) 32 : fTimeStep(time) 33 , fIndex(index) 33 , fIndex(index) 34 , fData(std::pair<std::unique_ptr<JumpingDat 34 , fData(std::pair<std::unique_ptr<JumpingData>, ReactionData*>(nullptr, 35 35 pReactionData)) 36 {} 36 {} 37 37 38 Event::Event(const G4double& time, const Index 38 Event::Event(const G4double& time, const Index& index, 39 std::unique_ptr<JumpingData>&& ju 39 std::unique_ptr<JumpingData>&& jumping) 40 : fTimeStep(time) 40 : fTimeStep(time) 41 , fIndex(index) 41 , fIndex(index) 42 , fData(std::pair<std::unique_ptr<JumpingDat 42 , fData(std::pair<std::unique_ptr<JumpingData>, ReactionData*>( 43 std::move(jumping), nullptr)) 43 std::move(jumping), nullptr)) 44 {} 44 {} 45 45 46 Event::~Event() = default; 46 Event::~Event() = default; 47 47 48 void Event::PrintEvent() const 48 void Event::PrintEvent() const 49 { 49 { 50 G4cout << "****PrintEvent::TimeStep : " << G 50 G4cout << "****PrintEvent::TimeStep : " << G4BestUnit(fTimeStep, "Time") 51 << " index : " << fIndex << " action 51 << " index : " << fIndex << " action : "; 52 if(std::get<0>(fData) == nullptr) 52 if(std::get<0>(fData) == nullptr) 53 { 53 { 54 G4cout << std::get<1>(fData)->GetReactant1 54 G4cout << std::get<1>(fData)->GetReactant1()->GetName() << " + " 55 << std::get<1>(fData)->GetReactant2 55 << std::get<1>(fData)->GetReactant2()->GetName() << " -> " 56 << std::get<1>(fData)->GetProducts( 56 << std::get<1>(fData)->GetProducts()->size() << G4endl; 57 } 57 } 58 else 58 else 59 { 59 { 60 G4cout << std::get<0>(fData)->first->GetNa 60 G4cout << std::get<0>(fData)->first->GetName() << " jumping to " 61 << std::get<0>(fData)->second << G4 61 << std::get<0>(fData)->second << G4endl; 62 } 62 } 63 } 63 } 64 64 65 G4bool comparatorEventSet::operator()(std::uni 65 G4bool comparatorEventSet::operator()(std::unique_ptr<Event> const& rhs, 66 std::uni 66 std::unique_ptr<Event> const& lhs) const 67 { 67 { 68 return rhs->GetTime() < lhs->GetTime(); 68 return rhs->GetTime() < lhs->GetTime(); 69 } 69 } 70 70 71 G4DNAEventSet::G4DNAEventSet() 71 G4DNAEventSet::G4DNAEventSet() 72 : fEventSet(comparatorEventSet()) 72 : fEventSet(comparatorEventSet()) 73 {} 73 {} 74 74 75 void G4DNAEventSet::CreateEvent(const G4double 75 void G4DNAEventSet::CreateEvent(const G4double& time, const Index& index, 76 Event::Reactio 76 Event::ReactionData* pReactionData) 77 { 77 { 78 auto pEvent = std::make_unique<Event>(time, 78 auto pEvent = std::make_unique<Event>(time, index, pReactionData); 79 AddEvent(std::move(pEvent)); 79 AddEvent(std::move(pEvent)); 80 } 80 } 81 81 82 void G4DNAEventSet::CreateEvent(const G4double 82 void G4DNAEventSet::CreateEvent(const G4double& time, const Index& index, 83 std::unique_pt 83 std::unique_ptr<Event::JumpingData> jum) 84 { 84 { 85 auto pEvent = std::make_unique<Event>(time, 85 auto pEvent = std::make_unique<Event>(time, index, std::move(jum)); 86 AddEvent(std::move(pEvent)); 86 AddEvent(std::move(pEvent)); 87 } 87 } 88 88 89 void G4DNAEventSet::RemoveEventOfVoxel(const I 89 void G4DNAEventSet::RemoveEventOfVoxel(const Index& key) 90 { 90 { 91 auto it = fEventMap.find(key); 91 auto it = fEventMap.find(key); 92 if(it != fEventMap.end()) 92 if(it != fEventMap.end()) 93 { 93 { 94 fEventSet.erase(it->second); 94 fEventSet.erase(it->second); 95 fEventMap.erase(it); 95 fEventMap.erase(it); 96 } 96 } 97 } 97 } 98 98 99 void G4DNAEventSet::RemoveEvent(EventSet::iter 99 void G4DNAEventSet::RemoveEvent(EventSet::iterator iter) 100 { 100 { 101 auto index = (*iter)->GetIndex(); 101 auto index = (*iter)->GetIndex(); 102 RemoveEventOfVoxel(index); 102 RemoveEventOfVoxel(index); 103 } 103 } 104 104 105 void G4DNAEventSet::AddEvent(std::unique_ptr<E 105 void G4DNAEventSet::AddEvent(std::unique_ptr<Event> pEvent) 106 { 106 { 107 // idea is no 2 events in one key (or index) 107 // idea is no 2 events in one key (or index) 108 auto key = pEvent->GetIndex(); 108 auto key = pEvent->GetIndex(); 109 RemoveEventOfVoxel(key); 109 RemoveEventOfVoxel(key); 110 auto it = fEventSet.emplace(std::move 110 auto it = fEventSet.emplace(std::move(pEvent)); 111 fEventMap[key] = std::get<0>(it); 111 fEventMap[key] = std::get<0>(it); 112 } 112 } 113 113 114 //____________________________________________ 114 //_____________________________________________________________________________________ 115 115 116 G4DNAEventSet::~G4DNAEventSet() { RemoveEventS 116 G4DNAEventSet::~G4DNAEventSet() { RemoveEventSet(); } 117 117 118 [[maybe_unused]] void G4DNAEventSet::PrintEven 118 [[maybe_unused]] void G4DNAEventSet::PrintEventSet() 119 { 119 { 120 G4cout<<G4endl; 120 G4cout<<G4endl; 121 G4cout << "********************************* 121 G4cout << "*****************************************************" << G4endl; 122 G4cout << "G4DNAEventSet::PrintEventSet() of 122 G4cout << "G4DNAEventSet::PrintEventSet() of : "<< this->size()<<" events "<< G4endl; 123 for(const auto& it : fEventSet) 123 for(const auto& it : fEventSet) 124 { 124 { 125 (*it).PrintEvent(); 125 (*it).PrintEvent(); 126 } 126 } 127 G4cout << "End PrintEventSet()" << G4endl; 127 G4cout << "End PrintEventSet()" << G4endl; 128 G4cout << "********************************* 128 G4cout << "*****************************************************" << G4endl; 129 G4cout << G4endl; 129 G4cout << G4endl; 130 } 130 }