13 <input>PartGunChargedStableGenSim.root</input>
14 <contact>Marcel Hohmann, mhohmann@student.unimelb.edu.au</contact>
15 <description> This steering script fully reconstructs particle gun events for a set of charged stable particles,
16and dumps ECL PID validation info in an ntuple and a set of histograms (one file for each pdgId).</description>
21from reconstruction
import add_reconstruction
22from ROOT
import Belle2
26for idx
in range(len(Belle2.Const.chargedStableSet)):
27 particle = Belle2.Const.chargedStableSet.at(idx)
29 if particle == Belle2.Const.deuteron:
31 pdgId = particle.getPDGCode()
32 chargedStableList.extend([pdgId, -pdgId])
36mergeChargeOfPdgIds = [
37 Belle2.Const.pion.getPDGCode(),
38 Belle2.Const.kaon.getPDGCode(),
39 Belle2.Const.proton.getPDGCode()
42basf2.set_random_seed(1337)
45main = basf2.create_path()
48inputFileName =
"../PartGunChargedStableGenSim.root"
49main.add_module(
"RootInput", inputFileName=inputFileName)
52main.add_module(
"Gearbox")
54main.add_module(
"Geometry")
57add_reconstruction(main)
61 "ECLChargedPIDDataAnalysisValidation",
62 outputFileName=
"ECLChargedPid",
63 inputPdgIdList=chargedStableList,
64 mergeChargeOfPdgIds=mergeChargeOfPdgIds)
67main.add_module(
"ProgressBar")
71main.add_module(
'Progress')