+++ /dev/null
-# $Id$
-# Flugg tag $Name$
-
-# Id: GNUmakefile,v 2.1 1998/07/02 15:26:52 gunter Exp
-# --------------------------------------------------------------
-# GNUmakefile for examples module. Gabriele Cosmo, 06/04/98
-# --------------------------------------------------------------
-
-name := mainAlAuAl
-G4TARGET := $(name)
-G4EXLIB := true
-
-ifndef FLUGGINSTALL
- FLUGGINSTALL = ../../..
-endif
-
-.PHONY: all
-all: lib bin
-
-include $(FLUGGINSTALL)/config/binmake.gmk
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#ifndef MyDetectorConstruction_h
-#define MyDetectorConstruction_h 1
-
-#include "G4VUserDetectorConstruction.hh"
-#include "globals.hh"
-
-class MyDetectorConstruction : public G4VUserDetectorConstruction
-{
- public:
- MyDetectorConstruction();
- ~MyDetectorConstruction();
- G4VPhysicalVolume* Construct();
-
- private:
- G4double expHall_x;
- G4double expHall_y;
- G4double expHall_z;
-
- G4double myBox_x;
- G4double myBox_y;
- G4double myBox_zA;
- G4double myBox_zB;
- G4double myBox_zC;
- G4double myBox_zD;
- G4double myBox_zE;
-};
-
-#endif
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#include "FGeometryInit.hh"
-#include "MyDetectorConstruction.hh"
-
-#define flukam flukam_
-
-extern "C" void flukam(const G4int & GeoFlag);
-
-int main() {
-
- FGeometryInit* theFGeometryInit = FGeometryInit::GetInstance();
-
- theFGeometryInit
- ->setDetConstruction(new MyDetectorConstruction());
-
-//flag for geometry:
-// 1 for GEANT4
-// 0 for FLUKA
-// 2 for Rubia
- const G4int flag = 1;
-
-//call fortran
- flukam(flag);
-
-//end
- return 0;
-}
-
-
-
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#include "MyDetectorConstruction.hh"
-
-#include "G4Material.hh"
-#include "G4MaterialTable.hh"
-#include "G4Element.hh"
-#include "G4ElementTable.hh"
-#include "G4Box.hh"
-#include "G4Sphere.hh"
-#include "G4LogicalVolume.hh"
-#include "G4LogicalVolumeStore.hh"
-#include "G4RotationMatrix.hh"
-#include "G4ThreeVector.hh"
-#include "G4Transform3D.hh"
-#include "G4PVPlacement.hh"
-
-MyDetectorConstruction::MyDetectorConstruction()
-{
- expHall_x = 10.*cm;
- expHall_y = 10.*cm;
- expHall_z = 10.*cm;
-
- myBox_x = 10.*cm;
- myBox_y = 10.*cm;
- myBox_zA = 5.*cm;
- myBox_zB = 0.00841665*cm;
- myBox_zC = 0.0010873*cm;
- myBox_zD = 0.07640375*cm;
- myBox_zE = 4.9140923*cm;
-
-}
-
-MyDetectorConstruction::~MyDetectorConstruction()
-{;}
-
-G4VPhysicalVolume* MyDetectorConstruction::Construct()
-{
-//===================================volumes
- G4cout << "MyDetectorConstruction::Construct start" << G4endl;
-
-
- //------------------------------ experimental hall
- G4Box * experimantalHall_box
- = new G4Box("expHall_b",expHall_x,expHall_y,expHall_z);
- G4LogicalVolume * experimantalHall_log
- = new G4LogicalVolume(experimantalHall_box,0,"expHall_L",0,0,0);
- G4VPhysicalVolume * experimantalHall_phys
- = new G4PVPlacement(0,G4ThreeVector(),"expHall_P",
- experimantalHall_log,0,false,0);
-
- //------------------------------ test box2
- G4Box * test_box2
- = new G4Box("myBox_b2",myBox_x,myBox_y,myBox_zA);
- G4LogicalVolume * test_log2
- = new G4LogicalVolume(test_box2,0,"myBox_L2",0,0,0);
- G4VPhysicalVolume * test_phys2
- = new G4PVPlacement(0,G4ThreeVector(0.,0.,-5.*cm),"myBox_P2",
- test_log2,experimantalHall_phys,false,0);
-
- //------------------------------ test box3
- G4Box * test_box3
- = new G4Box("myBox_b3",myBox_x,myBox_y,myBox_zB);
- G4LogicalVolume * test_log3
- = new G4LogicalVolume(test_box3,0,"myBox_L3",0,0,0);
- G4VPhysicalVolume * test_phys3
- = new G4PVPlacement(0,G4ThreeVector(0.,0.,0.00841665*cm),"myBox_P3",
- test_log3,experimantalHall_phys,false,0);
-
- //------------------------------ test box4
- G4Box * test_box4
- = new G4Box("myBox_b4",myBox_x,myBox_y,myBox_zC);
- G4LogicalVolume * test_log4
- = new G4LogicalVolume(test_box4,0,"myBox_L4",0,0,0);
- G4VPhysicalVolume * test_phys4
- = new G4PVPlacement(0,G4ThreeVector(0.,0.,0.0179206*cm),"myBox_P4",
- test_log4,experimantalHall_phys,false,0);
-
- //------------------------------ test box5
- G4Box * test_box5
- = new G4Box("myBox_b5",myBox_x,myBox_y,myBox_zD);
- G4LogicalVolume * test_log5
- = new G4LogicalVolume(test_box5,0,"myBox_L5",0,0,0);
- G4VPhysicalVolume * test_phys5
- = new G4PVPlacement(0,G4ThreeVector(0.,0.,0.09541165*cm),"myBox_P5",
- test_log5,experimantalHall_phys,false,0);
-
- //------------------------------ test box6
- G4Box * test_box6
- = new G4Box("myBox_b6",myBox_x,myBox_y,myBox_zE);
- G4LogicalVolume * test_log6
- = new G4LogicalVolume(test_box6,0,"myBox_L6",0,0,0);
- G4VPhysicalVolume * test_phys6
- = new G4PVPlacement(0,G4ThreeVector(0.,0.,5.0859077*cm),"myBox_P6",
- test_log6,experimantalHall_phys,false,0);
-
- //------------------------------------------------------------------
- G4int numLVVol = G4int(G4LogicalVolumeStore::GetInstance()->size());
-
- if (test_log6->GetMaterial())
- G4cout << "material is defined " << G4endl;
- else
- G4cout << "material is NOT defined " << G4endl;
-
- return experimantalHall_phys;
-}
+++ /dev/null
-# $Id$
-# Flugg tag $Name$
-
-# Id: GNUmakefile,v 2.1 1998/07/02 15:26:52 gunter Exp
-# --------------------------------------------------------------
-# GNUmakefile for examples module. Gabriele Cosmo, 06/04/98.
-# --------------------------------------------------------------
-
-name := mainBias
-G4TARGET := $(name)
-G4EXLIB := true
-
-ifndef FLUGGINSTALL
- FLUGGINSTALL = ../../..
-endif
-
-#the target of the first rule of the first makefile is the default goal
-#so "make" will do lib and bin targets.
-.PHONY: all
-all: lib bin
-
-include $(FLUGGINSTALL)/config/binmake.gmk
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#ifndef MyDetectorConstruction_h
-#define MyDetectorConstruction_h 1
-
-#include "G4VUserDetectorConstruction.hh"
-#include "globals.hh"
-
-class MyDetectorConstruction : public G4VUserDetectorConstruction
-{
- public:
- MyDetectorConstruction();
- ~MyDetectorConstruction();
- G4VPhysicalVolume* Construct();
-
- private:
- G4double expHall_rad;
- G4double expHall_z;
-
- G4double tar_rad;
- G4double tar_z;
-
- G4double litCil_rad;
- G4double litCil_z;
-
- G4double bigCil_rad;
- G4double bigCil_z;
-};
-
-#endif
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#include "FGeometryInit.hh"
-#include "MyDetectorConstruction.hh"
-
-#define flukam flukam_
-
-extern "C" void flukam(const G4int & GeoFlag);
-
-int main() {
-
- FGeometryInit* theFGeometryInit = FGeometryInit::GetInstance();
-
- theFGeometryInit
- ->setDetConstruction(new MyDetectorConstruction());
-
-//flag for geometry:
-// 1 for GEANT4
-// 0 for FLUKA
-// 2 for Rubia
- const G4int flag = 1;
-
-//call fortran
- flukam(flag);
-
-//end
- return 0;
-}
-
-
-
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#include "MyDetectorConstruction.hh"
-
-#include "G4Box.hh"
-#include "G4Tubs.hh"
-#include "G4Sphere.hh"
-#include "G4LogicalVolume.hh"
-#include "G4RotationMatrix.hh"
-#include "G4ThreeVector.hh"
-#include "G4Transform3D.hh"
-#include "G4PVPlacement.hh"
-
-MyDetectorConstruction::MyDetectorConstruction()
-{
- expHall_rad = 1000.*cm;
- expHall_z = 1000.*cm;
-
- tar_rad = 3.*cm;
- tar_z = 10.*cm;
-
- litCil_rad = 1000.*cm;
- litCil_z = 10.*cm;
-
- bigCil_rad = 1000.*cm;
- bigCil_z = 250.*cm;
-}
-
-MyDetectorConstruction::~MyDetectorConstruction()
-{;}
-
-G4VPhysicalVolume* MyDetectorConstruction::Construct()
-{
-//=======================volumes
-//-------------------- experimental hall
- G4Tubs * experimentalHall_tub
- = new G4Tubs("expHall_S",0.*cm,expHall_rad,expHall_z,0.*deg,360.*deg);
- G4LogicalVolume * experimentalHall_log
- = new G4LogicalVolume(experimentalHall_tub,0,"expHall_L",0,0,0);
- G4VPhysicalVolume * experimentalHall_phys
- = new G4PVPlacement(0,G4ThreeVector(),"expHall_P",
- experimentalHall_log,0,false,0);
-
- //------------------------------ big cylinder
- G4Tubs * bigCil_tub
- = new G4Tubs("bigCil_S",0.*cm,bigCil_rad,bigCil_z,0.*deg,360.*deg);
- G4LogicalVolume * bigCil_log
- = new G4LogicalVolume(bigCil_tub,0,"bigCil_L",0,0,0);
- G4VPhysicalVolume * bigCil_phys
- = new G4PVPlacement(0,G4ThreeVector(0.*cm,0.*cm,750.*cm),"bigCil_P",
- bigCil_log,experimentalHall_phys,false,0);
-
- //------------------------------ target
- G4Tubs * target_tub
- = new G4Tubs("tar_S",0.*cm,tar_rad,tar_z,0.*deg,360.*deg);
- G4LogicalVolume * target_log
- = new G4LogicalVolume(target_tub,0,"tar_L",0,0,0);
- G4VPhysicalVolume * target_phys
- = new G4PVPlacement(0,G4ThreeVector(0.*cm,0.*cm,-10.*cm),"tar_P",
- target_log,experimentalHall_phys,false,0);
-
- //------------------------------ little cylinders
- G4Tubs * litCil_tub
- = new G4Tubs("litCil_S",0.*cm,litCil_rad,litCil_z,0.*deg,360.*deg);
- G4LogicalVolume * litCil_log
- = new G4LogicalVolume(litCil_tub,0,"litCil_L",0,0,0);
-
- /* G4VPhysicalVolume * litCil_phy1
- = new G4PVPlacement(0,G4ThreeVector(0.*cm,0.*cm,10.*cm),
- "litCil_P1",litCil_log,experimentalHall_phys,false,0);
- G4VPhysicalVolume * litCil_phy2
- = new G4PVPlacement(0,G4ThreeVector(0.*cm,0.*cm,30.*cm),
- "litCil_P2",litCil_log,experimentalHall_phys,false,0);
- G4VPhysicalVolume * litCil_phy3
- = new G4PVPlacement(0,G4ThreeVector(0.*cm,0.*cm,50.*cm),
- "litCil_P3",litCil_log,experimentalHall_phys,false,0);
- G4VPhysicalVolume * litCil_phy4
- = new G4PVPlacement(0,G4ThreeVector(0.*cm,0.*cm,70.*cm),
- "litCil_P4",litCil_log,experimentalHall_phys,false,0);
- */
-
- for(G4int i=0;i<25;i++)
- {
- cout<<"Cilindretto num. "<<i<<endl;
- new G4PVPlacement(0,G4ThreeVector(0.*cm,0.*cm,(1+2*i)*10.*cm),"litCil_P",
- litCil_log,experimentalHall_phys,false,i);
- }
-
- //------------------------------------------------------------------
- return experimentalHall_phys;
-}
+++ /dev/null
-# $Id$
-# --------------------------------------------------------------
-# GNUmakefile for examples module. Gabriele Cosmo, 06/04/98.
-# --------------------------------------------------------------
-
-ifndef FLUGGINSTALL
- FLUGGINSTALL = ../..
-endif
-
-include $(FLUGGINSTALL)/config/architecture.gmk
-
-SUBDIRS = AlAuAl BiasEx MagSphereNotRep T36flugg
-
-.PHONY : all clean clean_libs
-
-all:
- @for dir in $(SUBDIRS); do (cd $$dir; $(MAKE)); done
-
-clean:
- @for dir in $(SUBDIRS); do (cd $$dir; $(MAKE) clean); done
-
-clean_libs:
- @for dir in $(SUBDIRS); do (cd $$dir; $(MAKE) clean_libs); done
+++ /dev/null
-# $Id$
-# Flugg tag $Name$
-
-# Id: GNUmakefile,v 2.1 1998/07/02 15:26:52 gunter Exp
-# --------------------------------------------------------------
-# GNUmakefile for examples module. Gabriele Cosmo, 06/04/98.
-# --------------------------------------------------------------
-
-name := mainSphereNotRep
-G4TARGET := $(name)
-G4EXLIB := true
-
-ifndef FLUGGINSTALL
- FLUGGINSTALL = ../../..
-endif
-
-#the target of the first rule of the first makefile is the default goal
-#so "make" will do lib and bin targets.
-.PHONY: all
-all: lib bin
-
-include $(FLUGGINSTALL)/config/binmake.gmk
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#ifndef ExN04Field_H
-#define ExN04Field_H 1
-
-#include "globals.hh"
-#include "G4MagneticField.hh"
-
-class ExN04Field : public G4MagneticField
-{
- public:
- ExN04Field();
- ~ExN04Field();
-
- void GetFieldValue( const double Point[3],
- double *Bfield ) const;
-
- private:
- G4double Bz;
- G4double rmax_sq;
- G4double zmax;
-};
-
-#endif
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-//
-// MyDetectorConstruction.hh, 2000/02/11 for flugg
-// Sara Vanini
-//
-//
-
-
-#ifndef MyDetectorConstruction_h
-#define MyDetectorConstruction_h 1
-
-#include "G4VUserDetectorConstruction.hh"
-#include "globals.hh"
-
-class G4Box;
-class G4Tubs;
-class G4Sphere;
-class G4VSolid;
-class G4LogicalVolume;
-class G4VPhysicalVolume;
-class G4Material;
-class G4UniformMagField;
-class G4IntersectionSolid;
-
-
-class MyDetectorConstruction : public G4VUserDetectorConstruction
-{
- public:
-
- MyDetectorConstruction();
- ~MyDetectorConstruction();
-
- public:
-
- G4VPhysicalVolume* Construct();
- // SetMagField(G4double fieldValue);
-
- private:
-
- G4Material* TubMaterial;
- G4double TubRad;
-
- G4Material* SphereMaterial;
- G4double SphereRad;
-
- G4int NbOfLayers,NbOfTubs,NbOfSpheres;
- G4double LayerThickness;
-
- G4double DetSizeX,DetSizeY,DetSizeZ;
-
- G4Material* defaultMaterial, *LayerMaterial;
- G4double WorldSize;
-
- G4Box* solidWorld; //pointer to the solid World
- G4LogicalVolume* logicWorld; //pointer to the logical World
- G4VPhysicalVolume* physiWorld; //pointer to the physical World
-
- G4Box* solidDet; //pointer to the solid Det
- G4LogicalVolume* logicDet; //pointer to the logical Det
- G4VPhysicalVolume* physiDet; //pointer to the physical Det
-
- G4Box* solidLayer; //pointer to the solid Layer
- G4LogicalVolume* logicLayer; //pointer to the logical Layer
- G4VPhysicalVolume* physiLayer1,* physiLayer2,* physiLayer3,
- * physiLayer4,* physiLayer5,* physiLayer6;
- //pointer to the physical Layer
-
- G4Box* solidTubLayer; //pointer to the solid Layer
- G4LogicalVolume* logicTubLayer; //pointer to the logical Layer
- G4VPhysicalVolume* physiTubLayer; //pointer to the physical Layer
-
- G4Sphere* solidSphere; //pointer to the solid Sphere
- G4LogicalVolume* logicSphere; //pointer to the logical Sphere
- G4VPhysicalVolume* physiSphere; //pointer to the physical Sphere
-
- G4Tubs* solidTub; //pointer to the solid Tub
- G4LogicalVolume* logicTub; //pointer to the logical Tub
- G4VPhysicalVolume* physiTub; //pointer to the physical Tub
-
- G4Box* solidBoxIntersTub;
- G4VSolid* solidTubSeg; //pointer to the solid Tub Segment
- G4LogicalVolume* logicTubSeg; //pointer to the logical Tub Segment
- G4VPhysicalVolume* physiTubSeg; //pointer to the physical Tub Segment
-
- G4Box* solidBoxIntersSph;
- G4VSolid* solidSphSeg1,*solidSphSeg2,*solidSphSeg3,*solidSphSeg4;
- //pointers to the solid Sphere Segment
- G4LogicalVolume* logicSphSeg1,*logicSphSeg2,*logicSphSeg3,*logicSphSeg4;
- //pointer to the logical Sphere Segment
- G4VPhysicalVolume* physiSphSeg1,*physiSphSeg2,*physiSphSeg3,*physiSphSeg4;
- //pointer to the physical Sphere Segment
-
-
-
-
- G4UniformMagField* magField; //pointer to the magnetic field
-
-
- private:
-
- void DefineMaterials();
- void ComputeDetectorParameters();
- G4VPhysicalVolume* ConstructDetector();
-};
-
-
-
-inline void MyDetectorConstruction::ComputeDetectorParameters()
-{
- // Compute derived parameters of the calorimeter
- LayerThickness = 2*SphereRad + 2*TubRad;
- DetSizeY = LayerThickness * NbOfLayers;
-}
-
-#endif
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#include "FGeometryInit.hh"
-#include "MyDetectorConstruction.hh"
-
-#define flukam flukam_
-
-extern "C" void flukam(const G4int & GeoFlag);
-
-int main() {
-
- FGeometryInit* theFGeometryInit = FGeometryInit::GetInstance();
-
- theFGeometryInit
- ->setDetConstruction(new MyDetectorConstruction());
-
-//flag for geometry:
-// 1 for GEANT4
-// 0 for FLUKA
-// 2 for Rubbia
- const G4int flag = 1;
-
-//call fortran
- flukam(flag);
-
-//end
- return 0;
-}
-
-
-
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#include "ExN04Field.hh"
-
-ExN04Field::ExN04Field()
-{
- Bz = 3.0*tesla;
- rmax_sq = sqr(50.*cm);
- zmax = 100.*cm;
-}
-
-ExN04Field::~ExN04Field()
-{;}
-
-void ExN04Field::GetFieldValue(const double Point[3],double *Bfield) const
-{
- Bfield[0] = 0.;
- Bfield[1] = 0.;
- if(abs(Point[2])<zmax && (sqr(Point[0])+sqr(Point[1]))<rmax_sq)
- { Bfield[2] = Bz; }
- else
- { Bfield[2] = 0.; }
-}
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-//
-// Example with Sphere and Tub layers for testing Magnetic Field in FLUGG
-// Sara Vanini, 11/02/00. Not replicated volumes!
-//
-
-
-#include "MyDetectorConstruction.hh"
-#include "G4Material.hh"
-#include "G4Tubs.hh"
-#include "G4Box.hh"
-#include "G4Sphere.hh"
-#include "G4LogicalVolume.hh"
-#include "G4PVPlacement.hh"
-#include "G4PVReplica.hh"
-#include "G4UniformMagField.hh"
-#include "G4FieldManager.hh"
-#include "G4VisAttributes.hh"
-#include "G4Colour.hh"
-#include "G4ios.hh"
-#include "G4IntersectionSolid.hh"
-
-
-MyDetectorConstruction::MyDetectorConstruction()
-:solidWorld(NULL),logicWorld(NULL),physiWorld(NULL),
- solidDet(NULL),logicDet(NULL),physiDet(NULL),
- solidLayer(NULL),logicLayer(NULL),physiLayer1(NULL),
- physiLayer2(NULL),physiLayer3(NULL),physiLayer4(NULL),
- physiLayer5(NULL),physiLayer6(NULL),
- solidTub(NULL),logicTub(NULL),physiTub(NULL),
- solidSphere(NULL),logicSphere(NULL),physiSphere(NULL),
- TubMaterial(NULL),SphereMaterial(NULL),defaultMaterial(NULL),
- magField(NULL)
-{
- // default parameter values of the calorimeter
- WorldSize = 200.*cm;
- TubRad = 3.*mm;
- SphereRad = 5.*mm;
- NbOfLayers = 6;
- DetSizeX = 5.1*cm;
- DetSizeZ = 8.*cm;
- NbOfTubs = 13;
- NbOfSpheres = 40;
-}
-
-
-
-MyDetectorConstruction::~MyDetectorConstruction()
-{}
-
-
-
-G4VPhysicalVolume* MyDetectorConstruction::Construct()
-{
- DefineMaterials();
- return ConstructDetector();
-}
-
-
-
-void MyDetectorConstruction::DefineMaterials()
-{
- //This function illustrates the possible ways to define materials
-
-G4String name, symbol; //a=mass of a mole;
-G4double a, z, density; //z=mean number of protons;
-G4int iz, n; //iz=number of protons in an isotope;
- // n=number of nucleons in an isotope;
-
-G4int ncomponents, natoms;
-G4double abundance, fractionmass;
-G4double temperature, pressure;
-
-//
-// define Elements
-//
-
-a = 1.01*g/mole;
-G4Element* H = new G4Element(name="Hydrogen",symbol="H" , z= 1., a);
-
-a = 12.01*g/mole;
-G4Element* C = new G4Element(name="Carbon" ,symbol="C" , z= 6., a);
-
-a = 14.01*g/mole;
-G4Element* N = new G4Element(name="Nitrogen",symbol="N" , z= 7., a);
-
-a = 16.00*g/mole;
-G4Element* O = new G4Element(name="Oxygen" ,symbol="O" , z= 8., a);
-
-a = 28.09*g/mole;
-G4Element* Si = new G4Element(name="Silicon",symbol="Si" , z= 14., a);
-
-a = 55.85*g/mole;
-G4Element* Fe = new G4Element(name="Iron" ,symbol="Fe", z=26., a);
-
-//
-// define an Element from isotopes, by relative abundance
-//
-
-G4Isotope* U5 = new G4Isotope(name="U235", iz=92, n=235, a=235.01*g/mole);
-G4Isotope* U8 = new G4Isotope(name="U238", iz=92, n=238, a=238.03*g/mole);
-
-G4Element* U = new G4Element(name="enriched Uranium", symbol="U", ncomponents=2);
-U->AddIsotope(U5, abundance= 90.*perCent);
-U->AddIsotope(U8, abundance= 10.*perCent);
-
-//
-// define simple materials
-//
-
-density = 2.700*g/cm3;
-a = 26.98*g/mole;
-G4Material* Al = new G4Material(name="Aluminium", z=13., a, density);
-
-density = 1.390*g/cm3;
-a = 39.95*g/mole;
-G4Material* lAr = new G4Material(name="liquidArgon", z=18., a, density);
-
-density = 8.960*g/cm3;
-a = 63.55*g/mole;
-G4Material* Cu = new G4Material(name="Copper" , z=29., a, density);
-
-density = 11.35*g/cm3;
-a = 207.19*g/mole;
-G4Material* Pb = new G4Material(name="Lead" , z=82., a, density);
-
-//
-// define a material from elements. case 1: chemical molecule
-//
-
-density = 1.000*g/cm3;
-G4Material* H2O = new G4Material(name="Water", density, ncomponents=2);
-H2O->AddElement(H, natoms=2);
-H2O->AddElement(O, natoms=1);
-
-density = 1.032*g/cm3;
-G4Material* Sci = new G4Material(name="Scintillator", density, ncomponents=2);
-Sci->AddElement(C, natoms=9);
-Sci->AddElement(H, natoms=10);
-
-density = 2.200*g/cm3;
-G4Material* SiO2 = new G4Material(name="quartz", density, ncomponents=2);
-SiO2->AddElement(Si, natoms=1);
-SiO2->AddElement(O , natoms=2);
-
-//
-// define a material from elements. case 2: mixture by fractional mass
-//
-
-density = 1.290*mg/cm3;
-G4Material* Air = new G4Material(name="Air" , density, ncomponents=2);
-Air->AddElement(N, fractionmass=0.7);
-Air->AddElement(O, fractionmass=0.3);
-
-//
-// define a material from elements and/or others materials (mixture of mixtures)
-//
-
-density = 0.200*g/cm3;
-G4Material* Aerog = new G4Material(name="Aerogel", density, ncomponents=3);
-Aerog->AddMaterial(SiO2, fractionmass=62.5*perCent);
-Aerog->AddMaterial(H2O , fractionmass=37.4*perCent);
-Aerog->AddElement (C , fractionmass= 0.1*perCent);
-
-//
-// examples of gas in non STP conditions
-//
-
-density = 27.*mg/cm3;
-pressure = 50.*atmosphere;
-temperature = 325.*kelvin;
-G4Material* CO2 = new G4Material(name="CarbonicGas", density, ncomponents=2,
- kStateGas,temperature,pressure);
-CO2->AddElement(C, natoms=1);
-CO2->AddElement(O, natoms=2);
-
-density = 0.3*mg/cm3;
-pressure = 2.*atmosphere;
-temperature = 500.*kelvin;
-G4Material* steam = new G4Material(name="WaterSteam", density, ncomponents=1,
- kStateGas,temperature,pressure);
-steam->AddMaterial(H2O, fractionmass=1.);
-
-//
-// examples of vacuum
-//
-
-density = universe_mean_density; //from PhysicalConstants.h
-pressure = 3.e-18*pascal;
-temperature = 2.73*kelvin;
-new G4Material(name="Galactic", z=1., a=1.01*g/mole, density,
- kStateGas,temperature,pressure);
-
-density = 1.e-5*g/cm3;
-pressure = 2.e-2*bar;
-temperature = STP_Temperature; //from PhysicalConstants.h
-G4Material* beam = new G4Material(name="Beam", density, ncomponents=1,
- kStateGas,temperature,pressure);
-beam->AddMaterial(Air, fractionmass=1.);
-
-//G4cout << *(G4Material::GetMaterialTable()) << endl;
-
- //default materials of the calorimeter
- TubMaterial = Al;
- SphereMaterial = Sci;
- LayerMaterial = Pb;
- defaultMaterial = Air;
-}
-
-
-
-G4VPhysicalVolume* MyDetectorConstruction::ConstructDetector()
-{
- // complete the Detector parameters definition
- ComputeDetectorParameters();
-
- //
- // World
- //
- solidWorld = new G4Box("World", //its name
- WorldSize/2,WorldSize/2,WorldSize/2); //its size
-
- logicWorld = new G4LogicalVolume(solidWorld, //its solid
- defaultMaterial, //its material
- "World"); //its name
-
- physiWorld = new G4PVPlacement(0, //no rotation
- G4ThreeVector(), //at (0,0,0)
- "World", //its name
- logicWorld, //its logical volume
- NULL, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
- //
- // Detector
- //
- solidDet=NULL; logicDet=NULL; physiDet=NULL;
- solidLayer=NULL; logicLayer=NULL; physiLayer1=NULL;
- physiLayer2=NULL; physiLayer3=NULL; physiLayer4=NULL;
- physiLayer5=NULL; physiLayer6=NULL;
- solidTubLayer=NULL; logicTubLayer=NULL; physiTubLayer=NULL;
-
- solidDet = new G4Box("Detector", //its name
- DetSizeX/2,DetSizeY/2,DetSizeZ/2);//size
-
- logicDet = new G4LogicalVolume(solidDet, //its solid
- defaultMaterial, //its material
- "Detector"); //its name
-
- physiDet = new G4PVPlacement(0, //no rotation
- G4ThreeVector(2.55*cm,-0.3*cm,4.0*cm),
- "Detector", //its name
- logicDet, //its logical volume
- physiWorld, //its mother volume
- false, //no boolean operation
- 0); //copy number
- //
- // Layer
- //
- solidLayer = new G4Box("Layer", //its name
- DetSizeX/2,LayerThickness/2,DetSizeZ/2); //size
-
- logicLayer = new G4LogicalVolume(solidLayer, //its solid
- LayerMaterial, //its material
- "Layer"); //its name
- /*
- if (NbOfLayers > 1)
- physiLayer = new G4PVReplica("Layer", //its name
- logicLayer, //its logical volume
- physiDet, //its mother
- kYAxis, //axis of replication
- NbOfLayers, //number of replica
- LayerThickness); //witdth of replica
- else
- physiLayer = new G4PVPlacement(0, //no rotation
- G4ThreeVector(), //at (0,0,0)
- "Layer", //its name
- logicLayer, //its logical volume
- physiDet, //its mother volume
- false, //no boolean operation
- 0); //copy number
- */
-
- physiLayer1 = new G4PVPlacement(0, //no rotation
- G4ThreeVector(0,
- -DetSizeY/2+LayerThickness/2,0),
- "Layer 1", //its name
- logicLayer, //its logical volume
- physiDet, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
- physiLayer2 = new G4PVPlacement(0, //no rotation
- G4ThreeVector(0,
- -DetSizeY/2+3*LayerThickness/2,0),
- "Layer 2", //its name
- logicLayer, //its logical volume
- physiDet, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
- physiLayer3 = new G4PVPlacement(0, //no rotation
- G4ThreeVector(0,
- -LayerThickness/2,0),
- "Layer 3", //its name
- logicLayer, //its logical volume
- physiDet, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
- physiLayer4 = new G4PVPlacement(0, //no rotation
- G4ThreeVector(0,
- +LayerThickness/2,0),
- "Layer 4", //its name
- logicLayer, //its logical volume
- physiDet, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
- physiLayer5 = new G4PVPlacement(0, //no rotation
- G4ThreeVector(0,
- +3*LayerThickness/2,0),
- "Layer 5", //its name
- logicLayer, //its logical volume
- physiDet, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
- physiLayer6 = new G4PVPlacement(0, //no rotation
- G4ThreeVector(0,
- +5*LayerThickness/2,0),
- "Layer 6", //its name
- logicLayer, //its logical volume
- physiDet, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
-
-
-
-
-
- /*
- //
- // Tub-Layer
- //
- solidTubLayer = new G4Box("Tub-Layer", //its name
- TubRad*NbOfTubs,TubRad,DetSizeX/2); //size
-
- logicTubLayer = new G4LogicalVolume(solidTubLayer, //its solid
- LayerMaterial, //its material
- "Tub Layer"); //its name
-
- G4RotationMatrix * rm = new G4RotationMatrix();
- G4double phi = 90*deg;
- rm->rotateY(phi);
-
- physiTubLayer = new G4PVPlacement(rm, //rotation
- G4ThreeVector(0,SphereRad,-0.1*cm),
- logicTubLayer, //its logical volume
- "Tub Layer", //its name
- logicLayer, //its mother volume
- false, //no boolean operation
- 0); //copy number
- */
- //
- // Tubs
- //
- solidTub=NULL; logicTub=NULL; physiTub=NULL;solidBoxIntersTub=NULL;
- solidTubSeg=NULL; logicTubSeg=NULL; physiTubSeg=NULL;
-
- if (TubRad > 0.)
- { solidTub = new G4Tubs("Tub", //its name
- 0,
- TubRad,
- DetSizeX/2,
- 0.*deg,
- 360.*deg);
-
- logicTub = new G4LogicalVolume(solidTub, //its solid
- TubMaterial, //its material
- "Tub"); //its name
-
- G4RotationMatrix * rm = new G4RotationMatrix();
- G4double phi = 90*deg;
- rm->rotateY(phi);
-
- for(int tubCopy=0; tubCopy<NbOfTubs; tubCopy++)
- {
- physiTub = new G4PVPlacement(rm, //rotation
- G4ThreeVector(0,SphereRad,
- -TubRad*NbOfTubs+TubRad*(1+2*tubCopy)-0.1*cm),
- logicTub, //its logical volume
- "Tub", //its name
- logicLayer, //its mother volume
- false, //no boolean operation
- tubCopy); //copy number
-
- }
-
-
- /*
- physiTub = new G4PVReplica("Tub", //its name
- logicTub, //its logical volume
- physiTubLayer, //its mother
- kXAxis, //axis of replication
- NbOfTubs, //number of replica
- 2*TubRad); //witdth of replica
- */
-
- //Tub segments
- solidBoxIntersTub = new G4Box("Tub segment", //its name
- TubRad,
- TubRad,
- DetSizeX/2);
-
- solidTubSeg = new G4IntersectionSolid("Tub segment",
- solidTub,
- solidBoxIntersTub,
- 0,
- G4ThreeVector(-0.4*cm,0,0));
-
- logicTubSeg = new G4LogicalVolume(solidTubSeg, //its solid
- TubMaterial, //its material
- "Tub segment"); //its name
-
- physiTubSeg = new G4PVPlacement(rm, //rotation
- G4ThreeVector(0,SphereRad,4.1*cm),
- logicTubSeg, //its logical volume
- "Tub segment", //its name
- logicLayer, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
- }
-
- //
- // Spheres
- //
- solidSphere=NULL; logicSphere=NULL; physiSphere=NULL;
-
- if (SphereRad > 0.)
- { solidSphere = new G4Sphere("Sphere",
- 0*cm,SphereRad,
- 0,360*deg,
- 0,180*deg);
-
- logicSphere = new G4LogicalVolume(solidSphere,
- SphereMaterial,
- "Sphere");
-
- G4double serie, element, Xposition, Yposition, Zposition;
- for(int copyNo=0;copyNo<40;copyNo++)
- {
- serie = int(copyNo/9);
- element = copyNo - serie * 9 + 1;
- Yposition = -TubRad;
-
- if (element<5)
- {
- Xposition = -2.55*cm + 2*SphereRad * element;
- Zposition = -4.0*cm + SphereRad + sqrt(3)*SphereRad*2*serie;
- }
- else
- {
- Xposition = -2.55*cm + SphereRad + (2*SphereRad)*(element-5);
- Zposition = -4.0*cm + SphereRad + sqrt(3)*SphereRad*(2*serie+1);
- }
-
- physiSphere = new G4PVPlacement(0, //rotation
- G4ThreeVector(Xposition,Yposition,Zposition),
- logicSphere, //its logical volume
- "Sphere", //its name
- logicLayer, //its mother volume
- false, //no boolean operation
- copyNo); //copy number
-
- }
-
- //Sphere segments
- solidBoxIntersSph = new G4Box("Sphere segment solid", //its name
- SphereRad,
- SphereRad,
- SphereRad);
-
-
- solidSphSeg1 = new G4IntersectionSolid("Sphere segment 1",
- solidSphere,
- solidBoxIntersSph,
- 0,
- G4ThreeVector(SphereRad,0,0));
-
- logicSphSeg1 = new G4LogicalVolume(solidSphSeg1,
- SphereMaterial,
- "Sphere segment 1");
-
- solidSphSeg2 = new G4IntersectionSolid("Sphere segment 2",
- solidSphere,
- solidBoxIntersSph,
- 0,
- G4ThreeVector(-0.4*cm,0,0));
-
- logicSphSeg2 = new G4LogicalVolume(solidSphSeg2,
- SphereMaterial,
- "Sphere segment 2");
-
- solidSphSeg3 = new G4IntersectionSolid("Sphere segment 3",
- solidSphere,
- solidBoxIntersSph,
- 0,
- G4ThreeVector(-0.9*cm,0,0));
-
- logicSphSeg3 = new G4LogicalVolume(solidSphSeg3,
- SphereMaterial,
- "Sphere segment 3");
-
- solidSphSeg4 = new G4IntersectionSolid("Sphere segment 4",
- solidSphere,
- solidBoxIntersSph,
- 0,
- G4ThreeVector(0,0,
- 8.0*cm-(9*sqrt(3)+2)*SphereRad));
-
- logicSphSeg4 = new G4LogicalVolume(solidSphSeg4,
- SphereMaterial,
- "Sphere segment 4");
-
-
- for(int s1=0;s1<5;s1++)
- {
- physiSphSeg1 = new G4PVPlacement(0, //no rotation
- G4ThreeVector(-DetSizeX/2,-TubRad,
- -4.0*cm+SphereRad+sqrt(3)*SphereRad*2*s1),
- logicSphSeg1, //its logical volume
- "Sphere segments 1", //its name
- logicLayer, //its mother volume
- false, //no boolean operation
- s1); //copy number
-
- physiSphSeg2 = new G4PVPlacement(0, //no rotation
- G4ThreeVector(2.45*cm,-TubRad,
- -4.0*cm+SphereRad+sqrt(3)*SphereRad*2*s1),
- logicSphSeg2, //its logical volume
- "Sphere segments 2", //its name
- logicLayer, //its mother volume
- false, //no boolean operation
- s1); //copy number
-
- physiSphSeg4 = new G4PVPlacement(0, //no rotation
- G4ThreeVector(
- -2.55*cm+SphereRad+2*SphereRad*s1,
- -TubRad,
- -4.0*cm+SphereRad+sqrt(3)*SphereRad*9),
- logicSphSeg4, //its logical volume
- "Sphere segments 4", //its name
- logicLayer, //its mother volume
- false, //no boolean operation
- s1); //copy number
- }
-
- for(int s2=0;s2<4;s2++)
- {
- physiSphSeg3 = new G4PVPlacement(0, //no rotation
- G4ThreeVector(2.95*cm,
- -TubRad,
- -4.0*cm+SphereRad+sqrt(3)*SphereRad*(2*s2+1)),
- logicSphSeg3, //its logical volume
- "Sphere segments 3", //its name
- logicLayer, //its mother volume
- false, //no boolean operation
- s2); //copy number
- }
- }
-
-
- //
- // Visualization attributes
- //
- logicWorld->SetVisAttributes (G4VisAttributes::Invisible);
- G4VisAttributes* simpleBoxVisAtt= new G4VisAttributes(G4Colour(1.0,1.0,1.0));
- simpleBoxVisAtt->SetVisibility(true);
- logicDet->SetVisAttributes(simpleBoxVisAtt);
-
- //
- //always return the physical World
- //
-
- return physiWorld;
-}
-
-
-
-/*
-void MyDetectorConstruction::SetMagField(G4double fieldValue)
-{
- //apply a global uniform magnetic field along Z axis
- FGeometryInit * ptrGeoInit = FGeometryInit::GetInstance();
- G4FieldManager* fieldMgr = ptrGeoInit->getFieldManager();
-
- if(magField) delete magField; //delete the existing magn field
-
- if(fieldValue!=0.) // create a new one if non nul
- { magField = new G4UniformMagField(G4ThreeVector(0.,0.,fieldValue));
- fieldMgr->SetDetectorField(magField);
- fieldMgr->CreateChordFinder(magField);
- } else {
- magField = NULL;
- fieldMgr->SetDetectorField(magField);
- }
- }
-
- */
+++ /dev/null
-# $Id$
-# Flugg tag $Name$
-
-# Id: GNUmakefile,v 2.1 1998/07/02 15:26:52 gunter Exp
-# --------------------------------------------------------------
-# GNUmakefile for examples module. Gabriele Cosmo, 06/04/98.
-# --------------------------------------------------------------
-
-name := mainT36
-G4TARGET := $(name)
-G4EXLIB := true
-
-ifndef FLUGGINSTALL
- FLUGGINSTALL = ../../..
-endif
-
-#the target of the first rule of the first makefile is the default goal
-#so "make" will do lib and bin targets.
-.PHONY: all
-all: lib bin
-
-include $(FLUGGINSTALL)/config/binmake.gmk
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-// MyDetectorConstruction.hh, 11/XII/1998 - Sara Vanini
-// GEANT4 tag $Name$
-//
-//
-
-#ifndef MyDetectorConstruction_h
-#define MyDetectorConstruction_h 1
-
-#include "G4VUserDetectorConstruction.hh"
-#include "globals.hh"
-
-class G4Box;
-class G4Tubs;
-class G4LogicalVolume;
-class G4VPhysicalVolume;
-class G4Material;
-class G4UniformMagField;
-//class T36EMCalorimeterSD;
-//class T36HADCalorimeterSD;
-
-class MyDetectorConstruction : public G4VUserDetectorConstruction
-{
- public:
-
- MyDetectorConstruction();
- ~MyDetectorConstruction();
- G4VPhysicalVolume* Construct();
- void SetMagField(G4double);
- void PrintCalorParameters();
-
-//for N02 (1+staies for the last HAD absorber):
- G4int GetNbOfLayers() {return (1+NbOfEMLayers+NbOfHADLayers);};
-
-
- G4int GetNbOfEMLayers() {return NbOfEMLayers;};
- G4int GetNbOfHADLayers() {return NbOfHADLayers;};
- G4double GetWorldSizeX() {return WorldSizeX;};
- G4double GetWorldSizeY() {return WorldSizeY;};
- G4double GetWorldSizeZ() {return WorldSizeZ;};
- G4double GetCalorSizeY() {return ModuleSizeY;};
- G4double GetCalorSizeZ() {return CalorSizeZ;};
- const G4VPhysicalVolume* GetEMphysiAbs() {return EMphysiAbsorber;};
- const G4VPhysicalVolume* GetEMphysiGap() {return EMphysiGap;};
- const G4VPhysicalVolume* GetHADphysiAbs() {return HADphysiAbsorber;};
- const G4VPhysicalVolume* GetlastHADAbsorber() {return lastHADphysiAbsorber;};
- const G4VPhysicalVolume* GetHADphysiGap() {return HADphysiGap;};
-
- const G4VPhysicalVolume* GetEMphysiAbMedShield() {return EMphysiAbMedShield;};
-
- const G4VPhysicalVolume* GetEMphysiSc() {return EMphysiSc;};
-
- const G4VPhysicalVolume* GetHADphysiAbMedShield() {return HADphysiAbMedShield;};
- const G4VPhysicalVolume* GetHADphysiAbExShield() {return HADphysiAbExShield;};
- const G4VPhysicalVolume* GetHADphysiSc() {return HADphysiSc;};
-
- const G4VPhysicalVolume* GetlastHADphysiAbs() {return lastHADphysiAbsorber;};
-
-
- private:
-
- G4double EMModuleThickness,EMModuleSizeZ;
- G4double HADModuleThickness,HADModuleSizeZ;
- G4Material* AbsorberMaterial,*GapMaterial,*defaultMaterial;
- G4Material* RodMaterial,*ScMaterial,*SpacerMaterial;
- G4Material* WLSMaterial,*frontPlateMaterial,*MedAbsMaterial,*ExAbsMaterial;
- G4Material* vacuum,*WLSMat,*Al,*MedLead,*ExLead;
- G4double AbsorberThickness,AbMedShieldThickness,AbExShieldThickness;
-
- G4double GapThickness;
-
- G4int NbOfEMLayers,NbOfHADLayers,NbOfModules;
- G4double LayerThickness;
-
- G4double ModuleSizeX,ModuleSizeY,ModuleSizeZ,CalorSizeZ;
- G4double EMWLSSizeX,EMWLSSizeZ,RodRad;
- G4double HADWLSSizeX,HADWLSSizeZ;
- G4double ScThickness,ScSizeY,AlSizeX,SpaSizeX;
- G4double SpacerThickness,SpacerSizeY;
-
- G4double WorldSizeX,WorldSizeY,WorldSizeZ;
-
- G4Box* solidWorld,*solidFrontVac; //pointer to the solid World
- G4LogicalVolume* logicWorld,*logicFrontVac; //pointer to the logical World
- G4VPhysicalVolume* physiWorld,*physiFrontVac; //pointer to the physical World
-
- G4Box* solidCalor,*solidModule,*solidEMWLS,*solidHADWLS,*solidSpace;
- G4LogicalVolume* logicCalor,*logicModule,*logicEMWLS,*logicHADWLS,*logicSpace;
- G4VPhysicalVolume* physiCalor,*physiModule,*physiSpace;
-
- G4Box* EMsolidLayer,*HADsolidLayer,*solidAl,*EMsolidModule,*HADsolidModule;
- G4LogicalVolume* EMlogicLayer,*HADlogicLayer,*logicAl,*EMlogicModule,*HADlogicModule;
- G4VPhysicalVolume* EMphysiLayer,*HADphysiLayer,*physiAl,*EMphysiModule,*HADphysiModule;
-
- G4Box* EMsolidAbsorber,*HADsolidAbsorber,*lastHADsolidAbsorber;
- G4LogicalVolume* EMlogicAbsorber,*HADlogicAbsorber,*lastHADlogicAbsorber;
- G4VPhysicalVolume* EMphysiAbsorber,*HADphysiAbsorber,*lastHADphysiAbsorber;
-
- G4Box* EMsolidAbMedShield,*HADsolidAbMedShield,*EMsolidAbExShield,*HADsolidAbExShield;
- G4LogicalVolume* EMlogicAbMedShield,*HADlogicAbMedShield,*EMlogicAbExShield,*HADlogicAbExShield;
- G4VPhysicalVolume* EMphysiAbMedShield,*HADphysiAbMedShield,*lastHADphysiAbMedShield;
- G4VPhysicalVolume* EMphysiAbExShield,*HADphysiAbExShield,*lastHADphysiAbExShield;
-
- G4Box* EMsolidGap,*HADsolidGap,*EMsolidSc,*HADsolidSc;
- G4LogicalVolume* EMlogicGap,*HADlogicGap,*EMlogicSc,*HADlogicSc;
- G4VPhysicalVolume* EMphysiGap,*HADphysiGap,*EMphysiSc,*HADphysiSc;
-
- G4Box* EMsolidSpacer,*HADsolidSpacer;
- G4Tubs* EMsolidRod,*HADsolidRod;
- G4LogicalVolume* EMlogicRod,*HADlogicRod,*EMlogicSpacer,*HADlogicSpacer;
-
- G4UniformMagField* magField; //pointer to the magnetic field
-
- //pointer to the sensitive detectors - EM and HAD
- // T36EMCalorimeterSD* EMcalorimeterSD;
- // T36HADCalorimeterSD* HADcalorimeterSD;
-
-//for N02
-// ExN02CalorimeterSD* calorimeterSD;
-
- private:
-
- void DefineMaterials();
- void ComputeCalorParameters();
- G4VPhysicalVolume* ConstructCalorimeter();
-};
-
-
-
-inline void MyDetectorConstruction::ComputeCalorParameters()
-{
- // Compute derived parameters of the calorimeter
-
- LayerThickness = AbsorberThickness + GapThickness;
- EMModuleThickness = NbOfEMLayers*LayerThickness;
- HADModuleThickness = NbOfHADLayers*LayerThickness;
-
- EMWLSSizeX = EMModuleThickness + HADModuleThickness + AbsorberThickness;
- HADWLSSizeX = HADModuleThickness + AbsorberThickness;
-
- ModuleSizeZ = 2*EMWLSSizeZ + EMModuleSizeZ;
- CalorSizeZ = NbOfModules*ModuleSizeZ;
- ModuleSizeX = AlSizeX + (NbOfEMLayers+NbOfHADLayers)*LayerThickness
- + AbsorberThickness + SpaSizeX;
-}
-
-#endif
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-//
-// Parameterisation for t36 EM and HAD layers
-//
-
-#ifndef T36LayerParam_H
-#define T36LayerParam_H 1
-
-#include "globals.hh"
-#include "G4VPVParameterisation.hh"
-
-class G4VPhysicalVolume;
-class G4Box;
-
-class T36LayerParam : public G4VPVParameterisation
-{
- public:
- T36LayerParam(G4double LayerThickness, G4double NbOfLayers);
- ~T36LayerParam();
- void ComputeTransformation
- (const G4int copyNo,G4VPhysicalVolume *physVol) const;
-
- private:
-
- G4double fLayerThickness;
- G4double fNbOfLayers;
-};
-
-#endif
-
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-//
-// Parameterisation for t36 module
-//
-
-#ifndef T36ModuleParam_H
-#define T36ModuleParam_H 1
-
-#include "globals.hh"
-#include "G4VPVParameterisation.hh"
-
-class G4VPhysicalVolume;
-class G4Box;
-
-class T36ModuleParam : public G4VPVParameterisation
-{
- public:
- T36ModuleParam(G4double ModuleSizeZ, G4double NbOfModules);
- ~T36ModuleParam();
- void ComputeTransformation
- (const G4int copyNo,G4VPhysicalVolume *physVol) const;
-
- private:
-
- G4double fModuleSizeZ;
- G4double fNbOfModules;
-};
-
-#endif
-
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#include "FGeometryInit.hh"
-#include "MyDetectorConstruction.hh"
-
-#define flukam flukam_
-
-extern "C" void flukam(const G4int & GeoFlag);
-
-int main() {
-
- FGeometryInit* theFGeometryInit = FGeometryInit::GetInstance();
-
- theFGeometryInit
- ->setDetConstruction(new MyDetectorConstruction());
-
-//flag for geometry:
-// 1 for GEANT4
-// 0 for FLUKA
-// 2 for Rubia
- const G4int flag = 1;
-
-//call fortran
- flukam(flag);
-
-//end
- return 0;
-}
-
-
-
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-//
-// MyDetectorConstruction.hh, 3/III/1998 - Sara Vanini
-// T36 calorimeter with parametric volumes insted of replicans
-//
-//
-
-//#include "T36EMCalorimeterSD.hh"
-//#include "T36HADCalorimeterSD.hh"
-#include "MyDetectorConstruction.hh"
-//#include "T36ModuleParam.hh"
-#include "T36LayerParam.hh"
-
-#include "G4Material.hh"
-#include "G4Box.hh"
-#include "G4Tubs.hh"
-#include "G4LogicalVolume.hh"
-#include "G4PVPlacement.hh"
-#include "G4PVReplica.hh"
-#include "G4PVParameterised.hh"
-#include "G4VPVParameterisation.hh"
-//#include "G4SDManager.hh"
-//#include "G4RunManager.hh"
-#include "G4UniformMagField.hh"
-#include "G4FieldManager.hh"
-#include "G4TransportationManager.hh"
-#include "G4VisAttributes.hh"
-#include "G4Colour.hh"
-#include "G4ios.hh"
-
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
-
-
-MyDetectorConstruction::MyDetectorConstruction()
-:solidWorld(NULL),logicWorld(NULL),physiWorld(NULL),
- solidCalor(NULL),logicCalor(NULL),physiCalor(NULL),
- solidModule(NULL),logicModule(NULL),physiModule(NULL),
- EMsolidModule(NULL),EMlogicModule(NULL),EMphysiModule(NULL),
- HADsolidModule(NULL),HADlogicModule(NULL),HADphysiModule(NULL),
- solidFrontVac(NULL),logicFrontVac(NULL),physiFrontVac(NULL),
-
- solidEMWLS(NULL),logicEMWLS(NULL),solidHADWLS(NULL),logicHADWLS(NULL),
- solidAl(NULL),logicAl(NULL),physiAl(NULL),
- solidSpace(NULL),logicSpace(NULL),physiSpace(NULL),
- EMsolidLayer(NULL),EMlogicLayer(NULL),EMphysiLayer(NULL),
- HADsolidLayer(NULL),HADlogicLayer(NULL),HADphysiLayer(NULL),
-
- EMsolidAbsorber(NULL),EMlogicAbsorber(NULL),EMphysiAbsorber(NULL),
- HADsolidAbsorber(NULL),HADlogicAbsorber(NULL),HADphysiAbsorber(NULL),
- lastHADsolidAbsorber(NULL),lastHADlogicAbsorber(NULL),lastHADphysiAbsorber(NULL),
- lastHADphysiAbMedShield(NULL),lastHADphysiAbExShield(NULL),
- EMsolidAbMedShield(NULL),EMlogicAbMedShield(NULL),EMphysiAbMedShield(NULL),
- EMsolidAbExShield(NULL),EMlogicAbExShield(NULL),EMphysiAbExShield(NULL),
- HADsolidAbMedShield(NULL),HADlogicAbMedShield(NULL),HADphysiAbMedShield(NULL),
- HADsolidAbExShield(NULL),HADlogicAbExShield(NULL),HADphysiAbExShield(NULL),
-
-
- EMsolidGap(NULL),EMlogicGap(NULL),EMphysiGap(NULL),
- HADsolidGap(NULL),HADlogicGap(NULL),HADphysiGap(NULL),
- EMsolidSc(NULL),EMlogicSc(NULL),EMphysiSc(NULL),
- HADsolidSc(NULL),HADlogicSc(NULL),HADphysiSc(NULL),
- AbsorberMaterial(NULL),GapMaterial(NULL),ScMaterial(NULL),WLSMaterial(NULL),
- RodMaterial(NULL),defaultMaterial(NULL),SpacerMaterial(NULL),
- frontPlateMaterial(NULL),MedAbsMaterial(NULL),ExAbsMaterial(NULL),
- EMsolidRod(NULL),HADsolidRod(NULL),EMlogicRod(NULL),HADlogicRod(NULL),
- EMsolidSpacer(NULL),HADsolidSpacer(NULL),EMlogicSpacer(NULL),HADlogicSpacer(NULL)
-
-{
- WorldSizeX = 4000*mm;
- WorldSizeY = 4000*mm;
- WorldSizeZ = 4000*mm;
-
- // default parameter values of the calorimeter
- // EMWLSSizeX = 1103.5*mm; //copre tutta la torre
- EMWLSSizeZ = 5*mm;
- // HADWLSSizeX = 887.5*mm; //come input fluka: dal layer 17 per tutta la torre
- HADWLSSizeZ = 3.5*mm;
- ScSizeY = 218*mm;
- ScThickness = 2.5*mm;
- RodRad = 1*mm;
- AlSizeX = 20*mm;
- SpaSizeX = 1*mm;
- EMModuleSizeZ = 218*mm;
- HADModuleSizeZ = 211*mm;
- SpacerThickness = 3.5*mm;
- SpacerSizeY = 21*mm;
-
-
- AbsorberThickness = 10.*mm;
- AbMedShieldThickness = 8.*mm;
- AbExShieldThickness = 1.*mm;
- GapThickness = 3.5*mm;
- NbOfEMLayers = 16;
- NbOfHADLayers = 65;
- EMModuleSizeZ = 218*mm;
- HADModuleSizeZ = 211*mm;
-
- ModuleSizeY = 700*mm;
- NbOfModules = 3;
-}
-
-
-MyDetectorConstruction::~MyDetectorConstruction()
-{}
-
-
-G4VPhysicalVolume* MyDetectorConstruction::Construct()
-{
- DefineMaterials();
- return ConstructCalorimeter();
-}
-
-
-void MyDetectorConstruction::DefineMaterials()
-{
- //This function illustrates the possible ways to define materials
-
-G4String name, symbol; //a=mass of a mole;
-G4double a, z, density; //z=mean number of protons;
-G4int iz, n; //iz=number of protons in an isotope;
- // n=number of nucleons in an isotope;
-
-G4int ncomponents, natoms;
-G4double abundance, fractionmass;
-G4double temperature, pressure;
-
-//
-// define Elements
-//
-
-a = 1.0079*g/mole;
-G4Element* H = new G4Element(name="Hydrogen",symbol="H" , z= 1., a);
-
-a = 12.01*g/mole;
-G4Element* C = new G4Element(name="Carbon" ,symbol="C" , z= 6., a);
-
-a = 14.007*g/mole;
-G4Element* N = new G4Element(name="Nitrogen",symbol="N" , z= 7., a);
-
-a = 15.999*g/mole;
-G4Element* O = new G4Element(name="Oxygen" ,symbol="O" , z= 8., a);
-
-a = 39.948*g/mole;
-G4Element* Ar = new G4Element(name="Argon" ,symbol="Ar" , z= 18., a);
-
-a = 207.19*g/mole;
-G4Element* Pb = new G4Element(name="Lead" ,symbol="Pb" , z= 82., a);
-
-a = 121.75*g/mole;
-G4Element* Sb = new G4Element(name="Antimony" ,symbol="Sb" , z= 51., a);
-
-
-//
-// define simple materials
-//
-
-density = universe_mean_density; //from PhysicalConstants.h
-pressure = 3.e-18*pascal;
-temperature = 2.73*kelvin;
-G4Material* vacuum = new G4Material(name="Galactic", z=1., a=1.01*g/mole, density,
- kStateGas,temperature,pressure);
-
-density = 2.700*g/cm3;
-a = 26.982*g/mole;
-G4Material* Al = new G4Material(name="Aluminium", z=13., a, density);
-
-// il numero di atomi e` stato normalizzato al numero di atomi della molecola!
-// posso usare la formula chimica minimale, quindi moltiplicare per un
-// fattore comune fino ad ottenere interi.
-
-density = 1.170*g/cm3;
-G4Material* WLSMat = new G4Material(name="PMMAWLS", density, ncomponents=3);
-WLSMat->AddElement(H, natoms=100);
-WLSMat->AddElement(C, natoms=57);
-WLSMat->AddElement(O, natoms=57);
-
-//N.B. le frazioni di massa sono state normalizzate ad 1 rispetto all`input di fluka!
-density = 11.3*g/cm3;
-G4Material* ExLead = new G4Material(name="LeadSB", density, ncomponents=2);
-ExLead->AddElement(Pb, fractionmass=0.96);
-ExLead->AddElement(Sb, fractionmass=0.04);
-
-density = 11.3*g/cm3;
-G4Material* MedLead = new G4Material(name="LeadSB0", density, ncomponents=2);
-MedLead->AddElement(Pb, fractionmass=0.96);
-MedLead->AddElement(Sb, fractionmass=0.04);
-
-density = 0.001225*g/cm3;
-G4Material* Air = new G4Material(name="Air", density, ncomponents=3);
-Air->AddElement(N, fractionmass=0.7555795);
-Air->AddElement(O, fractionmass=0.23158806);
-Air->AddElement(Ar, fractionmass=0.012832444);
-
-density = 1.044*g/cm3;
-G4Material* SciMat = new G4Material(name="SCSN38", density, ncomponents=2);
-SciMat->AddElement(H, natoms=1);
-SciMat->AddElement(C, natoms=1);
-
-
-//G4cout << *(G4Material::GetMaterialTable()) << endl;
-
- //default materials of the calorimeter
- GapMaterial = Air;
- defaultMaterial = Air;
- SpacerMaterial = Al;
- RodMaterial = SciMat;
- ScMaterial = SciMat;
-
- WLSMaterial = WLSMat;
- frontPlateMaterial = Al;
- MedAbsMaterial = MedLead;
- ExAbsMaterial = ExLead;
-}
-
-
-
-
-G4VPhysicalVolume* MyDetectorConstruction::ConstructCalorimeter()
-{
- // complete the Calor parameters definition
- ComputeCalorParameters();
-
- //
- // World
- //
- solidWorld = new G4Box("World", //its name
- WorldSizeX/2,WorldSizeY/2,WorldSizeZ/2); //its size
-
- logicWorld = new G4LogicalVolume(solidWorld, //its solid
- defaultMaterial, //vacuum
- "World"); //its name
-
- physiWorld = new G4PVPlacement(0, //no rotation
- G4ThreeVector(), //at (0,0,0)
- "World", //its name
- logicWorld, //its logical volume
- NULL, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
- //
- // Calorimeter
- //
- solidCalor = new G4Box("Calorimeter", //its name
- ModuleSizeX/2,ModuleSizeY/2,CalorSizeZ/2);//size
-
- logicCalor = new G4LogicalVolume(solidCalor, //its solid
- defaultMaterial, //its material
- "Calorimeter"); //its name
-
- physiCalor = new G4PVPlacement(0, //no rotation
- G4ThreeVector(ModuleSizeX/2-AlSizeX-SpaSizeX,0.,0.),
- "Calorimeter", //its name
- logicCalor, //its logical volume
- physiWorld, //its mother volume
- false, //no boolean operation
- 0); //copy number
- //
- // Module
- //
- solidModule = new G4Box("Module", //its name
- ModuleSizeX/2,ModuleSizeY/2,ModuleSizeZ/2);//size
-
- logicModule = new G4LogicalVolume(solidModule, //its solid
- defaultMaterial, //its material
- "Module"); //its name
-
-
- for(int iMod=0;iMod<NbOfModules;iMod++)
- {
- new G4PVPlacement(0,
- G4ThreeVector(0.,0.,
- -ModuleSizeZ/2*(NbOfModules-1)+iMod*ModuleSizeZ),
- "Module",
- logicModule,
- physiCalor,
- false,
- iMod);
- }
-
-
- //
- // EMWLS
- //
- solidEMWLS = new G4Box("EMWLS", //its name
- EMWLSSizeX/2,ModuleSizeY/2,EMWLSSizeZ/2);//size
-
- logicEMWLS = new G4LogicalVolume(solidEMWLS, //its solid
- WLSMaterial, //its material
- "EMWLS"); //its name
- //right and left side of module
- for(int t=0; t<2; t++)
- {
- int b=2*t-1;
- new G4PVPlacement(0, //no rotation
- G4ThreeVector(-ModuleSizeX/2+EMWLSSizeX/2+AlSizeX+SpaSizeX,
- 0.,b*(-ModuleSizeZ/2+EMWLSSizeZ/2)), //
- logicEMWLS, //its logical volume
- "EMWLS", //its name
- logicModule, //its mother volume
- false, //no boolean operation
- t); //copy number
- }
-
- //
- // HADWLS
- //
- solidHADWLS = new G4Box("HADWLS", //its name
- HADWLSSizeX/2,ModuleSizeY/2,HADWLSSizeZ/2);//size
-
- logicHADWLS = new G4LogicalVolume(solidHADWLS, //its solid
- WLSMaterial, //its material
- "HADWLS"); //its name
- //left and right side of module
- for(int u=0; u<2; u++)
- {
- int c=2*u-1;
- new G4PVPlacement(0, //no rotation
- G4ThreeVector(ModuleSizeX/2-HADWLSSizeX/2,
- 0.,c*(ModuleSizeZ/2-EMWLSSizeZ-HADWLSSizeZ/2)),
- logicHADWLS, //its logical volume
- "HADWLS", //its name
- logicModule, //its mother volume
- false, //no boolean operation
- u); //copy number
- }
-
-
- //
- // Al front plate
- //
- solidAl = new G4Box("AlFrontPlate", //its name
- AlSizeX/2,ModuleSizeY/2,ModuleSizeZ/2);//size
-
- logicAl = new G4LogicalVolume(solidAl, //its solid
- frontPlateMaterial, //its material
- "AlFrontPlate"); //its name
-
- physiAl = new G4PVPlacement(0, //no rotation
- G4ThreeVector(-ModuleSizeX/2+AlSizeX/2,0.,0.),
- logicAl, //its logical volume
- "AlFrontPlate", //its name
- logicModule, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
- //
- // front space with vacuum
- //
- solidSpace = new G4Box("VacFrontSpace", //its name
- SpaSizeX/2,ModuleSizeY/2,ModuleSizeZ/2);//size
-
- logicSpace = new G4LogicalVolume(solidSpace, //its solid
- defaultMaterial, //its material
- "VacFrontSpace"); //its name
-
- physiSpace = new G4PVPlacement(0, //no rotation
- G4ThreeVector(-ModuleSizeX/2+AlSizeX+SpaSizeX/2,0.,0.),
- logicSpace, //its logical vo
- "VacFrontSpace", //its name
- logicModule, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
-
-
-
-
- // *******************************
- // ******** EMModule *******
- // *******************************
-
- EMsolidModule = new G4Box("EMModule", //its name
- EMModuleThickness/2,ModuleSizeY/2,EMModuleSizeZ/2);//size
-
- EMlogicModule = new G4LogicalVolume(EMsolidModule, //its solid
- defaultMaterial, //its material
- "EMModule"); //its name
-
- EMphysiModule = new G4PVPlacement(0, //no rotation
- G4ThreeVector(-ModuleSizeX/2+AlSizeX+SpaSizeX+EMModuleThickness/2,
- 0.,0.),
- EMlogicModule, //its logical volume
- "EMModule", //its name
- logicModule, //its mother volume
- false, //no boolean operation
- 0); //copy number
- //
- // EMLayer
- //
- EMsolidLayer = new G4Box("EMLayer", //its name
- LayerThickness/2,ModuleSizeY/2,EMModuleSizeZ/2); //size
-
- EMlogicLayer = new G4LogicalVolume(EMsolidLayer, //its solid
- defaultMaterial, //its material
- "EMLayer"); //its name
-
- G4VPVParameterisation * EMLayersParam =
- new T36LayerParam(LayerThickness,NbOfEMLayers);
-
- EMphysiLayer = new G4PVParameterised("EMLayer", //its name
- EMlogicLayer, //its logical volume
- EMphysiModule, //its mother
- kXAxis, //axis of param
- NbOfEMLayers, //number of param
- EMLayersParam); //param
-
- //
- // Pb EMAbsorber
- //
- EMsolidAbsorber = new G4Box("EMAbsorber", //its name
- AbsorberThickness/2,ModuleSizeY/2,EMModuleSizeZ/2);
-
- EMlogicAbsorber = new G4LogicalVolume(EMsolidAbsorber, //its solid
- defaultMaterial, //its material
- "EMAbsorber"); //its name
-
- EMphysiAbsorber = new G4PVPlacement(0, //no rotation
- G4ThreeVector(-GapThickness/2,0.,0.), //its position
- "EMAbsorber", //its name
- EMlogicAbsorber, //its logical volume
- EMphysiLayer, //its mother
- false, //no boulean operat
- 0); //copy number
-
- //
- // Pb EMAbsorber Medium Shield
- //
- EMsolidAbMedShield = new G4Box("EMAbMedShield", //its name
- AbMedShieldThickness/2,ModuleSizeY/2,EMModuleSizeZ/2);
-
- EMlogicAbMedShield = new G4LogicalVolume(EMsolidAbMedShield, //its solid
- MedAbsMaterial, //its material
- "EMAbMedShield"); //its name
-
- EMphysiAbMedShield = new G4PVPlacement(0, //no rotation
- G4ThreeVector(0.,0.,0.), //its position
- "EMAbMedShield", //its name
- EMlogicAbMedShield, //its logical volume
- EMphysiAbsorber, //its mother
- false, //no boulean operat
- 0); //copy number
-
- //
- // 2 Pb EMAbsorber External Shields
- //
- EMsolidAbExShield = new G4Box("EMAbExShield", //its name
- AbExShieldThickness/2,ModuleSizeY/2,EMModuleSizeZ/2);
-
- EMlogicAbExShield = new G4LogicalVolume(EMsolidAbExShield, //its solid
- ExAbsMaterial, //its material
- "EMAbExShield"); //its name
-
-
- for(int s=0; s<2; s++)
- {
- EMphysiAbExShield = new G4PVPlacement(0, //no rotation
- G4ThreeVector((2*s-1)*(AbsorberThickness/2-AbExShieldThickness/2),0.,0.), //its position
- "EMAbExShield", //its name
- EMlogicAbExShield, //its logical volume
- EMphysiAbsorber, //its mother
- false, //no boulean operat
- s); //copy number
-
- }
-
-
- //
- // EMGap
- //
- EMsolidGap = new G4Box("EMGap",
- GapThickness/2,ModuleSizeY/2,EMModuleSizeZ/2);
-
- EMlogicGap = new G4LogicalVolume(EMsolidGap,
- GapMaterial,
- "EMGap");
-
- EMphysiGap = new G4PVPlacement(0, //no rotation
- G4ThreeVector(AbsorberThickness/2,0.,0.), //its position
- "EMGap", //its name
- EMlogicGap, //its logical volume
- EMphysiLayer, //its mother
- false, //no boulean operat
- 0); //copy number
-
-
- //
- // 3 EMScPlates
- //
- EMsolidSc = new G4Box("EMSc",
- ScThickness/2,ScSizeY/2,EMModuleSizeZ/2);
-
- EMlogicSc = new G4LogicalVolume(EMsolidSc,
- ScMaterial,
- "EMSc");
-
-
- for(int i=0; i<+3; i++)
- {
- EMphysiSc = new G4PVPlacement(0, //no rotation
- G4ThreeVector(0.,(i-1)*(ScSizeY+2*RodRad),0.), //its position
- "EMSc", //its name
- EMlogicSc, //its logical volume
- EMphysiGap, //its mother
- false, //no boulean operat
- i);
- }
-
-
- //
- // EMRod
- //
- EMsolidRod = new G4Tubs("EMRod",
- 0.,RodRad,EMModuleSizeZ/2,0*deg,360*deg);
-
- EMlogicRod = new G4LogicalVolume(EMsolidRod,
- RodMaterial,
- "EMRod");
-
- for(int h=0; h<2; h++)
- {
- new G4PVPlacement(0, //no rotation
- G4ThreeVector(0.,(2*h-1)*(ScSizeY/2+RodRad),0.), //its position
- "EMRod", //its name
- EMlogicRod, //its logical volume
- EMphysiGap, //its mother
- false, //no boulean operat
- h);
- }
-
- //
- // EMSpacer
- //
- EMsolidSpacer = new G4Box("EMSpacer",
- SpacerThickness/2,SpacerSizeY/2,EMModuleSizeZ/2);
-
- EMlogicSpacer = new G4LogicalVolume(EMsolidSpacer,
- SpacerMaterial,
- "EMSpacer");
-
- for(int k=0; k<2; k++)
- {
- int v=2*k-1;
- new G4PVPlacement(0, //no rotation
- G4ThreeVector(0.,v*(3*ScSizeY/2+2*RodRad+SpacerSizeY/2),0.), //its position
- "EMSpacer", //its name
- EMlogicSpacer, //its logical volume
- EMphysiGap, //its mother
- false, //no boulean operat
- k);
- }
-
-
-
- // *******************************
- // ******** HADModule *******
- // *******************************
-
- HADsolidModule = new G4Box("HADModule", //its name
- HADModuleThickness/2,ModuleSizeY/2,HADModuleSizeZ/2);//size
-
- HADlogicModule = new G4LogicalVolume(HADsolidModule, //its solid
- defaultMaterial, //its material
- "HADModule"); //its name
-
- HADphysiModule = new G4PVPlacement(0, //no rotation
- G4ThreeVector(ModuleSizeX/2-HADModuleThickness/2-
- AbsorberThickness,0.,0.), //at (0,0,0)
- HADlogicModule, //its logical volume
- "HADModule", //its name
- logicModule, //its mother volume
- false, //no boolean operation
- 0); //copy number
-
- //
- // HADLayer
- //
- HADsolidLayer = new G4Box("HADLayer", //its name
- LayerThickness/2,ModuleSizeY/2,HADModuleSizeZ/2); //size
-
- HADlogicLayer = new G4LogicalVolume(HADsolidLayer, //its solid
- defaultMaterial, //its material
- "HADLayer"); //its name
-
- G4VPVParameterisation * HADLayersParam =
- new T36LayerParam(LayerThickness,NbOfHADLayers);
-
- HADphysiLayer = new G4PVParameterised("HADLayer", //its name
- HADlogicLayer, //its logical volume
- HADphysiModule, //its mother
- kXAxis, //axis of param
- NbOfHADLayers, //number of param
- HADLayersParam); //param
- //
- // Pb HADAbsorber
- //
- HADsolidAbsorber = new G4Box("HADAbsorber", //its name
- AbsorberThickness/2,ModuleSizeY/2,HADModuleSizeZ/2);
-
- HADlogicAbsorber = new G4LogicalVolume(HADsolidAbsorber, //its solid
- defaultMaterial, //its material
- "HADAbsorber"); //its name
-
- HADphysiAbsorber = new G4PVPlacement(0, //no rotation
- G4ThreeVector(-GapThickness/2,0.,0.), //its position
- "HADAbsorber", //its name
- HADlogicAbsorber, //its logical volume
- HADphysiLayer, //its mother
- false, //no boulean operat
- 0); //copy number
-
-
- //
- // Pb HADAbsorber Medium Shield
- //
- HADsolidAbMedShield = new G4Box("HADAbMedShield", //its name
- AbMedShieldThickness/2,ModuleSizeY/2,HADModuleSizeZ/2);
-
- HADlogicAbMedShield = new G4LogicalVolume(HADsolidAbMedShield, //its solid
- MedAbsMaterial, //its material
- "HADAbMedShield"); //its name
-
- HADphysiAbMedShield = new G4PVPlacement(0, //no rotation
- G4ThreeVector(0.,0.,0.), //its position
- "HADAbMedShield", //its name
- HADlogicAbMedShield, //its logical volume
- HADphysiAbsorber, //its mother
- false, //no boulean operat
- 0); //copy number
-
- //
- // 2 Pb HADAbsorber External Shields
- //
- HADsolidAbExShield = new G4Box("HADAbExShield", //its name
- AbExShieldThickness/2,ModuleSizeY/2,HADModuleSizeZ/2);
-
- HADlogicAbExShield = new G4LogicalVolume(HADsolidAbExShield, //its solid
- ExAbsMaterial, //its material
- "HADAbExShield"); //its name
-
- for(int hs=0; hs<2; hs++)
- {
- HADphysiAbExShield = new G4PVPlacement(0, //no rotation
- G4ThreeVector((2*hs-1)*(AbsorberThickness/2-AbExShieldThickness/2),0.,0.), //its position
- "HADAbExShield", //its name
- HADlogicAbExShield, //its logical volume
- HADphysiAbsorber, //its mother
- false, //no boulean operat
- hs); //copy number
-
- }
-
- //
- // last Pb HADAbsorber
- //
-
-
- lastHADsolidAbsorber = new G4Box("lastHADAbsorber", //its name
- AbsorberThickness/2,ModuleSizeY/2,HADModuleSizeZ/2);
-
- lastHADlogicAbsorber = new G4LogicalVolume(lastHADsolidAbsorber, //its solid
- defaultMaterial, //its material
- "lastHADAbsorber"); //its name
-
- lastHADphysiAbsorber = new G4PVPlacement(0, //no rotation
- G4ThreeVector(ModuleSizeX/2-AbsorberThickness/2,0.,0.), //its position
- lastHADlogicAbsorber, //its logical volume
- "lastHADAbsorber", //its name
- logicModule, //its mother
- false, //no boulean operat
- 0); //copy number
-
- //
- // last Pb HADAbsorber Medium Shield
- //
- lastHADphysiAbMedShield = new G4PVPlacement(0, //no rotation
- G4ThreeVector(0.,0.,0.), //its position
- "lastHADAbMedShield", //its name
- HADlogicAbMedShield, //its logical volume
- lastHADphysiAbsorber, //its mother
- false, //no boulean operat
- 0); //copy number
-
- //
- // 2 last Pb HADAbsorber External Shields
- //
- for(int ls=0; ls<2; ls++)
- {
- lastHADphysiAbMedShield = new G4PVPlacement(0, //no rotation
- G4ThreeVector((2*ls-1)*(AbsorberThickness/2-AbExShieldThickness/2),0.,0.), //its position
- "lastHADAbExShield", //its name
- HADlogicAbExShield, //its logical volume
- lastHADphysiAbsorber, //its mother
- false, //no boulean operat
- s); //copy number
-
- }
-
-
- //
- // HADGap
- //
- HADsolidGap = new G4Box("HADGap",
- GapThickness/2,ModuleSizeY/2,HADModuleSizeZ/2);
-
- HADlogicGap = new G4LogicalVolume(HADsolidGap,
- GapMaterial,
- "HADGap");
-
- HADphysiGap = new G4PVPlacement(0, //no rotation
- G4ThreeVector(AbsorberThickness/2,0.,0.), //its position
- "HADGap", //its name
- HADlogicGap, //its logical volume
- HADphysiLayer, //its mother
- false, //no boulean operat
- 0); //copy number
-
-
- //
- // 3 HADScPlates
- //
- HADsolidSc = new G4Box("HADSc",
- ScThickness/2,ScSizeY/2,HADModuleSizeZ/2);
-
- HADlogicSc = new G4LogicalVolume(HADsolidSc,
- ScMaterial,
- "HADSc");
-
-
- for(int m=0; m<+3; m++)
- {
- HADphysiSc = new G4PVPlacement(0, //no rotation
- G4ThreeVector(0.,(m-1)*(ScSizeY+2*RodRad),0.), //its position
- "HADSc", //its name
- HADlogicSc, //its logical volume
- HADphysiGap, //its mother
- false, //no boulean operat
- m);
- }
-
-
- //
- // HADRod
- //
- HADsolidRod = new G4Tubs("HADRod",
- 0.,RodRad,HADModuleSizeZ/2,0*deg,360*deg);
-
- HADlogicRod = new G4LogicalVolume(HADsolidRod,
- RodMaterial,
- "HADRod");
-
- for(int q=0; q<2; q++)
- {
- new G4PVPlacement(0, //no rotation
- G4ThreeVector(0.,(2*q-1)*(ScSizeY/2+RodRad),0.), //its position
- "HADRod", //its name
- HADlogicRod, //its logical volume
- HADphysiGap, //its mother
- false, //no boulean operat
- q);
- }
-
- //
- // HADSpacer
- //
- HADsolidSpacer = new G4Box("HADSpacer",
- SpacerThickness/2,SpacerSizeY/2,HADModuleSizeZ/2);
-
- HADlogicSpacer = new G4LogicalVolume(HADsolidSpacer,
- SpacerMaterial,
- "HADSpacer");
-
- for(int z=0; z<2; z++)
- {
- int y=2*z-1;
- new G4PVPlacement(0, //no rotation
- G4ThreeVector(0.,y*(3*ScSizeY/2+2*RodRad+SpacerSizeY/2),0.), //its position
- "HADSpacer", //its name
- HADlogicSpacer, //its logical volume
- HADphysiGap, //its mother
- false, //no boulean operat
- z);
- }
- /*
- //
- // Sensitive EM and HAD Detectors: Absorber and Gap
- //
- G4SDManager* SDman = G4SDManager::GetSDMpointer();
- SDman->SetVerboseLevel(1);
-
- //for N02: N.B. non basta dichiarare la "madre" volume sensibile, bisogna
- // dichiarare anche ogni figlio!!!
-
- if(!EMcalorimeterSD)
- {
- G4String EMcaloSDname = "/mydet/EMcalorimeter";
- EMcalorimeterSD = new T36EMCalorimeterSD(EMcaloSDname,this);
- SDman->AddNewDetector( EMcalorimeterSD );
- }
-
- if (EMlogicAbsorber)
- EMlogicAbsorber->SetSensitiveDetector(EMcalorimeterSD);
- if (EMlogicAbExShield)
- EMlogicAbExShield->SetSensitiveDetector(EMcalorimeterSD);
- if (EMlogicAbMedShield)
- EMlogicAbMedShield->SetSensitiveDetector(EMcalorimeterSD);
-
- if (EMlogicGap)
- EMlogicGap->SetSensitiveDetector(EMcalorimeterSD);
- if (EMlogicSc)
- EMlogicSc->SetSensitiveDetector(EMcalorimeterSD);
- if (EMlogicRod)
- EMlogicRod->SetSensitiveDetector(EMcalorimeterSD);
- if (EMlogicSpacer)
- EMlogicSpacer->SetSensitiveDetector(EMcalorimeterSD);
-
- if(!HADcalorimeterSD)
- {
- G4String HADcaloSDname = "/mydet/HADcalorimeter";
- HADcalorimeterSD = new T36HADCalorimeterSD(HADcaloSDname,this);
- SDman->AddNewDetector( HADcalorimeterSD );
- }
-
- if (HADlogicAbsorber)
- HADlogicAbsorber->SetSensitiveDetector(HADcalorimeterSD);
- if (HADlogicAbExShield)
- HADlogicAbExShield->SetSensitiveDetector(HADcalorimeterSD);
- if (HADlogicAbMedShield)
- HADlogicAbMedShield->SetSensitiveDetector(HADcalorimeterSD);
-
- if (HADlogicGap)
- HADlogicGap->SetSensitiveDetector(HADcalorimeterSD);
- if (HADlogicSc)
- HADlogicSc->SetSensitiveDetector(HADcalorimeterSD);
- if (HADlogicRod)
- HADlogicRod->SetSensitiveDetector(HADcalorimeterSD);
- if (HADlogicSpacer)
- HADlogicSpacer->SetSensitiveDetector(HADcalorimeterSD);
-
- if (lastHADlogicAbsorber)
- lastHADlogicAbsorber->SetSensitiveDetector(HADcalorimeterSD);
- */
-
- /*
-
- //
- // Visualization attributes
- //
-
- logicWorld->SetVisAttributes (G4VisAttributes::Invisible);
- G4VisAttributes* simpleBoxVisAtt= new G4VisAttributes(G4Colour(1.0,1.0,1.0));
- simpleBoxVisAtt->SetVisibility(true);
- logicCalor->SetVisAttributes(simpleBoxVisAtt);
-
- */
- //
- //always return the physical World
- //
- return physiWorld;
-}
-
-void MyDetectorConstruction::PrintCalorParameters()
-{
- G4cout << "\n The calorimeter is made of " << NbOfEMLayers << "EM layers and "
- << NbOfHADLayers << "HAD layers of: [ "
- << AbsorberThickness/mm << "mm of " << MedAbsMaterial->GetName()
- << " + "
- << GapThickness/mm << "mm of " << GapMaterial->GetName() << " ] "
- << endl;
-}
-
-void MyDetectorConstruction::SetMagField(G4double fieldValue)
-{
- //apply a global uniform magnetic field along Z axis
- G4FieldManager* fieldMgr
- = G4TransportationManager::GetTransportationManager()->GetFieldManager();
-
- if(magField) delete magField; //delete the existing magn field
-
- if(fieldValue!=0.) // create a new one if non nul
- { magField = new G4UniformMagField(G4ThreeVector(0.,0.,fieldValue));
- fieldMgr->SetDetectorField(magField);
- fieldMgr->CreateChordFinder(magField);
- } else {
- magField = NULL;
- fieldMgr->SetDetectorField(magField);
- }
-}
-
-
-
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-//
-// T36LayerParam.cc, 3/III/99, Sara Vanini
-// parametrisation for t36 layers
-//
-
-#include "T36LayerParam.hh"
-
-#include "G4VPhysicalVolume.hh"
-#include "G4ThreeVector.hh"
-#include "G4Box.hh"
-
-T36LayerParam::T36LayerParam(G4double LayerThickness, G4double NbOfLayers)
-{
- fLayerThickness = LayerThickness;
- fNbOfLayers = NbOfLayers;
-}
-
-T36LayerParam::~T36LayerParam()
-{}
-
-void T36LayerParam::ComputeTransformation
-(const G4int copyNo,G4VPhysicalVolume *physVol) const
-{
- static int counter = 0;
- counter +=1;
- G4double Xposition= - fLayerThickness/2*(fNbOfLayers-1) + copyNo*fLayerThickness;
- G4ThreeVector origin(Xposition,0,0);
-
- physVol->SetTranslation(origin);
-}
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-//
-// T36ModuleParam.cc, 3/III/99, Sara Vanini
-// parametrisation for t36 modules
-//
-
-#include "T36ModuleParam.hh"
-
-#include "G4VPhysicalVolume.hh"
-#include "G4ThreeVector.hh"
-#include "G4Box.hh"
-#include "G4RotationMatrix.hh"
-
-
-T36ModuleParam::T36ModuleParam(G4double ModuleSizeZ, G4double NbOfModules)
-{
- fModuleSizeZ = ModuleSizeZ;
- fNbOfModules = NbOfModules;
-}
-
-T36ModuleParam::~T36ModuleParam()
-{}
-
-void T36ModuleParam::ComputeTransformation
-(const G4int copyNo,G4VPhysicalVolume *physVol) const
-{
- G4double Zposition= - fModuleSizeZ/2*(fNbOfModules-1) + copyNo*fModuleSizeZ;
- G4ThreeVector origin(0,0,Zposition);
- physVol->SetTranslation(origin);
-}
+++ /dev/null
-# --------------------------------------------------------------
-# GNUmakefile for examples module. Sara Vanini, 13/III/00
-# --------------------------------------------------------------
-
-name := mainFLUGG
-G4TARGET := $(name)
-G4EXLIB := true
-
-ifndef FLUGGINSTALL
- FLUGGINSTALL = ../../..
-endif
-
-.PHONY: all
-all: lib bin
-
-include $(FLUGGINSTALL)/config/binmake.gmk
+++ /dev/null
-// Define G4GEOMETRY_DEBUG for debugging information on cout
-
-#include "FGeometryInit.hh"
-#include "MyDetectorConstruction.hh"
-
-#define flukam flukam_
-
-extern "C" void flukam(const G4int & GeoFlag);
-
-int main() {
-
- FGeometryInit* theFGeometryInit = FGeometryInit::GetInstance();
-
- theFGeometryInit
- ->setDetConstruction(new MyDetectorConstruction());
-
-//flag for geometry:
-// 1 for GEANT4
-// 0 for FLUKA
-// 2 for Rubia
- const G4int flag = 1;
-
-//call fortran
- flukam(flag);
-
-//end
- return 0;
-}
-
-
-
-
+++ /dev/null
-# $Id$
-# Flugg tag $Name$
-
-# Id: GNUmakefile,v 2.1 1998/07/02 15:26:52 gunter Exp
-# --------------------------------------------------------------
-# GNUmakefile for examples module. Gabriele Cosmo, 06/04/98
-# --------------------------------------------------------------
-
-name := test
-G4TARGET := $(name)
-G4EXLIB := true
-
-ifndef FLUGGINSTALL
- FLUGGINSTALL = ../../..
-endif
-
-.PHONY: all
-all: lib bin
-
-include $(FLUGGINSTALL)/config/binmake.gmk
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#ifndef MyDetectorConstruction_h
-#define MyDetectorConstruction_h 1
-
-#include "G4VUserDetectorConstruction.hh"
-#include "G4RotationMatrix.hh"
-#include "globals.hh"
-
-class G4VPhysicalVolume;
-class G4Material;
-
-enum Shape {
-
- // CSG solids
- kBox, kTubs,kCons, kTorus, kTrd, kTrap, kPara, kSphere,
-
- // specific solids
- kPolyhedra, kPolycone, kEllipticalTube, kHype
-};
-
-
-class MyDetectorConstruction : public G4VUserDetectorConstruction
-{
- public:
- MyDetectorConstruction(Shape shape, G4int number,G4bool rotate);
- virtual ~MyDetectorConstruction();
-
- // methods
- virtual G4VPhysicalVolume* Construct();
-
- private:
- MyDetectorConstruction() {}
-
- // methods
- void CreateAir();
- void CreateWorld();
-
- // CSG solids
- void CreateBox();
- void CreateTubs();
- void CreateCons();
- void CreateTorus();
- void CreateTrd();
- void CreateTrap();
- void CreatePara();
- void CreateSphere();
-
- // specific solids
- void CreatePolyhedra();
- void CreatePolycone();
- void CreateEllipticalTube();
- void CreateHype();
-
- // data members
- G4Material* fAir;
- G4VPhysicalVolume* fWorld;
- G4RotationMatrix* fRotation;
-
- Shape fShape;
- G4int fNumber;
- G4bool fRotate;
-
- const G4double fWorldX;
- const G4double fWorldY;
- const G4double fWorldZ;
-
- const G4double fX;
- const G4double fY;
- const G4double fZ;
-
- const G4double fR;
- const G4double fDR;
- const G4double fSphi;
- const G4double fDphi;
- const G4double fFullPhi;
-
- const G4double fAlpha;
- const G4double fTheta;
- const G4double fPhi;
-
- const G4int fNofZPlanes;
- const G4int fNofSides;
-};
-
-#endif
-
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#include "MyDetectorConstruction.hh"
-
-#include "G4Material.hh"
-#include "G4MaterialTable.hh"
-#include "G4Element.hh"
-#include "G4ElementTable.hh"
-#include "G4Box.hh"
-#include "G4Tubs.hh"
-#include "G4Cons.hh"
-#include "G4Torus.hh"
-#include "G4Trd.hh"
-#include "G4Trap.hh"
-#include "G4Para.hh"
-#include "G4Sphere.hh"
-#include "G4Polyhedra.hh"
-#include "G4Polycone.hh"
-#include "G4EllipticalTube.hh"
-#include "G4Hype.hh"
-#include "G4LogicalVolume.hh"
-#include "G4LogicalVolumeStore.hh"
-#include "G4RotationMatrix.hh"
-#include "G4ThreeVector.hh"
-#include "G4Transform3D.hh"
-#include "G4PVPlacement.hh"
-
-#include "globals.hh"
-
-//_____________________________________________________________________________
-MyDetectorConstruction::MyDetectorConstruction(Shape shape,
- G4int number,
- G4bool rotate)
- : fAir(0),
- fWorld(0),
-
- fShape(shape),
- fNumber(number),
- fRotate(rotate),
-
- fWorldX(1.*m),
- fWorldY(1.*m),
- fWorldZ(1.*m),
-
- fX(5.*cm),
- fY(10.*cm),
- fZ(15.*cm),
-
- fR (5.*cm),
- fDR(5.*cm),
- fSphi (0.*deg),
- fDphi (315.*deg),
- fFullPhi(360.*deg),
-
- fAlpha(45.*deg),
- fTheta(30.*deg),
- fPhi (15.*deg),
-
- fNofZPlanes(4),
- fNofSides(8)
-{
- fRotation = new G4RotationMatrix();
- fRotation->rotateY(5.*deg);
-}
-
-//_____________________________________________________________________________
-MyDetectorConstruction::~MyDetectorConstruction()
-{;}
-
-//
-// private methods
-//
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreateAir()
-{
- G4double a, iz, z, density;
- G4String name, symbol;
- G4double temperature, pressure;
- G4int nel;
-
- //Air
- a = 14.01*g/mole;
- G4Element* elN = new G4Element(name="Nitrogen", symbol="N", iz=7., a);
- a = 16.00*g/mole;
- G4Element* elO = new G4Element(name="Oxigen", symbol="O", iz=8., a);
- density = 1.29*mg/cm3;
- fAir = new G4Material(name="Air", density, nel=2);
- fAir->AddElement(elN, .7);
- fAir->AddElement(elO, .3);
-}
-
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreateWorld()
-{
-
- G4Box* box
- = new G4Box("world_s", fWorldX, fWorldY, fWorldZ);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(box, fAir, "world_l", 0, 0, 0);
-
- fWorld
- = new G4PVPlacement(0, G4ThreeVector(), "world_p", log, 0, false, 0);
-
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreateBox()
-{
-
- // big box
-
- G4Box* solid_env
- = new G4Box("box_env_s", fX, fY, fZ*fNumber);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "box_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "box_env_p", log_env, fWorld, false, 0);
-
- // small daughter boxes
-
- if (fNumber>1) {
- G4Box* box
- = new G4Box("box_s", fX, fY, fZ);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(box, fAir, "box_l", 0, 0, 0);
-
- for (G4int i=0; i<fNumber; i++) {
-
- // position
- G4double zpos = -fNumber*fZ + fZ + i*2*fZ;
- G4ThreeVector position(0., 0., zpos);
-
- new G4PVPlacement(0, position, log, "box_p", log_env, false, 0);
- }
- }
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreateTubs()
-{
-
- // big tub segment
-
- G4Tubs* solid_env
- = new G4Tubs("tubs_env_s", fR, fR+fDR*fNumber, fZ, fSphi, fDphi);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "tubs_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "tubs_env_p", log_env, fWorld, false, 0);
-
- // small daughter tub segments
-
- if (fNumber>1) {
- for (G4int i=0; i<fNumber; i++) {
- G4Tubs* solid
- = new G4Tubs("tubs_s", fR+i*fDR, fR+(i+1)*fDR, fZ, fSphi, fDphi);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(solid, fAir, "tubs_l", 0, 0, 0);
-
- // position
- G4ThreeVector position(0., 0., 0.);
-
- new G4PVPlacement(0, position, log, "tubs_p", log_env, false, 0);
- }
- }
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreateTorus()
-{
-
- // big tourus segment
-
- G4Torus* solid_env
- = new G4Torus("torus_env_s",
- fR, fR+fDR, fR+5*fDR, fSphi, fDphi);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "torus_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "torus_env_p", log_env, fWorld, false, 0);
-
- // small daughter torus segments
- if (fNumber>1) {
- for (G4int i=0; i<fNumber; i++) {
- G4Torus* solid
- = new G4Torus("torus_s",
- fR, fR+fDR, fR+5*fDR,
- fSphi+i*(fDphi/fNumber), fDphi/fNumber);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(solid, fAir, "torus_l", 0, 0, 0);
-
- // position
- G4ThreeVector position(0., 0., 0.);
-
- new G4PVPlacement(0, position, log, "torus_p", log_env, false, 0);
- }
- }
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreateCons()
-{
-
- // big cons segment
-
- G4Cons* solid_env
- = new G4Cons("cons_env_s",
- fR, fR+fDR*fNumber, fR+fDR, fR+fDR+fDR*fNumber,
- fZ, fSphi, fDphi);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "cons_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "cons_env_p", log_env, fWorld, false, 0);
-
- // small daughter cons segments
-
- if (fNumber>1) {
- for (G4int i=0; i<fNumber; i++) {
- G4Cons* solid
- = new G4Cons("cons_s",
- fR+i*fDR, fR+(i+1)*fDR, fR+fDR+i*fDR, fR+fDR+(i+1)*fDR,
- fZ, fSphi, fDphi);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(solid, fAir, "cons_l", 0, 0, 0);
-
- // position
- G4ThreeVector position(0., 0., 0.);
-
- new G4PVPlacement(0, position, log, "cons_p", log_env, false, 0);
- }
- }
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreateTrd()
-{
-
- // big trd segment
-
- G4Trd* solid_env
- = new G4Trd("trd_env_s", 2*fX, 4*fX, 2*fY, 4*fY, fZ*fNumber);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "trd_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "trd_env_p", log_env, fWorld, false, 0);
-
- // small daughter trd segments
-
- G4Trd* solid
- = new G4Trd("trd_s", fX, 2*fX, fY, 2*fY, fZ);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(solid, fAir, "trd_l", 0, 0, 0);
-
- if (fNumber>1) {
- for (G4int i=0; i<fNumber; i++) {
- // position
- G4double zpos = -fNumber*fZ + fZ + i*2*fZ;
- G4ThreeVector position(0., 0., zpos);
-
- new G4PVPlacement(0, position, log, "trd_p", log_env, false, 0);
- }
- }
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreateTrap()
-{
- // big trap
-
- G4double theta = fTheta;
- G4double phi = fPhi;
- if (fNumber>1) {
- theta = 0.;
- phi = 0.;
- }
-
- G4Trap* solid_env
- = new G4Trap("trap_env_s",
- fZ*fNumber, theta, phi,
- fY, fX, 2*fX, 5.0*deg, fY, fX, 2*fX, 5.0*deg);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "trap_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "trap_env_p", log_env, fWorld, false, 0);
-
- // small daughter trap segments
-
- G4Trap* solid
- = new G4Trap("trap_env_s",
- fZ, theta, phi,
- fY, fX, 2*fX, 5.0*deg, fY, fX, 2*fX, 5.0*deg);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(solid, fAir, "trap_l", 0, 0, 0);
-
- if (fNumber>1) {
- for (G4int i=0; i<fNumber; i++) {
- // position
- G4double zpos = -fNumber*fZ + fZ + i*2*fZ;
- G4ThreeVector position(0., 0., zpos);
-
- new G4PVPlacement(0, position, log, "trap_p", log_env, false, 0);
- }
- }
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreatePara()
-{
-
- G4double theta = fTheta;
- G4double phi = fPhi;
- if (fNumber>1) {
- theta = 0.;
- phi = 0.;
- }
-
- // big para segment
-
- G4Para* solid_env
- = new G4Para("para_env_s", fX, fY, fZ*fNumber, fAlpha, theta, phi);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "para_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "para_env_p", log_env, fWorld, false, 0);
-
- // small daughter para segments
-
- G4Para* solid
- = new G4Para("para_s", fX, fY, fZ, fAlpha, theta, phi);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(solid, fAir, "para_l", 0, 0, 0);
-
- if (fNumber>1) {
- for (G4int i=0; i<fNumber; i++) {
- // position
- G4double zpos = -fNumber*fZ + fZ + i*2*fZ;
- G4ThreeVector position(0., 0., zpos);
-
- new G4PVPlacement(0, position, log, "para_p", log_env, false, 0);
- }
- }
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreateSphere()
-{
-
- // big sphere segment
-
- G4Sphere* solid_env
- = new G4Sphere("sphere_env_s",
- fR, fR+fDR*fNumber, fSphi, fDphi, fSphi, fDphi/2.);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "sphere_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "sphere_env_p", log_env, fWorld, false, 0);
-
- // small daughter sphere segments
- if (fNumber>1) {
- for (G4int i=0; i<fNumber; i++) {
- G4Sphere* solid
- = new G4Sphere("sphere_s",
- fR+i*fDR, fR+(i+1)*fDR, fSphi, fDphi, fSphi, fDphi/2.);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(solid, fAir, "sphere_l", 0, 0, 0);
-
- // position
- G4ThreeVector position(0., 0., 0.);
-
- new G4PVPlacement(0, position, log, "sphere_p", log_env, false, 0);
- }
- }
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreatePolyhedra()
-{
- // big polyhedra
-
- G4double* z_env = new G4double[fNofZPlanes];
- G4double* rmin_env = new G4double[fNofZPlanes];
- G4double* rmax_env = new G4double[fNofZPlanes];
-
- for (G4int i=0; i< fNofZPlanes; i++) {
- z_env[i] = -fNofZPlanes*fZ + fZ + i*2*fZ;
- rmin_env[i] = fR + (i+1) * fDR + (2*(i%2)-1) * (i+1) * fDR/2.;
- rmax_env[i] = fR + (i+2) * fDR + (2*(i%2)-1) * (i+1) * fDR/2.;
- }
-
- G4Polyhedra* solid_env
- = new G4Polyhedra("polyhedra_env_s", fSphi, fFullPhi,
- fNofSides, fNofZPlanes, z_env, rmin_env, rmax_env);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "polyhedra_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "polyhedra_env_p", log_env, fWorld, false, 0);
-
- // small daughter polyhedra segments
-
- if (fNumber>1) {
- for (G4int j=0; j<fNumber; j++) {
-
- G4double* z = new G4double[fNofZPlanes];
- G4double* rmin = new G4double[fNofZPlanes];
- G4double* rmax = new G4double[fNofZPlanes];
-
- for (G4int k=0; k< fNofZPlanes; k++) {
- G4double dr = (rmax_env[k] - rmin_env[k])/fNumber;
- z[k] = z_env[k];
- rmin[k] = rmin_env[k] + j*dr;
- rmax[k] = rmin_env[k] + (j+1)*dr;
- }
-
- G4Polyhedra* solid
- = new G4Polyhedra("polyhedra_s", fSphi, fFullPhi,
- fNofSides, fNofZPlanes, z, rmin, rmax);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(solid, fAir, "polyhedra_l", 0, 0, 0);
-
- // position
- G4ThreeVector position(0., 0., 0.);
- new G4PVPlacement(0, position, log, "polyhedra_p", log_env, false, 0);
-
- delete [] z;
- delete [] rmin;
- delete [] rmax;
- }
- }
-
- delete [] z_env;
- delete [] rmin_env;
- delete [] rmax_env;
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreatePolycone()
-{
- // big polycone
-
- G4double* z_env = new G4double[fNofZPlanes];
- G4double* rmin_env = new G4double[fNofZPlanes];
- G4double* rmax_env = new G4double[fNofZPlanes];
-
- for (G4int i=0; i< fNofZPlanes; i++) {
- z_env[i] = -fNofZPlanes*fZ + fZ + i*2*fZ;
- rmin_env[i] = fR + (i+1) * fDR + (2*(i%2)-1) * (i+1) * fDR/2.;
- rmax_env[i] = fR + (i+2) * fDR + (2*(i%2)-1) * (i+1) * fDR/2.;
- }
-
- G4Polycone* solid_env
- = new G4Polycone("polycone_env_s", fSphi, fFullPhi,
- fNofZPlanes, z_env, rmin_env, rmax_env);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "polycone_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "polycone_env_p", log_env, fWorld, false, 0);
-
- // small daughter polycone segments
-
- if (fNumber>1) {
- for (G4int j=0; j<fNumber; j++) {
-
- G4double* z = new G4double[fNofZPlanes];
- G4double* rmin = new G4double[fNofZPlanes];
- G4double* rmax = new G4double[fNofZPlanes];
-
- for (G4int k=0; k< fNofZPlanes; k++) {
- G4double dr = (rmax_env[k] - rmin_env[k])/fNumber;
- z[k] = z_env[k];
- rmin[k] = rmin_env[k] + j*dr;
- rmax[k] = rmin_env[k] + (j+1)*dr;
- }
-
- G4Polycone* solid
- = new G4Polycone("polycone_s", fSphi, fFullPhi,
- fNofZPlanes, z, rmin, rmax);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(solid, fAir, "polycone_l", 0, 0, 0);
-
- // position
- G4ThreeVector position(0., 0., 0.);
- new G4PVPlacement(0, position, log, "polycone_p", log_env, false, 0);
-
- delete [] z;
- delete [] rmin;
- delete [] rmax;
- }
- }
-
- delete [] z_env;
- delete [] rmin_env;
- delete [] rmax_env;
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreateEllipticalTube()
-{
-
- // big elliptical tube
- // visualization not available
-
- G4EllipticalTube* solid_env
- = new G4EllipticalTube("eltube_env_s",
- fR + 2.*fDR/5., fR + 4.*fDR/5., fZ);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "eltube_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "eltube_env_p", log_env, fWorld, false, 0);
-
- // smaller daughter elliptical tube
- if (fNumber>1) {
- G4EllipticalTube* solid
- = new G4EllipticalTube("eltube_s", fR, fR + fDR/5., fZ);
-
- G4LogicalVolume* log
- = new G4LogicalVolume(solid, fAir, "eltube_l", 0, 0, 0);
-
- // position
- G4ThreeVector position(0., 0., 0.);
-
- new G4PVPlacement(0, position, log, "eltube_p", log_env, false, 0);
- }
-
- if (fNumber>2) {
- // not supported cases
- G4cout << "Elliptical tube:" << G4endl;
- G4cout << "Only one daugher is implemented in the test." << G4endl;
- }
-}
-
-//_____________________________________________________________________________
-void MyDetectorConstruction::CreateHype()
-{
-
- // hyperbolic volume
- // visualization not available
-
- G4Hype* solid_env
- = new G4Hype("hype_env_s",
- fR, fR + fDR, fPhi, fTheta, fZ);
-
- G4LogicalVolume* log_env
- = new G4LogicalVolume(solid_env, fAir, "hype_env_l", 0, 0, 0);
-
- // rotation
- G4RotationMatrix* rotation = 0;
- if (fRotate) rotation = fRotation;
-
- new G4PVPlacement(rotation, G4ThreeVector(0.,0.,0.),
- "hype_env_p", log_env, fWorld, false, 0);
- if (fNumber>1) {
- // not supported cases
- G4cout << "Hype:" << G4endl;
- G4cout << "No daughers are implemented in the test." << G4endl;
- }
-}
-
-//
-// public methods
-//
-
-//_____________________________________________________________________________
-G4VPhysicalVolume* MyDetectorConstruction::Construct()
-{
- CreateAir();
-
- CreateWorld();
-
- switch (fShape) {
-
- // CSG solids
- case kBox:
- CreateBox();
- break;
- case kTubs:
- CreateTubs();
- break;
- case kCons:
- CreateCons();
- break;
- case kTorus:
- CreateTorus();
- break;
- case kTrd:
- CreateTrd();
- break;
- case kTrap:
- CreateTrap();
- break;
- case kPara:
- CreatePara();
- break;
- case kSphere:
- CreateSphere();
- break;
-
- // specific solids
- case kPolyhedra:
- CreatePolyhedra();
- break;
- case kPolycone:
- CreatePolycone();
- break;
- case kEllipticalTube:
- CreateEllipticalTube();
- break;
- case kHype:
- CreateHype();
- break;
- }
-
- return fWorld;
-}
+++ /dev/null
-// $Id$
-// Flugg tag $Name$
-
-#include "FGeometryInit.hh"
-#include "MyDetectorConstruction.hh"
-
-#define flukam flukam_
-
-extern "C" void flukam(const G4int & GeoFlag);
-
-int main() {
-
- FGeometryInit* theFGeometryInit = FGeometryInit::GetInstance();
-
- // Test cases:
- // MyDetectorConstruction(Shape shape, G4int number, G4bool rotate);
- //
- // Shapes:
- // kBox, kTubs,kCons, kTorus, kTrd, kTrap, kPara, kSphere,
- // kPolyhedra, kPolycone, kEllipticalTube, kHype
- //
- // Number: in the interval < 1, 5 >
- // =1: one volume without daughters
- // >1: a volume with embedded daughters
-
- MyDetectorConstruction* detector
- = new MyDetectorConstruction(kPara, 5, true);
-
- theFGeometryInit->setDetConstruction(detector);
-
-//flag for geometry:
-// 1 for GEANT4
-// 0 for FLUKA
-// 2 for Rubia
- const G4int flag = 1;
-
-//call fortran
- flukam(flag);
-
-//end
- return 0;
-}
-
-
-
-
+++ /dev/null
-* $Id$
-* Flugg tag $Name$
-*
-GLOBAL -1.
-TITLE
-ARLW500: 500 MeV p on Cu + 5 m - concrete shield: AF**2+GRS February 1994
-******************* BEAM, EVENT, DISCARD & THRESHOLDS *******************
-*== Set up BEAM parameters ============================================
-*needs arlw.pemf and fluarlw.o
-*
-DEFAULTS NEW-DEFA
-BEAM -0.5 PROTON
-*== BE CAREFUL: the beam z-starting point is supposed to be 1 m before ==
-*== the "effective" center of the target!!!! This is used by the sco- ==
-*== ring routines to compute angles properly. ==
-*DISCARD 3.0 4.0
-BEAMPOS 0.0 0.0 -110.0
-*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
-*
-*== Make sure we get the new EVENT generator ==========================
-EVENTYPE 1.0 EVAP
-*
-* /C.O.group/ NAA group/ NAA fact./ from / to /in step of
-*
-LOW-BIAS 73.0 72.0 0.95 1.0 50.0
-*
-GEOBEGIN COMBINAT
-GEOEND
-*
-* =============== Material to region correspondence ==============
-*
-* Material assignments
-** Mat. # 1 (black hole) in regions 2 and 29 no magfld
-* Mat. # 1 (black hole) in regions 1 to 4 no magfld
-ASSIGNMAT 1.0 2.0 29.0 27.0
-*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
-** Mat. # 12 (copper) in region 5 no magfld
-* Mat. # 12 (copper) in region 3 no magfld
-ASSIGNMAT 12.0 3.0
-** Mat. # 2 (vacuum) in region 6 no magfld
-* Mat. # 2 (vacuum) in region 1 no magfld
-ASSIGNMAT 2.0 1.0
-** Mat. # 25 (concrete) in regions 7-31 no magfld
-* Mat. # 25 (concrete) in regions 4-28 no magfld
-ASSIGNMAT 25.0 4.0 28.0
-*
-*
-******************* MATERIALS *******************************************
-*== New Material Definitions ==========================================
-*
-* Correctly define water, air and concrete and
-* define elements Na,Si,Ar,K,Ca which are not predefined and are needed
-* as components of concrete and air: Na=26,Si=27,Ar=28,K=29 and Ca=27
-*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
-MATERIAL 1.0 1.0079 0.089 3.0 1.0HYDROGEN
-MATERIAL 6.0 12.01 2.25 6.0 CARBON
-MATERIAL 7.0 14.007 0.001251 7.0 NITROGEN
-MATERIAL 8.0 15.999 0.001429 8.0 OXYGEN
-MATERIAL 12.0 24.305 1.738 9.0 MAGNESIU
-MATERIAL 13.0 26.982 2.70 10.0 ALUMINUM
-MATERIAL 26.0 55.847 7.87 11.0 IRON
-MATERIAL 29.0 63.546 8.96 12.0 COPPER
-MATERIAL 0.0 0.0 1.0 20.0 WATER
-MATERIAL 0.0 0.0 0.001205 24.0 AIR
-MATERIAL 0.0 0.0 2.35 25.0 CONCRETE
-MATERIAL 11.0 22.99 0.971 26.0 SODIUM
-MATERIAL 14.0 28.09 2.33 27.0 SILICON
-MATERIAL 18.0 39.99 0.001780 28.0 ARGON
-MATERIAL 19.0 39.10 0.86 29.0 POTASSIU
-MATERIAL 20.0 40.08 1.55 30.0 CALCIUM
-*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
-* Define material 20 as a compound (Water) H and O
-* overriding the predefined material
-COMPOUND 2.0 3.0 1.0 8.0 WATER
-* Define material 24 as a compound (Air at NTP density) N, O and Ar
-* overriding the predefined material
-COMPOUND -0.755 7.0 -0.232 8.0 -0.013 28.0AIR
-* Define material 25 as a compound (Concrete) H,O,Na,Mg,Al,Si,C,K,Ca,Fe
-* overriding the predefined material
-COMPOUND -0.01 3.0 -0.529107 8.0 -0.016 26.0CONCRETE
-COMPOUND -0.002 9.0 -0.033872 10.0 -0.337021 27.0CONCRETE
-COMPOUND -0.001 6.0 -0.013 29.0 -0.044 30.0CONCRETE
-COMPOUND -0.014 11.0 CONCRETE
-*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
-*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
-*
-*== Set up the correspondence between the FLUKA materials and the ======
-*== low-energy neutron module (MORSE) materials ========================
-* / FL92 Ind/ 1st ID / 2nd ID / 3rd ID /
-* LOW-MAT 31.0 293.0 CARBON
-LOW-MAT 28.0 293.0 ARGON
-*
-*== Switch on correct electromagnetic transport =========================
-*
-EMF
-EMFCUT -0.001 0.0001 1.0 50.0
-EMFRAY 1.0 1.0 50.0
-*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
-*
-*************** BIASSING FOR PRODUCTION AND TRANSPORT *******************
-*
-WW-THRES 100.0 100.0 1.0 39.0
-WW-THRES 1.0 1.0 40.0 40.0
-*
-WW-PROFI 1.0 1.0 10.0 1.0 1.0
-WW-PROFI 1.4 11.0 28.0 1.0 1.0
-WW-PROFI 1.96 29.0 40.0 1.0 1.0
-WW-PROFI 2.74 41.0 52.0 1.0 1.0
-WW-PROFI 3.84 53.0 64.0 1.0 1.0
-WW-PROFI 4.61 65.0 71.0 1.0 1.0
-WW-PROFI 5.53 72.0 72.0 1.0 1.0
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-BIASING 0.0 0.0 8.962E-02 1.0 50.0
-BIASING 0.0 0.0 8.962E-02 1.0
-BIASING 0.0 0.0 8.962E-02 4.0
-BIASING 0.0 0.0 1.452E-01 5.0
-BIASING 0.0 0.0 2.352E-01 6.0
-BIASING 0.0 0.0 3.810E-01 7.0
-BIASING 0.0 0.0 6.173E-01 8.0
-BIASING 0.0 0.0 1.000E+00 9.0
-BIASING 0.0 0.0 1.620E+00 10.0
-BIASING 0.0 0.0 2.624E+00 11.0
-BIASING 0.0 0.0 4.252E+00 12.0
-BIASING 0.0 0.0 6.887E+00 13.0
-BIASING 0.0 0.0 1.116E+01 14.0
-BIASING 0.0 0.0 1.808E+01 15.0
-BIASING 0.0 0.0 2.928E+01 16.0
-BIASING 0.0 0.0 4.744E+01 17.0
-BIASING 0.0 0.0 7.685E+01 18.0
-BIASING 0.0 0.0 1.245E+02 19.0
-BIASING 0.0 0.0 2.017E+02 20.0
-BIASING 0.0 0.0 3.267E+02 21.0
-BIASING 0.0 0.0 5.293E+02 22.0
-BIASING 0.0 0.0 8.575E+02 23.0
-BIASING 0.0 0.0 1.389E+03 24.0
-BIASING 0.0 0.0 2.250E+03 25.0
-BIASING 0.0 0.0 3.645E+03 26.0
-BIASING 0.0 0.0 5.906E+03 27.0
-BIASING 0.0 0.0 9.567E+03 28.0
-BIASING 0.0 0.0 9.567E+03 2.0
-*
-WW-FACTO 2.000E-01 8.000E+00 0.0 1.0 50.0 1.0
-WW-FACTO 2.000E-01 8.000E+00 0.0 1.0 1.0
-WW-FACTO 2.000E-01 8.000E+00 0.0 4.0 1.0
-WW-FACTO 2.000E-01 4.938E+00 0.0 5.0 1.0
-WW-FACTO 1.905E-01 3.048E+00 0.0 6.0 1.0
-WW-FACTO 1.176E-01 1.882E+00 0.0 7.0 1.0
-WW-FACTO 7.260E-02 1.162E+00 0.0 8.0 1.0
-WW-FACTO 4.481E-02 7.170E-01 0.0 9.0 1.0
-WW-FACTO 2.766E-02 4.426E-01 0.0 10.0 1.0
-WW-FACTO 1.708E-02 2.732E-01 0.0 11.0 1.0
-WW-FACTO 1.054E-02 1.686E-01 0.0 12.0 1.0
-WW-FACTO 6.506E-03 1.041E-01 0.0 13.0 1.0
-WW-FACTO 4.016E-03 6.426E-02 0.0 14.0 1.0
-WW-FACTO 2.479E-03 3.967E-02 0.0 15.0 1.0
-WW-FACTO 1.530E-03 2.449E-02 0.0 16.0 1.0
-WW-FACTO 9.447E-04 1.511E-02 0.0 17.0 1.0
-WW-FACTO 5.831E-04 9.330E-03 0.0 18.0 1.0
-WW-FACTO 3.600E-04 5.759E-03 0.0 19.0 1.0
-WW-FACTO 2.222E-04 3.555E-03 0.0 20.0 1.0
-WW-FACTO 1.372E-04 2.194E-03 0.0 21.0 1.0
-WW-FACTO 8.466E-05 1.355E-03 0.0 22.0 1.0
-WW-FACTO 5.226E-05 8.362E-04 0.0 23.0 1.0
-WW-FACTO 3.226E-05 5.162E-04 0.0 24.0 1.0
-WW-FACTO 1.991E-05 3.186E-04 0.0 25.0 1.0
-WW-FACTO 1.229E-05 1.967E-04 0.0 26.0 1.0
-WW-FACTO 7.588E-06 1.214E-04 0.0 27.0 1.0
-WW-FACTO 4.684E-06 7.494E-05 0.0 28.0 1.0
-WW-FACTO 4.684E-06 7.494E-05 0.0 2.0 1.0
-*
-* Set multiple scattering on for energies above 30 MeV
-*ACCURACY -0.03 -0.03
-******************* OUTPUT AND SCORING **********************************
-* Minimum output for fluxes, maximum output for densities
-*
-OUTLEVEL 1.0 7.0 1.0
-*
-* Score all-stars
-*
-SCORE 201.0 208.0
-*
-* Boundary crossings
-*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
-EXTRAWEI 1.0
-*
-USRBDX +99.0 1.0 -40.0 4.0 5.0 1.0 P1 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 4.0 5.0 1.0 P1 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 4.0 5.0 1.0 P1 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 4.0 5.0 1.0 P1 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 4.0 5.0 1.0 P1 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 4.0 5.0 1.0 C1 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 4.0 5.0 1.0 C1 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 4.0 5.0 1.0 C1 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 4.0 5.0 1.0 C1 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 4.0 5.0 1.0 C1 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 4.0 5.0 1.0 G1 0-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 4.0 5.0 1.0 G1 5-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 4.0 5.0 1.0 G1 10-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 4.0 5.0 1.0 G1 15-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 4.0 5.0 1.0 G1 20-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 4.0 5.0 1.0 N1 0-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 4.0 5.0 1.0 N1 5-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 4.0 5.0 1.0 N1 10-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 4.0 5.0 1.0 N1 15-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 4.0 5.0 1.0 N1 20-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 5.0 6.0 1.0 P2 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 5.0 6.0 1.0 P2 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 5.0 6.0 1.0 P2 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 5.0 6.0 1.0 P2 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 5.0 6.0 1.0 P2 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 5.0 6.0 1.0 C2 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 5.0 6.0 1.0 C2 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 5.0 6.0 1.0 C2 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 5.0 6.0 1.0 C2 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 5.0 6.0 1.0 C2 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 5.0 6.0 1.0 G2 0-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 5.0 6.0 1.0 G2 5-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 5.0 6.0 1.0 G2 10-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 5.0 6.0 1.0 G2 15-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 5.0 6.0 1.0 G2 20-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 5.0 6.0 1.0 N2 0-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 5.0 6.0 1.0 N2 5-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 5.0 6.0 1.0 N2 10-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 5.0 6.0 1.0 N2 15-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 5.0 6.0 1.0 N2 20-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 6.0 7.0 1.0 P3 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 6.0 7.0 1.0 P3 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 6.0 7.0 1.0 P3 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 6.0 7.0 1.0 P3 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 6.0 7.0 1.0 P3 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 6.0 7.0 1.0 C3 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 6.0 7.0 1.0 C3 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 6.0 7.0 1.0 C3 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 6.0 7.0 1.0 C3 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 6.0 7.0 1.0 C3 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 6.0 7.0 1.0 G3 0-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 6.0 7.0 1.0 G3 5-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 6.0 7.0 1.0 G3 10-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 6.0 7.0 1.0 G3 15-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 6.0 7.0 1.0 G3 20-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 6.0 7.0 1.0 N3 0-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 6.0 7.0 1.0 N3 5-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 6.0 7.0 1.0 N3 10-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 6.0 7.0 1.0 N3 15-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 6.0 7.0 1.0 N3 20-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 7.0 8.0 1.0 P4 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 7.0 8.0 1.0 P4 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 7.0 8.0 1.0 P4 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 7.0 8.0 1.0 P4 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 7.0 8.0 1.0 P4 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 7.0 8.0 1.0 C4 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 7.0 8.0 1.0 C4 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 7.0 8.0 1.0 C4 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 7.0 8.0 1.0 C4 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 7.0 8.0 1.0 C4 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 7.0 8.0 1.0 G4 0-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 7.0 8.0 1.0 G4 5-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 7.0 8.0 1.0 G4 10-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 7.0 8.0 1.0 G4 15-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 7.0 8.0 1.0 G4 20-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 7.0 8.0 1.0 N4 0-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 7.0 8.0 1.0 N4 5-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 7.0 8.0 1.0 N4 10-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 7.0 8.0 1.0 N4 15-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 7.0 8.0 1.0 N4 20-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 8.0 9.0 1.0 P5 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 8.0 9.0 1.0 P5 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 8.0 9.0 1.0 P5 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 8.0 9.0 1.0 P5 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 8.0 9.0 1.0 P5 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 8.0 9.0 1.0 C5 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 8.0 9.0 1.0 C5 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 8.0 9.0 1.0 C5 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 8.0 9.0 1.0 C5 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 8.0 9.0 1.0 C5 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 8.0 9.0 1.0 G5 0-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 8.0 9.0 1.0 G5 5-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 8.0 9.0 1.0 G5 10-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 8.0 9.0 1.0 G5 15-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 8.0 9.0 1.0 G5 20-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 8.0 9.0 1.0 N5 0-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 8.0 9.0 1.0 N5 5-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 8.0 9.0 1.0 N5 10-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 8.0 9.0 1.0 N5 15-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 8.0 9.0 1.0 N5 20-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 9.0 10.0 1.0 P6 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 9.0 10.0 1.0 P6 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 9.0 10.0 1.0 P6 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 9.0 10.0 1.0 P6 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 9.0 10.0 1.0 P6 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 9.0 10.0 1.0 C6 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 9.0 10.0 1.0 C6 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 9.0 10.0 1.0 C6 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 9.0 10.0 1.0 C6 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 9.0 10.0 1.0 C6 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 9.0 10.0 1.0 G6 0-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 9.0 10.0 1.0 G6 5-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 9.0 10.0 1.0 G6 10-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 9.0 10.0 1.0 G6 15-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 9.0 10.0 1.0 G6 20-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 9.0 10.0 1.0 N6 0-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 9.0 10.0 1.0 N6 5-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 9.0 10.0 1.0 N6 10-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 9.0 10.0 1.0 N6 15-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 9.0 10.0 1.0 N6 20-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 10.0 11.0 1.0 P7 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 10.0 11.0 1.0 P7 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 10.0 11.0 1.0 P7 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 10.0 11.0 1.0 P7 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 10.0 11.0 1.0 P7 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 10.0 11.0 1.0 C7 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 10.0 11.0 1.0 C7 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 10.0 11.0 1.0 C7 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 10.0 11.0 1.0 C7 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 10.0 11.0 1.0 C7 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 10.0 11.0 1.0 G7 0-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 10.0 11.0 1.0 G7 5-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 10.0 11.0 1.0 G7 10-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 10.0 11.0 1.0 G7 15-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 10.0 11.0 1.0 G7 20-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 10.0 11.0 1.0 N7 0-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 10.0 11.0 1.0 N7 5-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 10.0 11.0 1.0 N7 10-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 10.0 11.0 1.0 N7 15-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 10.0 11.0 1.0 N7 20-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 11.0 12.0 1.0 P8 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 11.0 12.0 1.0 P8 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 11.0 12.0 1.0 P8 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 11.0 12.0 1.0 P8 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 11.0 12.0 1.0 P8 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 11.0 12.0 1.0 C8 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 11.0 12.0 1.0 C8 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 11.0 12.0 1.0 C8 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 11.0 12.0 1.0 C8 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 11.0 12.0 1.0 C8 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 11.0 12.0 1.0 G8 0-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 11.0 12.0 1.0 G8 5-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 11.0 12.0 1.0 G8 10-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 11.0 12.0 1.0 G8 15-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 11.0 12.0 1.0 G8 20-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 11.0 12.0 1.0 N8 0-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 11.0 12.0 1.0 N8 5-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 11.0 12.0 1.0 N8 10-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 11.0 12.0 1.0 N8 15-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 11.0 12.0 1.0 N8 20-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 12.0 13.0 1.0 P9 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 12.0 13.0 1.0 P9 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 12.0 13.0 1.0 P9 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 12.0 13.0 1.0 P9 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 12.0 13.0 1.0 P9 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 12.0 13.0 1.0 C9 0-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 12.0 13.0 1.0 C9 5-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 12.0 13.0 1.0 C9 10-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 12.0 13.0 1.0 C9 15-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 12.0 13.0 1.0 C9 20-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 12.0 13.0 1.0 G9 0-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 12.0 13.0 1.0 G9 5-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 12.0 13.0 1.0 G9 10-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 12.0 13.0 1.0 G9 15-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 12.0 13.0 1.0 G9 20-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 12.0 13.0 1.0 N9 0-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 12.0 13.0 1.0 N9 5-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 12.0 13.0 1.0 N9 10-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 12.0 13.0 1.0 N9 15-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 12.0 13.0 1.0 N9 20-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 13.0 14.0 1.0 P100-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 13.0 14.0 1.0 P105-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 13.0 14.0 1.0 P1010-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 13.0 14.0 1.0 P1015-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 13.0 14.0 1.0 P1020-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 13.0 14.0 1.0 C100-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 13.0 14.0 1.0 C105-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 13.0 14.0 1.0 C1010-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 13.0 14.0 1.0 C1015-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 13.0 14.0 1.0 C1020-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 13.0 14.0 1.0 G100-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 13.0 14.0 1.0 G105-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 13.0 14.0 1.0 G1010-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 13.0 14.0 1.0 G1015-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 13.0 14.0 1.0 G1020-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 13.0 14.0 1.0 N100-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 13.0 14.0 1.0 N105-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 13.0 14.0 1.0 N1010-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 13.0 14.0 1.0 N1015-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 13.0 14.0 1.0 N1020-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 14.0 15.0 1.0 P110-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 14.0 15.0 1.0 P115-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 14.0 15.0 1.0 P1110-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 14.0 15.0 1.0 P1115-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 14.0 15.0 1.0 P1120-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 14.0 15.0 1.0 C110-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 14.0 15.0 1.0 C115-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 14.0 15.0 1.0 C1110-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 14.0 15.0 1.0 C1115-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 14.0 15.0 1.0 C1120-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 14.0 15.0 1.0 G110-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 14.0 15.0 1.0 G115-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 14.0 15.0 1.0 G1110-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 14.0 15.0 1.0 G1115-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 14.0 15.0 1.0 G1120-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 14.0 15.0 1.0 N110-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 14.0 15.0 1.0 N115-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 14.0 15.0 1.0 N1110-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 14.0 15.0 1.0 N1115-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 14.0 15.0 1.0 N1120-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 15.0 16.0 1.0 P120-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 15.0 16.0 1.0 P125-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 15.0 16.0 1.0 P1210-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 15.0 16.0 1.0 P1215-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 15.0 16.0 1.0 P1220-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 15.0 16.0 1.0 C120-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 15.0 16.0 1.0 C125-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 15.0 16.0 1.0 C1210-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 15.0 16.0 1.0 C1215-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 15.0 16.0 1.0 C1220-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 15.0 16.0 1.0 G120-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 15.0 16.0 1.0 G125-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 15.0 16.0 1.0 G1210-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 15.0 16.0 1.0 G1215-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 15.0 16.0 1.0 G1220-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 15.0 16.0 1.0 N120-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 15.0 16.0 1.0 N125-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 15.0 16.0 1.0 N1210-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 15.0 16.0 1.0 N1215-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 15.0 16.0 1.0 N1220-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 16.0 17.0 1.0 P130-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 16.0 17.0 1.0 P135-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 16.0 17.0 1.0 P1310-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 16.0 17.0 1.0 P1315-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 16.0 17.0 1.0 P1320-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 16.0 17.0 1.0 C130-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 16.0 17.0 1.0 C135-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 16.0 17.0 1.0 C1310-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 16.0 17.0 1.0 C1315-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 16.0 17.0 1.0 C1320-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 16.0 17.0 1.0 G130-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 16.0 17.0 1.0 G135-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 16.0 17.0 1.0 G1310-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 16.0 17.0 1.0 G1315-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 16.0 17.0 1.0 G1320-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 16.0 17.0 1.0 N130-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 16.0 17.0 1.0 N135-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 16.0 17.0 1.0 N1310-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 16.0 17.0 1.0 N1315-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 16.0 17.0 1.0 N1320-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 17.0 18.0 1.0 P140-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 17.0 18.0 1.0 P145-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 17.0 18.0 1.0 P1410-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 17.0 18.0 1.0 P1415-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 17.0 18.0 1.0 P1420-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 17.0 18.0 1.0 C140-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 17.0 18.0 1.0 C145-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 17.0 18.0 1.0 C1410-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 17.0 18.0 1.0 C1415-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 17.0 18.0 1.0 C1420-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 17.0 18.0 1.0 G140-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 17.0 18.0 1.0 G145-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 17.0 18.0 1.0 G1410-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 17.0 18.0 1.0 G1415-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 17.0 18.0 1.0 G1420-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 17.0 18.0 1.0 N140-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 17.0 18.0 1.0 N145-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 17.0 18.0 1.0 N1410-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 17.0 18.0 1.0 N1415-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 17.0 18.0 1.0 N1420-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 18.0 19.0 1.0 P150-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 18.0 19.0 1.0 P155-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 18.0 19.0 1.0 P1510-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 18.0 19.0 1.0 P1515-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 18.0 19.0 1.0 P1520-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 18.0 19.0 1.0 C150-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 18.0 19.0 1.0 C155-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 18.0 19.0 1.0 C1510-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 18.0 19.0 1.0 C1515-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 18.0 19.0 1.0 C1520-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 18.0 19.0 1.0 G150-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 18.0 19.0 1.0 G155-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 18.0 19.0 1.0 G1510-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 18.0 19.0 1.0 G1515-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 18.0 19.0 1.0 G1520-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 18.0 19.0 1.0 N150-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 18.0 19.0 1.0 N155-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 18.0 19.0 1.0 N1510-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 18.0 19.0 1.0 N1515-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 18.0 19.0 1.0 N1520-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 19.0 20.0 1.0 P160-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 19.0 20.0 1.0 P165-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 19.0 20.0 1.0 P1610-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 19.0 20.0 1.0 P1615-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 19.0 20.0 1.0 P1620-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 19.0 20.0 1.0 C160-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 19.0 20.0 1.0 C165-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 19.0 20.0 1.0 C1610-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 19.0 20.0 1.0 C1615-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 19.0 20.0 1.0 C1620-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 19.0 20.0 1.0 G160-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 19.0 20.0 1.0 G165-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 19.0 20.0 1.0 G1610-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 19.0 20.0 1.0 G1615-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 19.0 20.0 1.0 G1620-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 19.0 20.0 1.0 N160-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 19.0 20.0 1.0 N165-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 19.0 20.0 1.0 N1610-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 19.0 20.0 1.0 N1615-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 19.0 20.0 1.0 N1620-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 23.0 24.0 1.0 P170-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 20.0 21.0 1.0 P175-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 20.0 21.0 1.0 P1710-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 20.0 21.0 1.0 P1715-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 20.0 21.0 1.0 P1720-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 20.0 21.0 1.0 C170-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 20.0 21.0 1.0 C175-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 20.0 21.0 1.0 C1710-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 20.0 21.0 1.0 C1715-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 20.0 21.0 1.0 C1720-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 20.0 21.0 1.0 G170-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 20.0 21.0 1.0 G175-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 20.0 21.0 1.0 G1710-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 20.0 21.0 1.0 G1715-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 20.0 21.0 1.0 G1720-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 20.0 21.0 1.0 N170-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 20.0 21.0 1.0 N175-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 20.0 21.0 1.0 N1710-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 20.0 21.0 1.0 N1715-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 20.0 21.0 1.0 N1720-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 21.0 22.0 1.0 P180-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 21.0 22.0 1.0 P185-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 21.0 22.0 1.0 P1810-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 21.0 22.0 1.0 P1815-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 21.0 22.0 1.0 P1820-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 21.0 22.0 1.0 C180-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 21.0 22.0 1.0 C185-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 21.0 22.0 1.0 C1810-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 21.0 22.0 1.0 C1815-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 21.0 22.0 1.0 C1820-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 21.0 22.0 1.0 G180-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 21.0 22.0 1.0 G185-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 21.0 22.0 1.0 G1810-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 21.0 22.0 1.0 G1815-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 21.0 22.0 1.0 G1820-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 21.0 22.0 1.0 N180-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 21.0 22.0 1.0 N185-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 21.0 22.0 1.0 N1810-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 21.0 22.0 1.0 N1815-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 21.0 22.0 1.0 N1820-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 22.0 23.0 1.0 P190-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 22.0 23.0 1.0 P195-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 22.0 23.0 1.0 P1910-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 22.0 23.0 1.0 P1915-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 22.0 23.0 1.0 P1920-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 22.0 23.0 1.0 C190-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 22.0 23.0 1.0 C195-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 22.0 23.0 1.0 C1910-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 22.0 23.0 1.0 C1915-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 22.0 23.0 1.0 C1920-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 22.0 23.0 1.0 G190-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 22.0 23.0 1.0 G195-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 22.0 23.0 1.0 G1910-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 22.0 23.0 1.0 G1915-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 22.0 23.0 1.0 G1920-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 22.0 23.0 1.0 N190-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 22.0 23.0 1.0 N195-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 22.0 23.0 1.0 N1910-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 22.0 23.0 1.0 N1915-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 22.0 23.0 1.0 N1920-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 23.0 24.0 1.0 P200-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 23.0 24.0 1.0 P205-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 23.0 24.0 1.0 P2010-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 23.0 24.0 1.0 P2015-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 23.0 24.0 1.0 P2020-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 23.0 24.0 1.0 C200-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 23.0 24.0 1.0 C205-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 23.0 24.0 1.0 C2010-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 23.0 24.0 1.0 C2015-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 23.0 24.0 1.0 C2020-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 23.0 24.0 1.0 G200-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 23.0 24.0 1.0 G205-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 23.0 24.0 1.0 G2010-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 23.0 24.0 1.0 G2015-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 23.0 24.0 1.0 G2020-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 23.0 24.0 1.0 N200-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 23.0 24.0 1.0 N205-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 23.0 24.0 1.0 N2010-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 23.0 24.0 1.0 N2015-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 23.0 24.0 1.0 N2020-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 24.0 25.0 1.0 P210-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 24.0 25.0 1.0 P215-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 24.0 25.0 1.0 P2110-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 24.0 25.0 1.0 P2115-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 24.0 25.0 1.0 P2120-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 24.0 25.0 1.0 C210-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 24.0 25.0 1.0 C215-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 24.0 25.0 1.0 C2110-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 24.0 25.0 1.0 C2115-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 24.0 25.0 1.0 C2120-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 24.0 25.0 1.0 G210-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 24.0 25.0 1.0 G215-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 24.0 25.0 1.0 G2110-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 24.0 25.0 1.0 G2115-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 24.0 25.0 1.0 G2120-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 24.0 25.0 1.0 N210-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 24.0 25.0 1.0 N215-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 24.0 25.0 1.0 N2110-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 24.0 25.0 1.0 N2115-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 24.0 25.0 1.0 N2120-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 25.0 26.0 1.0 P220-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 25.0 26.0 1.0 P225-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 25.0 26.0 1.0 P2210-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 25.0 26.0 1.0 P2215-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 25.0 26.0 1.0 P2220-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 25.0 26.0 1.0 C220-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 25.0 26.0 1.0 C225-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 25.0 26.0 1.0 C2210-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 25.0 26.0 1.0 C2215-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 25.0 26.0 1.0 C2220-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 25.0 26.0 1.0 G220-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 25.0 26.0 1.0 G225-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 25.0 26.0 1.0 G2210-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 25.0 26.0 1.0 G2215-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 25.0 26.0 1.0 G2220-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 25.0 26.0 1.0 N220-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 25.0 26.0 1.0 N225-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 25.0 26.0 1.0 N2210-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 25.0 26.0 1.0 N2215-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 25.0 26.0 1.0 N2220-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 26.0 27.0 1.0 P230-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 26.0 27.0 1.0 P235-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 26.0 27.0 1.0 P2310-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 26.0 27.0 1.0 P2315-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 26.0 27.0 1.0 P2320-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 26.0 27.0 1.0 C230-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 26.0 27.0 1.0 C235-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 26.0 27.0 1.0 C2310-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 26.0 27.0 1.0 C2315-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 26.0 27.0 1.0 C2320-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 26.0 27.0 1.0 G230-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 26.0 27.0 1.0 G235-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 26.0 27.0 1.0 G2310-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 26.0 27.0 1.0 G2315-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 26.0 27.0 1.0 G2320-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 26.0 27.0 1.0 N230-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 26.0 27.0 1.0 N235-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 26.0 27.0 1.0 N2310-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 26.0 27.0 1.0 N2315-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 26.0 27.0 1.0 N2320-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 27.0 28.0 1.0 P240-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 27.0 28.0 1.0 P245-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 27.0 28.0 1.0 P2410-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 27.0 28.0 1.0 P2415-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 27.0 28.0 1.0 P2420-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 27.0 28.0 1.0 C240-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 27.0 28.0 1.0 C245-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 27.0 28.0 1.0 C2410-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 27.0 28.0 1.0 C2415-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 27.0 28.0 1.0 C2420-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 27.0 28.0 1.0 G240-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 27.0 28.0 1.0 G245-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 27.0 28.0 1.0 G2410-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 27.0 28.0 1.0 G2415-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 27.0 28.0 1.0 G2420-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 27.0 28.0 1.0 N240-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 27.0 28.0 1.0 N245-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 27.0 28.0 1.0 N2410-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 27.0 28.0 1.0 N2415-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 27.0 28.0 1.0 N2420-30dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 1.0 -40.0 28.0 2.0 1.0 P250-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 28.0 2.0 1.0 P255-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 28.0 2.0 1.0 P2510-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 28.0 2.0 1.0 P2515-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 1.0 -40.0 28.0 2.0 1.0 P2520-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 28.0 2.0 1.0 C250-5deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 28.0 2.0 1.0 C255-10deg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 28.0 2.0 1.0 C2510-15dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 28.0 2.0 1.0 C2515-20dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 209.0 -41.0 28.0 2.0 1.0 C2520-30dg
-USRBDX 1.0 0.01 20.0 &
-USRBDX +99.0 7.0 -42.0 28.0 2.0 1.0 G250-5deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 28.0 2.0 1.0 G255-10deg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 28.0 2.0 1.0 G2510-15dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 28.0 2.0 1.0 G2515-20dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 7.0 -42.0 28.0 2.0 1.0 G2520-30dg
-USRBDX 1.0 0.0001 40.0 &
-USRBDX +99.0 8.0 -43.0 28.0 2.0 1.0 N250-5deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 28.0 2.0 1.0 N255-10deg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 28.0 2.0 1.0 N2510-15dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 28.0 2.0 1.0 N2515-20dg
-USRBDX 1.0 1.E-14140.0 &
-USRBDX +99.0 8.0 -43.0 28.0 2.0 1.0 N2520-30dg
-USRBDX 1.0 1.E-14140.0 &
-*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
-*************** START SIMULATION ************************************
-* Read seeds from unit 1
-DUMPTHEM 101 MGDRAW
-RANDOMIZE 1.0
-*== Start the calculation ==============================================
-*
-* /No. protons/No. stars/Print time/Inc.beam//default interactive tlim./
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-START 1.0 99999999.
-STOP
+++ /dev/null
-* $Id$
-* Flugg tag $Name$
-*
-GLOBAL 0.0 -1.0
-TITLE
- 1 GeV p for testing magnetic field
-*needs t36.pemf
-DEFAULTS CALORIME
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-BEAM 1. 1. 1. PROTON
-BEAMPOS 2.0 -4. -1. 0.0 0.0
-* Switch on the new generator/heavies transported /option 13 selected
-EVENTYPE EVAP
-* Switch off all low energy neutron biasing
-LOW-BIAS 73.0 1.0 82.0
-* Switch on detailed multiple scattering for charged hadrons down to threshold
-ACCURACY 1.0 1.0
-* Pair production and bremss. down to 100 keV
-PAIRBREM +3.0 0.0 0.0001 3.0 50.0
-* Geometry
-GEOBEGIN COMBINAT
-GEOEND
-* *** Time cut-off ***
-TIME-CUT 2000.0 1.0 82.0
-* *** Materials ***
-MATERIAL 1.0 1.0079 0.089 3.0 1.0 HYDROGEN
-MATERIAL 6.0 12.01 2.26 6.0 CARBON
-MATERIAL 7.0 14.007 0.001251 7.0 NITROGEN
-MATERIAL 8.0 15.999 0.001429 8.0 OXYGEN
-MATERIAL 18.0 39.948 1.78 9.0 ARGON
-MATERIAL 13.0 29.982 2.70 10.0 ALUMINUM
-MATERIAL 26.0 55.85 7.87 17.0 IRON
-MATERIAL 26.0 55.85 7.87 18.0 IRON0
-MATERIAL 82.0 207.19 11.35 19.0 LEAD
-MATERIAL 51.0 121.75 6.62 20.0 ANTIMONY
-MATERIAL 11.3 21.0 LEADSB0
-COMPOUND -0.96 19.0 -0.04 20.0 LEADSB0
-* Be sure to be consistent with PEMF for e/pi purposes:
-MAT-PROP 804.96 21.0
-STERNHEI 6.159511 0.508 3.0000 0.246848 3.0000 0.0 21.0
-MATERIAL 11.3 22.0 LEADSB
-COMPOUND -0.96 19.0 -0.04 20.0 LEADSB
-* Be sure to be consistent with PEMF for e/pi purposes:
-MAT-PROP 804.96 22.0
-STERNHEI 6.159511 0.508 3.0000 0.246848 3.0000 0.0 22.0
-MATERIAL 1.044 23.0 SCSN38
-* Be sure to be consistent with PEMF for e/pi purposes:
-MAT-PROP 64.7 23.0
-STERNHEI 3.1997 0.1464 2.4855 0.16101 3.2393 0.0 23.0
-* Set the transport accuracy high in Scintillator:
-FLUKAFIX 0.05 23.0
-*
-COMPOUND +1.0 3.0 +1.0 6.0 SCSN38
-MATERIAL 0.001225 24.0 AIR
-COMPOUND -.9256E-03 7.0 -.2837E-03 8.0 -.01572E-3 9.0 AIR
-MATERIAL 1.17 25.0 PMMAWLS
-COMPOUND +1.0 3.0 +0.57 6.0 +0.57 8.0 PMMAWLS
-* Set the correspondence for low energy neutrons *
-LOW-MAT 3.0 1.0 -2.0 293.0 HYDROGEN
-LOW-MAT 6.0 6.0 -2.0 293.0 CARBON
-LOW-MAT 7.0 7.0 14.0 293.0 NITROGEN
-LOW-MAT 8.0 8.0 16.0 293.0 OXYGEN
-LOW-MAT 9.0 18.0 -2.0 293.0 ARGON
-LOW-MAT 10.0 13.0 -2.0 293.0 ALUMINUM
-LOW-MAT 17.0 26.0 -2.0 293.0 IRON
-LOW-MAT 18.0 26.0 -2.0 293.0 IRON
-LOW-MAT 19.0 82.0 -2.0 293.0 LEAD
-LOW-MAT 20.0 51.0 -2.0 293.0 ANTIMONY
-* Set the transport accuracy high in Scintillator:
-FLUKAFIX 0.05 23.0
-*
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-* Define region materials, cut-off's and step sizes *
-* void around
-ASSIGNMAT 2.0 1.0 2.0
-* layers with lead
-ASSIGNMAT 22.0 3.0 8.0 1.0 1.0
-* spheres with scintillator
-ASSIGNMAT 23.0 23.0 81.0 1.0 1.0
-* tubs with alluminum
-ASSIGNMAT 10.0 9.0 22.0 1.0 1.0
-* black hole
-ASSIGNMAT 1.0 82.0
-EMFCUT 0.0001 0.00001 1. 81.
-*
-* *** production
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-* Hadron and muon thresholds: *
-PART-THR -0.003 1.0 39.0 1.0
-PART-THR -0.0196 8.0 8.0 1.0
-PART-THR -0.05 9.0 9.0 1.0
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-* raggio di curvatura: R[m] = p[Gev/c] / (0.3*B[Tesla])
-MGNFIELD 20. .01 .05 60.0 0.0 0.0
-*
-SCORE 208.0 211.0 210.0
-*DISCARD 3.0
-*
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-* output on un-formatted file for cartesian binning
-EVENTBIN -0.1 208. -25. 5.1 4.5 8.0 CartesBin
-EVENTBIN 0. -5.1 0. 25.5 48. 40. &
-* output on un-formatted file for region by region binning
-EVENTBIN 8.0 208. -26. 81. RegionBin
-EVENTBIN 3. 1. &
-*
-OUTLEVEL 1.0 7.0 1.0
-*OUTLEVEL 1.0 7.0
-DUMPTHEM 101. MGDRAW
-*
-RANDOMIZ 1.0
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-START 1.0 99999999. 4000.0 0.0 1.0
-STOP
-
-
-
-
-
-
-
+++ /dev/null
-* $Id$
-* Flugg tag $Name$
-*
-GLOBAL -1.
-TITLE
- 10 GeV pi- on simple geometry for testing bias histories
-*needs argon.pemf
-DEFAULTS CALORIME
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-BEAM 10.0 PION-
-BEAMPOS -1.4 0.5 0.0 1.0
-* Switch on the new generator/heavies transported /option 13 selected
-EVENTYPE EVAP
-* Switch off all low energy neutron biasing
-LOW-BIAS 73.0 1.0 6.0
-* Switch on detailed multiple scattering for charged hadrons down to threshold
-ACCURACY 1.0 1.0
-* Pair production and bremss. down to 100 keV
-PAIRBREM +3.0 0.0 0.0001 3.0 50.0
-* Geometry
-GEOBEGIN COMBINAT
-GEOEND
-* *** Time cut-off ***
-TIME-CUT 2000.0 1.0 6.0
-* *** Materials ***
-MATERIAL 1.0 1.0079 0.089 3.0 1.0 HYDROGEN
-MATERIAL 6.0 12.01 2.26 6.0 CARBON
-MATERIAL 7.0 14.007 0.001251 7.0 NITROGEN
-MATERIAL 8.0 15.999 0.001429 8.0 OXYGEN
-MATERIAL 18.0 39.948 1.78 9.0 ARGON
-MATERIAL 13.0 29.982 2.70 10.0 ALUMINUM
-MATERIAL 26.0 55.85 7.87 17.0 IRON
-MATERIAL 26.0 55.85 7.87 18.0 IRON0
-MATERIAL 82.0 207.19 11.35 19.0 LEAD
-MATERIAL 51.0 121.75 6.62 20.0 ANTIMONY
-MATERIAL 11.3 21.0 LEADSB0
-COMPOUND -0.96 19.0 -0.04 20.0 LEADSB0
-* Be sure to be consistent with PEMF for e/pi purposes:
-MAT-PROP 804.96 21.0
-STERNHEI 6.159511 0.508 3.0000 0.246848 3.0000 0.0 21.0
-MATERIAL 11.3 22.0 LEADSB
-COMPOUND -0.96 19.0 -0.04 20.0 LEADSB
-* Be sure to be consistent with PEMF for e/pi purposes:
-MAT-PROP 804.96 22.0
-STERNHEI 6.159511 0.508 3.0000 0.246848 3.0000 0.0 22.0
-MATERIAL 1.044 23.0 SCSN38
-* Be sure to be consistent with PEMF for e/pi purposes:
-MAT-PROP 64.7 23.0
-STERNHEI 3.1997 0.1464 2.4855 0.16101 3.2393 0.0 23.0
-* Set the transport accuracy high in Scintillator:
-FLUKAFIX 0.05 23.0
-*
-COMPOUND +1.0 3.0 +1.0 6.0 SCSN38
-MATERIAL 0.001225 24.0 AIR
-COMPOUND -.9256E-03 7.0 -.2837E-03 8.0 -.01572E-3 9.0 AIR
-MATERIAL 1.17 25.0 PMMAWLS
-COMPOUND +1.0 3.0 +0.57 6.0 +0.57 8.0 PMMAWLS
-* Set the correspondence for low energy neutrons *
-LOW-MAT 3.0 1.0 -2.0 293.0 HYDROGEN
-LOW-MAT 6.0 6.0 -2.0 293.0 CARBON
-LOW-MAT 7.0 7.0 14.0 293.0 NITROGEN
-LOW-MAT 8.0 8.0 16.0 293.0 OXYGEN
-LOW-MAT 9.0 18.0 -2.0 293.0 ARGON
-LOW-MAT 10.0 13.0 -2.0 293.0 ALUMINUM
-LOW-MAT 17.0 26.0 -2.0 293.0 IRON
-LOW-MAT 18.0 26.0 -2.0 293.0 IRON
-LOW-MAT 19.0 82.0 -2.0 293.0 LEAD
-LOW-MAT 20.0 51.0 -2.0 293.0 ANTIMONY
-* Set the transport accuracy high in Scintillator:
-FLUKAFIX 0.05 23.0
-*
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-* Define region materials, cut-off's and step sizes *
-ASSIGNMAT 24.0 2.0 3.0 0.0 0.0 0.0
-ASSIGNMAT 22.0 4.0 0.0 0.0 0.0 0.0
-ASSIGNMAT 23.0 5.0 0.0 0.0 0.0 0.0
-ASSIGNMAT 1.0 6.0 0.0 0.0 0.0 0.0
-ASSIGNMAT 1.0 1.0 0.0 0.0 0.0 0.0
-STEPSIZE 10.0 1.0 6.0
-STEPSIZE 1.0 4.0 5.0
-STEPSIZE 1.0 6.0
-EMFCUT -0.100E-03+0.010E-03 6.0
-DELTARAY 0.00001 1.05 6.0
-*
-MULSOPT -3.0 -2.0 1.0 50.0
-* *** Switch on accurate treatment of photoelectric interactions and X-ray
-* *** production
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-EMFFLUO 1.0 1.0 6.0
-* *** Switch on Rayleigh (coherent) scattering for all regions
-EMFRAY 1.0 1.0 6.0
-* Hadron and muon thresholds: *
-PART-THR 0.075 1.0 5.0 1.0
-*PART-THR -0.05 9.0 9.0
-*PART-THR -0.0196 8.0 8.0 1.0
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-EMFFIX 3.0 0.1 5.0 0.1
-EMFFIX 4.0 0.1
-*** bias for testing lattice histories ***
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-BIASING 0.0 0.0 8.962E-02 1.0 6.0
-BIASING 0.0 0.0 8.962E-02 3.0
-BIASING 0.0 0.5 8.962E-02 4.0
-BIASING 0.0 2.0 1.452E-01 5.0
-LAM-BIAS 0.0 -0.1 1.0
-EMF-BIAS 2000.0 10.0 10.0 4.0 5.0 1.0 LPBEMF
-LOW-BIAS 0.0 1.0 0.9 4.0 5.0
-*
-SCORE 208.0 211.0 210.0
-EXTRAWEI 1.0
-EVENTDAT -54.0 BIASP10EV
-*
-OUTLEVEL 1.0 7.0 1.0
-*DUMPTHEM 101. MGDRAW
-*
-RANDOMIZ 1.0
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-START 1.0 99999999. 4000.0 0.0 1.0
-STOP
+++ /dev/null
-* $Id$
-* Flugg tag $Name$
-*
-GLOBAL -1.
-TITLE
-T36ap10: 10 GeV pi- on TEST36A for FLUKA+GEANT4 simulation
-*needs t36.pemf
-DEFAULTS CALORIME
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-BEAM 10.0 PION-
-BEAMPOS -10.0 0.0 0.0 1.0
-* Switch on the new generator/heavies transported /option 13 selected
-EVENTYPE EVAP
-* Switch on low-energy neutron transport
-*LOW-NEUT 72.0 22.0 0.0196 0.0
-* Switch off all low energy neutron biasing
-LOW-BIAS 73.0 1.0 44.0
-* Switch on detailed multiple scattering for charged hadrons down to threshold
-ACCURACY 1.0 1.0
-* Pair production and bremss. down to 100 keV
-PAIRBREM +3.0 0.0 0.0001 3.0 50.0
-* Geometry
-GEOBEGIN COMBINAT
-GEOEND
-* *** Time cut-off ***
-TIME-CUT 2000.0 1.0 43.0
-* *** Materials ***
-MATERIAL 1.0 1.0079 0.089 3.0 1.0 HYDROGEN
-MATERIAL 6.0 12.01 2.26 6.0 CARBON
-MATERIAL 7.0 14.007 0.001251 7.0 NITROGEN
-MATERIAL 8.0 15.999 0.001429 8.0 OXYGEN
-MATERIAL 18.0 39.948 1.78 9.0 ARGON
-MATERIAL 13.0 29.982 2.70 10.0 ALUMINUM
-MATERIAL 26.0 55.85 7.87 17.0 IRON
-MATERIAL 26.0 55.85 7.87 18.0 IRON0
-MATERIAL 82.0 207.19 11.35 19.0 LEAD
-MATERIAL 51.0 121.75 6.62 20.0 ANTIMONY
-MATERIAL 11.3 21.0 LEADSB0
-COMPOUND -0.96 19.0 -0.04 20.0 LEADSB0
-* Be sure to be consistent with PEMF for e/pi purposes:
-MAT-PROP 804.96 21.0
-STERNHEI 6.159511 0.508 3.0000 0.246848 3.0000 0.0 21.0
-MATERIAL 11.3 22.0 LEADSB
-COMPOUND -0.96 19.0 -0.04 20.0 LEADSB
-* Be sure to be consistent with PEMF for e/pi purposes:
-MAT-PROP 804.96 22.0
-STERNHEI 6.159511 0.508 3.0000 0.246848 3.0000 0.0 22.0
-MATERIAL 1.044 23.0 SCSN38
-* Be sure to be consistent with PEMF for e/pi purposes:
-MAT-PROP 64.7 23.0
-STERNHEI 3.1997 0.1464 2.4855 0.16101 3.2393 0.0 23.0
-* Set the transport accuracy high in Scintillator:
-FLUKAFIX 0.05 23.0
-*
-COMPOUND +1.0 3.0 +1.0 6.0 SCSN38
-MATERIAL 0.001225 24.0 AIR
-COMPOUND -.9256E-03 7.0 -.2837E-03 8.0 -.01572E-3 9.0 AIR
-MATERIAL 1.17 25.0 PMMAWLS
-*COMPOUND +1.0 3.0 +0.57 6.0 +0.57 8.0 PMMAWLS
-COMPOUND -5.94e-02 10.0 -4.03e-01 11.0 -5.37e-01 12.0PMMAWLS
-* Set the correspondence for low energy neutrons *
-LOW-MAT 3.0 1.0 -2.0 293.0 HYDROGEN
-LOW-MAT 6.0 6.0 -2.0 293.0 CARBON
-LOW-MAT 7.0 7.0 14.0 293.0 NITROGEN
-LOW-MAT 8.0 8.0 16.0 293.0 OXYGEN
-LOW-MAT 9.0 18.0 -2.0 293.0 ARGON
-LOW-MAT 10.0 13.0 -2.0 293.0 ALUMINUM
-LOW-MAT 17.0 26.0 -2.0 293.0 IRON
-LOW-MAT 18.0 26.0 -2.0 293.0 IRON
-LOW-MAT 19.0 82.0 -2.0 293.0 LEAD
-LOW-MAT 20.0 51.0 -2.0 293.0 ANTIMONY
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-* Define region materials, cut-off's and step sizes *
-STEPSIZE 10.0 1.0 44.0
-* External void: *
-ASSIGNMAT 1.0 44.0
-* Dummy region 2 with black hole: *
-ASSIGNMAT 1.0 2.0
-* Surrounding vacuum: *
-ASSIGNMAT 2.0 1.0
-STEPSIZE 10.0 1.0
-* Module front vacuum: *
-ASSIGNMAT 2.0 11.0
-STEPSIZE 10.0 11.0
-* Aluminium front plate: *
-ASSIGNMAT 10.0 10.0
-STEPSIZE 0.666 10.0
-EMFCUT +1.000E-03+0.020E-03 8.0
-DELTARAY 0.0001 10.0
-* WLS: *
-ASSIGNMAT 25.0 6.0 9.0
-STEPSIZE 0.5 6.0 9.0
-EMFCUT -0.100E-03+0.010E-03 6.0 9.0
-DELTARAY 0.00005 25.0
-* Lead 1st and last EM layers: *
-ASSIGNMAT 22.0 16.0 17.0
-STEPSIZE 0.035 16.0 17.0
-EMFCUT -0.100E-03+0.010E-03 16.0 17.0
-DELTARAY 0.00005 22.0
-* Lead 1st and last HAD layers: *
-ASSIGNMAT 22.0 30.0 31.0
-STEPSIZE 0.035 30.0 31.0
-EMFCUT -0.100E-03+0.010E-03 30.0 31.0
-DELTARAY 0.00005 22.0
-* Lead 1st and last HAD last layer: *
-ASSIGNMAT 22.0 34.0 35.0
-STEPSIZE 0.035 34.0 35.0
-EMFCUT -0.100E-03+0.010E-03 34.0 35.0
-DELTARAY 0.00005 22.0
-* Lead median EM layers: *
-ASSIGNMAT 21.0 15.0
-STEPSIZE 0.25 15.0
-EMFCUT +1.000E-03+0.020E-03 15.0
-DELTARAY 0.0002 21.0
-* Lead median HAD layers: *
-ASSIGNMAT 21.0 29.0
-STEPSIZE 0.25 29.0
-EMFCUT +1.000E-03+0.020E-03 29.0
-DELTARAY 0.0002 21.0
-* Lead median HAD last layer: *
-ASSIGNMAT 21.0 33.0
-STEPSIZE 0.25 33.0
-EMFCUT +1.000E-03+0.020E-03 33.0
-DELTARAY 0.0002 21.0
-* EM Air gaps: *
-ASSIGNMAT 24.0 18.0
-STEPSIZE 0.025 18.0
-EMFCUT -0.100E-03+0.010E-03 18.0
-DELTARAY 0.00005 24.0
-* HAD Air gaps: *
-ASSIGNMAT 24.0 36.0
-STEPSIZE 0.025 36.0
-EMFCUT -0.100E-03+0.010E-03 36.0
-DELTARAY 0.00005 24.0
-* PVC EM rods: *
-ASSIGNMAT 23.0 22.0 23.0
-STEPSIZE 0.05 22.0 23.0
-EMFCUT -0.100E-03+0.010E-03 22.0 23.0
-* PVC HAD rods: *
-ASSIGNMAT 23.0 40.0 41.0
-STEPSIZE 0.05 40.0 41.0
-EMFCUT -0.100E-03+0.010E-03 40.0 41.0
-* Al EM spacers: *
-ASSIGNMAT 10.0 24.0 25.0
-STEPSIZE 1.0 24.0 25.0
-EMFCUT -0.250E-03+0.010E-03 24.0 25.0
-* Al HAD spacers: *
-ASSIGNMAT 10.0 42.0 43.0
-STEPSIZE 1.0 42.0 43.0
-EMFCUT -0.250E-03+0.010E-03 42.0 43.0
-* SCSN38 EM towers: *
-ASSIGNMAT 23.0 19.0 21.0
-STEPSIZE 0.05 19.0 21.0
-EMFCUT -0.100E-03+0.010E-03 19.0 21.0
-DELTARAY 0.00001 1.05 23.0
-* SCSN38 HAD towers: *
-ASSIGNMAT 23.0 37.0 39.0
-STEPSIZE 0.05 37.0 39.0
-EMFCUT -0.100E-03+0.010E-03 37.0 39.0
-DELTARAY 0.00001 1.05 23.0
-*
-MULSOPT -3.0 -2.0 1.0 50.0
-* *** Switch on accurate treatment of photoelectric interactions and X-ray
-* *** production
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-EMFFLUO 1.0 1.0 44.0
-* *** Switch on Rayleigh (coherent) scattering for all regions
-EMFRAY 1.0 1.0 44.0
-* Hadron and muon thresholds: *
-PART-THR 0.075 1.0 39.0 1.0
-PART-THR -0.05 9.0 9.0
-PART-THR -0.0196 8.0 8.0 1.0
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-EMFFIX 10.0 0.1 21.0 0.1 22.0 0.1
-EMFFIX 23.0 0.1 24.0 0.1 25.0 0.1
-SCORE 208.0 211.0 210.0
-*DISCARD 3.0 4.0 7.0
-EXTRAWEI 1.0
-EVENTDAT -54.0 T36AP10EV
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-*EVENTBIN +8.0 208.0 -26.0 438.0 0.1 1.0 200nsquenc
-*EVENTBIN 421.0 -.1 -1.0 1.0 1. 1.0 &
-*EVENTBIN +8.0 208.0 -26.0 438.0 0.1 1.0 200nsnoq
-*EVENTBIN 421.0 -.1 -1.0 1.0 1. 1.0 &
-*EVENTBIN +8.0 208.0 -26.0 438.0 0.1 1.0 50nsquench
-*EVENTBIN 421.0 -.1 -1.0 1.0 1. 1.0 &
-*EVENTBIN +8.0 208.0 -26.0 438.0 0.1 1.0 2 usquench
-*EVENTBIN 421.0 -.1 -1.0 1.0 1. 1.0 &
-*EVENTBIN +8.0 208.0 -26.0 438.0 0.1 1.0 2 usnoq
-*EVENTBIN 421.0 -.1 -1.0 1.0 1. 1.0 &
-*EVENTBIN +8.0 208.0 -26.0 438.0 0.1 1.0 2 usphesta
-*EVENTBIN 421.0 -.1 -1.0 1.0 1. 1.0 &
-* output on formatted file
-EVENTBIN +8.0 208.0 28.0 438. 0.1 1.0 prova
-EVENTBIN 421.0 -.1 -1.0 1.0 1. 1.0 &
-* Define time gates and quenching factors for the detectors *
-TCQUENCH 200.E-09 0.0085 1.0 1.0
-TCQUENCH 200.E-09 2.0 2.0
-TCQUENCH 50.E-09 0.0085 0.0 3.0 3.0
-TCQUENCH 0.0085 0.0 4.0 4.0
-TCQUENCH 5.0 5.0
-TCQUENCH 0.0085 0.0 6.0 6.0
-*
-*OUTLEVEL 1.0 7.0 1.0
-RANDOMIZ 1.0
-*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
-START 1.0 99999999. 4000.0 0.0 1.0
-STOP
+++ /dev/null
-* $Id$
-* Flugg tag $Name$
-*
-GLOBAL 0.0 -1.0
-TITLE
-ALAUAL: 1 MeV e- on Au,0th Born,Ae,Ap=10,Ue,Up=10keV,fudgem=0,mcs,single=1
-DEFAULTS EM-CASCA PRECISION
-*needs wa_50.pemf
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-BEAM -0.001 0.0 0.0 ELECTRON
-BEAMPOS -5.0
-*
-*
-STEPSIZE 10.00 1.0 3.0
-MULSOPT -0.0 +0.0 1.0 50.0
-*MULSOPT 0.0 +0.0 +0.0 1.0 1.0 1.0 GLOBAL
-*MULSOPT 0.0 0.0 1.0 1.0 100000.0GLOBEMF
-SCORE 208.0
-OUTLEVEL 1.0 7.0 1.0
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-*
-*
-USRBIN -11.0 208.0 +95.0 0.3 0.0168333 1st Al
-USRBIN +0.0 +0.0 30.0 20.0 &
-USRBIN -11.0 208.0 +95.0 0.3 0.0190079 Au
-USRBIN +0.0 0.0168333 30.0 20.0 &
-USRBIN -11.0 208.0 +95.0 0.3 0.1718154 2nd Al
-USRBIN +0.0 0.0190079 30.0 120.0 &
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-DUMPTHEM 101. MGDRAW
-RANDOMIZ 1.0
-GEOBEGIN COMBINAT
-GEOEND
-MATERIAL 13.0 26.982 2.70 ALUMINUB
-MATERIAL 79.0 196.97 19.3 GOLD
-ASSIGNMAT 26.0 3.0 5.0 2.0
-ASSIGNMAT 27.0 4.0 4.0
-ASSIGNMAT 2.0 2.0 6.0 4.0
-ASSIGNMAT 1.0 1.0 7.0 6.0
-SCORE 208.0
-OUTLEVEL 1.0 7.0 1.0
-EMFFIX 26.0 0.08 27.0 0.08
-EMFFLUO 1.0 26.0 27.0
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-*USRBDX 111.0 3. -197.0 1.0 2.0 1.0 AlAuAngdis
-*USRBDX 1.02E-03 1.E-05 1.0 250.0 &
-*USRBDX 111.0 3. -197.0 2.0 3.0 1.0 AuAlAngdis
-*USRBDX 1.02E-03 1.E-05 1.0 250.0 &
-*
-*USRBIN -11.0 213.0 -196.0 0.3 0.0168333 1st Al
-*USRBIN +0.0 +0.0 30.0 20.0 &
-*USRBIN -11.0 213.0 -196.0 0.3 0.0190079 Au
-*USRBIN +0.0 0.0168333 30.0 20.0 &
-*USRBIN -11.0 213.0 -196.0 0.3 0.1718154 2nd Al
-*USRBIN +0.0 0.0190079 30.0 120.0 &
-*
-USRBIN -11.0 208.0 +95.0 0.3 0.0168333 1st Al
-USRBIN +0.0 +0.0 30.0 20.0 &
-USRBIN -11.0 208.0 +95.0 0.3 0.0190079 Au
-USRBIN +0.0 0.0168333 30.0 20.0 &
-USRBIN -11.0 208.0 +95.0 0.3 0.1718154 2nd Al
-USRBIN +0.0 0.0190079 30.0 120.0 &
-*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
-DUMPTHEM 101. MGDRAW
-RANDOMIZ 1.0
-START 100.0 999999999. 3000.
-STOP
+++ /dev/null
- MEDIUM=COPPER ,STERNCID=CU
- ELEM,RHO= 8.9600E+00,NE= 1
- ASYM=CU,Z=29.,A= 63.540,PZ= 1.00000E+00,RHOZ= 6.35400E+01
- 1.43545E+00 1.01100E+00 1.00000E-01 1.00015E+04 1.00010E+04
- *** NEW PEGS4 ***
- 0.610999064504728 1.00000000000000 841.000000000000
- 875.731305716962 88244456.3125273
- T T T T
- 15 32 91 105 17 0 561 561 1122 -1 593 594 1188
- 0.500000000000000 1.77827941003892 1.43545098353905
- 0.999000000000000 0.995000000000000 0.990000000000000
- 0.980000000000000 0.970000000000000 0.950000000000000
- 0.900000000000000 0.800000000000000 0.700000000000000
- 0.600000000000000 0.500000000000000 0.400000000000000
- 0.300000000000000 0.200000000000000
- 0.849189398442293 3.417046638488989E-004 1.374659777482253E-003
- 1.736294554051905E-003 3.534837070399302E-003 3.621849464053674E-003
- 7.528103104656042E-003 2.070716451928813E-002 5.136536987858970E-002
- 6.887449456940599E-002 9.365060297461308E-002 0.129739806221673
- 0.185579882045522 0.281134629598709 1.06355337508085
- 2.270935646809018E-004 9.147194976861120E-004 1.157983063356273E-003
- 2.366899515169069E-003 2.438485490347939E-003 5.110037500969093E-003
- 1.433604474061183E-002 3.698763001280055E-002 5.161252150260809E-002
- 7.227761792431542E-002 0.102600936843072 0.149883435358615
- 0.230916863908821 0.392723106157794 1.48329121175571
- 1.469710958807230E-004 5.938007407741466E-004 7.558374401992331E-004
- 1.558806111600911E-003 1.624948966410559E-003 3.466253179612586E-003
- 1.010488221242234E-002 2.763408133904894E-002 4.041105703872714E-002
- 5.839336918950563E-002 8.486739273520906E-002 0.126270253828059
- 0.197372071296751 0.339649517555693 0.590441969025817
- 2.09846182070850 1.156266955472039E-004 4.701551746537645E-004
- 6.050318664446227E-004 1.268181102212379E-003 1.347400984974657E-003
- 2.952364349625539E-003 9.021889553170867E-003 2.595796976106264E-002
- 3.896840856051469E-002 5.688455699831263E-002 8.289070665439385E-002
- 0.123097794308485 0.191846026324652 0.328997313945778
- 0.569177646128132 0.664860748300540 2.90250404457438
- 1.024442140764854E-004 4.200325530249554E-004 5.487971869190434E-004
- 1.182527691478965E-003 1.303459112214553E-003 2.982421315261351E-003
- 9.632893329399128E-003 2.885535040422161E-002 4.347933072915289E-002
- 6.257031851422155E-002 8.977610264484445E-002 0.131641405028440
- 0.202512753195387 0.341933028052638 0.580733569068956
- 0.669874645683716 0.734954965850430 3.82383785729629
- 9.539528345949113E-005 3.973794698778879E-004 5.336263600625372E-004
- 1.203617129944837E-003 1.399417040810306E-003 3.361125867683427E-003
- 1.140249954689701E-002 3.490381241942736E-002 5.171072946134427E-002
- 7.250725150214520E-002 0.101444074372393 0.145463178432099
- 0.219933505076904 0.365407394835647 0.606927107579612
- 0.685855036999265 0.743097672487427 0.778195033431291
- 4.77078281126595 9.098741816799116E-005 3.932703750371220E-004
- 5.570433717748470E-004 1.330691792530365E-003 1.613606852183603E-003
- 4.019457247313805E-003 1.414693786039926E-002 4.332698104940064E-002
- 6.217543595317791E-002 8.433329127362829E-002 0.114913402414051
- 0.161182361416185 0.238375129378176 0.386787930916409
- 0.627249512043885 0.697399685659242 0.749811551972454
- 0.782190517791951 0.800885016479982 5.71961570807765
- 8.858137606242726E-005 4.023067724507411E-004 6.084477588358875E-004
- 1.554600823328526E-003 1.979305088229608E-003 5.105901654503405E-003
- 1.818621125588336E-002 5.364325140857994E-002 7.347930997717561E-002
- 9.648062464785716E-002 0.127801939636711 0.176067456366905
- 0.254364137975870 0.402410660389435 0.641696222430587
- 0.705896105973484 0.755405699808674 0.786241076180617
- 0.804084885271972 0.814118983280487 6.61930486947425
- 8.838898277956635E-005 4.390136663249890E-004 7.282738275729352E-004
- 1.978905302768468E-003 2.571960964546340E-003 6.601444018286727E-003
- 2.287768178938359E-002 6.464564272526702E-002 8.497337931600504E-002
- 0.107918207372754 0.139219818698239 0.187260403588952
- 0.265518771246845 0.413396093401334 0.649583354019129
- 0.708169336244291 0.754517857579366 0.783572239418159
- 0.800414442215355 0.809888462359895 0.814941192736994
- 7.50225392331694 9.031014303484156E-005 5.045983680864370E-004
- 9.196521038862946E-004 2.614268734673931E-003 3.412625538025753E-003
- 8.587280608058089E-003 2.834139914104486E-002 7.540567147530604E-002
- 9.519395072986815E-002 0.117479516066840 0.147945319876034
- 0.194355307395947 0.272107360759077 0.420785873097472
- 0.654402350683423 0.709111925738535 0.753959647669214
- 0.782311049717823 0.798783804245821 0.808056373841334
- 0.813002427567950 0.814883209815486 8.36718336617754
- 9.452390588905237E-005 6.110814692623195E-004 1.209346430451149E-003
- 3.488571170722645E-003 4.453330650471174E-003 1.079895191281690E-002
- 3.381797859893438E-002 8.529656933205576E-002 0.103432460685645
- 0.124407986521690 0.154273048644801 0.200126603576783
- 0.277210312358769 0.425337699767445 0.658281114909872
- 0.710989573809871 0.754128417784140 0.781394885458524
- 0.797235547234056 0.806150274563570 0.810897547878908
- 0.812676881139201 0.810870658373656 9.19908185067994
- 1.048200102744689E-004 8.067235225175920E-004 1.626133636664516E-003
- 4.526025299438454E-003 5.572441479429311E-003 1.300805397702403E-002
- 3.883064217426056E-002 9.305928973257950E-002 0.109332998245010
- 0.129343538887382 0.158528853761778 0.204038061129440
- 0.281092109206013 0.429503620551089 0.662427627471569
- 0.712965357901367 0.753954236645052 0.779805366313455
- 0.794814339234128 0.803254614789722 0.807736895703925
- 0.809375153353538 0.807515286066143 0.797859661588136
- 9.93908024490647 1.226832380132651E-004 1.088030233953903E-003
- 2.138778965315737E-003 5.647382460098194E-003 6.673068662572995E-003
- 1.502082981503695E-002 4.295712718383193E-002 9.873141717171409E-002
- 0.113266446192085 0.132563551105796 0.161242821464703
- 0.206530288591658 0.283526122483506 0.432099514483037
- 0.665187455597803 0.713419005896692 0.751133523853253
- 0.774707974950275 0.788359641168807 0.796023992597159
- 0.800071112865452 0.801475457961209 0.799519560628746
- 0.789594353528373 0.757980103807373 10.5337616466434
- 1.499163718521550E-004 1.436781644247437E-003 2.693630017065568E-003
- 6.752073155552322E-003 7.653704084451293E-003 1.667462103567931E-002
- 4.599237362648133E-002 0.102446762677853 0.115657265742125
- 0.134393374884743 0.162794128378716 0.207915640300753
- 0.284758539495507 0.433124650947478 0.666459976144736
- 0.712718499921479 0.746555046646125 0.767346072297662
- 0.779325567490272 0.786032008313119 0.789541385550467
- 0.790656041449228 0.788204919333822 0.777376201413177
- 0.745221770672704 0.661880695048058 10.9596728105599
- 1.907219004113227E-004 1.855037040574427E-003 3.228075626018395E-003
- 7.692816648890472E-003 8.451734181505454E-003 1.792219893976180E-002
- 4.806265340016003E-002 0.104640407772906 0.117170004523048
- 0.135543355137737 0.163757455987703 0.208764330546628
- 0.285540480860605 0.433942449814374 0.667603502676306
- 0.712391699603643 0.743521745273253 0.762378840128340
- 0.773197445577125 0.779239938139424 0.782375401199902
- 0.782901884506977 0.779377413045056 0.767195624474139
- 0.733068585784848 0.651279542320367 0.488379465450155
- 11.1737840435805 2.447712800070303E-004 2.300772494849215E-003
- 3.719553831393843E-003 8.453169928106267E-003 9.007614531561878E-003
- 1.873356407548842E-002 4.935618207993655E-002 0.105957539711240
- 0.118063763687667 0.136222697600904 0.164309684542900
- 0.209273687013396 0.285964630366082 0.434388263539496
- 0.668222900178746 0.711929955747499 0.741332729780717
- 0.758970603000464 0.769060413659503 0.774684801668762
- 0.777200940059667 0.776606624096164 0.771481898520790
- 0.756992837989485 0.723825394103351 0.641533400381960
- 0.481089316709196 0.274856333001130 11.2165567530948
- 3.094275639678021E-004 2.728216165337644E-003 4.127488086111215E-003
- 9.014646380617381E-003 9.386532417820779E-003 1.927353694566818E-002
- 5.010702468902749E-002 0.106851361298563 0.118595500735770
- 0.136559112832996 0.164630342322504 0.209572764443509
- 0.286226828325535 0.434491947264056 0.668289646179944
- 0.711625660681953 0.740105199151283 0.757070494696247
- 0.766754628071282 0.771767942928041 0.773156054273353
- 0.770896653163239 0.763400913023420 0.749304912019702
- 0.711564230636179 0.624133451063270 0.468034413015250
- 0.271078628592389 0.117499196127710 11.1346267996325
- 3.816070957974347E-004 3.099520961025463E-003 4.426896986350157E-003
- 9.378881218470718E-003 9.629007451596366E-003 1.962096732399266E-002
- 5.049379726534314E-002 0.107412968956004 0.118887809240190
- 0.136696859535308 0.164802426310120 0.209723734562842
- 0.286384318599569 0.434426340278942 0.668110988633480
- 0.711467907755606 0.739441161918315 0.756014218117045
- 0.765094368703316 0.768969361151186 0.768675512692838
- 0.764024898089184 0.756697310876820 0.736112434188013
- 0.686531493639630 0.592010121085243 0.443931619806827
- 0.263577839196415 0.116481304170183 4.212112382279570E-002
- 10.9291189794085 4.564679333627508E-004 3.394468373902054E-003
- 4.623522841367416E-003 9.595481631978420E-003 9.764464929756925E-003
- 1.980085009046782E-002 5.076218409998651E-002 0.107581503109898
- 0.119029062561486 0.136880617522192 0.164855159260860
- 0.209746296946316 0.286340176697661 0.434467895304547
- 0.668259481401640 0.710943892294436 0.738102965922842
- 0.753776711573978 0.761521552500092 0.763617094710255
- 0.760925817252490 0.756283068417209 0.741807880405808
- 0.706541255716384 0.641825010334459 0.542188734365535
- 0.406562127905258 0.249168268246134 0.114620648144962
- 4.180075564615672E-002 1.387556326711029E-002 10.6667143293398
- 5.959822195272089E-004 3.606988834220608E-003 4.771598480201761E-003
- 9.771082091731871E-003 9.887939640181486E-003 1.999955323176958E-002
- 5.111114971480737E-002 0.108197208987232 0.119637363008827
- 0.137570266808845 0.165632130426384 0.210675440007075
- 0.287640268815787 0.436530588985946 0.671567379088255
- 0.714046646382775 0.740119814711873 0.754190876646184
- 0.760001850956293 0.759521494955472 0.756702709462926
- 0.744697646632946 0.714408724151974 0.660459115541665
- 0.582233337329666 0.481649189876793 0.361162626317997
- 0.228911384990201 0.111467420190292 4.161215791750126E-002
- 1.388025313347238E-002 4.454139800966516E-003
- -0.494305802109391 -0.470511152532346 -0.461224865905246
- -0.489597608488607 -0.516281952519905 -0.719332305250697
- -0.726556871204797 -0.121814774109748 -1.854596302268764E-002
- 5.542869250736178E-002 6.162763639126848E-002 4.196080083416576E-002
- 5.690648903709493E-002 0.100218602630435 -0.718323873230665
- -0.728510526242299 -0.722281095929796 -0.817658064595345
- -0.990484486086475 -1.15304240154705 -0.682105770761454
- -8.450484669191266E-002 1.003126414693042E-002 4.709420716816049E-002
- 5.447530041339482E-002 5.232437182123938E-002 5.463578176089638E-002
- 8.461731874405207E-002 0.135856742320561 -1.78361163837314
- -2.05424263799056 -2.22864887624572 -1.90265537255494
- -1.93583687628890 -1.44734443755627 -0.242222565396714
- 6.251128481121931E-002 4.769507159033101E-002 5.255718693963455E-002
- 5.870385486061108E-002 5.806317820591803E-002 5.746718442494692E-002
- 7.461747554490936E-002 0.126726201751702 0.142926314653302
- 0.180319339654099 0.225769162967997 0.538681372708018
- 0.176658896778531 -1.09895135566803 -0.553448789338112
- 0.310224921580050 0.287438753974524 0.138448352372563
- 0.105520285542344 9.420271547399407E-002 6.975795680198912E-002
- 6.378713666986446E-002 8.013498382205995E-002 0.122098867577380
- 0.137488805328502 0.160045087064814 -2.83508281941452
- -2.83643536749677 -3.79785144265855 -2.34755381299702
- -8.946638691751330E-002 0.893235325875784 0.936656249276755
- 0.595661709387612 0.316834481526184 0.163437575542153
- 0.106437398689035 9.207915069893489E-002 9.021194915880107E-002
- 9.941914333319306E-002 0.119468617725282 0.129976737754960
- 0.156743234568082 0.136522377928976 -3.73846579357945
- -3.69340845555231 -5.38151989232416 -2.66389424942828
- 1.34471763558930 2.23243855889524 1.39492669759537
- 0.937267710235826 0.490819168947819 0.236876576840910
- 0.129774951954739 8.300000998659308E-002 0.101333887491919
- 0.134967555355432 0.131518985274160 0.127328940628901
- 0.150170431039058 0.129936170355823 0.105429333458434
- 1.76273740612081 1.74134437561071 1.50870920419975
- 1.70909653529404 2.37465947269156 2.37812944586064
- 1.83588857555202 1.28000396780873 0.654401782105346
- 0.265209233657569 0.136932021921130 0.109359519934209
- 0.125204415608367 0.147814161926591 0.124532076110145
- 0.115992304621932 0.144601915688082 0.125652115887888
- 0.101908576494529 8.210886177505250E-002 3.07837197551526
- 3.12230663025336 3.05983246896855 3.12553138401831
- 3.25758736971195 3.11178499258743 2.56796657993217
- 1.60807687922357 0.724439683860962 0.296387804578199
- 6.111328082772458E-002 0.176132325517374 0.139543812573732
- 0.126499164353643 0.113044566679011 0.105495414710891
- 0.140624708869682 0.122933275023937 9.982938388367715E-002
- 8.043716702811693E-002 7.013782720035837E-002 7.63445835087088
- 7.50343385504622 7.15917446212399 6.74176595146814
- 5.35678116540714 4.07839226680043 2.94499173854266
- 1.78976456717262 0.815578537343986 0.289599811130219
- 3.635040163003180E-002 0.147564435400332 0.126111585094187
- 0.120857624914477 0.100555447679907 9.688164641557115E-002
- 0.136872836101129 0.120197684532868 9.763521186682117E-002
- 7.867639789962702E-002 6.867806766567207E-002 7.868446321544685E-002
- 11.6456463389307 11.3807941577023 10.5360657694945
- 9.90969494835700 7.46041664100540 5.23815220134980
- 3.43903624559408 1.85414598847768 0.832168629027338
- 0.271854549748802 1.781338087987258E-003 1.084016043627587E-003
- 8.546565879954672E-002 0.131883370073553 8.766643837738650E-002
- 8.655442747194717E-002 0.134543127881904 0.118981924999175
- 9.677825104510007E-002 7.801464235089235E-002 6.820002749392100E-002
- 7.803542552851178E-002 0.122859695233419 16.5573089460091
- 16.2051066992038 14.6938713706297 13.9214481456702
- 9.62443517166695 6.06855123636534 3.61351767088383
- 1.89927342266403 0.790924974014026 3.965986939006751E-002
- -0.173176484956848 -8.697317010396692E-002 4.897517393380278E-002
- 0.117333690661442 8.368400997490091E-002 8.609971876795552E-002
- 0.132766561483349 0.117263688334628 9.528008985058771E-002
- 7.682835040298698E-002 6.730111153458301E-002 7.737751530781513E-002
- 0.122187908622631 0.221297266386429 36.8431865732043
- 32.6017654592268 21.7786446860927 16.1684666790132
- 11.1474774101385 6.33054133339316 3.76622388242660
- 1.84461940763767 0.640506925302072 -3.552366771802097E-002
- -7.232830512118751E-002 -0.150477706561463 2.848600685821762E-002
- 0.114645777621235 8.644018087215975E-002 8.734892867599889E-002
- 0.130661156362011 0.114999345838127 9.337631717207133E-002
- 7.524188159656728E-002 6.615377367745198E-002 7.668571461632470E-002
- 0.121707825560116 0.220659663777446 -0.250542014608208
- 64.6878584883385 50.9029858407569 27.6937952402876
- 17.9846686006653 12.0325451829444 6.37142121890582
- 3.72949054725140 1.70090742716229 0.364761257527520
- 5.865019289760022E-002 -2.289655329982688E-002 -0.116992398327483
- 1.296541498229797E-002 0.113397169367209 9.737170725761535E-002
- 9.459650831394656E-002 0.127896001867479 0.111374375521172
- 9.012231198275890E-002 7.261374254520794E-002 6.430585096740292E-002
- 7.570384167959718E-002 0.121251573593287 -0.173397564272414
- -0.241897438659797 -0.277561839515280 99.2309413761823
- 68.3362710721001 31.5108601038243 19.9616924080261
- 12.3224841555360 6.15342152787925 3.48547905347723
- 1.47381551273788 0.136735068296946 0.110402032444553
- 6.771806409013054E-004 -9.949998568753195E-002 3.361418654800277E-003
- 0.108522172516394 0.111915424416138 0.107294169168757
- 0.124772692447458 0.106662917671979 8.580012121904380E-002
- 6.912760450699351E-002 6.190276056788604E-002 7.450232321864664E-002
- 0.137039396524594 -0.168985630984819 -0.240702985193956
- -0.271271756276407 -0.147163423312249 143.368848395971
- 86.7414402153656 33.8075925386188 19.6258844142982
- 12.7271873538992 5.45425418826646 3.19583147230973
- 1.09999175780579 -0.151471225148918 0.138450152249820
- 1.601759460205103E-002 -8.903071211515096E-002 -6.023469867846961E-003
- 0.107601626303331 0.122338140150097 0.115762886142086
- 0.122358233678549 0.103071708147245 8.249458464262924E-002
- 6.659279848327361E-002 6.010942169249330E-002 9.573252915383225E-002
- 0.136855331933094 -0.173420975463381 -0.223465138934253
- -0.273932290026083 -0.147597002644987 0.204302762094146
- 193.613162158656 101.376847528222 34.5066656556216
- 19.4508384024496 11.9646639626493 4.27626748467415
- 2.85304456722490 -0.374330442902627 -0.166253411907084
- 0.155316649684429 2.504160181843863E-002 -8.528784273579662E-002
- -1.735053143050077E-002 0.108027001013475 0.129192313638637
- 0.120552090423606 0.120674381931539 0.100648972494923
- 8.038888498415039E-002 6.485640180272426E-002 8.388542594285082E-002
- 9.687210261125938E-002 -0.127755688757693 -0.140547750019623
- -0.239442105277533 -0.265596814477817 -0.145963634651575
- 0.206599849720181 0.578683900352557 249.267296692751
- 111.043721316724 33.6844971182258 18.6698307176740
- 10.2235145206149 3.65296673862650 2.27485011140157
- 6.874068427431788E-002 -0.164186640778939 0.166042189057088
- 2.716307097768371E-002 -8.168309359065400E-002 -1.595102296173622E-002
- 0.105377650906351 0.131768711641732 0.124093622660808
- 0.119836857732318 9.929398958918192E-002 7.912782628590995E-002
- 8.702612447028175E-002 8.528380392817894E-002 0.125055988858081
- -7.032434231189182E-002 -0.179530321849191 -0.236333055182177
- -0.238824583943203 -0.135567523712235 0.169457793922654
- 0.589437924738425 0.833944918608770 312.442981979559
- 115.203485502087 31.3034727182859 16.7672710639169
- 8.12326620786712 3.83120921101380 -0.990250879374859
- 1.116514523395962E-002 -0.156035841154398 0.171667067099937
- 2.625865436329348E-002 -7.841813696455786E-002 -6.898829095440717E-003
- 0.101262369882722 0.132251297780601 0.127062861386335
- 0.119470627023733 9.848408317339995E-002 9.845714983924851E-002
- 8.816211436529688E-002 0.115603496768573 2.372771486428774E-002
- -0.142035336580027 -0.193130175119275 -0.210623533742721
- -0.196708345147935 -0.116316169251848 0.126656992001878
- 0.578573475063938 0.820831885966650 0.946950153903023
- 386.576694827933 113.429332903571 27.2580233810707
- 13.9405706423064 6.86616754473051 -2.35619325200601
- -0.697502823293261 -3.945027080058403E-002 -0.168205070393473
- 0.174889574622271 3.293538099679617E-002 -7.914275800703328E-002
- -1.731591766315214E-002 0.104081023690476 0.136322813344965
- 0.128618298775218 0.118563028920428 0.114622121663888
- 9.833876416058843E-002 0.117016050176856 1.636222013311262E-002
- -0.127382466774912 -0.170300076073484 -0.177894958491776
- -0.166550786683758 -0.145506037232891 -8.980468804399233E-002
- 8.776288391179135E-002 0.489400274046596 0.807496473635281
- 0.936230268167007 0.982723477597790 378.241376538593
- 100.815901249414 23.0082507716080 11.6570757410691
- 5.21819917404764 2.44735631163897 0.191343777743010
- -4.903536712791148E-002 -0.171740214272502 0.174991598115037
- 3.525221669275762E-002 -7.587930083703887E-002 -1.683031115807512E-002
- 0.105051439463215 0.139910618058070 0.131917646001462
- 0.133170678291659 0.113079886985310 0.124098176838317
- 8.757321836153823E-003 0.145782912183915 -0.163050980201789
- -0.162601833820956 -0.139733540270987 -0.112464513611217
- -9.171772409911902E-002 -5.860374914094959E-002 5.983160799826208E-002
- 0.383740623131072 0.820528828800634 0.931520267848478
- 0.982722975501078 0.994487574375447
- 1.02628030961715 1.03572933531917 1.04785638640801
- 1.07318365780219 1.10007642110578 1.15906849135782
- 1.33126099553276 1.77497042693800 2.39298793623706
- 3.28127946932280 4.59776177598211 6.72296859840004
- 10.5668625846004 18.9281778670557 0.453560188258482
- 0.454814831685763 0.459917146897453 0.466489567551909
- 0.480359255897494 0.495135028200734 0.527331376020316
- 0.622779422630411 0.870351352228302 1.21384962595784
- 1.70859591736945 2.44960523728416 3.65130701525833
- 5.82603440232494 10.5967701280954 22.6543628558898
- 0.164966921723971 0.165625327801312 0.168302783473291
- 0.171743428819740 0.178953230325109 0.186626851656812
- 0.203512653348098 0.252519787665081 0.375252101417614
- 0.543245389740821 0.785708605889912 1.15007859057137
- 1.74245897930740 2.81653625875180 5.18052986239215
- 11.2105920053246 22.6051223134551 7.288791507291369E-002
- 7.336976261174800E-002 7.530928558068856E-002 7.775586933185522E-002
- 8.267759771656265E-002 8.780176969021522E-002 9.908257730495798E-002
- 0.129711705336532 0.201530650650596 0.296379258789887
- 0.431453821211419 0.631381828248050 0.953900382930376
- 1.53730288714519 2.81663683236063 6.06540803293764
- 12.2323344924165 23.4179582181358 3.625519024572351E-002
- 3.661511941258479E-002 3.808716841657836E-002 4.000234567263755E-002
- 4.416926894521190E-002 4.857085884022008E-002 5.751097460596717E-002
- 7.956845040426089E-002 0.126826502756179 0.184982979120863
- 0.264726196576797 0.382166210333704 0.570345688008142
- 0.906540752446142 1.63315194862084 3.45237950878987
- 6.90760777247858 13.1764299600799 24.3630377519607
- 1.892415131134013E-002 1.923457003896833E-002 2.051331801361547E-002
- 2.219707713148001E-002 2.600890878299288E-002 2.996050146100037E-002
- 3.715921895444897E-002 5.386670472588886E-002 8.587565135259928E-002
- 0.122198752388007 0.170433756221695 0.240138192116325
- 0.351199093096838 0.549738644064802 0.972848560946663
- 2.00581143594605 3.94781556205136 7.46102609198211
- 13.7271114602293 24.8741810009724 1.006766600840621E-002
- 1.039859855581764E-002 1.171460812604292E-002 1.334313134482623E-002
- 1.657690730946941E-002 1.969538181805881E-002 2.542011403748115E-002
- 3.790642979727399E-002 5.915227809559716E-002 8.123459975278047E-002
- 0.109957288502795 0.151325367321192 0.216661093520064
- 0.331496130394186 0.572728175857727 1.15409281994616
- 2.24680637551892 4.22364096817398 7.74947519188266
- 14.0218398881312 25.1717464879599 5.449031935326175E-003
- 5.755395784263570E-003 6.961337329060755E-003 8.423379660167346E-003
- 1.120899859355198E-002 1.379478244120266E-002 1.838395446912395E-002
- 2.717874940976198E-002 4.010683539088371E-002 5.309394408196858E-002
- 6.978034258716714E-002 9.353825172537299E-002 0.132226168502247
- 0.196007910561165 0.332419358262716 0.659487218071759
- 1.27516619306397 2.38955165548796 4.37731065935576
- 7.91355637941927 14.1998244893436 25.3630989465803
- 2.986294032347544E-003 3.299828064474898E-003 4.490574112306929E-003
- 5.847178355596403E-003 8.153153861067620E-003 1.009142273370644E-002
- 1.326651770103546E-002 1.890815339616391E-002 2.660901575433134E-002
- 3.393419027322064E-002 4.329600380315657E-002 5.665922552679339E-002
- 7.825849522919258E-002 0.114264620540141 0.190600152512316
- 0.373354828172741 0.717464158176733 1.34036139966710
- 2.45145696066793 4.42812086358568 7.94201489305239
- 14.1822917679636 25.2314790456353 1.654738365174442E-003
- 1.955991544266701E-003 3.065365069144015E-003 4.260179865945153E-003
- 6.098752694506272E-003 7.494759590151635E-003 9.563598212653410E-003
- 1.281635423419653E-002 1.705249024069920E-002 2.107648529860774E-002
- 2.617000330213447E-002 3.351186875635646E-002 4.520854027755988E-002
- 6.575674212841603E-002 0.108557890748526 0.210652370846629
- 0.403452908241218 0.752775805067930 1.37597932327927
- 2.48470521147284 4.45569798471265 7.95603443314374
- 14.1542320812862 25.0293870041599 9.212717671642668E-004
- 1.202637862980467E-003 2.201080476538474E-003 3.200512669255721E-003
- 4.535771986339263E-003 5.428145431205384E-003 6.629131207798660E-003
- 8.399626201774859E-003 1.061382241954353E-002 1.268503095897745E-002
- 1.545460853515760E-002 1.953010727036823E-002 2.603358868145745E-002
- 3.751598044674180E-002 6.153690981203183E-002 0.118962471242990
- 0.227160927449610 0.423060336618701 0.772509469698665
- 1.39419321243313 2.49936629792453 4.46209308785537
- 7.93770760638419 14.0366253819856 24.4376927107242
- 5.155381987418017E-004 8.048663434315972E-004 1.683652506230009E-003
- 2.385777401193673E-003 3.249556001492471E-003 3.759251845415908E-003
- 4.413028243913369E-003 5.301034445028307E-003 6.385390785207515E-003
- 7.472422043121398E-003 8.971975365045922E-003 1.123290761246483E-002
- 1.488280519150187E-002 2.134749092350524E-002 3.488327910470466E-002
- 6.720900086829987E-002 0.127878613834769 0.237595319450774
- 0.433270322598770 0.781371506916852 1.40019154759785
- 2.49918783280148 4.44534030643853 7.86062198612490
- 13.6862287601452 22.7637077901232 2.889937550598433E-004
- 5.739946485798593E-004 1.277628549747239E-003 1.725647094372254E-003
- 2.227794710701443E-003 2.493261777346559E-003 2.818628225845784E-003
- 3.235666183865906E-003 3.752254877661209E-003 4.326819252615997E-003
- 5.146464612168209E-003 6.406184958768156E-003 8.455085551254616E-003
- 1.209283429399955E-002 1.971306249541293E-002 3.790976874002359E-002
- 7.178397082169330E-002 0.132890209389178 0.241823516853437
- 0.435598507863246 0.780068281622680 1.39183155753375
- 2.47518140637059 4.37640881214727 7.61023433112830
- 12.6559142792022 18.7031846271414 1.621816667027945E-004
- 4.236751737137879E-004 9.279446311453316E-004 1.197670281332155E-003
- 1.460857817351064E-003 1.586181658221620E-003 1.733806223855962E-003
- 1.919832872130996E-003 2.166780441466498E-003 2.472824344582404E-003
- 2.926611425333023E-003 3.631224257266256E-003 4.780457585737200E-003
- 6.823672053880687E-003 1.110285639771495E-002 2.135216477398128E-002
- 4.021628911561038E-002 7.412517158269895E-002 0.134537090680887
- 0.241988601277924 0.432999312521096 0.772225402710129
- 1.37295201250248 2.42397696806919 4.21219606417304
- 6.98826539667547 10.3109157408904 12.0196986725139
- 9.109708626749603E-005 3.201669642644581E-004 6.469310207842131E-004
- 7.870493664669141E-004 9.202454232290915E-004 9.729229837268646E-004
- 1.036589195126390E-003 1.114761035923955E-003 1.237932277000830E-003
- 1.404573608438075E-003 1.657294133217070E-003 2.051816106786847E-003
- 2.697243943113811E-003 3.845955139970538E-003 6.253416032474429E-003
- 1.202364462075342E-002 2.256098730995579E-002 4.145555627455164E-002
- 7.510361740653050E-002 0.134947108529340 0.241326214717627
- 0.430250245948474 0.763770910742107 1.34726779154874
- 2.33321805580351 3.86457583567012 5.70804856944837
- 6.64845346417834 5.69581409349461 5.121893634368665E-005
- 2.384928438212142E-004 4.313608252454213E-004 4.993733439525764E-004
- 5.580060016807216E-004 5.777601734075651E-004 6.050631710999452E-004
- 6.392028207469376E-004 7.024607049428123E-004 7.946659760347489E-004
- 9.356096264386916E-004 1.157078601295230E-003 1.519893493448731E-003
- 2.165002313926786E-003 3.520068634819039E-003 6.765086909988793E-003
- 1.266297713730843E-002 2.322273552856853E-002 4.202298703872537E-002
- 7.545791642605074E-002 0.134892098903654 0.240115728584909
- 0.425845293662378 0.748343259103744 1.29498545964159
- 2.14417217406410 3.15817036578428 3.69142558079176
- 3.16455702744749 2.13862643244894 2.880646583692442E-005
- 1.714199082710521E-004 2.756186767319492E-004 3.063571904275478E-004
- 3.292906238251384E-004 3.369726460802589E-004 3.477885639376647E-004
- 3.637119409156026E-004 3.974795100962286E-004 4.480969425669003E-004
- 5.272769630013910E-004 6.517972306880800E-004 8.556647504060090E-004
- 1.218328346070395E-003 1.979275293103630E-003 3.804866255303710E-003
- 7.113230750536667E-003 1.303103253489854E-002 2.356517130122146E-002
- 4.229885438912966E-002 7.549681865472709E-002 0.134259491298690
- 0.237183281119439 0.416627842052052 0.719378849523866
- 1.17851936648114 1.72266976719615 2.03346642587577
- 1.76215293455739 1.19689876780424 0.718198412703960
- 1.619776777419549E-005 1.182296386955726E-004 1.695384054428048E-004
- 1.822504045095869E-004 1.905907695019356E-004 1.940758164873096E-004
- 1.979613441695457E-004 2.058289023361006E-004 2.244251403848445E-004
- 2.521960695378520E-004 2.968504332088356E-004 3.668798509405220E-004
- 4.814119716560596E-004 6.854783351426832E-004 1.112308839778837E-003
- 2.139707425570724E-003 3.997799867910346E-003 7.319202782213715E-003
- 1.323087105394888E-002 2.371171994464270E-002 4.228050244040247E-002
- 7.489420174472702E-002 0.132264637811407 0.231646717596495
- 0.394175485343282 0.632976544083363 0.914941018265144
- 1.09593376471096 0.979769958017444 0.669035464694675
- 0.403099140317204 0.231402115189048 8.638767663991825E-006
- 7.826537487936863E-005 1.011857629531652E-004 1.058867031413109E-004
- 1.090966203181663E-004 1.103601212086103E-004 1.122482355720289E-004
- 1.160662667207323E-004 1.262320877574119E-004 1.421252717577717E-004
- 1.670239880883371E-004 2.063263027432630E-004 2.707571732649823E-004
- 3.853389855605248E-004 6.257573379725212E-004 1.203077108115907E-003
- 2.245147910067375E-003 4.106756982831484E-003 7.409923058826972E-003
- 1.326445248961913E-002 2.355662255341170E-002 4.171134010821154E-002
- 7.342644971037034E-002 0.126548250973767 0.209955977266937
- 0.328580804350521 0.468939232146495 0.571219212025101
- 0.535151461759327 0.372248325720426 0.225661242680377
- 0.129692848652285 7.341553696061404E-002 1.771765698797410E-005
- 4.990314847390114E-005 5.929627969284019E-005 6.106609363657351E-005
- 6.221520584765141E-005 6.280210983471284E-005 6.364763620146623E-005
- 6.568325377127339E-005 7.134388060167624E-005 8.033564682937378E-005
- 9.438626244696898E-005 1.165489665726852E-004 1.529193494251748E-004
- 2.177151024270758E-004 3.535591547286173E-004 6.798873315368561E-004
- 1.267230991778807E-003 2.312510587803473E-003 4.166306781572246E-003
- 7.426604625458251E-003 1.318185199562715E-002 2.326353173903793E-002
- 4.027882535965311E-002 6.755842650446583E-002 0.108618960951261
- 0.165422206308709 0.233159274936030 0.288751008259468
- 0.285606227760049 0.208197396896506 0.126826380337919
- 7.302093871347584E-002 4.133514041149848E-002 2.329335753165929E-002
- 0.454338340939134 0.000000000000000E+000 0.000000000000000E+000
- 1.225724997345843E-007 5.878819436061741E-006 2.744348369898673E-005
- 1.883443277070690E-004 1.778585026215318E-003 1.153007388832872E-002
- 2.496406797711890E-002 4.368944706886852E-002 7.108753940066155E-002
- 0.113819316402540 0.187247521972059 0.687655438706215
- 0.000000000000000E+000 0.000000000000000E+000 2.269101906358851E-007
- 8.111786799447408E-006 3.296502702501841E-005 2.015794082270257E-004
- 1.702553704161646E-003 1.025875028716568E-002 2.167838142673472E-002
- 3.770344576686584E-002 6.134382528532480E-002 9.854982644095579E-002
- 0.162981296377584 0.293194476285180 1.12904789109976
- 0.000000000000000E+000 0.000000000000000E+000 6.847719063798405E-007
- 1.591839718571880E-005 5.278480175923643E-005 2.732788236575117E-004
- 1.959652850404258E-003 1.053384711075480E-002 2.119039345920095E-002
- 3.596158088389166E-002 5.769089830030873E-002 9.195158230940953E-002
- 0.151474007182676 0.271858782691193 0.486084479517408
- 1.75464585856300 0.000000000000000E+000 0.000000000000000E+000