ATLASソフトウエア講習会2016


ATLASでは主に3つのイベントディスプレイ

がありますが、ここではAtlantisとVP1の使い方を説明します。


Atlantisの使い方

.thumb_fig_14.png

VP1の使い方

.thumb_fig_33.png

ATLAS実験の公式のイベントディスプレイは
https://twiki.cern.ch/twiki/bin/view/AtlasPublic/EventDisplayPublicResults
にあります。

https://twiki.cern.ch/twiki/bin/viewauth/Atlas/OfficialEventDisplayGuidelines


特定のイベントのファイルを取得する方法

例として、
https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/CONFNOTES/ATLAS-CONF-2016-045/
Auxiliary Figure Figure 5で使われているRun 302393 Event 3804660240のイベントに注目します。

AODファイルを取って来る方法

http://ami.in2p3.fr/
のDataset BrowserのSimple Searchで

data16_13TeV.%302393%physics_Main.merge.AOD.%

を検索して、公式に作られたAODがどのタグを使っているか調べます。
このランがPeriod C4であることが表示され、
https://twiki.cern.ch/twiki/bin/view/AtlasProtected/DataMCForAnalysis#Data_and_MC_processing_campaigns
から、
https://twiki.cern.ch/twiki/bin/view/AtlasProtected/DataMCForAnalysis2015#Preferred_data_sets_data16_13_Te
を見ると、Physics Container data16_13TeV.periodC.physics_Main.PhysCont.AOD.t0pro20_v01を使えば良いことがわかります。
Dataset BrowserのSimple Searchで

data16_13TeV.periodC.physics_Main.PhysCont.AOD.t0pro20_v01

を検索し、Contained_datasetsをクリックすると、データセット
data16_13TeV.00302393.physics_Main.merge.AOD.f711_m1620 を使えば良いことがわかります。
http://ami.in2p3.fr/
のAMI-TagsのShow/Editで
でf711を調べると、AtlasProduction 20.7.6.7を使っていることがわかります。

Pandaの設定をします。

$ setupATLAS
$ lsetup 'asetup AtlasProduction 20.7.6.7 here'
$ lsetup emi
$ voms-proxy-init -voms atlas 
$ lsetup panda

Run番号とEvent番号を含むテキストファイルを作ります。

$ echo 302393 3804660240 > Run302393_Event3804660240.txt

prunで該当のイベントを含むAODファイルを抜き出し、acmd.pyコマンドを走らせ、該当のイベントのみ抜き出し、保存します。

$ prun --exec "acmd.py filter-files -s Run302393_Event3804660240.txt -o AOD.Run302393_Event3804660240.pool.root %IN" --outputs AOD.Run302393_Event3804660240.pool.root --athenaTag AtlasProduction,20.7.6.7  --eventPickDataType AOD --eventPickStreamName physics_Main --eventPickAmiTag f711 --eventPickEvtList Run302393_Event3804660240.txt --noBuild --outDS user.oda.AOD.Run302393_Event3804660240/

userで始まる出力のデータセット名は適宜変えて下さい。
ジョブが終われば、メールが来るので、rucio getで出力データセットを取って来ます。

$ setupATLAS
$ lsetup emi
$ voms-proxy-init -voms atlas
$ lsetup rucio
$ rucio get user.oda:user.oda.AOD.Run302393_Event3804660240_AOD.Run302393_Event3804660240.pool.root/

複数のイベントを1ジョブで取って来ることも可能なはずですが、1ジョブ1イベントずつにした方が無難です。

RAWファイルを取って来る方法

AODの場合と同様に、

$ prun --exec "acmd.py filter-files -s Run302393_Event3804660240.txt -o RAW.Run302393_Event3804660240.data %IN" --outputs RAW.Run302393_Event3804660240.data --athenaTag AtlasProduction,20.7.6.7  --eventPickDataType RAW --eventPickStreamName physics_Main --eventPickEvtList Run302393_Event3804660240.txt --noBuild --outDS user.oda.RAW.Run302393_Event3804660240/

を実行します。AMIタグは指定する必要ありません。

