add inlude files

include/ActionInitialization.hh

+//
+// ********************************************************************
+// * License and Disclaimer                                           *
+// *                                                                  *
+// * The  Geant4 software  is  copyright of the Copyright Holders  of *
+// * the Geant4 Collaboration.  It is provided  under  the terms  and *
+// * conditions of the Geant4 Software License,  included in the file *
+// * LICENSE and available at  http://cern.ch/geant4/license .  These *
+// * include a list of copyright holders.                             *
+// *                                                                  *
+// * Neither the authors of this software system, nor their employing *
+// * institutes,nor the agencies providing financial support for this *
+// * work  make  any representation or  warranty, express or implied, *
+// * regarding  this  software system or assume any liability for its *
+// * use.  Please see the license in the file  LICENSE  and URL above *
+// * for the full disclaimer and the limitation of liability.         *
+// *                                                                  *
+// * This  code  implementation is the result of  the  scientific and *
+// * technical work of the GEANT4 collaboration.                      *
+// * By using,  copying,  modifying or  distributing the software (or *
+// * any work based  on the software)  you  agree  to acknowledge its *
+// * use  in  resulting  scientific  publications,  and indicate your *
+// * acceptance of all terms of the Geant4 Software license.          *
+// ********************************************************************
+//
+// $Id: ActionInitialization.hh 76485 2013-11-11 10:54:47Z gcosmo $
+//
+/// \file ActionInitialization.hh
+/// \brief Definition of the ActionInitialization class
+
+#ifndef ActionInitialization_h
+#define ActionInitialization_h 1
+
+#include "G4VUserActionInitialization.hh"
+
+class DetectorConstruction;
+
+/// Action initialization class.
+///
+
+class ActionInitialization : public G4VUserActionInitialization
+{
+  public:
+    ActionInitialization(DetectorConstruction*);
+    //ActionInitialization();
+    virtual ~ActionInitialization();
+
+    virtual void BuildForMaster() const;
+    virtual void Build() const;
+    
+    
+  private:
+    DetectorConstruction* fDetector;
+
+};
+
+#endif

include/CollimatorConstruction.hh

+#ifndef COLLIMATORCONSTRUCTION_H
+#define COLLIMATORCONSTRUCTION_H
+
+#include "globals.hh"
+#include "G4VUserDetectorConstruction.hh"
+#include "G4SystemOfUnits.hh"
+#include "G4LogicalVolume.hh"
+
+class CollimatorConstruction
+{
+  public:
+    G4VSolid* ConstructCollimator(G4double c_thick, G4double c_coating_thick);
+    G4VSolid* ConstructCollimatorCoating(G4double c_thick, G4double c_coating_thick);
+    G4VSolid* ConstructSmallCollimator(G4double c_thick, G4double c_coating_thick);
+    G4VSolid* ConstructSAMBACollimator(G4double c_thick);
+
+    G4double GetCollimatorSizeX()  {return halfCollimatorSizeX;};
+    G4double GetCollimatorSizeY()  {return halfCollimatorSizeY;};
+    G4double GetCollimatorSizeZ()  {return halfCollimatorSizeZ;};
+
+  private:
+
+    //Initial block dimensions
+    G4double halfCollimatorSizeX ;
+    G4double halfCollimatorSizeY ;
+    G4double halfCollimatorSizeZ ;
+
+};
+#endif

include/DetectorConstruction.hh

