MC Generation for ttH

20 July 2015

Generating POWHEG+Herwig++

Jobs fail with
IncludeError: include file PowhegControl/PowhegControl_ttH_Common.py can not be found

--> Wrong release version. With 19.2.3.7, there is no PowhegControl_ttH_Common.py in $AtlasArea/InstallArea/jobOptions/PowhegControl (/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.3/AtlasProduction/19.2.3.7/InstallArea/jobOptions/PowhegControl). It is present in 19.2.4.1. --> Note: Shouldn't I be using release 20.x for MC15? But that does not seem to have PowhegControl at all.

3 June 2015

Running on Grid

aMC@NLO & Herwig++

Using the LHE inputs used for the latest official ttH MC (group.phys-gener:group.phys-gener.aMCatNLOHerwigpp.v2_2_1.341270.ttH_inclusive_13TeV.mc15_v2/). This submits 40 jobs with 3125 events each (the odd number because the LHE inputs have 12500 events).

asetup 19.2.3.7,here
pathena --cmtConfig=x86_64-slc6-gcc47-opt --individualOutDS --inDS=group.phys-gener:group.phys-gener.aMCatNLOHerwigpp.v2_2_1.341270.ttH_inclusive_13TeV.mc15_v2/ --nFiles=10 --nSkipFiles=2 --trf "Generate_tf.py --ecmEnergy=13000 --runNumber=341270 --maxEvents=3125 --randomSeed=%RNDM:110 --jobConfig=MC15.341270.aMcAtNloHerwigppEvtGen_UEEE5_CTEQ6L1_CT10ME_ttH125_inc_semil.py --inputGeneratorFile=%IN --outputEVNTFile=%OUT.evgen.pool.root --firstEvent=%SKIPEVENTS" --outDS user.jmansour.MC15.aMcAtNloHerwigpp.test.00014 --disableAutoRetry --noBuild --nEventsPerJob=3125 --nEventsPerFile=12500

POWHEG & Pythia6

Submits 20 jobs with 5000 events each (100,000 events total). Note that if we use more events per job, they will crash since there is a memory leak somewhere. I believe the limit is around 6000 events, depending on the memory of the grid nodes.

asetup 19.2.4.1,here
pathena --individualOutDS --trf "Generate_tf.py --ecmEnergy 13000 --evgenJobOpts=MC15JobOpts-00-00-65_v2.tar.gz --runNumber 500001 --maxEvents 5000 --jobConfig MC15.500001.PowhegPythiaEvtGen_P2012_ttH_semilep_jason.py --outputEVNTFile test.EVNT.pool.root --steering=afterburn --randomSeed=%RNDM:210 --outputEVNTFile=%OUT.evgen.pool.root" --outDS user.jmansour.MC15.PowhegPythia.test.00028 --nJobs=20 --nEventsPerJob=5000 --disableAutoRetry --noBuild

If you get a "NameError", there was a typo in skeleton.ABtoEVGEN.py. Copy /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.3/AtlasProduction/19.2.3.7/InstallArea/jobOptions/EvgenJobTransforms/skeleton.ABtoEVGEN.py to InstallArea /jobOptions/EvgenJobTransforms/skeleton.ABtoEVGEN.py (relative to where you did asetup and are running pathena from) and fix the typo (a spurious "l" at the beginning of a line).

JasonDhiaMansour - 2015-06-03

27 May 2015

Generating EVNT:

setupATLAS
asetup 19.2.X.Y-VAL,r6,hereGenerate_tf.py --ecmEnergy 13000 --runNumber 500001 \
 --firstEvent 1 --maxEvents 5000 --randomSeed 678 \
  --jobConfig MC15.500001.PowhegPythiaEvtGen_P2012_ttH_semilep_T0.py \
  --outputEVNTFile test.EVNT.pool.root --steering=afterburn

MC15.500001.PowhegPythiaEvtGen_P2012_ttH_semilep_T0.py

