Generate 1 to 10 GeV electrons

tracking/options/tracker_reconstruction.py

@@ -118,7 +118,7 @@ vtxcluster_trk_init = TrackParamVertexClusterInit("vtxcluster_trk_init",
 
 # Tracking algorithms
 trk_find_alg = TrackFindingAlgorithm("trk_find_alg",
-        inputSourceLinks="BarrelTrackSourceLinks",
+        inputSourceLinks="BarrelSourceLinks",
         inputInitialTrackParameters= "InitTrackParams",#"InitTrackParamsFromClusters", 
         outputTrajectories="trajectories",
         OutputLevel=DEBUG)
@@ -129,7 +129,7 @@ parts_from_fit = ParticlesFromTrackFit("parts_from_fit",
         OutputLevel=DEBUG)
 
 trk_find_alg1 = TrackFindingAlgorithm("trk_find_alg1",
-        inputSourceLinks="BarrelTrackSourceLinks",
+        inputSourceLinks="BarrelSourceLinks",
         inputInitialTrackParameters= "InitTrackParamsFromClusters", 
         outputTrajectories="trajectories1",
         OutputLevel=DEBUG)

tracking/scripts/gen_central_electrons.cxx

@@ -39,7 +39,7 @@ void gen_central_electrons(int n_events = 100,
         FourVector(0.0, 0.0, 0.0, 0.938), 2212, 4);
 
     // Define momentum
-    Double_t p     = r1->Uniform(9.0, 10.0);
+    Double_t p     = r1->Uniform(1.0, 10.0);
     Double_t phi   = r1->Uniform(0.0, 2.0 * M_PI);
     Double_t costh = r1->Uniform(cos_theta_min, cos_theta_max);
     Double_t th    = std::acos(costh);