[807] | 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 | |
---|