WIP: Crystal Calorimeter Simulation

• Added two variables: tracker_region_zmax and tracker_region_rmax into topside_defs.xml. ddsim/npsim is registering a plugin which needs those variables defined.
• Added bash script to run simulation.
• Added python script

options/example_crystal.py

+from Gaudi.Configuration import *
+  
+from GaudiKernel.DataObjectHandleBase import DataObjectHandleBase
+from Configurables import ApplicationMgr, EICDataSvc, PodioOutput, GeoSvc
+from GaudiKernel import SystemOfUnits as units
+
+geo_service  = GeoSvc("GeoSvc", detectors=["topside.xml"])
+podioevent   = EICDataSvc("EventDataSvc", inputs=["sim_input.root"], OutputLevel=DEBUG)
+
+from Configurables import PodioInput
+from Configurables import Jug__Digi__CrystalEndcapsDigi as CrystalEndcapsDigi
+from Configurables import Jug__Reco__CrystalEndcapsReco as CrystalEndcapsReco
+from Configurables import Jug__Reco__CalorimeterIslandCluster as IslandCluster
+from Configurables import Jug__Reco__ClusterRecoCoG as RecoCoG
+
+podioinput = PodioInput("PodioReader", collections=["mcparticles","CrystalEcalHits"], OutputLevel=DEBUG)
+emcaldigi = CrystalEndcapsDigi(inputHitCollection="CrystalEcalHits", outputHitCollection="RawDigiEcalHits")
+emcalreco = CrystalEndcapsReco(inputHitCollection="RawDigiEcalHits", outputHitCollection="RecoEcalHits",
+                               minModuleEdep=1.0*units.MeV)
+emcalcluster = IslandCluster(inputHitCollection="RecoEcalHits", outputClusterCollection="EcalClusters",
+                             minClusterCenterEdep=30*units.MeV, groupRange=2.0)
+clusterreco = RecoCoG(clusterCollection="EcalClusters", logWeightBase=4.2, moduleDimZName="CrystalBox_z_length")
+
+
+out = PodioOutput("out", filename="rec_output.root")
+out.outputCommands = ["keep EcalClusters"]
+
+ApplicationMgr(
+    TopAlg = [podioinput, emcaldigi, emcalreco, emcalcluster, clusterreco, out],
+    EvtSel = 'NONE',
+    EvtMax   = 10000,
+    ExtSvc = [podioevent],
+    OutputLevel=ERROR
+ )