Fixing empty plots

modified: emcal_electrons.sh modified: options/crystal_calorimeter_reco.py

ecal/options/crystal_calorimeter_reco.py

@@ -10,9 +10,22 @@ if "JUGGLER_DETECTOR" in os.environ :
   detector_name = str(os.environ["JUGGLER_DETECTOR"])
 
 # todo add checks
-input_sim_file  = str(os.environ["JUGGLER_SIM_FILE"])
-output_rec_file = str(os.environ["JUGGLER_REC_FILE"])
-n_events = str(os.environ["JUGGLER_N_EVENTS"])
+input_sim_file = "jug_input.root"
+if "JUGGLER_SIM_FILE" in os.environ :
+  input_sim_file  = str(os.environ["JUGGLER_SIM_FILE"])
+else :
+    print(" ERROR : JUGGLER_SIM_FILE not set" )
+
+output_rec_file = "jug_rec.root"
+if "JUGGLER_REC_FILE" in os.environ :
+  output_rec_file = str(os.environ["JUGGLER_REC_FILE"])
+else :
+    print(" ERROR : JUGGLER_REC_FILE not set" )
+
+
+n_events = 100
+if "JUGGLER_N_EVENTS" in os.environ :
+  n_events = str(os.environ["JUGGLER_N_EVENTS"])
 
 geo_service  = GeoSvc("GeoSvc", detectors=["{}.xml".format(detector_name)])
 podioevent   = EICDataSvc("EventDataSvc", inputs=[input_sim_file], OutputLevel=DEBUG)
@@ -31,15 +44,28 @@ from Configurables import Jug__Reco__ClusterRecoCoG as RecoCoG
 podioinput = PodioInput("PodioReader", collections=["mcparticles","CrystalEcalHits"], OutputLevel=DEBUG)
 
 ## copiers to get around input --> output copy bug. Note the "2" appended to the output collection.
-copier = MCCopier("MCCopier", inputCollection="mcparticles", outputCollection="mcparticles2",OutputLevel=DEBUG) 
+copier    = MCCopier("MCCopier", inputCollection="mcparticles", outputCollection="mcparticles2",OutputLevel=DEBUG) 
 calcopier = CalCopier("CalCopier", inputCollection="CrystalEcalHits", outputCollection="CrystalEcalHits2",OutputLevel=DEBUG) 
 
-emcaldigi = CrystalEndcapsDigi("ecal_digi", inputHitCollection="CrystalEcalHits", outputHitCollection="RawDigiEcalHits")
-emcalreco = CrystalEndcapsReco("ecal_reco", inputHitCollection="RawDigiEcalHits", outputHitCollection="RecoEcalHits",
-                               minModuleEdep=1.0*units.MeV)
-emcalcluster = IslandCluster("emcal_cluster", inputHitCollection="RecoEcalHits", outputClusterCollection="EcalClusters",
-                             minClusterCenterEdep=30*units.MeV, groupRange=2.0)
-clusterreco = RecoCoG("cluster_reco", clusterCollection="EcalClusters", logWeightBase=4.2, moduleDimZName="CrystalBox_z_length")
+emcaldigi = CrystalEndcapsDigi("ecal_digi", 
+        inputHitCollection="CrystalEcalHits", 
+        outputHitCollection="RawDigiEcalHits")
+
+emcalreco = CrystalEndcapsReco("ecal_reco", 
+        inputHitCollection="RawDigiEcalHits", 
+        outputHitCollection="RecoEcalHits",
+        minModuleEdep=0.00001*units.MeV)
+
+emcalcluster = IslandCluster("emcal_cluster", 
+        inputHitCollection="RecoEcalHits", 
+        outputClusterCollection="EcalClusters",
+        minClusterCenterEdep=0.00001*units.MeV, 
+        groupRange=2.0)
+
+clusterreco = RecoCoG("cluster_reco", 
+        clusterCollection="EcalClusters", 
+        logWeightBase=4.2, 
+        moduleDimZName="CrystalBox_z_length")
 
 out = PodioOutput("out", filename=output_rec_file)