########################################################### # NA62 MACRO TEMPLATE FOR THE INITIALIZATION PHASE # ########################################################### # # ¢ # Placeholders (these are set on a job-to-job basis): # MCSEED1 # MCRUNNUMBER # NUMBEROFEVENTS # MCOUTFILENAME # ## Output file ############################################ /output/fileName MCOUTFILENAME ## Beam Type ############################################## # /beam/SetBeam TYPE # Type can be # turtle    - NA62 beam # gps       - G4 Generic particle source; useful for tests # matrixgun - Used in the LAV matrices generation ########################################################### /beam/SetBeam turtle ## Matrixgun option ####################################### # /MatrixGun/particlePerBin N #   Set the number of particle generated in a event #   (Every event change matrix bin: see setBin command) # /MatrixGun/channelId N #   Select in witch LAV block generate the events # /MatrixGun/setBin N #   Select in witch matrix bin will generate the events #   If without argument (or negative argument) loop over bins # /MatrixGun/setRandomBin #   Select a random matrix bin to every generation # /MatrixGun/setFixedBin #   Use single bin for all generation # /MatrixGun/matrixLoad filename #   Select the file path of matrix to load. It could be also an #   empty matrix: the only relevant information is the binning ################################################################# ## Beam options ################################################# # # /beam/fileName: name of the datacard file for fast MC description. #                 Turtle is the external beam generator. #                 The datacard uses the standard turtle-type syntax. #                 The datacard repository is Beam/datacard. #                 Pre-defined datacards: # #         k12hika+_cedar_xx    : G4 starts in front of the CEDAR (default). #         k12hika+_gtk1_xx     : G4 starts in front of the GTK1 and Turtle #                                simulates the CEDAR material. #         k12hika+_bluetube_xx : G4 starts after the GTK3 and Turtle #                                simulates the CEDAR and GTK material #                                (for testing purposes). # # The suffix _xx of the datacard name identifies the type of the beam # particle. The available types are: #    kp kaon+ #    pp pion+ #    pr proton ########################################################################## /beam/fileName Beam/datacard/k12hika+_cedar_kp #################################################################################### # Control muon halo simulation # If /beam/Halo is On, ONLY the halo component before GTK3 is generated # or optionally, one can propagate through the setup the 525 MHz pi+ # component in the beam after GTK3. # The input binary file contains a structure with the kinematics of muons # from K+/- and pi+/- decays, or pions, produced by the HALO program # The user should specify the type of the particle (mu+, mu- or pi+) # Make sure that beamOn is larger than the number of particles in the bin file # # The following binary files are available in # /castor/cern.ch/user/b/balev/na62mc/halo # The file should be copied locally, since it is not directly readable from CASTOR # Note: the rate (in kHz) is equal to the number of particles in each file # # FILE                          PARTICLES    Description # ---------------------------   ----------   ------------------------ # muons_upstream_gtk3_pip.bin    99143 mu+   mu+ from pi+ before GTK3 # muons_upstream_gtk3_kap.bin    11791 mu+   mu+ from K+  before GTK3 # muons_upstream_gtk3_pim.bin    21374 mu-   mu- from pi- before GTK3 # muons_upstream_gtk3_kam.bin     3034 mu-   mu- from K-  before GTK3 # pions_downstream_gtk3.bin     526561 pi+   pi+ in beam after GTK3 #################################################################################### /beam/Halo Off /beam/HaloFile muons_upstream_gtk3_pip.bin /beam/HaloType mu+ ## Decay options ######################################################## # # /decay/force: # # 0   the beam particle decays according to the particle type # 1   the beam particle is a K+ and decays according to the specified #     decay type (default, see /decay/type, mandatory in this case) # # /decay/type: # # Hadronic modes # #  0  K+ --> pi+ nu nu #  1  K+ --> pi+ pi0 #  2  K+ --> pi+ pi0 gamma (IB) # # 10  K+ --> pi+ pi+ pi- #                 /decay/radcor:  0 not applied (default) #                                 1 Coulomb factor #                                 2 Coulomb factor + PHOTOS # # 11  K+ --> pi+ pi0 pi0, empirical parameterization, PLB686 (2010) 101 #                 /decay/radcor:  0 not applied (default) #                                 1 applied (PHOTOS) # # Leptonic modes Kl2(g) # # 20  K+ --> e+ nu       (2 body) # 21  K+ --> e+ nu gamma (IB) as in RK definition # 22  K+ --> e+ nu gamma (SD+) # 23  K+ --> e+ nu gamma (SD-) # 24  K+ --> e+ nu gamma (INT+) # 25  K+ --> e+ nu gamma (INT-) # # 30  K+ --> mu+ nu       (2 body) # 31  K+ --> mu+ nu gamma (IB) as in RK definition # 32  K+ --> mu+ nu gamma (SD+) # 33  K+ --> mu+ nu gamma (SD-) # 34  K+ --> mu+ nu gamma (INT+) # 35  K+ --> mu+ nu gamma (INT-) # # Semileptonic modes Kl3(g) # # 40  K+ --> pi0 e+ nu  (Ke3) #                 /decay/radcor:  0 not applied (default) #                                 1 applied (PHOTOS+Ginsberg) # 41  K+ --> pi0 mu+ nu (Kmu3) #                 /decay/radcor:  same as for mode 40 # 42  K+ --> pi0 e+ nu gamma  (from KLOE library) # 43  K+ --> pi0 mu+ nu gamma (from KLOE library) # # Semileptonic modes Kl4(g) # # 60  K+ --> pi+ pi- e+ nu (Ke4) #                 /decay/radcor:  0 not applied (default) #                                 1 applied (PHOTOS) # 61  K+ --> pi0 pi0 e+ nu (Ke4) #                 /decay/radcor:  same as for mode 60 # # The Ferrara 1-loop ChPT generators: not validated # # 70  K+ --> pi+ pi- e+ nu (Ke4) # 71  K+ --> pi0 pi0 e+ nu (Ke4) # 72  K+ --> pi+ pi- mu+ nu (Kmu4) # 73  K+ --> pi0 pi0 mu+ nu (Kmu4) # # Hadronic modes with photons # # 80  K+ --> pi+ pi0 gamma (IB), Egamma > 10 MeV in kaon frame # 81  K+ --> pi+ pi0 gamma (SD) # 82  K+ --> pi+ pi0 gamma (INT) # 83  K+ --> pi+ gamma gamma # 84  K+ --> pi+ gamma e+ e- # 85  K+ --> pi+ gamma mu+ mu- # # Four-lepton modes # # 100 K+ --> e+  nu e+  e- # 101 K+ --> mu+ nu e+  e- # 102 K+ --> e+  nu mu+ mu- # 103 K+ --> mu+ nu mu+ mu- # # FCNC, LFV and LNV modes # # 120 K+ --> pi+ e+ e- #                 /decay/radcor:  0 not applied (default) #                                 1 applied (PHOTOS) # 121 K+ --> pi+ mu+ mu- # # LFV and LNV modes: flat distribution over phase space # # 130 K+ --> pi+ mu+ e- # 131 K+ --> pi+ mu- e+ # 132 K+ --> pi- mu+ e+ # 133 K+ --> pi- e+  e+ # 134 K+ --> pi- mu+ mu+ # # /decay/pizeroDecay: governs the simulation of decays of pi0's from the kaon decay. #                     Effective only of there are pi0's in the final state. #                     The number of digits should be equal to the number of pi0's. #                     Each digit encodes the simulation of the decay of a pi0 #                     in the folowing way: 0: pi0 passed to Geant, #                     1: pi0-->gg; 2: pi0-->ge+e- (Dalitz), #                     3: pi0-->e+e-, 4: pi0-->e+e-e+e- (double Dalitz). #                     Default = 0. # # /decay/zMin, /decay/zMax: z coordinate in mm of the beginning and end of #                           the decay region. These values are effective #                           in case of forced decay only. ########################################################################## /decay/force 1 /decay/type 0 /decay/radcor 0 /decay/pizeroDecay 0 /decay/zMin 104000. /decay/zMax 180000. ## Detector options ###################################################### # All SubDetectors are enabled by default. # /Detector/EnableSubDetector , /Detector/DisableSubDetector # let enable/disable them individually # /Detector/UpdateGeometry must be issued after any change ########################################################################### #/Detector/DisableSubDetector CHANTI #/Detector/DisableSubDetector Cedar #/Detector/DisableSubDetector CHOD #/Detector/DisableSubDetector GigaTracker #/Detector/DisableSubDetector IRC #/Detector/DisableSubDetector LAV #/Detector/DisableSubDetector LKr #/Detector/DisableSubDetector RICH #/Detector/DisableSubDetector Spectrometer #/Detector/DisableSubDetector SAC #/Detector/DisableSubDetector MUV1 #/Detector/DisableSubDetector MUV2 #/Detector/DisableSubDetector MUV3 ### Cedar options #################################################### # # Enable the Cherenkov effect in the Cedar: true or false (default) # ###################################################################### /Detector/Cedar/EnableCherenkovEffect true ## LAV options ############################################################################ # # ### Select which parts of the LAV should be used ### # All stations/vessels/layers/bananas are enabled by default # # # Stations #   - goes from 0 (LAV station 1) to 11 (LAV station 12) # /Detector/LAV/EnableAllStations # /Detector/LAV/DisableAllStations # /Detector/LAV/EnableStation # /Detector/LAV/DisableStation # # Vessels #   - as for Stations # /Detector/LAV/EnableAllVessels # /Detector/LAV/DisableAllVessels # /Detector/LAV/EnableVessel # /Detector/LAV/DisableVessel # # Layers #   - as for Stations #   - goes from 0 to 4 for stations 1-8, from 0-3 for stations 9-12 # /Detector/LAV/EnableAllLayers # /Detector/LAV/DisableAllLayers # /Detector/LAV/EnableLayer # /Detector/LAV/DisableLayer # # Bananas #   - as for Stations #   - as for Layers #   - : 0-7 for stations 1-5, 0-11 for stations 6-8, #          0-14 for stations 9-11, 0-15 for station 12 # /Detector/LAV/EnableAllBananas # /Detector/LAV/DisableAllBananas # /Detector/LAV/EnableBanana # /Detector/LAV/DisableBanana # # Example: use only LAV station 3 (=2) # /Detector/LAV/DisableAllStations # /Detector/LAV/EnableStation 2 # # Example: use only layer 4 of station 3 (=2), do not use steel vessel # /Detector/LAV/DisableAllStations # /Detector/LAV/DisableAllLayers # /Detector/LAV/EnableStation 2 # /Detector/LAV/EnableLayer 2 4 # /Detector/LAV/DisableVessel 2 # # ### Select how the PbGl block should be simulated #   - can be: #     standard - approximate optical photons tracking using map files (fast, default) #     optical  - do full tracking of optical photons (slow) # /Detector/LAV/SelectBlockSimulation # # ### Select matrices for standard simulation ### #/Detector/LAV/EfficiencyMatrix ./LAV/LAVEff.txt #/Detector/LAV/TimeMatrix ./LAV/LAVDelay.txt ############################################################################################# ### Enable RICH fast simulation ######################################### # # Full simulation enabled by default: # /Detector/UpdateGeometry must be issued after any change # # Example: use Fast Simulation for the RICH: # /Detector/RICH/EnableFastSimulation # #   = false - standard G4 simulation (slow, default) #           true  - No refractive index in neon, optical photons generation and transport (fast) # ################################################################### /Detector/RICH/EnableFastSimulation false ### Enable MUV1 full optical simulation ######################################### # # Full optical simulation disabled by default: # #   = false - No photons are produced, saving information for parametrization  (default) #           true  - Full optical simulation (very slow, not recommended) # ################################################################### /Detector/MUV1/EnableFullOpticalSimulation false /Detector/UpdateGeometry ### User beam (for testing purposes) ################################################### #/beam/SetBeam gps ##/gps/particle e+ ##/gps/position 0 0 215 m ##/gps/momentumAmp 10.0 GeV ##/gps/direction 0 0 1 ## 1D accelerator beam ## ##/gps/verbose 1 #/gps/particle e+ #/gps/pos/type Beam ## ## the incident surface is in the x-y plane #/gps/pos/rot1 1 0 0 #/gps/pos/rot2 0 1 0 ## ## the beam spot is centered at few cm upstream the RICH and is ## of 1d gaussian shape with a 1mm central plateau #/gps/pos/shape Circle #/gps/pos/centre  0 0.2 215 m #/gps/pos/radius 1. mm #/gps/pos/sigma_r .1 mm ## ## the beam is travelling along the z-axis with 20 microrad dispersion #/gps/ang/rot1 0 1 0 #/gps/ang/rot2 1 0 0 #/gps/ang/type beam1d #/gps/ang/sigma_r 20e-6 rad ## ## the beam energy is in gaussian profile #/gps/ene/type Gauss #/gps/ene/mono 10.0 GeV #/gps/ene/sigma 100. MeV # ##/tracking/verbose 1 ########################################################################################## ## Handling of KineParts ################################################################# ## The beam particle and two generations of its decay daughters are always stored ## The following commands are listed in order of priority of execution ## 1. Save all steps or tracks if 1 (default 0) ##    If each step is saved, then the track is not saved /kine/SaveAllSteps 0 /kine/SaveAllTracks 0 ## 2. Tracks which always should be kept ##    - by particle name ##    - by creator process ##    - by end process ##    By default no tracks are saved by those commands #/kine/ParticleTypeToSave pi+ #/kine/CreatorProcessToSave Decay #/kine/EndProcessToSave Decay ## 3. Tracks which always should be rejected ##    - by particle name ##    - by creator process ##    - by end process ##    By default optical photons, products of pi,K inelastic scattering ##    and bremsstrahlung, and scintillating particles are not stored /kine/ParticleTypeToReject opticalphoton /kine/CreatorProcessToReject PionPlusInelastic /kine/CreatorProcessToReject PionMinusInelastic /kine/CreatorProcessToReject KaonPlusInelastic /kine/CreatorProcessToReject eBrem /kine/CreatorProcessToReject muIoni /kine/EndProcessToReject Scintillation ## 4. For the rest of the tracks save if ##    - energy is above MinimumEnergy (in MeV) ##    - interaction level is less or equal to MaximumInteractionLevel ##    Interaction level is 0 for the beam and increments by 1 for each secondary particle ##    The default energy cut is 0 (all particle stored) ##    The default interaction level is 2 /kine/MinimumEnergy 0. /kine/MaximumInteractionLevel 2 ## Print information for saved tracks for each event if 1 (default 0) /kine/SaveTrackVerbose 0 ## Random Generator ###################################################### /random/seedDecay MCSEED1 ## Run Number ###################################################### /run/number MCRUNNUMBER ## Number of events ###################################################### /run/beamOn NUMBEROFEVENTS