+//
+// ********************************************************************
+// * License and Disclaimer                                           *
+// *                                                                  *
+// * The  Geant4 software  is  copyright of the Copyright Holders  of *
+// * the Geant4 Collaboration.  It is provided  under  the terms  and *
+// * conditions of the Geant4 Software License,  included in the file *
+// * LICENSE and available at  http://cern.ch/geant4/license .  These *
+// * include a list of copyright holders.                             *
+// *                                                                  *
+// * Neither the authors of this software system, nor their employing *
+// * institutes,nor the agencies providing financial support for this *
+// * work  make  any representation or  warranty, express or implied, *
+// * regarding  this  software system or assume any liability for its *
+// * use.  Please see the license in the file  LICENSE  and URL above *
+// * for the full disclaimer and the limitation of liability.         *
+// *                                                                  *
+// * This  code  implementation is the result of  the  scientific and *
+// * technical work of the GEANT4 collaboration.                      *
+// * By using,  copying,  modifying or  distributing the software (or *
+// * any work based  on the software)  you  agree  to acknowledge its *
+// * use  in  resulting  scientific  publications,  and indicate your *
+// * acceptance of all terms of the Geant4 Software license.          *
+// ********************************************************************
+//
+/// \file OpNovice/include/OpNoviceDetectorConstruction.hh
+/// \brief Definition of the OpNoviceDetectorConstruction class
+//
+//
+//
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+#ifndef DetectorConstruction_h
+#define DetectorConstruction_h 1
+
+#include "globals.hh"
+#include "G4VUserDetectorConstruction.hh"
+#include "DetectorMessenger.hh"
+#include "G4SystemOfUnits.hh"
+#include "G4MaterialPropertiesTable.hh"
+#include "SoleilMaterials.hh"
+#include "G4GDMLParser.hh"
+class DetectorMessenger;
+class G4LogicalVolume;
+class G4Material;
+class SoleilMaterials;
+class G4Box;
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+class DetectorConstruction : public G4VUserDetectorConstruction
+{
+  public:
+    DetectorConstruction();
+    virtual ~DetectorConstruction();
+    static G4VPhysicalVolume *
+    GetPhysicalVolumeByName(const G4String &name);
+
+  public:
+    virtual G4VPhysicalVolume* Construct();
+    void UpdateGeometry();
+    //static G4VPhysicalVolume* GetPhysicalVolumeByName(const G4String &name);
+    void SetSize  (G4double);
+    void SetCollimatorMaterial(G4String);
+    void SetTargetMaterial(G4String);
+    void SetGeContainerMaterial(G4String);
+    void SetGeContainerMaterialCoating(G4String);
+    void SetSampleMaterial(G4String);
+    void SetOutputDirectory(G4String);
+    void SetWorldMaterial(G4String);
+    void SetDetectorType(G4int);
+    void SetGeDetectorLength(G4double);
+    void SetGeDetectorThickness(G4double);
+    void SetGeDetectorWidth(G4double);
+    void SetContactThickness(G4double);
+    void SetSampleLength(G4double);
+    void SetSampleThickness(G4double);
+    void SetSampleWidth(G4double);
+    void SetBeWindowRadius(G4double);
+    void SetDetectorCollimatorX(G4double);
+    void SetCollimatorThickness(G4double);
+    void SetDistanceCollimatorDetector(G4double);
+    void SetNumberOfTargetSamples(G4int);
+    void SetVolName(G4ThreeVector);
+    void SetDetectorName(G4String);
+    void SetSetupName(G4String);
+    void SetDataType(G4String);
+
+
+
+  public:
+    //const G4VPhysicalVolume* GetWorld() {return physicPenSampleBox;};
+    const G4VPhysicalVolume* GetWorld() {return physicWorldBox;};
+    G4int GetDetectorType(){return fDetectorType;};
+    G4double GetDetectorCollimatorX(){return fDetectorCollimatorX;};
+    G4double GetDetectorCollimatorY(){return fDetectorCollimatorY;};
+    G4double GetSourceContainerY(){return fSourceContainerY;};
+    G4String GetDetectorName(){return fDetectorName;};
+    G4String GetDetectorOutputDataDirectory(){return data_output_directory;};
+    G4String GetVolName(){return fVolName;};
+
+    G4double GetCollimatorThickness()  {return CollimatorThickness;};
+    G4double GetContactThickness()  {return ContactThickness;};
+    G4double GetGeDetectorLength()  {return detectorInnerRadius*2.0;};
+    G4double GetGeDetectorThickness()  {return detectorThickness*2.0;};
+    G4double GetGeDetectorWidth()  {return detectorHeight*2.0;};
+    G4double GetSampleLength()  {return sampleInnerRadius*2.0;};
+    G4double GetSampleThickness()  {return sampleThickness*2.0;};
+    G4double GetSampleWidth()  {return sampleHeight*2.0;};
+    G4double GetBeWindowRadius()  {return beWindowRadius;};
+    G4double GetSamplePositionZ()  {return zPositionSample;};
+    G4double GetDistanceCollimatorDetector()  {return distanceCollimatorDetector;};
+    G4int GetNumberOfTargetSamples()  {return nSamples;};
+    G4String GetSetupName(){return fSetupName;};
+    G4String GetDataType(){return fDataType;};
+    G4Material*        GetWorldMaterial()   {return fWorldMaterial;};
+    G4Material*        GetTargetMaterial()   {return fTargetMaterial;};
+    G4Material*        GetCollimatorMaterial()   {return fCollimatorMaterial;};
+    G4ThreeVector* fSourceVector;
+
+    void               DefineMaterials();
+
+  private:
+
+    G4double halfSizeDarkBoxX;
+    G4double halfSizeDarkBoxY;
+    G4double halfSizeDarkBoxZ;
+
+    G4double detectorInnerRadius;
+    G4double detectorThickness;
+    G4double detectorHeight;
+    G4double CollimatorThickness;
+    G4double ContactThickness;
+    G4int nSamples;
+    G4double sampleInnerRadius;
+    G4double sampleThickness;
+    G4double sampleHeight;
+    G4double beWindowRadius;
+    G4double zPositionSample;
+
+
+    G4double fSiliconPlate_h;
+    G4double fHolderWidth;
+
+    G4Material* fWorldMaterial;
+    G4Material* fTargetMaterial;
+    G4Material* fCollimatorMaterial;
+    G4Material* fSampleMaterial;
+    G4String fSetupName;
+    G4double halfSourceContainerThickness;
+    G4double halfCollimatorThickness;
+    G4double distanceCollimatorDetector;
+    G4Box* fWorldBox;
+    G4Tubs* geDetectorTube;
+    G4Box* AlContactBox;
+    G4Tubs* sampleTube;
+
+    G4LogicalVolume* logicCollimatorLEAPSCoating;
+    G4LogicalVolume* logicCollimatorLEAPS;
+    G4LogicalVolume* logicCollimatorLEAPSSmall;
+    G4LogicalVolume* logicCollimatorSAMBA;
+    G4LogicalVolume* logicSample;
+    G4LogicalVolume* geLogicVolume;
+    G4LogicalVolume* logicContactVolume;
+    G4LogicalVolume* logicBeWindow;
+    G4LogicalVolume* logicMetalWindow;
+    G4LogicalVolume* logicMetalWindowBack;
+    G4LogicalVolume* logicMetalTube;
+    G4LogicalVolume* logicWorldBox;
+    G4LogicalVolume* logicGeContainer;
+    G4LogicalVolume* logicGeContainerCoating;
+    G4LogicalVolume* logicSupportPins;
+    G4LogicalVolume* logicGeCapot;
+
+    G4VSolid* solidBeWindow;
+    G4VSolid* solidMetalWindow;
+    G4VSolid* solidMetalWindowBack;
+    G4VSolid* solidMetalTube;
+    G4VSolid* solidGeContainer;
+    G4VSolid* solidGeContainerCoating;
+    G4VSolid* solidSupportPins;
+    G4VSolid* solidGeCapotPrimitive;
+    G4VSolid* solidGeCapot;
+    G4VSolid* solidBoxCapotInner;
+    
+
+    G4VPhysicalVolume* physicWorldBox;
+    G4VPhysicalVolume* physicPenSampleBox;
+    G4VPhysicalVolume* physicSourceContainer;
+    G4VPhysicalVolume* physicCollimator;
+
+
+
+    G4Material* PVTMaterial;
+    G4Material* PenMaterial;
+    G4Material* materialBialkali;
+    G4Material* materialSi;
+    G4Material* materialBeWindow;
+    G4Material* materialGe;
+    G4Material* materialAir;
+    G4Material* fVacuum;
+    G4Material* materialTriggerFoilEJ212;
+    G4Material* Pstyrene;
+    G4Material* fGlass;
+    G4Material* materialPMMA;
+    G4Material* materialGreaseEJ550;
+    G4Material* materialTeflon;
+    G4Material* materialVikuiti;
+    G4Material* materialPolyethylene;
+    G4Material* materialTitanium;
+    G4Material* materialCollimatorCoating;
+    G4Material* materialContacts;
+    G4Material* materialSample;
+    G4Material* materialMetalTube;
+    G4Material* materialSupportPins;
+    G4Material* materialGeContainer;
+    G4Material* materialGeContainerCoating;
+
+    G4int fDetectorType;
+    G4double fDetectorCollimatorX;
+    G4double fDetectorCollimatorY;
+    G4double fSourceContainerY;
+    G4String fDetectorName;
+    G4String fVolName;
+    G4String data_output_directory;
+    G4String fDataType;
+
+    SoleilMaterials* materialConstruction;
+    DetectorMessenger* fDetectorMessenger;
+};
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+#endif /*DetectorConstruction_h*/

