| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 //
28 /// \file ScanDamage.hh
29 /// \brief Definition of the ScanDamage class
30
31 #ifndef ScanDamage_h
32 #define ScanDamage_h
33
34 #include <map>
35 #include <vector>
36 #include <string>
37 #include <tuple>
38 #include <set>
39 #include "Damage.hh"
40 #include <filesystem>
41 namespace fs = std::filesystem;
42 class TFile;
43 using ullint = unsigned long long int;
44 //....oooOO0OOooo........oooOO0OOooo........oo
45
46 struct VoxelData
47 {
48 VoxelData(int chromo, int domain, ullint f
49 {
50 fChromosome = chromo;
51 fDomain = domain;
52 fFirstBpCopyNum = firstBpCN;
53 }
54
55 ~VoxelData() {}
56
57 int fChromosome{0};
58 int fDomain{0};
59 ullint fFirstBpCopyNum{0};
60 };
61
62 //....oooOO0OOooo........oooOO0OOooo........oo
63
64 typedef std::vector<std::vector<ullint> > Tabl
65
66 //....oooOO0OOooo........oooOO0OOooo........oo
67
68 class Trier {
69 public:
70 bool operator()(const std::vector<ullint>&
71 {
72 bool bb = false;
73 if(a[0] < b[0]) bb = true;
74 return bb;
75 }
76 };
77
78 //....oooOO0OOooo........oooOO0OOooo........oo
79
80 class ScanDamage
81 {
82 public:
83 ScanDamage();
84 ~ScanDamage() = default;
85 std::map<unsigned int,std::map<unsigned in
86 void SetThresholdEnergy(double e) {fThresh
87 void SetProbabilityForIndirectSBSelection(
88 double GetThresholdEnergy() {return fThres
89 double GetProbabilityForIndirectSBSelectio
90 std::map<int, Table> GetMergedSBData() {re
91 double GetEdepSumInNucleus() {return fEdep
92 double GetTotalNbBpPlacedInGeo() {return f
93 double GetTotalNbHistonePlacedInGeo() {ret
94 double GetNucleusVolume() {return fNucleus
95 double GetNucleusMassDensity() {return fNu
96 double GetNucleusMass() {return fNucleusMa
97 std::map<int,ullint> GetChromosomeBpSizesM
98 void SkipScanningIndirectDamage() {fSkipSc
99 bool SkippedScanningIndirectDamage() {retu
100 private:
101 void ScanDamageFromPhys();
102 void ScanDamageFromChem();
103 void RetrieveVoxelBp();
104 void FillVoxelData();
105 void AnaPhysRootFile(const std::string fil
106 void AnaChemRootFile(fs::directory_entry e
107 void AnaPhysRootTree1(TFile*);
108 void AnaPhysRootTree2(TFile*);
109 void SortPhysTableWithSelection();
110 void SortChemTableWithSelection();
111 void ReadCellandVoxelDefFilePaths();
112 void MergeDamageFromPhysChem();
113 std::tuple<unsigned int, unsigned int> Get
114 double fThresholdEnergy{17.5};//eV
115 std::string fCellDefFilePath{""};
116 std::set<std::string> fVoxelDefFilesList;
117 std::map<std::string, int> fBpPerVoxel;
118 std::vector<VoxelData> fVoxels;
119 std::map<int, Table> fphysTables, fphysSle
120 std::map<unsigned int,std::map<unsigned in
121 double fEdepSumInNucleus{0}; //eV
122 int corruptedFiles = 0;
123 double fTotalNbBpPlacedInGeo{0};
124 double fTotalNbHistonePlacedInGeo{0};
125 double fNucleusVolume{0};
126 double fNucleusMassDensity{0};
127 double fNucleusMass{0};
128 double fProbabilityForIndirectSB{0.4};
129 std::map<int,ullint> fChromosomeBpMap; //S
130 bool fSkipScanningIndirectDamage{false};
131 };
132
133 //....oooOO0OOooo........oooOO0OOooo........oo
134
135 #endif
136
137
138