Running the BaBar Monte Carlo

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

Informations

• Publié par: Shean
• Nombre de lectures: 16
• Langue: English

Extrait

Ł Ł Ł

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