include/DetectorMessenger.hh

+#ifndef DetectorMessenger_h
+#define DetectorMessenger_h 1
+
+#include "G4UImessenger.hh"
+#include "globals.hh"
+
+class DetectorConstruction;
+class G4UIdirectory;
+class G4UIcmdWithAString;
+class G4UIcmdWithADoubleAndUnit;
+class G4UIcmdWithADouble;
+class G4UIcmdWithoutParameter;
+class G4UIcmdWithAnInteger;
+class G4UIcmdWithABool;
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+class DetectorMessenger: public G4UImessenger
+{
+  public:
+
+    DetectorMessenger(DetectorConstruction* );
+   ~DetectorMessenger();
+
+    virtual void SetNewValue(G4UIcommand*, G4String);
+
+  private:
+
+    DetectorConstruction*      fDetector;
+
+
+    G4UIdirectory*             fGe2SOSDir;
+    G4UIdirectory*             fDetDir;
+
+    G4UIcmdWithAString*         commandSetWorldMaterial;
+    G4UIcmdWithAnInteger*       commandSetDetectorType;
+    G4UIcmdWithAnInteger*       commandSetNumberOfTargetSamples;
+    G4UIcmdWithADoubleAndUnit*  commandSetGeDetectorLength;
+    G4UIcmdWithADoubleAndUnit*  commandSetGeDetectorThickness;
+    G4UIcmdWithADoubleAndUnit*  commandSetGeDetectorWidth;
+    G4UIcmdWithADoubleAndUnit*  commandSetContactThickness;
+    G4UIcmdWithADoubleAndUnit*  commandSetBeWindowRadius;
+    G4UIcmdWithADoubleAndUnit*  commandSetCollimatorThickness;
+    G4UIcmdWithADoubleAndUnit*  commandSetDistanceCollimatorDetector;
+    G4UIcmdWithAString*         commandSetCollimatorMaterial;
+    G4UIcmdWithAString*         commandSetTargetMaterial;
+    G4UIcmdWithAString*         commandSetGeContainerMaterial;
+    G4UIcmdWithAString*         commandSetGeContainerMaterialCoating;
+    G4UIcmdWithAString*         commandSetDetectorName;
+    G4UIcmdWithAString*         commandSetSetupName;
+    G4UIcmdWithAString*         commandSetDataType;
+    G4UIcmdWithAString*         commandSetOutputDirectory;
+};
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+#endif

include/EventAction.hh

+#ifndef EventAction_h
+#define EventAction_h 1
+
+#include "G4UserEventAction.hh"
+#include "globals.hh"
+#include "G4ThreeVector.hh"
+
+class DetectorConstruction;
+class RunAction;
+
+/// Event action class
+///
+
+class EventAction : public G4UserEventAction
+{
+  public:
+    EventAction(DetectorConstruction* det, RunAction* runAction);
+    virtual ~EventAction();
+
+    virtual void BeginOfEventAction(const G4Event* event);
+    virtual void EndOfEventAction(const G4Event* event);
+    virtual void AddInfo(G4double xF, G4double yF, G4double zF, G4double depE);
+
+
+  private:
+
+    DetectorConstruction*   fDetector;
+    RunAction* 	fRunAction;
+
+    int64_t eventNumber;
+    G4double xFirstPen;
+    G4double yFirstPen;
+    G4double zFirstPen;
+    G4double depositedEnergyPENStackedSample1;
+     
+
+};
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+#endif

include/PhysicsList.hh

+//
+// ********************************************************************
+// * License and Disclaimer                                           *
+// *                                                                  *
+// * The  Geant4 software  is  copyright of the Copyright Holders  of *
+// * the Geant4 Collaboration.  It is provided  under  the terms  and *
+// * conditions of the Geant4 Software License,  included in the file *
+// * LICENSE and available at  http://cern.ch/geant4/license .  These *
+// * include a list of copyright holders.                             *
+// *                                                                  *
+// * Neither the authors of this software system, nor their employing *
+// * institutes,nor the agencies providing financial support for this *
+// * work  make  any representation or  warranty, express or implied, *
+// * regarding  this  software system or assume any liability for its *
+// * use.  Please see the license in the file  LICENSE  and URL above *
+// * for the full disclaimer and the limitation of liability.         *
+// *                                                                  *
+// * This  code  implementation is the result of  the  scientific and *
+// * technical work of the GEANT4 collaboration.                      *
+// * By using,  copying,  modifying or  distributing the software (or *
+// * any work based  on the software)  you  agree  to acknowledge its *
+// * use  in  resulting  scientific  publications,  and indicate your *
+// * acceptance of all terms of the Geant4 Software license.          *
+// ********************************************************************
+//
+/// \file OpNovice/include/OpNovicePhysicsList.hh
+/// \brief Definition of the OpNovicePhysicsList class
+//
+//
+//
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+#ifndef PhysicsList_h
+#define PhysicsList_h 1
+
+#include "globals.hh"
+#include "G4VUserPhysicsList.hh"
+
+class PhysicsListMessenger;
+
+class G4Cerenkov;
+class G4Scintillation;
+class G4OpAbsorption;
+class G4OpRayleigh;
+class G4OpMieHG;
+class G4OpBoundaryProcess;
+
+class G4VPhysicsConstructor;
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+class PhysicsList : public G4VUserPhysicsList
+{
+  public:
+
+    PhysicsList();
+    virtual ~PhysicsList();
+
+  public:
+
+    virtual void ConstructParticle();
+    virtual void ConstructProcess();
+
+    virtual void SetCuts();
+
+    //these methods Construct physics processes and register them
+    void ConstructDecay();
+    void ConstructEM();
+    void ConstructOp();
+
+    //for the Messenger
+    void SetVerbose(G4int);
+    void SetNbOfPhotonsCerenkov(G4int);
+
+  private:
+    //added by Luis
+    G4VPhysicsConstructor*  emPhysicsList;
+    G4int VerboseLevel;
+    G4int OpVerbLevel;
+
+    G4double cutForGamma;
+    G4double cutForElectron;
+    G4double cutForPositron;
+    //end of adds
+    static G4ThreadLocal G4int fVerboseLevel;
+    static G4ThreadLocal G4int fMaxNumPhotonStep;
+
+    static G4ThreadLocal G4Cerenkov* fCerenkovProcess;
+    static G4ThreadLocal G4Scintillation* fScintillationProcess;
+    static G4ThreadLocal G4OpAbsorption* fAbsorptionProcess;
+    static G4ThreadLocal G4OpRayleigh* fRayleighScatteringProcess;
+    static G4ThreadLocal G4OpMieHG* fMieHGScatteringProcess;
+    static G4ThreadLocal G4OpBoundaryProcess* fBoundaryProcess;
+};
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+#endif /* PhysicsList_h */

