Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/digits_hits/scorer/src/G4PSCylinderSurfaceCurrent.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 /digits_hits/scorer/src/G4PSCylinderSurfaceCurrent.cc (Version 11.3.0) and /digits_hits/scorer/src/G4PSCylinderSurfaceCurrent.cc (Version 9.0.p2)


  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 // G4PSCylinderSurfaceCurrent                     
 29 #include "G4PSCylinderSurfaceCurrent.hh"          
 30                                                   
 31 #include "G4SystemOfUnits.hh"                     
 32 #include "G4StepStatus.hh"                        
 33 #include "G4Track.hh"                             
 34 #include "G4VSolid.hh"                            
 35 #include "G4VPhysicalVolume.hh"                   
 36 #include "G4VPVParameterisation.hh"               
 37 #include "G4UnitsTable.hh"                        
 38 #include "G4GeometryTolerance.hh"                 
 39 #include "G4VScoreHistFiller.hh"                  
 40                                                   
 41 //////////////////////////////////////////////    
 42 // (Description)                                  
 43 //   This is a primitive scorer class for scor    
 44 //  Current version assumes only for G4Tubs sh    
 45 //                                                
 46 // Surface is defined at the inner surface of     
 47 // Direction                   R    R+dR          
 48 //   0  IN || OUT            ->|<-  |             
 49 //   1  IN                   ->|    |             
 50 //   2  OUT                    |<-  |             
 51 //                                                
 52 // Created: 2007-03-21  Tsukasa ASO               
 53 // 2010-07-22   Introduce Unit specification.     
 54 // 2020-10-06   Use G4VPrimitivePlotter and fi    
 55 //              vs. surface current * track we    
 56 //                                                
 57 //////////////////////////////////////////////    
 58                                                   
 59 G4PSCylinderSurfaceCurrent::G4PSCylinderSurfac    
 60                                                   
 61                                                   
 62   : G4PSCylinderSurfaceCurrent(name, direction    
 63 {}                                                
 64                                                   
 65 G4PSCylinderSurfaceCurrent::G4PSCylinderSurfac    
 66                                                   
 67                                                   
 68                                                   
 69   : G4VPrimitivePlotter(name, depth)              
 70   , HCID(-1)                                      
 71   , fDirection(direction)                         
 72   , EvtMap(nullptr)                               
 73   , weighted(true)                                
 74   , divideByArea(true)                            
 75 {                                                 
 76   DefineUnitAndCategory();                        
 77   SetUnit(unit);                                  
 78 }                                                 
 79                                                   
 80 G4bool G4PSCylinderSurfaceCurrent::ProcessHits    
 81                                                   
 82 {                                                 
 83   G4StepPoint* preStep = aStep->GetPreStepPoin    
 84   G4VSolid* solid      = ComputeCurrentSolid(a    
 85   auto  tubsSolid    = static_cast<G4Tubs*>(so    
 86                                                   
 87   G4int dirFlag = IsSelectedSurface(aStep, tub    
 88   // G4cout << " pos " << preStep->GetPosition    
 89   if(dirFlag > 0)                                 
 90   {                                               
 91     if(fDirection == fCurrent_InOut || fDirect    
 92     {                                             
 93       G4TouchableHandle theTouchable = preStep    
 94       //                                          
 95       G4double current = 1.0;                     
 96       if(weighted)                                
 97         current = preStep->GetWeight();  // Cu    
 98       //                                          
 99       if(divideByArea)                            
100       {                                           
101         G4double square = 2. * tubsSolid->GetZ    
102                           tubsSolid->GetInnerR    
103                           tubsSolid->GetDeltaP    
104         current = current / square;  // Curren    
105       }                                           
106                                                   
107       G4int index = GetIndex(aStep);              
108       EvtMap->add(index, current);                
109                                                   
110       if(!hitIDMap.empty() && hitIDMap.find(in    
111       {                                           
112         auto filler = G4VScoreHistFiller::Inst    
113         if(filler == nullptr)                     
114         {                                         
115           G4Exception("G4PSCylinderSurfaceCurr    
116                       JustWarning,                
117                       "G4TScoreHistFiller is n    
118                       "not filled.");             
119         }                                         
120         else                                      
121         {                                         
122           filler->FillH1(hitIDMap[index], preS    
123         }                                         
124       }                                           
125     }                                             
126   }                                               
127                                                   
128   return true;                                    
129 }                                                 
130                                                   
131 G4int G4PSCylinderSurfaceCurrent::IsSelectedSu    
132                                                   
133 {                                                 
134   G4TouchableHandle theTouchable =                
135     aStep->GetPreStepPoint()->GetTouchableHand    
136   G4double kCarTolerance =                        
137     G4GeometryTolerance::GetInstance()->GetSur    
138                                                   
139   if(aStep->GetPreStepPoint()->GetStepStatus()    
140   {                                               
141     // Entering Geometry                          
142     G4ThreeVector stppos1 = aStep->GetPreStepP    
143     G4ThreeVector localpos1 =                     
144       theTouchable->GetHistory()->GetTopTransf    
145     if(std::fabs(localpos1.z()) > tubsSolid->G    
146       return -1;                                  
147     G4double localR2 =                            
148       localpos1.x() * localpos1.x() + localpos    
149     G4double InsideRadius = tubsSolid->GetInne    
150     if(localR2 >                                  
151          (InsideRadius - kCarTolerance) * (Ins    
152        localR2 <                                  
153          (InsideRadius + kCarTolerance) * (Ins    
154     {                                             
155       return fCurrent_In;                         
156     }                                             
157   }                                               
158                                                   
159   if(aStep->GetPostStepPoint()->GetStepStatus(    
160   {                                               
161     // Exiting Geometry                           
162     G4ThreeVector stppos2 = aStep->GetPostStep    
163     G4ThreeVector localpos2 =                     
164       theTouchable->GetHistory()->GetTopTransf    
165     if(std::fabs(localpos2.z()) > tubsSolid->G    
166       return -1;                                  
167     G4double localR2 =                            
168       localpos2.x() * localpos2.x() + localpos    
169     G4double InsideRadius = tubsSolid->GetInne    
170     if(localR2 >                                  
171          (InsideRadius - kCarTolerance) * (Ins    
172        localR2 <                                  
173          (InsideRadius + kCarTolerance) * (Ins    
174     {                                             
175       return fCurrent_Out;                        
176     }                                             
177   }                                               
178                                                   
179   return -1;                                      
180 }                                                 
181                                                   
182 void G4PSCylinderSurfaceCurrent::Initialize(G4    
183 {                                                 
184   EvtMap = new G4THitsMap<G4double>(detector->    
185   if(HCID < 0)                                    
186     HCID = GetCollectionID(0);                    
187   HCE->AddHitsCollection(HCID, (G4VHitsCollect    
188 }                                                 
189                                                   
190 void G4PSCylinderSurfaceCurrent::clear() { Evt    
191                                                   
192 void G4PSCylinderSurfaceCurrent::PrintAll()       
193 {                                                 
194   G4cout << " MultiFunctionalDet  " << detecto    
195   G4cout << " PrimitiveScorer " << GetName() <    
196   G4cout << " Number of entries " << EvtMap->e    
197   for(const auto& [copy, current] : *(EvtMap->    
198   {                                               
199     G4cout << "  copy no.: " << copy << "  cur    
200     if(divideByArea)                              
201     {                                             
202       G4cout << *(current) / GetUnitValue() <<    
203     }                                             
204     else                                          
205     {                                             
206       G4cout << *(current) << " [tracks]";        
207     }                                             
208     G4cout << G4endl;                             
209   }                                               
210 }                                                 
211                                                   
212 void G4PSCylinderSurfaceCurrent::SetUnit(const    
213 {                                                 
214   if(divideByArea)                                
215   {                                               
216     CheckAndSetUnit(unit, "Per Unit Surface");    
217   }                                               
218   else                                            
219   {                                               
220     if(unit.empty())                              
221     {                                             
222       unitName  = unit;                           
223       unitValue = 1.0;                            
224     }                                             
225     else                                          
226     {                                             
227       G4String msg = "Invalid unit [" + unit +    
228                      GetUnit() + "] ) for " +     
229       G4Exception("G4PSCylinderSurfaceCurrent:    
230                   JustWarning, msg);              
231     }                                             
232   }                                               
233 }                                                 
234                                                   
235 void G4PSCylinderSurfaceCurrent::DefineUnitAnd    
236 {                                                 
237   // Per Unit Surface                             
238   new G4UnitDefinition("percentimeter2", "perc    
239                        (1. / cm2));               
240   new G4UnitDefinition("permillimeter2", "perm    
241                        (1. / mm2));               
242   new G4UnitDefinition("permeter2", "perm2", "    
243 }                                                 
244