#services that need to be reinit ServiceMgr.MultipleEventLoopMgr.ToBeReinitialized = ["StoreGateSvc", "DetectorStore"] # # I/O stuff # ServiceMgr.EventSelector.InputCollections = [ "root://eosatlas//eos/atlas/atlascerngroupdisk/proj-sit/digitization/RTT/mc10/mc10_7TeV.005008.CavernInput.merge.HITS.e4_e607_s951_s952_tid170551_00/HITS.170551._000011.pool.root.1" ] # Set up for seeking. ServiceMgr.EventSelector.CollectionType = "SeekableROOT" include( "AthenaPoolCnvSvc/WriteAthenaPool_jobOptions.py" ) from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream( "Stream1" ) Stream1.OutputFile = "MultiplePassFile.root"; # ** mandatory parameter ** // The output file name Stream1.ForceRead=TRUE #force read of output data objs Stream1.ItemList=["EventInfo#*", "McEventCollection#*"] #-------------------------------------------------------------- #--- Monitoring and Debug printouts #-------------------------------------------------------------- MessageSvc.OutputLevel = INFO Stream1.OutputLevel = DEBUG #StoreGateSvc.Dump=1 #StoreGateSvc.OutputLevel=DEBUG from AthenaCommon.AppMgr import theAuditorSvc from GaudiAud.GaudiAudConf import ChronoAuditor theAuditorSvc += ChronoAuditor()
if numThreads >= 2: InDetSiSPSeededTrackFinder.Cardinality = numThreads topSequence += InDetSiSPSeededTrackFinder # Print algorithms if doPrint: printfunc(topSequence) # Set the number of events to be processed theApp.EvtMax = EvtMax # Output file if doDump: from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream outStream = AthenaPoolOutputStream( "OutStream", "SiSPSeededTracksStandaloneFromESD.pool.root") outStream.ItemList = [ "xAOD::EventInfo#EventInfo", "xAOD::EventAuxInfo#EventInfoAux." ] outStream.ItemList += ["TrackCollection#" + TracksLocation] #-------------------------------------------------------------- # Set output lvl (VERBOSE, DEBUG, INFO, WARNING, ERROR, FATAL) #-------------------------------------------------------------- ServiceMgr.MessageSvc.OutputLevel = INFO ServiceMgr.MessageSvc.Format = "% F%50W%S%7W%R%T %0W%M" if numThreads >= 2: from SCT_ConditionsAlgorithms.SCTCondAlgCardinality import sctCondAlgCardinality sctCondAlgCardinality.set(numThreads)
ServiceMgr.MultipleEventLoopMgr.ToBeReinitialized = [ "StoreGateSvc", "DetectorStore" ] # # I/O stuff # ServiceMgr.EventSelector.InputCollections = [ "root://eosatlas.cern.ch//eos/atlas/atlascerngroupdisk/proj-sit/digitization/RTT/mc10/mc10_7TeV.005008.CavernInput.merge.HITS.e4_e607_s951_s952_tid170551_00/HITS.170551._000011.pool.root.1" ] include("AthenaPoolCnvSvc/WriteAthenaPool_jobOptions.py") from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("Stream1") Stream1.OutputFile = "MultiplePassFile.root" # ** mandatory parameter ** // The output file name Stream1.ItemList = ["EventInfo#*", "McEventCollection#*"] #-------------------------------------------------------------- #--- Monitoring and Debug printouts #-------------------------------------------------------------- MessageSvc.OutputLevel = INFO Stream1.OutputLevel = DEBUG #StoreGateSvc.Dump=1 #StoreGateSvc.OutputLevel=DEBUG from AthenaCommon.AppMgr import theAuditorSvc from GaudiAud.GaudiAudConf import ChronoAuditor theAuditorSvc += ChronoAuditor()
"root://eosatlas//eos/atlas/atlascerngroupdisk/proj-sit/digitization/RTT/mc10/mc10_7TeV.105001.pythia_minbias.merge.HITS.e577_s932_s952_tid170554_00/HITS.170554._000034.pool.root.1", "root://eosatlas//eos/atlas/atlascerngroupdisk/proj-sit/digitization/RTT/mc10/mc10_7TeV.105001.pythia_minbias.merge.HITS.e577_s932_s952_tid170554_00/HITS.170554._000043.pool.root.1" ] ServiceMgr.EventMixer.OutputRunNumber=54321 ServiceMgr.EventMixer.EventNumbers=[1, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 47] #-------------------------------------------------------------- #--- Output stream configuration #-------------------------------------------------------------- include( "AthenaPoolCnvSvc/WriteAthenaPool_jobOptions.py" ) from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream( "Stream1" ) Stream1.OutputFile = "MixedFile.root"; # ** mandatory parameter ** // The output file name Stream1.ForceRead=TRUE; #force read of output data objs Stream1.ItemList=["MergedEventInfo#*", "EventInfo#*", "McEventCollection#*", "SiHitCollection#SCT_Hits", "SiHitCollection/PixelHits"]; #SYNTAX ERROR: / not # #-------------------------------------------------------------- #--- Monitoring and Debug printouts #-------------------------------------------------------------- MessageSvc.OutputLevel = INFO ServiceMgr.EventMixer.OutputLevel = VERBOSE MessageSvc.setVerbose += [ "MixingEventSelector::TriggerList" ] Stream1.OutputLevel = DEBUG #StoreGateSvc.Dump=1 #StoreGateSvc.OutputLevel=DEBUG from AthenaCommon.AppMgr import theAuditorSvc
def hits_persistency(): """ HITS POOL file persistency """ from G4AtlasApps.SimFlags import simFlags from AthenaCommon.DetFlags import DetFlags from AthenaCommon.AthenaCommonFlags import athenaCommonFlags from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream ## Not yet understood, but need to treat StreamHITS as alg in Hive. ## Seems to also work fine outside of Hive, but to be extra safe I'm ## only changing it in Hive. from AthenaCommon.ConcurrencyFlags import jobproperties as concurrencyProps if concurrencyProps.ConcurrencyFlags.NumThreads() > 0: as_alg = True else: as_alg = False ## NB. Two-arg constructor is needed, since otherwise metadata writing fails! stream1 = AthenaPoolOutputStream("StreamHITS", athenaCommonFlags.PoolHitsOutput(), asAlg=as_alg) ## Write geometry tag info - move to main method #import EventInfoMgt.EventInfoMgtInit ## EventInfo & TruthEvent always written by default stream1.ForceRead = True stream1.ItemList = [ "EventInfo#*", "McEventCollection#TruthEvent", "JetCollection#*" ] ## If we are running quasi-stable particle simulation, include the original event record if hasattr( simFlags, 'IncludeParentsInG4Event' ) and simFlags.IncludeParentsInG4Event.statusOn and simFlags.IncludeParentsInG4Event( ): stream1.ItemList += ["McEventCollection#GEN_EVENT"] stream1.ItemList += ["xAOD::JetContainer#*", "xAOD::JetAuxContainer#*"] ## Make stream aware of aborted events stream1.AcceptAlgs = ["G4AtlasAlg"] ## Detectors ## Inner Detector if DetFlags.ID_on(): stream1.ItemList += [ "SiHitCollection#*", "TRTUncompressedHitCollection#*", "TrackRecordCollection#CaloEntryLayer" ] ## Calo if DetFlags.Calo_on(): stream1.ItemList += [ "CaloCalibrationHitContainer#*", "LArHitContainer#*", "TileHitVector#*", #"SimpleScintillatorHitCollection#*", "TrackRecordCollection#MuonEntryLayer" ] ## Muon if DetFlags.Muon_on(): stream1.ItemList += [ "RPCSimHitCollection#*", "TGCSimHitCollection#*", "CSCSimHitCollection#*", "MDTSimHitCollection#*", "TrackRecordCollection#MuonExitLayer" ] if hasattr(simFlags, 'SimulateNewSmallWheel'): if simFlags.SimulateNewSmallWheel(): stream1.ItemList += ["GenericMuonSimHitCollection#*"] ## Lucid if DetFlags.Lucid_on(): stream1.ItemList += ["LUCID_SimHitCollection#*"] ## FwdRegion if DetFlags.FwdRegion_on(): stream1.ItemList += ["SimulationHitCollection#*"] ## ZDC if DetFlags.ZDC_on(): stream1.ItemList += [ "ZDC_SimPixelHit_Collection#*", "ZDC_SimStripHit_Collection#*" ] ## ALFA if DetFlags.ALFA_on(): stream1.ItemList += [ "ALFA_HitCollection#*", "ALFA_ODHitCollection#*" ] ## AFP if DetFlags.AFP_on(): stream1.ItemList += [ "AFP_TDSimHitCollection#*", "AFP_SIDSimHitCollection#*" ] ### Ancillary scintillators #stream1.ItemList += ["ScintillatorHitCollection#*"] ## TimingAlg stream1.ItemList += ["RecoTimingObj#EVNTtoHITS_timings"] ## Add cosmics and test beam configuration hit persistency if required cf. geom tag layout = simFlags.SimLayout.get_Value() if "tb" not in layout: from AthenaCommon.BeamFlags import jobproperties if jobproperties.Beam.beamType() == 'cosmics' or \ (hasattr(simFlags, "WriteTR") and simFlags.WriteTR.statusOn) or \ (hasattr(simFlags, "ReadTR") and simFlags.ReadTR.statusOn): stream1.ItemList += [ "TrackRecordCollection#CosmicRecord", "TrackRecordCollection#CosmicPerigee" ] else: ## CTB-specific if layout.startswith("ctb"): if simFlags.LArFarUpstreamMaterial.statusOn and simFlags.LArFarUpstreamMaterial.get_Value( ): stream1.ItemList.append( "TrackRecordCollection#LArFarUpstreamMaterialExitLayer" ) ## Persistency of test-beam layout if layout.startswith('ctb') or layout.startswith('tb_Tile2000_'): stream1.ItemList += ["TBElementContainer#*"]
#-------------------------------------------------------------- # EventLoop #-------------------------------------------------------------- from AthenaCommon.ConcurrencyFlags import jobproperties as jp nThreads = jp.ConcurrencyFlags.NumThreads() if nThreads > 0: EventLoop = Service("AthenaHiveEventLoopMgr") else: EventLoop = Service("AthenaEventLoopMgr") EventLoop.UseSecondaryEventNumber = True svcMgr += EventLoop #-------------------------------------------------------------- # Output options #-------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream = AthenaPoolOutputStream( "Stream" ) Stream.OutputFile = "OutputRDO.root" Stream.ItemList = [ 'xAOD::EventInfo#EventInfo', 'xAOD::EventAuxInfo#EventInfoAux.' ] #-------------------------------------------------------------- # Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL ) #-------------------------------------------------------------- svcMgr.MessageSvc = Service( "MessageSvc" ) svcMgr.MessageSvc.OutputLevel = WARNING svcMgr.MessageSvc.debugLimit = 100000 # No stats printout include( "AthenaPoolTest/NoStats_jobOptions.py" )
# Copyright (C) 2002-2019 CERN for the benefit of the ATLAS collaboration include.block('EventOverlayJobTransforms/SignalOutputItemList_jobOptions.py') from AthenaCommon.DetFlags import DetFlags from OverlayCommonAlgs.OverlayFlags import overlayFlags # TODO: some event info update will be needed # The output - signal from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream signalStream = AthenaPoolOutputStream('StreamRDO_SGNL', overlayFlags.PoolSignalRDOOutput(), asAlg=True) signalStream.ItemList = [] # Event info signalStream.ItemList += [ 'xAOD::EventInfo#' + overlayFlags.sigPrefix() + 'EventInfo', 'xAOD::EventAuxInfo#' + overlayFlags.sigPrefix() + 'EventInfoAux.' ] signalStream.ItemList += [ 'RecoTimingObj#' + overlayFlags.sigPrefix() + 'EVNTtoHITS_timings' ] if DetFlags.overlay.Truth_on(): signalStream.ItemList += [ 'McEventCollection#' + overlayFlags.sigPrefix() + 'TruthEvent' ]
from Starlight_i.Starlight_iConf import Starlight_i genSeq += Starlight_i("Starlight") genSeq.Starlight.McEventKey = "GEN_EVENT" evgenConfig.generators += ["Starlight"] ## Extra stream persistency ## 2101 == EventInfo, 133273 == MCTruth (HepMC) from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream _evgenstream = AthenaPoolOutputStream("StreamEVGEN") _evgenstream.ItemList = ["2101#*", "133273#GEN_EVENT"] del _evgenstream
raise RuntimeError("No outputHitsFile provided.") from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream try: Stream1 = AthenaPoolOutputStream( "StreamHITS", athenaCommonFlags.PoolHitsOutput() ) except: Stream1 = AthenaPoolOutputStream( "StreamHITS", "DidNotSetOutputName.root" ) # The next line is an example on how to exclude clid's if they are causing a problem #Stream1.ExcludeList = ['6421#*'] #-------------------------------------------------------------- # Specify collections for output HIT files, as not all are required. #-------------------------------------------------------------- #Truth Stream1.ItemList=["EventInfo#*", "McEventCollection#TruthEvent", # mc truth (hepmc) "JetCollection#*", # Truth jets reconstructed (optionally) during evgen "TrackRecordCollection#MuonEntryLayer"] # others not used in pileup # "TrackRecordCollection#MuonExitLayer", # not used in pileup # "TrackRecordCollection#CaloEntryLayer"] # not used in pileup # Deal with "new" truth jet collections properly from PyJobTransforms.trfUtils import releaseIsOlderThan if releaseIsOlderThan(20,0): #Hack to maintain compatibility of G4AtlasApps trunk with #19.2.X.Y after EDM changes in release 20.0.0. Stream1.ItemList += ["xAOD::JetContainer_v1#*", "xAOD::JetAuxContainer_v1#*"] else: Stream1.ItemList += ["xAOD::JetContainer#*", "xAOD::JetAuxContainer#*"]
svcMgr.ProxyProviderSvc.OutputLevel = DEBUG svcMgr.AthenaPoolAddressProviderSvcPrimary.OutputLevel = DEBUG svcMgr.AthenaPoolAddressProviderSvcSecondary.OutputLevel = DEBUG svcMgr.DoubleEventSelector.OutputLevel = DEBUG #-------------------------------------------------------------- # Output options #-------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream # TODO: noTag=True to avoid warning, needs EventInfo overlay implemented Stream1 = AthenaPoolOutputStream("Stream1", asAlg=True, noTag=True) Stream1.OutputLevel = INFO Stream1.OutputFile = "OutputRDO.root" # List of DO's to write out Stream1.ItemList = [] Stream1.ItemList += ["McEventCollection#TruthEvent"] Stream1.ItemList += ["RecoTimingObj#EVNTtoHITS_timings"] #-------------------------------------------------------------- # Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL ) #-------------------------------------------------------------- svcMgr.MessageSvc.debugLimit = 100000 # No stats printout include("AthenaPoolTest/NoStats_jobOptions.py") #============================================================== # # End of job options file #
fullItemList += [ "TBScintillatorRawCont#*" ] fullItemList += [ "TBTDC#*" ] fullItemList += [ "TBTDCRawCont#*" ] fullItemList += [ "TBTrack#*" ] fullItemList += [ "TBTriggerPatternUnit#*" ] fullItemList += [ "TBADCRawCont#*" ] if doTile: fullItemList += [ "TileBeamElemContainer#*" ] #adding digits for the moment fullItemList += [ "LArDigitContainer#*" ] fullItemList += [ "TileDigitsContainer#*" ] # if not doSim: Stream1.ItemList = fullItemList # else: # ToolSvc.Stream1.ItemList = fullItemList #----- # Tricks for conditions #----- #conddb.iovdbsvc.dbConnection="impl=cool;techno=sqlite;schema=sqlite200/ALLP200.db;X:TBDP200" #conddb.iovdbsvc.dbConnection="" conddb.addFolder("LAR","/LAR/BadChannels/BadChannels<tag>LARBadChannelsBadChannels-MC-empty</tag>") if not doSim: conddb.addFolder("LAR","/LAR/BadChannels/MissingFEBs<tag>LARBadChannelsMissingFEBs-empty</tag>") else: conddb.iovdbsvc.Folders += ["<db>impl=cool;techno=logical;schema=COOLONL_LAR;X:TBDP200</db> /LAR/BadChannels/MissingFEBs<tag>LARBadChannelsMissingFEBs-empty</tag>"] #conddb.addFolder("CALO","/CALO/Identifier/CaloTTOnOffIdMapAtlas")
#-------------------------------------------------------------- from Digitization.DigitizationWriteMetaData import writeDigitizationMetadata writeDigitizationMetadata() #-------------------------------------------------------------- # Pool Output (Change this to use a different file) #-------------------------------------------------------------- if DetFlags.writeRDOPool.any_on(): from AthenaCommon.AthenaCommonFlags import athenaCommonFlags from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream streamRDO = AthenaPoolOutputStream( "StreamRDO", athenaCommonFlags.PoolRDOOutput.get_Value(), asAlg=True) streamRDO.ForceRead = True from Digitization.DigiOutput import getStreamRDO_ItemList streamRDO.ItemList = getStreamRDO_ItemList(logDigitization_flags) streamRDO.AcceptAlgs += [digitizationFlags.digiSteeringConf.get_Value()] streamRDO.OutputFile = athenaCommonFlags.PoolRDOOutput() if athenaCommonFlags.UseLZMA(): ServiceMgr.AthenaPoolCnvSvc.PoolAttributes += [ "DatabaseName = '" + athenaCommonFlags.PoolRDOOutput() + "'; COMPRESSION_ALGORITHM = '2'" ] ServiceMgr.AthenaPoolCnvSvc.PoolAttributes += [ "DatabaseName = '" + athenaCommonFlags.PoolRDOOutput() + "'; COMPRESSION_LEVEL = '1'" ] else: ServiceMgr.AthenaPoolCnvSvc.PoolAttributes += [ "DatabaseName = '" + athenaCommonFlags.PoolRDOOutput() + "'; COMPRESSION_ALGORITHM = '1'"
] ServiceMgr.EventMixer.OutputRunNumber = 54321 ServiceMgr.EventMixer.EventNumbers = [ 1, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 47 ] #-------------------------------------------------------------- #--- Output stream configuration #-------------------------------------------------------------- include("AthenaPoolCnvSvc/WriteAthenaPool_jobOptions.py") from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("Stream1") Stream1.OutputFile = "MixedFile.root" # ** mandatory parameter ** // The output file name Stream1.ItemList = [ "MergedEventInfo#*", "EventInfo#*", "McEventCollection#*", "SiHitCollection#SCT_Hits", "SiHitCollection/PixelHits" ] #SYNTAX ERROR: / not # #-------------------------------------------------------------- #--- Monitoring and Debug printouts #-------------------------------------------------------------- MessageSvc.OutputLevel = INFO ServiceMgr.EventMixer.OutputLevel = VERBOSE MessageSvc.setVerbose += ["MixingEventSelector::TriggerList"] Stream1.OutputLevel = DEBUG #StoreGateSvc.Dump=1 #StoreGateSvc.OutputLevel=DEBUG from AthenaCommon.AppMgr import theAuditorSvc
topSequence.ttbarFilter.jetsName = 'Kt4LCTopoJets' # use new LC topo jets #topSequence.ttbarFilter.maxDeltaRcx = 100.4 # effectively switch off # cluster thinning print topSequence #============================================================== # # End of job options file # ############################################################### #-------------------------------------------------------------- #--- Secondary Write portion ----- Don't change it !!! #-------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream StreamDPD = AthenaPoolOutputStream("StreamDPD") StreamDPD.ItemList = ['EventInfo#*', 'TrackRecordCollection#*'] StreamDPD.ItemList += ['ElectronContainer#ElectronAODCollection'] StreamDPD.ItemList += ['egDetailContainer#egDetailAOD'] StreamDPD.ItemList += ['MissingET#MET_RefFinal'] StreamDPD.ItemList += ['MissingET#MET_LocHadObj'] StreamDPD.ItemList += ['CaloClusterContainer#CaloCalTopoCluster'] StreamDPD.ItemList += ['Analysis::MuonContainer#StacoMuonCollection'] StreamDPD.ItemList += ['JetCollection#Kt4LCTopoJets'] StreamDPD.ItemList += ['JetCollection#Kt4TruthJets'] #StreamDPD.ItemList += ['JetCollection#Cone4H1TowerJets'] StreamDPD.ItemList += ['JetKeyDescriptor#JetKeyMap'] StreamDPD.ItemList += ['Rec::TrackParticleContainer#TrackParticleCandidate'] StreamDPD.ItemList += ['Rec::TrackParticleContainer#StacoTrackParticles'] StreamDPD.ItemList += ['TruthParticleContainer#SpclMC'] StreamDPD.ItemList += ['McEventCollection#GEN_AOD']
######## # Read the McEventCollection # if not 'DUMP' in dir(): DUMP = True pass if DUMP: job += CfgMgr.DumpMC("ReadGenEvent", McEventKey=MCEVENTKEY) #--------------------------------------------------------------- # Pool Persistency #--------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream outStream = AthenaPoolOutputStream("OutStream") outStream.ItemList = [ "EventInfo#McEventInfo", "McEventCollection#" + MCEVENTKEY ] if not "OUTPUT" in dir(): OUTPUT = "mc.event.pool" outStream.OutputFile = OUTPUT #-------------------------------------------------------------- # Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL ) #-------------------------------------------------------------- #svcMgr.MessageSvc.OutputLevel = DEBUG #============================================================== # # End of job options file #
theApp.CreateSvc += [svcMgr.THistSvc.getFullName()] if not 'TUPLEFILENAME' in dir(): TUPLEFILENAME = 'mcaod.ttbar.root' pass svcMgr.THistSvc.Output += [ "mcaod DATAFILE='%s' OPT='RECREATE'" % TUPLEFILENAME ] #-------------------------------------------------------------- # Output options #-------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream outStream = AthenaPoolOutputStream("StreamEvGen") outStream.ItemList = [ "EventInfo#*", "McEventCollection#MyGEN_EVENT", "TruthParticleContainer#*", ] OUTPUT = os.path.join(os.path.dirname(INPUT[0]), "filtered.%s" % os.path.basename(INPUT[0])) import os outStream.OutputFile = OUTPUT #-------------------------------------------------------------- # Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL ) #-------------------------------------------------------------- svcMgr.MessageSvc.defaultLimit = 4000000 #svcMgr.MessageSvc.OutputLevel = ERROR #==============================================================
from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream # Stream1 Stream1 = AthenaPoolOutputStream( "Stream1", noTag=True ) Stream1.OutputFile = "SimplePoolFile1.root" # List of DO's to write out Stream1.ItemList += fullItemList # Stream2 Stream2 = AthenaPoolOutputStream( "Stream2", noTag=True ) Stream2.OutputFile = "SimplePoolFile2.root" # List of DO's to write out Stream2.ItemList += fullItemList # Remove dumm_A from Stream1 Stream1.ItemList = [ o for o in Stream1.ItemList if not o in 'dummy_A#*' ] Stream1.ItemList = [ o for o in Stream1.ItemList if not o in 'IAthenaPoolTestCollection#*' ] # Remove dumm_E from Stream2 Stream2.ItemList = [ o for o in Stream2.ItemList if not o in 'dummy_E#*' ] Stream2.ItemList = [ o for o in Stream2.ItemList if not o in 'AthenaPoolTestMatrix#*' ] printfunc ("Stream1.ItemList ",Stream1.ItemList) printfunc ("Stream2.ItemList ",Stream2.ItemList) #-------------------------------------------------------------- # Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL ) #-------------------------------------------------------------- svcMgr.MessageSvc.OutputLevel = WARNING svcMgr.MessageSvc.debugLimit = 100000
#Explicitly specify the output file catalog from PoolSvc.PoolSvcConf import PoolSvc svcMgr += PoolSvc() svcMgr.PoolSvc.WriteCatalog = "xmlcatalog_file:SplittableData.xml" from AthenaPoolCnvSvc.AthenaPoolCnvSvcConf import AthenaPoolCnvSvc svcMgr += AthenaPoolCnvSvc() svcMgr.AthenaPoolCnvSvc.CommitInterval = 10 include("AthenaPoolMultiTest/ExampleStreamConfig.py") from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream DataStream = AthenaPoolOutputStream("DataStream", "AthenaPoolMultiTest_Splittable0.root", True) DataStream.ItemList = exampleItemList # comes from ExampleStreamConfig DataStream.ItemList += ["FauxTriggerMap#*"] # add item not in StreamConfig DataStream.MetadataItemList += exampleMetadataList DataStream.AcceptAlgs = ["PassAllFilter"] DataStream.RequireAlgs = ["PassAllFilter"] DataStream.VetoAlgs = ["PassNoneFilter"] #-------------------------------------------------------------- # Consequently, should have 3 output files # Numbers 3, 4, from PassAll and 1 from Passnone #-------------------------------------------------------------- # # Now configure output collection # #-------------------------------------------------------------- from RegistrationServices.RegistrationServicesConf import RegistrationStream
filAlg=ttbarFilterAlgorithm("ttbarFilter") # to change properties of the filter instantiate it like this: # filAlg=ttbarFilterAlgorithm("ttbarFilter",JetMinEt=40*GeV) topSequence += ConfigurablePyAlgorithm("ttbarFilter") #============================================================== # # End of job options file # ############################################################### #-------------------------------------------------------------- #--- Secondary Write portion ----- Don't change it !!! #-------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream StreamDPD = AthenaPoolOutputStream("StreamDPD") StreamDPD.ItemList = ['EventInfo#*', 'TrackRecordCollection#*'] StreamDPD.ItemList += ['ElectronContainer#ElectronAODCollection'] StreamDPD.ItemList += ['egDetailContainer#egDetailAOD'] StreamDPD.ItemList += ['MissingET#MET_RefFinal'] StreamDPD.ItemList += ['MissingET#MET_LocHadObj'] StreamDPD.ItemList += ['CaloClusterContainer#CaloCalTopoCluster'] StreamDPD.ItemList += ['Analysis::MuonContainer#StacoMuonCollection'] StreamDPD.ItemList += ['JetCollection#Kt4LCTopoJets'] StreamDPD.ItemList += ['JetKeyDescriptor#JetKeyMap'] StreamDPD.ItemList += ['JetCollection#Kt4TruthJets'] StreamDPD.ItemList += ['Rec::TrackParticleContainer#TrackParticleCandidate'] StreamDPD.ItemList += ['Rec::TrackParticleContainer#StacoTrackParticles'] StreamDPD.ItemList += ['TruthParticleContainer#SpclMC'] StreamDPD.ItemList += ['McEventCollection#GEN_AOD'] StreamDPD.OutputFile = "SkimmedThin.AOD.pool.root"
theApp.EvtMax = EVTMAX #-------------------------------------------------------------- # Private Application Configuration options #-------------------------------------------------------------- ## schedule a noop algorithm: it just print where it is from PerfMonTests.PerfMonTestsConf import PerfMonTest__NoopAlg topSequence += PerfMonTest__NoopAlg("NoopAlg") #--------------------------------------------------------------- # Pool Persistency #--------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream outStream = AthenaPoolOutputStream("OutStream") outStream.ItemList = ["EventInfo#McEventInfo"] if not 'OUTPUT' in dir(): OUTPUT = "my.data.pool" pass outStream.OutputFile = OUTPUT ## somehow better configure the AthenaPoolCnvSvc for our small ## persistency output job svcMgr.AthenaPoolCnvSvc.CommitInterval = 1000 svcMgr.MessageSvc.OutputLevel = ERROR #============================================================== # # End of job options file #
esdList += CfgItemList("TrigEsd", items=["JetCollection#HLT_AutoKey*"], allowWildCard=True) assert (esdList.TrigEsd() == ["JetCollection#HLT_AutoKey*"]) msg.info("Should see a warning about attempt at adding invalid element...") foo = CfgItemList("Foo", items=["Bla#Blah*"]) assert (foo() == []) del foo ## FIXME: AthenaPoolOutputStream is bringing theApp into our namespace... msg.info("Creating an AthenaPoolOutputStream...") from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream StreamESD = AthenaPoolOutputStream("StreamESD") StreamESD.ItemList = esdList() msg.info("ESD item list: %r", esdList()) msg.info("removing TrigEsd item list") del esdList.TrigEsd StreamESD.ItemList = esdList() msg.info("ESD item list: %r", StreamESD.ItemList) caught = False try: esdList += AthenaPoolOutputStream("FooFoo") except TypeError as err: caught = True msg.info(err) msg.info("Error has been caught, good")
# Copyright (C) 2002-2019 CERN for the benefit of the ATLAS collaboration include.block('EventOverlayJobTransforms/OverlayOutput_jobOptions.py') from AthenaCommon.DetFlags import DetFlags from AthenaCommon.AthenaCommonFlags import athenaCommonFlags from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream from OverlayCommonAlgs.OverlayFlags import overlayFlags outStream = AthenaPoolOutputStream('StreamRDO', athenaCommonFlags.PoolRDOOutput(), asAlg=True) outStream.ItemList = [] # Event info outStream.ItemList += [ 'xAOD::EventInfo#EventInfo', 'xAOD::EventAuxInfo#EventInfoAux.' ] # Timings outStream.ItemList += ['RecoTimingObj#EVNTtoHITS_timings'] if not overlayFlags.isDataOverlay(): outStream.ItemList += ['RecoTimingObj#HITStoRDO_timings'] # Truth if DetFlags.overlay.Truth_on(): outStream.ItemList += ['McEventCollection#TruthEvent'] if 'TrackRecordCollection' in overlayFlags.optionalContainerMap(): for collection in overlayFlags.optionalContainerMap( )['TrackRecordCollection']:
if not 'ALGMODE' in dir(): ALGMODE = 'cpp' pass if ALGMODE == 'py': from McParticleTests.Lib import PyMcAodSymLinkTests as McAodSymLinkTests else: from McParticleTests.McParticleTestsConf import McAodSymLinkTests topSequence += McAodSymLinkTests(OutputLevel=DEBUG) #--------------------------------------------------------------- # Pool Persistency #--------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream outStream = AthenaPoolOutputStream("OutStream") outStream.ItemList = ["EventInfo#*"] outStream.ItemList += ["McEventCollection#GEN_EVENT"] outStream.ItemList += ["TruthParticleContainer#SpclMC"] if 'OUTPUT' not in dir(): OUTPUT = "mc.aod.symlinktests.pool" pass outStream.OutputFile = OUTPUT outStream.EvtConversionSvc = "AthenaPoolCnvSvc" outStream.ForceRead = True #force read of output data objs svcMgr.MessageSvc.defaultLimit = 4000000 svcMgr.MessageSvc.OutputLevel = ERROR #============================================================== # # End of job options file
from Starlight_i.Starlight_iConf import Starlight_i genSeq += Starlight_i("Starlight") genSeq.Starlight.McEventKey = "GEN_EVENT" evgenConfig.generators += ["Starlight"] ## Extra stream persistency ## 2101 == EventInfo, 133273 == MCTruth (HepMC) from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream _evgenstream = AthenaPoolOutputStream("StreamEVGEN") _evgenstream.ItemList = ["2101#*","133273#GEN_EVENT"] del _evgenstream
include( "GeneratorObjectsAthenaPool/GeneratorObjectsAthenaPool_joboptions.py" ) include( "LArAthenaPool/LArAthenaPool_joboptions.py" ) include( "G4SimAthenaPOOL/G4SimAthenaPOOL_joboptions.py" ) # Get a handle to the ServiceManager from AthenaCommon.AppMgr import ServiceMgr as svcMgr # Get a handle to the ApplicationManager from AthenaCommon.AppMgr import theApp # # Pool output # from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream( "Stream1" ) Stream1.OutputFile = "PoolOutputFileName Stream1.ItemList=["EventInfo#*", "McEventCollection#TruthEvent"] # mc truth (hepmc) Stream1.ItemList+=["LArHitContainer#*"] Stream1.ForceRead = True # # Pool input # svcMgr.EventSelector.InputCollections = PoolInputCollection # # the Tile, LAr and Calo detector description package # from AthenaCommon.GlobalFlags import jobproperties jobproperties.Global.DetDescrVersion=Geometry from AtlasGeoModel import SetGeometryVersion