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 ------------------------------------------------------------

(10)Reading AOD with Athena (2008.5.11卒業)

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の プロットが出ない。