140 likes | 236 Views
Generation of W+jets with cut on W pt. Alexander Solodkov IHEP, Protvino. Software versions. Following recommendation of Philippe Calfayan I tried to run modified version of Alpgen with cut on boson pt
E N D
Generation of W+jets with cut on W pt Alexander Solodkov IHEP, Protvino Sanya.Solodkov@cern.ch
Software versions • Following recommendation of Philippe Calfayan I tried to run modified version of Alpgen with cut on boson pt • Last Athena release – 17.7.0 contains latest versions of all Alpgen packages, it means that production can be easily done on the grid • External/AlpGen-00-00-06 uses Alpgen v2.1.4.p5 • Generators/AlpGen_i-00-00-43 • Generators/AlpGenControl-00-00-22 Sanya.Solodkov@cern.ch
Test with Z tau tau Np1 • For first test with cut on Z pt the following datasets were produced (5000 events each) • 147122.ZtautauNp1 - no pt cut • 147154.ZtautauNp1 - 0 < pt < 40 • 147186.ZtautauNp1 - 40 < pt < 70 • 147218.ZtautauNp1 - 70 < pt < 140 • 147250.ZtautauNp1 - 140 < pt < 280 • 147282.ZtautauNp1 - 280 < pt < 500 • 147314.ZtautauNp1 - 500 < pt Sanya.Solodkov@cern.ch
Test with Z tau tau Np1 • For pt > 500 GeV job failed with error: • TestHepMC ERROR EFFICIENCY IS TOO LOW! 97.9624% found, but at least: 98% required AthenaEventLoopMgr ERROR Error in Algorithm Finalize severity=ERROR severity=ERROR • Had to rerun job with extra parameter • postExec=topAlg.TestHepMC.EffFailThreshold=0.975 • Next few slides – attempt to reproduce plot shown in this presentation: • https://indico.cern.ch/getFile.py/access?contribId=2&resId=0&materialId=slides&confId=254873 Sanya.Solodkov@cern.ch
Test with Z tau tau Np1 (no sum) Sanya.Solodkov@cern.ch
Test with Z tau tau zoomed Sanya.Solodkov@cern.ch
Test with Z tau tau stacked hist Sanya.Solodkov@cern.ch
Test with Z tau tau stacked hist Sanya.Solodkov@cern.ch
Test with W+jets • W+njets without pt cut and with pt>150 cut were produced for 3 possible decay modes and number of jets 0 <= n <=5 – 36 datasets in total • 6 jobs failed • WenuNp5incl • WtaunuNp5incl • WenuNp0_pt150 • WenuNp5incl_pt150 • WmunuNp0_pt150 • WtaunuNp0_pt150 Sanya.Solodkov@cern.ch
Test with W+jets • In failed jobs errors are always similar • At line 108 of file ../src/initalpgen.F (unit = 71, file = 'inparmAlpGen.dat') • Fortran runtime error: End of file • At line 2064 of file ../alplib/alpgen.f (unit = 11, file = 'alpout.par') • Fortran runtime error: End of file • At line 3614 of file ../alplib/alpgen.f (unit = 11, file = 'alpout.par') • Fortran runtime error: End of file • Log files are available at http://cern.ch/solodkov/alpgen • Asked Philippe if he knows what’s wrong Sanya.Solodkov@cern.ch
W e nu + 1,2,3,4 jetss Sanya.Solodkov@cern.ch
W mu nu + 1,2,3,4 jetss Sanya.Solodkov@cern.ch
W tau nu + 1,2,3,4 jetss Sanya.Solodkov@cern.ch
Conclusions • Event generation with cut on boson pt seems to be working • WlnuNp1 – WlnuNp4 are fine • Still need to understand the reason of failures in some jobs for Np0 and Np5inclusive • I guess next step – to produce more “official” datasets • How to proceed ? • Private production or request to production team ? • Which release should be used ? • Output from 17.7.0 might be incompatible with 17.2.X.Y • How many events are needed ? Sanya.Solodkov@cern.ch