EventDisplay_RPVLL

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" を選択。
"All studies" を選択した場合は、その次の画面で、Analysisボタンを押す。
Muonsなどを選択する。

次のイベントに移動しようとすると、クラッシュしてしまうことが多い。

特定のイベント選びたい時は-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