from Gaudi.Configuration import * from Configurables import ApplicationMgr, FCCDataSvc, PodioOutput podioevent = FCCDataSvc("EventDataSvc", input="output_hcalSim_e50GeV_eta036_1events.root") # reads HepMC text file and write the HepMC::GenEvent to the data service from Configurables import PodioInput podioinput = PodioInput("PodioReader", collections=["HCalHits"], OutputLevel=DEBUG) from Configurables import GeoSvc geoservice = GeoSvc("GeoSvc", detectors=[ 'file:Detector/DetFCChhBaseline1/compact/FCChh_DectEmptyMaster.xml', 'file:Detector/DetFCChhHCalTile/compact/FCChh_HCalBarrel_TileCal.xml'], OutputLevel = INFO) from Configurables import CreateCaloCells createcells = CreateCaloCells("CreateCaloCells", doCellCalibration = False, addCellNoise = False, filterCellNoise = False, OutputLevel = DEBUG) createcells.hits.Path="HCalHits" createcells.cells.Path="HCalCells" out = PodioOutput("out", filename="output_HCalCells_digitisation_noNoise.root", OutputLevel = DEBUG) out.outputCommands = ["keep *"] #CPU information from Configurables import AuditorSvc, ChronoAuditor chra = ChronoAuditor() audsvc = AuditorSvc()
hcalExtBarrelReadoutName = "HCalExtBarrelReadout" hcalBarrelReadoutPhiEtaName = "BarHCal_Readout_phieta" hcalExtBarrelReadoutPhiEtaName = "ExtBarHCal_Readout_phieta" hcalEndcapReadoutName = "HECPhiEta" hcalFwdReadoutName = "HFwdPhiEta" # Number of events num_events = 3 from Gaudi.Configuration import * from Configurables import ApplicationMgr, FCCDataSvc, PodioOutput podioevent = FCCDataSvc("EventDataSvc", input="output_fullCalo_SimAndDigi_e50GeV_"+str(num_events)+"events.root") # reads HepMC text file and write the HepMC::GenEvent to the data service from Configurables import PodioInput podioinput = PodioInput("PodioReader", collections = [ecalBarrelCellsName, ecalEndcapCellsName, ecalFwdCellsName, hcalBarrelCellsName, hcalExtBarrelCellsName, hcalEndcapCellsName, hcalFwdCellsName]) from Configurables import GeoSvc geoservice = GeoSvc("GeoSvc", detectors=[ 'file:Detector/DetFCChhBaseline1/compact/FCChh_DectEmptyMaster.xml', 'file:Detector/DetFCChhTrackerTkLayout/compact/Tracker.xml', 'file:Detector/DetFCChhECalInclined/compact/FCChh_ECalBarrel_withCryostat.xml', 'file:Detector/DetFCChhHCalTile/compact/FCChh_HCalBarrel_TileCal.xml', 'file:Detector/DetFCChhHCalTile/compact/FCChh_HCalExtendedBarrel_TileCal.xml', 'file:Detector/DetFCChhCalDiscs/compact/Endcaps_coneCryo.xml', 'file:Detector/DetFCChhCalDiscs/compact/Forward_coneCryo.xml'], OutputLevel = WARNING) ecalBarrelNoisePath = "/eos/project/f/fccsw-web/testsamples/elecNoise_ecalBarrel_50Ohm_traces2_2shieldWidth.root" ecalEndcapNoisePath = "/eos/project/f/fccsw-web/testsamples/elecNoise_emec_50Ohm_2shieldWidth_6layers.root"
from Gaudi.Configuration import * from Configurables import ApplicationMgr, FCCDataSvc, PodioOutput podioevent = FCCDataSvc("EventDataSvc", input="output_ecalSim_e50GeV_eta0_10events.root") # reads HepMC text file and write the HepMC::GenEvent to the data service from Configurables import PodioInput podioinput = PodioInput("PodioReader", collections=["ECalHits", "ECalPositionedHits"], OutputLevel=DEBUG) from Configurables import GeoSvc geoservice = GeoSvc( "GeoSvc", detectors=[ 'file:Detector/DetFCChhBaseline1/compact/FCChh_DectEmptyMaster.xml', 'file:Detector/DetFCChhECalSimple/compact/FCChh_ECalBarrel_Mockup.xml' ], OutputLevel=INFO) # common ECAL specific information # readout name ecalReadoutName = "ECalHitsPhiEta" # active material identifier name ecalIdentifierName = "active_layer" # active material volume name ecalVolumeName = "LAr_sensitive" # ECAL bitfield names & values ecalFieldNames = ["system", "ECAL_Cryo", "bath", "EM_barrel"]
import os from Gaudi.Configuration import * import GaudiKernel.SystemOfUnits as units from Configurables import ApplicationMgr, FCCDataSvc, PodioOutput podioevent = FCCDataSvc( "EventDataSvc", input= "root://eospublic.cern.ch//eos/experiment/fcc/ee/tutorial/fccee_idea_pgun.root" ) from Configurables import PodioInput, ReadTestConsumer podioinput = PodioInput("PodioReader", collections=["positionedHits_DCH"], OutputLevel=DEBUG) # Parses the given xml file from Configurables import GeoSvc geoservice = GeoSvc("GeoSvc", detectors=[ os.environ.get("FCC_DETECTORS", "") + '/Detector/DetFCCeeIDEA/compact/FCCee_DectMaster.xml', ]) from Configurables import CreateDCHHits createhits = CreateDCHHits("CreateDCHHits", readoutName="DriftChamberCollection", EdepCut=100 * 1e-9, DCACut=0.8, OutputLevel=INFO)
from Gaudi.Configuration import * from Configurables import ApplicationMgr, FCCDataSvc, PodioOutput podioevent = FCCDataSvc("EventDataSvc") # reads HepMC text file and write the HepMC::GenEvent to the data service from Configurables import PodioInput, ReadTestConsumer podioinput = PodioInput( "PodioReader", filename="output.root", collections=["allGenVertices", "allGenParticles", "hits"], OutputLevel=DEBUG) checker = ReadTestConsumer() out = PodioOutput("out", filename="out2.root", OutputLevel=DEBUG) out.outputCommands = ["keep *"] ApplicationMgr( TopAlg=[podioinput, checker, out], EvtSel='NONE', EvtMax=3, ExtSvc=[podioevent], )
import argparse recoparser = argparse.ArgumentParser() recoparser.add_argument('-r','--rootfile', type=str, default=None, help='rootfile') recoargs, _ = recoparser.parse_known_args() rootfile = recoargs.rootfile print "==================================" print "rootfile = ", rootfile print "==================================" ## read in generated particles from ROOT via podio podioevent = FCCDataSvc("EventDataSvc", input=rootfile) #"Sim/SimPapas/data/ee_ZH_Zmumu_Hbb.root") from Configurables import PodioInput, ReadTestConsumer podioinput = PodioInput("PodioReader", collections=["GenVertex", "GenParticle"], OutputLevel=DEBUG) from CMS_detector_cfg import detservice from papas_cfg import papasalg_pdebug #output fcc particles to root from Configurables import PodioOutput out = PodioOutput("out", OutputLevel=INFO) out.outputCommands = ["keep *"] from Configurables import ApplicationMgr ApplicationMgr( ## all algorithms should be put here TopAlg=[podioinput, papasalg_pdebug, out], EvtSel='NONE',
from Gaudi.Configuration import * from Configurables import ApplicationMgr, FCCDataSvc, PodioOutput # v01 production - min. bias events podioevent = FCCDataSvc( "EventDataSvc", input="/eos/project/f/fccsw-web/testsamples/caloCells_minBiasEvent.root") # reads HepMC text file and write the HepMC::GenEvent to the data service from Configurables import PodioInput podioinput = PodioInput("in", collections=["ECalBarrelCells"]) from Configurables import GeoSvc geoservice = GeoSvc( "GeoSvc", detectors=[ 'file:Detector/DetFCChhBaseline1/compact/FCChh_DectEmptyMaster.xml', 'file:Detector/DetFCChhECalInclined/compact/FCChh_ECalBarrel_withCryostat.xml' ], OutputLevel=WARNING) # Pileup in ECal Barrel # readout name ecalBarrelReadoutNamePhiEta = "ECalBarrelPhiEta" # geometry tool from Configurables import TubeLayerPhiEtaCaloTool ecalBarrelGeometry = TubeLayerPhiEtaCaloTool( "EcalBarrelGeo",
if bFieldOff: pileupNoisePath = "/afs/cern.ch/work/c/cneubuse/public/FCChh/noBfield/noiseBarrel_mu" + str( puEvents) + ".root" inputPileupNoisePerCell = "/afs/cern.ch/work/c/cneubuse/public/FCChh/cellNoise_map_segHcal_noiseLevelElectronicsPileup_mu" + str( puEvents) + ".root" ############################################################################################################## ####### INPUT ############# ############################################################################################################## # reads HepMC text file and write the HepMC::GenEvent to the data service from Configurables import ApplicationMgr, FCCDataSvc, PodioInput, PodioOutput podioevent = FCCDataSvc("EventDataSvc", input=input_name) podioinput = PodioInput("PodioReader", collections=inputCollections, OutputLevel=DEBUG) print "Reading collections: ", inputCollections ############################################################################################################## ####### CELL POSITIONS TOOLS ############# ############################################################################################################## #Configure tools for calo cell positions from Configurables import CellPositionsECalBarrelTool, CellPositionsHCalBarrelNoSegTool, CellPositionsHCalBarrelTool, CellPositionsCaloDiscsTool, CellPositionsTailCatcherTool ECalBcells = CellPositionsECalBarrelTool("CellPositionsECalBarrel", readoutName=ecalBarrelReadoutName, OutputLevel=INFO) EMECcells = CellPositionsCaloDiscsTool("CellPositionsEMEC", readoutName=ecalEndcapReadoutName,