Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/composite_calorimeter/src/CCaloOrganization.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/advanced/composite_calorimeter/src/CCaloOrganization.cc (Version 11.3.0) and /examples/advanced/composite_calorimeter/src/CCaloOrganization.cc (Version 6.1)


  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