(Do not edit this file, it was generated from beget.tex.) This section describes the data cards used to control BEGET. The data cards are read in with the FFREAD package from the CERN library. The defaults are set in subroutine BEGFFKEY and BEGINIT. Data Card Function --------- ---------------------------------------------------------------- ============Generator Controls============================================== GEN Choose what to generate CP Select CP violation and/or \bo-\bobar\ mixing BACK Activate beam background from TURTLE files BACF Set the location of a TURTLE file list COSM Control cosmic ray generator ============Decay Controls================================================== B0B0 Selectively deactivate \bo\ and/or \bobar\ decay modes JPSI Control \J\ decay JETD Specify a file to read in JETSET decay information ============Beam Controls=================================================== BEAM Set the beam energies SBEAM Set the beam energy spread VERT Set the primary vertex location SVERT Set the primary vertex spread IVERT Set the primary vertex smearing method ============General Controls================================================ GENP Print event information on the terminal JETO Print the decay table for the listed particles RAN Control the random number seeds EVNT Set the number of events to generate REJE Control user event rejection Default Data Card Value Default Function --------- ----- ----------------------------------------------------- ======================General Controls===================================== GEN 2 Generate \upsbb, decay using JETSET CP 0 No CP violation or \bo-\bobar\ mixing BACK 0 do not add TURTLE background BACF '' no default TURTLE file list COSM (see below) default BaBar geometry used if cosmics are generated ======================Decay Controls======================================= B0B0 1 1 Turn on all \bo\ decays JPSI -1 Turn on all \J\ decays, unless B0B0 is not 1 1 (see below) JETD '' Do not read any external decay file ======================Beam Controls======================================= BEAM 9.0 -3.109 $p_z(e^-)$ = 9.0 GeV/c, $p_z(e^+)$ = -3.109 GeV/c SBEAM 0. 0. zero energy spread on beam (no smearing) VERT 0. 0. 0. event origin is (0,0,0) SVERT 0. 0. 0. no vertex smearing IVERT 0 0 0 use gaussians (sigmas set by SVERT) ======================General Controls==================================== GENP 0 0 Do not print the events on the terminal. JETO 0*20 no print out of decay tables RAN 1 0 use default seed #1 = 9876 54321 EVNT 1 set Nbeg_events = 1 REJE 0 accept all events **************************************************************************** GEN IMODE KF x1 x2 x3 x4 x5 x6 Choose what to generate **************************************************************************** GEN 2 0 0. 0. 0. 0. 0. 0. Defaults OVERVIEW +--------------------------------------------------------------+ |-100 NULL Physics Generator | | -99 call BEGUSER | +--------------------------------------------------------------+ | Single Particle (Pt) | | -13 kf Ptmin Ptmax cos(themin) cos(themax) phimin phimax | | -12 kf Ptmin Ptmax themin themax phimin phimax | | -11 kf Pt cos(theta) phi | | -10 kf Pt theta phi | +--------------------------------------------------------------+ | Single Particle (P) | | -3 kf Pmin Pmax cos(themin) cos(themax) phimin phimax | | -2 kf Pmin Pmax themin themax phimin phimax | | -1 kf P cos(theta) phi | | 0 kf P theta phi | +--------------------------------------------------------------+ | 1 Iopt Upsilon(4S) | | 2 Iopt Upsilon(4S) -> B0 B0~ only | | 3 Iopt Upsilon(4S) -> B+ B- only | +--------------------------------------------------------------+ | 4 kf e+e- -> qq~ (JETSET) | | 5 kf e+e- -> gamma -> "onium" (JETSET call luonia) | +--------------------------------------------------------------+ | 10 tau+ tau- KORALB | | 11 beam-gas EPC | | 12 cosmic rays HemiCosm | +--------------------------------------------------------------+ DETAILS ------- IMODE=-100 Do not call any physics generators. Useful if you want to look at TURTLE background by itself. IMODE= -99 Call subroutine BEGUSER to generate the event. IMODE=-10 to -13 same as IMODE 0 to -3 except: Pt (momentum transverse to beam pipe) instead of P momentum in GeV/c IMODE=-3 same as IMODE=-2 except: x3 x4 = cos(The_min) cos(The_max) Angles in degrees, momentum in GeV/c IMODE=-2 same as IMODE=0 except: x1 x2 = Pmin Pmax x3 x4 = The_min The_max x5 x6 = Phi_min Phi_max A random value for P, theta and phi is chosen for each event in the ranges specified by x1-x6 Angles in degrees, momentum in GeV/c IMODE=-1 same as IMODE=0 except: x1 x2 x3 = P cos(theta) Phi Angles in degrees, Momentum in GeV/c IMODE=0 User defined particle and four-vector. The particle is allowed to decay using JETSET. KF = JETSET code to identify particle x1 x2 x3 = P Theta Phi x4 x5 x6 = ignored P = momentum of particle in lab frame (GeV/c) Theta, Phi = direction of mom. vector in lab frame (degrees) IMODE=1 Generate Upsilon(4S) at rest in cm frame. Decay using JETSET and boost to lab frame. Angular distribution of decay: Iopt = 0 (default) flat in cos(theta) Iopt = 1 sin^2(theta) IMODE=2 Same as IMODE=1 with Upsilon(4S) -> B0 B0~ decay only. IMODE=3 Same as IMODE=1 with Upsilon(4S) -> B+ B- decay only. IMODE=4 Generate e+e- -> q q~ with JETSET, decay with JETSET. Get Ecm and boost from BEAM parameters. KF=0 Generate all quark pairs up to heaviest allowed by Ecm KF=n Generate only quark pairs defined by n n=1,2,3,4,5 = d,u,s,c,b IMODE=5 Generate e+e- -> gamma -> "onium" -> ggg or gg gamma -> parton shower with JETSET, decay with JETSET. Get Ecm and boost from BEAM parameters. KF=0 Generate g g g events only KF=n Generate mixture of ggg and gg gamma determined by Q^2 of n n=1,2,3,4,5 = d,u,s,c,b IMODE=10 Generate tau+ tau- using KORALB IMODE=11 Generate Beam-Gas interactions EPC library calculates the cross-sections IMODE=12 cosmic-ray generator, see data card COSM **************************************************************************** CP Ichoice Par1 Par2 Par3 Par4 IPar Select CP violation/B0 mixing **************************************************************************** CP 0 Defaults Ichoice=0 no CP violation or B0 mixing Ichoice=1 use BBKING version of CP/mixing routine ASMXCP_BBKING Ichoice=2 use ASLUND version of CP/mixing routine ASMXCP_ASLUND Defaults for Ichoice=1 and 2 are different CP 1 xd xs sincp coscp icp .7 8. .5 0. 1 Defaults CP 2 sin_2alpha sin_2beta sin_2gamma coscp .5 .8 .5 0. Defaults **************************************************************************** BACK imode ix1 ix2 rmean Control beam background from TURTLE files **************************************************************************** BACK 0 0 0 0. Defaults Add off-momentum beam particles and bremsstrahlung photons to the current event. The current implementation samples particles from a set of TURTLE files previously prepared. If no GEN data card is present then no physics generator will be called, i.e. the event will contain only background. imode .ne.0 ==> turn on background = 1 end job at end of background input file = 2 rewind file if EOF reached and reuse same background events ix1,ix2 range of beam crossings over which to generate backgrounds. (0 corresponds to trigger beam crossing.) rmean = 0.0 then use probability per beam crossing data contained in above file to calculate the number of events. > 0.0 then input a fixed rmean number of events spread over the beam crossing range ix1,ix2. Use prob/crossing data to determine RELATIVE sampling frequency from files. < 0.0 then use a Poisson with mean -rmean to determine the number of events to spread over over the selected beam crossing range. Use prob/crossing data to determine RELATIVE sampling frequency from files. **************************************************************************** BACf 'filename' Location of TURTLE file list **************************************************************************** BACF '' Defaults filename = Full pathname to a file containing a list of TURTLE files and some information about out the files were generated. For example file $BFROOT/data/ap75.list which contains: 2 1.47000e-02 25 /nfs/juno/u5/data/heb.dat 4.34200e-02 20 /nfs/juno/u5/data/leb.dat The first number in the file is the number of files of TURTLE rays in the list. In this sample there are two files of turtle rays (one for high energy beam (heb) and one for low (leb)). The first column is the probability per crossing for the events in this file and the second column is the number of events contained in the file. The format of the ray files differ in some important respects from the previous OBJEGS input files and should not be considered compatible. **************************************************************************** Control cosmic ray generator COSM FFuse Pmode Pmin Pmax Geant DetLen DetRad DetOff **************************************************************************** COSM TRUE 3 1.0 1000.0 TRUE 870.0 420.0 0.0 Defaults FFuse Use FFREAD card (not local defaults in BegHemiCosm) Pmode Momentum spectrum type Pmode = 0 Constant = Pmin Pmode = 1 1/p^2 Pmode = 2 From measured values Pmode = 3 DarInv from theory, with cos(theta) dependence Pmin Minimum p (GeV/c) Pmax Maximum p (GeV/c) Geant Use GEANT detector geometry (if present) instead of the following DetLen Detector cylinder length (cm) DetRad Detector cylinder radius (cm) DetOff Detector cylinder offset (cm) example: 1/p^2 spectrum: COSM TRUE 1 p spectrum from files and use a large detector: COSM TRUE 2 1.0 1000.0 FALSE 12000.0 500.0 0.0 # See the file src/hemicosm.F for more details, or run one-event jobs and read the printouts. **************************************************************************** B0B0 IB0 IB0~ Selectively deactivate B0 and/or B0~ decay modes **************************************************************************** B0B0 1 1 Defaults IB0 = controls B0 decay IB0~= controls B0~ decay processed by subroutine BDECAY example: generates: GEN 2 Upsilon(4S) -> B0 B0~ B0B0 11 3 | +---> mu- nu_mu~ X +------> Psi K0L +------> e+e- or mu+mu- For decay modes 11,12, and 13 the Psi decay modes are switched to e+e- and mu+mu- only. The user can override this with the JPSI data card (see below) IB0=0 turn off all decays 1 allow all decays 2 only B0 -> e nu_e X 3 -> mu nu_mu X 4 -> tau nu_tau X 5 unused 6 unused 7 unused 8 only B0 -> e X or mu X or tau X 9 -> e X or mu X 10 unused 11 only B0 -> Psi K0L \ 12 -> Psi K0S | Psi -> e+e- or mu+mu- 13 -> Psi K*0 / (see JPSI data card) 14 -> D+ D- 15 -> D*+ D*- 16 -> pi+ pi- 17 -> rho+/- pi-/+ 18 -> a1+/- pi-/+ **************************************************************************** JPSI n Control J/Psi decay modes **************************************************************************** JPSI -1 Defaults n=-1 all decays allowed unless B0B0 selected decay 11,12,13 then only e+e- and mu+mu- allowed 0 turn off all decays 1 all decays allowed 2 Psi -> e+e- only 3 mu+mu- only 9 e+e- and mu+mu- only **************************************************************************** JETD 'filename' Specify a file to read in JETSET decay information **************************************************************************** JETD '' Defaults WARNING!!!!! WARNING!!!!! You must use single quotes around the filename WARNING!!!!! e.g. JETD 'myfile.dat' WARNING!!!!! filename = Relative or absolute pathname that contains JETSET decay information a la LUUPDA. This info overrides the built-in decays. The file must contain the full decay table for a given particle but does not need to redefine every particle. ASLUND and BBKING users will recognize this as .newdat The LUCOMP mapping in the stand-alone and GEANT version is KF LUCOMP(KF) ----- ---------- 70533 404 Upsilon(4S) 80533 405 Upsilon(5S) nXXXXX 410+abs(n) user definable particles See Section 2.3.1 concerning LUCOMP in ASLUND **************************************************************************** BEAM Pz(electron) Pz(positron) Ecm (GeV) Set the beam energies **************************************************************************** BEAM 9.0 -3.1 0. Defaults Pz_electron, Pz_positron z component of momentum of electron and positron beam (GeV). Px = Py = 0. Ecm center-of-mass energy in (GeV) You may set EITHER the two beam momenta OR one beam momentum and the center-of-mass energy. Set the unused quantity to zero and it will automatically be filled in with the correct value. Note that the sign of the beam momenta is important. **************************************************************************** SBEAM Sigma(electron) Sigma(positron) (GeV) Set the beam energy spread **************************************************************************** SBEAM 0. 0. Defaults The beam smearing is handled as in the native ASLUND generator, which smears the nominal beam energy by a gaussian for every event. For each event the center-of-mass is calculated and the mass of the Upsilon(4S) is set equal to this energy. The Upsilon(4S) is then generated with zero momentum. The boost of the final state is also based on the smeared energies on an event by event basis. Warning, If the center-of-mass energy for a smeared event happens to be below the mass threshold for B0 production an empty event will be generated. The actual threshold in JETSET is slightly higher than twice the B0 mass. **************************************************************************** VERT x y z (mm) Set the primary vertex location **************************************************************************** VERT 0. 0. 0. Defaults The given coordinates, in millimeters, sets the primary vertex position **************************************************************************** SVERT dx dy dz (mm) Set the primary vertex spread **************************************************************************** SVERT 0. 0. 0. Defaults **************************************************************************** IVERT Ix Iy Iz Set the primary vertex smearing method **************************************************************************** IVERT 0 0 0 Defaults 0=smear with gaussian (sigma given by SVERT) 1=smear with flat distribution (-delta,delta), (delta given by SVERT **************************************************************************** GENP IHEP ILUND IGPKINE IGPVERT Print event information on the terminal **************************************************************************** GENP 0 0 0 0 Defaults IHEP >0 call heplist(IHEP) after every event to list particles ILUND>0 call lulist(ILUND) after every event to list particles IGPKINE>0 call GPKINE(0) after every event to list GEANT particles IGPVERT>0 call GPVERT(0) after every event to list GEANT vertices IGPKINE and IGPVERT work only in GEANT programs All four printouts can be activated simultaneously **************************************************************************** JETO KF1 KF2 ... KF Print the decay table for the listed particles **************************************************************************** JETO 0 0 ... 0 Defaults For any non-zero JETSET particle code the decay table for that particle will be listed at the end of the initialization pass through the JETSET generator. example: JETO 511 443 # list the decay table for B0 and J/PSI **************************************************************************** RAN Iseed1 Iseed2 Control the random number seeds **************************************************************************** RAN 1 0 Defaults This data card initializes the random number seed for the BEGET random number generator (RANECU). RAN behaves exactly like the RNDM data card of GEANT. The preselected random number seeds (see below) result in distinct sequences 1 billion numbers apart. Not all integer seeds result in a good or random sequence, so the authors recommend using the preselected ones. Iseed2=0 Iseed1 can be a number between 1 and 100 that causes a preselected seed to initialize the random number generator. Iseed2<>0 The random number generator will be initialized with seeds Iseed1 and Iseed2. **************************************************************************** EVNT Nbeg_events Set the number of events to generate **************************************************************************** EVNT 1 Defaults This data card merely sets a variable (Nbeg_events) in the BEGET common block. It is up to the user to check it and call the BEGET subroutine the appropriate number of times. ASLUND and BBSIM users should not use this card since event counting is already built in to those programs. **************************************************************************** REJE Ireje Activate user event rejection function **************************************************************************** REJE 0 Defaults This data card, in conjuction with the use function BEGUREJECT, allows the user to generate specific events by running a generator and then deciding whether the event is interesting. This is most useful in the context of sending events to GEANT, where simulation time is quite long. Ireje = 1 Call user function BEGUREJECT. If the functions returns true then generate another event call BEGUREJECT again. The event is not passed on to the main program until BEGUREJECT returns false. Ireje = 0 Do not reject any events.