include/PrimaryGeneratorAction.hh

+#ifndef PrimaryGeneratorAction_h
+#define PrimaryGeneratorAction_h 1
+
+#include "G4VUserPrimaryGeneratorAction.hh"
+#include "PrimaryGeneratorMessenger.hh"
+#include "globals.hh"
+
+class G4ParticleGun;
+class G4Event;
+class PrimaryGeneratorMessenger;
+class DetectorConstruction;
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+class PrimaryGeneratorAction : public G4VUserPrimaryGeneratorAction
+{
+
+public:
+	PrimaryGeneratorAction(DetectorConstruction*);
+	~PrimaryGeneratorAction();
+
+// public:
+// 	PrimaryGeneratorAction();
+// 	virtual ~PrimaryGeneratorAction();
+
+public:
+	virtual void GeneratePrimaries(G4Event*);
+	void  SetSourceType(G4int newType);
+	void  SetSourceDirectionType(G4int newType);
+	void  SetSourceGeometry(G4int newType);
+	void  SetSourceEnergy(G4double newEnergy);
+	void  SetParticleName(G4int Z, G4int A, G4double excitEnergy);
+	void  SetSourcePositionX(G4double newValue);
+	void  SetSourcePositionY(G4double newValue);
+	void  SetSourcePositionZ(G4double newValue);
+	void  SetSourceDiameter(G4double newValue);
+	void  SetSourcePolarizationDegree(G4double newValue);
+	void  SetSourcePolarizationAngle(G4double newValue);
+	G4int GetSourceType(void){return fSourceType;};
+	G4int GetSourceGeometry(void){return fSourceGeometry;};
+	G4double GetSourceEnergy(void){return fSourceEnergy;};
+	G4String GetParticleName(void){return fParticleName;};
+	G4double GetSourcePositionX(void){return fPositionX;};
+	G4double GetSourcePositionY(void){return fPositionY;};
+	G4double GetSourcePositionZ(void){return fPositionZ;};
+	G4double GetSourceDiameter(void){return fDiameter;};
+	G4double GetSourcePolarizationAngle(void){return fTheta_polar;};
+	G4double GetSourcePolarizationDegree(void){return fPolarization_degree;};
+private:
+	G4double fPositionX;
+	G4double fPositionY;
+	G4double fPositionZ;
+	G4double fDiameter;
+	G4String fParticleName;
+	PrimaryGeneratorMessenger* fPrimaryMessenger;
+	G4ParticleGun*  fParticleGun;
+	G4int           fSourceType;
+	G4int           fSourceDirectionType;
+	G4int           fSourceGeometry;
+	G4double	fSourceEnergy;
+	G4double	fTheta_polar;
+	G4double	fPolarization_degree;
+	DetectorConstruction*      fDetector;
+	G4double		fPhotonWavelength;
+	PrimaryGeneratorMessenger* fGunMessenger;
+	G4ThreeVector fPoint;
+};
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+#endif /*PrimaryGeneratorAction_h*/

include/PrimaryGeneratorMessenger.hh

+#ifndef PrimaryGeneratorMessenger_h
+#define PrimaryGeneratorMessenger_h 1
+
+#include "G4UImessenger.hh"
+#include "globals.hh"
+
+class PrimaryGeneratorAction;
+class G4UIdirectory;
+class G4UIcmdWithADouble;
+class G4UIcmdWithADoubleAndUnit;
+class G4UIcmdWithAnInteger;
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+class PrimaryGeneratorMessenger: public G4UImessenger
+{
+  public:
+    PrimaryGeneratorMessenger(PrimaryGeneratorAction* );
+    virtual ~PrimaryGeneratorMessenger();
+
+    virtual void SetNewValue(G4UIcommand*, G4String);
+
+  private:
+    PrimaryGeneratorAction* 		fAction;
+    G4UIdirectory*                  fGunDir;
+    G4UIcmdWithAnInteger*			fSourceType;
+    G4UIcmdWithAnInteger*			fSourceDirectionType;
+    G4UIcmdWithAnInteger*			fSourceGeometry;
+    G4UIcmdWithADoubleAndUnit* fSourceEnergy;
+    G4UIcmdWithADouble* fSourcePolarizationAngle;
+    G4UIcmdWithADouble* fSourcePolarizationDegree;
+    G4UIcmdWithADoubleAndUnit* fSourcePositionX;
+    G4UIcmdWithADoubleAndUnit* fSourcePositionY;
+    G4UIcmdWithADoubleAndUnit* fSourcePositionZ;
+    G4UIcmdWithADoubleAndUnit* fSourceDiameter;
+};
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+#endif

include/RunAction.hh

+//
+// ********************************************************************
+// * License and Disclaimer                                           *
+// *                                                                  *
+// * The  Geant4 software  is  copyright of the Copyright Holders  of *
+// * the Geant4 Collaboration.  It is provided  under  the terms  and *
+// * conditions of the Geant4 Software License,  included in the file *
+// * LICENSE and available at  http://cern.ch/geant4/license .  These *
+// * include a list of copyright holders.                             *
+// *                                                                  *
+// * Neither the authors of this software system, nor their employing *
+// * institutes,nor the agencies providing financial support for this *
+// * work  make  any representation or  warranty, express or implied, *
+// * regarding  this  software system or assume any liability for its *
+// * use.  Please see the license in the file  LICENSE  and URL above *
+// * for the full disclaimer and the limitation of liability.         *
+// *                                                                  *
+// * This  code  implementation is the result of  the  scientific and *
+// * technical work of the GEANT4 collaboration.                      *
+// * By using,  copying,  modifying or  distributing the software (or *
+// * any work based  on the software)  you  agree  to acknowledge its *
+// * use  in  resulting  scientific  publications,  and indicate your *
+// * acceptance of all terms of the Geant4 Software license.          *
+// ********************************************************************
+//
+/// \file OpNovice/include/OpNoviceRunAction.hh
+/// \brief Definition of the OpNoviceRunAction class
+//
+//
+//
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+#ifndef RunAction_h
+#define RunAction_h 1
+
+#include "globals.hh"
+#include "G4UserRunAction.hh"
+#include "RunActionMessenger.hh"
+#include "G4AnalysisManager.hh"
+
+class DetectorConstruction;
+class PrimaryGeneratorAction;
+class RunActionMessenger;
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+class G4Timer;
+class G4Run;
+
+class RunAction : public G4UserRunAction
+{
+  public:
+    RunAction(DetectorConstruction* det, PrimaryGeneratorAction* prim=0);
+    ~RunAction();
+
+  public:
+    virtual void BeginOfRunAction(const G4Run* aRun);
+    virtual void EndOfRunAction(const G4Run* aRun);
+    void SetFileName(G4String);
+
+  private:
+    G4AnalysisManager* fMan;
+    DetectorConstruction*   fDetector;
+    PrimaryGeneratorAction* fPrimary;
+    G4Timer* fTimer;
+    G4String fFileName;
+
+    RunActionMessenger* fRunActionMessenger;
+};
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+#endif /*RunAction_h*/

