Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 ////////////////////////////////////////////// 23 /////////////////////////////////////////////////////////////////////////////// 27 // File: CCalSeppingAction.cc 24 // File: CCalSeppingAction.cc 28 // Description: Study profiling during the ste 25 // Description: Study profiling during the steps 29 ////////////////////////////////////////////// 26 /////////////////////////////////////////////////////////////////////////////// 30 #include <iostream> << 31 << 32 #include "CCalSteppingAction.hh" 27 #include "CCalSteppingAction.hh" 33 #include "CCalRunAction.hh" << 34 28 35 #include "G4AnalysisManager.hh" << 36 #include "G4SystemOfUnits.hh" << 37 #include "G4SDManager.hh" 29 #include "G4SDManager.hh" 38 #include "G4StepPoint.hh" 30 #include "G4StepPoint.hh" >> 31 #include <iostream> 39 #include "G4ThreeVector.hh" 32 #include "G4ThreeVector.hh" 40 #include "G4RunManager.hh" << 33 #include <math.h> >> 34 >> 35 #ifdef G4ANALYSIS_USE >> 36 #include "CCalAnalysis.hh" >> 37 #endif >> 38 >> 39 >> 40 CCalSteppingAction::CCalSteppingAction(){ >> 41 >> 42 #ifdef G4ANALYSIS_USE >> 43 CCalAnalysis* analysis = CCalAnalysis::getInstance(); >> 44 timeHistoMaxBin=analysis->maxbin(); >> 45 #else >> 46 timeHistoMaxBin=1; >> 47 #endif >> 48 >> 49 int i; >> 50 for (i=0; i<200; i++) {timeDeposit[i] = 0.;} >> 51 for (i=0; i<70; i++) {LateralProfile[i] = 0.;} 41 52 42 CCalSteppingAction::CCalSteppingAction() << 43 { << 44 timeHistoMaxBin = 200; << 45 for (G4int i=0; i<200; i++) {timeDeposit[i] << 46 for (G4int i=0; i<70; i++) {LateralProfile[ << 47 } 53 } 48 54 >> 55 49 CCalSteppingAction::~CCalSteppingAction(){ 56 CCalSteppingAction::~CCalSteppingAction(){ >> 57 G4cout << "CCalSteppingAction deleted" << G4endl; 50 } 58 } 51 59 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 60 67 G4float fde = (G4float)(de/CLHEP::GeV); << 61 void CCalSteppingAction::UserSteppingAction(const G4Step* aStep){ 68 //G4cout << " TSliceID= " << TSliceID << " << 69 62 70 timeDeposit[TSliceID] += fde; << 63 G4StepPoint* PostStepPoint= aStep->GetPostStepPoint(); >> 64 G4StepPoint* PreStepPoint= aStep->GetPreStepPoint(); >> 65 int TSliceID; >> 66 >> 67 TSliceID = static_cast<int>( (PostStepPoint->GetGlobalTime() ) / nanosecond); >> 68 TSliceID = TSliceID<timeHistoMaxBin ? TSliceID : timeHistoMaxBin-1; >> 69 timeDeposit[TSliceID] += aStep->GetTotalEnergyDeposit() / GeV; 71 70 72 G4ThreeVector HitPoint = 0.5*(PostStepPoint- 71 G4ThreeVector HitPoint = 0.5*(PostStepPoint->GetPosition()+ 73 PreStepPoint-> << 72 PreStepPoint->GetPosition()); 74 // Because the beam axis has been defined as 73 // Because the beam axis has been defined as the x-axis, 75 // the lateral displacement is given in term 74 // the lateral displacement is given in terms of the y and z positions. 76 G4double perp = << 75 double perp = std::sqrt(HitPoint.y()*HitPoint.y()+HitPoint.z()*HitPoint.z()); 77 std::sqrt(HitPoint.y()*HitPoint.y()+HitPoi << 76 int radialPosition = std::min(69,int(perp/cm)); 78 G4int ir = (G4int)perp; << 77 LateralProfile[radialPosition] += aStep->GetTotalEnergyDeposit() / GeV; 79 //G4cout << " perp= " << perp << " ir= " << << 78 80 G4int radialPosition = std::max(0,std::min(6 << 81 LateralProfile[radialPosition] += fde; << 82 //G4cout << " done " << G4endl; << 83 } 79 } 84 80 >> 81 85 void CCalSteppingAction::endOfEvent(){ 82 void CCalSteppingAction::endOfEvent(){ 86 83 87 G4AnalysisManager* man = G4AnalysisManager:: << 84 #ifdef G4ANALYSIS_USE 88 #ifdef debug << 85 CCalAnalysis* analysis = CCalAnalysis::getInstance(); 89 G4double totalFilledProfileHcal = 0.0; << 86 analysis->InsertLateralProfile(LateralProfile); 90 #endif << 87 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 88 #endif 119 for (G4int j=0; j<timeHistoMaxBin; ++j) { << 89 120 man->FillH1(IDTimeHist+j,timeDeposit[j]); << 90 int i=0; 121 #ifdef debug << 91 for (i=0; i<70; i++){LateralProfile[i] = 0.;} 122 G4cout << "Fill Time slice histo " << j << << 92 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 << 140 G4cout << " --- End of event --- " << G4endl << 141 93 142 } 94 } 143 95