Geant4 Cross Reference

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


  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 //    *    PurgMagSteppingAction.cc     *          31 //    *    PurgMagSteppingAction.cc     *
 32 //    *                                 *          32 //    *                                 *
 33 //    ***********************************          33 //    ***********************************
 34 //                                                 34 //
                                                   >>  35 // $Id: PurgMagSteppingAction.cc,v 1.5 2006/06/29 16:06:21 gunter Exp $
                                                   >>  36 // GEANT4 tag $Name: geant4-09-01-patch-02 $
 35 //                                                 37 //
 36                                                    38 
 37 //....oooOO0OOooo........oooOO0OOooo........oo     39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 38                                                    40 
 39 #include "PurgMagSteppingAction.hh"                41 #include "PurgMagSteppingAction.hh"
 40 #include "PurgMagRunAction.hh"                     42 #include "PurgMagRunAction.hh"
 41 #include "PurgMagDetectorConstruction.hh"          43 #include "PurgMagDetectorConstruction.hh"
 42                                                    44 
 43 #include "G4SteppingManager.hh"                    45 #include "G4SteppingManager.hh"
 44 #include "G4Electron.hh"                       << 
 45 #include "G4Gamma.hh"                          << 
 46 #include "G4Positron.hh"                       << 
 47 #include "G4VTouchable.hh"                         46 #include "G4VTouchable.hh"
 48 #include "G4VPhysicalVolume.hh"                    47 #include "G4VPhysicalVolume.hh"
 49 #include "G4AnalysisManager.hh"                << 
 50                                                    48 
 51 #include "G4SystemOfUnits.hh"                  <<  49 #ifdef G4ANALYSIS_USE
                                                   >>  50 #include"PurgMagAnalysisManager.hh"
                                                   >>  51 #endif
 52                                                    52 
 53 //....oooOO0OOooo........oooOO0OOooo........oo     53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 54                                                    54 
 55 PurgMagSteppingAction::PurgMagSteppingAction(c <<  55 PurgMagSteppingAction::PurgMagSteppingAction(PurgMagRunAction* run,PurgMagDetectorConstruction* det)
 56                PurgMagDetectorConstruction* de <<  56 :PurgMagRun(run),Detector(det)
 57 :Detector(det)                                 << 
 58 { }                                                57 { }
 59                                                    58 
 60 //....oooOO0OOooo........oooOO0OOooo........oo     59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 61                                                    60 
 62 PurgMagSteppingAction::~PurgMagSteppingAction(     61 PurgMagSteppingAction::~PurgMagSteppingAction()
 63 { }                                                62 { }
 64                                                    63 
 65 //....oooOO0OOooo........oooOO0OOooo........oo     64 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 66                                                    65 
 67 void PurgMagSteppingAction::UserSteppingAction     66 void PurgMagSteppingAction::UserSteppingAction(const G4Step* aStep)
 68                                                    67   
 69 {                                                  68 { 
 70   G4AnalysisManager* man = G4AnalysisManager:: <<  69 #ifdef G4ANALYSIS_USE
 71                                                << 
 72   //Collection at SSD in N-tuples. Electrons a     70   //Collection at SSD in N-tuples. Electrons and photons separated
 73   //Prestep point in World, next volume Measur     71   //Prestep point in World, next volume MeasureVolume, process transportation
                                                   >>  72  PurgMagAnalysisManager* analysis = PurgMagAnalysisManager::getInstance();  
 74   if ((aStep->GetPreStepPoint()->GetPhysicalVo     73   if ((aStep->GetPreStepPoint()->GetPhysicalVolume() == Detector->GetWorld())&&
 75       (aStep->GetTrack()->GetNextVolume() == D     74       (aStep->GetTrack()->GetNextVolume() == Detector->GetMeasureVolume())&&
 76       //(aStep->GetTrack()->GetMomentumDirecti     75       //(aStep->GetTrack()->GetMomentumDirection().z()>0.)&& // only particles with positive momentum
 77       (aStep->GetPostStepPoint()->GetProcessDe     76       (aStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName() == "Transportation"))
 78     {                                              77     {
 79       G4double gx, gy, gz, ge, gpx, gpy, gpz,      78       G4double gx, gy, gz, ge, gpx, gpy, gpz, ex, ey, ez, ee; 
 80       G4double epx, epy, epz, px, py, pz, pe,  <<  79       G4double epx, epy, epz, px, py, pz, pe, ppx, ppy, ppz;
                                                   >>  80    
 81                                                    81      
 82       // Electrons                                 82       // Electrons 
 83       if(aStep->GetTrack()->GetDynamicParticle <<  83       if(aStep->GetTrack()->GetDynamicParticle()->GetDefinition()->GetParticleName()
 84    == G4Electron::Definition())                <<  84       == "e-")
 85   {//Position                                      85   {//Position 
 86          ex = (aStep->GetTrack()->GetPosition(     86          ex = (aStep->GetTrack()->GetPosition().x())/cm;
 87          ey = (aStep->GetTrack()->GetPosition(     87          ey = (aStep->GetTrack()->GetPosition().y())/cm; 
 88          ez = (aStep->GetTrack()->GetPosition(     88          ez = (aStep->GetTrack()->GetPosition().z())/cm;
 89    // Energy                                       89    // Energy
 90    ee = (aStep->GetTrack()->GetKineticEnergy()     90    ee = (aStep->GetTrack()->GetKineticEnergy())/MeV;
 91    // Momentum                                     91    // Momentum
 92    epx = aStep->GetTrack()->GetMomentum().x();     92    epx = aStep->GetTrack()->GetMomentum().x();
 93    epy = aStep->GetTrack()->GetMomentum().y();     93    epy = aStep->GetTrack()->GetMomentum().y();
 94    epz = aStep->GetTrack()->GetMomentum().z();     94    epz = aStep->GetTrack()->GetMomentum().z();
 95                                                    95     
 96    // Fill N-tuple electrons  ( id = 1)        <<  96    // Fill N-tuple electrons
 97    man->FillNtupleDColumn(1,0, ex);            <<  97    analysis->fill_Tuple_Electrons(ex, ey, ez, ee, epx, epy, epz);
 98    man->FillNtupleDColumn(1,1, ey);            << 
 99    man->FillNtupleDColumn(1,2, ez);            << 
100    man->FillNtupleDColumn(1,3, ee);            << 
101    man->FillNtupleDColumn(1,4, epx);           << 
102    man->FillNtupleDColumn(1,5, epy);           << 
103    man->FillNtupleDColumn(1,6, epz);           << 
104    man->AddNtupleRow(1);                       << 
105       }                                            98       }
106                                                    99 
107       // Photons                                  100       // Photons
108       if (aStep->GetTrack()->GetDynamicParticl << 101       if (aStep->GetTrack()->GetDynamicParticle()->GetDefinition()->
109     G4Gamma::Definition())                     << 102           GetParticleName() == "gamma")
110   {                                               103   {
111                                                   104     
112     // Position                                   105     // Position
113           gx = (aStep->GetTrack()->GetPosition    106           gx = (aStep->GetTrack()->GetPosition().x())/cm;
114     gy = (aStep->GetTrack()->GetPosition().y()    107     gy = (aStep->GetTrack()->GetPosition().y())/cm;
115     gz = (aStep->GetTrack()->GetPosition().z()    108     gz = (aStep->GetTrack()->GetPosition().z())/cm;
116                                                   109     
117     // Energy                                     110     // Energy
118     ge = (aStep->GetTrack()->GetKineticEnergy(    111     ge = (aStep->GetTrack()->GetKineticEnergy())/MeV;
119                                                   112     
120     // Momentum                                   113     // Momentum
121     gpx = aStep->GetTrack()->GetMomentum().x()    114     gpx = aStep->GetTrack()->GetMomentum().x();
122     gpy = aStep->GetTrack()->GetMomentum().y()    115     gpy = aStep->GetTrack()->GetMomentum().y();
123     gpz = aStep->GetTrack()->GetMomentum().z()    116     gpz = aStep->GetTrack()->GetMomentum().z();
124                                                   117 
125     // Fill N-tuple photons (id=2)             << 118     // Fill N-tuple photons
126     man->FillNtupleDColumn(2,0, gx);           << 119     analysis->fill_Tuple_Gamma(gx, gy, gz, ge, gpx, gpy, gpz);
127     man->FillNtupleDColumn(2,1, gy);           << 
128     man->FillNtupleDColumn(2,2, gz);           << 
129     man->FillNtupleDColumn(2,3, ge);           << 
130     man->FillNtupleDColumn(2,4, gpx);          << 
131     man->FillNtupleDColumn(2,5, gpy);          << 
132     man->FillNtupleDColumn(2,6, gpz);          << 
133     man->AddNtupleRow(2);                      << 
134   }                                               120   }
135                                                   121 
136                                                   122 
137       // Positrons                                123       // Positrons
138       if (aStep->GetTrack()->GetDynamicParticl << 124       if (aStep->GetTrack()->GetDynamicParticle()->GetDefinition()->GetParticleName() == "e+")
139     G4Positron::Definition())                  << 
140   {                                               125   {
141                                                   126 
142     // Position                                   127     // Position
143           px = (aStep->GetTrack()->GetPosition    128           px = (aStep->GetTrack()->GetPosition().x())/cm;
144     py = (aStep->GetTrack()->GetPosition().y()    129     py = (aStep->GetTrack()->GetPosition().y())/cm;
145     pz = (aStep->GetTrack()->GetPosition().z()    130     pz = (aStep->GetTrack()->GetPosition().z())/cm;
146                                                   131     
147     // Energy                                     132     // Energy
148     pe = (aStep->GetTrack()->GetKineticEnergy(    133     pe = (aStep->GetTrack()->GetKineticEnergy())/MeV;
149                                                   134     
150     // Momentum                                   135     // Momentum
151     ppx = aStep->GetTrack()->GetMomentum().x()    136     ppx = aStep->GetTrack()->GetMomentum().x();
152     ppy = aStep->GetTrack()->GetMomentum().y()    137     ppy = aStep->GetTrack()->GetMomentum().y();
153     ppz = aStep->GetTrack()->GetMomentum().z()    138     ppz = aStep->GetTrack()->GetMomentum().z();
154                                                   139     
155     // Fill Ntuple positrons ( id = 3)         << 140     // Fill Ntuple positrons
156     man->FillNtupleDColumn(3,0, px);           << 141     analysis->fill_Tuple_Positrons(px, py, pz, pe, ppx, ppy, ppz);
157     man->FillNtupleDColumn(3,1, py);           << 
158     man->FillNtupleDColumn(3,2, pz);           << 
159     man->FillNtupleDColumn(3,3, pe);           << 
160     man->FillNtupleDColumn(3,4, ppx);          << 
161     man->FillNtupleDColumn(3,5, ppy);          << 
162     man->FillNtupleDColumn(3,6, ppz);          << 
163     man->AddNtupleRow(3);                      << 
164   }                                               142   }
165     }                                             143     }
                                                   >> 144 #endif
166 }                                                 145 }
167 //....oooOO0OOooo........oooOO0OOooo........oo    146 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
168                                                   147