Permissions:

find . -type d -exec fs setacl -dir {} -acl system:anyuser rl \;

Grid:

General MC production

  1. https://svnweb.cern.ch/trac/atlasphys/browser/Physics/Generators/ExampleJobOptions [example job options]
  2. https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/RunningGenerators
  3. https://twiki.cern.ch/twiki/bin/view/AtlasProtected/HowToTestProdJobs [validation: test JOs]
  4. https://indico.cern.ch/event/273749/session/12/contribution/43/material/slides/0.pdf [Ulla: MC production for Run II - Higgs] {SLIDES}
  5. https://indico.cern.ch/event/338847/contribution/5/material/slides/0.pdf [Rei: Higgs (N)NLO MC&Tools for RUN-2] {SLIDES}
  6. http://home.thep.lu.se/~torbjorn/pythia81html/ParticleProperties.html [Status code and Particle properties]
  7. https://cds.cern.ch/record/1753200/files/ATL-COM-SOFT-2014-036.pdf [Monte Carlo Truth Task Force Report]
  8. PDF lhaglue numbering: https://lhapdf.hepforge.org/pdfsets.html
  9. Validation rules for MCgenerators ATLAS: https://twiki.cern.ch/twiki/bin/view/AtlasProtected/GeneratorValidationRules
  10. Matrix Element validation tool - valME: https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/MatrixElementValidation
  11. General validation rules in ATLAS: https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/GeneratorValidationRules
  12. LoopTools: http://www.feynarts.de/looptools/
  13. Dump HepMC format from EVNT file:

setupATLAS

asetup 20.0.0.3,AtlasProduction,here

cmt co Event/xAOD/xAODTruthCnv

cd Event/xAOD/xAODTruthCnv/cmt

cmt config

gmake

cd ../

mkdir run

cd run

