Running the BaBar Monte Carlo

Running the BaBar Monte Carlo

English
18 Pages
Read
Download
Downloading requires you to have access to the YouScribe library
Learn all about the services we offer

Description

Running the BaBar Monte CarloEvent GeneratorsAbi SofferColorado StateBaBar Tutorial11/6/2000• The big picture• Basics of running the generators• Changing things:– Decay modes– Beam parameters– Filtering for interesting eventsThe Big PictureGenFwkInt Bbsim SimApp BearEvtGen Geant3 Detector Reco(jetset) Interactions Responsebkqedkoralbetc.GeneratorsQA BbsimModuleApp AppGenFwkInt Bbsim SimApp BearEvtGen Geant3 Detector Reco(jetset) Interactions Responsebkqedkoralbetc.xdr file objy objyuse & loseBasics of Running BaBar Softwarehttp://www.slac.stanford.edu/BFROOT/www/Physics/Tools/generators/GenFwkInt/GeneratorsQA.htmlCreate a test-release directory:> newrel -s $MYWORK -t 9.8.0 980-genscratch disk: releasedir name$BFROOT/work/#/id> cd 980-gen> srtpath (Set the needed env. variables)> addpkg workdir (Get a work directory)> gmake workdir.setup (Set up workdir soft links)Generating Basics> cd workdir> sp3analboot (Read beam params from conditions DB)> GeneratorsQAApp> source PARENT/GeneratorsQA/EvtGen.tcl> ev begin -nev 1000Koralb,> exitBkqed,# ofetc.eventsPrintout for first 10 (or any other number of) events :StdHepPrint:StdHep Track info for event 1 :StdHepPrint:Trk# Stat Id Dtr1 DtrN Mom1 MomN Px Py Pz E Vx Vy VzStdHepPrint: 1 3 e- 3 3 0 0 -0.1651 -0.0010 8.9870 8.9886 0.1040 0.3304 -0.6490StdHepPrint: 2 3 e+ ...

Subjects

Informations

