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,