cp ../share/* .

Open the file called

HepMCTruthReader_jobOptions.py

and replace “evnt.pool.root” with the name of your EVNT file.

athena HepMCTruthReader_jobOptions.py > output.txt

Rivet

Sherpa

mg5_aMC: madgraph+MC@NLO [lxplus]

The output of MadGraph is LHE File, while the input for a Rivet Analysis is a hepmc file. So the first step is to install the following converter: lhef2hepmc.tar.gz After unpacking, you can install it via the command

make HEPMC_PREFIX=/home/student/local/common

Then you can run the code

./lhef2hepmc INPUT.LHE OUTPUT.HEPMC

Then you can run RIVET on that file

rivet -a MC_MCNET2014 OUTPUT.HEPMC rivet-mkhtml Rivet.yoda

Powheg / Powheg OTF

Setup and run Powheg OTF on lxplus (tested with athena version 19.1.1.5 -see setup below-)

POWHEG-BOX-V2 release: 2819

mkdir testarea/
cd testarea/

Download packages:

  •   MC14JobOptions-00-00-33
  •   POWHEG-BOX-V2
  •   PowhegControl-00-00-12
  •   External/Powheg-00-00-10
 svn checkout [--revision n] --username anonymous --password anonymous svn://powhegbox.mib.infn.it/trunk/POWHEG-BOX-V2
 svn list [--revision n] --username anonymous --password anonymous svn://powhegbox.mib.infn.it/trunk/User-Processes-V2
 svn co [--revision n] --username anonymous --password anonymous svn://powhegbox.mib.infn.it/trunk/User-Processes-V2/XXX
 

In additional: download POWHEG-BOX-V2 processes and compile them inside POWHEG-BOX

NOTE: to compile HVJ change the FC compiler to gfortran inside the Makefile.

(1) MC14JobOptions

  • modify JO in share/
  • modify gencontrol/PowhegControl_postInclude.py to add relevant processes
  • modify/add gencontrol/PowhegControl_<process>.py to add/modify relevant processes
(2) POWHEG-BOX-V2: download and compile
(3) PowhegControl -00-00-12: compile

(4) External/Powheg: compile

###################

setupATLAS 
asetup AtlasProduction,19.1.1.5,slc6,gcc47,64,here # <=== ATHENA VERSION 19.1.1.5

export POWHEGPATH=/afs/<a target="_blank" href="http://cern.ch/work/c/cpandini/testarea/mcgen_powhegOTF/POWHEG-BOX-V2">cern.ch/work/c/cpandini/testarea/mcgen_powhegOTF/POWHEG-BOX-V2</a> 
export LHAPATH=/afs/<a target="_blank" href="http://cern.ch/sw/lcg/external/MCGenerators_hepmc2.06.08/lhapdf/5.8.9/share/PDFsets">cern.ch/sw/lcg/external/MCGenerators_hepmc2.06.08/lhapdf/5.8.9/share/PDFsets</a>

export FASTJET_CONFIG=/afs/<a target="_blank" href="http://cern.ch/sw/lcg/external/fastjet/3.0.3/x86_64-slc6-gcc46-opt">cern.ch/sw/lcg/external/fastjet/3.0.3/x86_64-slc6-gcc46-opt</a>
export PATH=$FASTJET_CONFIG/bin:$PATH
export FASTJETPATH=$FASTJET_CONFIG

# pdfs                                                                                                                                                                                                       
export PDFS=/afs/<a target="_blank" href="http://cern.ch/sw/lcg/external/MCGenerators_hepmc2.06.08/lhapdf/5.8.9/x86_64-slc6-gcc46-opt">cern.ch/sw/lcg/external/MCGenerators_hepmc2.06.08/lhapdf/5.8.9/x86_64-slc6-gcc46-opt</a>
export PDFLIB=$PDFS/lib
export LD_LIBRARY_PATH=$PDFLIB:$LD_LIBRARY_PATH
export PATH=$PATH:$PDFS/bin
###################

Generate_tf.py --ecmEnergy=13000. --runNumber=189900 --firstEvent=1 --maxEvents=5000 --randomSeed=54298752 --jobConfig=<a target="_blank" href="http://MC14.189900.PowhegPythia8_AU2CT10_ggZH125_eebb.py">MC14.189900.PowhegPythia8_AU2CT10_ggZH125_eebb.py</a> --outputEVNTFile=testcarlo.189900_powheg_highstat5k.EVNT.pool.root
###################

Powheg+Pythia matching

Powheg input card parameters: powheg.input

Parameters to test ggZH with POWHEG-BOX-V2 OTF (see above):

  • ncall1 400000
  • itmx1 4 [Olivier(1) - commented ~70 seconds]
  • ncall2 5000000 [Olivier(40000) - much lower - test]
  • itmx2 5 [Olivier(1) - test]
  • foldx 1
  • foldcsi 1 [alternate name for foldx]
  • foldy 1
  • foldphi 1
  • nubound 1000000 [Olivier 2M - test]
  • xupbound 2 [only in Olivier input card]
  • ixmax 1
  • icsimax 1
  • iymax 1
  • charmthr 1.5
  • bottomthr 5.0
  • charmthrpdf 1.5
  • lhans1/lhans2 11000 [set CT10nlo - Olivier/old using 10800 CT10]
  • withnegweights 1 [only Olivier - allow negative weights]
  • facscfact/renscfact 1.0 [as default]
  • pdfreweight 1 [allow pdf reweighting]
  • hwidth 0.00403 [in GeV -> modify Olivier's inputcard]
  • vdecaymode/hdecaymode [set in JO] --> hardcoded to generate ee,mm,tautau,veve,vmvm,vtauvtau separately
  • storeinfo_rwgt 1 [store info for reweighting]
  • doublefst 1 [?]
  • nohad 1 [?]
  • LOevents 1 [ok for ggZH, change for qqVH]
  • bornktmin [default 0d0 - generation cut - minimum kt in underlying Born]
  • bornsuppfact [default 0d0 - mass parameter for Born suppression factor - if < 0 suppfact = 1]
  • manyseeds + other parameters for parallel run -> not set for now [ to be tested ]

**Diboson**

Meetings

Useful old meetings [indico links]

C++ useful references:

Unfolding:

EFT effective field theories:

My Links

My Personal Preferences

  • Show tool-tip topic info on mouse-over of WikiWord links, on or off:
    • Set LINKTOOLTIPINFO = off

  • Preference for the editor, default is the WYSIWYG editor. The options are raw, wysiwyg:
    • Set EDITMETHOD = wysiwyg

  • More preferences TWiki has system wide preferences settings defined in TWikiPreferences. You can customize preferences settings to your needs: To overload a system setting, (1) do a "raw view" on TWikiPreferences, (2) copy a Set VARIABLE = value bullet, (3) do a "raw edit" of your user profile page, (4) add the bullet to the bullet list above, and (5) customize the value as needed. Make sure the settings render as real bullets (in "raw edit", a bullet requires 3 or 6 spaces before the asterisk).

Related Topics

UserForm
First Name Carlo Enrico
Last Name Pandini
Titles

Email carlo.enrico.pandini@cernNOSPAMPLEASE.ch
Telephone

Mobile

Skype ID

Department

Organization

URL

Location

Region

Country Italy
Status Update

Edit personal data
Edit | Attach | Watch | Print version | History: r28 < r27 < r26 < r25 < r24 | Backlinks | Raw View | Raw edit | More topic actions
Topic revision: r28 - 2015-01-29 - CarloEnricoPandini
 
    • 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