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

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

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;}
    }
Edit | Attach | Watch | Print version | History: r21 < r20 < r19 < r18 < r17 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r21 - 2017-07-24 - RustemOspanov
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback