http://epp.phys.kyushu-u.ac.jp/~oda/pukiwiki/index.php?EventDisplay2016 既に作られたDAOD_RPVLLファイルを使ってVP1を実行する場合 †ファイルをダウンロードする †setupATLAS lsetup emi voms-proxy-init -voms atlas lsetup rucio rucio get --nrandom 1 data16_13TeV.00311481.physics_Main.recon.DAOD_RPVLL.f758_m1616_r8669 Athenaの設定 †asetup AtlasProduction 20.7.8.7 # r8669 uses AtlasProduction 20.7.8.7 VP1の実行 †vp1 data16_13TeV.00311481.physics_Main.recon.DAOD_RPVLL.f758_m1616_r8669/DAOD_RPVLL.* "All studies" か "Analysis Object (AOD) only studies" を選択。 次のイベントに移動しようとすると、クラッシュしてしまうことが多い。 特定のイベント選びたい時は-evtfilterオプションが使える。例えば、 vp1 -evtfilter=311481:742611918 data16_13TeV.00311481.physics_Main.recon.DAOD_RPVLL.f758_m1616_r8669/DAOD_RPVLL.09788089._000727.pool.root.1 selects run 311481 event 742611918. -evtfilter=FILTER : Only process the events in the FILTER list. Format is a comma separated list of event numbers, where each event number can optionally be preceeded by a run number and a semicolon. Prepend the list with an "%" to instead exclude the events. For instance select a few nice events like: --evtfilter=7,13,5300:26,54 Or exclude a few bad events like: --evtfilter=%7,13,5300:26,54 It is also possible to specify one or more files with run/evt numbers specified on each line: --evtfilter=%7,13,5300:26,54,myevents1.txt Atlantisを実行する場合 †見たいイベントのRAWファイルを取得する †setupATLAS lsetup 'asetup AtlasProduction 20.7.8.7 here' lsetup emi voms-proxy-init -voms atlas lsetup panda 見たいイベントのラン番号とイベント番号を含むテキストファイルを作る。 echo 311481 742611918 > Run311481_Event742611918.txt それを用いてグリッドジョブを投げる。 prun --exec "acmd.py filter-files -s Run311481_Event742611918.txt -o RAW.Run311481_Event742611918.data %IN" --outputs RAW.Run311481_Event742611918.data --athenaTag AtlasProduction,20.7.8.7 --eventPickDataType RAW --eventPickStreamName physics_Main --eventPickEvtList Run311481_Event742611918.txt --noBuild --outDS user.oda.RAW.Run311481_Event742611918/ RAWファイルはディスク上にないことが多いので、少し時間が掛かる。 lsetup rucio rucio get user.oda.RAW.Run311481_Event742611918_RAW.Run311481_Event742611918.data/ 複数のイベントをまとめて抜き出したい時 †以下のようにラン番号とイベント番号をスペースで区切ったRun_Event.listという名前のファイルを作る。 303079 412182676 310341 433397286 310370 1504194326 310370 433397286 310405 1515131235 310863 915987167 310969 1420586256 311170 1605219278 311244 1279674068 311287 151741500 311473 1156910711 311481 1179856748 311481 1423908103 311481 1515115319 311481 1560163113 311481 1638352153 311481 1805419211 311481 1903118424 311481 2070604365 311481 574987370 311481 742611918 311481 968441316 311481 972014010 311481 999969045 さらに、以下の内容のSubmit.pyというファイルを作る import os for line in open('Run_Event.list', 'r'): itemList = line[:-1].split(' ') run=itemList[0] event=itemList[1] os.system("echo \""+run+" "+event+"\">Run_Event.txt") os.system("prun --exec \"acmd.py filter-files -s Run_Event.txt -o RAW.Run"+run+"_Event"+event+".data %IN\" --outputs RAW.Run"+run+"_Event"+event+".data --athenaTag AtlasProduction,20.7.8.7 --eventPickDataType RAW --eventPickStreamName physics_Main --eventPickEvtList Run_Event.txt --noBuild --outDS user.oda.RAW.HNL.Run."+run+"_Event"+event+"/") (出力データセット名を変えないと連続してジョブを投げられない。) setupATLAS lsetup 'asetup AtlasProduction 20.7.8.7 here' lsetup emi voms-proxy-init -voms atlas lsetup panda を行ってから、 python Submit.py を実行する。 Jive XMLファイルの生成 †r8669のコマンドのpreExecに'RAWtoESD:jp.Rec.doJiveXML.set_Value_and_Lock(True)'を追加する。 setupATLAS asetup AtlasProduction 20.7.8.7 Reco_tf.py --AMI r8669 --inputBSFile user.oda.RAW.Run311481_Event742611918_RAW.Run311481_Event742611918.data/user.oda.11593022._000001.RAW.Run311481_Event742611918.data --outputDAOD_RPVLLFile DAOD_RPVLL.pool.root --maxEvents 10 --preExec 'RAWtoESD:from InDetRecExample.InDetJobProperties import InDetFlags;InDetFlags.doLargeD0.set_Value_and_Lock(True);InDetFlags.keepAdditionalHitsOnTrackParticle.set_Value_and_Lock(True);' 'ESDtoDPD:from MuonRecExample.MuonRecFlags import muonRecFlags;muonRecFlags.prdToxAOD=True;' 'RAWtoESD:jp.Rec.doJiveXML.set_Value_and_Lock(True)' を実行すると、JiveXML_311481_742611918.xmlというファイルができる。 Atlantisの実行 †setupATLAS lsetup atlantis atlantis |