topSequence = AlgSequence() from Pythia_i.Pythia_iConf import Pythia topSequence += Pythia() #------------------------------------------------------------- # define your physics process here # Declare the Pythia to be the Generator to be used topSequence.Pythia.PythiaCommand += ["pysubs msel 6"] #selection of the subprocess gg -> ttbar topSequence.Pythia.PythiaCommand += ["pysubs msub 81 1"] topSequence.Pythia.PythiaCommand += ["pysubs msub 82 1"] #------------------------------------------------------------- #------------------------------------------------------------- # Pool Persistency #------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream # EventInfo Converters include("EventAthenaPool/EventAthenaPool_joboptions.py") PoolSvc = Service("PoolSvc") Stream1 = AthenaPoolOutputStream("Stream1") Stream1.ItemList += ["EventInfo#*"] Stream1.ItemList += ["McEventCollection#*"] Stream1.ForceRead = TRUE Stream1.OutputFile = "McEvent.root" #PoolSvc.Output = "McEvent.root" #==============================================================
#theApp.Dlls += [ "RootHistCnv" ] # Change the following line to "ROOT" for ROOT persistency theApp.HistogramPersistency = "ROOT" #-------------------------------------------------------------- # NTuple output file #-------------------------------------------------------------- #RootHistSvc = Service( "RootHistSvc" ) #HbookHistSvc.NPAWC = 1500000 #HistogramPersistencySvc = Service( "HistogramPersistencySvc" ) #HistogramPersistencySvc.OutputFile = "histo.root" svcMgr.NTupleSvc.Output = [ "FILE1 DATAFILE='pythiaB.root' OPT='NEW' TYP='ROOT'" ] ############################################################### # Add POOL persistency from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream include( "GeneratorObjectsAthenaPool/GeneratorObjectsAthenaPool_joboptions.py" ) Stream1 = AthenaPoolOutputStream( "Stream1" ) Stream1.OutputFile = "pythiaB.pool.root" # 2101 = EventInfo # 133273 = MCTruth (HepMC) Stream1.ItemList += [ "2101#*", "133273#*" ] Stream1.AcceptAlgs = [ "BSignalFilter" ] #============================================================== # # End of job options file # ###############################################################
0.01, 0.01, 56.0 ] # mm x, y, z. Default taken from: Simulation/G4Sim/MCTruth/MCTruth/VertexPositioner.hh EventBoost.GaussianMean += [0.0, 0.0, 0.0] # mm x, y, z EventBoost.DoFlatVertexSmearing = False EventBoost.FlatSmearingBoundaryMin += [-0.01, -0.01, -0.01] #mm x, y, z EventBoost.FlatSmearingBoundaryMax += [0.01, 0.01, 0.01] #mm x, y, z EventBoost.OutputLevel = 3 ######### Streams & persistency ######### ################################################## # Setup output stream. # ################################################## from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream myOutStream = AthenaPoolOutputStream("myOutStream") myOutStream.EvtConversionSvc = "AthenaPoolCnvSvc" myOutStream.OutputFile = "boosted.pool.root" myOutStream.ItemList += ["EventInfo#*", "McEventCollection#*"] myOutStream.ForceRead = TRUE ################################################## # Do consistency checks # ################################################## #from TruthExamples.TruthExamplesConf import TestHepMC #job += TestHepMC() #job.TestHepMC.EnergyDifference = 0.001*GeV # tolerance on energy conservation #job.TestHepMC.OutputLevel = WARNING ################################################## # Make CBNT # ##################################################
#-------------------------------------------------------------- # JobOptions for the loading of the AthenaSealSvc #-------------------------------------------------------------- # Check the dictionary in memory for completeness #include( "AthenaSealSvc/AthenaSealSvc_joboptions.py" ) svcMgr.AthenaSealSvc.CheckDictionary = True #-------------------------------------------------------------- # Output options #-------------------------------------------------------------- # Stream's output file from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("Stream1") Stream1.OutputFile = "SimplePoolFile.root" # List of DO's to write out Stream1.ItemList += ["CaloCellContainer#*"] Stream1.ItemList += ["EventInfo#*"] #-------------------------------------------------------------- # Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL ) #-------------------------------------------------------------- svcMgr.MessageSvc.OutputLevel = WARNING svcMgr.MessageSvc.debugLimit = 100000 svcMgr.MessageSvc.errorLimit = 100000 #svcMgr.ClassIDSvc.OutputLevel = DEBUG LArCellContFakeWriter.OutputLevel = DEBUG AthenaEventLoopMgr = Service("AthenaEventLoopMgr") AthenaEventLoopMgr.OutputLevel = INFO
#-------------------------------------------------------------- # 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 # ###############################################################
myAtRndmGenSvc.EventReseeding = False ServiceMgr += myAtRndmGenSvc ## Add an action from G4AtlasApps.SimFlags import simFlags simFlags.OptionalUserActionList.addAction( 'G4UA::MaterialStepRecorderTool'['Run', 'Event', 'Step']) ############### The Material hit collection ################## from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream # --- check dictionary ServiceMgr.AthenaSealSvc.CheckDictionary = True # --- commit interval (test) ServiceMgr.AthenaPoolCnvSvc.OutputLevel = DEBUG ServiceMgr.AthenaPoolCnvSvc.CommitInterval = 10 MaterialStream = AthenaPoolOutputStream('MaterialStream') MaterialStream.OutputFile = "MaterialStepCollection.root" MaterialStream.ItemList += ['Trk::MaterialStepCollection#*'] ############################################################## include("G4AtlasApps/G4Atlas.flat.configuration.py") ## Populate alg sequence from AthenaCommon.CfgGetter import getAlgorithm topSeq += getAlgorithm("G4AtlasAlg", tryDefaultConfigurable=True) include("InDetSLHC_Example/postInclude.SLHC_Setup.py") #--- End jobOptions.GeantinoMapping.py file ------------------------------
from AthenaServices.AthenaServicesConf import ToyNextPassFilterTool ToolSvc += ToyNextPassFilterTool(NPasses=3) #services that need to be reinit 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()
from AthenaCommon.AppMgr import ServiceMgr if not hasattr(ServiceMgr, 'PileUpMergeSvc'): from PileUpTools.PileUpToolsConf import PileUpMergeSvc ServiceMgr += PileUpMergeSvc() # Pool input (Change this to use a different file) from AthenaCommon.AppMgr import ServiceMgr if not hasattr(ServiceMgr, 'EventSelector'): import AthenaPoolCnvSvc.ReadAthenaPool ServiceMgr.EventSelector.InputCollections = PoolHitsInput # Pool Output from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("Stream1") Stream1.OutputFile = PoolHITsOutput Stream1.ItemList+=["EventInfo#*"]; Stream1.ItemList+=["LArHitContainer#*"]; # the LAr and Calo detector description package include( "CaloDetMgrDetDescrCnv/CaloDetMgrDetDescrCnv_joboptions.py" ) include( "LArDetMgrDetDescrCnv/LArDetMgrDetDescrCnv_joboptions.py" ) # Pool Converters include( "LArAthenaPool/LArAthenaPool_joboptions.py" ) #algorithm from AthenaCommon.AlgSequence import AlgSequence topSequence = AlgSequence() from LArDigitization.LArDigitizationConf import LArHitMerger
#-------------------------------------------------------------- # AthenaRoot Persistency #-------------------------------------------------------------- include('GeneratorObjectsAthenaPool/GeneratorObjectsAthenaPool_joboptions.py') from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("Stream1") # 2101 == EventInfo # 133273 == MC Truth (HepMC) # 54790518 == HijigEventParams Stream1.ItemList += ["2101#*"] Stream1.ItemList += ["133273#*"] Stream1.ItemList += ["54790518#*"] Stream1.OutputFile = "hijing.flow.p_pb.pool.root" #============================================================== # # End of job options file # ###############################################################
topSequence += ObjectThinner ####### If the data contains the Truth, then do this, uncomment this lines # Make Truth Jet on the Fly if there is truth in this data #from JetRec.JetGetters import * #from JetRec.JetRecFlags import jetFlags #jetFlags.inputFileType = 'GEN' #make_StandardJetGetter('Cone',0.4,'Truth').jetAlgorithmHandle() # Slim the McEventCollection #include ( "DPDUtils/TruthSlimming.py" ) #################### End the treatment of the Truth ######## from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream StreamDPD = AthenaPoolOutputStream("StreamDPD") StreamDPD.OutputFile = "D2PD.pool.root" #Metadata and luminosity block StreamDPD.ItemList += ["IOVMetaDataContainer#*"] StreamDPD.ItemList += ["LumiBlockCollection#*"] #Event Level StreamDPD.ItemList += ['EventInfo#*'] StreamDPD.ItemList += ['PileUpEventInfo#*'] StreamDPD.ItemList += ['EventStreamInfo#*'] StreamDPD.ItemList += ['MergedEventInfo#*'] # Inner Detector TrackParticles - # for Inner detector TrackParticles associated to Electrons/Photons/Muons/TauJets StreamDPD.ItemList += ['Rec::TrackParticleContainer#TrackParticleCandidate']
FTKDumpCond = FTKDumpCondAlgo("FTKDumpCondAlgo" , OutputLevel=INFO) FTKDumpCond.IBLMode = runArgs.IBLMode FTKDumpCond.DumpBadModules = False FTKDumpCond.DumpGlobalToLocalMap = False FTKDumpCond.DumpModulePositions = True alg += FTKDumpCond else : alg += FTKRoadFinder alg += FTKTrackFitter if not FTKTrackFitter.doTrackFile : from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream( "Stream1" ) Stream1.OutputFile = "FTKSimulationOutput.pool.root" Stream1.ItemList += [ "FTKAthTrackContainer#*", "FTKTrackFitterStats#*", "TrigInDetTrackCollection#*" ] ## Post-include if hasattr(runArgs,"postInclude"): for fragment in runArgs.postInclude: include(fragment) ## Post-exec if hasattr(runArgs,"postExec"): ftkLog.info("transform post-exec") for cmd in runArgs.postExec:
AANTupleStream.ExistDataHeader = False AANTupleStream.OutputLevel = VERBOSE HistogramPersistencySvc = Service("HistogramPersistencySvc") HistogramPersistencySvc.OutputFile = "histo.root" #--------------------------------------------------------------- # Pool Persistency #--------------------------------------------------------------- if not doNotWritePOOL: from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream include( "GeneratorObjectsAthenaPool/GeneratorObjectsAthenaPool_joboptions.py") Stream1 = AthenaPoolOutputStream("Stream1") Stream1.OutputFile = "athena.evgen.root" # 2101 = EventInfo # 133273 = MCTruth (HepMC) Stream1.ItemList += ["2101#*", "133273#*"] Stream1.AcceptAlgs = ["BSignalFilter"] include("AthenaSealSvc/AthenaSealSvc_joboptions.py") AthenaSealSvc.CheckDictionary = TRUE #============================================================== # # End of job options file # ###############################################################
#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()
#-------------------------------------------------------------- # Output options # LAr Hit # Converters: include("EventAthenaPool/EventAthenaPool_joboptions.py") #include( "AFP_EventAthenaPool/AFP_EventAthenaPool_joboptions.py" ) # ItemList: include("EventAthenaPool/EventAthenaPoolItemList_joboptions.py") #include( "AFP_EventAthenaPool/AFP_EventAthenaPoolItemList_joboptions.py" ) print fullItemList from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("Stream1") Stream1.OutputFile = strOutputRootFile # List of DO's to write out - we don't use fullItemList, just list a # few we want #Stream1.ItemList += fullItemList Stream1.ItemList += ["EventInfo#*", "McEventCollection#TruthEvent"] # mc truth (hepmc) #Stream1.ItemList+= ["AFP_TDDigiCollection#*"] Stream1.ItemList += ["AFP_TDLocRecoEvCollection#*"] #Stream1.ItemList+= ["AFP_SiDigiCollection#*"] Stream1.ItemList += ["AFP_SIDLocRecoEvCollection#*"] #-------------------------------------------------------------- # Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL ) MessageSvc = Service("MessageSvc") MessageSvc.OutputLevel = 3
"zLimits -4000. 4000.", # required for shape=cylinder # "ptLimits 0. 0.", # lower and upper limits. The upper limit is optional # "phiLimits 0. 0.", # lower and upper limits. The upper limit is optional # "etaLimits 0. 0.", # lower and upper limits. The upper limit is optional "rLimits 650. 1950." ] # required for shape=cylinder # "weightLimits 0. 0."] # lower and upper limits. The upper limit is optional topSequence += BeamHaloGenerator print BeamHaloGenerator # Required when monitoring plots are enabled. theApp.HistogramPersistency = "ROOT" if not hasattr(ServiceMgr, 'THistSvc'): from GaudiSvc.GaudiSvcConf import THistSvc ServiceMgr += THistSvc() ServiceMgr.THistSvc.Output += [ "BeamHalo DATAFILE='BeamHalo.root' OPT='RECREATE'" ] # For debugging from TruthExamples.TruthExamplesConf import DumpMC topSequence += DumpMC() from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("StreamEVGEN") Stream1.OutputFile = "evgen.pool.root" Stream1.ItemList += ['EventInfo#*', 'McEventCollection#*']
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" #StreamDPD.OutputFile = outFileName StreamDPD.ForceRead=TRUE; #force read of output data objs StreamDPD.AcceptAlgs=["ttbarFilter"] ChronoStatSvc = Service( "ChronoStatSvc" ) ChronoStatSvc.ChronoDestinationCout = True ChronoStatSvc.PrintUserTime = True ChronoStatSvc.PrintSystemTime = True ChronoStatSvc.PrintEllapsedTime = True AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" ) AthenaPoolCnvSvc.UseDetailChronoStat = True MessageSvc = Service("MessageSvc")
if not hasattr(ServiceMgr, 'EventSelector'): import AthenaPoolCnvSvc.ReadAthenaPool ServiceMgr.EventSelector.InputCollections = PoolHitsInput ServiceMgr.EventSelector.SkipEvents = SkipEvents if DetFlags.writeRDOPool.any_on(): # Pool Output from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("Stream1") Stream1.ItemList += ["EventInfo#*"] Stream1.ItemList += ["TBEventInfo#*"] Stream1.ItemList += ["McEventCollection#*"] Stream1.OutputFile = PoolRDOOutput #--- Initial checks ------------------------------------------- #print job configuration DetFlags.Print() #check job configuration if not (DetFlags.geometry.any_on()): AthError("this digitization job needs some DetFlags.geometry On") if not (DetFlags.digitize.any_on()): AthError("this *digitization* job needs some DetFlags.digitize On") if not (DetFlags.writeRDOPool.any_on()): log.warning("this digitization job will not write any RDO object") #--- GeoModel stuff ------------------------------------------- from AthenaCommon.GlobalFlags import globalflags
# Jan24,06 turns ON decays charm and bottom but not pi0, lambda, ... "ihpr2 12 2", # Turns ON retaining of particle history - truth information: "ihpr2 21 1" ] #--------------------------------------------------------------- # Pool Persistency #--------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream #--- Stream1 --- Stream1 = AthenaPoolOutputStream("Stream1") theApp.Dlls += ["GeneratorObjectsAthenaPoolPoolCnv"] PoolSvc = Service("PoolSvc") # 2101 == EventInfo # 133273 == MCTruth (HepMC) # 54790518 == HijigEventParams Stream1.ItemList += ["2101#*"] Stream1.ItemList += ["133273#*"] Stream1.ItemList += ["54790518#*"] Stream1.OutputFile = "hijing.pool.root" #============================================================== # # End of job options file # ###############################################################
from AthenaCommon.AppMgr import ServiceMgr include("AthenaCommon/Atlas_Gen.UnixStandardJob.py") include("PartPropSvc/PartPropSvc.py") # AthenaPool converters include("EventAthenaPool/EventAthenaPool_joboptions.py") include("GeneratorObjectsAthenaPool/GeneratorObjectsAthenaPool_joboptions.py") #include( "AthenaPoolCnvSvc/WriteAthenaPool_jobOptions.py" ) from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("StreamEVGEN") Stream1.ItemList += ["EventInfo#*", "McEventCollection#*"] Stream1.OutputFile = "mc13.007903.cavernbg_sf05_25ns.evgen._0001.pool.root" from CavernBkgGenerator.CavernBkgGeneratorConf import MuonBackgroundConverter topSequence += MuonBackgroundConverter("MuonBackgroundConverter") MuonCavern = topSequence.MuonBackgroundConverter # The following is needed to load the Athena Random # Number Generation Service. from AthenaServices.AthenaServicesConf import AtRndmGenSvc ServiceMgr += AtRndmGenSvc("AtRndmGenSvc", Seeds=["SINGLE 2040160768 790101"]) # AtRndmGenSvc.ReadFromFile = true; ServiceMgr.AtRndmGenSvc.OutputLevel = DEBUG
#-------------------------------------------------------------- # 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" )
## Add an action from G4AtlasApps.SimFlags import simFlags simFlags.OptionalUserActionList.addAction('G4UA::MaterialStepRecorderTool') # suppress the enormous amount of MC output from TruthExamples.TruthExamplesConf import DumpMC DumpMC.VerboseOutput = False ############### The Material hit collection ################## from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream # --- commit interval (test) ServiceMgr.AthenaPoolCnvSvc.OutputLevel = DEBUG MaterialStream = AthenaPoolOutputStream('MaterialStream') MaterialStream.OutputFile = outPath + "MaterStepFile.root" MaterialStream.ItemList += ['MaterialStepVector#*'] ############################################################## include("G4AtlasApps/G4Atlas.flat.configuration.py") ## Populate alg sequence from AthenaCommon.CfgGetter import getAlgorithm topSeq += getAlgorithm("G4AtlasAlg", tryDefaultConfigurable=True) from InDetIBL_Example.SLHC_Setup import SLHC_Setup SLHC_Setup = SLHC_Setup() #--- End jobOptions.GeantinoMapping.py file ------------------------------
include("InDetRawDataByteStreamCnv/PixelWriteBS_jobOptions.py") # ------------------------------------------------------------ # persistency # ------------------------------------------------------------ if doWriteESD or doWriteAOD or ('doCopyRDO' in dir() and doCopyRDO): # --- load setup from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream # --- check dictionary ServiceMgr.AthenaSealSvc.CheckDictionary = True if doWriteESD: # --- create stream StreamESD = AthenaPoolOutputStream("StreamESD") StreamESD.OutputFile = "InDetRecESD.root" # --- save MC collection if truth turned on if InDetFlags.doTruth(): StreamESD.ItemList += ["McEventCollection#*"] # ---- load list of objects include("InDetRecExample/WriteInDetESD.py") if doWriteAOD: # --- create stream StreamAOD = AthenaPoolOutputStream("StreamAOD") StreamAOD.OutputFile = "InDetRecAOD.root" # --- save MC collection if truth turned on if InDetFlags.doTruth(): StreamAOD.ItemList += ["McEventCollection#*"] # --- load list of objects include("InDetRecExample/WriteInDetAOD.py")
job.Pythia.PythiaCommand += [ "grav 1 2", "grav 2 1111", "grav 3 7000", "grav 4 3500", "grav 5 350", "grav 6 3500", "grav 7 2.5E+1", "grav 8 0" ] job.Pythia.PythiaCommand += [ "pyinit pylistf 1", "pyinit dumpr 0 3", "pystat 1" ] #---------------------------------------------------------------- from TruthExamples.TruthExamplesConf import DumpMC #job += DumpMC() #--------------------------------------------------------------- # Ntuple service output #--------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream #from AthenaServices.AthenaServicesConf import AthenaOutputStream from PoolSvc.PoolSvcConf import PoolSvc from AthenaPoolCnvSvc.AthenaPoolCnvSvcConf import AthenaPoolCnvSvc Stream1 = AthenaPoolOutputStream("StreamEVGEN") Stream1.ItemList += ["EventInfo#*", "McEventCollection#*"] Stream1.OutputFile = "Pythia_Graviton.pool.root" #============================================================== # # End of job options file # ###############################################################
job.Pythia.PythiaCommand += [ "pydat3 mdme 89 1 0", #enu "pydat3 mdme 90 1 0" ] #munu from TruthExamples.TruthExamplesConf import PrintMC job += PrintMC() job.PrintMC.FirstEvent = 1 job.PrintMC.LastEvent = 3 #--------------------------------------------------------------- # Pool Persistency #--------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("StreamEVGEN") Stream1.OutputFile = "matchig.pool.root" Stream1.ItemList += ['EventInfo#*', 'McEventCollection#*'] #--------------------------------------------------------------- # Ntuple service output #--------------------------------------------------------------- #include( 'CBNT_Athena/CBNT_AthenaAware_jobOptions.py' ) #include( 'CBNT_Athena/CBNT_EventInfo_jobOptions.py' ) #include( 'RecExCommon/CBNT_Truth_jobOptions.py' ) #theApp.CreateSvc += ['AthenaSealSvc'] #ServiceMgr += CfgMgr.AthenaSealSvc() #--------------------------------------------------------------- # AANT persistency
## get a handle on the ServiceManager from AthenaCommon.AppMgr import ServiceMgr as svcMgr #Explicitly specify the output file catalog from PoolSvc.PoolSvcConf import PoolSvc svcMgr += PoolSvc() svcMgr.PoolSvc.WriteCatalog = "xmlcatalog_file:Catalog1.xml" from AthenaPoolCnvSvc.AthenaPoolCnvSvcConf import AthenaPoolCnvSvc svcMgr += AthenaPoolCnvSvc() svcMgr.AthenaPoolCnvSvc.CommitInterval = 10 from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("Stream1") Stream1.OutputFile = "HLT.root" # ** mandatory parameter ** // The output file name Stream1.ItemList += ["TrigDec::TrigDecision#TrigDecision"] Stream1.ItemList += ["TrigConf::HLTAODConfigData#AODConfig"] Stream1.ItemList += ["TrigConf::Lvl1AODConfigData#AODConfig"] Stream1.ItemList += ["TrigConf::Lvl1AODPrescaleConfigData#AODConfig"] theApp.Dlls += ["TrigEventAthenaPoolPoolCnv"] #-------------------------------------------------------------- # Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL) #-------------------------------------------------------------- svcMgr.MessageSvc.OutputLevel = 3 svcMgr.PoolSvc.OutputLevel = 2 svcMgr.AthenaPoolCnvSvc.OutputLevel = 2 Stream1.OutputLevel = 2
#services that need to be reinit 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()
job += CopyMcEventCollection() #-------------------------------------------------------------- # Output options #-------------------------------------------------------------- from StoreGate.StoreGateConf import StoreGateSvc ServiceMgr += StoreGateSvc("StoreGateSvc") ServiceMgr.StoreGateSvc.Dump = True from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream outStream = AthenaPoolOutputStream("OverlayOutputStream") outStream.Store = ServiceMgr.StoreGateSvc # Stream's output file outStream.OutputFile = "InDetOverlayOut.pool.root" #FIXME: ### FIXME: this does not work because I and O stores are separate. #FIXME: # Copy everything from the input and must force reading of all input #FIXME: # objects #FIXME: #outStream.ForceRead = True #FIXME: #outStream.TakeItemsFromInput = True # Need to set this to False if an explicit list is to be used # Can use an explicit list to avoid broken classes. Requires TakeItemsFromInput=False. outStream.ItemList = [] # Define output stream outStream.ItemList += ["EventInfo#*"] #FIXME: not InDet-specific outStream.ItemList += ["PileUpEventInfo#*"] #FIXME: not InDet-specific outStream.ItemList += ["McEventCollection#*"] #FIXME: not InDet-specific
job.Pyquench.nuclearA = 208 job.Pyquench.minBiasYN = False job.Pyquench.bimpact = 0. job.Pyquench.ptmin = 100. job.Pyquench.sqrtSNN = 5520. # Number of events to be processed (default is 10) theApp.EvtMax = 50 #--------------------------------------------------------------- # Pool Persistency #--------------------------------------------------------------- include("GeneratorObjectsAthenaPool/GeneratorObjectsAthenaPool_joboptions.py") include("AthenaSealSvc/AthenaSealSvc_joboptions.py") AthenaSealSvc.CheckDictionary = TRUE from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("StreamEVGEN") #Stream1.ForceRead=TRUE Stream1.OutputFile = "r1040b0pyquench.root" # 2101 = EventInfo # 133273 = MCTruth (HepMC) Stream1.ItemList += ["2101#*", "133273#*"] #============================================================== # # End of job options file # ###############################################################
include("EventOverlayJobTransforms/MuonMcSignal_jobOptions.py") include("EventOverlayJobTransforms/LArMcSignal_jobOptions.py") include("EventOverlayJobTransforms/TileMcSignal_jobOptions.py") include("EventOverlayJobTransforms/L1Signal_jobOptions.py") from OverlayCommonAlgs.OverlayCommonAlgsConf import UpdateEventInfo job += UpdateEventInfo() job.UpdateEventInfo.InfoType = "McEventInfo" # The output - signal from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream signalStream = AthenaPoolOutputStream("StreamRDO_MC") signalStream.OutputFile = SignalCollection signalStream.Store = "TemporaryStore" signalStream.ItemList += [ "EventInfo#*", "EventStreamInfo#*", "PileUpEventInfo#*" ] if DetFlags.overlay.Truth_on(): signalStream.ItemList += ["McEventCollection#*", "TrackRecordCollection#*"] if DetFlags.overlay.CSC_on(): signalStream.ItemList += ["CscSimDataCollection#CSC_SDO"] if DetFlags.overlay.MDT_on(): signalStream.ItemList += ["MuonSimDataCollection#MDT_SDO"] if DetFlags.overlay.RPC_on(): signalStream.ItemList += ["MuonSimDataCollection#RPC_SDO"] if DetFlags.overlay.TGC_on():
topAlg += PrintMC() #-------------------------------------------------------------- # Event related parameters #-------------------------------------------------------------- # Number of events to be processed (default is 10) theApp.EvtMax = 10 #-------------------------------------------------------------- # Algorithms Private Options #-------------------------------------------------------------- Pythia = topAlg.Pythia Pythia.PythiaCommand = ["pyinit user madgraph"] from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream( "StreamEVGEN" ) Stream1.OutputFile = "baur.pool.root" Stream1.ItemList += [ 'EventInfo#*', 'McEventCollection#*' ] ############################################################### # ############################################################### #from AthenaCommon.AppMgr import ServiceMgr #ServiceMgr.EventSelector.RunNumber = runArgs.runNumber #ServiceMgr.EventSelector.FirstEvent = runArgs.firstEvent #---------------------------------------------------------------- # Ntuple service output #---------------------------------------------------------------- #--------------------------------------------------------------- # Ntuple service output
topAlg = AlgSequence("TopAlg") #from Pythia_i.Pythia_iConf import Pythia include("MC09JobOptions/MC9_Pythia_Common.py") Pythia.PythiaCommand = [ "pyinit user protos", "pydat1 parj 90 20000.", "pydat3 mdcy 15 1 0" ] # ... Tauola include("MC09JobOptions/MC9_Tauola_Fragment.py") # ... Photos include("MC09JobOptions/MC9_Photos_Fragment.py") #from TruthExamples.TruthExamplesConf import DumpMC #topAlg += DumpMC() #--------------------------------------------------------------- # Pool Persistency #--------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream("StreamEVGEN") Stream1.OutputFile = "protos.pool.root" Stream1.ItemList += ['EventInfo#*', 'McEventCollection#*'] # # End of job options file # ###############################################################
ServiceMgr += myAtRndmGenSvc ## add the material step recording action SimFlags.OptionalUserActionList.addAction('G4UA::MaterialStepRecorderTool', ['Run', 'Event', 'Step']) ############### The Material hit collection ################## from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream # --- check dictionary ServiceMgr.AthenaSealSvc.CheckDictionary = True # --- commit interval (test) ServiceMgr.AthenaPoolCnvSvc.OutputLevel = DEBUG ServiceMgr.AthenaPoolCnvSvc.CommitInterval = 10 MaterialStream = AthenaPoolOutputStream('MaterialStream') MaterialStream.OutputFile = "MaterialStepFile.root" MaterialStream.ItemList += ['Trk::MaterialStepCollection#*'] ############################################################## include("G4AtlasApps/G4Atlas.flat.configuration.py") # Add the beam effects algorithm from AthenaCommon.CfgGetter import getAlgorithm topSeq += getAlgorithm("BeamEffectsAlg", tryDefaultConfigurable=True) ## Populate alg sequence topSeq += getAlgorithm("G4AtlasAlg", tryDefaultConfigurable=True) #--- End jobOptions.GeantinoMapping.py file ------------------------------
"root://eosatlas//eos/atlas/atlascerngroupdisk/proj-sit/digitization/RTT/mc10/mc10_7TeV.005008.CavernInput.merge.HITS.e4_e607_s951_s952_tid170551_00/HITS.170551._000111.pool.root.1" ] ServiceMgr += getConfigurable("EventSelectorAthenaPool")("MinBiasSelector" ) ServiceMgr.MinBiasSelector.InputCollections = [ "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