Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/tracking/src/G4RichTrajectoryPoint.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 /tracking/src/G4RichTrajectoryPoint.cc (Version 11.3.0) and /tracking/src/G4RichTrajectoryPoint.cc (Version 9.5)


  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 // G4RichTrajectoryPoint class implementation  <<  26 //
                                                   >>  27 // $Id: G4RichTrajectoryPoint.cc,v 1.6 2010-02-22 21:26:56 allison Exp $
                                                   >>  28 // GEANT4 tag $Name: not supported by cvs2svn $
                                                   >>  29 //
                                                   >>  30 //
                                                   >>  31 // ---------------------------------------------------------------
                                                   >>  32 //
                                                   >>  33 // G4RichTrajectoryPoint.cc
 27 //                                                 34 //
 28 // Contact:                                        35 // Contact:
 29 //   Questions and comments on G4TrajectoryPoi     36 //   Questions and comments on G4TrajectoryPoint, on which this is based,
 30 //   should be sent to                             37 //   should be sent to
 31 //     Katsuya Amako  (e-mail: Katsuya.Amako@k     38 //     Katsuya Amako  (e-mail: Katsuya.Amako@kek.jp)
 32 //     Makoto  Asai   (e-mail: asai@slac.stanf <<  39 //     Makoto  Asai   (e-mail: asai@kekvax.kek.jp)
 33 //     Takashi Sasaki (e-mail: Takashi.Sasaki@     40 //     Takashi Sasaki (e-mail: Takashi.Sasaki@kek.jp)
 34 //   and on the extended code to:                  41 //   and on the extended code to:
 35 //     John Allison   (e-mail: John.Allison@ma     42 //     John Allison   (e-mail: John.Allison@manchester.ac.uk)
 36 //     Joseph Perl    (e-mail: perl@slac.stanf     43 //     Joseph Perl    (e-mail: perl@slac.stanford.edu)
 37 // ------------------------------------------- <<  44 //
                                                   >>  45 // ---------------------------------------------------------------
 38                                                    46 
 39 #include "G4RichTrajectoryPoint.hh"                47 #include "G4RichTrajectoryPoint.hh"
 40                                                    48 
 41 #include "G4AttDef.hh"                         <<  49 #include "G4Track.hh"
                                                   >>  50 #include "G4Step.hh"
                                                   >>  51 #include "G4VProcess.hh"
                                                   >>  52 
 42 #include "G4AttDefStore.hh"                        53 #include "G4AttDefStore.hh"
                                                   >>  54 #include "G4AttDef.hh"
 43 #include "G4AttValue.hh"                           55 #include "G4AttValue.hh"
 44 #include "G4Step.hh"                           << 
 45 #include "G4Track.hh"                          << 
 46 #include "G4UIcommand.hh"                      << 
 47 #include "G4UnitsTable.hh"                         56 #include "G4UnitsTable.hh"
 48 #include "G4VProcess.hh"                       << 
 49                                                    57 
 50 // #define G4ATTDEBUG                          <<  58 //#define G4ATTDEBUG
 51 #ifdef G4ATTDEBUG                                  59 #ifdef G4ATTDEBUG
 52 #  include "G4AttCheck.hh"                     <<  60 #include "G4AttCheck.hh"
 53 #endif                                             61 #endif
 54                                                    62 
 55 #include <sstream>                                 63 #include <sstream>
 56                                                    64 
 57 G4Allocator<G4RichTrajectoryPoint>*& aRichTraj <<  65 G4Allocator<G4RichTrajectoryPoint> aRichTrajectoryPointAllocator;
 58 {                                              << 
 59   G4ThreadLocalStatic G4Allocator<G4RichTrajec << 
 60   return _instance;                            << 
 61 }                                              << 
 62                                                    66 
 63 G4RichTrajectoryPoint::G4RichTrajectoryPoint() <<  67 G4RichTrajectoryPoint::G4RichTrajectoryPoint():
                                                   >>  68   fpAuxiliaryPointVector(0),
                                                   >>  69   fTotEDep(0.),
                                                   >>  70   fRemainingEnergy(0.),
                                                   >>  71   fpProcess(0),
                                                   >>  72   fPreStepPointStatus(fUndefined),
                                                   >>  73   fPostStepPointStatus(fUndefined),
                                                   >>  74   fPreStepPointGlobalTime(0),
                                                   >>  75   fPostStepPointGlobalTime(0),
                                                   >>  76   fPreStepPointWeight(1.),
                                                   >>  77   fPostStepPointWeight(1.)
                                                   >>  78 {}
 64                                                    79 
 65 G4RichTrajectoryPoint::G4RichTrajectoryPoint(c <<  80 G4RichTrajectoryPoint::G4RichTrajectoryPoint(const G4Track* aTrack):
 66   : fPosition(aTrack->GetPosition()),          <<  81   G4TrajectoryPoint(aTrack->GetPosition()),
 67     fPreStepPointGlobalTime(aTrack->GetGlobalT <<  82   fpAuxiliaryPointVector(0),
 68     fPostStepPointGlobalTime(aTrack->GetGlobal <<  83   fTotEDep(0.),
 69     fpPreStepPointVolume(aTrack->GetTouchableH <<  84   fRemainingEnergy(aTrack->GetKineticEnergy()),
 70     fpPostStepPointVolume(aTrack->GetNextTouch <<  85   fpProcess(0),
 71     fPreStepPointWeight(aTrack->GetWeight()),  <<  86   fPreStepPointStatus(fUndefined),
 72     fPostStepPointWeight(aTrack->GetWeight())  <<  87   fPostStepPointStatus(fUndefined),
                                                   >>  88   fPreStepPointGlobalTime(aTrack->GetGlobalTime()),
                                                   >>  89   fPostStepPointGlobalTime(aTrack->GetGlobalTime()),
                                                   >>  90   fpPreStepPointVolume(aTrack->GetTouchableHandle()),
                                                   >>  91   fpPostStepPointVolume(aTrack->GetNextTouchableHandle()),
                                                   >>  92   fPreStepPointWeight(aTrack->GetWeight()),
                                                   >>  93   fPostStepPointWeight(aTrack->GetWeight())
 73 {}                                                 94 {}
 74                                                    95 
 75 G4RichTrajectoryPoint::G4RichTrajectoryPoint(c <<  96 G4RichTrajectoryPoint::G4RichTrajectoryPoint(const G4Step* aStep):
 76   : fPosition(aStep->GetPostStepPoint()->GetPo <<  97   G4TrajectoryPoint(aStep->GetPostStepPoint()->GetPosition()),
 77     fpAuxiliaryPointVector(aStep->GetPointerTo <<  98   fpAuxiliaryPointVector(aStep->GetPointerToVectorOfAuxiliaryPoints()),
 78     fTotEDep(aStep->GetTotalEnergyDeposit())   <<  99   fTotEDep(aStep->GetTotalEnergyDeposit())
 79 {                                                 100 {
 80   G4StepPoint* preStepPoint = aStep->GetPreSte    101   G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
 81   G4StepPoint* postStepPoint = aStep->GetPostS    102   G4StepPoint* postStepPoint = aStep->GetPostStepPoint();
 82   if (aStep->GetTrack()->GetCurrentStepNumber( << 103   if (aStep->GetTrack()->GetCurrentStepNumber() <= 0) {  // First step
 83   {                                            << 
 84     fRemainingEnergy = aStep->GetTrack()->GetK    104     fRemainingEnergy = aStep->GetTrack()->GetKineticEnergy();
 85   }                                            << 105   } else {
 86   else {                                       << 
 87     fRemainingEnergy = preStepPoint->GetKineti    106     fRemainingEnergy = preStepPoint->GetKineticEnergy() - fTotEDep;
 88   }                                               107   }
 89   fpProcess = postStepPoint->GetProcessDefined    108   fpProcess = postStepPoint->GetProcessDefinedStep();
 90   fPreStepPointStatus = preStepPoint->GetStepS    109   fPreStepPointStatus = preStepPoint->GetStepStatus();
 91   fPostStepPointStatus = postStepPoint->GetSte    110   fPostStepPointStatus = postStepPoint->GetStepStatus();
 92   fPreStepPointGlobalTime = preStepPoint->GetG    111   fPreStepPointGlobalTime = preStepPoint->GetGlobalTime();
 93   fPostStepPointGlobalTime = postStepPoint->Ge    112   fPostStepPointGlobalTime = postStepPoint->GetGlobalTime();
 94   fpPreStepPointVolume = preStepPoint->GetTouc    113   fpPreStepPointVolume = preStepPoint->GetTouchableHandle();
 95   fpPostStepPointVolume = postStepPoint->GetTo    114   fpPostStepPointVolume = postStepPoint->GetTouchableHandle();
 96   fPreStepPointWeight = preStepPoint->GetWeigh    115   fPreStepPointWeight = preStepPoint->GetWeight();
 97   fPostStepPointWeight = postStepPoint->GetWei    116   fPostStepPointWeight = postStepPoint->GetWeight();
                                                   >> 117 
                                                   >> 118   /*
                                                   >> 119   G4cout << "fpAuxiliaryPointVector "
                                                   >> 120    << (void*) fpAuxiliaryPointVector;
                                                   >> 121   G4cout << ": ";
                                                   >> 122   if (fpAuxiliaryPointVector) {
                                                   >> 123     G4cout << "size: " << fpAuxiliaryPointVector->size();
                                                   >> 124     for (size_t i = 0; i < fpAuxiliaryPointVector->size(); ++i)
                                                   >> 125       G4cout << "\n  " << (*fpAuxiliaryPointVector)[i];
                                                   >> 126   } else {
                                                   >> 127     G4cout << "non-existent";
                                                   >> 128   }
                                                   >> 129   G4cout << G4endl;
                                                   >> 130 
                                                   >> 131   static const G4Step* lastStep = 0;
                                                   >> 132   if (aStep && aStep == lastStep) {
                                                   >> 133     G4cout << "********* aStep is same as last" << G4endl;
                                                   >> 134   }
                                                   >> 135   lastStep = aStep;
                                                   >> 136 
                                                   >> 137   static std::vector<G4ThreeVector>*  lastAuxiliaryPointVector = 0;
                                                   >> 138   if (fpAuxiliaryPointVector &&
                                                   >> 139       fpAuxiliaryPointVector == lastAuxiliaryPointVector) {
                                                   >> 140     G4cout << "********* fpAuxiliaryPointVector is same as last" << G4endl;
                                                   >> 141   }
                                                   >> 142   lastAuxiliaryPointVector = fpAuxiliaryPointVector;
                                                   >> 143   */
 98 }                                                 144 }
 99                                                   145 
100 // This is certainly wrong. the fpAuxiliaryPoi << 146 G4RichTrajectoryPoint::G4RichTrajectoryPoint
101 // and must be deep copied. Note also the copy << 147 (const G4RichTrajectoryPoint &right):
102 // However, G4RichTrajectory requires a copy c << 148   G4TrajectoryPoint(right),
103 G4RichTrajectoryPoint::G4RichTrajectoryPoint(c << 149   fpAuxiliaryPointVector(right.fpAuxiliaryPointVector),
                                                   >> 150   fTotEDep(right.fTotEDep),
                                                   >> 151   fRemainingEnergy(right.fRemainingEnergy),
                                                   >> 152   fpProcess(right.fpProcess),
                                                   >> 153   fPreStepPointStatus(right.fPreStepPointStatus),
                                                   >> 154   fPostStepPointStatus(right.fPostStepPointStatus),
                                                   >> 155   fPreStepPointGlobalTime(right.fPreStepPointGlobalTime),
                                                   >> 156   fPostStepPointGlobalTime(right.fPostStepPointGlobalTime),
                                                   >> 157   fpPreStepPointVolume(right.fpPreStepPointVolume),
                                                   >> 158   fpPostStepPointVolume(right.fpPostStepPointVolume),
                                                   >> 159   fPreStepPointWeight(right.fPreStepPointWeight),
                                                   >> 160   fPostStepPointWeight(right.fPostStepPointWeight)
                                                   >> 161 {}
104                                                   162 
105 G4RichTrajectoryPoint::~G4RichTrajectoryPoint( << 163 G4RichTrajectoryPoint::~G4RichTrajectoryPoint()
                                                   >> 164 {
                                                   >> 165   if(fpAuxiliaryPointVector) {
                                                   >> 166     /*
                                                   >> 167     G4cout << "Deleting fpAuxiliaryPointVector at "
                                                   >> 168      << (void*) fpAuxiliaryPointVector
                                                   >> 169      << G4endl;
                                                   >> 170     */
                                                   >> 171     delete fpAuxiliaryPointVector;
                                                   >> 172   }
                                                   >> 173 }
106                                                   174 
107 const std::map<G4String, G4AttDef>* G4RichTraj << 175 const std::map<G4String,G4AttDef>*
                                                   >> 176 G4RichTrajectoryPoint::GetAttDefs() const
108 {                                                 177 {
109   G4bool isNew;                                   178   G4bool isNew;
110   std::map<G4String, G4AttDef>* store = G4AttD << 179   std::map<G4String,G4AttDef>* store
                                                   >> 180     = G4AttDefStore::GetInstance("G4RichTrajectoryPoint",isNew);
111   if (isNew) {                                    181   if (isNew) {
                                                   >> 182 
                                                   >> 183     // Copy base class att defs...
                                                   >> 184     *store = *(G4TrajectoryPoint::GetAttDefs());
                                                   >> 185 
112     G4String ID;                                  186     G4String ID;
113                                                   187 
114     ID = "Pos";                                << 
115     (*store)[ID] = G4AttDef(ID, "Position", "P << 
116     ID = "Aux";                                   188     ID = "Aux";
117     (*store)[ID] =                             << 189     (*store)[ID] = G4AttDef(ID, "Auxiliary Point Position",
118       G4AttDef(ID, "Auxiliary Point Position", << 190           "Physics","G4BestUnit","G4ThreeVector");
119     ID = "TED";                                   191     ID = "TED";
120     (*store)[ID] = G4AttDef(ID, "Total Energy  << 192     (*store)[ID] = G4AttDef(ID,"Total Energy Deposit",
                                                   >> 193           "Physics","G4BestUnit","G4double");
121     ID = "RE";                                    194     ID = "RE";
122     (*store)[ID] = G4AttDef(ID, "Remaining Ene << 195     (*store)[ID] = G4AttDef(ID,"Remaining Energy",
                                                   >> 196           "Physics","G4BestUnit","G4double");
123     ID = "PDS";                                   197     ID = "PDS";
124     (*store)[ID] = G4AttDef(ID, "Process Defin << 198     (*store)[ID] = G4AttDef(ID,"Process Defined Step",
                                                   >> 199           "Physics","","G4String");
125     ID = "PTDS";                                  200     ID = "PTDS";
126     (*store)[ID] = G4AttDef(ID, "Process Type  << 201     (*store)[ID] = G4AttDef(ID,"Process Type Defined Step",
                                                   >> 202           "Physics","","G4String");
127     ID = "PreStatus";                             203     ID = "PreStatus";
128     (*store)[ID] = G4AttDef(ID, "Pre-step-poin << 204     (*store)[ID] = G4AttDef(ID,"Pre-step-point status",
                                                   >> 205           "Physics","","G4String");
129     ID = "PostStatus";                            206     ID = "PostStatus";
130     (*store)[ID] = G4AttDef(ID, "Post-step-poi << 207     (*store)[ID] = G4AttDef(ID,"Post-step-point status",
                                                   >> 208           "Physics","","G4String");
131     ID = "PreT";                                  209     ID = "PreT";
132     (*store)[ID] = G4AttDef(ID, "Pre-step-poin << 210     (*store)[ID] = G4AttDef(ID,"Pre-step-point global time",
                                                   >> 211           "Physics","G4BestUnit","G4double");
133     ID = "PostT";                                 212     ID = "PostT";
134     (*store)[ID] = G4AttDef(ID, "Post-step-poi << 213     (*store)[ID] = G4AttDef(ID,"Post-step-point global time",
                                                   >> 214           "Physics","G4BestUnit","G4double");
135     ID = "PreVPath";                              215     ID = "PreVPath";
136     (*store)[ID] = G4AttDef(ID, "Pre-step Volu << 216     (*store)[ID] = G4AttDef(ID,"Pre-step Volume Path",
                                                   >> 217                             "Physics","","G4String");
137     ID = "PostVPath";                             218     ID = "PostVPath";
138     (*store)[ID] = G4AttDef(ID, "Post-step Vol << 219     (*store)[ID] = G4AttDef(ID,"Post-step Volume Path",
                                                   >> 220                             "Physics","","G4String");
139     ID = "PreW";                                  221     ID = "PreW";
140     (*store)[ID] = G4AttDef(ID, "Pre-step-poin << 222     (*store)[ID] = G4AttDef(ID,"Pre-step-point weight",
                                                   >> 223           "Physics","","G4double");
141     ID = "PostW";                                 224     ID = "PostW";
142     (*store)[ID] = G4AttDef(ID, "Post-step-poi << 225     (*store)[ID] = G4AttDef(ID,"Post-step-point weight",
                                                   >> 226           "Physics","","G4double");
143   }                                               227   }
144   return store;                                   228   return store;
145 }                                                 229 }
146                                                   230 
147 static G4String Status(G4StepStatus stps)      << 231 static G4String Status(G4StepStatus s)
148 {                                                 232 {
149   G4String status;                                233   G4String status;
150   switch (stps) {                              << 234   switch (s) {
151     case fWorldBoundary:                       << 235   case fWorldBoundary:         status = "fWorldBoundary"; break;
152       status = "fWorldBoundary";               << 236   case fGeomBoundary:          status = "fGeomBoundary"; break;
153       break;                                   << 237   case fAtRestDoItProc:        status = "fAtRestDoItProc"; break;
154     case fGeomBoundary:                        << 238   case fAlongStepDoItProc:     status = "fAlongStepDoItProc"; break;
155       status = "fGeomBoundary";                << 239   case fPostStepDoItProc:      status = "fPostStepDoItProc"; break;
156       break;                                   << 240   case fUserDefinedLimit:      status = "fUserDefinedLimit"; break;
157     case fAtRestDoItProc:                      << 241   case fExclusivelyForcedProc: status = "fExclusivelyForcedProc"; break;
158       status = "fAtRestDoItProc";              << 242   case fUndefined:             status = "fUndefined"; break;
159       break;                                   << 243   default:                     status = "Not recognised"; break;
160     case fAlongStepDoItProc:                   << 
161       status = "fAlongStepDoItProc";           << 
162       break;                                   << 
163     case fPostStepDoItProc:                    << 
164       status = "fPostStepDoItProc";            << 
165       break;                                   << 
166     case fUserDefinedLimit:                    << 
167       status = "fUserDefinedLimit";            << 
168       break;                                   << 
169     case fExclusivelyForcedProc:               << 
170       status = "fExclusivelyForcedProc";       << 
171       break;                                   << 
172     case fUndefined:                           << 
173       status = "fUndefined";                   << 
174       break;                                   << 
175     default:                                   << 
176       status = "Not recognised";               << 
177       break;                                   << 
178   }                                               244   }
179   return status;                                  245   return status;
180 }                                                 246 }
181                                                   247 
182 static G4String Path(const G4TouchableHandle&     248 static G4String Path(const G4TouchableHandle& th)
183 {                                                 249 {
184   std::ostringstream oss;                         250   std::ostringstream oss;
185   G4int depth = th->GetHistoryDepth();            251   G4int depth = th->GetHistoryDepth();
186   for (G4int i = depth; i >= 0; --i) {            252   for (G4int i = depth; i >= 0; --i) {
187     oss << th->GetVolume(i)->GetName() << ':'  << 253     oss << th->GetVolume(i)->GetName()
                                                   >> 254   << ':' << th->GetCopyNumber(i);
188     if (i != 0) oss << '/';                       255     if (i != 0) oss << '/';
189   }                                               256   }
190   return oss.str();                               257   return oss.str();
191 }                                                 258 }
192                                                   259 
193 std::vector<G4AttValue>* G4RichTrajectoryPoint    260 std::vector<G4AttValue>* G4RichTrajectoryPoint::CreateAttValues() const
194 {                                                 261 {
195   // Create base class att values...              262   // Create base class att values...
                                                   >> 263   std::vector<G4AttValue>* values = G4TrajectoryPoint::CreateAttValues();
196                                                   264 
197   auto values = new std::vector<G4AttValue>;   << 265   if (fpAuxiliaryPointVector) {
198   values->push_back(G4AttValue("Pos", G4BestUn << 266     std::vector<G4ThreeVector>::iterator iAux;
199                                                << 267     for (iAux = fpAuxiliaryPointVector->begin();
200 #ifdef G4ATTDEBUG                              << 268    iAux != fpAuxiliaryPointVector->end(); ++iAux) {
201   G4cout << G4AttCheck(values, GetAttDefs());  << 269       values->push_back(G4AttValue("Aux",G4BestUnit(*iAux,"Length"),""));
202 #endif                                         << 
203                                                << 
204   if (fpAuxiliaryPointVector != nullptr) {     << 
205     for (const auto& iAux : *fpAuxiliaryPointV << 
206       values->push_back(G4AttValue("Aux", G4Be << 
207     }                                             270     }
208   }                                               271   }
209                                                   272 
210   values->push_back(G4AttValue("TED", G4BestUn << 273   values->push_back(G4AttValue("TED",G4BestUnit(fTotEDep,"Energy"),""));
211   values->push_back(G4AttValue("RE", G4BestUni << 
212                                                   274 
213   if (fpProcess != nullptr) {                  << 275   values->push_back(G4AttValue("RE",G4BestUnit(fRemainingEnergy,"Energy"),""));
214     values->push_back(G4AttValue("PDS", fpProc << 276 
215     values->push_back(                         << 277   if (fpProcess) {
216       G4AttValue("PTDS", G4VProcess::GetProces << 278     values->push_back
217   }                                            << 279       (G4AttValue("PDS",fpProcess->GetProcessName(),""));
218   else {                                       << 280     values->push_back
219     values->push_back(G4AttValue("PDS", "None" << 281       (G4AttValue
220     values->push_back(G4AttValue("PTDS", "None << 282        ("PTDS",G4VProcess::GetProcessTypeName(fpProcess->GetProcessType()),
                                                   >> 283   ""));
                                                   >> 284   } else {
                                                   >> 285     values->push_back(G4AttValue("PDS","None",""));
                                                   >> 286     values->push_back(G4AttValue("PTDS","None",""));
221   }                                               287   }
222                                                   288 
223   values->push_back(G4AttValue("PreStatus", St << 289   values->push_back
                                                   >> 290     (G4AttValue("PreStatus",Status(fPreStepPointStatus),""));
224                                                   291 
225   values->push_back(G4AttValue("PostStatus", S << 292   values->push_back
                                                   >> 293     (G4AttValue("PostStatus",Status(fPostStepPointStatus),""));
226                                                   294 
227   values->push_back(G4AttValue("PreT", G4BestU << 295   values->push_back
                                                   >> 296     (G4AttValue("PreT",G4BestUnit(fPreStepPointGlobalTime,"Time"),""));
228                                                   297 
229   values->push_back(G4AttValue("PostT", G4Best << 298   values->push_back
                                                   >> 299     (G4AttValue("PostT",G4BestUnit(fPostStepPointGlobalTime,"Time"),""));
230                                                   300 
231   if (fpPreStepPointVolume && (fpPreStepPointV << 301   if (fpPreStepPointVolume && fpPreStepPointVolume->GetVolume()) {
232     values->push_back(G4AttValue("PreVPath", P << 302     values->push_back(G4AttValue("PreVPath",Path(fpPreStepPointVolume),""));
233   }                                            << 303   } else {
234   else {                                       << 304     values->push_back(G4AttValue("PreVPath","None",""));
235     values->push_back(G4AttValue("PreVPath", " << 
236   }                                               305   }
237                                                   306 
238   if (fpPostStepPointVolume && (fpPostStepPoin << 307   if (fpPostStepPointVolume && fpPostStepPointVolume->GetVolume()) {
239     values->push_back(G4AttValue("PostVPath",  << 308     values->push_back(G4AttValue("PostVPath",Path(fpPostStepPointVolume),""));
240   }                                            << 309   } else {
241   else {                                       << 310     values->push_back(G4AttValue("PostVPath","None",""));
242     values->push_back(G4AttValue("PostVPath",  << 
243   }                                               311   }
244                                                   312 
245   std::ostringstream oss1;                     << 313   {
246   oss1 << fPreStepPointWeight;                 << 314     std::ostringstream oss;
247   values->push_back(G4AttValue("PreW", oss1.st << 315     oss << fPreStepPointWeight;
                                                   >> 316     values->push_back
                                                   >> 317       (G4AttValue("PreW",oss.str(),""));
                                                   >> 318   }
248                                                   319 
249   std::ostringstream oss2;                     << 320   {
250   oss2 << fPostStepPointWeight;                << 321     std::ostringstream oss;
251   values->push_back(G4AttValue("PostW", oss2.s << 322     oss << fPostStepPointWeight;
                                                   >> 323     values->push_back
                                                   >> 324       (G4AttValue("PostW",oss.str(),""));
                                                   >> 325   }
252                                                   326 
253 #ifdef G4ATTDEBUG                                 327 #ifdef G4ATTDEBUG
254   G4cout << G4AttCheck(values, GetAttDefs());  << 328   G4cout << G4AttCheck(values,GetAttDefs());
255 #endif                                            329 #endif
256                                                   330 
257   return values;                                  331   return values;
258 }                                                 332 }
259                                                   333