include ('PowhegControl_tth_semilep_jason.py')

PowhegControl_tth_semilep_jason.py

from AthenaCommon.Logging import logging  # loads logger
msg = logging.getLogger( "JobOptions" )

evgenConfig.description    = "POWHEG+Pythia6 ttH production with Powheg hdamp equal top mass, TESTING, Perugia 2012 tune, at least one lepton filter, had top"
evgenConfig.keywords    = [ 'SM', 'top', 'lepton', 'Higgs']
evgenConfig.contact     = ['jason.mansour@cern.ch', 'lily.asquith@cern.ch']
# MC15: everything goes inside this if statement
if runArgs.trfSubstepName == 'generate' :
   include('PowhegControl/PowhegControl_ttH_Common.py')
   PowhegConfig.topdecaymode = 22222
   # bbbar
   PowhegConfig.hdecaymode = 5
   PowhegConfig.hdamp        = 172.5
   #msg.info("nEvents %s", PowhegConfig.nEvents)
   PowhegConfig.nEvents    *= 3
   #PowhegConfig.bornsuppfact = supp
   #msg.info("nEvents %s", PowhegConfig.nEvents)
   #msg.info(PowhegConfig.configurable_parameters)
   #import sys
   #sys.exit(1)
   PowhegConfig.generateRunCard()
   PowhegConfig.generateEvents()
   include('MC15JobOptions/PowhegPythia_Perugia2012_Common.py')
   include('MC15JobOptions/Pythia_Tauola.py')
   include('MC15JobOptions/Pythia_Photos.py')
   include('MC15JobOptions/TTbarWToLeptonFilter.py')
   filtSeq.TTbarWToLeptonFilter.NumLeptons = -1
   filtSeq.TTbarWToLeptonFilter.Ptcut = 0.
   include('MC15JobOptions/BoostedHadTopAndTopPair.py')
   # filtSeq.BoostedHadTopAndTopPair.tHadPtMin  = top_pt_min*1000
   # filtSeq.BoostedHadTopAndTopPair.tHadPtMax  = top_pt_max*1000
#  Run EvtGen as afterburner
include('MC15JobOptions/Pythia_Powheg_EvtGen.py')

This might fail with

18:34:11   Failed to read LHEF event information in USELHEF,
18:34:11   assume end of file has been reached.
18:34:11 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
18:34:11 !!! LHEF TERMINATES NORMALLY: NO MORE EVENTS IN FILE !!!
18:34:11 !!! PLEASE IGNORE ANY ATHENA ERROR MESSAGES LIKE !!!
18:34:11 !!! AthenaEventLoopMgr  ERROR Terminating event processing loop due to errors !!!
18:34:11 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
18:34:11 Pythia              FATAL ../src/GenModule.cxx:51 (StatusCode GenModule::execute()): code 0: this->callGenerator()

This is because not enough LHE events are generated. By default it generates numEvents*1.1 events, and subsequent steps filter them out. If the efficiency is too low, we get this error message. Solution: Add something like

PowhegConfig.nEvents    *= 3

to the jobOptions.

Convert EVNT to XAOD:

setupATLAS
asetup AtlasDerivation,20.1.4.2,gcc48,here
Reco_tf.py --inputEVNTFile test.EVNT.pool.root --outputDAODFile test.pool.root --reductionConf TRUTH1

Checking XAOD:

cd /afs/cern.ch/user/j/jmansour/truthbrowser
# first time, start a new RootCore project:
# rcSetup Base,2.3.10
# rc make_skeleton TruthBrowser
source rcSetup.sh
# rc find_packages
# rc compile
TruthBrowser ~/topMCGen/gen3/DAOD_TRUTH1.test.pool.root

Make Rivet plots:

...

-- JasonDhiaMansour - 2015-07-06

Edit | Attach | Watch | Print version | History: r5 < r4 < r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r5 - 2015-08-10 - JasonDhiaMansour
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

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