| 1 | # $Id: tiaraSpecifications.py,v 1.3 2003/06/20 12:41:07 dressel Exp $
|
|---|
| 2 | # -------------------------------------------------------------------
|
|---|
| 3 | # GEANT4 tag $Name: $
|
|---|
| 4 | # -------------------------------------------------------------------
|
|---|
| 5 | #
|
|---|
| 6 | import CLHEP
|
|---|
| 7 | import G4Kernel
|
|---|
| 8 |
|
|---|
| 9 | ColWidth = {"43":40*CLHEP.cm,
|
|---|
| 10 | "68":80*CLHEP.cm}
|
|---|
| 11 |
|
|---|
| 12 | BinEdgesScinti = {"43":range(4,45),
|
|---|
| 13 | "68":range(6,45) + range(46,71,2)}
|
|---|
| 14 |
|
|---|
| 15 | BinEdgesBonner = {"43": [4.500E+07,
|
|---|
| 16 | 3.500E+07,
|
|---|
| 17 | 2.750E+07,
|
|---|
| 18 | 2.250E+07,
|
|---|
| 19 | 1.750E+07,
|
|---|
| 20 | 1.350E+07,
|
|---|
| 21 | 1.000E+07,
|
|---|
| 22 | 6.700E+06,
|
|---|
| 23 | 4.490E+06,
|
|---|
| 24 | 3.010E+06,
|
|---|
| 25 | 2.020E+06,
|
|---|
| 26 | 1.350E+06,
|
|---|
| 27 | 9.070E+05,
|
|---|
| 28 | 4.980E+05,
|
|---|
| 29 | 2.240E+05,
|
|---|
| 30 | 8.650E+04,
|
|---|
| 31 | 1.500E+04,
|
|---|
| 32 | 3.350E+03,
|
|---|
| 33 | 4.540E+02,
|
|---|
| 34 | 2.260E+01,
|
|---|
| 35 | 5.040E+00,
|
|---|
| 36 | 1.120E+00,
|
|---|
| 37 | 4.140E-01,
|
|---|
| 38 | 1.000E-04],
|
|---|
| 39 | "68": [ 8.000E+07,
|
|---|
| 40 | 6.500E+07,
|
|---|
| 41 | 5.500E+07,
|
|---|
| 42 | 4.500E+07,
|
|---|
| 43 | 3.500E+07,
|
|---|
| 44 | 2.750E+07,
|
|---|
| 45 | 2.250E+07,
|
|---|
| 46 | 1.750E+07,
|
|---|
| 47 | 1.350E+07,
|
|---|
| 48 | 1.000E+07,
|
|---|
| 49 | 6.700E+06,
|
|---|
| 50 | 4.490E+06,
|
|---|
| 51 | 3.010E+06,
|
|---|
| 52 | 2.020E+06,
|
|---|
| 53 | 1.350E+06,
|
|---|
| 54 | 9.070E+05,
|
|---|
| 55 | 4.980E+05,
|
|---|
| 56 | 2.240E+05,
|
|---|
| 57 | 8.650E+04,
|
|---|
| 58 | 1.500E+04,
|
|---|
| 59 | 3.350E+03,
|
|---|
| 60 | 4.540E+02,
|
|---|
| 61 | 2.260E+01,
|
|---|
| 62 | 5.040E+00,
|
|---|
| 63 | 1.120E+00,
|
|---|
| 64 | 4.140E-01,
|
|---|
| 65 | 1.000E-04]}
|
|---|
| 66 |
|
|---|
| 67 |
|
|---|
| 68 | tallyBinEdges = {}
|
|---|
| 69 | tallyBinEdges["43"] = [0,10,35,45]
|
|---|
| 70 | tallyBinEdges["68"] = [0,10,60,70]
|
|---|
| 71 |
|
|---|
| 72 | sourceTallyEdges = {}
|
|---|
| 73 | sourceTallyEdges["43"] = [0,36.3, 45.5,50]
|
|---|
| 74 | sourceTallyEdges["68"] = [0,60.8, 72.5,80]
|
|---|
| 75 |
|
|---|
| 76 | class Experiment(object):
|
|---|
| 77 | def __init__(self,
|
|---|
| 78 | energy,
|
|---|
| 79 | minNeutronEnergyCut,
|
|---|
| 80 | particleCut,
|
|---|
| 81 | shieldWidth,
|
|---|
| 82 | shieldMaterial):
|
|---|
| 83 | self.energy = "%(energy)d" % vars()
|
|---|
| 84 | self.minNeutronEnergyCut = minNeutronEnergyCut
|
|---|
| 85 | self.particleCut = particleCut
|
|---|
| 86 | self.shieldWidth = shieldWidth
|
|---|
| 87 | self.shieldMaterial = shieldMaterial
|
|---|
| 88 | self.binEdgesBonner = BinEdgesBonner[self.energy]
|
|---|
| 89 | self.binEdgesScinti = BinEdgesScinti[self.energy]
|
|---|
| 90 | self.colWidth = 0
|
|---|
| 91 | if self.shieldWidth <= 50.0*CLHEP.cm:
|
|---|
| 92 | self.colWidth = ColWidth[self.energy]
|
|---|
| 93 |
|
|---|
| 94 |
|
|---|
| 95 | class Specifications(object):
|
|---|
| 96 | def __init__(self,dimensions, experiment, materials):
|
|---|
| 97 | self.dimensions = dimensions
|
|---|
| 98 | self.experiment = experiment
|
|---|
| 99 | self.materials = materials
|
|---|
| 100 |
|
|---|