include/RunActionMessenger.hh

+// #ifndef RunActionMessenger_h
+// #define RunActionMessenger_h 1
+//
+// #include "G4UImessenger.hh"
+// #include "globals.hh"
+//
+// class RunAction;
+// class G4UIdirectory;
+// class G4UIcmdWithAString;
+// class G4UIcmdWithADoubleAndUnit;
+// class G4UIcmdWithADouble;
+// class G4UIcmdWithoutParameter;
+// class G4UIcmdWithAnInteger;
+//
+// class RunActionMessenger: public G4UImessenger
+// {
+//   public:
+//
+//     RunActionMessenger(G4UserRunAction );
+//    ~RunActionMessenger();
+//
+//     virtual void SetNewValue(G4UIcommand*, G4String);
+//
+//   private:
+//
+//     G4Run*      fRunAction;
+//
+//     G4UIdirectory*             fPENDir;
+//     G4UIdirectory*             fRunDir;
+//     G4UIcmdWithAString*        fRunNameCMD;
+// };
+//
+// //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+//
+// #endif

include/SiliconPlateConstruction.hh

+#ifndef SILICONPLATECONSTUCTION_H
+#define SILICONPLATECONSTUCTION_H
+
+#include "globals.hh"
+#include "G4VUserDetectorConstruction.hh"
+#include "DetectorMessenger.hh"
+#include "G4SystemOfUnits.hh"
+
+
+class SiliconPlateConstruction
+{
+  public:
+    G4VSolid* ConstructPlate();
+
+  private:
+    G4double fSiliconPlate_h;
+    G4double fHolderWidth;
+};
+#endif

include/SoleilMaterials.hh

+/**
+ * @file PenMaterials.hh
+ * @author: (modified by) Luis Manzanillas
+ * @date 2020, Max Planck Institute for Physics
+ */
+
+
+#ifndef SoleilMaterials_H
+#define SoleilMaterials_H
+
+#include "globals.hh"
+
+
+class SoleilMaterials {
+public:
+  SoleilMaterials();
+  ~SoleilMaterials();
+  
+  void Construct();
+  
+private:
+  G4double lightYieldAntracene;
+    
+};
+
+#endif

include/SteppingAction.hh

+//
+// ********************************************************************
+// * License and Disclaimer                                           *
+// *                                                                  *
+// * The  Geant4 software  is  copyright of the Copyright Holders  of *
+// * the Geant4 Collaboration.  It is provided  under  the terms  and *
+// * conditions of the Geant4 Software License,  included in the file *
+// * LICENSE and available at  http://cern.ch/geant4/license .  These *
+// * include a list of copyright holders.                             *
+// *                                                                  *
+// * Neither the authors of this software system, nor their employing *
+// * institutes,nor the agencies providing financial support for this *
+// * work  make  any representation or  warranty, express or implied, *
+// * regarding  this  software system or assume any liability for its *
+// * use.  Please see the license in the file  LICENSE  and URL above *
+// * for the full disclaimer and the limitation of liability.         *
+// *                                                                  *
+// * This  code  implementation is the result of  the  scientific and *
+// * technical work of the GEANT4 collaboration.                      *
+// * By using,  copying,  modifying or  distributing the software (or *
+// * any work based  on the software)  you  agree  to acknowledge its *
+// * use  in  resulting  scientific  publications,  and indicate your *
+// * acceptance of all terms of the Geant4 Software license.          *
+// ********************************************************************
+//
+// $Id: SteppingAction.hh 68752 2013-04-05 10:23:47Z gcosmo $
+//
+/// \file optical/LXe/include/LXeSteppingAction.hh
+/// \brief Definition of the LXeSteppingAction class
+//
+#ifndef SteppingAction_H
+#define SteppingACtion_H 1
+
+#include "globals.hh"
+#include "G4UserSteppingAction.hh"
+
+#include "G4OpBoundaryProcess.hh"
+
+
+class EventAction;
+
+class SteppingAction : public G4UserSteppingAction
+{
+  public:
+
+    SteppingAction(EventAction* eventAction);
+    virtual ~SteppingAction();
+    virtual void UserSteppingAction(const G4Step*);
+
+    void SetOneStepPrimaries(G4bool b){fOneStepPrimaries=b;}
+    G4bool GetOneStepPrimaries(){return fOneStepPrimaries;}
+
+  private:
+
+    G4bool fOneStepPrimaries;
+
+    G4OpBoundaryProcessStatus fExpectedNextStatus;
+
+    EventAction* fEventAction;
+};
+
+#endif

input_files/.~lock.sipmQE.csv#

+,hayward,pclg-16.mppmu.mpg.de,13.08.2018 12:44,file:///home/iwsatlas1/hayward/.config/libreoffice/4;
\ No newline at end of file

input_files/BC408.csv

