Introduction

• This page describes how to run a complete analysis from DxAOD using Physics Analysis Ntuple framework (PhysicsAnp) - work in progress

References for xAOD tools

• SoftwareTutorialxAODAnalysisInAthena - xAOD analysis with Athena
• Event
  • ExtendedPileupReweighting
  • AsgAnalysisInterfaces/AsgAnalysisInterfaces
  • AnalysisMetadata
  • TopAnalysis/trunk/util/top-xaod.cxx
  • AthAnalysisBase - TWiki for xAOD with athena (including metadata)
  • AssociationUtils
• Electrons
  • ElectronEfficiencyRun2
• Muons
  • MCPAnalysisGuidelinesMC15
  • Muon trigger TWiki
• Jets
  • JetEtmissRecommendations2016
• Trigger
  • https://twiki.cern.ch/twiki/bin/view/Atlas/TriggerAnalysisTools
  • Trigger/TrigAnalysis/TrigDecisionTool/TrigDecisionTool/TrigDecisionTool.h

Grid production

PhysicsAth code

• Set up and compile code:

   $ mkdir -p ~/testarea/AtlasAth
   $ cd ~/testarea/AtlasAth
   $ svn co svn+ssh://svn.cern.ch/reps/manc/PhysicsNtuple/PhysicsAth/trunk PhysicsNtuple/PhysicsAth
   $ ln -s PhysicsNtuple/PhysicsAth/macros/setup/setup_atlas_release_cvmfs.sh .
   $ source setup_atlas_release_cvmfs.sh
   $ cd PhysicsNtuple/PhysicsAth/cmt
   $ cmt config; make -j4

• Run local job:

   $ athena share/procReadEGAM.py -c "inputDir='data16_13TeV.00300655.physics_Main.merge.DAOD_EGAM1.f708_m1606_p2667';outFile='out.root';EvtMax=1000"

• Submit grid job:

   $ localSetupPandaClient
   $ pathena $TestArea/PhysicsNtuple/PhysicsAth/share/procReadEGAM.py --inDS=mc15_13TeV.410000.PowhegPythiaEvtGen_P2012_ttbar_hdamp172p5_nonallhad.merge.DAOD_EGAM1.e3698_s2608_s2183_r7725_r7676_p2666 --outDS=user.rustem.mc15_13TeV.410000.PowhegPythiaEvtGen_P2012_ttbar_hdamp172p5_nonallhad.EGAM1.r7725_egam_prodv03 --destSE=UKI-NORTHGRID-MAN-HEP_LOCALGROUPDISK

Grid commands:

• file_list_rucio.py
  • prepare list of datasets for grid production
• subGrid.py
  • prepare shell scripts for submitting grid pathena jobs
  • prepare text files with names of output datasets

Useful links

• ttH Run 2 framework based on TopPhys
• Introduction to the TopxAOD framework, aka Top WG framework
• TopPhys xAOD packages
• xAOD data classes
• ATLAS code browser

Useful datasets

mc15_13TeV.343365.aMcAtNloPythia8EvtGen_A14_NNPDF23_NNPDF30ME_ttH125_dilep.merge.AOD.e4706_s2726_r7772_r7676
mc15_13TeV.343366.aMcAtNloPythia8EvtGen_A14_NNPDF23_NNPDF30ME_ttH125_semilep.merge.AOD.e4706_s2726_r7772_r7676
mc15_13TeV.343367.aMcAtNloPythia8EvtGen_A14_NNPDF23_NNPDF30ME_ttH125_allhad.merge.AOD.e4706_s2726_r7772_r7676

mc15_13TeV.410155.aMcAtNloPythia8EvtGen_MEN30NLO_A14N23LO_ttW.merge.AOD.e5070_s2726_r7772_r7676
mc15_13TeV.410156.aMcAtNloPythia8EvtGen_MEN30NLO_A14N23LO_ttZnunu.merge.AOD.e5070_s2726_r7772_r7676
mc15_13TeV.410157.aMcAtNloPythia8EvtGen_MEN30NLO_A14N23LO_ttZqq.merge.AOD.e5070_s2726_r7772_r7676

mc15_13TeV.410218.aMcAtNloPythia8EvtGen_MEN30NLO_A14N23LO_ttee.merge.AOD.e5070_s2726_r7772_r7676
mc15_13TeV.410219.aMcAtNloPythia8EvtGen_MEN30NLO_A14N23LO_ttmumu.merge.AOD.e5070_s2726_r7772_r7676
mc15_13TeV.410220.aMcAtNloPythia8EvtGen_MEN30NLO_A14N23LO_tttautau.merge.AOD.e5070_s2726_r7772_r7676

Truth jet association

bool isHS = false;
if((jet->isAvailable< ElementLink<xAOD::JetContainer> >("GhostTruthAssociationLink") &&
        jet->auxdata< ElementLink<xAOD::JetContainer> >("GhostTruthAssociationLink").isValid()) )
    {
        const xAOD::Jet* tjet = * jet->auxdata< ElementLink<xAOD::JetContainer> >("GhostTruthAssociationLink");
        isHS = tjet->pt()>10000;
        if(isHS) { std::cout << "found HS jet " << std::endl;}
    }