def setUp(self): config = {} config['layers'] = 444 config['thickness_layer'] = 40.0 config['thickness_bar'] = 10.0 config['density_iron'] = 7.87 config['density_scint'] = 1.06 self.max_size = config['layers'] * config['thickness_layer'] / 2 self.dist = GA.composite_z(config)
physics_list = G4.G4physicslists.QGSP_BERT() gRunManager.SetUserInitialization(physics_list) physics_list.SetDefaultCutValue(1.0 * mm) physics_list.SetCutsWithDefault() # # Generator actions # pos = [] width = rc['width'] * rc['bars'] pos.append(uniform(loc=-width / 2, scale=width)) pos.append(uniform(loc=-width / 2, scale=width)) pos.append(GeneratorAction.composite_z(rc)) #config['vertex'] mom = [0, 0, config['energy_MeV']] pid = config['pid'] if GeneratorAction.is_neutrino_code(config['pid']): if config['distribution'] == 'point': raise NotImplementedError assert config['distribution'] == 'electron' or\ config['distribution'] == 'muon' or\ config['distribution'] == 'flat' log.warning('Energy argument ignored... assuming 3.8 GeV muons')