+499.9975281,2637.96921438,0.0735794201,1.638
+499.0131531,2637.96921438,0.0735794201,1.638
+497.9906311,2608.99345188,0.0971923978,1.638
+497.0058289,2622.29994345,0.0971923978,1.638
+495.9829102,2613.79187337,0.0971923978,1.638
+494.9976807,2603.73826131,0.0971923978,1.638
+494.0122681,2603.73826131,0.0971923978,1.638
+492.9887085,2590.49569101,0.0971923978,1.638
+492.0028687,2598.38399004,0.1271069389,1.638
+491.0168457,2598.38399004,0.1271069389,1.638
+489.9926453,2577.05542506,0.1271069389,1.638
+489.0061951,2564.05714059,0.1271069389,1.638
+487.9815979,2548.347806775,0.1271069389,1.638
+486.994751,2547.17915955,0.1271069389,1.638
+486.0077209,2539.136196315,0.1271069389,1.638
+484.9825134,2558.170144095,0.1605650934,1.638
+483.9950867,2542.053145245,0.1605650934,1.638
+483.0074158,2535.851649495,0.1605650934,1.638
+481.9815369,2546.09103633,0.1605650934,1.638
+480.9935608,2542.91445021,0.1605650934,1.638
+480.0053406,2534.316192555,0.1962429487,1.638
+479.0167847,2534.316192555,0.1962429487,1.638
+477.9901733,2537.71558458,0.1962429487,1.638
+477.0013428,2527.87653999,0.1962429487,1.638
+476.0123291,2541.11598456,0.1962429487,1.638
+474.9850769,2533.591619055,0.2321267072,1.638
+473.995636,2527.809791055,0.2321267072,1.638
+473.0060425,2508.79228482,0.2321267072,1.638
+472.0162354,2539.898498895,0.2321267072,1.638
+470.9881287,2519.72970048,0.2653402458,1.638
+469.9979553,2508.25898682,0.2653402458,1.638
+469.0075684,2510.4745092,0.2653402458,1.638
+468.0169067,2497.209343245,0.2653402458,1.638
+466.9880676,2483.85275016,0.3009036907,1.638
+465.9971313,2485.87373862,0.3009036907,1.638
+465.0059204,2485.87373862,0.342032518,1.638
+464.0145874,2482.69067997,0.342032518,1.638
+462.9849243,2457.825832335,0.342032518,1.638
+461.9931641,2442.67083426,0.3835639848,1.638
+461.0012512,2451.069351945,0.3835639848,1.638
+460.0091858,2416.89066903,0.3835639848,1.638
+459.0168152,2408.9391483,0.4144530728,1.638
+457.9862061,2406.85883574,0.4144530728,1.638
+456.9935608,2384.49520812,0.4501835674,1.638
+456.0006714,2382.02038893,0.4501835674,1.638
+455.0076294,2375.125178595,0.4501835674,1.58
+454.0144043,2359.919995455,0.4833418929,1.58
+452.9827271,2335.293207525,0.4833418929,1.58
+451.9891052,2335.293207525,0.5181897534,1.58
+450.9953613,2337.74506704,0.5181897534,1.58
+450.0013733,2329.04335257,0.5181897534,1.58
+449.0071716,2305.46831358,0.5576980796,1.58
+448.0128784,2306.587293105,0.5576980796,1.58
+447.0183411,2309.6601039,0.5576980796,1.58
+445.9853516,2281.91255616,0.5935497891,1.58
+444.990448,2281.91255616,0.5935497891,1.58
+443.9953918,2269.258074225,0.5935497891,1.58
+443.0001221,2252.6670906,0.6336132432,1.58
+442.0047302,2237.651696235,0.6722671765,1.58
+441.0090332,2219.98972629,0.7102842552,1.58
+440.0133057,2219.099689905,0.7102842552,1.58
+439.017334,2195.708355705,0.7493359978,1.58
+437.9828491,2156.42220642,0.8212886163,1.58
+436.9865112,2138.893565205,0.8929768043,1.58
+435.9900208,2118.61300251,0.9188273262,1.58
+434.9933167,2118.61300251,0.9555453103,1.58
+433.9964294,2089.30505946,0.9555453103,1.58
+432.9994202,2065.22979345,1,1.58
+432.0021362,2034.325722525,1,1.58
+431.0047607,2003.324659395,1,1.58
+430.0071411,1965.261281505,1,1.58
+429.0094604,1930.182494175,1,1.58
+428.0115051,1911.199187295,1,1.58
+427.0133972,1873.339853985,0.9748127801,1.58
+426.0151062,1845.794583555,0.9748127801,1.58
+425.0166626,1808.42680845,0.9383478032,1.58
+424.0180664,1760.687495655,0.9383478032,1.58
+422.980896,1684.109136045,0.9045387669,1.58
+421.9819031,1649.05928523,0.9045387669,1.58
+420.9827881,1608.124432425,0.8586336183,1.58
+419.983429,1569.42019701,0.8586336183,1.58
+418.9839783,1513.166508165,0.8137608066,1.58
+417.984314,1471.78243713,0.8137608066,1.58
+416.9845276,1420.110859155,0.7421656982,1.58
+415.9844971,1365.21822354,0.7421656982,1.58
+414.9843445,1310.214549885,0.7085935918,1.58
+413.9840393,1251.404446035,0.6604066717,1.58
+412.9834595,1188.74112192,0.6109184358,1.58
+411.9828491,1129.27714089,0.5709067043,1.58
+410.9819946,1059.51222561,0.5709067043,1.58
+409.9810181,990.855348855,0.5304646915,1.58
+409.0183411,990.855348855,0.4884942698,1.58
+408.0169983,921.709204005,0.4339902618,1.58
+407.0155029,846.175757805,0.3928507868,1.58
+406.0138245,768.784033815,0.3928507868,1.58
+405.0119934,696.447057525,0.3580601843,1.58
+404.0100098,622.65129687,0.3235818214,1.58
+403.0078735,548.89030866,0.2783967701,1.58
+402.0055542,480.56503656,0.2394436776,1.58
+401.0030518,358.826066295,0.2394436776,1.58
+400.0004272,307.247589945,0.2394436776,1.58
+399,0,0.2,1.58

input_files/OldPen.csv