$ rucio get user.oda.RAW.Run302393_Event3804660240_RAW.Run302393_Event3804660240.data/

RAWファイルをESDとAODファイルに変換する方法

ESDファイルは保存されていないので、検出器の詳しい情報を見たい場合、RAWファイルに戻る必要があります。

実データを処理する時は、データベースのアクセスの関係上login.icepp.jpだと非常に時間が掛るので、lxplus.cern.chで行った方が良いです。

https://ami.in2p3.fr/app?subapp=amiTags_show&userdata=f711
でわかるf711タグの設定にならい、

$ setupATLAS
$ lsetup 'asetup AtlasProduction 20.7.6.7'

とAtlasProduction 20.7.6.7を設定します。

login.icepp.jpを使う時はFroNTierのプロキシを追加します。

$ export FRONTIER_SERVER=$FRONTIER_SERVER"(proxyurl=http://conddb-px02.icepp.jp:3128)"

RAWファイルのシンボリックリンクを作成します。

$ ln -s user.oda.RAW.Run302393_Event3804660240_RAW.Run302393_Event3804660240.data/user.oda.10284977._000001.RAW.Run302393_Event3804660240.data RAW.Run302393_Event3804660240.data

そして、

$ Reco_tf.py --AMI f711 --inputBSFile RAW.Run302393_Event3804660240.data --outputESDFile ESD.Run302393_Event3804660240.pool.root 

を実行します。

--outputAODFile AOD.Run302393_Event3804660240.pool.root

を付け加えれば、AODファイルも同時に作れます。

Jive XMLファイルを作る

$ GetTfCommand.py --AMI f711

を実行すると、

PyJobTransforms.<module> 2016-12-21 14:22:55,672 INFO logging set in /cvmfs /atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc49-opt/20.7.6/AtlasProduction/20.7.6.7/InstallArea/share/bin/GetTfCommand.py

Information about tag f711:
This is a T0 tag.
This tag consists of 1 transform command(s).
Transform commands follow below.
Input and output file names (if present) are only suggestions.

asetup AtlasProduction,20.7.6.7
Reco_tf.py --conditionsTag all:CONDBR2-BLKPA-2016-14 --ignoreErrors 'False' --autoConfiguration='everything' --maxEvents '-1' --AMITag 'f711' --postExec  'e2d:from AthenaCommon.AppMgr import ServiceMgr;import MuonRPC_Cabling.MuonRPC_CablingConfig;ServiceMgr.MuonRPC_CablingSvc.RPCMapfromCool=False;ServiceMgr.MuonRPC_CablingSvc.CorrFileName="LVL1confAtlasRUN2_ver016.corr";ServiceMgr.MuonRPC_CablingSvc.ConfFileName="LVL1confAtlasRUN2_ver016.data";' --preExec  'all:jobproperties.Beam.bunchSpacing.set_Value_and_Lock(25);'  'r2e:from InDetPrepRawDataToxAOD.SCTxAODJobProperties import SCTxAODFlags;SCTxAODFlags.Prescale.set_Value_and_Lock(50);' --geometryVersion all:ATLAS-R2-2015-04-00-00

Input file arguments:
 --inputBSFile '/afs/cern.ch/atlas/offline/test/data11_7TeV.00191920.physics_JetTauEtmiss.merge.RAW._lb0257._SFO-9._0001.1.10evts'

