diff --git a/Reconstruction/egamma/egammaValidation/scripts/EgammaARTmonitoring_plotsMaker.py b/Reconstruction/egamma/egammaValidation/scripts/EgammaARTmonitoring_plotsMaker.py index 4e310673feea3f8e7c7a8c6649593033886034db..fbf065d2f9e4d95b4d2005d641dd299501a745cd 100755 --- a/Reconstruction/egamma/egammaValidation/scripts/EgammaARTmonitoring_plotsMaker.py +++ b/Reconstruction/egamma/egammaValidation/scripts/EgammaARTmonitoring_plotsMaker.py @@ -13,8 +13,8 @@ def GetKeyNames(f0, dir = ""): f0.cd(dir) return [key.GetName() for key in gDirectory.GetListOfKeys()] -f1 = TFile(sys.argv[1]) -f2 = TFile(sys.argv[2]) +f_baseline = TFile(sys.argv[1]) +f_nightly = TFile(sys.argv[2]) particleType = sys.argv[3] @@ -22,8 +22,8 @@ fO = TFile("BN_ComparisonPlots_"+particleType+".hist.root", "RECREATE") gROOT.SetBatch(kTRUE) gStyle.SetOptStat(0) -for folder in GetKeyNames(f1): - for histo in GetKeyNames(f1,folder): +for folder in GetKeyNames(f_nightly): + for histo in GetKeyNames(f_nightly,folder): c1 = TCanvas() @@ -48,21 +48,21 @@ for folder in GetKeyNames(f1): c1.Update() - h_Base = f1.Get(folder+'/'+histo) + h_Base = f_baseline.Get(folder+'/'+histo) h_Base.SetLineColor(4) h_Base.SetLineWidth(2) h_Base.GetXaxis().SetLabelOffset(1.20) h_Base.GetYaxis().SetTitleSize(0.045) h_Base.GetYaxis().SetTitleOffset(0.95) - h_Night = f2.Get(folder+'/'+histo) + h_Night = f_nightly.Get(folder+'/'+histo) h_Night.SetMarkerStyle(8) h_Night.SetMarkerSize(0.5) mainPad.cd() - h_Base.Draw() - h_Night.Draw("same p") + if not "2D" in histo: h_Base.Draw() + h_Night.Draw("same p" if not "2D" in histo else "colz") c1.Update() @@ -73,7 +73,7 @@ for folder in GetKeyNames(f1): leg.SetNColumns(2) leg.SetFillStyle(0) leg.SetBorderSize(0) - leg.AddEntry(h_Base , "Baseline", "l") + if not "2D" in histo: leg.AddEntry(h_Base , "Baseline", "l") leg.AddEntry(h_Night, "Nightly" , "p") leg.Draw() @@ -81,26 +81,28 @@ for folder in GetKeyNames(f1): ratioPad.cd() - h1clone = h_Night.Clone() - h1clone.Sumw2() - h1clone.SetStats(0) - h1clone.Divide(h_Base) - h1clone.SetMarkerColor(1) - h1clone.SetMarkerStyle(20) - if "Efficiency" in folder: - h1clone.GetYaxis().SetRangeUser(h1clone.GetMinimum()*0.7,h1clone.GetMaximum()*1.3) - gStyle.SetOptStat(0) - h1clone.GetXaxis().SetLabelSize(0.10) - h1clone.GetXaxis().SetTitleSize(0.17) - h1clone.GetYaxis().SetLabelSize(0.10) - h1clone.GetYaxis().SetRangeUser(0.75, 1.25) - h1clone.GetYaxis().SetTitle("Ratio") - h1clone.GetYaxis().CenterTitle(1) - h1clone.GetYaxis().SetTitleSize(0.15) - h1clone.GetYaxis().SetTitleOffset(0.3) - h1clone.GetYaxis().SetNdivisions(505) - - h1clone.Draw("p") + if not "2D" in histo: + h1clone = h_Night.Clone() + h1clone.Sumw2() + h1clone.SetStats(0) + h1clone.Divide(h_Base) + h1clone.SetMarkerColor(1) + h1clone.SetMarkerStyle(20) + if "Efficiency" in folder: + h1clone.GetYaxis().SetRangeUser(h1clone.GetMinimum()*0.7,h1clone.GetMaximum()*1.3) + gStyle.SetOptStat(0) + h1clone.GetXaxis().SetLabelSize(0.10) + h1clone.GetXaxis().SetTitleSize(0.17) + h1clone.GetYaxis().SetLabelSize(0.10) + h1clone.GetYaxis().SetRangeUser(0.75, 1.25) + h1clone.GetYaxis().SetTitle("Ratio") + h1clone.GetYaxis().CenterTitle(1) + h1clone.GetYaxis().SetTitleSize(0.15) + h1clone.GetYaxis().SetTitleOffset(0.3) + h1clone.GetYaxis().SetNdivisions(505) + + h1clone.Draw("p") + c1.Update() c1.SaveAs(folder+'_'+var_name+".png" ) diff --git a/Reconstruction/egamma/egammaValidation/share/egamma_art_checker_joboptions.py b/Reconstruction/egamma/egammaValidation/share/egamma_art_checker_joboptions.py index 147b1bb99bc042dfd7d9004904d445b2e192f5ed..624bf5650be1b963f2a5d746999bc91d2e4a0b5c 100755 --- a/Reconstruction/egamma/egammaValidation/share/egamma_art_checker_joboptions.py +++ b/Reconstruction/egamma/egammaValidation/share/egamma_art_checker_joboptions.py @@ -6,8 +6,7 @@ from AthenaCommon.AppMgr import theApp from AthenaCommon.AppMgr import ServiceMgr as svcMgr from AthenaCommon.AlgSequence import AlgSequence from GaudiSvc.GaudiSvcConf import THistSvc -import AthenaRootComps.ReadAthenaxAODHybrid - +import AthenaPoolCnvSvc.ReadAthenaPool from MCTruthClassifier.MCTruthClassifierConf import MCTruthClassifier MCClassifier = MCTruthClassifier( name = "MCTruthClassifier", diff --git a/Reconstruction/egamma/egammaValidation/src/ShowerShapesHistograms.cxx b/Reconstruction/egamma/egammaValidation/src/ShowerShapesHistograms.cxx index 1bedb98af63c4feadeeb1ff2cad69dffecd2ca3f..fb35ad1543fa80133698a71f8beb7d01bb6f8023 100644 --- a/Reconstruction/egamma/egammaValidation/src/ShowerShapesHistograms.cxx +++ b/Reconstruction/egamma/egammaValidation/src/ShowerShapesHistograms.cxx @@ -22,9 +22,16 @@ StatusCode ShowerShapesHistograms::initializePlots() { histoMap["fside"] = (new TH1D(Form("%s_%s",m_name.c_str(),"fside" ), ";f_{side}; f_{side} Events" , 350, 0. , 3.5 )); histoMap["wtots1"] = (new TH1D(Form("%s_%s",m_name.c_str(),"wtots1" ), ";w_{s, tot}; w_{s, tot} Events" , 100, 0. , 10. )); histoMap["ws3"] = (new TH1D(Form("%s_%s",m_name.c_str(),"ws3" ), ";w_{s, 3}; w_{s, 3} Events" , 100, 0. , 1. )); - histoMap["lateral"] = (new TH1D(Form("%s_%s",m_name.c_str(),"lateral"), ";Lateral of seed; Events", 10, 0, 1)); + histoMap["lateral"] = (new TH1D(Form("%s_%s",m_name.c_str(),"lateral"), ";Lateral of seed; Events", 50, 0, 1)); + histoMap["second_R"] = (new TH1D(Form("%s_%s",m_name.c_str(),"second_R"), ";Second R; Events", 150,0. , 15000.)); + histoMap["EMFrac"] = (new TH1D(Form("%s_%s",m_name.c_str(),"EMFrac"), ";EMFrac; Events", 50, 0, 1)); + histo2DMap["lateral_second_R_2D"] = (new TH2D(Form("%s_%s",m_name.c_str(),"lateral_second_R_2D"), ";Lateral of seed; Second R", 50, 0, 1, 150,0.,15000.)); ATH_CHECK(m_rootHistSvc->regHist(m_folder+"lateral", histoMap["lateral"])); + ATH_CHECK(m_rootHistSvc->regHist(m_folder+"second_R", histoMap["second_R"])); + ATH_CHECK(m_rootHistSvc->regHist(m_folder+"EMFrac", histoMap["EMFrac"])); + ATH_CHECK(m_rootHistSvc->regHist(m_folder+"lateral_second_R_2D", histo2DMap["lateral_second_R_2D"])); + ATH_CHECK(m_rootHistSvc->regHist(m_folder+"rhad", histoMap["rhad"])); ATH_CHECK(m_rootHistSvc->regHist(m_folder+"hadleak", histoMap["hadleak"])); ATH_CHECK(m_rootHistSvc->regHist(m_folder+"reta", histoMap["reta"])); @@ -56,13 +63,24 @@ void ShowerShapesHistograms::fill(const xAOD::Egamma& egamma) { } double lateral(0.); - + double emfrac(0.); + double second_r(0.); const std::vector<const xAOD::CaloCluster*> topoclusters = xAOD::EgammaHelpers::getAssociatedTopoClusters(egamma.caloCluster()); if (topoclusters.at(0)->retrieveMoment(xAOD::CaloCluster::LATERAL,lateral)) { histoMap["lateral"]->Fill(lateral); } + if (topoclusters.at(0)->retrieveMoment(xAOD::CaloCluster::ENG_FRAC_EM ,emfrac)) { + histoMap["EMFrac"]->Fill(emfrac); + } + + if (topoclusters.at(0)->retrieveMoment(xAOD::CaloCluster::SECOND_R,second_r)) { + histoMap["second_R"]->Fill(second_r); + } + + histo2DMap["lateral_second_R_2D"]->Fill(lateral,second_r); + if(egamma.showerShapeValue(rhad, xAOD::EgammaParameters::Rhad)) histoMap["rhad"]->Fill(rhad); if(egamma.showerShapeValue(Reta, xAOD::EgammaParameters::Reta)) histoMap["reta"]->Fill(Reta); if(egamma.showerShapeValue(Rphi, xAOD::EgammaParameters::Rphi)) histoMap["rphi"]->Fill(Rphi); diff --git a/Reconstruction/egamma/egammaValidation/src/ShowerShapesHistograms.h b/Reconstruction/egamma/egammaValidation/src/ShowerShapesHistograms.h index c06bfdd6e887b708fd2c5f3c171c22ab1e80061b..f2bf921be3f5f09eb696c21e94d2536e6c919cda 100644 --- a/Reconstruction/egamma/egammaValidation/src/ShowerShapesHistograms.h +++ b/Reconstruction/egamma/egammaValidation/src/ShowerShapesHistograms.h @@ -11,6 +11,7 @@ #include "EgammaAnalysisInterfaces/IAsgElectronLikelihoodTool.h" #include "TH1.h" +#include "TH2.h" namespace egammaMonitoring{ @@ -29,6 +30,8 @@ namespace egammaMonitoring{ m_rootHistSvc(rootHistSvc) {} std::map<std::string, TH1D* > histoMap; + std::map<std::string, TH2D* > histo2DMap; + StatusCode initializePlots(); void fill(const xAOD::Egamma& egamma); diff --git a/Reconstruction/egamma/egammaValidation/test/test_electron.sh b/Reconstruction/egamma/egammaValidation/test/test_electron.sh index 1a507d4ecf280719ef9af959ef842f61f5eb3a10..593375b162dc16b3e5bf4000fbc7c4301216dfa5 100755 --- a/Reconstruction/egamma/egammaValidation/test/test_electron.sh +++ b/Reconstruction/egamma/egammaValidation/test/test_electron.sh @@ -29,14 +29,13 @@ case $ArtProcess in unset ATHENA_NUM_PROC unset ATHENA_PROC_NUMBER - AODMerge_tf.py --inputAODFile=art_core_*/Nightly_AOD_electron.pool.root --outputAOD_MRGFile=Nightly_AOD_electron.pool.root + AODMerge_tf.py --inputAODFile=art_core_*/Nightly_AOD_electron.pool.root --outputAOD_MRGFile=Nightly_AOD_electron.pool.root --postExec "all:CILMergeAOD.removeItem(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList\");CILMergeAOD.add(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList.ENG_CALIB_TOT.ENG_CALIB_OUT_L.ENG_CALIB_EMB0.ENG_CALIB_EME0.ENG_CALIB_TILEG3.ENG_CALIB_DEAD_TOT.ENG_CALIB_DEAD_EMB0.ENG_CALIB_DEAD_TILE0.ENG_CALIB_DEAD_TILEG3.ENG_CALIB_DEAD_EME0.ENG_CALIB_DEAD_HEC0.ENG_CALIB_DEAD_FCAL.ENG_CALIB_DEAD_LEAKAGE.ENG_CALIB_DEAD_UNCLASS.ENG_CALIB_FRAC_EM.ENG_CALIB_FRAC_HAD.ENG_CALIB_FRAC_REST.LATERAL.ENG_FRAC_EM\");StreamAOD.ItemList=CILMergeAOD()" echo "art-result: $? AODMerge" set +e checkFile.py Nightly_AOD_electron.pool.root > checkFile_Nightly_electron.txt - checkxAOD.py Nightly_AOD_electron.pool.root > checkxAOD_Nightly_electron.txt echo "art-result: $? checks_files" @@ -65,7 +64,7 @@ case $ArtProcess in unset ATHENA_NUM_PROC unset ATHENA_PROC_NUMBER - Reco_tf.py --inputRDOFile=$x --outputAODFile=Nightly_AOD_electron.pool.root --maxEvents=2000 --autoConfiguration="everything" --preExec="from ParticleBuilderOptions.AODFlags import AODFlags; AODFlags.ThinGeantTruth.set_Value_and_Lock(False); rec.doTrigger=False; rec.doMuon=False; rec.doTau=False ; from RecExConfig.RecAlgsFlags import recAlgs; recAlgs.doMuonSpShower=False ; rec.doBTagging=False ; recAlgs.doEFlow=False ; recAlgs.doEFlowJet=False ; recAlgs.doMissingET=False ; recAlgs.doMissingETSig=False ; from JetRec.JetRecFlags import jetFlags ; jetFlags.Enabled=False ; " + Reco_tf.py --inputRDOFile=$x --outputAODFile=Nightly_AOD_electron.pool.root --maxEvents=2000 --autoConfiguration="everything" --preExec="from ParticleBuilderOptions.AODFlags import AODFlags; AODFlags.ThinGeantTruth.set_Value_and_Lock(False); rec.doTrigger=False; rec.doMuon=False; rec.doTau=False ; from RecExConfig.RecAlgsFlags import recAlgs; recAlgs.doMuonSpShower=False ; rec.doBTagging=False ; recAlgs.doEFlow=False ; recAlgs.doEFlowJet=False ; recAlgs.doMissingET=False ; recAlgs.doMissingETSig=False ; from JetRec.JetRecFlags import jetFlags ; jetFlags.Enabled=False ; " --postExec "ESDtoAOD:CILMergeAOD.removeItem(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList\");CILMergeAOD.add(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.LATERAL.ENG_FRAC_EM.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList.ENG_CALIB_TOT.ENG_CALIB_OUT_L.ENG_CALIB_EMB0.ENG_CALIB_EME0.ENG_CALIB_TILEG3.ENG_CALIB_DEAD_TOT.ENG_CALIB_DEAD_EMB0.ENG_CALIB_DEAD_TILE0.ENG_CALIB_DEAD_TILEG3.ENG_CALIB_DEAD_EME0.ENG_CALIB_DEAD_HEC0.ENG_CALIB_DEAD_FCAL.ENG_CALIB_DEAD_LEAKAGE.ENG_CALIB_DEAD_UNCLASS.ENG_CALIB_FRAC_EM.ENG_CALIB_FRAC_HAD.ENG_CALIB_FRAC_REST\");StreamAOD.ItemList=CILMergeAOD()" "POOLMergeAthenaMPAOD0:CILMergeAOD.removeItem(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList\");CILMergeAOD.add(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.LATERAL.ENG_FRAC_EM.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList.ENG_CALIB_TOT.ENG_CALIB_OUT_L.ENG_CALIB_EMB0.ENG_CALIB_EME0.ENG_CALIB_TILEG3.ENG_CALIB_DEAD_TOT.ENG_CALIB_DEAD_EMB0.ENG_CALIB_DEAD_TILE0.ENG_CALIB_DEAD_TILEG3.ENG_CALIB_DEAD_EME0.ENG_CALIB_DEAD_HEC0.ENG_CALIB_DEAD_FCAL.ENG_CALIB_DEAD_LEAKAGE.ENG_CALIB_DEAD_UNCLASS.ENG_CALIB_FRAC_EM.ENG_CALIB_FRAC_HAD.ENG_CALIB_FRAC_REST\");StreamAOD.ItemList=CILMergeAOD()" echo "art-result: $? reconstruction" diff --git a/Reconstruction/egamma/egammaValidation/test/test_gamma.sh b/Reconstruction/egamma/egammaValidation/test/test_gamma.sh index da84d21e05cfa06011463546a4c7e0a2020beab0..260fb0e9cdc6e2e3f3686e287499e7e4e4496aac 100755 --- a/Reconstruction/egamma/egammaValidation/test/test_gamma.sh +++ b/Reconstruction/egamma/egammaValidation/test/test_gamma.sh @@ -11,6 +11,7 @@ # art-output: *.hist.root # art-output: *.txt # art-output: *.png +# art-output: log.* echo "ArtProcess: $ArtProcess" @@ -29,14 +30,13 @@ case $ArtProcess in unset ATHENA_NUM_PROC unset ATHENA_PROC_NUMBER - AODMerge_tf.py --inputAODFile=art_core_*/Nightly_AOD_gamma.pool.root --outputAOD_MRGFile=Nightly_AOD_gamma.pool.root + AODMerge_tf.py --inputAODFile=art_core_*/Nightly_AOD_gamma.pool.root --outputAOD_MRGFile=Nightly_AOD_gamma.pool.root --postExec "all:CILMergeAOD.removeItem(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList\");CILMergeAOD.add(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList.ENG_CALIB_TOT.ENG_CALIB_OUT_L.ENG_CALIB_EMB0.ENG_CALIB_EME0.ENG_CALIB_TILEG3.ENG_CALIB_DEAD_TOT.ENG_CALIB_DEAD_EMB0.ENG_CALIB_DEAD_TILE0.ENG_CALIB_DEAD_TILEG3.ENG_CALIB_DEAD_EME0.ENG_CALIB_DEAD_HEC0.ENG_CALIB_DEAD_FCAL.ENG_CALIB_DEAD_LEAKAGE.ENG_CALIB_DEAD_UNCLASS.ENG_CALIB_FRAC_EM.ENG_CALIB_FRAC_HAD.ENG_CALIB_FRAC_REST.LATERAL.ENG_FRAC_EM\");StreamAOD.ItemList=CILMergeAOD()" echo "art-result: $? AODMerge" set +e checkFile.py Nightly_AOD_gamma.pool.root > checkFile_Nightly_gamma.txt - checkxAOD.py Nightly_AOD_gamma.pool.root > checkxAOD_Nightly_gamma.txt echo "art-result: $? checks_files" @@ -65,7 +65,7 @@ case $ArtProcess in unset ATHENA_NUM_PROC unset ATHENA_PROC_NUMBER - Reco_tf.py --inputRDOFile=$x --outputAODFile=Nightly_AOD_gamma.pool.root --maxEvents=2000 --autoConfiguration="everything" --preExec="from ParticleBuilderOptions.AODFlags import AODFlags; AODFlags.ThinGeantTruth.set_Value_and_Lock(False); rec.doTrigger=False; rec.doMuon=False; rec.doTau=False ; from RecExConfig.RecAlgsFlags import recAlgs; recAlgs.doMuonSpShower=False ; rec.doBTagging=False ; recAlgs.doEFlow=False ; recAlgs.doEFlowJet=False ; recAlgs.doMissingET=False ; recAlgs.doMissingETSig=False ; from JetRec.JetRecFlags import jetFlags ; jetFlags.Enabled=False ; " + Reco_tf.py --inputRDOFile=$x --outputAODFile=Nightly_AOD_gamma.pool.root --maxEvents=2000 --autoConfiguration="everything" --preExec="from ParticleBuilderOptions.AODFlags import AODFlags; AODFlags.ThinGeantTruth.set_Value_and_Lock(False); rec.doTrigger=False; rec.doMuon=False; rec.doTau=False ; from RecExConfig.RecAlgsFlags import recAlgs; recAlgs.doMuonSpShower=False ; rec.doBTagging=False ; recAlgs.doEFlow=False ; recAlgs.doEFlowJet=False ; recAlgs.doMissingET=False ; recAlgs.doMissingETSig=False ; from JetRec.JetRecFlags import jetFlags ; jetFlags.Enabled=False ; " --postExec "ESDtoAOD:CILMergeAOD.removeItem(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList\");CILMergeAOD.add(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.LATERAL.ENG_FRAC_EM.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList.ENG_CALIB_TOT.ENG_CALIB_OUT_L.ENG_CALIB_EMB0.ENG_CALIB_EME0.ENG_CALIB_TILEG3.ENG_CALIB_DEAD_TOT.ENG_CALIB_DEAD_EMB0.ENG_CALIB_DEAD_TILE0.ENG_CALIB_DEAD_TILEG3.ENG_CALIB_DEAD_EME0.ENG_CALIB_DEAD_HEC0.ENG_CALIB_DEAD_FCAL.ENG_CALIB_DEAD_LEAKAGE.ENG_CALIB_DEAD_UNCLASS.ENG_CALIB_FRAC_EM.ENG_CALIB_FRAC_HAD.ENG_CALIB_FRAC_REST\");StreamAOD.ItemList=CILMergeAOD()" "POOLMergeAthenaMPAOD0:CILMergeAOD.removeItem(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList\");CILMergeAOD.add(\"xAOD::CaloClusterAuxContainer#CaloCalTopoClustersAux.SECOND_R.LATERAL.ENG_FRAC_EM.SECOND_LAMBDA.CENTER_MAG.CENTER_LAMBDA.FIRST_ENG_DENS.ENG_FRAC_MAX.ISOLATION.ENG_BAD_CELLS.N_BAD_CELLS.BADLARQ_FRAC.ENG_POS.AVG_LAR_Q.AVG_TILE_Q.EM_PROBABILITY.BadChannelList.ENG_CALIB_TOT.ENG_CALIB_OUT_L.ENG_CALIB_EMB0.ENG_CALIB_EME0.ENG_CALIB_TILEG3.ENG_CALIB_DEAD_TOT.ENG_CALIB_DEAD_EMB0.ENG_CALIB_DEAD_TILE0.ENG_CALIB_DEAD_TILEG3.ENG_CALIB_DEAD_EME0.ENG_CALIB_DEAD_HEC0.ENG_CALIB_DEAD_FCAL.ENG_CALIB_DEAD_LEAKAGE.ENG_CALIB_DEAD_UNCLASS.ENG_CALIB_FRAC_EM.ENG_CALIB_FRAC_HAD.ENG_CALIB_FRAC_REST\");StreamAOD.ItemList=CILMergeAOD()" echo "art-result: $? reconstruction"