+499.9975281,4.636854764,0.2311307823,1.638
+499.0131531,4.6223330011,0.2457378723,1.638
+497.9906311,4.6100458996,0.2671333553,1.638
+497.0058289,4.5999604972,0.2884423357,1.638
+495.9829102,4.5975784689,0.2993551854,1.638
+494.9976807,4.5902685335,0.3033910257,1.638
+494.0122681,4.5852167759,0.2892926152,1.638
+492.9887085,4.582127571,0.283545021,1.638
+492.0028687,4.5557538935,0.2819296273,1.638
+491.0168457,4.5543465332,0.2913973663,1.638
+489.9926453,4.5584481127,0.3023784924,1.638
+489.0061951,4.5541708295,0.3105717202,1.638
+487.9815979,4.5291191218,0.3357357563,1.638
+486.994751,4.524959898,0.3525517441,1.638
+486.0077209,4.5271542737,0.3689279544,1.638
+484.9825134,4.5186733443,0.3705296386,1.638
+483.9950867,4.492494371,0.3703527709,1.638
+483.0074158,4.4854299898,0.3748594628,1.638
+481.9815369,4.4933655324,0.3796494486,1.638
+480.9935608,4.4864501946,0.3924045953,1.638
+480.0053406,4.4804352639,0.3991128139,1.638
+479.0167847,4.4682122057,0.4121093745,1.638
+477.9901733,4.4619346818,0.4337728496,1.638
+477.0013428,4.453645325,0.4433133535,1.638
+476.0123291,4.4465876174,0.4510739933,1.638
+474.9850769,4.4359539988,0.4663618901,1.638
+473.995636,4.4237721498,0.4792405375,1.638
+473.0060425,4.4193593877,0.4851572702,1.638
+472.0162354,4.4045167479,0.4982230044,1.638
+470.9881287,4.3918101186,0.5244333082,1.638
+469.9979553,4.3852738355,0.5398743273,1.638
+469.0075684,4.3736873361,0.5485711491,1.638
+468.0169067,4.3623110664,0.5685419113,1.638
+466.9880676,4.3523993185,0.5764699301,1.638
+465.9971313,4.3356369873,0.6013593871,1.638
+465.0059204,4.3277635726,0.6069619849,1.638
+464.0145874,4.3133631277,0.6235914663,1.638
+462.9849243,4.3037170976,0.6408068741,1.638
+461.9931641,4.292715356,0.650616896,1.638
+461.0012512,4.2763422145,0.6762077441,1.638
+460.0091858,4.2563041167,0.6809739877,1.638
+459.0168152,4.2489215944,0.6980632626,1.638
+457.9862061,4.2461454552,0.7150664118,1.638
+456.9935608,4.2310963511,0.7245845296,1.638
+456.0006714,4.2178040764,0.7448038375,1.638
+455.0076294,4.2127418728,0.7496966674,1.638
+454.0144043,4.1927083531,0.7632892145,1.638
+452.9827271,4.1822418954,0.7815125954,1.638
+451.9891052,4.1709369102,0.7926547732,1.638
+450.9953613,4.1620948871,0.811922833,1.638
+450.0013733,4.1505857694,0.8165901094,1.638
+449.0071716,4.1397999242,0.8307459985,1.638
+448.0128784,4.1253951712,0.8472740767,1.638
+447.0183411,4.1062430914,0.857508622,1.638
+445.9853516,4.0936594226,0.8748593092,1.638
+444.990448,4.0829953022,0.8862892872,1.638
+443.9953918,4.0719910974,0.8925319275,1.638
+443.0001221,4.0592772985,0.9087250521,1.638
+442.0047302,4.0442170045,0.9148842647,1.638
+441.0090332,4.0294931026,0.9299703337,1.638
+440.0133057,4.0117320118,0.9348231219,1.638
+439.017334,4.0029342496,0.9442864625,1.638
+437.9828491,3.988445202,0.9563531829,1.638
+436.9865112,3.9714091296,0.9597094543,1.638
+435.9900208,3.9566052514,0.9704347349,1.638
+434.9933167,3.9364411035,0.9771234369,1.638
+433.9964294,3.9203283288,0.9804235187,1.638
+432.9994202,3.9068945006,0.9876752109,1.638
+432.0021362,3.8870645155,0.9901002146,1.638
+431.0047607,3.8677316852,0.9946755974,1.638
+430.0071411,3.8469201861,0.9939415363,1.638
+429.0094604,3.8282046513,0.9970100722,1.638
+428.0115051,3.8063759167,0.9961355224,1.638
+427.0133972,3.7809849749,0.9945080935,1.638
+426.0151062,3.7528954341,0.9934814392,1.638
+425.0166626,3.7213677103,0.9931296809,1.638
+424.0180664,3.695200467,0.9903432261,1.638
+422.980896,3.6630365015,0.9873671972,1.638
+421.9819031,3.6308458335,0.9853279914,1.638
+420.9827881,3.5916351855,0.9772126257,1.638
+419.983429,3.5526919937,0.9723126705,1.638
+418.9839783,3.5116679452,0.9684620259,1.638
+417.984314,3.4708162428,0.958219664,1.638
+416.9845276,3.4274689573,0.9475563932,1.638
+415.9844971,3.3783436091,0.9288847642,1.638
+414.9843445,3.3225901698,0.9175698245,1.638
+413.9840393,3.2672260846,0.9076564669,1.638
+412.9834595,3.2028559014,0.8867623724,1.638
+411.9828491,3.1284923121,0.8721513615,1.638
+410.9819946,3.043444648,0.84578891,1.638
+409.9810181,2.9512627115,0.8251778327,1.638
+409.0183411,2.8646068668,0.8016966096,1.638
+408.0169983,2.7562238017,0.7923277282,1.638
+407.0155029,2.63086472,0.7589677593,1.638
+406.0138245,2.4966808386,0.7466640726,1.638
+405.0119934,2.3512113107,0.7060971712,1.638
+404.0100098,2.2035180178,0.6973018006,1.638
+403.0078735,2.0385444504,0.6509535283,1.638
+402.0055542,1.8620816153,0.6358511177,1.638
+401.0030518,1.6825841994,0.5992219981,1.638
+400.0004272,1.5017908466,0.5856299569,1.638
+399,0,0.57,1.638

input_files/PEN_ABS_04Degree.csv

