| 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 // This example is provided by the Geant4-DNA collaboration
27 // dnadamage3 example is derived from the chem6 example
28 // chem6 example authors: W. G. Shin and S. Incerti (CENBG, France)
29 //
30 // Any report or published results obtained using the Geant4-DNA software
31 // shall cite the following Geant4-DNA collaboration publication:
32 // Med. Phys. 37 (2010) 4692-4708
33 // J. Comput. Phys. 274 (2014) 841-882
34 // The Geant4-DNA web site is available at http://geant4-dna.org
35 //
36 // Authors: J. Naoki D. Kondo (UCSF, US)
37 // J. Ramos-Mendez and B. Faddegon (UCSF, US)
38 //
39 /// \file ScoreStrandBreaks.hh
40 /// \brief Definition of the ScoreStrandBreaks class
41
42 #ifndef DNADAMAGE2_ScoreSB_h
43 #define DNADAMAGE2_ScoreSB_h 1
44
45 #include "DetectorConstruction.hh"
46 #include "MoleculeInserter.hh"
47
48 #include "G4DNAIRT.hh"
49 #include "G4THitsMap.hh"
50 #include "G4UIcmdWithADoubleAndUnit.hh"
51 #include "G4UIcmdWithAString.hh"
52 #include "G4UIcmdWithAnInteger.hh"
53 #include "G4UImessenger.hh"
54 #include "G4VPrimitiveScorer.hh"
55
56 class G4VAnalysisManager;
57 class G4MolecularConfiguration;
58
59 typedef std::map<G4int, std::map<G4int, std::map<G4int, std::map<G4int, G4double>>>>
60 EnergyDepositMap;
61 typedef std::map<G4int, std::vector<std::vector<G4int>>> DamageMap;
62
63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
64
65 class ScoreStrandBreaks : public G4VPrimitiveScorer, public G4UImessenger
66 {
67 public:
68 ScoreStrandBreaks(G4String name, DetectorConstruction*, G4double*);
69 ~ScoreStrandBreaks() override;
70
71 void WriteWithAnalysisManager(G4VAnalysisManager*);
72 DamageMap GetDamageMap() { return fDirectDamageMap; }
73 void ClearDamageMap() { fDirectDamageMap.clear(); }
74
75 void Initialize(G4HCofThisEvent*) override;
76 void EndOfEvent(G4HCofThisEvent*) override;
77 void Clear();
78 void DrawAll() override;
79 void PrintAll() override;
80 void AbsorbResultsFromWorkerScorer(G4VPrimitiveScorer*);
81 void OutputAndClear(G4double, G4double);
82 void SetNewValue(G4UIcommand*, G4String) override;
83
84 void ASCII(G4double, G4double);
85
86 private:
87 G4int fnbOfEvents = 0;
88 G4bool fDNAInserted = false;
89 G4double fBreakEnergy = 17.5 * eV;
90 G4double fEnergyDeposit = 0;
91 G4double* fRadius = nullptr;
92
93 G4String fOutputName = "DirectDamageInfo";
94 G4String fOutputType = "ascii";
95 DamageMap fDirectDamageMap;
96 DamageMap fIndirectDamageMap;
97 EnergyDepositMap fEnergyDepositMap;
98
99 G4UIcmdWithAString* fpOutputFileUI = nullptr;
100 G4UIcmdWithAString* fpOutputTypeUI = nullptr;
101 G4UIcmdWithADoubleAndUnit* fpBreakEnergyUI = nullptr;
102
103 DetectorConstruction* fpDetector = nullptr;
104 std::map<G4int, G4double> fDoseArray;
105
106 MoleculeInserter* fpGun = nullptr;
107
108 protected:
109 G4bool ProcessHits(G4Step*, G4TouchableHistory*) override;
110 };
111
112 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
113
114 #endif
115