Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/composite_calorimeter/src/CCalSteppingAction.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/CCalSteppingAction.cc (Version 11.3.0) and /examples/advanced/composite_calorimeter/src/CCalSteppingAction.cc (Version 10.0)


  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: CCalSeppingAction.cc                      27 // File: CCalSeppingAction.cc
 28 // Description: Study profiling during the ste     28 // Description: Study profiling during the steps
 29 //////////////////////////////////////////////     29 ///////////////////////////////////////////////////////////////////////////////
 30 #include <iostream>                                30 #include <iostream>
                                                   >>  31 #include <cmath>
 31                                                    32 
 32 #include "CCalSteppingAction.hh"                   33 #include "CCalSteppingAction.hh"
 33 #include "CCalRunAction.hh"                    << 
 34                                                    34 
 35 #include "G4AnalysisManager.hh"                << 
 36 #include "G4SystemOfUnits.hh"                      35 #include "G4SystemOfUnits.hh"
 37 #include "G4SDManager.hh"                          36 #include "G4SDManager.hh"
 38 #include "G4StepPoint.hh"                          37 #include "G4StepPoint.hh"
 39 #include "G4ThreeVector.hh"                        38 #include "G4ThreeVector.hh"
 40 #include "G4RunManager.hh"                     << 
 41                                                    39 
 42 CCalSteppingAction::CCalSteppingAction()       <<  40 #ifdef G4ANALYSIS_USE  
 43 {                                              <<  41 #include "CCalAnalysis.hh"
 44   timeHistoMaxBin = 200;                       <<  42 #endif
 45   for (G4int i=0; i<200; i++) {timeDeposit[i]  <<  43 
 46   for (G4int i=0; i<70;  i++) {LateralProfile[ <<  44 
                                                   >>  45 CCalSteppingAction::CCalSteppingAction(){
                                                   >>  46 
                                                   >>  47 #ifdef G4ANALYSIS_USE  
                                                   >>  48   CCalAnalysis* analysis = CCalAnalysis::getInstance();
                                                   >>  49   timeHistoMaxBin=analysis->maxbin();
                                                   >>  50 #else
                                                   >>  51   timeHistoMaxBin=1;
                                                   >>  52 #endif
                                                   >>  53 
                                                   >>  54   int i; 
                                                   >>  55   for (i=0; i<200; i++) {timeDeposit[i] = 0.;}
                                                   >>  56   for (i=0; i<70;  i++) {LateralProfile[i] = 0.;}
                                                   >>  57 
 47 }                                                  58 }
 48                                                    59 
                                                   >>  60 
 49 CCalSteppingAction::~CCalSteppingAction(){         61 CCalSteppingAction::~CCalSteppingAction(){
                                                   >>  62   G4cout << "CCalSteppingAction deleted" << G4endl;
 50 }                                                  63 }
 51                                                    64   
 52 void CCalSteppingAction::UserSteppingAction(co << 
 53 {                                              << 
 54   G4double de = aStep->GetTotalEnergyDeposit() << 
 55   if(de < CLHEP::eV) { return; }               << 
 56                                                << 
 57   const G4StepPoint*  PostStepPoint= aStep->Ge << 
 58   const G4StepPoint*  PreStepPoint= aStep->Get << 
 59   G4double time =                              << 
 60     (PostStepPoint) ? PostStepPoint->GetGlobal << 
 61                                                << 
 62   G4int it = (G4int)time;                      << 
 63   it = std::min(it, 10000);                    << 
 64   //G4cout << "## time= " << time << " it= " < << 
 65   G4int TSliceID = std::max(0,std::min(it,time << 
 66                                                    65 
 67   G4float fde = (G4float)(de/CLHEP::GeV);      <<  66 void CCalSteppingAction::UserSteppingAction(const G4Step* aStep){
 68   //G4cout << "  TSliceID= " << TSliceID << "  << 
 69                                                    67 
 70   timeDeposit[TSliceID] += fde;                <<  68   G4StepPoint*  PostStepPoint= aStep->GetPostStepPoint(); 
                                                   >>  69   G4StepPoint*  PreStepPoint= aStep->GetPreStepPoint(); 
                                                   >>  70   int TSliceID;
                                                   >>  71 
                                                   >>  72   if ( PostStepPoint->GetGlobalTime() / nanosecond > 1.0E9 ) TSliceID = 999999999;
                                                   >>  73   else TSliceID = static_cast<int>( PostStepPoint->GetGlobalTime() / nanosecond );
                                                   >>  74   TSliceID = TSliceID<timeHistoMaxBin ? TSliceID : timeHistoMaxBin-1;
                                                   >>  75   timeDeposit[TSliceID] += aStep->GetTotalEnergyDeposit() / GeV;
 71                                                    76 
 72   G4ThreeVector HitPoint = 0.5*(PostStepPoint-     77   G4ThreeVector HitPoint = 0.5*(PostStepPoint->GetPosition()+
 73                                 PreStepPoint-> <<  78         PreStepPoint->GetPosition()); 
 74   // Because the beam axis has been defined as     79   // Because the beam axis has been defined as the x-axis, 
 75   // the lateral displacement is given in term     80   // the lateral displacement is given in terms of the y and z positions. 
 76   G4double perp =                              <<  81   double perp = std::sqrt(HitPoint.y()*HitPoint.y()+HitPoint.z()*HitPoint.z());
 77     std::sqrt(HitPoint.y()*HitPoint.y()+HitPoi <<  82   int radialPosition = std::min(69,int(perp/cm));
 78   G4int ir = (G4int)perp;                      <<  83   LateralProfile[radialPosition] += aStep->GetTotalEnergyDeposit() / GeV;
 79   //G4cout << "  perp= " << perp << " ir= " << <<  84   
 80   G4int radialPosition = std::max(0,std::min(6 << 
 81   LateralProfile[radialPosition] += fde;       << 
 82   //G4cout << "  done " << G4endl;             << 
 83 }                                                  85 }
 84                                                    86 
                                                   >>  87 
 85 void CCalSteppingAction::endOfEvent(){             88 void CCalSteppingAction::endOfEvent(){
 86                                                    89 
 87   G4AnalysisManager* man = G4AnalysisManager:: <<  90 #ifdef G4ANALYSIS_USE  
 88 #ifdef debug                                   <<  91   CCalAnalysis* analysis = CCalAnalysis::getInstance();
 89   G4double totalFilledProfileHcal = 0.0;       <<  92   analysis->InsertLateralProfile(LateralProfile);  
 90 #endif                                         <<  93   analysis->InsertTime(timeDeposit); 
 91   static G4int IDlateralProfile = -1;          << 
 92   if (IDlateralProfile < 0) {                  << 
 93     IDlateralProfile = man->GetH1Id("h500");   << 
 94   }                                            << 
 95   for (G4int i=0; i<70; ++i) {                 << 
 96     man->FillH1(IDlateralProfile+i,LateralProf << 
 97 #ifdef debug                                   << 
 98     G4cout << "Fill Profile Hcal histo " << i  << 
 99     totalFilledProfileHcal += LateralProfile[i << 
100 #endif                                         << 
101   }                                            << 
102                                                << 
103 #ifdef debug                                   << 
104   G4cout << "CCalAnalysis::InsertLateralProfil << 
105          << " histo " << totalFilledProfileHca << 
106 #endif                                         << 
107                                                << 
108   static G4int IDTimeHist = -1;                << 
109   if (IDTimeHist < 0) {                        << 
110     IDTimeHist = man->GetH1Id("h300");         << 
111   }                                            << 
112   static G4int IDTimeProfile = -1;             << 
113   if (IDTimeProfile < 0) {                     << 
114     IDTimeProfile = man->GetH1Id("h901");      << 
115   }                                            << 
116 #ifdef debug                                   << 
117   G4double totalFilledTimeProfile = 0.0;       << 
118 #endif                                             94 #endif
119   for (G4int j=0; j<timeHistoMaxBin; ++j) {    <<  95   
120     man->FillH1(IDTimeHist+j,timeDeposit[j]);  <<  96   int i=0;
121 #ifdef debug                                   <<  97   for (i=0; i<70; i++){LateralProfile[i] = 0.;}
122     G4cout << "Fill Time slice histo " << j << <<  98   for (i=0; i<200; i++){timeDeposit[i] = 0.;}
123     totalFilledTimeProfile += timeDeposit[j];  << 
124 #endif                                         << 
125                                                << 
126     G4double t = j + 0.5;                      << 
127     man->FillH1(IDTimeProfile+1,t,timeDeposit[ << 
128 #ifdef debug                                   << 
129       G4cout << "Fill Time profile histo 1 wit << 
130 #endif                                         << 
131   }                                            << 
132 #ifdef debug                                   << 
133   G4cout << "CCalAnalysis::InsertTime: Total f << 
134          << totalFilledTimeProfile << G4endl;  << 
135 #endif                                         << 
136                                                << 
137   for (G4int i=0; i<70; i++){LateralProfile[i] << 
138   for (G4int i=0; i<200; i++){timeDeposit[i] = << 
139                                                    99 
140   G4cout << " --- End of event --- " << G4endl    100   G4cout << " --- End of event --- " << G4endl;
141                                                   101   
142 }                                                 102 }  
143                                                   103