* ****************************************************************** * General Macro for Monte Carlo Production and GEANT simulation * * for HERA-B Physics * * * * Author: T. Lohse * * Date: 1999-2001 * * * * List of modifications: * * - December 27, 2002 (V. Saveliev) * * adopted new MC naming convention for paths * * - May 1, 2003 (U. Husemann): * * added new open charm channels: 240xx 241xx 242xx 243xx 244xx * * - July 30, 2003 (U. Husemann): * * o updated open charm channels (240xx 241xx 242xx 243xx 244xx): * * -> allow c.c. channels with same process ID (A. Gorisek) * * o added inelastic channels 1260x (K0S), 1270x (Lambda) * * o added inelastic channels 1880x, 1890x, 1900x (phi, K*) * * o added new group of multistrange baryons (Xi, Omega) * * 8010x, 8020x, 8030x, 8040x (works with ARTE >= ARTE-04-01-r5)* * o added new NRQCD charmonium channels (J/psi, Psi(2S)) * * 4220x, 4300x, 4310x, 4320x * * o added commands to adjust beam size and position * * (works with ARTE >= ARTE-04-01-r5) * * o updated decay lists according to J. Spenglers new TDEC table * * * * ****************************************************************** * * *** Script for variables *** * online = on arte_version = 0401r3 geo_vers_1 = '02' geo_vers_2 = '1205' detector = 1 magnet = on wiremask = '00010000' * interactive = on * e_beam = 920. phys_proc = 41000 statistics = 500 * production_site = desy mc_outdir = ./ run_number = '00001' file_number = '00001' * geanformat = 'G' geomformat = 'G' recoformat = 'G' * if [online]=on then interactive = off statistics = 1000000 production_site = desy phys_proc = $ENV(MCP) geo_vers_1 = $ENV(MCY) geo_vers_2 = $ENV(MCV) wiremask = $ENV(MCW) endif * *############################################################################# * detector Defines the detector geometry to be used: * ======== * 0 = generator level; * 1 = full geometry; * * for detector>0 targets are switched on/off according to: * * wiremask Defines unused target wires --> abcdefgh * ======== * a => I inner 1=in 0=out * b => I outer 1=in 0=out * c => I upper 1=in 0=out * d => I lower 1=in 0=out * e => II inner 1=in 0=out * f => II outer 1=in 0=out * g => II upper 1=in 0=out * h => II lower 1=in 0=out * * Note: In addition (!) only target of correct material type * are considered (last digit of physproc) * *############################################################################# * *... physproc Defines the physics process to be generated: * *---------------------------------------------------------------------- * * Overview: group 1: Normal Inelastic Proton Nucleus Interactions * --------- group 2: Open Charm Production * group 3: Open Beauty Production (incl. CP Channels) * group 4: Charmonium Production * group 5: Bottomonium Production * group 6: Rare Decays of Charm Hadrons * group 7: Rare Decays of Beauty Hadrons * * note: physproc -> -physproc for charge conjugated special decays * of heavy quark pairs * note: In the following 'g g' represents any partonic * initial state relevant for heavy quark antiquark production * *---------------------------------------------------------------------- * group 1: Normal Inelastic Proton Nucleus Interactions * -------- * 1000j: no pT-cut, unweighted events j=0: p-C * j=1: p-Al * j=2: p-Ti * j=3: p-W * j=9: p-any * 1010j: pT-cut at 2 GeV, unweighted events * 1020j: pT-cut at 3 GeV, unweighted events * 1110j: pT-threshold at 2 GeV, weighted events * 1120j: pT-threshold at 4 GeV, weighted events * 1130j: pT-threshold at 6 GeV, weighted events * 1200j: pN -> n(rho0) X -> n(mu+ mu-) X * 1210j: pN -> n(rho0) X -> n(mu+ mu- pi0) X * 1220j: pN -> n(omega) X -> n(mu+ mu-) X * 1230j: pN -> n(omega) X -> n(mu+ mu- pi0) X * 1240j: pN -> n(phi) X -> n(mu+ mu-) X * 1250j: pN -> n(phi) X -> n(mu+ mu- pi0) X * 1260j: pN -> n(K0S) X * 1270j: pN -> n(Lambda/Lambda~) X * 1510j: Drell-Yan q\bar{q} -> e+e- * 1520j: Drell-Yan q\bar{q} -> mu+mu- * 1530j: Drell-Yan NLO: q\bar{q}->Gamma* g->e+e-g / qg->Gamma* q->e+ e- q * 1540j: Drell-Yan NLO: q\bar{q}->Gamma* g->e+e-g / qg->Gamma* q->mu+mu-q * 1880j: pN -> n(phi) X -> n(K+ K-) X * 1890j: pN-> n(K*0) X -> n(K+ pi-) X * 1900j: pN-> n(K*0~) X -> n(K- pi+) X * * ---------------------------------------------------------------------- * group 2: Open Charm Production * -------- * 200kj: g g -> c cbar j=0: p-C * j=1: p-Al * j=2: p-Ti * j=3: p-W * j=9: p-any * fragmentation modes: * k=1: only particle containing the heavy quark * k=2: half string fragmentation from heavy quarks * k=3: independent fragmentation of heavy quarks * k=4: full string fragmentation from heavy quarks * k=5: all fragments from pN strings kept * 201kj: g g -> c cbar -> mu X * 202kj: g g -> c cbar -> e X * 203kj: g g -> c cbar -> mu+ numu mu- numu~ X * 204kj: g g -> c cbar -> e+ nue e- nue~ X * 205kj: as 203kj with pT(mu1)>1.0 and pT(mu2)>1.0 * 206kj: as 204kj with pT(e1)>1.0 and pT(e2)>1.0 * 210kj: g g -> c cbar -> chadr~ D*+ X -> l+ nu D0 pi+ X * -> l+ nu pi+ l- nu~ X * 211kj: as 210kj with pT(lept1)>1.0 and pT(lept2)>1.0 * 212kj: as 210kj with mass(lept1,lept2)>2.5 * 213kj as 210kj with muons in final state * 214kj as 213kj with pT(lept1)>1.0 and pT(lept2)>1.0 * 215kj: as 213kj with mass(lept1,lept2)>2.5 * 220kj: g g -> c cbar -> chadr~ D*+ X -> l+ nu D0 pi+ X * -> l+ nu pi+ pi+ K- X * 221kj: as 220kj, pT(lept)>1.0GeV * 240kj: g g -> c cbar -> D0 X -> K- pi+ X * 241kj: g g -> c cbar -> D+ X -> K- pi+ pi+ * (incl. resonant decays K*~ -> K- pi+, K1430_0*~ -> K- pi+) * 242kj: g g -> c cbar -> D_s+ X -> K+ K- pi+ * (incl. resonant decays phi -> K+ K- , K*~ -> K- pi+, * K1430_0*~ -> K- pi+, f_0 -> K+ K-) * 243kj: g g -> c cbar -> D*+ X -> D0 pi+ X -> K- pi+ pi+ * 244kj: g g -> c cbar -> D+ X -> K_S pi+ (K_S decay in GEANT) * * ---------------------------------------------------------------------- * group 3: Open Beauty Production (incl. CP Channels) * -------- * * 300kj: g g -> b bbar j=0: p-C * j=1: p-Al * j=2: p-Ti * j=3: p-W * j=9: p-any * fragmentation modes: * k=1: only particle containing the heavy quark * k=2: half string fragmentation from heavy quarks * k=3: independent fragmentation of heavy quarks * k=4: full string fragmentation from heavy quarks * k=5: all fragments from pN strings kept * 301kj: g g -> b bbar -> c X * 302kj: g g -> b bbar -> u X * 303kj: g g -> b bbar -> mu X * 304kj: g g -> b bbar -> e X * 305kj: g g -> b bbar -> J/Psi X -> mu+ mu- X * 306kj: g g -> b bbar -> J/Psi X -> e+ e- X * 307kj: g g -> b bbar -> B0 X * 308kj: g g -> b bbar -> B0s X * 309kj: g g -> b bbar -> Lambda_b X * 310kj: g g -> b bbar -> B0 bhadr X -> l nu l nu X * 311kj: g g -> b bbar -> B**+ X -> B0 pi+ X -> J/Psi K0S pi+ X * -> mu+ mu- K0S pi+ X * 312kj: g g -> b bbar -> B**+ X -> B0 pi+ X -> J/Psi K0S pi+ X * -> e+ e- K0S pi+ X * 320kj: g g -> b bbar -> B0 X -> J/Psi K0S X -> mu+ mu- K0S X * 321kj: g g -> b bbar -> B0 X -> J/Psi K0S X -> e+ e- K0S X * 322kj: g g -> b bbar -> B0 X -> J/Psi K*0 X -> mu+ mu- K0S pi0 X * 323kj: g g -> b bbar -> B0 X -> J/Psi K*0 X -> e+ e- K0S pi0 X * 324kj: g g -> b bbar -> B0 X -> J/Psi K*0 X -> mu+ mu- K+ pi- X * or B0 -> B0bar -> J/Psi K*0bar -> mu+ mu- K- pi+ * 325kj: g g -> b bbar -> B0 X -> J/Psi K*0 X -> e+ e- K+ pi- X * or B0 -> B0bar -> J/Psi K*0bar -> e+ e- K- pi+ * 326kj: g g -> b bbar -> B0 X -> J/Psi rho0 X -> mu+ mu- pi+ pi- X * 327kj: g g -> b bbar -> B0 X -> J/Psi rho0 X -> e+ e- pi+ pi- X * 328kj: g g -> b bbar -> B0 X -> chi_c1 K0S X -> mu+ mu- gamma K0S X * 329kj: g g -> b bbar -> B0 X -> chi_c1 K0S X -> e+ e- gamma K0S X * 340kj: g g -> b bbar -> B0 X -> pi+ pi- X * 341kj: g g -> b bbar -> B0 X -> K+ pi- X or B0 -> B0bar -> K- pi+ * 342kj: g g -> b bbar -> B0 X -> a1+- pi-+ X -> pi+ pi- pi+ pi- X * 360kj: g g -> b bbar -> B0s X -> Ds pi X * 361kj: g g -> b bbar -> B0s X -> Ds K X * 362kj: g g -> b bbar -> B0s X -> K+ K- X * 363kj: g g -> b bbar -> B0s X -> J/Psi Phi X -> mu+ mu- K+ K- X * 364kj: g g -> b bbar -> B0s X -> J/Psi Phi X -> e+ e- K+ K- X * 370kj: g g -> b bbar -> bhadr B0 X -> l nu D* l nu X * -> l nu D0 pi l nu X * 371kj: g g -> b bbar -> bhadr B0s X -> l nu Ds l nu X * * ---------------------------------------------------------------------- * group 4: Charmonium Production * -------- * * 4000j: g g -> J/Psi X, weighted events j=0: p-C * j=1: p-Al * j=2: p-Ti * j=3: p-W * j=9: p-any * 4010j: g g -> (J/Psi or chi_c) X, weighted events * 4100j: g g -> J/Psi X -> mu+ mu- X, weighted events * 4110j: g g -> J/Psi X -> e+ e- X, weighted events * 4120j: g g -> chi_c1 X -> J/Psi gamma X -> mu+ mu- gamma X, * weighted events * 4130j: g g -> chi_c1 X -> J/Psi gamma X -> e+ e- gamma X, * weighted events * 4200j: g g (q qbar) -> charmonium via Color Singlet processes * 4210j: g g (q qbar) -> charmonium via NRQCD, i.e. CSM+COM channels * 4220j: g g (q qbar) -> J/Psi X -> mu+ mu- X (NRQCD) * feed-down from Psi(2S), chi_c0,1,2 * Note: The following channels work only with ARTE >= ARTE-04-01-r5 * To make them work with older versions, replace 'LUND PSI2S' * by 'LUND COM -psi2S' (very inefficient!) * 4300j: g g (q qbar) -> Psi(2S) X -> mu+ mu- X (NRQCD) * 4310j: g g (q qbar) -> Psi(2S) X -> J/Psi pi+ pi- X * -> mu+ mu- pi+ pi- X (NRQCD) * 4320j: g g (q qbar) -> Psi(2S) X -> Jpsi eta X * -> mu+ mu- gamma gamma (NRQCD) * ---------------------------------------------------------------------- * group 5: Bottomonium Production * -------- * * 5000j: g g -> Upsilon(1S) X, weighted events j=0: p-C * j=1: p-Al * j=2: p-Ti * j=3: p-W * j=9: p-any * 5010j: g g -> (Upsilon(1S) or chi_b(1P)) X, weighted events * 5100j: g g -> Upsilon(1S) X -> mu+ mu- X, weighted events * 5110j: g g -> Upsilon(1S) X -> e+ e- X, weighted events * 5120j: g g -> chi_b(1P) X -> Upsilon(1S) gamma X -> mu+ mu- gamma X, * weighted events * 5130j: g g -> chi_b(1P) X -> Upsilon(1S) gamma X -> e+ e- gamma X, * weighted events * * ---------------------------------------------------------------------- * group 6: Rare Decays of Charm Hadrons * -------- * * 600kj: g g -> D0 X -> mu+ mu- X j=0: p-C * j=1: p-Al * j=2: p-Ti * j=3: p-W * j=9: p-any * fragmentation modes: * k=1: only particle containing the heavy quark * k=2: half string fragmentation from heavy quarks * k=3: independent fragmentation of heavy quarks * k=4: full string fragmentation from heavy quarks * k=5: all fragments from pN strings kept * 601kj: g g -> D0 X -> mu+ mu- gamma X * 602kj: g g -> D0 X -> e+ e- gamma X * 610kj: g g -> D0 X -> gamma gamma X * *---------------------------------------------------------------------- * group 7: Rare Decays of Beauty Hadrons * -------- * * 700kj: g g -> B0 X -> mu+ mu- X j=0: p-C * j=1: p-Al * j=2: p-Ti * j=3: p-W * j=9: p-any * fragmentation modes: * k=1: only particle containing the heavy quark * k=2: half string fragmentation from heavy quarks * k=3: independent fragmentation of heavy quarks * k=4: full string fragmentation from heavy quarks * k=5: all fragments from pN strings kept * 701kj: g g -> B0s X -> mu+ mu- X * 702kj: g g -> B0 X -> mu+ mu- gamma X * 703kj: g g -> B0 X -> e+ e- gamma X * 704kj: g g -> B0s X -> mu+ mu- gamma X * 705kj: g g -> B0s X -> e+ e- gamma X * 710kj: g g -> B0 X -> gamma gamma X * 711kj: g g -> B0s X -> gamma gamma X * 720kj: g g -> B0 X -> K0S mu+ mu- X * 721kj: g g -> B0 X -> K0S e+ e- X * 722kj: g g -> B+ X -> K+ mu+ mu- X * 723kj: g g -> B+ X -> K+ e+ e- X * 724kj: g g -> B0 X -> K*0 mu+ mu- X -> K+ pi- mu+ mu- X * or B0 -> B0bar -> K*0bar mu+ mu- -> K- pi+ mu+ mu- * 725kj: g g -> B0 X -> K*0 e+ e- X -> K+ pi- e+ e- X * or B0 -> B0bar -> K*0bar e+ e- -> K- pi+ e+ e- * 726kj: g g -> B+ X -> K*+ mu+ mu- X -> K0S pi+ mu+ mu- X * 727kj: g g -> B+ X -> K*+ e+ e- X -> K0S pi+ e+ e- X * 730kj: g g -> B0s X -> phi gamma X -> K+ K- gamma X * 740kj: g g -> Lambda_b X -> Lambda gamma X * 741kj: g g -> Lambda_b X -> Lambda mu+ mu- X * 742kj: g g -> Lambda_b X -> Lambda e+ e- X * * ---------------------------------------------------------------------- * group 8: Multistrange Baryons * -------- * * Note: The requested decay channels will only work with ARTE versions * >= ARTE-04-01-r5 (decay simulated in GEANT, not in physics generator!) * * j=0: p-C * j=1: p-Al * j=2: p-Ti * j=3: p-W * j=9: p-any * 8010j: g g -> Xi- X -> Lambda pi- X * 8020j: g g -> Xi*0(1530) X -> Xi- pi+ X -> Lambda pi- pi+ X * 8030j: g g -> Xi*-(1530) X -> Xi- pi0 X -> Lambda pi- pi0 X * 8040j: g g -> Omega- X -> further decays according to PDG * *############################################################################# * *... MC production site * site = '00' if ([production_site]=amsterdam .or. [production_site]=1) then site = '01' elseif ([production_site]=austin .or. [production_site]=2) then site = '02' elseif ([production_site]=beijing .or. [production_site]=3) then site = '03' elseif ([production_site]=berlin .or. [production_site]=4) then site = '04' elseif ([production_site]=bologna .or. [production_site]=5) then site = '05' elseif ([production_site]=dortmund .or. [production_site]=6) then site = '06' elseif ([production_site]=northwestern .or. [production_site]=7) then site = '07' elseif ([production_site]=glasgow .or. [production_site]=8) then site = '08' elseif ([production_site]=desy .or. [production_site]=9) then site = '09' elseif ([production_site]=uni_hamburg .or. [production_site]=10) then site = '10' elseif ([production_site]=mpi_heidelb .or. [production_site]=11) then site = '11' elseif ([production_site]=uni_heidelb .or. [production_site]=12) then site = '12' elseif ([production_site]=helsinki .or. [production_site]=13) then site = '13' elseif ([production_site]=kiev .or. [production_site]=14) then site = '14' elseif ([production_site]=nbi .or. [production_site]=15) then site = '15' elseif ([production_site]=ljubljana .or. [production_site]=16) then site = '16' elseif ([production_site]=ucla .or. [production_site]=17) then site = '17' elseif ([production_site]=lund .or. [production_site]=18) then site = '18' elseif ([production_site]=itep_moscow .or. [production_site]=19) then site = '19' elseif ([production_site]=mephi_moscow .or. [production_site]=20) then site = '20' elseif ([production_site]=munich .or. [production_site]=21) then site = '21' elseif ([production_site]=yale .or. [production_site]=22) then site = '22' elseif ([production_site]=princeton .or. [production_site]=23) then site = '23' elseif ([production_site]=roma .or. [production_site]=24) then site = '24' elseif ([production_site]=siegen .or. [production_site]=25) then site = '25' elseif ([production_site]=brookhaven .or. [production_site]=26) then site = '26' elseif ([production_site]=zeuthen .or. [production_site]=27) then site = '27' else message 'unknown production site!!!!!' endif * * *** Physics Channels Generator *** * *... Generator with multyplicity tuning for HARD Processes * /GENERATOR/REWEIGHT 3 * *... Phys Processes * charge = [phys_proc] phys_proc = abs([phys_proc]) itarg = mod([phys_proc],10) * iphys = ([phys_proc]-[itarg])/10 group = ([iphys]-mod([iphys],1000))/1000 fragm = mod([iphys],10) iphys = ([iphys]-[fragm])/10 if ([group]=2 .or. [group]=3 .or. [group]=6 .or. [group]=7) then QQ [fragm] endif * goto [iphys] * 100: LUND inel goto 9999 101: INEL 'PTC' 2. LUND inel goto 9999 102: INEL 'PTC' 3. LUND inel goto 9999 111: INEL 'PTT' 2. LUND inel goto 9999 112: INEL 'PTT' 4. LUND inel goto 9999 113: INEL 'PTT' 6. LUND inel goto 9999 120: LUND inel particle rho. BEGdecay rho. decay 1 0 1. mu+ mu- ENDdecay rho. goto 9999 121: LUND inel particle rho. BEGdecay rho. decay 1 0 1. mu+ mu- pi. ENDdecay rho. goto 9999 122: LUND inel particle omega BEGdecay omega decay 1 0 1. mu+ mu- ENDdecay omega goto 9999 123: LUND inel particle omega BEGdecay omega decay 1 0 1. mu+ mu- pi. ENDdecay omega goto 9999 124: LUND inel particle phi BEGdecay phi decay 1 0 1. mu+ mu- ENDdecay phi goto 9999 125: LUND inel particle phi BEGdecay phi decay 1 0 1. mu+ mu- pi. ENDdecay phi goto 9999 126: LUND inel particle K_S goto 9999 127: if [charge]>0 then spec = Lam. else spec = Lam~. endif LUND inel particle [spec] goto 9999 130: LUND inel INEL 'PTS' 1. particle rho. BEGdecay rho. decay 1 0 1. mu+ mu- ENDdecay rho. goto 9999 131: LUND inel INEL 'PTS' 1. particle rho. BEGdecay rho. decay 1 0 1. mu+ mu- pi. ENDdecay rho. goto 9999 132: LUND inel INEL 'PTS' 1.2 particle omega BEGdecay omega decay 1 0 1. mu+ mu- ENDdecay omega goto 9999 133: LUND inel INEL 'PTS' 1.4 particle omega BEGdecay omega decay 1 0 1. mu+ mu- pi. ENDdecay omega goto 9999 134: LUND inel INEL 'PTS' 1. particle phi BEGdecay phi decay 1 0 1. mu+ mu- ENDdecay phi goto 9999 135: LUND inel INEL 'PTS' 1. particle phi BEGdecay phi decay 1 0 1. mu+ mu- pi. ENDdecay phi goto 9999 136: LUND inel INEL 'PTS' 1.4 LEPTON 1 0.6 particle omega BEGdecay omega decay 1 0 1. mu+ mu- pi. ENDdecay omega goto 9999 151: LUND 'DY' PARTICLE e+ e- goto 9999 152: LUND 'DY' PARTICLE mu+ mu- goto 9999 153: LUND 'NLODY' particle e+ e- goto 9999 154: LUND 'NLODY' particle mu+ mu- goto 9999 188: LUND inel particle phi BEGdecay phi decay 1 3 1. K+ K- ENDdecay phi goto 9999 189: LUND inel particle K*. BEGdecay K*. decay 1 3 1. K+ pi- ENDdecay K*. goto 9999 190: LUND inel particle K*~ BEGdecay K*~ decay 1 3 1. K- pi+ ENDdecay K*~ goto 9999 * 200: LUND cc goto 9999 201: particle mu+ mu- LUND cc goto 9999 202: particle e+ e- LUND cc goto 9999 203: particle mu+ LUND cc -mu- BEGdecay D+ decay 1 42 0.20300 mu+ nu_mu K~ decay 1 42 0.18900 mu+ nu_mu K*~ decay 1 42 0.01450 mu+ nu_mu K~ pi. decay 1 42 0.01450 mu+ nu_mu K- pi+ decay 1 42 0.03200 mu+ nu_mu K*~ pi. decay 1 42 0.03200 mu+ nu_mu K*- pi+ decay 1 42 0.00290 mu+ nu_mu pi. decay 1 42 0.00290 mu+ nu_mu eta decay 1 42 0.00290 mu+ nu_mu eta` decay 1 42 0.00290 mu+ nu_mu rho. decay 1 42 0.00290 mu+ nu_mu omega ENDdecay D+ BEGdecay D. decay 1 42 0.22080 mu+ nu_mu K- decay 1 42 0.17530 mu+ nu_mu K*- decay 1 42 0.01300 mu+ nu_mu K~ pi- decay 1 42 0.01300 mu+ nu_mu K- pi. decay 1 42 0.02600 mu+ nu_mu K*~ pi- decay 1 42 0.02600 mu+ nu_mu K*- pi. decay 1 42 0.01300 mu+ nu_mu pi- decay 1 42 0.01300 mu+ nu_mu rho- ENDdecay D. BEGdecay D_s+ decay 1 42 0.12500 mu+ nu_mu eta decay 1 42 0.12500 mu+ nu_mu eta` decay 1 42 0.18750 mu+ nu_mu phi decay 1 42 0.03125 mu+ nu_mu K+ K- decay 1 42 0.03125 mu+ nu_mu K. K~ ENDdecay D_s+ BEGdecay Lam_c+ decay 1 42 0.20000 mu+ nu_mu Lam. decay 1 42 0.05550 mu+ nu_mu Sig. decay 1 42 0.05550 mu+ nu_mu Sig*. decay 1 42 0.03330 mu+ nu_mu n decay 1 42 0.02220 mu+ nu_mu Delta. decay 1 42 0.06660 mu+ nu_mu p pi- decay 1 42 0.06660 mu+ nu_mu n pi. ENDdecay Lam_c+ BEGdecay c-hadron decay 1 42 0.50000 mu+ nu_mu squark specflav ENDdecay c-hadron goto 9999 204: particle e+ LUND cc -e- BEGdecay D+ decay 1 42 0.20300 e+ nu_e K~ decay 1 42 0.18900 e+ nu_e K*~ decay 1 42 0.01450 e+ nu_e K~ pi. decay 1 42 0.01450 e+ nu_e K- pi+ decay 1 42 0.03200 e+ nu_e K*~ pi. decay 1 42 0.03200 e+ nu_e K*- pi+ decay 1 42 0.00290 e+ nu_e pi. decay 1 42 0.00290 e+ nu_e eta decay 1 42 0.00290 e+ nu_e eta` decay 1 42 0.00290 e+ nu_e rho. decay 1 42 0.00290 e+ nu_e omega ENDdecay D+ BEGdecay D. decay 1 42 0.22080 e+ nu_e K- decay 1 42 0.17530 e+ nu_e K*- decay 1 42 0.01300 e+ nu_e K~ pi- decay 1 42 0.01300 e+ nu_e K- pi. decay 1 42 0.02600 e+ nu_e K*~ pi- decay 1 42 0.02600 e+ nu_e K*- pi. decay 1 42 0.01300 e+ nu_e pi- decay 1 42 0.01300 e+ nu_e rho- ENDdecay D. BEGdecay D_s+ decay 1 42 0.12500 e+ nu_e eta decay 1 42 0.12500 e+ nu_e eta` decay 1 42 0.18750 e+ nu_e phi decay 1 42 0.03125 e+ nu_e K+ K- decay 1 42 0.03125 e+ nu_e K. K~ ENDdecay D_s+ BEGdecay Lam_c+ decay 1 42 0.20000 e+ nu_e Lam. decay 1 42 0.05550 e+ nu_e Sig. decay 1 42 0.05550 e+ nu_e Sig*. decay 1 42 0.03330 e+ nu_e n decay 1 42 0.02220 e+ nu_e Delta. decay 1 42 0.06660 e+ nu_e p pi- decay 1 42 0.06660 e+ nu_e n pi. ENDdecay Lam_c+ BEGdecay c-hadron decay 1 42 0.50000 e+ nu_e squark specflav ENDdecay c-hadron goto 9999 205: LEPTON 2 1. 1. 0. iphys = 203 goto 203 206: LEPTON 2 1. 1. 0. iphys = 204 goto 204 210: if [charge]>0 then spec = D*+ else spec = D*- endif LUND cc -[spec] BEGdecay D*+ decay 0 3 1. D. pi+ goto D*+dec 210D*+ret: ENDdecay D*+ BEGdecay D+ decay 1 42 0.20300 e+ nu_e K~ decay 1 42 0.18900 e+ nu_e K*~ decay 1 42 0.01450 e+ nu_e K~ pi. decay 1 42 0.01450 e+ nu_e K- pi+ decay 1 42 0.03200 e+ nu_e K*~ pi. decay 1 42 0.03200 e+ nu_e K*- pi+ decay 1 42 0.00290 e+ nu_e pi. decay 1 42 0.00290 e+ nu_e eta decay 1 42 0.00290 e+ nu_e eta` decay 1 42 0.00290 e+ nu_e rho. decay 1 42 0.00290 e+ nu_e omega decay 1 42 0.20300 mu+ nu_mu K~ decay 1 42 0.18900 mu+ nu_mu K*~ decay 1 42 0.01450 mu+ nu_mu K~ pi. decay 1 42 0.01450 mu+ nu_mu K- pi+ decay 1 42 0.03200 mu+ nu_mu K*~ pi. decay 1 42 0.03200 mu+ nu_mu K*- pi+ decay 1 42 0.00290 mu+ nu_mu pi. decay 1 42 0.00290 mu+ nu_mu eta decay 1 42 0.00290 mu+ nu_mu eta` decay 1 42 0.00290 mu+ nu_mu rho. decay 1 42 0.00290 mu+ nu_mu omega ENDdecay D+ BEGdecay D. decay 1 42 0.22080 e+ nu_e K- decay 1 42 0.17530 e+ nu_e K*- decay 1 42 0.01300 e+ nu_e K~ pi- decay 1 42 0.01300 e+ nu_e K- pi. decay 1 42 0.02600 e+ nu_e K*~ pi- decay 1 42 0.02600 e+ nu_e K*- pi. decay 1 42 0.01300 e+ nu_e pi- decay 1 42 0.01300 e+ nu_e rho- decay 1 42 0.22080 mu+ nu_mu K- decay 1 42 0.17530 mu+ nu_mu K*- decay 1 42 0.01300 mu+ nu_mu K~ pi- decay 1 42 0.01300 mu+ nu_mu K- pi. decay 1 42 0.02600 mu+ nu_mu K*~ pi- decay 1 42 0.02600 mu+ nu_mu K*- pi. decay 1 42 0.01300 mu+ nu_mu pi- decay 1 42 0.01300 mu+ nu_mu rho- ENDdecay D. BEGdecay D_s+ decay 1 42 0.12500 e+ nu_e eta decay 1 42 0.12500 e+ nu_e eta` decay 1 42 0.18750 e+ nu_e phi decay 1 42 0.03125 e+ nu_e K+ K- decay 1 42 0.03125 e+ nu_e K. K~ decay 1 42 0.12500 mu+ nu_mu eta decay 1 42 0.12500 mu+ nu_mu eta` decay 1 42 0.18750 mu+ nu_mu phi decay 1 42 0.03125 mu+ nu_mu K+ K- decay 1 42 0.03125 mu+ nu_mu K. K~ ENDdecay D_s+ BEGdecay Lam_c+ decay 1 42 0.20000 e+ nu_e Lam. decay 1 42 0.05550 e+ nu_e Sig. decay 1 42 0.05550 e+ nu_e Sig*. decay 1 42 0.03330 e+ nu_e n decay 1 42 0.02220 e+ nu_e Delta. decay 1 42 0.06660 e+ nu_e p pi- decay 1 42 0.06660 e+ nu_e n pi. decay 1 42 0.20000 mu+ nu_mu Lam. decay 1 42 0.05550 mu+ nu_mu Sig. decay 1 42 0.05550 mu+ nu_mu Sig*. decay 1 42 0.03330 mu+ nu_mu n decay 1 42 0.02220 mu+ nu_mu Delta. decay 1 42 0.06660 mu+ nu_mu p pi- decay 1 42 0.06660 mu+ nu_mu n pi. ENDdecay Lam_c+ BEGdecay c-hadron decay 1 42 0.50000 e+ nu_e squark specflav decay 1 42 0.50000 mu+ nu_mu squark specflav ENDdecay c-hadron channel [spec] 1 0 goto 9999 211: LEPTON 2 1. 1. 0. iphys = 210 goto 210 212: LEPTON 2 0. 0. 2.5 iphys = 210 goto 210 213: if [charge]>0 then spec = D*+ else spec = D*- endif LUND cc -[spec] BEGdecay D*+ decay 0 3 1. D. pi+ goto D*+dec 213D*+ret: ENDdecay D*+ BEGdecay D+ decay 1 42 0.40600 mu+ nu_mu K~ decay 1 42 0.37800 mu+ nu_mu K*~ decay 1 42 0.02900 mu+ nu_mu K~ pi. decay 1 42 0.02900 mu+ nu_mu K- pi+ decay 1 42 0.06400 mu+ nu_mu K*~ pi. decay 1 42 0.06400 mu+ nu_mu K*- pi+ decay 1 42 0.00580 mu+ nu_mu pi. decay 1 42 0.00580 mu+ nu_mu eta decay 1 42 0.00580 mu+ nu_mu eta` decay 1 42 0.00580 mu+ nu_mu rho. decay 1 42 0.00580 mu+ nu_mu omega ENDdecay D+ BEGdecay D. decay 1 42 0.44160 mu+ nu_mu K- decay 1 42 0.35060 mu+ nu_mu K*- decay 1 42 0.02600 mu+ nu_mu K~ pi- decay 1 42 0.02600 mu+ nu_mu K- pi. decay 1 42 0.05200 mu+ nu_mu K*~ pi- decay 1 42 0.05200 mu+ nu_mu K*- pi. decay 1 42 0.02600 mu+ nu_mu pi- decay 1 42 0.02600 mu+ nu_mu rho- ENDdecay D. BEGdecay D_s+ decay 1 42 0.25000 mu+ nu_mu eta decay 1 42 0.25000 mu+ nu_mu eta` decay 1 42 0.37500 mu+ nu_mu phi decay 1 42 0.06250 mu+ nu_mu K+ K- decay 1 42 0.06250 mu+ nu_mu K. K~ ENDdecay D_s+ BEGdecay Lam_c+ decay 1 42 0.40000 mu+ nu_mu Lam. decay 1 42 0.11100 mu+ nu_mu Sig. decay 1 42 0.11150 mu+ nu_mu Sig*. decay 1 42 0.06660 mu+ nu_mu n decay 1 42 0.04440 mu+ nu_mu Delta. decay 1 42 0.13320 mu+ nu_mu p pi- decay 1 42 0.13320 mu+ nu_mu n pi. ENDdecay Lam_c+ BEGdecay c-hadron decay 1 42 1.00000 mu+ nu_mu squark specflav ENDdecay c-hadron channel [spec] 1 0 goto 9999 * 214: LEPTON 2 1. 1. 0. iphys = 213 goto 213 * 215: LEPTON 2 0. 0. 2.5 iphys = 213 goto 213 * 220: if [charge]>0 then spec = D*+ else spec = D*- endif LEPTON 1. 0. 0. 0. LUND cc -[spec] BEGdecay D*+ decay 0 3 1. D. pi+ goto D*+dec 220D*+ret: ENDdecay D*+ BEGdecay D. decay 0 0 0.03650 K- pi+ goto D0dec 220D0ret: ENDdecay D. channel [spec] 1 0 channel D. 1 0 goto 9999 221: if [charge]>0 then spec = D*+ else spec = D*- endif LEPTON 1. 1. 0. 0. LUND cc -[spec] BEGdecay D*+ decay 0 3 1. D. pi+ goto D*+dec 221D*+ret: ENDdecay D*+ BEGdecay D. decay 0 0 0.03650 K- pi+ goto D0dec 221D0ret: ENDdecay D. channel [spec] 1 0 channel D. 1 0 goto 9999 240: if [charge]>0 then spec = D. else spec = D~ endif LUND cc -[spec] BEGdecay D. decay 0 0 1. K- pi+ goto D0dec 240D0ret: ENDdecay D. channel [spec] 1 0 goto 9999 241: if [charge]>0 then spec = D+ else spec = D- endif LUND cc -[spec] BEGdecay D+ decay 0 0 0.0860 K- pi+ pi+ decay 0 0 0.0128 K*~ pi+ decay 0 0 0.0230 K1430_0*~ pi+ goto D+dec 241D+ret: ENDdecay D+ BEGdecay K*. decay 0 0 0.3300 K+ pi- ENDdecay K*. BEGdecay K1430_0*. decay 0 0 0.3300 K+ pi- ENDdecay K1430_0*. channel [spec] 3 0 channel K*~ 1 0 channel K1430_0*~ 1 0 goto 9999 242: if [charge]>0 then spec = D_s+ else spec = D_s- endif LUND cc -[spec] BEGdecay D_s+ decay 0 0 0.0090 K+ K- pi+ decay 0 0 0.0360 phi pi+ decay 0 0 0.0330 K+ K*~ decay 0 0 0.0049 f_0 pi+ decay 0 0 0.0070 K+ K1430_0*~ ENDdecay D_s+ BEGdecay phi decay 0 0 0.4920 K+ K- ENDdecay phi BEGdecay K*. decay 0 0 0.3300 K+ pi- ENDdecay K*. BEGdecay K1430_0*. decay 0 0 0.3300 K+ pi- ENDdecay K1430_0*. BEGdecay f_0 decay 0 0 0.3300 K+ K- ENDdecay f_0 channel [spec] 5 0 channel phi 1 0 channel K*~ 1 0 channel K1430_0*~ 1 0 channel f_0 1 0 goto 9999 243: if [charge]>0 then spec = D*+ else spec = D*- endif LUND cc -[spec] BEGdecay D*+ decay 0 0 0.6770 D. pi+ ENDdecay D*+ BEGdecay D. decay 0 0 0.03650 K- pi+ goto D0dec 243D0ret: ENDdecay D. channel [spec] 1 0 channel D. 1 0 goto 9999 244: if [charge]>0 then spec = D+ else spec = D- endif LUND cc -[spec] BEGdecay [spec] decay 0 0 0.0139 K_S pi+ goto D+dec 244D+ret: ENDdecay D+ channel [spec] 1 0 goto 9999 * 300: LUND bb goto 9999 301: if [charge]>0 then spec = b-hadron else spec = b-hadron~ endif LUND bb BEGdecay b-hadron decay 0 42 0.10700 nu_e~ e- cquark specflav decay 0 42 0.10700 nu_mu~ mu- cquark specflav decay 0 42 0.04000 nu_tau~ tau- cquark specflav decay 0 42 0.50800 uquark~ dquark cquark specflav decay 0 42 0.08100 uquark~ cquark dquark specflav decay 0 42 0.14200 cquark~ squark cquark specflav decay 0 42 0.01000 cquark~ cquark squark specflav decay 0 42 0.00500 cquark~ squark uquark specflav goto b-hadrondec 301b-hadronret: ENDdecay b-hadron BEGdecay B. decay 0 85 1. goto B0dec 301B0ret: ENDdecay B. BEGdecay B+ decay 0 85 1. goto B+dec 301B+ret: ENDdecay B+ BEGdecay B_s. decay 0 85 1. goto B0sdec 301B0sret: ENDdecay B_s. BEGdecay B_c+ decay 0 85 1. goto Bc+dec 301Bc+ret: ENDdecay B_c+ BEGdecay Lam_b. decay 0 85 1. goto Lam_bdec 301Lam_bret: ENDdecay Lam_b. channel [spec] 8 1 channel B. 1 0 channel B+ 1 0 channel B_s. 1 0 channel B_c+ 1 0 channel Lam_b. 1 0 goto 9999 302: if [charge]>0 then spec = b-hadron else spec = b-hadron~ endif LUND bb BEGdecay b-hadron decay 0 42 0.14000 nu_e~ e- uquark specflav decay 0 42 0.14000 nu_mu~ mu- uquark specflav decay 0 42 0.03300 nu_tau~ tau- uquark specflav decay 0 42 0.66700 uquark~ dquark uquark specflav goto b-hadrondec 302b-hadronret: ENDdecay b-hadron BEGdecay B. decay 0 85 1. goto B0dec 302B0ret: ENDdecay B. BEGdecay B+ decay 0 85 1. goto B+dec 302B+ret: ENDdecay B+ BEGdecay B_s. decay 0 85 1. goto B0sdec 302B0sret: ENDdecay B_s. BEGdecay B_c+ decay 0 85 1. goto Bc+dec 302Bc+ret: ENDdecay B_c+ BEGdecay Lam_b. decay 0 85 1. goto Lam_bdec 302Lam_bret: ENDdecay Lam_b. channel [spec] 4 1 channel B. 1 0 channel B+ 1 0 channel B_s. 1 0 channel B_c+ 1 0 channel Lam_b. 1 0 goto 9999 303: particle mu+ mu- LUND bb goto 9999 304: particle e+ e- LUND bb goto 9999 305: lund bb -Jpsi BEGdecay Jpsi decay 0 0 1. mu- mu+ goto Jpsidec 305Jpsiret: ENDdecay Jpsi channel Jpsi 1 0 goto 9999 306: lund bb -Jpsi BEGdecay Jpsi decay 0 0 1. e- e+ goto Jpsidec 306Jpsiret: ENDdecay Jpsi channel Jpsi 1 0 goto 9999 307: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] goto 9999 308: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] goto 9999 309: if [charge]>0 then spec = Lam_b. else spec = Lam_b~. endif LUND bb -[spec] goto 9999 310: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] LEPTON 2 0. 0. 0. goto 9999 311: if [charge]>0 then spec = B_2*+ else spec = B_2*- endif LUND bb -[spec] BEGdecay B_2*+ decay 0 0 1. B. pi+ goto B_2*+dec 311B_2*+ret: ENDdecay B_2*+ BEGdecay B. decay 0 0 1. Jpsi K_S goto B0dec 311B0ret: ENDdecay B. BEGdecay Jpsi decay 0 0 1. mu- mu+ goto Jpsidec 311Jpsiret: ENDdecay Jpsi channel [spec] 1 1 channel B. 1 0 channel Jpsi 1 0 goto 9999 312: if [charge]>0 then spec = B_2*+ else spec = B_2*- endif LUND bb -[spec] BEGdecay B_2*+ decay 0 0 1. B. pi+ goto B_2*+dec 312B_2*+ret: ENDdecay B_2*+ BEGdecay B. decay 0 0 1. Jpsi K_S goto B0dec 312B0ret: ENDdecay B. BEGdecay Jpsi decay 0 0 1. e- e+ goto Jpsidec 312Jpsiret: ENDdecay Jpsi channel [spec] 1 1 channel B. 1 0 channel Jpsi 1 0 goto 9999 320: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. Jpsi K_S goto B0dec 320B0ret: ENDdecay B. BEGdecay Jpsi decay 0 0 1. mu- mu+ goto Jpsidec 320Jpsiret: ENDdecay Jpsi channel [spec] 1 1 channel Jpsi 1 0 goto 9999 321: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. Jpsi K_S goto B0dec 321B0ret: ENDdecay B. BEGdecay Jpsi decay 0 0 1. e- e+ goto Jpsidec 321Jpsiret: ENDdecay Jpsi channel [spec] 1 1 channel Jpsi 1 0 goto 9999 322: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. Jpsi K*. goto B0dec 322B0ret: ENDdecay B. BEGdecay Jpsi decay 0 0 1. mu- mu+ goto Jpsidec 322Jpsiret: ENDdecay Jpsi BEGdecay K*. decay 0 3 1. K_S pi. goto Kstar0dec 322Kstar0ret: ENDdecay K*. channel [spec] 1 1 channel Jpsi 1 0 channel K*. 1 0 goto 9999 323: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. Jpsi K*. goto B0dec 323B0ret: ENDdecay B. BEGdecay Jpsi decay 0 0 1. e- e+ goto Jpsidec 323Jpsiret: ENDdecay Jpsi BEGdecay K*. decay 0 3 1. K_S pi. goto Kstar0dec 323Kstar0ret: ENDdecay K*. channel [spec] 1 1 channel Jpsi 1 0 channel K*. 1 0 goto 9999 324: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. Jpsi K*. goto B0dec 324B0ret: ENDdecay B. BEGdecay Jpsi decay 0 0 1. mu- mu+ goto Jpsidec 324Jpsiret: ENDdecay Jpsi BEGdecay K*. decay 0 3 1. K+ pi- goto Kstar0dec 324Kstar0ret: ENDdecay K*. channel [spec] 1 1 channel Jpsi 1 0 channel K*. 1 0 goto 9999 325: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. Jpsi K*. goto B0dec 325B0ret: ENDdecay B. BEGdecay Jpsi decay 0 0 1. e- e+ goto Jpsidec 325Jpsiret: ENDdecay Jpsi BEGdecay K*. decay 0 3 1. K+ pi- goto Kstar0dec 325Kstar0ret: ENDdecay K*. channel [spec] 1 1 channel Jpsi 1 0 channel K*. 1 0 goto 9999 326: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. Jpsi rho. goto B0dec 326B0ret: ENDdecay B. BEGdecay Jpsi decay 0 0 1. mu- mu+ goto Jpsidec 326Jpsiret: ENDdecay Jpsi BEGdecay rho. decay 0 3 1. pi+ pi- goto rho0dec 326rho0ret: ENDdecay rho. channel [spec] 1 1 channel Jpsi 1 0 channel rho. 1 0 goto 9999 327: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. Jpsi rho. goto B0dec 327B0ret: ENDdecay B. BEGdecay Jpsi decay 0 0 1. e- e+ goto Jpsidec 327Jpsiret: ENDdecay Jpsi BEGdecay rho. decay 0 3 1. pi+ pi- goto rho0dec 327rho0ret: ENDdecay rho. channel [spec] 1 1 channel Jpsi 1 0 channel rho. 1 0 goto 9999 328: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. chi_c1 K_S goto B0dec 328B0ret: ENDdecay B. BEGdecay chi_c1 decay 0 0 1. Jpsi gamma goto chi_c1dec 328chi_c1ret: ENDdecay chi_c1 BEGdecay Jpsi decay 0 0 1. mu- mu+ goto Jpsidec 328Jpsiret: ENDdecay Jpsi channel [spec] 1 1 channel chi_c1 1 0 channel Jpsi 1 0 goto 9999 329: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. chi_c1 K_S goto B0dec 329B0ret: ENDdecay B. BEGdecay chi_c1 decay 0 0 1. Jpsi gamma goto chi_c1dec 329chi_c1ret: ENDdecay chi_c1 BEGdecay Jpsi decay 0 0 1. e- e+ goto Jpsidec 329Jpsiret: ENDdecay Jpsi channel [spec] 1 1 channel chi_c1 1 0 channel Jpsi 1 0 goto 9999 340: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. pi+ pi- goto B0dec 340B0ret: ENDdecay B. channel [spec] 1 1 goto 9999 341: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. K+ pi- goto B0dec 341B0ret: ENDdecay B. channel [spec] 1 1 goto 9999 342: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 0.5 a_1+ pi- decay 0 0 0.5 a_1- pi+ goto B0dec 342B0ret: ENDdecay B. BEGdecay a_1+ decay 0 0 1. rho. pi+ goto a1+dec 342a1+ret: ENDdecay a_1+ BEGdecay rho. decay 0 3 1. pi+ pi- goto rho0dec 342rho0ret: ENDdecay rho. channel [spec] 2 1 channel a_1+ 1 0 channel rho. 1 0 goto 9999 360: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] BEGdecay B_s. decay 0 0 1. D_s- pi+ goto B0sdec 360B0sret: ENDdecay B_s. channel [spec] 1 1 goto 9999 361: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] BEGdecay B_s. decay 0 0 1. D_s- K+ goto B0sdec 361B0sret: ENDdecay B_s. channel [spec] 1 1 goto 9999 362: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] BEGdecay B_s. decay 0 0 1. K- K+ goto B0sdec 362B0sret: ENDdecay B_s. channel [spec] 1 1 goto 9999 363: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] BEGdecay B_s. decay 0 0 1. Jpsi phi goto B0sdec 363B0sret: ENDdecay B_s. BEGdecay Jpsi decay 0 0 1. mu- mu+ goto Jpsidec 363Jpsiret: ENDdecay Jpsi BEGdecay phi decay 0 3 1. K+ K- goto phidec 363phiret: ENDdecay phi channel [spec] 1 1 channel Jpsi 1 0 channel phi 1 0 goto 9999 364: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] BEGdecay B_s. decay 0 0 1. Jpsi phi goto B0sdec 364B0sret: ENDdecay B_s. BEGdecay Jpsi decay 0 0 1. e- e+ goto Jpsidec 364Jpsiret: ENDdecay Jpsi BEGdecay phi decay 0 3 1. K+ K- goto phidec 364phiret: ENDdecay phi channel [spec] 1 1 channel Jpsi 1 0 channel phi 1 0 goto 9999 370: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 42 0.5 nu_e e+ D*- decay 0 42 0.5 nu_mu mu+ D*- decay 1 42 0.09524 nu_e e+ D- decay 1 42 0.26190 nu_e e+ D*- decay 1 42 0.02381 nu_e e+ D_1- decay 1 42 0.02381 nu_e e+ D_0*- decay 1 42 0.03810 nu_e e+ D_1*- decay 1 42 0.05714 nu_e e+ D_2*- decay 1 42 0.09524 nu_mu mu+ D- decay 1 42 0.26190 nu_mu mu+ D*- decay 1 42 0.02381 nu_mu mu+ D_1- decay 1 42 0.02381 nu_mu mu+ D_0*- decay 1 42 0.03810 nu_mu mu+ D_1*- decay 1 42 0.05714 nu_mu mu+ D_2*- ENDdecay B. BEGdecay B+ decay 1 42 0.09524 nu_e e+ D~ decay 1 42 0.26190 nu_e e+ D*~ decay 1 42 0.02381 nu_e e+ D_1~ decay 1 42 0.02381 nu_e e+ D_0*~ decay 1 42 0.03810 nu_e e+ D_1*~ decay 1 42 0.05714 nu_e e+ D_2*~ decay 1 42 0.09524 nu_mu mu+ D~ decay 1 42 0.26190 nu_mu mu+ D*~ decay 1 42 0.02381 nu_mu mu+ D_1~ decay 1 42 0.02381 nu_mu mu+ D_0*~ decay 1 42 0.03810 nu_mu mu+ D_1*~ decay 1 42 0.05714 nu_mu mu+ D_2*~ ENDdecay B+ BEGdecay B_s. decay 1 42 0.09524 nu_e e+ D_s- decay 1 42 0.26190 nu_e e+ D_s*- decay 1 42 0.02381 nu_e e+ D_s1- decay 1 42 0.02381 nu_e e+ D_s0*- decay 1 42 0.03810 nu_e e+ D_s1*- decay 1 42 0.05714 nu_e e+ D_s2*- decay 1 42 0.09524 nu_mu mu+ D_s- decay 1 42 0.26190 nu_mu mu+ D_s*- decay 1 42 0.02381 nu_mu mu+ D_s1- decay 1 42 0.02381 nu_mu mu+ D_s0*- decay 1 42 0.03810 nu_mu mu+ D_s1*- decay 1 42 0.05714 nu_mu mu+ D_s2*- ENDdecay B_s. BEGdecay B_c+ decay 1 42 0.05660 nu_e e+ eta_c decay 1 42 0.16510 nu_e e+ Jpsi decay 1 42 0.05660 nu_mu mu+ eta_c decay 1 42 0.16510 nu_mu mu+ Jpsi decay 1 42 0.06600 e+ nu_e B_s. decay 1 42 0.19810 e+ nu_e B_s*. decay 1 42 0.06600 mu+ nu_mu B_s. decay 1 42 0.19810 mu+ nu_mu B_s*. decay 1 42 0.00470 e+ nu_e B. decay 1 42 0.00940 e+ nu_e B*. decay 1 42 0.00470 mu+ nu_mu B. decay 1 42 0.00940 mu+ nu_mu B*. ENDdecay B_c+ BEGdecay Lam_b. decay 1 42 0.5 nu_e~ e- Lam_c+ decay 1 42 0.5 nu_mu~ mu- Lam_c+ ENDdecay Lam_b. BEGdecay b-hadron decay 1 42 0.5 nu_e~ e- cquark specflav decay 1 42 0.5 nu_mu~ mu- cquark specflav ENDdecay b-hadron BEGdecay D*+ decay 0 3 1. D. pi+ goto D*+dec 370D*+ret: ENDdecay D*+ channel [spec] 2 1 channel D*+ 1 0 goto 9999 371: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] LEPTON 2 0. 0. 0. BEGdecay B_s. decay 0 42 0.13333 nu_e e+ D_s- decay 0 42 0.36667 nu_e e+ D_s*- decay 0 42 0.13333 nu_mu mu+ D_s- decay 0 42 0.36667 nu_mu mu+ D_s*- goto B0sdec 371B0sret: ENDdecay B_s. channel [spec] 4 1 goto 9999 372: if [charge]>0 then spec = B. else spec = B.~ endif LUND bb -[spec] BEGdecay B. decay 0 42 1. nu_mu mu+ D- goto B0dec 372B0ret: ENDdecay B. BEGdecay D+ decay 0 0 1. K- pi+ pi+ goto D+dec 372D+ret: ENDdecay D+ channel [spec] 1 1 channel D+ 1 0 goto 9999 373: if [charge]>0 then spec = B+ else spec = B- endif LUND bb -[spec] BEGdecay B+ decay 0 42 1. nu_mu mu+ D~ goto B+dec 373B+ret: ENDdecay B+ BEGdecay D. decay 0 0 0.03650 K- pi+ decay 0 0 0.00624 K- pi+ rho. decay 0 0 0.01800 K- pi+ pi+ pi- goto D0dec 373D0ret: ENDdecay D. BEGdecay rho. decay 0 3 1. pi+ pi- goto rho0dec 373rho0ret: ENDdecay rho. channel [spec] 1 0 channel D. 3 0 channel rho. 1 0 goto 9999 374: if [charge]>0 then spec = B. else spec = B.~ endif LUND bb -[spec] BEGdecay B. decay 0 42 1. nu_mu mu+ D*- goto B0dec 374B0ret: ENDdecay B. BEGdecay D*+ decay 0 3 1. D. pi+ goto D*+dec 374D*+ret: ENDdecay D*+ BEGdecay D. decay 0 0 0.03650 K- pi+ decay 0 0 0.00624 K- pi+ rho. decay 0 0 0.01800 K- pi+ pi+ pi- goto D0dec 374D0ret: ENDdecay D. BEGdecay rho. decay 0 3 1. pi+ pi- goto rho0dec 374rho0ret: ENDdecay rho. channel [spec] 1 1 channel D*+ 1 0 channel D. 3 0 channel rho. 1 0 goto 9999 * 400: PSI REWEI LUND jpsi -Jpsi goto 9999 401: PSI REWEI PARTICLE Jpsi chi_c0 chi_c1 chi_c2 LUND jpsi goto 9999 410: PSI REWEI LUND jpsi -Jpsi BEGdecay Jpsi decay 0 0 1. mu- mu+ goto Jpsidec 410Jpsiret: ENDdecay Jpsi channel Jpsi 1 0 goto 9999 411: PSI REWEI LUND jpsi -Jpsi BEGdecay Jpsi decay 0 0 1. e- e+ goto Jpsidec 411Jpsiret: ENDdecay Jpsi channel Jpsi 1 0 goto 9999 412: PSI REWEI LUND jpsi -chi_c1 BEGdecay chi_c1 decay 0 0 1. Jpsi gamma goto chi_c1dec 412chi_c1ret: ENDdecay chi_c1 BEGdecay Jpsi decay 0 0 1. mu- mu+ goto Jpsidec 412Jpsiret: ENDdecay Jpsi channel chi_c1 1 0 channel Jpsi 1 0 goto 9999 413: PSI REWEI LUND jpsi -chi_c1 BEGdecay chi_c1 decay 0 0 1. Jpsi gamma goto chi_c1dec 413chi_c1ret: ENDdecay chi_c1 BEGdecay Jpsi decay 0 0 1. e- e+ goto Jpsidec 413Jpsiret: ENDdecay Jpsi channel chi_c1 1 0 channel Jpsi 1 0 goto 9999 420: LUND CSM goto 9999 421: LUND COM goto 9999 422: LUND COM -Jpsi BEGdecay Jpsi decay 0 0 1. mu- mu+ goto Jpsidec 422Jpsiret: ENDdecay Jpsi channel Jpsi 1 0 goto 9999 430: LUND PSI2S BEGdecay psi2S decay 1 0 1. mu- mu+ ENDdecay psi2S goto 9999 431: LUND PSI2S BEGdecay psi2S decay 1 0 1. Jpsi pi+ pi- ENDdecay psi2S BEGdecay Jpsi decay 1 0 1. mu- mu+ ENDdecay Jpsi goto 9999 432: LUND PSI2S BEGdecay psi2S decay 1 0 1. Jpsi eta ENDdecay psi2S BEGdecay Jpsi decay 1 0 1. mu- mu+ ENDdecay Jpsi BEGdecay eta decay 1 0 1. gamma gamma ENDdecay eta goto 9999 * 500: UPSILON REWEI LUND upsi -Upsilon goto 9999 501: UPSILON REWEI PARTICLE Upsilon chi_b0 chi_b1 chi_b2 LUND upsi goto 9999 510: UPSILON REWEI LUND upsi -Upsilon BEGdecay Upsilon decay 0 0 1. mu- mu+ goto Upsilondec 510Upsilonret: ENDdecay Upsilon channel Upsilon 1 0 goto 9999 511: UPSILON REWEI LUND upsi -Upsilon BEGdecay Upsilon decay 0 0 1. e- e+ goto Upsilondec 511Upsilonret: ENDdecay Upsilon channel Upsilon 1 0 goto 9999 512: UPSILON REWEI LUND upsi -chi_b1 BEGdecay chi_b1 decay 0 0 1. Upsilon gamma goto chi_b1dec 512chi_b1ret: ENDdecay chi_b1 BEGdecay Upsilon decay 0 0 1. mu- mu+ goto Upsilondec 512Upsilonret: ENDdecay Upsilon channel chi_b1 1 0 channel Upsilon 1 0 goto 9999 513: UPSILON REWEI LUND upsi -chi_b1 BEGdecay chi_b1 decay 0 0 1. Upsilon gamma goto chi_b1dec 513chi_b1ret: ENDdecay chi_b1 BEGdecay Upsilon decay 0 0 1. e- e+ goto Upsilondec 513Upsilonret: ENDdecay Upsilon channel chi_b1 1 0 channel Upsilon 1 0 goto 9999 * 600: if [charge]>0 then spec = D. else spec = D~ endif LUND cc -[spec] BEGdecay D. decay 0 0 1. mu+ mu- goto D0dec 600D0ret: ENDdecay D. channel [spec] 1 0 goto 9999 601: if [charge]>0 then spec = D. else spec = D~ endif LUND cc -[spec] BEGdecay D. decay 0 0 1. mu+ mu- gamma goto D0dec 601D0ret: ENDdecay D. channel [spec] 1 0 goto 9999 602: if [charge]>0 then spec = D. else spec = D~ endif LUND cc -[spec] BEGdecay D. decay 0 0 1. e+ e- gamma goto D0dec 602D0ret: ENDdecay D. channel [spec] 1 0 goto 9999 610: if [charge]>0 then spec = D. else spec = D~ endif LUND cc -[spec] BEGdecay D. decay 0 0 1. gamma gamma goto D0dec 610D0ret: ENDdecay D. channel [spec] 1 0 goto 9999 * * 700: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. mu+ mu- goto B0dec 700B0ret: ENDdecay B. channel [spec] 1 1 goto 9999 701: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] BEGdecay B_s. decay 0 0 1. mu+ mu- goto B0sdec 701B0sret: ENDdecay B_s. channel [spec] 1 1 goto 9999 702: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. mu+ mu- gamma goto B0dec 702B0ret: ENDdecay B. channel [spec] 1 1 goto 9999 703: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. e+ e- gamma goto B0dec 703B0ret: ENDdecay B. channel [spec] 1 1 goto 9999 704: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] BEGdecay B_s. decay 0 0 1. mu+ mu- gamma goto B0sdec 704B0sret: ENDdecay B_s. channel [spec] 1 1 goto 9999 705: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] BEGdecay B_s. decay 0 0 1. e+ e- gamma goto B0sdec 705B0sret: ENDdecay B_s. channel [spec] 1 1 goto 9999 710: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. gamma gamma goto B0dec 710B0ret: ENDdecay B. channel [spec] 1 1 goto 9999 711: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] BEGdecay B_s. decay 0 0 1. gamma gamma goto B0sdec 711B0sret: ENDdecay B_s. channel [spec] 1 1 goto 9999 720: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. K_S mu+ mu- goto B0dec 720B0ret: ENDdecay B. channel [spec] 1 1 goto 9999 721: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. K_S e+ e- goto B0dec 721B0ret: ENDdecay B. channel [spec] 1 1 goto 9999 722: if [charge]>0 then spec = B+ else spec = B- endif LUND bb -[spec] BEGdecay B+ decay 0 0 1. K+ mu+ mu- goto B+dec 722B+ret: ENDdecay B+ channel [spec] 1 0 goto 9999 723: if [charge]>0 then spec = B+ else spec = B- endif LUND bb -[spec] BEGdecay B+ decay 0 0 1. K+ e+ e- goto B+dec 723B+ret: ENDdecay B+ channel [spec] 1 0 goto 9999 724: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. K*. mu+ mu- goto B0dec 724B0ret: ENDdecay B. BEGdecay K*. decay 0 3 1. K+ pi- goto K0*dec 724K0*ret: ENDdecay K*. channel [spec] 1 1 channel K*. 1 0 goto 9999 725: if [charge]>0 then spec = B. else spec = B~ endif LUND bb -[spec] BEGdecay B. decay 0 0 1. K*. e+ e- goto B0dec 725B0ret: ENDdecay B. BEGdecay K*. decay 0 3 1. K+ pi- goto K0*dec 725K0*ret: ENDdecay K*. channel [spec] 1 1 channel K*. 1 0 goto 9999 726: if [charge]>0 then spec = B+ else spec = B- endif LUND bb -[spec] BEGdecay B+ decay 0 0 1. K*+ mu+ mu- goto B+dec 726B+ret: ENDdecay B+ BEGdecay K*+ goto K*+dec decay 0 3 1. K_S pi+ 726K*+ret: ENDdecay K*+ channel [spec] 1 0 channel K*+ 1 0 goto 9999 727: if [charge]>0 then spec = B+ else spec = B- endif LUND bb -[spec] BEGdecay B+ decay 0 0 1. K*+ e+ e- goto B+dec 727B+ret: ENDdecay B+ BEGdecay K*+ goto K*+dec decay 0 3 1. K_S pi+ 727K*+ret: ENDdecay K*+ channel [spec] 1 0 channel K*+ 1 0 goto 9999 730: if [charge]>0 then spec = B_s. else spec = B_s~ endif LUND bb -[spec] BEGdecay B_s. decay 0 0 1. phi gamma goto B0sdec 730B0sret: ENDdecay B_s. BEGdecay phi decay 0 3 1. K+ K- goto phidec 730phiret: ENDdecay phi channel [spec] 1 1 channel phi 1 0 goto 9999 740: if [charge]>0 then spec = Lam_b. else spec = Lam_b~. endif LUND bb -[spec] BEGdecay Lam_b. decay 0 0 1. Lam. gamma goto Lam_bdec 740Lam_bret: ENDdecay Lam_b. channel [spec] 1 0 goto 9999 741: if [charge]>0 then spec = Lam_b. else spec = Lam_b~. endif LUND bb -[spec] BEGdecay Lam_b. decay 0 0 1. Lam. mu+ mu- goto Lam_bdec 741Lam_bret: ENDdecay Lam_b. channel [spec] 1 0 goto 9999 742: if [charge]>0 then spec = Lam_b. else spec = Lam_b~. endif LUND bb -[spec] BEGdecay Lam_b. decay 0 0 1. Lam. e+ e- goto Lam_bdec 742Lam_bret: ENDdecay Lam_b. channel [spec] 1 0 goto 9999 * 801: if [charge]>0 then spec = Xi- else spec = Xi~+ endif LUND inel particle [spec] BEGdecay Xi- decay 1 0 1. Lam. pi- ENDdecay Xi- BEGdecay Lam. decay 1 0 1. p pi- ENDdecay Lam. channel [spec] 1 0 goto 9999 802: if [charge]>0 then spec = Xi*. else spec = Xi*~. endif LUND inel particle [spec] BEGdecay Xi*. decay 1 0 1. Xi- pi+ ENDdecay Xi*. BEGdecay Xi- decay 1 0 1. Lam. pi- ENDdecay Xi- BEGdecay Lam. decay 1 0 1. p pi- ENDdecay Lam. channel [spec] 1 0 goto 9999 803: if [charge]>0 then spec = Xi*- else spec = Xi*~+ endif LUND inel particle [spec] BEGdecay Xi*- decay 1 0 1. Xi- pi. ENDdecay Xi*- BEGdecay Xi- decay 1 0 1. Lam. pi- ENDdecay Xi- BEGdecay Lam. decay 1 0 1. p pi- ENDdecay Lam. channel [spec] 1 0 goto 9999 804: if [charge]>0 then spec = Ome- else spec = Ome~+ endif LUND inel particle [spec] BEGdecay Ome- decay 1 0 0.678 Lam. K- decay 1 0 0.236 Xi. pi- decay 1 0 0.086 Xi- pi. ENDdecay Ome- BEGdecay Xi- decay 1 0 1. Lam. pi- ENDdecay Xi- BEGdecay Lam. decay 1 0 1. p pi- ENDdecay Lam. BEGdecay Xi. decay 1 0 1. Lam. pi. ENDdecay Xi. channel [spec] 1 0 goto 9999 9999: * * *** Output Files (DST,HBOOK,PS) *** * *... IOFILE HB_MC_GEN DST * if [online]=on then IOFILE DST_OUT /tmp/.shm_MC -OM IOFILE MINI_OUT /tmp/.shm_MC -OM else * if [detector]=0 then out_name = [mc_outdir]/p[phys_proc]/d[geo_vers_1].[geo_vers_2]/w[wiremask]/gen_0/run[site]_[run_number]/run[site]_[run_number]_file[file_number] elseif [detector]=1 then out_name = [mc_outdir]/p[phys_proc]/d[geo_vers_1].[geo_vers_2]/w[wiremask]/gen/run[site]_[run_number]/run[site]_[run_number]_file[file_number] endif IOFILE DST_out [out_name].dst -O[geanformat] *... IOFILE HBOOK Out IO/HBOOK p[phys_proc]_w[wiremask]_run[site]_[run_number].hbook -O *... IOFILE PRISM Out IOFILE PRISM p[phys_proc]_d[geo_vers_1].[geo_vers_2]_w[wiremask]_run[site]_[run_number].ps -IG * * compression * IO/ZIP DST_OUT LZO IO/ZIP MINI_OUT LZO * * *** Detector Geometry for GEANT *** * if [detector]=0 then goto detend elseif [detector]=1 then GEO_VERS = [geo_vers_1].[geo_vers_2] VERS [GEO_VERS] if [online]=on then GEAN READ else GEAN endif * geant/control/swit 6 1 * HBDIGI/DGLEV OTR 3 setenv OtrHitPreparation 1 * else message 'Detector Geometry not yet available!!!' goto detend endif * *... GEO_file out. * if [online]<>'on' then if [run] = 0001 then IOFILE GEOMOUT [out_name].geo -O[geomformat] MMWRIT GEOMOUT NGHD CLOSE GEOMOUT endif endif * detend: * *... Target Configuration * if [itarg]=0 then TARGET 'C' if [detector]>0 then wire 'mat' 'C' endif elseif [itarg]=1 then TARGET 'Al' if [detector]>0 then wire 'mat' 'Al' endif elseif [itarg]=2 then TARGET 'Ti' if [detector]>0 then wire 'mat' 'Ti' endif elseif [itarg]=3 then TARGET 'W' if [detector]>0 then wire 'mat' 'W' endif elseif [itarg]=9 then if [detector]=0 then message 'Target material has to be specified on generator level!!!' message '************* Use Ti as default **************' TARGET 'Ti' endif else message 'Invalid target material, Ti used as default!!!' TARGET 'Ti' endif * if ([detector]>0) then * wire '1i' 1. wire '1o' 1. wire '1t' 1. wire '1b' 1. wire '2i' 1. wire '2o' 1. wire '2t' 1. wire '2b' 1. * IIb = mod([wiremask],10) wiremask = ([wiremask] - [IIb])/10 IIt = mod([wiremask],10) wiremask = ([wiremask] - [IIt])/10 IIo = mod([wiremask],10) wiremask = ([wiremask] - [IIo])/10 IIi = mod([wiremask],10) wiremask = ([wiremask] - [IIi])/10 Ib = mod([wiremask],10) wiremask = ([wiremask] - [Ib])/10 It = mod([wiremask],10) wiremask = ([wiremask] - [It])/10 Io = mod([wiremask],10) wiremask = ([wiremask] - [Io])/10 Ii = mod([wiremask],10) * if [Ii]=0 then wire '1i' 0. endif if [Io]=0 then wire '1o' 0. endif if [It]=0 then wire '1t' 0. endif if [Ib]=0 then wire '1b' 0. endif if [IIi]=0 then wire '2i' 0. endif if [IIo]=0 then wire '2o' 0. endif if [IIt]=0 then wire '2t' 0. endif if [IIb]=0 then wire '2b' 0. endif endif * * ... Change in the beam size (command available since ARTE-04-01-r5). * * The default beam width: pbeam pos 1i 0.0445 0. pbeam pos 1b 0.0445 0. pbeam pos 2i 0.0445 0. pbeam pos 2b 0.0445 0. * * Narrow beam component: * pbeam pos 1i 0.0200 0. * pbeam pos 1b 0.0160 0. * pbeam pos 2i 0.0200 0. * pbeam pos 2b 0.0160 0. * * *... Batch/Interactive Job * if [interactive]<>'on' then RUN [statistics] EXIT endif * goto done * *... Decay channels for the B0 * B0dec: decay 1 42 0.02000 nu_e e+ D- decay 1 42 0.05500 nu_e e+ D*- decay 1 42 0.00500 nu_e e+ D_1- decay 1 42 0.00500 nu_e e+ D_0*- decay 1 42 0.00800 nu_e e+ D_1*- decay 1 42 0.01200 nu_e e+ D_2*- decay 1 42 0.02000 nu_mu mu+ D- decay 1 42 0.05500 nu_mu mu+ D*- decay 1 42 0.00500 nu_mu mu+ D_1- decay 1 42 0.00500 nu_mu mu+ D_0*- decay 1 42 0.00800 nu_mu mu+ D_1*- decay 1 42 0.01200 nu_mu mu+ D_2*- decay 1 42 0.00700 nu_tau tau+ D- decay 1 42 0.02100 nu_tau tau+ D*- decay 1 0 0.00350 D- pi+ decay 1 0 0.01100 D- rho+ decay 1 0 0.00550 D- a_1+ decay 1 0 0.00420 D*- pi+ decay 1 0 0.00900 D*- rho+ decay 1 0 0.01800 D*- a_1+ decay 1 0 0.01500 D- D_s+ decay 1 0 0.01850 D- D_s*+ decay 1 0 0.01350 D*- D_s+ decay 1 0 0.02500 D*- D_s*+ decay 1 0 0.00040 eta_c K. decay 1 0 0.00070 eta_c K*. decay 1 0 0.00080 Jpsi K. decay 1 0 0.00140 Jpsi K*. decay 1 0 0.00190 chi_c1 K. decay 1 0 0.00250 chi_c1 K*. decay 1 0 0.00001 pi+ pi- decay 1 0 0.00003 rho+ pi- decay 1 0 0.00001 pi+ rho- decay 1 0 0.00002 rho+ rho- decay 1 32 0.02000 cd_0~ ud_0 decay 1 32 0.04500 cd_1~ ud_1 decay 1 48 0.38013 uquark dquark~ cquark~ dquark decay 1 48 0.02100 uquark squark~ cquark~ dquark decay 1 48 0.01800 uquark cquark~ dquark~ dquark decay 1 48 0.00090 uquark cquark~ squark~ dquark decay 1 13 0.13000 cquark squark~ cquark~ dquark decay 1 13 0.00600 cquark dquark~ cquark~ dquark decay 1 48 0.00600 uquark dquark~ uquark~ dquark decay 1 48 0.00400 cquark squark~ uquark~ dquark goto [iphys]B0ret * *... Decay channels for the J/Psi * Jpsidec: decay 1 0 0.05930 e- e+ decay 1 0 0.05880 mu- mu+ decay 1 12 0.88190 rndmflav rndmflav~ goto [iphys]Jpsiret * *... Decay channels for the chi_c1 * chi_c1dec: decay 1 0 0.31600 Jpsi gamma decay 1 12 0.68400 rndmflav rndmflav~ goto [iphys]chi_c1ret * *... Decay channels for the Upsilon * Upsilondec: decay 1 0 0.02380 e- e+ decay 1 0 0.02480 mu- mu+ decay 1 0 0.02670 tau- tau+ decay 1 32 0.01500 dquark dquark~ decay 1 32 0.04500 uquark uquark~ decay 1 32 0.01500 squark squark~ decay 1 32 0.04500 cquark cquark~ decay 1 4 0.77570 gluon gluon gluon decay 1 4 0.02900 gamma gluon gluon goto [iphys]Upsilonret * *... Decay channels for the chi_b1 * chi_b1dec: decay 1 0 0.35000 Upsilon gamma decay 1 32 0.65000 gluon gluon goto [iphys]chi_b1ret * *... Decay channels for the K*0 * Kstar0dec: decay 1 3 0.66500 K+ pi- decay 1 3 0.33300 K. pi. decay 1 0 0.00200 K. gamma goto [iphys]Kstar0ret * *... Decay channels for the B0s * B0sdec: decay 1 42 0.02000 nu_e e+ D_s- decay 1 42 0.05500 nu_e e+ D_s*- decay 1 42 0.00500 nu_e e+ D_s1- decay 1 42 0.00500 nu_e e+ D_s0*- decay 1 42 0.00800 nu_e e+ D_s1*- decay 1 42 0.01200 nu_e e+ D_s2*- decay 1 42 0.02000 nu_mu mu+ D_s- decay 1 42 0.05500 nu_mu mu+ D_s*- decay 1 42 0.00500 nu_mu mu+ D_s1- decay 1 42 0.00500 nu_mu mu+ D_s0*- decay 1 42 0.00800 nu_mu mu+ D_s1*- decay 1 42 0.01200 nu_mu mu+ D_s2*- decay 1 42 0.00800 nu_tau tau+ D_s- decay 1 42 0.02400 nu_tau tau+ D_s*- decay 1 0 0.00350 D_s- pi+ decay 1 0 0.01100 D_s- rho+ decay 1 0 0.00550 D_s- a_1+ decay 1 0 0.00420 D_s*- pi+ decay 1 0 0.00900 D_s*- rho+ decay 1 0 0.01800 D_s*- a_1+ decay 1 0 0.01500 D_s- D_s+ decay 1 0 0.01850 D_s- D_s*+ decay 1 0 0.01350 D_s*- D_s+ decay 1 0 0.02500 D_s*- D_s*+ decay 1 0 0.00020 eta_c eta decay 1 0 0.00020 eta_c eta` decay 1 0 0.00070 eta_c phi decay 1 0 0.00040 Jpsi eta decay 1 0 0.00040 Jpsi eta` decay 1 0 0.00140 Jpsi phi decay 1 0 0.00100 chi_c1 eta decay 1 0 0.00090 chi_c1 eta` decay 1 0 0.00250 chi_c1 phi decay 1 0 0.00001 pi+ K- decay 1 0 0.00003 rho+ K- decay 1 0 0.00001 pi+ K*- decay 1 0 0.00002 rho+ K*- decay 1 48 0.43702 uquark dquark~ cquark~ squark decay 1 13 0.08000 uquark cquark~ dquark~ squark decay 1 13 0.07000 cquark squark~ cquark~ squark decay 1 13 0.02000 cquark cquark~ squark~ squark decay 1 42 0.01500 uquark dquark~ uquark~ squark decay 1 42 0.00500 cquark squark~ uquark~ squark goto [iphys]B0sret * *... Decay channels for the phi * phidec: decay 1 3 0.49200 K+ K- decay 1 3 0.33700 K_L K_S decay 1 0 0.04400 rho- pi+ decay 1 0 0.04400 rho. pi. decay 1 0 0.04400 rho+ pi- decay 1 1 0.02300 pi+ pi- pi. decay 1 0 0.01300 gamma eta decay 1 0 0.00200 gamma omega decay 1 0 0.00100 pi. gamma goto [iphys]phiret * *... Decay channels for the B+ * B+dec: decay 1 42 0.02000 nu_e e+ D~ decay 1 42 0.05500 nu_e e+ D*~ decay 1 42 0.00500 nu_e e+ D_1~ decay 1 42 0.00500 nu_e e+ D_0*~ decay 1 42 0.00800 nu_e e+ D_1*~ decay 1 42 0.01200 nu_e e+ D_2*~ decay 1 42 0.02000 nu_mu mu+ D~ decay 1 42 0.05500 nu_mu mu+ D*~ decay 1 42 0.00500 nu_mu mu+ D_1~ decay 1 42 0.00500 nu_mu mu+ D_0*~ decay 1 42 0.00800 nu_mu mu+ D_1*~ decay 1 42 0.01200 nu_mu mu+ D_2*~ decay 1 42 0.00700 nu_tau tau+ D~ decay 1 42 0.02100 nu_tau tau+ D*~ decay 1 0 0.00350 D~ pi+ decay 1 0 0.01100 D~ rho+ decay 1 0 0.00550 D~ a_1+ decay 1 0 0.00420 D*~ pi+ decay 1 0 0.00900 D*~ rho+ decay 1 0 0.01800 D*~ a_1+ decay 1 0 0.01500 D~ D_s+ decay 1 0 0.01850 D~ D_s*+ decay 1 0 0.01350 D*~ D_s+ decay 1 0 0.02500 D*~ D_s*+ decay 1 0 0.00040 eta_c K+ decay 1 0 0.00070 eta_c K*+ decay 1 0 0.00080 Jpsi K+ decay 1 0 0.00140 Jpsi K*+ decay 1 0 0.00190 chi_c1 K+ decay 1 0 0.00250 chi_c1 K*+ decay 1 0 0.00001 pi+ pi. decay 1 0 0.00003 rho+ pi. decay 1 0 0.00001 pi+ rho. decay 1 0 0.00002 rho+ rho. decay 1 32 0.06000 cd_1~ uu_1 decay 1 48 0.44353 uquark dquark~ cquark~ uquark decay 1 48 0.02150 uquark squark~ cquark~ uquark decay 1 13 0.09000 cquark squark~ cquark~ uquark decay 1 13 0.00400 cquark dquark~ cquark~ uquark decay 1 48 0.00600 uquark dquark~ uquark~ uquark decay 1 48 0.00400 cquark squark~ uquark~ uquark goto [iphys]B+ret * *... Decay channels for the B_c+ * Bc+dec: decay 1 0 0.04700 nu_tau tau+ decay 1 11 0.12200 cquark squark~ decay 1 11 0.00600 cquark dquark~ decay 1 42 0.01200 nu_e e+ eta_c decay 1 42 0.03500 nu_e e+ Jpsi decay 1 42 0.01200 nu_mu mu+ eta_c decay 1 42 0.03500 nu_mu mu+ Jpsi decay 1 42 0.00300 nu_tau tau+ eta_c decay 1 42 0.00700 nu_tau tau+ Jpsi decay 1 42 0.15000 uquark dquark~ cquark~ cquark decay 1 42 0.03700 uquark cquark~ dquark~ cquark decay 1 42 0.00800 uquark squark~ cquark~ cquark decay 1 42 0.00200 uquark cquark~ squark~ cquark decay 1 42 0.05000 cquark squark~ cquark~ cquark decay 1 42 0.01500 cquark cquark~ squark~ cquark decay 1 42 0.00300 cquark dquark~ cquark~ cquark decay 1 42 0.00100 cquark cquark~ dquark~ cquark decay 1 42 0.01400 e+ nu_e B_s. decay 1 42 0.04200 e+ nu_e B_s*. decay 1 42 0.01400 mu+ nu_mu B_s. decay 1 42 0.04200 mu+ nu_mu B_s*. decay 1 42 0.24000 dquark~ uquark squark bquark~ decay 1 42 0.06500 dquark~ squark uquark bquark~ decay 1 42 0.01200 squark~ uquark squark bquark~ decay 1 42 0.00300 squark~ squark uquark bquark~ decay 1 42 0.00100 e+ nu_e B. decay 1 42 0.00200 e+ nu_e B*. decay 1 42 0.00100 mu+ nu_mu B. decay 1 42 0.00200 mu+ nu_mu B*. decay 1 42 0.01400 dquark~ uquark dquark bquark~ decay 1 42 0.00300 dquark~ dquark uquark bquark~ goto [iphys]Bc+ret * *... Decay channels for the Lambda_b * Lam_bdec: decay 1 42 0.10500 nu_e~ e- Lam_c+ decay 1 42 0.10500 nu_mu~ mu- Lam_c+ decay 1 42 0.04000 nu_tau~ tau- Lam_c+ decay 1 0 0.00770 Lam_c+ pi- decay 1 0 0.02000 Lam_c+ rho- decay 1 0 0.02350 Lam_c+ a_1- decay 1 0 0.02850 Lam_c+ D_s- decay 1 0 0.04350 Lam_c+ D_s*- decay 1 0 0.00110 eta_c Lam. decay 1 0 0.00220 Jpsi Lam. decay 1 0 0.00440 chi_c1 Lam. decay 1 48 0.42910 uquark~ dquark cquark ud_0 decay 1 13 0.08000 uquark~ cquark dquark ud_0 decay 1 13 0.07000 cquark~ squark cquark ud_0 decay 1 13 0.02000 cquark~ cquark squark ud_0 decay 1 42 0.01500 uquark~ dquark uquark ud_0 decay 1 42 0.00500 cquark~ squark uquark ud_0 goto [iphys]Lam_bret * *... Decay channels for the B_2*+ * B_2*+dec: decay 1 0 0.30000 B. pi+ decay 1 0 0.15000 B+ pi. decay 1 0 0.16000 B*. pi+ decay 1 0 0.08000 B*+ pi. decay 1 0 0.13000 B*. pi+ pi. decay 1 0 0.06000 B*+ pi+ pi- decay 1 0 0.08000 B. pi+ pi. decay 1 0 0.04000 B+ pi+ pi- goto [iphys]B_2*+ret * *... Decay channels for the generic b hadron * b-hadrondec: decay 1 42 0.10500 nu_e~ e- cquark specflav decay 1 42 0.10500 nu_mu~ mu- cquark specflav decay 1 42 0.04000 nu_tau~ tau- cquark specflav decay 1 42 0.50000 uquark~ dquark cquark specflav decay 1 42 0.08000 uquark~ cquark dquark specflav decay 1 42 0.14000 cquark~ squark cquark specflav decay 1 42 0.01000 cquark~ cquark squark specflav decay 1 42 0.01500 uquark~ dquark uquark specflav decay 1 42 0.00500 cquark~ squark uquark specflav goto [iphys]b-hadronret * *... Decay channels for the D0 * D0dec: decay 1 42 0.03400 e+ nu_e K- decay 1 42 0.02700 e+ nu_e K*- decay 1 42 0.00200 e+ nu_e K~ pi- decay 1 42 0.00200 e+ nu_e K- pi. decay 1 42 0.00400 e+ nu_e K*~ pi- decay 1 42 0.00400 e+ nu_e K*- pi. decay 1 42 0.00200 e+ nu_e pi- decay 1 42 0.00200 e+ nu_e rho- decay 1 42 0.03400 mu+ nu_mu K- decay 1 42 0.02700 mu+ nu_mu K*- decay 1 42 0.00200 mu+ nu_mu K~ pi- decay 1 42 0.00200 mu+ nu_mu K- pi. decay 1 42 0.00400 mu+ nu_mu K*~ pi- decay 1 42 0.00400 mu+ nu_mu K*- pi. decay 1 42 0.00200 mu+ nu_mu pi- decay 1 42 0.00200 mu+ nu_mu rho- decay 1 0 0.03650 K- pi+ decay 1 0 0.04500 K*- pi+ decay 1 0 0.07300 K- rho+ decay 1 0 0.06200 K*- rho+ decay 1 0 0.02100 K~ pi. decay 1 0 0.02100 K*~ pi. decay 1 0 0.02100 K*~ eta decay 1 0 0.00610 K~ rho. decay 1 0 0.01500 K*~ rho. decay 1 0 0.02500 K~ omega decay 1 0 0.00880 K~ phi decay 1 0 0.07400 K- a_1+ decay 1 0 0.01090 K1270_1- pi+ decay 1 0 0.00410 K- K+ decay 1 0 0.00200 K*- K+ decay 1 0 0.00350 K- K*+ decay 1 0 0.00110 K~ K. decay 1 0 0.00100 K*~ K. decay 1 0 0.00270 K*~ K*. decay 1 0 0.00160 pi+ pi- decay 1 0 0.00160 pi. pi. decay 1 0 0.00180 phi rho. decay 1 0 0.01100 K- pi+ pi. decay 1 0 0.00630 K- pi+ rho. decay 1 0 0.00520 K- K+ K~ decay 1 0 0.01800 K~ pi+ pi- decay 1 0 0.01600 K*~ pi+ pi- decay 1 0 0.00340 K- K. pi+ decay 1 0 0.00360 K*~ K+ pi- decay 1 0 0.00090 K_S K_S K_S decay 1 0 0.00060 phi pi+ pi- decay 1 0 0.01500 pi+ pi- pi. decay 1 0 0.09230 K- pi+ pi. pi. decay 1 0 0.01800 K- pi+ pi+ pi- decay 1 0 0.02200 K~ pi+ pi- pi. decay 1 0 0.00770 K*~ pi+ pi- pi. decay 1 0 0.00900 K~ K+ K- pi. decay 1 0 0.00750 pi+ pi+ pi- pi- decay 1 0 0.02400 K- pi+ pi+ pi- pi. decay 1 0 0.00850 K~ pi+ pi+ pi- pi- decay 1 0 0.06700 K~ pi+ pi- pi. pi. decay 1 0 0.05110 K~ rho. pi. pi. pi. decay 1 0 0.01700 pi+ pi+ pi- pi- pi. decay 1 0 0.00040 rho. pi+ pi+ pi- pi- decay 1 0 0.00280 K+ K- pi+ pi- pi. goto [iphys]D0ret * * Decay channels for the D+ * D+dec: decay 1 42 0.07000 e+ nu_e K~ decay 1 42 0.06500 e+ nu_e K*~ decay 1 42 0.00500 e+ nu_e K~ pi. decay 1 42 0.00500 e+ nu_e K- pi+ decay 1 42 0.01100 e+ nu_e K*~ pi. decay 1 42 0.01100 e+ nu_e K*- pi+ decay 1 42 0.00100 e+ nu_e pi. decay 1 42 0.00100 e+ nu_e eta decay 1 42 0.00100 e+ nu_e eta` decay 1 42 0.00100 e+ nu_e rho. decay 1 42 0.00100 e+ nu_e omega decay 1 42 0.07000 mu+ nu_mu K~ decay 1 42 0.06500 mu+ nu_mu K*~ decay 1 42 0.00500 mu+ nu_mu K~ pi. decay 1 42 0.00500 mu+ nu_mu K- pi+ decay 1 42 0.01100 mu+ nu_mu K*~ pi. decay 1 42 0.01100 mu+ nu_mu K*- pi+ decay 1 42 0.00100 mu+ nu_mu pi. decay 1 42 0.00100 mu+ nu_mu eta decay 1 42 0.00100 mu+ nu_mu eta` decay 1 42 0.00100 mu+ nu_mu rho. decay 1 42 0.00100 mu+ nu_mu omega decay 1 0 0.02600 K~ pi+ decay 1 0 0.01900 K*~ pi+ decay 1 0 0.06600 K~ rho+ decay 1 0 0.04100 K*~ rho+ decay 1 0 0.04500 K1400_1~ pi+ decay 1 0 0.07600 K~ a_1+ decay 1 0 0.00730 K~ K+ decay 1 0 0.00470 K*~ K+ decay 1 0 0.00470 K~ K*+ decay 1 0 0.02600 K*~ K*+ decay 1 0 0.00100 pi. pi+ decay 1 0 0.00060 pi. rho+ decay 1 0 0.00660 eta pi+ decay 1 0 0.00500 eta rho+ decay 1 0 0.00300 eta` pi+ decay 1 0 0.00300 eta` rho+ decay 1 0 0.00060 rho. pi+ decay 1 0 0.00060 rho. rho+ decay 1 0 0.00100 omega pi+ decay 1 0 0.00100 omega rho+ decay 1 0 0.00600 phi pi+ decay 1 0 0.00500 phi rho+ decay 1 0 0.01200 K~ pi+ pi. decay 1 0 0.00570 K*~ pi+ rho. decay 1 0 0.06700 K- pi+ pi+ decay 1 0 0.00800 K- rho+ pi+ decay 1 0 0.00220 pi+ pi+ pi- decay 1 0 0.02700 K~ K+ K~ decay 1 0 0.00400 K- K+ pi+ decay 1 0 0.01900 phi pi+ pi. decay 1 0 0.01200 K~ pi+ pi+ pi- decay 1 0 0.00200 K*~ pi+ pi+ pi- decay 1 0 0.00900 K- pi+ pi+ pi. decay 1 0 0.02180 pi+ pi+ pi- pi. decay 1 0 0.00100 K- pi+ pi+ pi+ pi- decay 1 0 0.02200 K- pi+ pi+ pi. pi. decay 1 0 0.08700 K~ pi+ pi+ pi- pi. decay 1 0 0.00100 K~ rho. pi+ pi+ pi- decay 1 0 0.00190 K- rho. pi+ pi+ pi. decay 1 0 0.00150 pi+ pi+ pi+ pi- pi- decay 1 0 0.00280 rho. pi+ pi+ pi- pi. goto [iphys]D+ret * *... Decay channels for the D*+ * D*+dec: decay 1 3 0.67700 D. pi+ decay 1 3 0.30700 D+ pi. decay 1 0 0.01600 D+ gamma goto [iphys]D*+ret * *... Decay channels for the K0* * K0*dec: decay 1 3 0.66500 K+ pi- decay 1 3 0.33300 K. pi. decay 1 0 0.00200 K. gamma goto [iphys]K0*ret * *... Decay channels for the K*+ * K*+dec: decay 1 3 0.66600 K. pi+ decay 1 3 0.33300 K+ pi. decay 1 0 0.00100 K+ gamma goto [iphys]K*+ret * *... Decay channels for the rho0 * rho0dec: decay 1 3 0.98908 pi+ pi- decay 1 0 0.00990 pi+ pi- gamma decay 1 0 0.00062 pi. gamma decay 1 0 0.00031 eta gamma decay 1 0 0.00005 mu- mu+ decay 1 0 0.00004 e- e+ goto [iphys]rho0ret * *... Decay channels for the a1+ * a1+dec: decay 1 0 0.50000 rho. pi+ decay 1 0 0.50000 rho+ pi. goto [iphys]a1+ret * done: ENDKUMAC