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