| 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 // File: CCaloOrganization.cc
28 // Description: Packing, unpacking and other r
29 // calorimetric numbering schema
30 //////////////////////////////////////////////
31 #include "CCaloOrganization.hh"
32
33 //#define debug
34
35 unsigned int CCaloOrganization::packindex(G4in
36 G4in
37 //So this is the actual encoding of the inde
38 //top 4 bits encode Detector type
39 //
40 // Should work for all calorimeter with no d
41
42 unsigned int idx=(det&15)<<28; //bits 28-3
43 idx+=(((z+1)/2)&1)<<20; //bits 20
44 idx+=(eta&1023)<<10; //bits 10-1
45 idx+=(phi&1023); //bits 0-9
46 #ifdef debug
47 G4cout << " ECAL packing " << det << " " <<
48 << " into " << idx << G4endl;
49 #endif
50 return idx;
51 }
52
53 unsigned int CCaloOrganization::packindex(G4in
54 G4in
55 //So this is the actual encoding of the inde
56 //top 4 bits encode Detector type
57 //next 4 bits encode depth information
58 // Should work for all calorimeter with no d
59
60 unsigned int idx=(det&15)<<28; //bits 28-31
61 idx+=(depth&15)<<24; //bits 24-27
62 idx+=(z&1)<<20; //bits 20
63 idx+=(eta&1023)<<10; //bits 10-19
64 idx+=(phi&1023); //bits 0-9
65 #ifdef debug
66 G4cout << " HCAL packing " << det << " " <<
67 << " " << phi << " into " << idx <<
68 #endif
69 return idx;
70 }
71
72
73 void CCaloOrganization::unpackindex(const unsi
74 G4int& eta
75 det = (idx>>28)&15;
76 z = (idx>>20)&1;
77 z = 2*z-1;
78 eta = (idx>>10)&1023;
79 phi = (idx&1023);
80 }
81
82
83 void CCaloOrganization::unpackindex(const unsi
84 G4int& dep
85 G4int& phi
86 det = (idx>>28)&15;
87 depth=(idx>>24)&15;
88 z = (idx>>20)&1;
89 eta = (idx>>10)&1023;
90 phi = (idx&1023);
91 }
92
93
94 G4int CCaloOrganization::getUnitWithMaxEnergy(
95
96 //look for max
97 G4int UnitWithMaxEnergy = 0;
98 G4float maxEnergy = 0.;
99
100 for(std::map<G4int,G4float,std::less<G4int>
101 iter != themap.end(); iter++){
102
103 if( maxEnergy < (*iter).second) {
104 maxEnergy = (*iter).second;
105 UnitWithMaxEnergy = (*iter).first;
106 }
107 }
108 #ifdef debug
109 G4cout << " *** max energy of " << maxEnergy
110 << UnitWithMaxEnergy;
111 G4int det,z,eta,phi;
112 unpackindex(UnitWithMaxEnergy, det, z, eta,
113 G4cout << " corresponding to z= " << z << "
114 << G4endl;
115 #endif
116 return UnitWithMaxEnergy;
117
118 }
119
120
121 G4float CCaloOrganization::energyInMatrix(G4in
122 G4int
123 std::m
124
125 G4int det,z,eta,phi;
126 this->unpackindex(crystalWithMaxEnergy, det,
127 #ifdef debug
128 G4int ncristals=0;
129 #endif
130 G4int goBackInEta = nCellInEta/2;
131 G4int goBackInPhi = nCellInPhi/2;
132 G4int startEta = eta-goBackInEta;
133 G4int startPhi = phi-goBackInPhi;
134
135 G4float totalEnergy = 0.;
136
137 for(G4int ieta=startEta; ieta<startEta+nCell
138 for(G4int iphi=startPhi; iphi<startPhi+nCe
139
140 G4int index = this->packindex(det,z,ieta
141 totalEnergy += themap[index];
142 #ifdef debug
143 ncristals+=1;
144 G4cout << "ieta - iphi - E = " << ieta <
145 << themap[index] << G4endl;
146 #endif
147 }
148 }
149
150 #ifdef debug
151 G4cout << "Energy in " << nCellInEta << " ce
152 << nCellInPhi << " cells in phi matri
153 << " for " << ncristals << " crystals
154 #endif
155 return totalEnergy;
156
157 }
158