Newer
Older
# ==========================================================================
# AIDA Detector description implementation
# --------------------------------------------------------------------------
# Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
# All rights reserved.
#
# For the licensing terms see $DD4hepINSTALL/LICENSE.
# For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
#
# ==========================================================================
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
from g4units import ns
def run():
import DigiTest
digi = DigiTest.Test(geometry=None)
attenuation = digi.attenuation
input = digi.input_action('DigiParallelActionSequence/READER')
#========================================================================================================
signal = input.adopt_action('DigiROOTInput/SignalReader', mask=0x0, input = [digi.next_input()])
digi.info('Created input.signal')
#========================================================================================================
digi.info('Creating spillover sequence for EARLIER bunch crossings.....')
#========================================================================================================
spillover = input.adopt_action('DigiSequentialActionSequence/Spillover-25')
evtreader = spillover.adopt_action('DigiROOTInput/Reader-25ns', mask=0x1, input = [digi.next_input()])
attenuate = spillover.adopt_action('DigiHitAttenuatorExp/Attenuator-25ns',
t0 = -25 * ns, masks = [evtreader.mask], containers = attenuation)
hist_drop = spillover.adopt_action('DigiHitHistoryDrop/Drop-25ns', masks = [evtreader.mask])
digi.info('Created input.spillover25')
#========================================================================================================
spillover = input.adopt_action('DigiSequentialActionSequence/Spillover-50')
evtreader = spillover.adopt_action('DigiROOTInput/Reader-50ns', mask=0x2, input = [digi.next_input()])
attenuate = spillover.adopt_action('DigiHitAttenuatorExp/Attenuator-50ns',
t0 = -50 * ns, masks = [evtreader.mask], containers = attenuation)
hist_drop = spillover.adopt_action('DigiHitHistoryDrop/Drop-50ns', masks = [evtreader.mask])
digi.info('Created input.spillover50')
#========================================================================================================
spillover = input.adopt_action('DigiSequentialActionSequence/Spillover-75')
evtreader = spillover.adopt_action('DigiROOTInput/Reader-75ns', mask=0x3, input = [digi.next_input()])
attenuate = spillover.adopt_action('DigiHitAttenuatorExp/Attenuator-75ns',
t0 = -75 * ns, masks = [evtreader.mask], containers = attenuation)
hist_drop = spillover.adopt_action('DigiHitHistoryDrop/Drop-75ns', masks = [evtreader.mask])
digi.info('Created input.spillover75')
#========================================================================================================
digi.info('Creating spillover sequence for LATER bunch crossings.....')
#========================================================================================================
spillover = input.adopt_action('DigiSequentialActionSequence/Spillover+25')
evtreader = spillover.adopt_action('DigiROOTInput/Reader+25ns', mask=0x4, input = [digi.next_input()])
attenuate = spillover.adopt_action('DigiHitAttenuatorExp/Attenuator+25ns',
t0 = 25 * ns, masks = [evtreader.mask], containers = attenuation)
hist_drop = spillover.adopt_action('DigiHitHistoryDrop/Drop+25ns', masks = [evtreader.mask])
digi.info('Created input.spillover25')
#========================================================================================================
spillover = input.adopt_action('DigiSequentialActionSequence/Spillover+50')
evtreader = spillover.adopt_action('DigiROOTInput/Reader+50ns', mask=0x5, input = [digi.next_input()])
attenuate = spillover.adopt_action('DigiHitAttenuatorExp/Attenuator+50ns',
t0 = 50 * ns, masks = [evtreader.mask], containers = attenuation)
hist_drop = spillover.adopt_action('DigiHitHistoryDrop/Drop_50ns', masks = [evtreader.mask])
digi.info('Created input.spillover50')
#========================================================================================================
spillover = input.adopt_action('DigiSequentialActionSequence/Spillover+75')
evtreader = spillover.adopt_action('DigiROOTInput/Reader+75ns', mask=0x6, input = [digi.next_input()])
attenuate = spillover.adopt_action('DigiHitAttenuatorExp/Attenuator+75ns',
t0 = 75 * ns, masks = [evtreader.mask], containers = attenuation)
hist_drop = spillover.adopt_action('DigiHitHistoryDrop/Drop+75ns', masks = [evtreader.mask])
digi.info('Created input.spillover75')
#========================================================================================================
event = digi.event_action('DigiSequentialActionSequence/EventAction')
combine = event.adopt_action('DigiContainerCombine/Combine', masks=[0x0, 0x1, 0x2, 0x3], deposit_mask=0xFEED)
combine.erase_combined = True # Not thread-safe! only do in SequentialActionSequence
evtdump = event.adopt_action('DigiStoreDump/StoreDump')
digi.info('Created event.dump')
#========================================================================================================
evt_todo = 5
evt_done = digi.run(num_events = evt_todo, num_threads = 5, parallel = 3)
if evt_done == evt_todo: result = "PASSED"
else: result = "FAILED"
digi.always('%s Test finished after processing %d events.'%(result, digi.events_done(),))
digi.always('Test done. Exiting')
if __name__ == '__main__':
run()