source: trunk/source/geometry/magneticfield/test/field06/src/MyExSteppingAction.cc @ 1347

Last change on this file since 1347 was 1347, checked in by garnier, 13 years ago

geant4 tag 9.4
File size: 8.6 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer                                           *
4// *                                                                  *
5// * The  Geant4 software  is  copyright of the Copyright Holders  of *
6// * the Geant4 Collaboration.  It is provided  under  the terms  and *
7// * conditions of the Geant4 Software License,  included in the file *
8// * LICENSE and available at  http://cern.ch/geant4/license .  These *
9// * include a list of copyright holders.                             *
10// *                                                                  *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work  make  any representation or  warranty, express or implied, *
14// * regarding  this  software system or assume any liability for its *
15// * use.  Please see the license in the file  LICENSE  and URL above *
16// * for the full disclaimer and the limitation of liability.         *
17// *                                                                  *
18// * This  code  implementation is the result of  the  scientific and *
19// * technical work of the GEANT4 collaboration.                      *
20// * By using,  copying,  modifying or  distributing the software (or *
21// * any work based  on the software)  you  agree  to acknowledge its *
22// * use  in  resulting  scientific  publications,  and indicate your *
23// * acceptance of all terms of the Geant4 Software license.          *
24// ********************************************************************
25//
26//
27// $Id: MyExSteppingAction.cc,v 1.1 2006/11/28 15:19:42 tnikitin Exp $
28// GEANT4 tag $Name: geant4-09-04-ref-00 $
29//
30//
31                                                                               
32//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......                                                                               
34#include "ExN01SteppingAction.hh"
35                                                                               
36#include "ExN01DetectorConstruction.hh"
37#include "ExN01EventAction.hh"
38                                                                               
39#include "G4Step.hh"
40 #include "G4UnitsTable.hh"
41                                                                               
42////#include "G4RunManager.hh"
43                                                                               
44//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
45ExN01SteppingAction::ExN01SteppingAction(ExN01DetectorConstruction* det,
46                                          ExN01EventAction* evt)
47:detector(det), event_action(evt)                                         
48{ }                                                                               
49                                                     
50//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......                                                                               
51ExN01SteppingAction::~ExN01SteppingAction()
52{ }
53                                                                               
54//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......                                                                               
55void ExN01SteppingAction::UserSteppingAction(const G4Step* aStep)
56{
57 // get volume of the current step
58  G4VPhysicalVolume* volume
59   = aStep->GetPreStepPoint()->GetTouchableHandle()->GetVolume();
60  //
61  //tanja
62  //get name of particle
63  G4double EndStepX=0.;
64  G4double EndStepY=0.;
65  G4double EndStepZ=0.;
66    G4String particleName = (aStep -> GetTrack() -> GetDynamicParticle()
67                           -> GetDefinition() -> GetParticleName());
68    if( particleName =="e-")
69    {
70      //      G4cout<<particleName<<G4endl;
71     EndStepX =
72       aStep->GetPreStepPoint()->GetPosition().x();
73   
74       EndStepY =
75         aStep->GetPreStepPoint()->GetPosition().y();
76       EndStepZ =
77         aStep->GetPreStepPoint()->GetPosition().z();
78       G4cout<<particleName<<"  "<<"PreStep"<<" "
79            << std::setw(6) << G4BestUnit(aStep->GetPreStepPoint()->GetPosition().x(),"Length")<<"  "
80            << std::setw(6) << G4BestUnit(aStep->GetPreStepPoint()->GetPosition().y(),"Length")<<"  "
81            << std::setw(6) << G4BestUnit(aStep->GetPreStepPoint()->GetPosition().z(),"Length")<<"  "
82              << std::setw(6) << G4BestUnit(std::sqrt(EndStepX*EndStepX+EndStepY*EndStepY),"Length")<<"  " 
83           
84              << std::setw(6) << G4BestUnit(aStep->GetStepLength(),"Length")<<"  "
85            << std::setw(10) <<aStep->GetPreStepPoint()->GetTouchableHandle()->GetVolume()->GetName()<<"  "
86            <<G4endl;
87       //G4cout<<particleName<<"  "<<"PostStep  "
88       //     << std::setw(6) << G4BestUnit(aStep->GetPostStepPoint()->GetPosition().x(),"Length")<<"  "
89       //     << std::setw(6) << G4BestUnit(aStep->GetPostStepPoint()->GetPosition().y(),"Length")<<"  "
90       //     << std::setw(6) << G4BestUnit(aStep->GetPostStepPoint()->GetPosition().z(),"Length")<<"  "
91       //       << std::setw(6) << G4BestUnit(aStep->GetStepLength(),"Length")<<"  "
92       //     
93       ////     <<G4endl;
94       //
95       // const G4int verboseLevel=1;
96       //G4cout.precision(6);
97       // G4cout.setf(ios_base::fixed,ios_base::floatfield);
98       //  G4cout << std::setw( 9) << "X(mm)" << " "
99       //     << std::setw( 9) << "Y(mm)" << " "
100       //     << std::setw( 9) << "Z(mm)" << " "
101       //     << std::setw( 9) << " N_x " << " "
102       //     << std::setw( 9) << " N_y " << " "
103       //     << std::setw( 9) << " N_z " << " "
104       //        //  << std::setw( 9) << " Delta|N|" << " "
105       //  //   << std::setw( 9) << " Delta(N_z) " << " "
106       //  //   << std::setw( 9) << "StepLen" << " "
107       //  //   << std::setw( 9) << "PhsStep" << " "
108       //     << std::setw( 9) << "Safety" << " "
109       //     << std::setw(18) << "NextVolume" << " "
110       //     << G4endl;
111   
112 
113   
114       //G4cout.precision(6);
115       //G4cout <<std::setw( 9) <<aStep->GetPreStepPoint()->GetPosition().x() << " "
116       //     << std::setw( 9) <<aStep->GetPreStepPoint()->GetPosition().y() << " "
117       //     << std::setw( 9) <<aStep->GetPreStepPoint()->GetPosition().z() << " "
118       //        //   << std::setw( 9) <<aStep->GetPreStepPoint()->GetVelocity.x() << " "
119       //  //   << std::setw( 9) << aStep->GetPreStepPoint()->GetVelocity.y() << " "
120       //  //    << std::setw( 9) <<aStep->GetPreStepPoint()->GetVelocity.z() << " ";
121       //        // G4cout.precision(2);
122         
123       ///       // << std::setw( 9) <<aStep->GetPreStepPoint()->GetSafety << " "
124       //    << std::setw(12) << aStep->GetPreStepPoint()->GetTouchableHandle()->GetVolume() << " "
125       //     <<G4endl;                                                                         
126
127
128
129     }
130      if( particleName =="chargedgeantino")
131    {
132      //      G4cout<<particleName<<G4endl;
133     EndStepX =
134       aStep->GetPreStepPoint()->GetPosition().x();
135   
136       EndStepY =
137         aStep->GetPreStepPoint()->GetPosition().y();
138       EndStepZ =
139         aStep->GetPreStepPoint()->GetPosition().z();
140       G4cout<<particleName<<"  "
141            << std::setw(6) << G4BestUnit(aStep->GetPreStepPoint()->GetPosition().x(),"Length")<<"  "
142            << std::setw(6) << G4BestUnit(aStep->GetPreStepPoint()->GetPosition().y(),"Length")<<"  "
143            << std::setw(6) << G4BestUnit(aStep->GetPreStepPoint()->GetPosition().z(),"Length")<<"  "
144            <<G4endl;
145     
146     }
147      if( particleName =="mu+")
148    {
149      //      G4cout<<particleName<<G4endl;
150     EndStepX =
151       aStep->GetPreStepPoint()->GetPosition().x();
152   
153       EndStepY =
154         aStep->GetPreStepPoint()->GetPosition().y();
155       EndStepZ =
156         aStep->GetPreStepPoint()->GetPosition().z();
157       G4cout<<particleName<<"  "
158            << std::setw(6) << G4BestUnit(aStep->GetPreStepPoint()->GetPosition().x(),"Length")<<"  "
159            << std::setw(6) << G4BestUnit(aStep->GetPreStepPoint()->GetPosition().y(),"Length")<<"  "
160            << std::setw(6) << G4BestUnit(aStep->GetPreStepPoint()->GetPosition().z(),"Length")<<"  "
161            <<G4endl;
162     
163     }
164  //end tanja
165  // collect energy and track length step by step
166
167  //end tanja
168  //example of saving random number seed of this event, under condition
169  //// if (condition) G4RunManager::GetRunManager()->rndmSaveThisEvent();
170}
171                                                                               
172//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Note: See TracBrowser for help on using the repository browser.