updated 2008.5.12
今までやったことを復習してみる：------------------------------------------------------- > source /afs/cern.ch/sw/contrib/CMT/v1r20p20070720/mgr/setup.sh cmthome> cmt config cmthome> source setup.sh -tag=13.0.30 <--これをいつもやること 13.0.30> cmt co -r UserAnalysis-00-10-12 PhysicsAnalysis/AnalysisCommon/UserAnalysis ...> cd cmt cmt> source setup.sh cmt> gmake run> get_files AthExHelloWorld_jobOptions.py run> athena.py AthExHelloWorld_jobOptions.py run> vi myjob run> bsub -q 1nh myjob ---------------------------------------------------------- > cd Control/AthenaExamples/AthExHelloWorld/src > cd ../cmt > cmt config > source setup.sh [or .csh] > gmake > athena.py AthExHelloWorld_jobOptions.py を実行した ------------------------------------------------------------- > cmt create MyNewPackage MyNewPackage-00-00-01 > cd MyNewPackage/cmt > cd ../src --> Create a file MyAlg.cxx: > mkdir components > cd components ---> MyNewPackage_entries.cxx c, MyNewPackage_load.cxxも作れ > mkdir MyNewPackage ---> cd MyNewPackage ----> Create a file MyAlg.h: > cd ..> mkdir share> cd share --->MyJobOptions.py: packageビルド > cd ../cmt ---> (> cmt config > source setup.sh > gmake)のセット！ > cd ~/testarea/13.0.30PhysicsAnalysis/AnalysisCommon/UserAnalysis/run > cp ../../../../MyNewPackage/share/MyJobOptions.py . > athena.py MyJobOptions.py ------------------------------------------------------------（６）More Advanced Setup Options （2008.5.11）
全くわからない。スキップする。（７）Running the Full Chain（2008.5.11）
from Generation to production of Analysis Object Data (AOD) using Atlfast wto produce AOD directly. Atlfast take input from any of the event generator, simulation, digitization, or Event Summary Data (ESD) files.（8）Generating Events with Athena （2008.5.11卒業）
See McGeneratorsForAtlas for further information. Pythiaを例にする。Pythia options はPythia Manualを見るか PythiaForAtlasページを見よ。Running Pythia
run> get_files PDGTABLE.MeV ーー＞内容はこれ DELLのcygwinページから $ scp -r kondo@lxplus.cern.ch:/afs/cern.ch/user/k/kondo/testarea/13.0.30/Physic sAnalysis/AnalysisCommon/UserAnalysis/run/PDGTABLE.MeV .PDGTABLE.MeV とやったらコピーで来た！！粒子のテーブルで質量・幅やその誤差もある。 run> get_files jobOptions.pythia.py############################################################### # # JOB OPTIONS FILE. THIS FILE IS TO DEMONSTRATE HOW THE USER CAN # CHANGE THE PYTHIA PARAMETERS. IT DOES NOT NECESARILY COMPLY TO # THE ATLAS DEFAULTS FOR PHYSICS EVENT GENERATION. #============================================================== #-------------------------------------------------------------- # Private Application Configuration option from AthenaCommon.AppMgr import ServiceMgr ServiceMgr.MessageSvc.OutputLevel = DEBUG #-------------------------------------------------------------- # Event related parameters #-------------------------------------------------------------- # Number of events to be processed (default is 10) theApp.EvtMax = 5 #-------------------------------------------------------------- # Algorithms Private Options #-------------------------------------------------------------- from AthenaServices.AthenaServicesConf import AtRndmGenSvc ServiceMgr += AtRndmGenSvc() ServiceMgr.AtRndmGenSvc.Seeds = ["PYTHIA 4789899 989240512", "PYTHIA_INIT 820021 2347532"] # ServiceManager.AtRndmGenSvc.ReadFromFile = true; from AthenaCommon.AlgSequence import AlgSequence job=AlgSequence() from Pythia_i.Pythia_iConf import Pythia job +=Pythia() job.Pythia.PythiaCommand = ["pysubs msel 13","pysubs ckin 3 18.","pypars mstp 43 2"] ここで MSEL = 13 ーーー＞Z0 + Jet Production CKIN(3) = 18. ーーー＞the minimum PT for hard processes in GeV MSTP(43) = 2 ーーー＞Z0/gamma* interference off (only Z0 included) job.Pythia.PythiaCommand += ["pypars mstp 51 19070", "pypars mstp 52 2", "pypars mstp 53 19070", "pypars mstp 54 2", "pypars mstp 55 19070", "pypars mstp 56 2"] from TruthExamples.TruthExamplesConf import DumpMC job += DumpMC() # # End of job options file # ###############################################################> athena.py jobOptions.pythia.py > athena_gen.out で５イベントがgenerateされた！ Now modify the Pythia options to generate Z0->e+e- events only. myGenOptions.pyを作った。#-------------------------------------------------------------- # Private Application Configuration option #-------------------------------------------------------------- from AthenaCommon.AppMgr import ServiceMgr ServiceMgr.MessageSvc.OutputLevel = INFO #-------------------------------------------------------------- # Event related parameters #-------------------------------------------------------------- # Number of events to be processed (default is 10) theApp.EvtMax = 10 #-------------------------------------------------------------- # Run Number #-------------------------------------------------------------- EventSelector = Service( "EventSelector" ) EventSelector.RunNumber = 999999 #-------------------------------------------------------------- # Algorithms Private Options #-------------------------------------------------------------- from AthenaServices.AthenaServicesConf import AtRndmGenSvc ServiceMgr += AtRndmGenSvc() ServiceMgr.AtRndmGenSvc.Seeds = ["PYTHIA 4789899 989240512", "PYTHIA_INIT 820021 2347532"] これがランダムNo.の種らしい。変えよ。 # from AthenaCommon.AlgSequence import AlgSequence job=AlgSequence() # from Pythia_i.Pythia_iConf import Pythia job +=Pythia() # Generate Z->ee job.Pythia.PythiaCommand = ["pysubs msel 0","pysubs msub 1 1", "pypars mstp 43 2","pydat3 mdme 174 1 0", MSTP(43) = 2 ーーー＞Z0/gamma* interference off (only Z0 included) "pydat3 mdme 175 1 0","pydat3 mdme 176 1 0", "pydat3 mdme 177 1 0","pydat3 mdme 178 1 0", "pydat3 mdme 179 1 0","pydat3 mdme 180 1 0", "pydat3 mdme 181 1 0","pydat3 mdme 182 1 1", ＜−−−e "pydat3 mdme 183 1 0","pydat3 mdme 184 1 0", muonならpydat3 mdme 184 1 1 "pydat3 mdme 185 1 0","pydat3 mdme 186 1 0", "pydat3 mdme 187 1 0"] from TruthExamples.TruthExamplesConf import DumpMC job += DumpMC() #--------------------------------------------------------------- # Pool Persistency #--------------------------------------------------------------- from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream Stream1 = AthenaPoolOutputStream( "Stream1" ) Stream1.OutputFile = "pythia.pool.root" ＜−−ここに出力される。 Stream1.ItemList += [ 'EventInfo#*', 'McEventCollection#*' ] ###############################################################You might also want to change the output level threshold (MessageSvc.OutputLevel) to cut down on the printout. To stop the dumping of events to the output remove the DumpMC from theApp.TopAlg. わからない！！ ＜＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝ > athena.py myGenOptions.py > athena_gen.out -rw-r--r-- 1 kondo zp 858721 May 11 19:16 athena_gen1.out -rw-r--r-- 1 kondo zp 242819 May 11 19:16 pythia.pool.root が確かにできた！！ Pythiaの出力例とその読解努力はPythia出力例.txtにある。 PDGのページにMC用の粒子番号つけ もあることがわかった。