Geant4  10.01.p02
ML2Convergence.cc
Go to the documentation of this file.
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 // The code was written by :
27 // ^Claudio Andenna claudio.andenna@ispesl.it, claudio.andenna@iss.infn.it
28 // *Barbara Caccia barbara.caccia@iss.it
29 // with the support of Pablo Cirrone (LNS, INFN Catania Italy)
30 // with the contribute of Alessandro Occhigrossi*
31 //
32 // ^INAIL DIPIA - ex ISPESL and INFN Roma, gruppo collegato Sanità, Italy
33 // *Istituto Superiore di Sanità and INFN Roma, gruppo collegato Sanità, Italy
34 // Viale Regina Elena 299, 00161 Roma (Italy)
35 // tel (39) 06 49902246
36 // fax (39) 06 49387075
37 //
38 // more information:
39 // http://g4advancedexamples.lngs.infn.it/Examples/medical-linac
40 //
41 //*******************************************************//
42 
43 
44 #include "ML2Convergence.hh"
45 
47 {}
48 
50  G4String FileExperimentalData, G4String FileExperimentalDataOut,
51  G4bool bComp, G4int maxNumEvents, G4int nRecycling, G4int maxLoops)
52  :ML2ExpVoxels(0)
53 {
54  nGeometry=0;
55  nMaxLoops=maxLoops;
57  bCompareExp=bComp;
59  if (bCompareExp){nMaxLoops=-1;};
60  fileExperimentalData=FileExperimentalData;
61 
62  // if the flag compareExp if true and the experimental data is given create the class CML2ExpVoxels
64  {
65  ML2ExpVoxels=new CML2ExpVoxels(bCompareExp, saveEvents, seed, FileExperimentalData, FileExperimentalDataOut);
66  if (!ML2ExpVoxels->loadData())
67  {
68  ML2ExpVoxels=0;
69  std::cout <<"I don't have any convergence criteria set, I'll do " << nMaxLoops << " loop(s) for each rotation"<< G4endl;
70  }
71  else
72  {
73  ML2ExpVoxels->setRecycling(nRecycling);
74  }
75  }
76  maxNumberOfEvents=maxNumEvents;
77 }
78 
80 {
81  if (ML2ExpVoxels!=0)
82  {delete ML2ExpVoxels;}
83 
84 }
85 void CML2Convergence::add(const G4Step* aStep)
86 {
87  // accumulate events in the CML2ExpVoxels class (if created)
88  if (ML2ExpVoxels!=0)
89  {
90  if (aStep->GetTotalEnergyDeposit()>0.)
91  {ML2ExpVoxels->add(aStep);}
92  }
93 }
95 {
96  G4bool bStopRun=false;
97  if (ML2ExpVoxels!=0) // true if the experimental data file exists and is used to check the convergence
98  {
99  bStopRun=convergenceCriteria();
100  return bStopRun;
101  }
102  else // true if no experiemental data file is used. In this case it runs "nMaxLoops" loops.
103  {
104  idCurrentLoop--;
105  if (idCurrentLoop==0)
106  {
107  bStopRun=true;
108  }
109  }
110  return bStopRun;
111 }
113 {
114  G4bool bStopRun=true;
115  G4int nEventsAccumulated=0;
116  if (bCompareExp)
117  {
118  nEventsAccumulated=ML2ExpVoxels->getMaxNumberOfEvents();
119  // It checks if the maximum number of events is reached at least in one voxel. Having more rotations the limits is incremented each rotation
121  {bStopRun = true; ML2ExpVoxels->resetNEventsInVoxels();}
122  else
123  {bStopRun = false;}
124  }
125  std::cout <<"\n ++++++++++++++++++++ \n";
126  std::cout <<"current geometry: " << nGeometry;
127  std::cout << "\nNumber of events accumulated in the current geometry:"<<
128  nEventsAccumulated<<"\nNumber of events to be accumulated:" <<
129  maxNumberOfEvents<< "\n -------------------------\n";
130  return bStopRun;
131 }
G4String fileExperimentalData
G4bool loadData()
Definition: ML2ExpVoxels.cc:77
void add(const G4Step *aStep)
int G4int
Definition: G4Types.hh:78
CML2ExpVoxels * ML2ExpVoxels
void resetNEventsInVoxels()
bool G4bool
Definition: G4Types.hh:79
Definition: G4Step.hh:76
G4bool convergenceCriteria()
G4double GetTotalEnergyDeposit() const
G4int getMaxNumberOfEvents()
void setRecycling(int recycling)
Definition: ML2ExpVoxels.hh:66
#define G4endl
Definition: G4ios.hh:61
void add(G4ThreeVector pos, G4double depEnergy, G4double density)