| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 ////////////////////////////////////////////// 26 ///////////////////////////////////////////////////////////////////////////////
27 // File: CCaloOrganization.cc 27 // File: CCaloOrganization.cc
28 // Description: Packing, unpacking and other r 28 // Description: Packing, unpacking and other related utilities for
29 // calorimetric numbering schema 29 // calorimetric numbering schema
30 ////////////////////////////////////////////// 30 ///////////////////////////////////////////////////////////////////////////////
31 #include "CCaloOrganization.hh" 31 #include "CCaloOrganization.hh"
32 32
33 //#define debug 33 //#define debug
34 34
35 unsigned int CCaloOrganization::packindex(G4in << 35 unsigned int CCaloOrganization::packindex(int det, int z, int eta,
36 G4in << 36 int phi) const {
37 //So this is the actual encoding of the inde 37 //So this is the actual encoding of the index:
38 //top 4 bits encode Detector type 38 //top 4 bits encode Detector type
39 // 39 //
40 // Should work for all calorimeter with no d 40 // Should work for all calorimeter with no depth information
41 41
42 unsigned int idx=(det&15)<<28; //bits 28-3 42 unsigned int idx=(det&15)<<28; //bits 28-31 (21-23 are free for now)
43 idx+=(((z+1)/2)&1)<<20; //bits 20 43 idx+=(((z+1)/2)&1)<<20; //bits 20
44 idx+=(eta&1023)<<10; //bits 10-1 44 idx+=(eta&1023)<<10; //bits 10-19
45 idx+=(phi&1023); //bits 0-9 45 idx+=(phi&1023); //bits 0-9
46 #ifdef debug 46 #ifdef debug
47 G4cout << " ECAL packing " << det << " " << 47 G4cout << " ECAL packing " << det << " " << z << " " << eta << " " << phi
48 << " into " << idx << G4endl; << 48 << " into " << idx << G4endl;
49 #endif 49 #endif
50 return idx; 50 return idx;
51 } 51 }
52 52
53 unsigned int CCaloOrganization::packindex(G4in << 53 unsigned int CCaloOrganization::packindex(int det, int depth, int z, int eta,
54 G4in << 54 int phi) const {
55 //So this is the actual encoding of the inde 55 //So this is the actual encoding of the index:
56 //top 4 bits encode Detector type 56 //top 4 bits encode Detector type
57 //next 4 bits encode depth information 57 //next 4 bits encode depth information
58 // Should work for all calorimeter with no d 58 // Should work for all calorimeter with no depth information
59 59
60 unsigned int idx=(det&15)<<28; //bits 28-31 60 unsigned int idx=(det&15)<<28; //bits 28-31 (21-23 are free for now)
61 idx+=(depth&15)<<24; //bits 24-27 61 idx+=(depth&15)<<24; //bits 24-27
62 idx+=(z&1)<<20; //bits 20 62 idx+=(z&1)<<20; //bits 20
63 idx+=(eta&1023)<<10; //bits 10-19 63 idx+=(eta&1023)<<10; //bits 10-19
64 idx+=(phi&1023); //bits 0-9 64 idx+=(phi&1023); //bits 0-9
65 #ifdef debug 65 #ifdef debug
66 G4cout << " HCAL packing " << det << " " << 66 G4cout << " HCAL packing " << det << " " << depth << " " << z << " " << eta
67 << " " << phi << " into " << idx << << 67 << " " << phi << " into " << idx << G4endl;
68 #endif 68 #endif
69 return idx; 69 return idx;
70 } 70 }
71 71
72 72
73 void CCaloOrganization::unpackindex(const unsi << 73 void CCaloOrganization::unpackindex(const unsigned int& idx, int& det, int& z,
74 G4int& eta << 74 int& eta, int& phi) const {
75 det = (idx>>28)&15; 75 det = (idx>>28)&15;
76 z = (idx>>20)&1; 76 z = (idx>>20)&1;
77 z = 2*z-1; 77 z = 2*z-1;
78 eta = (idx>>10)&1023; 78 eta = (idx>>10)&1023;
79 phi = (idx&1023); 79 phi = (idx&1023);
80 } 80 }
81 81
82 82
83 void CCaloOrganization::unpackindex(const unsi << 83 void CCaloOrganization::unpackindex(const unsigned int& idx, int& det,
84 G4int& dep << 84 int& depth, int& z, int& eta,
85 G4int& phi << 85 int& phi) const {
86 det = (idx>>28)&15; 86 det = (idx>>28)&15;
87 depth=(idx>>24)&15; 87 depth=(idx>>24)&15;
88 z = (idx>>20)&1; 88 z = (idx>>20)&1;
89 eta = (idx>>10)&1023; 89 eta = (idx>>10)&1023;
90 phi = (idx&1023); 90 phi = (idx&1023);
91 } 91 }
92 92
93 93
94 G4int CCaloOrganization::getUnitWithMaxEnergy( << 94 int CCaloOrganization::getUnitWithMaxEnergy(std::map<int,float,std::less<int> >& themap){
95 95
96 //look for max 96 //look for max
97 G4int UnitWithMaxEnergy = 0; << 97 int UnitWithMaxEnergy = 0;
98 G4float maxEnergy = 0.; << 98 float maxEnergy = 0.;
99 << 99
100 for(std::map<G4int,G4float,std::less<G4int> << 100 for(std::map<int,float,std::less<int> >::iterator iter = themap.begin();
101 iter != themap.end(); iter++){ 101 iter != themap.end(); iter++){
102 << 102
103 if( maxEnergy < (*iter).second) { << 103 if( maxEnergy < (*iter).second) {
104 maxEnergy = (*iter).second; << 104 maxEnergy = (*iter).second;
105 UnitWithMaxEnergy = (*iter).first; 105 UnitWithMaxEnergy = (*iter).first;
106 } << 106 }
107 } << 107 }
108 #ifdef debug 108 #ifdef debug
109 G4cout << " *** max energy of " << maxEnergy 109 G4cout << " *** max energy of " << maxEnergy << " MeV was found in Unit id "
110 << UnitWithMaxEnergy; << 110 << UnitWithMaxEnergy;
111 G4int det,z,eta,phi; << 111 int det,z,eta,phi;
112 unpackindex(UnitWithMaxEnergy, det, z, eta, 112 unpackindex(UnitWithMaxEnergy, det, z, eta, phi);
113 G4cout << " corresponding to z= " << z << " 113 G4cout << " corresponding to z= " << z << " eta= " << eta << " phi = " << phi
114 << G4endl; << 114 << G4endl;
115 #endif 115 #endif
116 return UnitWithMaxEnergy; 116 return UnitWithMaxEnergy;
117 117
118 } 118 }
119 119
120 120
121 G4float CCaloOrganization::energyInMatrix(G4in << 121 float CCaloOrganization::energyInMatrix(int nCellInEta, int nCellInPhi,
122 G4int << 122 int crystalWithMaxEnergy,
123 std::m << 123 std::map<int,float,std::less<int> >& themap){
124 124
125 G4int det,z,eta,phi; << 125 int det,z,eta,phi;
126 this->unpackindex(crystalWithMaxEnergy, det, 126 this->unpackindex(crystalWithMaxEnergy, det, z, eta, phi);
127 #ifdef debug << 127 int ncristals=0;
128 G4int ncristals=0; << 128
129 #endif << 129 int goBackInEta = nCellInEta/2;
130 G4int goBackInEta = nCellInEta/2; << 130 int goBackInPhi = nCellInPhi/2;
131 G4int goBackInPhi = nCellInPhi/2; << 131 int startEta = eta-goBackInEta;
132 G4int startEta = eta-goBackInEta; << 132 int startPhi = phi-goBackInPhi;
133 G4int startPhi = phi-goBackInPhi; <<
134 133
135 G4float totalEnergy = 0.; << 134 float totalEnergy = 0.;
136 135
137 for(G4int ieta=startEta; ieta<startEta+nCell << 136 for(int ieta=startEta; ieta<startEta+nCellInEta; ieta++){
138 for(G4int iphi=startPhi; iphi<startPhi+nCe << 137 for(int iphi=startPhi; iphi<startPhi+nCellInPhi; iphi++){
139 138
140 G4int index = this->packindex(det,z,ieta << 139 int index = this->packindex(det,z,ieta,iphi);
141 totalEnergy += themap[index]; 140 totalEnergy += themap[index];
142 #ifdef debug <<
143 ncristals+=1; 141 ncristals+=1;
>> 142 #ifdef debug
144 G4cout << "ieta - iphi - E = " << ieta < 143 G4cout << "ieta - iphi - E = " << ieta << " " << iphi << " "
145 << themap[index] << G4endl; << 144 << themap[index] << G4endl;
146 #endif 145 #endif
147 } 146 }
148 } 147 }
149 << 148
150 #ifdef debug 149 #ifdef debug
151 G4cout << "Energy in " << nCellInEta << " ce 150 G4cout << "Energy in " << nCellInEta << " cells in eta times "
152 << nCellInPhi << " cells in phi matri << 151 << nCellInPhi << " cells in phi matrix = " << totalEnergy
153 << " for " << ncristals << " crystals << 152 << " for " << ncristals << " crystals" << G4endl;
154 #endif 153 #endif
155 return totalEnergy; 154 return totalEnergy;
156 155
157 } 156 }
158 157