Published by
Reads 16
Language English
Report a problem
Ł Ł Ł
Running the BaBar Monte Carlo Event Generators
The big picture Basics of running the generators Changing things: – Decay modes – Beam parameters – Filtering for interesting events
Abi Soffer Colorado State BaBar Tutorial 11/6/2000
GenFwkInt
EvtGen (jetset)
bkqed
koralb
etc.
The Big Picture
Bbsim
Geant3 Interactions
SimApp
Detector Response
Bear
Reco
GeneratorsQA App
GenFwkInt
EvtGen (jetset)
bkqed
koralb
etc.
use & lose
BbsimModule App
Bbsim
Geant3 Interactions
SimApp
Detector Response
xdr file
Bear
Reco
objy
objy
Basics of Running BaBar Software
http://www.slac.stanford.edu/BFROOT/www/Physics/Tools/generators/G enFwkInt/GeneratorsQA.html Create a test-release directory: > newrel -s $MYWORK   -t 9.8.0   980-gen
scratch disk: release dir name $BFROOT/work/#/id
> cd 980-gen > srtpath <cr twice> (Set the needed env. variables) > addpkg workdir (Get a work directory) > gmake workdir.setup (Set up workdir soft links)
> > > > > >
Generating Basics
cd workdir sp3analboot (Read beam params from conditions DB) GeneratorsQAApp source PARENT/GeneratorsQA/ EvtGen .tcl ev begin -nev 1000 exit
# of events
Koralb, Bkqed, etc.
Printout for first 10 (or any other number of) events :
StdHepPrint:StdHep Track info for event 1 : StdHepPrint:Trk# Stat Id Dtr1 DtrN Mom1 MomN Px Py Pz E Vx Vy Vz StdHepPrint: 1 3 e- 3 3 0 0 -0.1651 -0.0010 8.9870 8.9886 0.1040 0.3304 -0.6490 StdHepPrint: 2 3 e+ 3 3 0 0 0.0581 0.0010 -3.1089 3.1094 0.1040 0.3304 -0.6490 StdHepPrint: 3 2 Upsilon(4S) 4 5 1 2 -0.1070 -0.0000 5.8782 12.0980 0.1040 0.3304 -0.6490 StdHepPrint: 4 2 B0 6 11 3 3 0.1446 -0.1457 3.0860 6.1186 0.1040 0.3304 -0.6490 StdHepPrint: 5 2 B0 24 27 3 3 -0.2516 0.1457 2.7922 5.9794 0.1040 0.3304 -0.6490 StdHepPrint: 6 2 rho0 12 13 4 4 0.0473 -0.0806 1.7868 1.9766 0.1070 0.3274 -0.5856 StdHepPrint: 7 1 K+ 0 0 4 4 0.2306 0.0134 0.0868 0.5519 0.1070 0.3274 -0.5856 StdHepPrint: 8 2 anti-K0 14 14 4 4 -0.0838 0.1379 0.2517 0.5806 0.1070 0.3274 -0.5856 StdHepPrint: 9 2 rho- 15 16 4 4 -0.1249 0.2401 0.4044 0.8257 0.1070 0.3274 -0.5856 StdHepPrint: 10 2 pi0 19 20 4 4 0.0347 0.0412 0.1547 0.2122 0.1070 0.3274 -0.5856 StdHepPrint: 11 2 anti-D0 21 23 4 4 0.0406 -0.4976 0.4016 1.9715 0.1070 0.3274 -0.5856 StdHepPrint: 12 1 pi+ 0 0 6 6 0.1607 0.0334 1.7467 1.7599 0.1070 0.3274 -0.5856
etc.
Upsilon(4S) | B0 -----------------------------------------------------B0 _  |_ rho0 K+ ---anti-K0 rho- ------pi0 anti-D0 |_ D*- ---------------------a_0+ -------------rho- pi+  |_ pi+ pi- |_ K_L0 |_ pi- pi0 |_ K+ e- anti-nu_e |_ anti-D0 pi- |_ eta pi+ |_ pi- pi0  |_ K*+ mu- anti-nu_mu |_ pi- pi+ pi0  |_ K+ pi0
Many histograms and ntuples:
PAW > hi/file 1 GfiMCTruth.hbook PAW > nt/lis
===> Directory :  3784 (N) collision Tuple  2 (N) collision generator  6 (N) GqaMCAnalysis Tuple  1 (N) collisions PAW > hi/lis
 ===> Directory :  987 (1) collision Pz  1340 (1) D+ Momentum in B0 decays  1341 (1) D- Momentum in B0 decays  1342 (1) D0 Momentum in B0 decays  1343 (1) DB Momentum in B0 decays  3783 (1) tau->K cos helicity  etc.
Using the Events
> cd ../ > addpkg GeneratorsQA (get package to make changes) See GeneratorsQA/GqaMCAnalysis.cc , where the ntuples & histograms are made:
HepAList<BtaCandidate>*  mcList ; getTmpAList (anEvent, mcList , _btaTruthList.value()); HepAListIterator<BtaCandidate>  iter (*( mcList )); BtaCandidate  cand ; while ( cand = iter ()){  // do things with cand : extract quantities, histogram, etc. }
> gmake lib GeneratorsQA.bin
Changing What™s Generated
> source PARENT/GeneratorsQA/ XxxXxx .tcl
Ł EvtGen (& Jetset) Ł Koralb Ł Bkqed Ł Bhwide Ł GamGam Ł Diag36 Ł SingleParticle  Ł TurtleRead Ł Hemicosm  Ł StdHepRead
hadronic Ν # Ν % m # m % e + e % 2-photon 4-lepton single particles Turtle background cosmics Reads ascii files
Changing Generator Parameters
> module talk GfiXxxXxx GfiXxxXxx>  help > exit
All commands of all generator modules are listed at http://www.slac.stanford.edu/BFROOT/www/Physics/Tools/generator s/GenFwkInt/GenFwkInt.html
Ł eg.: Generating continuum/BBbar: > module talk GfiEvtGen GfiEvtGen> GENERATE set continuum > exit
or Upsilon(4S)
Changing Decay Modes in EvtGen
eg.: B 0 # ϑ % B 0 | generic: > module talk GfiEvtGen GfiXxxXxx> UDECAY set BtoPiPi.dec > exit
BtoPiPi.dec contains (leaving out some details) Decay B0 0.7 Enddecay End
pi+ pi %
Br. Fraction
Final state
PHSP ;
Decay model
Hundreds of examples in package ProdDecayFiles or EvtGen/DECAY.DEC