PFRunESDtoAOD_TopoTowers_mc21_14TeV.py

Go to the documentation of this file.

# Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
 
if __name__=="__main__":
 
    from AthenaConfiguration.AllConfigFlags import initConfigFlags
    cfgFlags = initConfigFlags()
    cfgFlags.Concurrency.NumThreads=1
    cfgFlags.Exec.MaxEvents=100
    cfgFlags.Input.isMC=True
    cfgFlags.Input.Files = ["/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/PFlowTests/mc21_14TeV/mc21_14TeV.601229.PhPy8EG_A14_ttbar_hdamp258p75_SingleLep.recon.AOD.e8481_s4383_r15934/AOD.41490164._005514.pool.root.1"]
    cfgFlags.Output.AODFileName="output_AOD.root"
    cfgFlags.Output.doWriteAOD=True
    #We need to use EM scale topotowers for pflow
    cfgFlags.Calo.FwdTower.prepareLCW=False
    cfgFlags.PF.useTopoTowers = True
    #Global links don't currently work with topotowers (other CP objects do not support topotowers)
    cfgFlags.PF.useElPhotLinks = False
    cfgFlags.PF.useMuLinks = False
    #Auto configure works by reading the conditions tag  metadata from the input file.
    #Hence for files produced using pre-CREST conditions tags this choice is not
    #compatible with the latest releases which only support CREST conditions tags.
    #So we force it to use the current Run 4 tag
    #See https://its.cern.ch/jira/browse/ATLASRECTS-8434
    from AthenaConfiguration.TestDefaults import defaultConditionsTags
    cfgFlags.IOVDb.GlobalTag = defaultConditionsTags.RUN4_MC
    cfgFlags.fillFromArgs()
    cfgFlags.lock()
 
    from AthenaConfiguration.MainServicesConfig import MainServicesCfg
    cfg = MainServicesCfg(cfgFlags)
 
    from AthenaPoolCnvSvc.PoolReadConfig import PoolReadCfg
    cfg.merge(PoolReadCfg(cfgFlags))
 
    #Build EM scale combined topocluster + tower container
    from CaloRec.CaloFwdTopoTowerConfig import CaloFwdTopoTowerCfg
    cfg.merge(CaloFwdTopoTowerCfg(cfgFlags,CaloTopoClusterContainerKey="CaloTopoClusters",
                                           towerContainerKey="CaloFwdTopoTowers",
                                           TowerBuilderName="CaloTopoTowerFromClusterMaker",
                                           TowerMergerName="CaloTopoSignalMaker",
                                           signalContainerKey="CaloTopoSignal"))
 
    #Also rebuild LCW scale combined topocluster + tower container - ensures both containers are of same size, 
    #which is an assumption in the pflow code when dealing with LCW weights
    cfg.merge(CaloFwdTopoTowerCfg(cfgFlags,TowerBuilderName="CaloCalTopoTowerFromClusterMaker",
                                           TowerMergerName="CaloCalTopoSignalMaker",
                                           signalContainerKey="CaloCalTopoSignal"))
 
    from eflowRec.PFRun3Config import PFFullCfg
    cfg.merge(PFFullCfg(cfgFlags,runTauReco=True))
 
    from eflowRec.PFRun3Remaps import ListRemaps
 
    list_remaps=ListRemaps(cfg, 'AOD')
    for mapping in list_remaps:
        cfg.merge(mapping)    
 
    cfg.run()