+385.0	0.5	1.62
+386.0	0.5	1.62
+387.0	0.5	1.62
+388.0	0.5	1.62
+389.0	0.5	1.62
+390.0	0.5	1.62
+391.0	0.5	1.62
+392.0	0.5	1.62
+393.0	0.5	1.62
+394.0	0.5	1.62
+395.0	0.5	1.62
+396.0	0.4573644558630363	1.62
+397.0	1.0740610589844586	1.62
+398.0	1.7072657476709239	1.62
+399.0	2.3574204197498254	1.62
+400.0	3.0249788020208506	1.62
+401.0	5.810850499940214	1.62
+402.0	5.579036955957408	1.62
+403.0	5.720068393193435	1.62
+404.0	6.000858949903191	1.62
+405.0	6.596741701274864	1.62
+406.0	7.175266253980685	1.62
+407.0	7.95933024864374	1.62
+408.0	8.904735497859637	1.62
+409.0	9.938842468654386	1.62
+410.0	10.930710719290358	1.62
+411.0	12.076420634142739	1.62
+412.0	13.088500191983051	1.62
+413.0	14.073937593186496	1.62
+414.0	15.19892212221158	1.62
+415.0	16.410148390465178	1.62
+416.0	17.421524507796825	1.62
+417.0	18.558903488687392	1.62
+418.0	19.349832321574727	1.62
+419.0	20.33319679097421	1.62
+420.0	21.258241611700168	1.62
+421.0	22.325860714844275	1.62
+422.0	23.336578695023963	1.62
+423.0	24.316836551679213	1.62
+424.0	24.99858043377279	1.62
+425.0	26.01006713331937	1.62
+426.0	26.755446795536773	1.62
+427.0	27.78836619648352	1.62
+428.0	28.853946577575783	1.62
+429.0	29.56966111545351	1.62
+430.0	30.173023992543776	1.62
+431.0	30.721586946295748	1.62
+432.0	32.135896494411384	1.62
+433.0	32.493790532871536	1.62
+434.0	33.034732341917355	1.62
+435.0	33.86532230933362	1.62
+436.0	34.459610363833846	1.62
+437.0	34.69280348527001	1.62
+438.0	35.85388269112591	1.62
+439.0	35.95167177377591	1.62
+440.0	37.06222287702453	1.62
+441.0	37.22794456662142	1.62
+442.0	37.59877587353113	1.62
+443.0	38.0759977010301	1.62
+444.0	38.33171072415995	1.62
+445.0	38.71898466192937	1.62
+446.0	39.343661250776364	1.62
+447.0	39.93226464799785	1.62
+448.0	40.42670361856144	1.62
+449.0	40.4978132092906	1.62
+450.0	41.91846813427499	1.62
+451.0	41.83871368417419	1.62
+452.0	41.469010992836935	1.62
+453.0	42.6926885641174	1.62
+454.0	42.909709713069496	1.62
+455.0	42.99953683212048	1.62
+456.0	43.25962168173261	1.62
+457.0	43.26706371146058	1.62
+458.0	43.54122928989146	1.62
+459.0	43.921057466114604	1.62
+460.0	45.129720926771746	1.62
+461.0	44.784302167470706	1.62
+462.0	45.192724143272855	1.62
+463.0	45.14966483941289	1.62
+464.0	45.3341155984293	1.62
+465.0	45.533174544974074	1.62
+466.0	45.93394087712653	1.62
+467.0	44.610517384648574	1.62
+468.0	45.639446365161305	1.62
+469.0	46.30164173891022	1.62
+470.0	46.37084746394715	1.62
+471.0	47.008241198314664	1.62
+472.0	46.29983080407387	1.62
+473.0	46.97826137191687	1.62
+474.0	46.858942162968795	1.62
+475.0	46.8924161417752	1.62
+476.0	47.187913881227004	1.62
+477.0	47.161979191279116	1.62
+478.0	47.68610916936892	1.62
+479.0	47.351715585994064	1.62
+480.0	46.66151664248891	1.62
+481.0	47.768700672556875	1.62
+482.0	47.542047861107676	1.62
+483.0	47.57897444171281	1.62
+484.0	48.35256361807604	1.62
+485.0	47.834603829586214	1.62
+486.0	47.827742516571064	1.62
+487.0	48.10589665740626	1.62
+488.0	49.43302251276107	1.62
+489.0	47.91207453211282	1.62
+490.0	48.43073743230514	1.62
+491.0	47.87306559027086	1.62
+492.0	47.77260066680517	1.62
+493.0	47.719152960614906	1.62
+494.0	48.89081550039993	1.62
+495.0	47.75201450134951	1.62
+496.0	48.147452648384736	1.62
+497.0	49.11124452771704	1.62
+498.0	49.67915536871652	1.62
+499.0	49.04706451474215	1.62
+500.0	49.657850595239225	1.62
+501.0	48.428819017725715	1.62
+502.0	47.263254645398966	1.62
+503.0	48.87699167493794	1.62
+504.0	47.83048198679239	1.62
+505.0	47.71923864048888	1.62
+506.0	48.09293484119233	1.62
+507.0	49.28514442176318	1.62
+508.0	49.390024122313456	1.62
+509.0	49.608961186582235	1.62
+510.0	48.445674541083704	1.62
+511.0	48.617349084832895	1.62
+512.0	51.75626489653614	1.62
+513.0	51.08181268328694	1.62
+514.0	49.67331250799426	1.62
+515.0	49.53781833671384	1.62
+516.0	50.75417223389212	1.62
+517.0	50.27398553350497	1.62
+518.0	50.01431603661094	1.62
+519.0	50.754127814096684	1.62
+520.0	50.09099141646328	1.62
+521.0	51.897643563718816	1.62
+522.0	49.996162405382584	1.62
+523.0	51.30569590452687	1.62
+524.0	49.48312141897378	1.62
+525.0	50.627586227688326	1.62
+526.0	52.63105799375451	1.62
+527.0	53.931801006528644	1.62
+528.0	54.22016804490344	1.62
+529.0	52.63888751542547	1.62
+530.0	49.55946365805137	1.62
+531.0	54.004665217560955	1.62
+532.0	53.49736156388983	1.62
+533.0	54.909876898410246	1.62
+534.0	54.75938305376352	1.62
+535.0	54.61479799486388	1.62
+536.0	54.205351434181054	1.62
+537.0	55.239817822116514	1.62
+538.0	53.83712260936581	1.62
+539.0	55.57197133501288	1.62
+540.0	56.8525222673966	1.62
+541.0	58.725745467187686	1.62
+542.0	54.10466954279019	1.62
+543.0	53.06013988226705	1.62
+544.0	59.08128537230183	1.62
+545.0	56.03324782308435	1.62
+546.0	55.040673707008146	1.62
+547.0	59.15951545376455	1.62
+548.0	59.156402114072854	1.62
+549.0	54.40516796760547	1.62
+550.0	54.663638139833665	1.62
+551.0	54.710727326578606	1.62
+552.0	55.45471567848257	1.62
+553.0	52.94362220837098	1.62
+554.0	54.90272054042494	1.62
+555.0	57.86114827826741	1.62
+556.0	57.7433316118438	1.62
+557.0	56.71265418472376	1.62
+558.0	62.0913359539485	1.62
+559.0	57.71662210560112	1.62
+560.0	60.63191484712538	1.62
+561.0	68.84901866853032	1.62
+562.0	53.93076508530472	1.62
+563.0	58.691841156659336	1.62
+564.0	60.294554482233174	1.62
+565.0	56.28094078564472	1.62
+566.0	65.39158103616748	1.62
+567.0	55.59833878805142	1.62
+568.0	58.45115112464616	1.62
+569.0	66.55507125716012	1.62
+570.0	62.133406233141976	1.62
+571.0	57.925907883163305	1.62
+572.0	57.76034865323523	1.62
+573.0	63.546758360279924	1.62
+574.0	52.90101646470013	1.62
+575.0	66.50	1.62