Output file arguments:
 --outputDAOD_IDNCBFile 'myDAOD_IDNCB.pool.root' --outputDRAW_EGZFile 'myDRAW_EGZ.data' --outputDESDM_IDALIGNFile 'myDESDM_IDALIGN.pool.root' --outputDESDM_SGLELFile 'myDESDM_SGLEL.pool.root' --outputDESDM_SLTTMUFile 'myDESDM_SLTTMU.pool.root' --outputNTUP_MUONCALIBFile 'myNTUP_MUONCALIB.root' --outputDRAW_ZMUMUFile 'myDRAW_ZMUMU.data' --outputAODFile 'myAOD.pool.root' --outputDAOD_IDTIDEFile 'myDAOD_IDTIDE.pool.root' --outputDRAW_TAUMUHFile 'myDRAW_TAUMUH.data' --outputDESDM_EGAMMAFile 'myDESDM_EGAMMA.pool.root' --outputDESDM_MCPFile 'myDESDM_MCP.pool.root' --outputDESDM_CALJETFile 'myDESDM_CALJET.pool.root' --outputDESDM_PHOJETFile 'myDESDM_PHOJET.pool.root' --outputDESDM_TILEMUFile 'myDESDM_TILEMU.pool.root' --outputDRAW_RPVLLFile 'myDRAW_RPVLL.data' --outputDAOD_IDTRKVALIDFile 'myDAOD_IDTRKVALID.pool.root' --outputDAOD_SCTVALIDFile 'myDAOD_SCTVALID.pool.root' --outputDESDM_EXOTHIPFile 'myDESDM_EXOTHIP.pool.root' --outputESDFile 'myESD.pool.root' --outputHISTFile 'myHIST.root'

AMI outputs:
{'outputDAOD_IDNCBFile': {'dstype': 'DAOD_IDNCB', 'ifMatch': 'data[0-9][0-9]_(cos|1beam|.*eV|comm)(?!.*DRAW.*)(.*Background.*)'}, 'outputDESDM_CALJETFile': {'dstype': 'DESDM_CALJET', 'ifMatch': 'data[0-9][0-9]_.*eV(?!.*DRAW.*)(.*physics_Main.*)'}, 'outputDESDM_EGAMMAFile': {'dstype': 'DESDM_EGAMMA', 'ifMatch': 'data[0-9][0-9]_.*eV(?!.*DRAW.*)(.*physics_Main.*)'}, 'outputDESDM_SLTTMUFile': {'dstype': 'DESDM_SLTTMU', 'ifMatch': 'data[0-9][0-9]_.*eV(?!.*DRAW.*)(.*physics_Main.*)'}, 'outputDRAW_EGZFile': {'dstype': 'DRAW_EGZ', 'ifMatch': 'data[0-9][0-9]_.*eV(?!.*DRAW.*)(.*physics_Main\..*)'}, 'outputDESDM_IDALIGNFile': {'dstype': '!replace RAW ESDM', 'ifMatch': 'data[0-9][0-9]_(1beam|.*eV|comm)(.*physics_Main.*)(.*\.DRAW_EGZ\.*|.*DRAW_ZMUMU.*)'}, 'outputDESDM_SGLELFile': {'dstype': 'DESDM_SGLEL', 'ifMatch': 'data[0-9][0-9]_.*eV(?!.*DRAW.*)(.*physics_Main.*)'}, 'outputAODFile': {'dstype': '!replace RAW AOD', 'ifMatch': 'data[0-9][0-9]_(cos|1beam|.*eV|comm|hi)(?!.*DRAW_RPVLL.*)(?!.*Background.*)'}, 'outputNTUP_MUONCALIBFile': {'dstype': 'NTUP_MUONCALIB', 'ifMatch': 'data[0-9][0-9]_(1beam|.*eV|comm)(?!.*DRAW.*)(?!.*Background.*)'}, 'outputDRAW_ZMUMUFile': {'dstype': 'DRAW_ZMUMU', 'ifMatch': 'data[0-9][0-9]_.*eV(?!.*DRAW.*)(.*physics_Main\..*)'}, 'outputDAOD_IDTIDEFile': {'dstype': 'DAOD_IDTIDE', 'ifMatch': 'data[0-9][0-9]_.*eV(?!.*DRAW.*)(.*physics_Main.*)'}, 'outputDRAW_TAUMUHFile': {'dstype': 'DRAW_TAUMUH', 'ifMatch': 'data[0-9][0-9]_.*eV(?!.*DRAW.*)(.*physics_Main\..*)'}, 'outputDAOD_SCTVALIDFile': {'dstype': 'DAOD_SCTVALID', 'ifMatch': 'data[0-9][0-9]_(cos|1beam|.*eV|comm|hi)(?!.*DRAW.*)(.*express_express.*)'}, 'outputDESDM_MCPFile': {'dstype': 'DESDM_MCP', 'ifMatch': 'data[0-9][0-9]_(1beam|.*eV|comm)(?!.*DRAW.*)(.*physics_Main\..*|.*CosmicMuons.*)'}, 'outputDRAW_RPVLLFile': {'dstype': 'DRAW_RPVLL', 'ifMatch': 'data[0-9][0-9]_.*eV(?!.*DRAW.*)(.*physics_Main\..*)'}, 'outputDESDM_PHOJETFile': {'dstype': 'DESDM_PHOJET', 'ifMatch': 'data[0-9][0-9]_.*eV(?!.*DRAW.*)(.*physics_Main.*)'}, 'outputDESDM_TILEMUFile': {'dstype': 'DESDM_TILEMU', 'ifMatch': 'data[0-9][0-9]_.*eV(?!.*DRAW.*)(.*physics_Main.*)'}, 'outputDAOD_IDTRKVALIDFile': {'dstype': 'DAOD_IDTRKVALID', 'ifMatch': 'data[0-9][0-9]_(cos|1beam|.*eV|comm)(?!.*DRAW.*)(.*MinBias.*|.*IDCosmic.*)'}, 'outputDESDM_EXOTHIPFile': {'dstype': 'DESDM_EXOTHIP', 'ifMatch': 'data[0-9][0-9]_(1beam|.*eV|comm)(?!.*DRAW.*)(.*physics_Main\..*)'}, 'outputESDFile': {'dstype': '!replace RAW ESD', 'ifMatch': '(?!.*DRAW.*)(?!.(.*physics_Main\..*|.*Background.*|.*L1Topo.*))'}, 'outputHISTFile': {'dstype': 'HIST', 'ifMatch': '(?!.(.*DRAW.*|.*debugrec.*))'}}

Possible output data types: ['DAOD_IDNCB', 'DESDM_CALJET', '!replace RAW ESDM', 'DESDM_SGLEL', 'HIST', 'DRAW_ZMUMU', 'DESDM_PHOJET', 'DESDM_SLTTMU', 'DAOD_IDTIDE', 'DRAW_TAUMUH', 'DESDM_EGAMMA', 'DESDM_MCP', 'DRAW_EGZ', 'NTUP_MUONCALIB', 'DESDM_TILEMU', 'DRAW_RPVLL', 'DAOD_IDTRKVALID', 'DAOD_SCTVALID', 'DESDM_EXOTHIP', '!replace RAW ESD', '!replace RAW AOD']

と出て来ます。この--preExecに

'r2e:jp.Rec.doJiveXML.set_Value_and_Lock(True)'

を付け加えることで、ESDとJiveXMLファイルが同時に作れます。

$ Reco_tf.py --conditionsTag all:CONDBR2-BLKPA-2016-14 --ignoreErrors 'False' --autoConfiguration='everything' --maxEvents '-1' --AMITag 'f711' --postExec  'e2d:from AthenaCommon.AppMgr import ServiceMgr;import MuonRPC_Cabling.MuonRPC_CablingConfig;ServiceMgr.MuonRPC_CablingSvc.RPCMapfromCool=False;ServiceMgr.MuonRPC_CablingSvc.CorrFileName="LVL1confAtlasRUN2_ver016.corr";ServiceMgr.MuonRPC_CablingSvc.ConfFileName="LVL1confAtlasRUN2_ver016.data";' --preExec  'all:jobproperties.Beam.bunchSpacing.set_Value_and_Lock(25);'  'r2e:from InDetPrepRawDataToxAOD.SCTxAODJobProperties import SCTxAODFlags;SCTxAODFlags.Prescale.set_Value_and_Lock(50);' 'r2e:jp.Rec.doJiveXML.set_Value_and_Lock(True)' --geometryVersion all:ATLAS-R2-2015-04-00-00 --inputBSFile RAW.Run302393_Event3804660240.data --outputESDFile ESD.Run302393_Event3804660240.pool.root

トップ   編集 凍結解除 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2016-12-22 (木) 14:21:34 (1027d)