1-hello |
2-mod |
3-pack |
4-pythia |
5-ATLFAST |
6-aodread |
7-sim |
8-dig |
9-rec |
10-aod |
今までやったことを復習してみる:
----------------------------------------------------------
run> bsub -q 1nh myjob
----------------------------------------------------------
> athena.py AthExHelloWorld_jobOptions.py を実行した
-----------------------------------------------------------
> athena.py MyJobOptions.py
-----------------------------------------------------------
run> get_files PDGTABLE.MeV
run> get_files jobOptions.pythia.py
> athena.py jobOptions.pythia.py > athena_gen.out
> athena.py myGenOptions.py > athena_gen.out
-rw-r--r-- 1 kondo zp 242819 May 11 19:16 pythia.pool.root
------------------------------------------------------------
> get_files POOLtoAtlfasttoAOD.py
> get_files Atlfast_ReadPOOLFile.py
> athena.py POOLtoAtlfasttoAOD.py > athena_fast.out
> get_files POOLtoAtlfasttoCBNT.py
> athena.py POOLtoAtlfasttoCBNT.py > athena_fastCBNT.out
------------------------------------------------------------
Obtain the following file:
> get_files AnalysisSkeleton_topOptions.py
EventSelector.InputCollections = [ "AOD.pool.root" ]
この内容を見ると
# Author: Ketevi A. Assamagan
# BNL, June 12, 2004
# get a handle on the ServiceManager which holds all the services
from AthenaCommon.AppMgr import ServiceMgr
# Event selector
import AthenaPoolCnvSvc.ReadAthenaPool
#EventSelector.BackNavigation = True
# Particle Properties
from PartPropSvc.PartPropSvcConf import PartPropSvc
# the POOL converters
include( "ParticleBuilderOptions/ESD_PoolCnv_jobOptions.py" )
include( "ParticleBuilderOptions/AOD_PoolCnv_jobOptions.py")
include( "ParticleBuilderOptions/McAOD_PoolCnv_jobOptions.py")
include( "EventAthenaPool/EventAthenaPool_joboptions.py" )
# The AOD input file
ServiceMgr.EventSelector.InputCollections = [ "AOD.pool.root" ]
# Athena-Aware NTuple making Tools
CBNTAthenaAware = True
include ("CBNT_Athena/CBNT_AthenaAware_jobOptions.py")
include ("CBNT_Athena/CBNT_EventInfo_jobOptions.py")
# list of the algorithms to be executed at run time
from UserAnalysis.UserAnalysisConf import AnalysisSkeleton
CBNT_AthenaAware += AnalysisSkeleton()
############# The properties of the AnalysisSkeleton Algorithm
CBNT_AthenaAware.AnalysisSkeleton.MCParticleContainer = "SpclMC"
CBNT_AthenaAware.AnalysisSkeleton.ElectronContainer = "ElectronAODCollection"
CBNT_AthenaAware.AnalysisSkeleton.DeltaRMatchCut = 0.2
CBNT_AthenaAware.AnalysisSkeleton.MaxDeltaR = 0.9999
CBNT_AthenaAware.AnalysisSkeleton.ElectronEtCut = 10.0*GeV
CBNT_AthenaAware.AnalysisSkeleton.ElectronEtaCut = 2.5
CBNT_AthenaAware.AnalysisSkeleton.ElectronCone = 0.9
CBNT_AthenaAware.AnalysisSkeleton.OutputLevel = DEBUG
##########################################
# setup TTree registration Service
# save ROOT histograms and Tuple
from GaudiSvc.GaudiSvcConf import THistSvc
ServiceMgr += THistSvc()
ServiceMgr.THistSvc.Output = [ "AANT DATAFILE='AnalysisSkeleton.aan.root' OPT='R
ECREATE'" ]
from AnalysisTools.AnalysisToolsConf import AANTupleStream
topSequence += AANTupleStream()
topSequence.AANTupleStream.ExtraRefNames = [ "StreamESD","Stream1" ]
topSequence.AANTupleStream.OutputName = 'AnalysisSkeleton.aan.root'
topSequence.AANTupleStream.WriteInputDataHeader = True
topSequence.AANTupleStream.OutputLevel = WARNING
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
ServiceMgr.MessageSvc.OutputLevel = ERROR
# Number of Events to process
theApp.EvtMax = 500
###################### For interactive analysis
#include ("PyAnalysisCore/InitPyAnalysisCore.py")
###################### Detail time measurement and auditors
# Use auditors
#theApp.AuditAlgorithms=True
#theApp.Dlls += [ "GaudiAud" ]
from GaudiSvc.GaudiSvcConf import AuditorSvc
ServiceMgr.AuditorSvc.Auditors += [ "ChronoAuditor"]
AthenaPoolCnvSvc = Service("AthenaPoolCnvSvc")
AthenaPoolCnvSvc.UseDetailChronoStat = TRUE
#StoreGateSvc = Service( "StoreGateSvc" )
#StoreGateSvc.Dump = True
#MessageSvc.OutputLevel = DEBUG
> athena.py AnalysisSkeleton_topOptions.py > athena_analysis.out
-rwxr-xr-x 1 kondo zp 2846 May 11 20:34 AnalysisSkeleton_topOptions.py*
-rw-r--r-- 1 kondo zp 489 May 11 20:39 PoolFileCatalog.xml.BAK
-rw-r--r-- 1 kondo zp 11943 May 11 20:39 athena_analysis.out
-rw-r--r-- 1 kondo zp 11201 May 11 20:39 AnalysisSkeleton.aan.root
------------------Using Root -------------------------
> root
displayがセットされない! ーーー> ここでROOTのお勉強
root [0] b = new TBrowser();
Click on the run folder on the left (2nd down)
Double click on AnalysisSkeleton.aan.root
Click on ROOT files
Double click on AnalysisSkeleton.aan.root
Double click on Electron
Double click on each histogram in turn
開いたら電子がない???
100 Z->ee AOD events is available at CERN at:
/afs/cern.ch/user/l/lloyd/public/zee100_13.0.30_aod.pool.root
You can either read this directly on lxplus or copy it to your own space/institute. These events are not part of the official production and should only be used for testing.
scp kondo@lxplus.cern.ch:/afs/cern.ch/user/l/lloyd/public/zee100_13.0.30_aod.pool.root .
で持ってきてみたらいろいろeventや分布はある。
vector.m_momemtum.m_eなど見ると分布している。306 events, 250GeVまでtailあり。
10 eventはできている。20に増やすぞ!
でもelec_pt分布は0エントリーだ。RUNとかBCIDとか20が出ている。しかしElectronPtの
プロットが出ない。