Inclusive Kinematics JB DA

options/reconstruction.py

@@ -681,13 +681,25 @@ if 'acadia' in detector_version:
             outputHitCollection="MRICHRecHits")
     algorithms.append(mrich_reco)
 
-# Electron kinematics
-electron_incl_kin = InclusiveKinematicsElectron("electron_incl_kin",
+# Inclusive kinematics
+incl_kin_electron = InclusiveKinematicsElectron("incl_kin_electron",
         inputMCParticles="mcparticles",
         inputParticles="ReconstructedParticles",
         outputData="InclusiveKinematicsElectron"
 )
-algorithms.append(electron_incl_kin)
+algorithms.append(incl_kin_electron)
+incl_kin_jb = InclusiveKinematicsJ("incl_kin_jb",
+        inputMCParticles="mcparticles",
+        inputParticles="ReconstructedParticles",
+        outputData="InclusiveKinematicsJB"
+)
+algorithms.append(incl_kin_jb)
+incl_kin_da = InclusiveKinematicsDA("incl_kin_da",
+        inputMCParticles="mcparticles",
+        inputParticles="ReconstructedParticles",
+        outputData="InclusiveKinematicsDA"
+)
+algorithms.append(incl_kin_da)
 
 # Output
 podout = PodioOutput("out", filename=output_rec)