Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/dna/utils/src/G4DNAEventSet.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

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