Skip to content
Snippets Groups Projects
Commit ba60521b authored by Frank Winklmeier's avatar Frank Winklmeier
Browse files

Merge branch 'sweep_24.0_main_2024-04-24' into 'main'

2024-04-24: merge of 24.0 into main

Closes ATLASRECTS-7974

See merge request atlas/athena!70844
parents 4cee2046 844b5911
No related branches found
No related tags found
No related merge requests found
Showing
with 0 additions and 1167 deletions
......@@ -22,7 +22,6 @@ atlas_add_dictionary( CaloLocalHadCalibDict
LINK_LIBRARIES ${Boost_LIBRARIES} ${ROOT_LIBRARIES} ${CLHEP_LIBRARIES} AthenaBaseComps AthenaKernel CaloConditions CaloDetDescrLib CaloDmDetDescr CaloEvent CaloGeoHelpers CaloIdentifier CaloSimEvent CaloUtilsLib CxxUtils GaudiKernel GeneratorObjects PathResolver StoreGateLib xAODCaloEvent )
# Install files from the package:
atlas_install_joboptions( share/*.txt share/*.py )
atlas_install_runtime( share/CaloHadDMCoeff_init_v1.txt share/CaloHadDMCoeff_init_v2.txt )
atlas_install_python_modules( python/*.py POST_BUILD_CMD ${ATLAS_FLAKE8} )
atlas_install_scripts( scripts/CaloLCW_tf.py scripts/GetCaloLCWPerformance.py )
#
# set inFileName and outFileName to convert TProfile2D classification to pool
#
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
import IOVDbSvc.IOVDb
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
from CaloLocalHadCalib.CaloLocalHadCalibConf import CaloReadLCClassificationFile
ReadLCClassification = CaloReadLCClassificationFile("ReadLCClassification")
ReadLCClassification.LCClassificationFileName = inFileName
ReadLCClassification.ClassificationKey="EMFracClassify"
topSequence += ReadLCClassification
# Number of events to be processed (default is 10)
theApp.EvtMax = 1
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
from RegistrationServices.OutputConditionsAlg import OutputConditionsAlg
myOCA=OutputConditionsAlg("myOCA",outFileName)
# List of objects container type#key#foldername
myOCA.ObjectList = ["CaloLocalHadCoeff#EMFracClassify#/CALO/HadCalibration/CaloEMFrac"]
myOCA.WriteIOV = True
myOCA.IOVTagList = ["CaloEMFrac2-CSC05-BERT" ]
myOCA.OutputLevel= Lvl.DEBUG
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
MessageSvc = Service( "MessageSvc" )
MessageSvc.OutputLevel = Lvl.INFO
DetectorStore = svcMgr.DetectorStore
DetectorStore.Dump = TRUE
DetectorStore.OutputLevel = Lvl.INFO
IOVDbSvc = Service( "IOVDbSvc" )
IOVDbSvc.OutputLevel = Lvl.INFO
svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myEMFrac05BERT200.db;dbname=OFLP200"
#
# set inFileName and outFileName to convert TProfile2D classification to pool
#
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
import IOVDbSvc.IOVDb
if not 'inFileNames' in dir():
inFileNames = ["AntiKt4CB_JES_NumInv_JX.root", "AntiKt6CB_JES_NumInv_JX.root" ]
if not 'inJetCollectionNames' in dir():
inJetCollectionNames = ["AntiKt4LCTopoJets", "AntiKt4LCTopoJets" ]
if not 'outFileName' in dir():
outFileName = "AntiKt4a6CB_JES_40_NumInv_JX.pool.root"
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
from CaloLocalHadCalib.CaloLocalHadCalibConf import CaloReadLCJetEnergyScaleFile
ReadLCJetEnergyScale = CaloReadLCJetEnergyScaleFile("ReadLCJetEnergyScale")
ReadLCJetEnergyScale.LCJetEnergyScaleFileNames = inFileNames
ReadLCJetEnergyScale.LCJetEnergyScaleJetCollectionNames = inJetCollectionNames
ReadLCJetEnergyScale.CorrectionKey="JESCorrection"
topSequence += ReadLCJetEnergyScale
# Number of events to be processed (default is 10)
theApp.EvtMax = 1
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
from RegistrationServices.OutputConditionsAlg import OutputConditionsAlg
myOCA=OutputConditionsAlg("myOCA",outFileName)
# List of objects container type#key#foldername
myOCA.ObjectList = ["CaloLocalHadCoeff#JESCorrection#/CALO/HadCalibration2/CaloJetEnergyScale"]
myOCA.WriteIOV = True
myOCA.IOVTagList = ["CaloHadJESCorr2-GEO18-QGSP-BERT-MU00" ]
myOCA.OutputLevel= Lvl.DEBUG
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
MessageSvc = Service( "MessageSvc" )
MessageSvc.OutputLevel = Lvl.INFO
DetectorStore = svcMgr.DetectorStore
DetectorStore.Dump = TRUE
DetectorStore.OutputLevel = Lvl.INFO
IOVDbSvc = Service( "IOVDbSvc" )
IOVDbSvc.OutputLevel = Lvl.INFO
svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myJESGEO18BERTMU00200.db;dbname=OFLP200"
#
# set inFileName and outFileName to convert TProfile2D classification to pool
#
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
import IOVDbSvc.IOVDb
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
from CaloLocalHadCalib.CaloLocalHadCalibConf import CaloReadLCOutOfClusterFile
ReadLCOutOfCluster = CaloReadLCOutOfClusterFile("ReadLCOutOfCluster")
ReadLCOutOfCluster.LCOutOfClusterFileName = inFileName
ReadLCOutOfCluster.CorrectionKey="OOCCorrection"
topSequence += ReadLCOutOfCluster
# Number of events to be processed (default is 10)
theApp.EvtMax = 1
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
from RegistrationServices.OutputConditionsAlg import OutputConditionsAlg
myOCA=OutputConditionsAlg("myOCA",outFileName)
# List of objects container type#key#foldername
myOCA.ObjectList = ["CaloLocalHadCoeff#OOCCorrection#/CALO/HadCalibration/CaloOutOfCluster"]
myOCA.WriteIOV = True
myOCA.IOVTagList = ["CaloHadOOCCorr2-CSC05-BERT" ]
myOCA.OutputLevel= Lvl.DEBUG
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
MessageSvc = Service( "MessageSvc" )
MessageSvc.OutputLevel = Lvl.INFO
DetectorStore = svcMgr.DetectorStore
DetectorStore.Dump = TRUE
DetectorStore.OutputLevel = Lvl.INFO
IOVDbSvc = Service( "IOVDbSvc" )
IOVDbSvc.OutputLevel = Lvl.INFO
svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myOOC05BERT200.db;dbname=OFLP200"
#
# set inFileName and outFileName to convert TProfile2D classification to pool
#
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
import IOVDbSvc.IOVDb
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
from CaloLocalHadCalib.CaloLocalHadCalibConf import CaloReadLCOutOfClusterFile
ReadLCOutOfCluster = CaloReadLCOutOfClusterFile("ReadLCOutOfCluster")
ReadLCOutOfCluster.LCOutOfClusterFileName = inFileName
ReadLCOutOfCluster.CorrectionKey="OOCPi0Correction"
topSequence += ReadLCOutOfCluster
# Number of events to be processed (default is 10)
theApp.EvtMax = 1
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
from RegistrationServices.OutputConditionsAlg import OutputConditionsAlg
myOCA=OutputConditionsAlg("myOCA",outFileName)
# List of objects container type#key#foldername
myOCA.ObjectList = ["CaloLocalHadCoeff#OOCPi0Correction#/CALO/HadCalibration/CaloOutOfClusterPi0"]
myOCA.WriteIOV = True
myOCA.IOVTagList = ["CaloHadOOCCorrPi02-CSC05-BERT" ]
myOCA.OutputLevel= Lvl.DEBUG
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
MessageSvc = Service( "MessageSvc" )
MessageSvc.OutputLevel = Lvl.INFO
DetectorStore = svcMgr.DetectorStore
DetectorStore.Dump = TRUE
DetectorStore.OutputLevel = Lvl.INFO
IOVDbSvc = Service( "IOVDbSvc" )
IOVDbSvc.OutputLevel = Lvl.INFO
svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myOOCPi005BERT200.db;dbname=OFLP200"
#
# set inFileName and outFileName to convert TProfile2D weights to pool
#
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
import IOVDbSvc.IOVDb
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
from CaloLocalHadCalib.CaloLocalHadCalibConf import CaloReadLCWeightsFile
ReadLCWeights = CaloReadLCWeightsFile("ReadLCWeights")
ReadLCWeights.LCWeightFileName = inFileName
ReadLCWeights.CorrectionKey="H1ClusterCellWeights"
topSequence += ReadLCWeights
# Number of events to be processed (default is 10)
theApp.EvtMax = 1
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
from RegistrationServices.OutputConditionsAlg import OutputConditionsAlg
myOCA=OutputConditionsAlg("myOCA",outFileName)
# List of objects container type#key#foldername
myOCA.ObjectList = ["CaloLocalHadCoeff#H1ClusterCellWeights#/CALO/HadCalibration2/H1ClusterCellWeights"]
myOCA.WriteIOV = True
myOCA.IOVTagList = ["CaloH1CellWeights2-GEO08-QGSP-BERT-Log" ]
myOCA.OutputLevel= Lvl.DEBUG
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
MessageSvc = Service( "MessageSvc" )
MessageSvc.OutputLevel = Lvl.INFO
DetectorStore = svcMgr.DetectorStore
DetectorStore.Dump = TRUE
DetectorStore.OutputLevel = Lvl.INFO
iovdbsvc = Service( "IOVDbSvc" )
IOVDbSvc.OutputLevel = Lvl.INFO
svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myWeightGeo08BERTLog200.db;dbname=OFLP200"
#
# set inFileName and outFileName to convert TProfile2D classification to pool
#
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
import IOVDbSvc.IOVDb
if not 'inFileNameC' in dir():
inFileNameC="LCWeights_h6prod.100004.piplusminuszero_logE.qgspbert.jid000040-50_v3.Phi5._sum.classify.root"
#inFileNameC="LCWeights_h6prod.100004.piplusminuszero_logE.ftfpbert.jid000054-69_v3._sum.classify.root"
if not 'inFileNameW' in dir():
inFileNameW="LCWeights_h6prod.100004.piplusminus_logE.qgspbert.jid000040-50_v3.Phi5._sum.nonInv.weights.root"
#inFileNameW="LCWeights_h6prod.100004.piplusminus_logE.ftfpbert.jid000054-69_v3._sum.nonInv.weights.root"
if not 'inFileNameO' in dir():
inFileNameO="LCWeights_h6prod.100004.piplusminus_logE.qgspbert.jid000040-50_v3.Phi5._sum.ooc.root"
#inFileNameO="LCWeights_h6prod.100004.piplusminus_logE.ftfpbert.jid000054-69_v3._sum.ooc.root"
if not 'inFileNameOPi0' in dir():
inFileNameOPi0="LCWeights_h6prod.100004.pizero_logE.qgspbert.jid000040-50_v3.Phi5._sum.ooc.root"
#inFileNameOPi0="LCWeights_h6prod.100004.pizero_logE.ftfpbert.jid000054-69_v3._sum.ooc.root"
if not 'outFileName' in dir():
outFileName="LCWeights_h6prod.QGSPBERT.100004.vv08.jid000040-50_v3.Phi5.piplusminuszero_logE.classify_weights_ooc_oocPi0.pool.root"
#outFileName="LCWeights_h6prod.FTFPBERT.vv05.100004.piplusminuszero_logE.classify_weights_ooc_oocPi0.pool.root"
if not os.path.exists(inFileNameC):
exit("No file "+inFileNameW)
if not os.path.exists(inFileNameW):
exit("No file "+inFileNameW)
if not os.path.exists(inFileNameO):
exit("No file "+inFileNameO)
if not os.path.exists(inFileNameOPi0):
exit("No file "+inFileNameOPi0)
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
from CaloLocalHadCalib.CaloLocalHadCalibConf import CaloReadLCClassificationFile, CaloReadLCWeightsFile, CaloReadLCOutOfClusterFile
ReadLCClassification = CaloReadLCClassificationFile("ReadLCClassification")
ReadLCWeights = CaloReadLCWeightsFile("ReadLCWeights")
ReadLCOutOfCluster = CaloReadLCOutOfClusterFile("ReadLCOutOfCluster")
ReadLCOutOfClusterPi0 = CaloReadLCOutOfClusterFile("ReadLCOutOfClusterPi0")
ReadLCClassification.LCClassificationFileName = inFileNameC
ReadLCClassification.ClassificationKey="EMFracClassify"
topSequence += ReadLCClassification
ReadLCWeights.LCWeightFileName = inFileNameW
ReadLCWeights.CorrectionKey="H1ClusterCellWeights"
topSequence += ReadLCWeights
ReadLCOutOfCluster.LCOutOfClusterFileName = inFileNameO
ReadLCOutOfCluster.CorrectionKey="OOCCorrection"
topSequence += ReadLCOutOfCluster
ReadLCOutOfClusterPi0.LCOutOfClusterFileName = inFileNameOPi0
ReadLCOutOfClusterPi0.CorrectionKey="OOCPi0Correction"
topSequence += ReadLCOutOfClusterPi0
# Number of events to be processed (default is 10)
theApp.EvtMax = 1
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
from RegistrationServices.OutputConditionsAlg import OutputConditionsAlg
myOCA=OutputConditionsAlg("myOCA",outFileName)
# List of objects container type#key#foldername
myOCA.ObjectList = [
"CaloLocalHadCoeff#EMFracClassify#/CALO/HadCalibration2/CaloEMFrac",
"CaloLocalHadCoeff#H1ClusterCellWeights#/CALO/HadCalibration2/H1ClusterCellWeights",
"CaloLocalHadCoeff#OOCCorrection#/CALO/HadCalibration2/CaloOutOfCluster",
"CaloLocalHadCoeff#OOCPi0Correction#/CALO/HadCalibration2/CaloOutOfClusterPi0"
]
myOCA.WriteIOV = True
#myOCA.IOVTagList = [
#"CaloEMFrac2-H6-2004-01-FTFP-BERT",
#"CaloH1CellWeights2-H6-2004-01-FTFP-BERT",
#"CaloHadOOCCorr2-H6-2004-01-FTFP-BERT",
#"CaloHadOOCCorrPi02-H6-2004-01-FTFP-BERT"
#]
#myOCA.IOVTagList = [
#"CaloEMFrac2-H6-2004-01-QGSP-BERT",
#"CaloH1CellWeights2-H6-2004-01-QGSP-BERT",
#"CaloHadOOCCorr2-H6-2004-01-QGSP-BERT",
#"CaloHadOOCCorrPi02-H6-2004-01-QGSP-BERT"
#]
#myOCA.IOVTagList = [
#"CaloEMFrac2-GEO08-QGSP-BERT-NClus",
#"CaloH1CellWeights2-GEO08-QGSP-BERT-NClus",
#"CaloHadOOCCorr2-GEO08-QGSP-BERT-ISol-NClus",
#"CaloHadOOCCorrPi02-GEO08-QGSP-BERT-ISol-NClus"
#]
#myOCA.IOVTagList = [
# "CaloEMFrac2-GEO08-FTFP-BERT-NClus",
# "CaloH1CellWeights2-GEO08-FTFP-BERT-NClus",
# "CaloHadOOCCorr2-GEO08-FTFP-BERT-NClus",
# "CaloHadOOCCorrPi02-GEO08-FTFP-BERT-NClus"
# ]
#myOCA.IOVTagList = [
#"CaloEMFrac2-GEO16-QGSP-BERT-Side",
#"CaloH1CellWeights2-GEO16-QGSP-BERT-Side",
#"CaloHadOOCCorr2-GEO16-QGSP-BERT-Side",
#"CaloHadOOCCorrPi02-GEO16-QGSP-BERT-Side"
#]
myOCA.IOVTagList = [
"CaloEMFrac2-GEO180100-QGSP-BERT",
"CaloH1CellWeights2-GEO180100-QGSP-BERT",
"CaloHadOOCCorr2-GEO180100-QGSP-BERT",
"CaloHadOOCCorrPi02-GEO180100-QGSP-BERT"
]
myOCA.OutputLevel= Lvl.DEBUG
AthenaPoolCnvSvc.PoolAttributes += [ "STREAM_MEMBER_WISE = '0'" ]
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
MessageSvc = Service( "MessageSvc" )
MessageSvc.OutputLevel = Lvl.INFO
DetectorStore = svcMgr.DetectorStore
DetectorStore.Dump = TRUE
DetectorStore.OutputLevel = Lvl.INFO
IOVDbSvc = Service( "IOVDbSvc" )
IOVDbSvc.OutputLevel = Lvl.INFO
#svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myH60401FTFPBERT200.db;dbname=OFLP200"
#svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myH60401QGSPBERT200.db;dbname=OFLP200"
#svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myGeo08QGSPBERTNClusIsol200.db;dbname=OFLP200"
#svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myGeo08FTFPBERTNClus200.db;dbname=OFLP200"
#svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myGeo16QGSPBERT200.db;dbname=OFLP200"
svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myGeo18QGSPBERT200.db;dbname=OFLP200"
# $Id: GetLCClassification_jobOptions.py,v 1.1.1.1 2008-11-04 08:56:11 menke Exp $
#
# File: CaloLocalHadCalib/share/GetLCClassification_jobOptions.py
# Created: October 2008, S.Menke
# Purpose: Generate Classification Histos in a athena job
#
from CaloLocalHadCalib.CaloLocalHadCalibConf import GetLCClassification
from AthenaCommon.SystemOfUnits import deg, GeV, MeV
from AthenaCommon.AlgSequence import AlgSequence
from math import ( pi as m_pi, log10 as m_log10 )
GetLCC = GetLCClassification("GetLCC")
GetLCC.OutputFileName = outFileNameLCC
#
# Example how to set dimensions for classiciation:
#
#GetLCC.ClassificationDimensions = {
# 'side': ('side',-1.5,1.5,1),
# '|eta|': ('|eta|',0.,5.,25),
# 'phi': ('phi',-m_pi,m_pi,1),
# 'log10(E_clus (MeV))': ('log10(E_clus (MeV))',m_log10(200*MeV),m_log10(1*TeV),13),
# 'log10(<rho_cell (MeV/mm^3)>)-log10(E_clus (MeV))': ('log10(<rho_cell (MeV/mm^3)>)-log10(E_clus (MeV))',-9.0,-4.0,20),
# 'log10(lambda_clus (mm))': ('log10(lambda_clus (mm))',0.0,4.0,20)}
#
GetLCC.ClusterCollectionName = "CaloTopoClusters"
topSequence += GetLCC
# File: CaloLocalHadCalib/share/GetLCDeadMaterialTree_jobOptions.py
# Created: Febryary 2009, G. Pospelov
# Purpose: Generate special small DeadMaterialTree from ESD
from CaloLocalHadCalib.CaloLocalHadCalibConf import GetLCDeadMaterialTree
from AthenaCommon.SystemOfUnits import deg, GeV, MeV
from AthenaCommon.AlgSequence import AlgSequence
GetLCDM = GetLCDeadMaterialTree("GetLCDM")
GetLCDM.HadDMCoeffInitFile = "CaloHadDMCoeff_init_v2.txt"
GetLCDM.ClusterCollectionName = "CaloTopoClusters"
GetLCDM.ClusterCollectionNameCalib = "CaloCalTopoClusters"
GetLCDM.doSaveCalibClusInfo = False # to save additional info from collection with calibrated clusters
if 'outFileNameLCDM' in dir():
GetLCDM.OutputFileName = outFileNameLCDM
else:
GetLCDM.OutputFileName = "DeadMaterialTree.root"
topSequence += GetLCDM
# $Id: GetLCOutOfCluster_jobOptions.py,v 1.1.1.1 2008-11-04 08:56:11 menke Exp $
#
# File: CaloLocalHadCalib/share/GetLCOutOfCluster_jobOptions.py
# Created: October 2008, S.Menke
# Purpose: Generate OutOfCluster Histos in a athena job
#
from CaloLocalHadCalib.CaloLocalHadCalibConf import GetLCOutOfCluster
from AthenaCommon.SystemOfUnits import deg, GeV, MeV
from AthenaCommon.AlgSequence import AlgSequence
from math import ( pi as m_pi, log10 as m_log10 )
GetLCO = GetLCOutOfCluster("GetLCO")
GetLCO.OutputFileName = outFileNameLCO
#
# Example how to set dimensions for out-of-cluster corrections:
#
#GetLCO.OutOfClusterDimensions = {
# 'side': ('side',-1.5,1.5,1),
# '|eta|': ('|eta|',0.,5.,50),
# 'phi': ('phi',-m_pi,m_pi,1),
# 'log10(E_clus (MeV))': ('log10(E_clus (MeV))',m_log10(200*MeV),m_log10(1*TeV),22),
# 'log10(lambda_clus (mm))': ('log10(lambda_clus (mm))',0.0,4.0,20),
# 'weight': ('weight',0.,5.,1)}
#
GetLCO.ClusterCollectionName = "CaloTopoClusters"
topSequence += GetLCO
# Athena job to get performance plots for local hadronic calibration on single pions
# (runs with 'athena CaloLocalHadCalib/GetLCSinglePionsPerf.py' )
#outFileName="LCSinglePionsPerformance.root"
if not 'inFileName' in dir():
#inFileName=['/tmp/pospelov/data/ESD.170972._001733.pool.root.1']
#inFileName=['ESD.371531._000576.pool.root.1']
inFileName=['ESD.371530._000247.pool.root.1']
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
## load POOL support
import AthenaPoolCnvSvc.ReadAthenaPool
## general job configuration
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
#DetDescrVersion = 'ATLAS-CSC-05-00-00'
#DetDescrVersion = 'ATLAS-GEO-08-00-01'
#DetDescrVersion = 'ATLAS-GEO-16-00-00'
#DetDescrVersion = 'ATLAS-GEO-18-01-00'
DetDescrVersion = 'ATLAS-R2-2015-03-01-00'
include("RecExCond/AllDet_detDescr.py")
#from AthenaCommon.BeamFlags import jobproperties
#jobproperties.Beam.numberOfCollisions = 8.0
## configure object key store to recognize calo cells
from RecExConfig.ObjKeyStore import ObjKeyStore, objKeyStore
oks = ObjKeyStore()
oks.addStreamESD('CaloCellContainer', ['AllCalo'] )
## re-do topo clusters on EM scale
from CaloRec.CaloTopoClusterFlags import jobproperties
jobproperties.CaloTopoClusterFlags.doTopoClusterLocalCalib = True
jobproperties.CaloTopoClusterFlags.doCellWeightCalib = False
from CaloRec.CaloClusterTopoGetter import CaloClusterTopoGetter, addSnapshot
CaloClusterTopoGetter()
from CaloRec.CaloRecConf import CaloClusterSnapshot
TopoCalibMoments=topSequence.CaloTopoCluster.TopoCalibMoments
TopoCalibMoments.MomentsNames += ["ENG_CALIB_OUT_L"]
topSequence.CaloTopoCluster.TopoCalibMoments.MatchDmType = 1 # 1=loose, 2=medium (default), 3=tight
addSnapshot("OOCPi0Calib","CaloOOCTopoCluster")
addSnapshot("LocalCalib","CaloWTopoCluster")
theApp.EvtMax = -1
#--------------------------------------------------------------
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
#--------------------------------------------------------------
MessageSvc = svcMgr.MessageSvc
MessageSvc.OutputLevel = Lvl.INFO
MessageSvc.infoLimit = 100000
svcMgr.EventSelector.InputCollections = inFileName
include ("CaloLocalHadCalib/GetLCSinglePionsPerf_jobOptions.py")
lcPerf.useRecoEfficiency = True
print topSequence
ChronoStatSvc = Service( "ChronoStatSvc" )
ChronoStatSvc.ChronoDestinationCout = True
ChronoStatSvc.PrintUserTime = True
ChronoStatSvc.PrintSystemTime = True
ChronoStatSvc.PrintEllapsedTime = True
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
AthenaPoolCnvSvc.UseDetailChronoStat = True
#include("RecExCommon/RecoUtils.py")
svcMgr.StoreGateSvc.Dump=True
# This is the job options file to get performance plots for local hadronic
# calibration on single pions
from CaloLocalHadCalib.CaloLocalHadCalibConf import GetLCSinglePionsPerf
from AthenaCommon.SystemOfUnits import deg, GeV, MeV
from AthenaCommon.AlgSequence import AlgSequence
from math import pi
lcPerf = GetLCSinglePionsPerf("LocalHadPerformance")
# collection name to study cluster moments
lcPerf.ClusterBasicCollName = "CaloTopoClusters"
# collections names to study engReco wrt Truth after different correction steps
lcPerf.ClusterCollectionNames = ["CaloTopoClusters", "CaloWTopoCluster", "CaloOOCTopoCluster", "CaloCalTopoClusters"]
#lcPerf.ClusterCollectionNames = ["CaloTopoCluster_EMAll", "CaloWTopoCluster_EMAll", "CaloOOCTopoCluster_EMAll", "CaloCalTopoCluster_EMAll"]
#lcPerf.ClusterCollectionNames = ["CaloTopoCluster_HadAll", "CaloWTopoCluster_HadAll", "CaloOOCTopoCluster_HadAll", "CaloCalTopoCluster_HadAll"]
lcPerf.CalibrationHitContainerNames= ["LArCalibrationHitInactive","LArCalibrationHitActive","TileCalibrationCellHitCnt"]
lcPerf.DMCalibrationHitContainerNames= ["LArCalibrationHitDeadMaterial", "TileCalibrationDMHitCnt"]
lcPerf.DistanceCut = 1.5
lcPerf.doEngRecOverTruth = True
lcPerf.doEngTag = True
lcPerf.doEngRecSpect = True
lcPerf.doEngNoiseClus = True
lcPerf.doClusMoments = True
lcPerf.doRecoEfficiency = True
lcPerf.etamin = 0.0
lcPerf.etamax = 5.0
lcPerf.netabin = 25
lcPerf.phimin = -pi
lcPerf.phimax = pi
lcPerf.nphibin = 1
lcPerf.logenermin = 2.0
lcPerf.logenermax = 6.4
lcPerf.nlogenerbin = 22
lcPerf.isTestbeam = False
lcPerf.useGoodClus = False
lcPerf.usePionClustersOnly = False
lcPerf.useRecoEfficiency = False
if 'outFileName' in dir():
lcPerf.OutputFileName = outFileName
else:
lcPerf.OutputFileName = "LCSinglePionsPerformance.root"
topSequence += lcPerf
# $Id: GetLCWeights_jobOptions.py,v 1.2 2008-12-08 16:57:41 menke Exp $
#
# File: CaloLocalHadCalib/share/GetLCWeights_jobOptions.py
# Created: October 2008, S.Menke
# Purpose: Generate Weight Histos in a athena job
#
from CaloLocalHadCalib.CaloLocalHadCalibConf import GetLCWeights
from AthenaCommon.SystemOfUnits import deg, GeV, MeV, TeV
from AthenaCommon.AlgSequence import AlgSequence
from math import ( pi as m_pi, log10 as m_log10 )
GetLCW = GetLCWeights("GetLCW")
GetLCW.OutputFileName = outFileNameLCW
GetLCW.UseInversionMethod = True
#
# Example how to set dimensions for two samplings (here EMB1 & 2):
#
# GetLCW.SamplingDimensions = {
# 'EMB1:EMB1': ('EMB1',0.5,1.5,1),
# 'EMB1:side': ('side',-1.5,1.5,1),
# 'EMB1:|eta|': ('|eta|',0.,1.6,16),
# 'EMB1:phi': ('phi',-m_pi,m_pi,1),
# 'EMB1:log10(E_clus (MeV))': ('log10(E_clus (MeV))',m_log10(200*MeV),m_log10(1*TeV),22),
# 'EMB1:log10(rho_cell (MeV/mm^3))': ('log10(rho_cell (MeV/mm^3))',-7.0,1.0,20),
# 'EMB1:weight': ('weight',-2.0,3.0,1),
# 'EMB2:EMB2': ('EMB2',1.5,2.5,1),
# 'EMB2:side': ('side',-1.5,1.5,1),
# 'EMB2:|eta|': ('|eta|',0.,1.6,16),
# 'EMB2:phi': ('phi',-m_pi,m_pi,1),
# 'EMB2:log10(E_clus (MeV))': ('log10(E_clus (MeV))',m_log10(200*MeV),m_log10(1*TeV),22),
# 'EMB2:log10(rho_cell (MeV/mm^3))': ('log10(rho_cell (MeV/mm^3))',-7.0,1.0,20),
# 'EMB2:weight': ('weight',-2.0,3.0,1)}
#
GetLCW.ClusterCollectionName = "CaloTopoClusters"
GetLCW.CalibrationHitContainerNames = ["LArCalibrationHitInactive"
,"LArCalibrationHitActive"
,"TileCalibHitInactiveCell"
,"TileCalibHitActiveCell"]
topSequence += GetLCW
import user
import ROOT
import cppyy
import re
import math
import glob,os,sys
os.environ["DISPLAY"]=''
def calcWeights(infile, outfile):
fnew = ROOT.TFile(outfile,"RECREATE")
histolist = []
f = ROOT.TFile.Open(infile)
l = f.GetListOfKeys()
while l.GetEntries() > 0:
key = l.First()
if re.compile('^inv_weight').search(key.GetName(),0):
mw=re.compile('weight').search(key.GetName(),0)
ROOT.gROOT.cd()
hiw = ROOT.TProfile2D()
f.GetObject(key.GetName(),hiw)
ROOT.gROOT.cd()
print "create new histo for key", key.GetName()
ROOT.gROOT.cd()
hnn = ROOT.TProfile2D(hiw)
hnn.Reset()
hnn.SetName(mw.string[mw.start():])
hnn.SetTitle(mw.string[mw.start():])
hnn.SetYTitle("log10(#rho_{cell} (MeV/mm^{3}))")
hnn.SetZTitle("E_{tot}/E_{reco}");
hnn.SetDirectory(ROOT.gROOT)
histolist += [hnn]
xax=hiw.GetXaxis()
yax=hiw.GetYaxis()
sumw = hiw.GetSum()
ntot = hiw.GetEntries()
for ix in range(0,hiw.GetNbinsX()):
for iy in range(0,hiw.GetNbinsY()):
iBin = hiw.GetBin(ix+1,iy+1)
if hiw.GetBinContent(iBin) > 0:
leclus=xax.GetBinCenter(ix+1)
invratio = hiw.GetBinContent(iBin)
ratio = 1.0/invratio
inverr = hiw.GetBinError(iBin)
n = int(hiw.GetBinEntries(iBin)*ntot/sumw)
if n > 0:
err = inverr/invratio*ratio/math.sqrt(n)
if err < 1.0*ratio:
ymin = yax.GetBinLowEdge(iy+1)
ymax = yax.GetBinUpEdge(iy+1)
lw = math.log10(invratio)
for ie in range(0,n):
lhro = ymin + (ie+0.5)*(ymax-ymin)/n
dr = err
if ie%2 == 1:
dr = -dr
pass
hnn.Fill(leclus,lhro+lw,ratio+dr)
pass
pass
pass
pass
pass
pass
pass
l.RemoveFirst()
pass
f.Close()
f.Delete()
fnew.cd()
for x in histolist:
x.SetDirectory(fnew)
x.Write()
fnew.Close()
fnew.Delete()
#os.system('rm -rf '+tempdir)
###############################################################################
#
###############################################################################
if __name__ == '__main__':
if len(sys.argv) != 3:
exit("Usage 'python calcWeights.py infile outfile'")
else:
infile = sys.argv[1]
outfile = sys.argv[2]
if not os.path.exists(infile):
os.system("No such file "+infile)
calcWeights(infile, outfile)
# This is the job options file to create classifcation, weighting and
# out-of-cluster histos from ESD
###############################################################
#
# Job options file
#
# Authors: Sven Menke
#==============================================================
doClassification=True
doWeighting=True
doOutOfCluster=True
doDeadMaterial=True
if not 'ClassificationType' in dir():
ClassificationType="None"
#ClassificationType="ParticleID_HAD"
if not 'inFileName' in dir():
inFileName=['ESD.371530._000247.pool.root.1']
if not 'outFileNameLCC' in dir():
outFileNameLCC='classify.root'
#outFileNameLCC=ClassificationType+'_'+'classify.root'
if not 'outFileNameLCW' in dir():
outFileNameLCW='weights.root'
#outFileNameLCW=ClassificationType+'_'+'weights.root'
if not 'outFileNameLCO' in dir():
outFileNameLCO='ooc.root'
#outFileNameLCO=ClassificationType+'_'+'ooc.root'
if not 'outFileNameLCDM' in dir():
outFileNameLCDM='dmc.root'
#outFileNameLCDM=ClassificationType+'_'+'DeadMaterialTree.root'
#DetDescrVersion = 'ATLAS-GEO-08-00-01'
#DetDescrVersion = 'ATLAS-GEO-16-00-00'
DetDescrVersion = 'ATLAS-GEO-18-01-00'
include("RecExCond/AllDet_detDescr.py")
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
## load POOL support
import AthenaPoolCnvSvc.ReadAthenaPool
## general job configuration
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
#from AthenaCommon.BeamFlags import jobproperties
#jobproperties.Beam.numberOfCollisions = 8.0
## configure object key store to recognize calo cells
from RecExConfig.ObjKeyStore import ObjKeyStore, objKeyStore
oks = ObjKeyStore()
oks.addStreamESD('CaloCellContainer', ['AllCalo'] )
## re-do topo clusters on EM scale
from CaloRec.CaloTopoClusterFlags import jobproperties
jobproperties.CaloTopoClusterFlags.doTopoClusterLocalCalib = False
jobproperties.CaloTopoClusterFlags.doCellWeightCalib = False
from CaloRec.CaloClusterTopoGetter import CaloClusterTopoGetter
CaloClusterTopoGetter()
topSequence.CaloTopoCluster.TopoCalibMoments.MomentsNames += ["ENG_CALIB_OUT_L","ENG_CALIB_FRAC_EM"]
topSequence.CaloTopoCluster.TopoCalibMoments.MatchDmType = 1 # 1=loose, 2=medium (default), 3=tight
theApp.EvtMax = -1
#--------------------------------------------------------------
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
#--------------------------------------------------------------
MessageSvc = svcMgr.MessageSvc
MessageSvc.OutputLevel = Lvl.INFO
MessageSvc.infoLimit = 1000000
svcMgr.EventSelector.InputCollections = inFileName
if not 'doClassification' in dir():
doClassification=False
pass
if not 'doWeighting' in dir():
doWeighting=False
pass
if not 'doOutOfCluster' in dir():
doOutOfCluster=False
pass
if not 'doDeadMaterial' in dir():
doDeadMaterial=False
pass
if doClassification:
include ("CaloLocalHadCalib/GetLCClassification_jobOptions.py")
GetLCC.ClassificationType=ClassificationType
pass
if doWeighting:
include ("CaloLocalHadCalib/GetLCWeights_jobOptions.py")
GetLCW.ClassificationType=ClassificationType
pass
if doOutOfCluster:
include ("CaloLocalHadCalib/GetLCOutOfCluster_jobOptions.py")
GetLCO.ClassificationType=ClassificationType
pass
if doDeadMaterial:
include ("CaloLocalHadCalib/GetLCDeadMaterialTree_jobOptions.py")
pass
print topSequence
ChronoStatSvc = Service( "ChronoStatSvc" )
ChronoStatSvc.ChronoDestinationCout = True
ChronoStatSvc.PrintUserTime = True
ChronoStatSvc.PrintSystemTime = True
ChronoStatSvc.PrintEllapsedTime = True
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
AthenaPoolCnvSvc.UseDetailChronoStat = True
# This is the job options file to create classifcation histos from ESD
###############################################################
#
# Job options file
#
# Authors: Sven Menke
#==============================================================
DetDescrVersion = 'ATLAS-GEO-08-00-01'
include("RecExCond/AllDet_detDescr.py")
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
## load POOL support
import AthenaPoolCnvSvc.ReadAthenaPool
## general job configuration
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
## configure object key store to recognize calo cells
from RecExConfig.ObjKeyStore import ObjKeyStore, objKeyStore
oks = ObjKeyStore()
oks.addStreamESD('CaloCellContainer', ['AllCalo'] )
## re-do topo clusters on EM scale
from CaloRec.CaloTopoClusterFlags import jobproperties
jobproperties.CaloTopoClusterFlags.doTopoClusterLocalCalib = False
jobproperties.CaloTopoClusterFlags.doCellWeightCalib = False
from CaloRec.CaloClusterTopoGetter import CaloClusterTopoGetter
CaloClusterTopoGetter()
topSequence.CaloTopoCluster.TopoCalibMoments.MomentsNames += ["ENG_CALIB_OUT_L"]
theApp.EvtMax = -1
#--------------------------------------------------------------
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
#--------------------------------------------------------------
MessageSvc = svcMgr.MessageSvc
MessageSvc.OutputLevel = Lvl.INFO
MessageSvc.infoLimit = 1000000
svcMgr.EventSelector.InputCollections = inFileName
include ("CaloLocalHadCalib/GetLCClassification_jobOptions.py")
print topSequence
ChronoStatSvc = Service( "ChronoStatSvc" )
ChronoStatSvc.ChronoDestinationCout = True
ChronoStatSvc.PrintUserTime = True
ChronoStatSvc.PrintSystemTime = True
ChronoStatSvc.PrintEllapsedTime = True
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
AthenaPoolCnvSvc.UseDetailChronoStat = True
# File: CaloLocalHadCalib/share/GetLCDeadMaterial_jobOptions.py
# Created: Febryary 2009, G. Pospelov
# Purpose: Run postprocessing algorithms for the special DeadMaterialTree
#
#
# set inFileName and outFileName to convert TProfile2D classification to pool
#
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
import IOVDbSvc.IOVDb
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
###############################################################################
#
###############################################################################
#mysfx="ATLAS.QGSPBERT.MC11.SINGLEPI.20062011.v1"
mysfx="mc10_7TeV.107410-107412.r2424.piplusminuszero.LCWeightsDM.16652.v2" # 15.06.2011
from CaloLocalHadCalib.CaloLocalHadCalibConf import GetLCDeadMaterial
GetLCDM = GetLCDeadMaterial("GetLCDM")
GetLCDM.HadDMCoeffInputFile = "CaloHadDMCoeff_init_v2.txt"
#GetLCDM.HadDMCoeffInputFile = "CaloLocalHadCoeff_output_ATLAS.QGSPBERT.MC10.SINGLEPI.21102010.v8.txt"
GetLCDM.HadDMCoeffOutputFile = "CaloLocalHadCoeff_output_"+mysfx+".txt"
GetLCDM.ReportProfiles = "report_CaloLocalHadCoeff_profiles_"+mysfx+".ps"
GetLCDM.ReportMinimization = "report_CaloLocalHadCoeff_minim_"+mysfx+".ps"
GetLCDM.ReportCheck = "report_CaloLocalHadCoeff_check_"+mysfx+".ps"
GetLCDM.DoFit = True
GetLCDM.DoMinimization = True
GetLCDM.DoPool = True
GetLCDM.DoCheck = True
# new feature to run on dmtrees from di-jets
#GetLCDM.NormalizationTypeForFit="const"
#GetLCDM.NormalizationTypeForMinim="const"
#GetLCDM.ClassificationType = "particleid"
GetLCDM.CorrectionKey="HadDMCoeff2"
#GetLCDM.InputRootFiles = ["/ltmp/gdp/LCWeights/def1/user.pospelov.mc10_7TeV.107412.singlepiminus.LCWeights.16016.def1._sum.dmc.root"]
#GetLCDM.InputRootFiles += ["/ltmp/gdp/LCWeights/def1/user.pospelov.mc10_7TeV.107412.singlepizero.LCWeights.16016.def1._sum.dmc.root"]
#GetLCDM.InputRootFiles += ["/ltmp/gdp/LCWeights/def1/user.pospelov.mc10_7TeV.107412.singlepiplus.LCWeights.16016.def1._sum.dmc.root"]
#GetLCDM.InputRootFiles = ["/ltmp/gdp/LCWeights/def2/user.pospelov.mc10_7TeV.107412.singlepiminus.LCWeights.16016.def2._sum.dmc.root"]
#GetLCDM.InputRootFiles += ["/ltmp/gdp/LCWeights/def2/user.pospelov.mc10_7TeV.107412.singlepizero.LCWeights.16016.def2._sum.dmc.root"]
#GetLCDM.InputRootFiles += ["/ltmp/gdp/LCWeights/def2/user.pospelov.mc10_7TeV.107412.singlepiplus.LCWeights.16016.def2b._sum.dmc.root"]
GetLCDM.InputRootFiles = ["/ltmp/gdp/aaa/user.pospelov.mc10_7TeV.107410.singlepiplus.r2424.LCWeights.16652.v1._sum.dmc.root"]
GetLCDM.InputRootFiles += ["/ltmp/gdp/aaa/user.pospelov.mc10_7TeV.107411.singlepiminus.r2424.LCWeights.16652.v1._sum.dmc.root"]
GetLCDM.InputRootFiles += ["/ltmp/gdp/aaa/user.pospelov.mc10_7TeV.107412.singlepizero.r2424.LCWeights.16652.v1._sum.dmc.root"]
#GetLCDM.MaxEventsPerFile = 500000
topSequence += GetLCDM
# Number of events to be processed (default is 10)
theApp.EvtMax = 1
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
from RegistrationServices.OutputConditionsAlg import OutputConditionsAlg
#myOCA=OutputConditionsAlg("myOCA","CaloLocalHadCoeff.dmcorr."+mysfx+".pool.root")
myOCA=OutputConditionsAlg("myOCA",mysfx+".pool.root")
# List of objects container type#key#foldername
myOCA.ObjectList = ["CaloLocalHadCoeff#HadDMCoeff2#/CALO/HadCalibration2/CaloDMCorr2"]
myOCA.WriteIOV = True
#myOCA.IOVTagList = ["CaloHadDMCorr2-CSC05-BERT" ]
#myOCA.IOVTagList = ["CaloHadDMCorr2-MC09V-BERT" ]
#myOCA.IOVTagList = ["CaloHadDMCorr2-GEO08-QGSP-BERT" ]
#myOCA.IOVTagList = ["CaloHadDMCorr2-GEO08-FTFP-BERT" ]
#myOCA.IOVTagList = ["CaloHadDMCorr2-GEO16-QGSP-BERT" ]
#myOCA.IOVTagList = ["CaloHadDMCorr2-GEO16-QGSP-BERT-Side" ]
myOCA.IOVTagList = ["CaloHadDMCorr2-GEO180100-QGSP-BERT" ]
myOCA.OutputLevel= DEBUG
AthenaPoolCnvSvc.PoolAttributes += [ "STREAM_MEMBER_WISE = '0'" ]
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
MessageSvc = Service( "MessageSvc" )
MessageSvc.OutputLevel = Lvl.INFO
DetectorStore = svcMgr.DetectorStore
DetectorStore.Dump = TRUE
DetectorStore.OutputLevel = Lvl.INFO
IOVDbSvc = Service( "IOVDbSvc" )
IOVDbSvc.OutputLevel = Lvl.INFO
svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myGeo18DmQGSPBERT200.db;dbname=OFLP200"
#svcMgr.IOVDbSvc.dbConnection = "sqlite://X;schema=myGeo08DmFTFPBERT200.db;dbname=OFLP200"
# This is the job options file to create DeadMaterialTree from ESD
###############################################################
#
# Job options file
#
# Authors: Gennady Pospelov
#==============================================================
#
# athena.py -c 'inFileName=["ESD.pool.root"];outFileName="dmtreeX.root"' CaloLocalHadCalib/deadmaterialtree.py
#
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
#inFileName=["/ltmp/gdp/aaa/ESD.081901._001996.pool.root.1"]
#outFileNameLCDM="dmtreeX.root"
## load POOL support
import AthenaPoolCnvSvc.ReadAthenaPool
## general job configuration
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
#from AthenaCommon.DetFlags import DetFlags
#DetFlags.ID_setOff()
#DetFlags.Muon_setOff()
#DetFlags.LVL1_setOff()
#DetDescrVersion = 'ATLAS-CSC-05-00-00'
DetDescrVersion = 'ATLAS-GEO-08-00-01'
include("RecExCond/AllDet_detDescr.py")
## configure object key store to recognize calo cells
from RecExConfig.ObjKeyStore import ObjKeyStore, objKeyStore
oks = ObjKeyStore()
oks.addStreamESD('CaloCellContainer', ['AllCalo'] )
## re-do topo clusters on EM scale
from CaloRec.CaloTopoClusterFlags import jobproperties
jobproperties.CaloTopoClusterFlags.doTopoClusterLocalCalib = False
jobproperties.CaloTopoClusterFlags.doCellWeightCalib = False
from CaloRec.CaloClusterTopoGetter import CaloClusterTopoGetter
CaloClusterTopoGetter()
topSequence.CaloTopoCluster.TopoCalibMoments.MomentsNames += ["ENG_CALIB_OUT_L","ENG_CALIB_FRAC_EM"]
topSequence.CaloTopoCluster.TopoCalibMoments.MatchDmType = 1 # 1=loose, 2=medium (default), 3=tight
theApp.EvtMax = -1
#--------------------------------------------------------------
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
#--------------------------------------------------------------
MessageSvc = svcMgr.MessageSvc
MessageSvc.OutputLevel = Lvl.INFO
MessageSvc.infoLimit = 1000000
if 'inFileName' in dir():
svcMgr.EventSelector.InputCollections = inFileName
include ("CaloLocalHadCalib/GetLCDeadMaterialTree_jobOptions.py")
print topSequence
ChronoStatSvc = Service( "ChronoStatSvc" )
ChronoStatSvc.ChronoDestinationCout = True
ChronoStatSvc.PrintUserTime = True
ChronoStatSvc.PrintSystemTime = True
ChronoStatSvc.PrintEllapsedTime = True
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
AthenaPoolCnvSvc.UseDetailChronoStat = True
# This is the job options file to create out-of-cluster histos from ESD
###############################################################
#
# Job options file
#
# Authors: Sven Menke
#==============================================================
DetDescrVersion = 'ATLAS-GEO-08-00-01'
include("RecExCond/AllDet_detDescr.py")
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
## load POOL support
import AthenaPoolCnvSvc.ReadAthenaPool
## general job configuration
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
## configure object key store to recognize calo cells
from RecExConfig.ObjKeyStore import ObjKeyStore, objKeyStore
oks = ObjKeyStore()
oks.addStreamESD('CaloCellContainer', ['AllCalo'] )
## re-do topo clusters on EM scale
from CaloRec.CaloTopoClusterFlags import jobproperties
jobproperties.CaloTopoClusterFlags.doTopoClusterLocalCalib = False
jobproperties.CaloTopoClusterFlags.doCellWeightCalib = False
from CaloRec.CaloClusterTopoGetter import CaloClusterTopoGetter
CaloClusterTopoGetter()
topSequence.CaloTopoCluster.TopoCalibMoments.MomentsNames += ["ENG_CALIB_OUT_L"]
theApp.EvtMax = -1
#--------------------------------------------------------------
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
#--------------------------------------------------------------
MessageSvc = svcMgr.MessageSvc
MessageSvc.OutputLevel = Lvl.INFO
MessageSvc.infoLimit = 1000000
svcMgr.EventSelector.InputCollections = inFileName
include ("CaloLocalHadCalib/GetLCOutOfCluster_jobOptions.py")
print topSequence
ChronoStatSvc = Service( "ChronoStatSvc" )
ChronoStatSvc.ChronoDestinationCout = True
ChronoStatSvc.PrintUserTime = True
ChronoStatSvc.PrintSystemTime = True
ChronoStatSvc.PrintEllapsedTime = True
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
AthenaPoolCnvSvc.UseDetailChronoStat = True
# This is the job options file to create weight histos from ESD
###############################################################
#
# Job options file
#
# Authors: Sven Menke
#==============================================================
DetDescrVersion = 'ATLAS-GEO-08-00-01'
include("RecExCond/AllDet_detDescr.py")
import AthenaCommon.Constants as Lvl
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AppMgr import theApp
## load POOL support
import AthenaPoolCnvSvc.ReadAthenaPool
## general job configuration
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
## configure object key store to recognize calo cells
from RecExConfig.ObjKeyStore import ObjKeyStore, objKeyStore
oks = ObjKeyStore()
oks.addStreamESD('CaloCellContainer', ['AllCalo'] )
## re-do topo clusters on EM scale
from CaloRec.CaloTopoClusterFlags import jobproperties
jobproperties.CaloTopoClusterFlags.doTopoClusterLocalCalib = False
jobproperties.CaloTopoClusterFlags.doCellWeightCalib = False
from CaloRec.CaloClusterTopoGetter import CaloClusterTopoGetter
CaloClusterTopoGetter()
topSequence.CaloTopoCluster.TopoCalibMoments.MomentsNames += ["ENG_CALIB_OUT_L"]
theApp.EvtMax = -1
#--------------------------------------------------------------
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
#--------------------------------------------------------------
MessageSvc = svcMgr.MessageSvc
MessageSvc.OutputLevel = Lvl.INFO
MessageSvc.infoLimit = 1000000
svcMgr.EventSelector.InputCollections = inFileName
include ("CaloLocalHadCalib/GetLCWeights_jobOptions.py")
print topSequence
ChronoStatSvc = Service( "ChronoStatSvc" )
ChronoStatSvc.ChronoDestinationCout = True
ChronoStatSvc.PrintUserTime = True
ChronoStatSvc.PrintSystemTime = True
ChronoStatSvc.PrintEllapsedTime = True
AthenaPoolCnvSvc = Service( "AthenaPoolCnvSvc" )
AthenaPoolCnvSvc.UseDetailChronoStat = True
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment