![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
SLHCの物理。数が少いので、同種2レプトントリガーの使えるHー＞WWをふたつ捕まえるチャンネルを調べる。 とりあえずHー＞lnuqqのイベントを20kとって来てAnalysisSkeltonをうごかした。
http://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/EXOTICS/
WZの共鳴生成のpythia scriptがあるので、WorkBookの12.0.2のとおりの設定で、Generationをやってみる。(atx7)でもatlfastがうまく動かないので11.0.3に戻ってやる。粒子やジェットを作っておかないと質量再構成ができないので、generationのチェックができないから。(atx13)
11.0.3だと中身のあるAODファイルができた。でもZeeZmmOnAODExample_jobOptions.pyで解析してもWARNING No AOD electron container found in TDSとなって、Zが作れない。 最初はUserAnalysisでうまくいかず。AnalysisExamplesでやりなおしても変なエラーが出るので以下をした。
pool_insertFileToCatalog WZResonance.pool.root pool_insertFileToCatalog AOD_WZ_FastSim.pool.root
でもelectron containerがないのは同じ。rootでAODをみるとelectron containerがあるのによくわからない。
less CLIDDBout.txt /Electron 1341992142 ElectronContainer
Atlfastで作られたElectronContainerは当然のごとくKeyが違うので、JOのコンテナ名のKeyにAtlfastをつける。これ常識。
http://www.hep.ucl.ac.uk/atlas/atlfast/
のUsersGuide参照。これでうまく動くと思ったが、やっぱり駄目だった。結局atx5の10.0.4を試すことにする。atx5はおそい。atx13の11.0.3でAnalysisExamplesは失敗したが、UserAnalysisではやっぱりAtlfastを付けると成功した。??
ChiLのpythiaの変更はどうもathenaに組み込まれているようだ。12.0.2でmissalingnやパイルアップをすればユニークな発表が作れるかも知れないが時間がなくなってきた。
RecExCommonのRecExCommon_topOptions.pyにAODクラスの記述がある。PhysicsAnalysis/AnalysisCommon
include( "ParticleEventAthenaPool/AOD_PoolCnv_jobOptions.py") include( "ParticleEventAthenaPool/AOD_OutputList_jobOptions.py") # set output item list StreamAOD.ItemList += AOD_ItemList
AnalysisExamples ではZll、H4l、ttbarがある。ttbarをgenするにはget_files toprex.datが必要。Simulationはeta cut+-3, DeltaIntersection(Old value:x100)で400秒ぐらいかかる。
https://uimon.cern.ch/twiki/bin/view/Atlas/WorkBook
11.0.1でAOD作成をしようとしたが以下のエラーが出る。Release Statusに対処方法あり。チェックアウト中。対処法にしたがってやっても結果は同じ。11.0.2をインストールした。
history1202WBook.txtTH1.Print Name = TridimMEN2T, Entries= 2800, Total sum= 1
Title = ln(N) vs (E fraction)**0.7 vs Mass/(Mass+1)
NbinsX= 20, xmin= 0, xmax=1, NbinsY= 20, ymin= 0, ymax=1, NbinsZ= 7, zmin= 0, zmax=3.8
THistSvc.sysFinalize() FATAL Standard std::exception is caught
THistSvc.sysFinalize() ERROR St9bad_alloc
ServiceManager WARNING Finalization of service THistSvc failed
10.0.4で当面は演習する。GEN_EVENTのコメントアウトを忘れないこと。
http://wwwlhc.icepp.s.u-tokyo.ac.jp/~jtanaka/ATLAS/ana/tutorials/atlasj/tutorial/ana01.html
/atlas/atlasj/atlas_japan/pro/src/atlasj_atlfast/atlasj_atlfast/ATLASJ_Event.h
CBNT 解析環境のSLC３計算機への移植は、library をあわせたあと、src のconfigure.inをみて include が正しくセットされるようにする。詳しくは以下を参照。
http://kkuatx.hep.kyokyo-u.ac.jp/~ryuichi/pukiwiki/index.php?myMemo
Rome Workshop 用粒子生成プロダクションスクリプトを使い、更にatlfastを使ってCBNTを作った。ヒッグス質量分布を作ると以下のようになる。山下くんの作ったPPTファイル。
http://www.phys.ufl.edu/~jblender/herwig/herwig.html
これでmSUGRAを選択し、ROME SU3パラメーターをいれると読みやすいisajet fileと数字だけのherwig susy fileができる。ただし、isajet7.67となっている。
Event Generator をどうやって走らせているのか？
cmt co Generators/GeneratorOptionsRome
で調べられる。AtlasG4Aimのワークエリアでcheck outしてrome.004202.JimmyZmumu.py をそのままrunしてみた。あっというまに終了。
https://uimon.cern.ch/twiki/bin/view/Atlas/RomeListOfSamples
rome.004421.evgen.SU1_Jimmy_coann._000xx.pool.root with xx=01 to 40 (5K events per file) (Sven Vahsen) (5=) SU1 DC2 stau coannihilation model (isasugra 7.69? m0=70 geV m12 = 350 geV A0 = 0 tan beta = 10 sgn(mu) = + ) 100k 10K events corresponds to an integrated luminosity of 1480 inverse picobarns. Hinchliffe
SU1 Jimmy.JimmyCommand += ["susyfile mean.isawig"] SU3 Jimmy.JimmyCommand += ["susyfile dc1new.txt"]
以下のチュートリアルにはsusyfileの作り方はなかった。
http://stavrop.home.cern.ch/stavrop/Generators/tutorial/generators.html
こちらにはisawigというsusyfileのgeneratorヘのリンクが有る。
http://www-theory.lbl.gov/~ianh/monte/Generators/Herwig/
リンクは
http://www-thphys.physics.ox.ac.uk/users/PeterRichardson/HERWIG/isawig.html
http://atlassw1.phy.bnl.gov/lxr/source/atlas/PhysicsAnalysis/SUSYPhys/SUSYPhysUser/
