3 from beamparameters
import add_beamparameters
9 from background
import add_output
20 d = datetime.datetime.today()
21 print(d.strftime(
'job start: %Y-%m-%d %H:%M:%S\n'))
26 patha =
"/group/belle2/BGcampaigns/g4sim/" + folder
28 pathb = patha +
'/output/'
29 pathc = patha +
'/beast_hist/'
31 for i
in range(num * 10, (num + 1) * 10):
32 filenn = pathb + bgType +
'_' + accRing +
'_study-phase2-' + str(i) +
'.root'
36 histfile = pathc + bgType +
'_' + accRing +
'_' + str(num) +
'.root'
38 input_module = register_module(
'RootInput')
39 input_module.param(
'inputFileNames', fnames)
40 main.add_module(input_module)
43 histo = register_module(
'HistoManager')
44 histo.param(
'histoFileName', histfile)
45 main.add_module(histo)
47 gearbox = register_module(
'Gearbox')
48 gearbox.param(
'fileName',
'/geometry/Beast2_phase2.xml')
49 main.add_module(gearbox)
51 main.add_module(
"Progress")
56 MIP_to_PE = [27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27]
57 he3digi = register_module(
'He3Digitizer')
58 he3digi.param(
'conversionFactor', 0.303132019)
59 he3digi.param(
'useMCParticles',
False)
60 main.add_module(he3digi)
61 diadigi = register_module(
'BeamDigitizer')
62 diadigi.param(
'WorkFunction', 13.25)
63 diadigi.param(
'FanoFactor', 0.382)
64 main.add_module(diadigi)
65 pindigi = register_module(
'PinDigitizer')
66 pindigi.param(
'WorkFunction', 3.64)
67 pindigi.param(
'FanoFactor', 0.13)
68 main.add_module(pindigi)
69 clawsdigi = register_module(
'ClawsDigitizer')
70 clawsdigi.param(
'C_MIP_to_PE', MIP_to_PE)
71 main.add_module(clawsdigi)
72 qcssdigi = register_module(
'QcsmonitorDigitizer')
73 qcssdigi.param(
'ScintCell', 40)
74 qcssdigi.param(
'C_keV_to_MIP', 1629.827)
75 qcssdigi.param(
'C_MIP_to_PE', 15.0)
76 qcssdigi.param(
'MIPthres', 0.5)
77 main.add_module(qcssdigi)
78 fangsdigi = register_module(
'FANGSDigitizer')
79 main.add_module(fangsdigi)
80 tpcdigi = register_module(
'TpcDigitizer')
81 main.add_module(tpcdigi)
83 beamab_study = register_module(
'BeamabortStudy')
84 main.add_module(beamab_study)
86 claws_study = register_module(
'ClawsStudy')
87 main.add_module(claws_study)
89 fangs_study = register_module(
'FANGSStudy')
90 main.add_module(fangs_study)
92 tpc_study = register_module(
'MicrotpcStudy')
93 main.add_module(tpc_study)
95 he3_study = register_module(
'He3tubeStudy')
96 main.add_module(he3_study)
98 pin_study = register_module(
'PindiodeStudy')
101 qcs_study = register_module(
'QcsmonitorStudy')
106 print(
'Event Statistics:')
107 print(basf2.statistics)
109 d = datetime.datetime.today()
110 print(d.strftime(
'job finish: %Y-%m-%d %H:%M:%S\n'))