Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/purging_magnet/src/PurgMagSteppingVerbose.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/purging_magnet/src/PurgMagSteppingVerbose.cc (Version 11.3.0) and /examples/advanced/purging_magnet/src/PurgMagSteppingVerbose.cc (Version 9.6)


  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 // Code developed by:                              26 // Code developed by:
 27 //  S.Larsson                                      27 //  S.Larsson
 28 //                                                 28 //
 29 //    ************************************         29 //    ************************************
 30 //    *                                  *         30 //    *                                  *
 31 //    *    PurgMagSteppingVerbose.cc     *         31 //    *    PurgMagSteppingVerbose.cc     *
 32 //    *                                  *         32 //    *                                  *
 33 //    ************************************         33 //    ************************************
 34 //                                                 34 //
                                                   >>  35 // $Id$
 35 //                                                 36 //
 36                                                    37 
 37 #include "PurgMagSteppingVerbose.hh"               38 #include "PurgMagSteppingVerbose.hh"
 38 #include "G4SteppingManager.hh"                    39 #include "G4SteppingManager.hh"
 39                                                    40 
 40 #include "G4UnitsTable.hh"                         41 #include "G4UnitsTable.hh"
 41                                                    42 
 42 //....oooOO0OOooo........oooOO0OOooo........oo     43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 43                                                    44 
 44 PurgMagSteppingVerbose::PurgMagSteppingVerbose     45 PurgMagSteppingVerbose::PurgMagSteppingVerbose()
 45 {}                                                 46 {}
 46                                                    47 
 47 //....oooOO0OOooo........oooOO0OOooo........oo     48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 48                                                    49 
 49 PurgMagSteppingVerbose::~PurgMagSteppingVerbos     50 PurgMagSteppingVerbose::~PurgMagSteppingVerbose()
 50 {}                                                 51 {}
 51                                                    52 
 52 //....oooOO0OOooo........oooOO0OOooo........oo     53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 53                                                    54 
 54 void PurgMagSteppingVerbose::StepInfo()            55 void PurgMagSteppingVerbose::StepInfo()
 55 {                                                  56 {
 56   CopyState();                                     57   CopyState();
 57                                                    58   
 58   G4int prec = G4cout.precision(3);                59   G4int prec = G4cout.precision(3);
 59                                                    60   
 60   if( verboseLevel >= 1 ){                         61   if( verboseLevel >= 1 ){
 61     if( verboseLevel >= 4 ) VerboseTrack();        62     if( verboseLevel >= 4 ) VerboseTrack();
 62     if( verboseLevel >= 3 ){                       63     if( verboseLevel >= 3 ){
 63       G4cout << G4endl;                            64       G4cout << G4endl;    
 64       G4cout << std::setw( 5) << "#Step#"          65       G4cout << std::setw( 5) << "#Step#"     << " "
 65        << std::setw( 6) << "X"          << "       66        << std::setw( 6) << "X"          << "    "
 66        << std::setw( 6) << "Y"          << "       67        << std::setw( 6) << "Y"          << "    "  
 67        << std::setw( 6) << "Z"          << "       68        << std::setw( 6) << "Z"          << "    "
 68        << std::setw( 9) << "KineE"      << " "     69        << std::setw( 9) << "KineE"      << " "
 69        << std::setw( 9) << "dEStep"     << " "     70        << std::setw( 9) << "dEStep"     << " "  
 70        << std::setw(10) << "StepLeng"              71        << std::setw(10) << "StepLeng"     
 71        << std::setw(10) << "TrakLeng"              72        << std::setw(10) << "TrakLeng" 
 72        << std::setw(10) << "NextVolu"              73        << std::setw(10) << "NextVolu" 
 73        << std::setw(10) << "Process"   << G4en     74        << std::setw(10) << "Process"   << G4endl;           
 74     }                                              75     }
 75                                                    76     
 76     G4cout << std::setw( 5) << fTrack->GetCurr     77     G4cout << std::setw( 5) << fTrack->GetCurrentStepNumber() << " "
 77      << std::setw( 6) << G4BestUnit(fTrack->Ge     78      << std::setw( 6) << G4BestUnit(fTrack->GetPosition().x(),"Length")
 78      << std::setw( 6) << G4BestUnit(fTrack->Ge     79      << std::setw( 6) << G4BestUnit(fTrack->GetPosition().y(),"Length")
 79      << std::setw( 6) << G4BestUnit(fTrack->Ge     80      << std::setw( 6) << G4BestUnit(fTrack->GetPosition().z(),"Length")
 80      << std::setw( 6) << G4BestUnit(fTrack->Ge     81      << std::setw( 6) << G4BestUnit(fTrack->GetKineticEnergy(),"Energy")
 81      << std::setw( 6) << G4BestUnit(fStep->Get     82      << std::setw( 6) << G4BestUnit(fStep->GetTotalEnergyDeposit(),"Energy")
 82      << std::setw( 6) << G4BestUnit(fStep->Get     83      << std::setw( 6) << G4BestUnit(fStep->GetStepLength(),"Length")
 83      << std::setw( 6) << G4BestUnit(fTrack->Ge     84      << std::setw( 6) << G4BestUnit(fTrack->GetTrackLength(),"Length");
 84                                                    85     
 85     // if( fStepStatus != fWorldBoundary){         86     // if( fStepStatus != fWorldBoundary){ 
 86     if( fTrack->GetNextVolume() != 0 ) {           87     if( fTrack->GetNextVolume() != 0 ) { 
 87       G4cout << std::setw(10) << fTrack->GetNe     88       G4cout << std::setw(10) << fTrack->GetNextVolume()->GetName();
 88     } else {                                       89     } else {
 89       G4cout << std::setw(10) << "OutOfWorld";     90       G4cout << std::setw(10) << "OutOfWorld";
 90     }                                              91     }
 91                                                    92     
 92     if(fStep->GetPostStepPoint()->GetProcessDe     93     if(fStep->GetPostStepPoint()->GetProcessDefinedStep() != NULL){
 93       G4cout << std::setw(10) << fStep->GetPos     94       G4cout << std::setw(10) << fStep->GetPostStepPoint()->GetProcessDefinedStep()
 94   ->GetProcessName();                              95   ->GetProcessName();
 95     } else {                                       96     } else {
 96       G4cout << "User Limit";                      97       G4cout << "User Limit";
 97     }                                              98     }
 98                                                    99     
 99     G4cout << G4endl;                             100     G4cout << G4endl;
100                                                   101     
101     if( verboseLevel == 2 ){                      102     if( verboseLevel == 2 ){
102       G4int tN2ndariesTot = fN2ndariesAtRestDo    103       G4int tN2ndariesTot = fN2ndariesAtRestDoIt +
103                       fN2ndariesAlongStepDoIt     104                       fN2ndariesAlongStepDoIt +
104                       fN2ndariesPostStepDoIt;     105                       fN2ndariesPostStepDoIt;
105       if(tN2ndariesTot>0){                        106       if(tN2ndariesTot>0){
106   G4cout << "    :----- List of 2ndaries - "      107   G4cout << "    :----- List of 2ndaries - "
107          << "#SpawnInStep=" << std::setw(3) <<    108          << "#SpawnInStep=" << std::setw(3) << tN2ndariesTot 
108          << "(Rest="  << std::setw(2) << fN2nd    109          << "(Rest="  << std::setw(2) << fN2ndariesAtRestDoIt
109          << ",Along=" << std::setw(2) << fN2nd    110          << ",Along=" << std::setw(2) << fN2ndariesAlongStepDoIt
110          << ",Post="  << std::setw(2) << fN2nd    111          << ",Post="  << std::setw(2) << fN2ndariesPostStepDoIt
111          << "), "                                 112          << "), "
112          << "#SpawnTotal=" << std::setw(3) <<     113          << "#SpawnTotal=" << std::setw(3) << (*fSecondary).size()
113          << " ---------------"                    114          << " ---------------"
114          << G4endl;                               115          << G4endl;
115                                                   116   
116   for(size_t lp1=(*fSecondary).size()-tN2ndari    117   for(size_t lp1=(*fSecondary).size()-tN2ndariesTot; 
117       lp1<(*fSecondary).size(); lp1++){           118       lp1<(*fSecondary).size(); lp1++){
118     G4cout << "    : "                            119     G4cout << "    : "
119      << std::setw(6)                              120      << std::setw(6)
120      << G4BestUnit((*fSecondary)[lp1]->GetPosi    121      << G4BestUnit((*fSecondary)[lp1]->GetPosition().x(),"Length")
121      << std::setw(6)                              122      << std::setw(6)
122      << G4BestUnit((*fSecondary)[lp1]->GetPosi    123      << G4BestUnit((*fSecondary)[lp1]->GetPosition().y(),"Length")
123      << std::setw(6)                              124      << std::setw(6)
124      << G4BestUnit((*fSecondary)[lp1]->GetPosi    125      << G4BestUnit((*fSecondary)[lp1]->GetPosition().z(),"Length")
125      << std::setw(6)                              126      << std::setw(6)
126      << G4BestUnit((*fSecondary)[lp1]->GetKine    127      << G4BestUnit((*fSecondary)[lp1]->GetKineticEnergy(),"Energy")
127      << std::setw(10)                             128      << std::setw(10)
128      << (*fSecondary)[lp1]->GetDefinition()->G    129      << (*fSecondary)[lp1]->GetDefinition()->GetParticleName();
129     G4cout << G4endl;                             130     G4cout << G4endl;
130   }                                               131   }
131                                                   132               
132   G4cout << "    :----------------------------    133   G4cout << "    :-----------------------------"
133          << "---------------------------------    134          << "----------------------------------"
134          << "-- EndOf2ndaries Info -----------    135          << "-- EndOf2ndaries Info ---------------"
135          << G4endl;                               136          << G4endl;
136       }                                           137       }
137     }                                             138     }
138                                                   139     
139   }                                               140   }
140   G4cout.precision(prec);                         141   G4cout.precision(prec);
141   G4cout<< "exit PurgMagSteppingVerbose::StepI    142   G4cout<< "exit PurgMagSteppingVerbose::StepInfo   " <<G4endl;
142 }                                                 143 }
143                                                   144 
144 //....oooOO0OOooo........oooOO0OOooo........oo    145 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
145                                                   146 
146 void PurgMagSteppingVerbose::TrackingStarted()    147 void PurgMagSteppingVerbose::TrackingStarted()
147 {                                                 148 {
148   CopyState();                                    149   CopyState();
149 G4int prec = G4cout.precision(3);                 150 G4int prec = G4cout.precision(3);
150   if( verboseLevel > 0 ){                         151   if( verboseLevel > 0 ){
151                                                   152 
152     G4cout << std::setw( 5) << "Step#"      <<    153     G4cout << std::setw( 5) << "Step#"      << " "
153            << std::setw( 6) << "X"          <<    154            << std::setw( 6) << "X"          << "    "
154      << std::setw( 6) << "Y"          << "        155      << std::setw( 6) << "Y"          << "    "  
155      << std::setw( 6) << "Z"          << "        156      << std::setw( 6) << "Z"          << "    "
156      << std::setw( 9) << "KineE"      << " "      157      << std::setw( 9) << "KineE"      << " "
157      << std::setw( 9) << "dEStep"     << " "      158      << std::setw( 9) << "dEStep"     << " "  
158      << std::setw(10) << "StepLeng"               159      << std::setw(10) << "StepLeng"  
159      << std::setw(10) << "TrakLeng"               160      << std::setw(10) << "TrakLeng"
160      << std::setw(10) << "NextVolu"               161      << std::setw(10) << "NextVolu"
161      << std::setw(10) << "Process"    << G4end    162      << std::setw(10) << "Process"    << G4endl;       
162                                                   163 
163     G4cout << std::setw( 5) << fTrack->GetCurr    164     G4cout << std::setw( 5) << fTrack->GetCurrentStepNumber() << " "
164      << std::setw( 6) << G4BestUnit(fTrack->Ge    165      << std::setw( 6) << G4BestUnit(fTrack->GetPosition().x(),"Length")
165      << std::setw( 6) << G4BestUnit(fTrack->Ge    166      << std::setw( 6) << G4BestUnit(fTrack->GetPosition().y(),"Length")
166      << std::setw( 6) << G4BestUnit(fTrack->Ge    167      << std::setw( 6) << G4BestUnit(fTrack->GetPosition().z(),"Length")
167      << std::setw( 6) << G4BestUnit(fTrack->Ge    168      << std::setw( 6) << G4BestUnit(fTrack->GetKineticEnergy(),"Energy")
168      << std::setw( 6) << G4BestUnit(fStep->Get    169      << std::setw( 6) << G4BestUnit(fStep->GetTotalEnergyDeposit(),"Energy")
169      << std::setw( 6) << G4BestUnit(fStep->Get    170      << std::setw( 6) << G4BestUnit(fStep->GetStepLength(),"Length")
170      << std::setw( 6) << G4BestUnit(fTrack->Ge    171      << std::setw( 6) << G4BestUnit(fTrack->GetTrackLength(),"Length");
171                                                   172 
172     if(fTrack->GetNextVolume()){                  173     if(fTrack->GetNextVolume()){
173       G4cout << std::setw(10) << fTrack->GetNe    174       G4cout << std::setw(10) << fTrack->GetNextVolume()->GetName() << " ";
174     } else {                                      175     } else {
175       G4cout << std::setw(10) << "OutOfWorld"     176       G4cout << std::setw(10) << "OutOfWorld" << " ";
176     }                                             177     }
177     G4cout << std::setw(10) << "initStep" << G    178     G4cout << std::setw(10) << "initStep" << G4endl;
178   }                                               179   }
179   G4cout.precision(prec);                         180   G4cout.precision(prec);
180   G4cout<< "exit PurgMagSteppingVerbose::Track    181   G4cout<< "exit PurgMagSteppingVerbose::TrackingStarted()   " <<G4endl;
181 }                                                 182 }
182                                                   183