From cfa787fe0b1d1cbfde05b6876937a21b82048fd1 Mon Sep 17 00:00:00 2001 From: Voutsi Date: Wed, 9 Mar 2022 16:20:57 +0100 Subject: [PATCH] Update config_LST1_muons.cfg Cleaning up --- config_LST1_muons.cfg | 108 +++--------------------------------------- 1 file changed, 6 insertions(+), 102 deletions(-) diff --git a/config_LST1_muons.cfg b/config_LST1_muons.cfg index 8e776cf..6695607 100644 --- a/config_LST1_muons.cfg +++ b/config_LST1_muons.cfg @@ -24,84 +24,17 @@ SEED @seed4@ 323 0 // Seed for 4th random number sequence, t * [ Primary particle options ] * PRMPAR 6 // mu- -* #ifdef PRMPAR -* PRMPAR $(PRMPAR) -* #endif -* -* #ifndef ONAXIS -* #define DIFFUSE 1 -* #endif -* ERANGE 5.832 1E3 // Energy range of primary particle (in GeV): muons with 50% opening angle or more -ERANGE 8.418 1E3 // Energy range of primary particle (in GeV): muons with 80% opening angle or more -* #ifdef EMIN -* # ifdef EMAX -* ERANGE $(EMIN)E3 $(EMAX)E3 // Requires EMIN and EMAX in units of TeV. -* # endif -* #endif +* ERANGE 8.418 1E3 // Energy range of primary particle (in GeV): muons with 80% opening angle or more * ESLOPE -2.0 // Slope of primary energy spectrum (-2.0 is equal CPU time per decade) -* #ifdef ESLOPE -* ESLOPE $(ESLOPE) // Requires spectral slope ESLOPE (<0) -* #endif -* -NSHOW 1000 // number of showers to generate -* #ifdef NSHOW -* NSHOW $(NSHOW) // Requires NSHOW environment variable. -* #endif -* -THETAP 0. 0. // Range of zenith angles (degrees) -* #ifdef ZENITH_ANGLE -* THETAP $(ZENITH_ANGLE) $(ZENITH_ANGLE) // Range of zenith angles (degrees) -* #endif -* -PHIP 180. 180. // Range of azimuth angles (degree): primaries coming from North -* #ifdef FROM_SOUTH -PHIP 0. 0. // Range of azimuth angles (degree): primaries coming from South -* #else -* # ifdef FROM_EAST -PHIP 90. 90. // Range of azimuth angles (degree): primaries coming from East -* # else -* # ifdef FROM_WEST -PHIP 270. 270. // Range of azimuth angles (degree): primaries coming from West -* # else -PHIP 180. 180. // Range of azimuth angles (degree): primaries coming from North -* # endif -* # endif -* #endif -* -* #if defined(ALIGN_B_FIELD) -* Geomagnetic field is assumed to be aligned with geographic North (array x axis). -* # define AZM_FROM_SOUTH 0. -* # define AZM_FROM_EAST 90. -* # define AZM_FROM_NORTH 180. -* # define AZM_FROM_WEST 270. -* #else -* Geomagnetic field aligned as predicted for year 2020.0. -* # define AZM_FROM_SOUTH -5.3195 -* # define AZM_FROM_EAST 84.6805 -* # define AZM_FROM_NORTH 174.6805 -* # define AZM_FROM_WEST 264.6805 -* #endif -* -* #if defined(FROM_SOUTH) -* PHIP $(AZM_FROM_SOUTH) $(AZM_FROM_SOUTH) // CORSIKA azimuth angles (degree) for primaries coming from South. -* IACT setenv AZM 180 // Corresponding astronomical azimuth, from geographical North towards East. -* #elif defined(FROM_EAST) -* PHIP $(AZM_FROM_EAST) $(AZM_FROM_EAST) // CORSIKA azimuth angles (degree) for primaries coming from East. -* IACT setenv AZM 90 // Corresponding astronomical azimuth, from geographical North towards East. -* #elif defined(FROM_WEST) -* PHIP $(AZM_FROM_WEST) $(AZM_FROM_WEST) // CORSIKA azimuth angles (degree) for primaries coming from West. -* IACT setenv AZM 270 // Corresponding astronomical azimuth, from geographical North towards East. -* #else -* PHIP $(AZM_FROM_NORTH) $(AZM_FROM_NORTH) // CORSIKA azimuth angles (degree) for primaries coming from North. -* IACT setenv AZM 0 // Corresponding astronomical azimuth, from geographical North towards East. -* #endif -* -*# ifdef TEL_LST +* NSHOW 1000 // number of showers to generate +* +* THETAP 0. 0. // Range of zenith angles (degrees) +* PHIP 180. 180. // Range of azimuth angles (degree): primaries coming from North +* * VIEWCONE 0. 2.0 // Such that at half of maximum Cherenkov angle the ring still fits into FoV. (LST: actual FoV-diam=4.6 deg) * VIEWCONE 0. 1.1 // Such that at maximum Cherenkov angle the ring still fits into FoV. (LST: actual FoV-diam=4.6 deg) * VIEWCONE 0. 0.9 // Such that at maximum Cherenkov angle the ring still fits into FoV. (LST: required FoV-diam=4.2 deg) -*# endif * * Optionally override prepared demo run settings NSHOW @NShower@ @@ -132,8 +65,6 @@ ATMOSPHERE 99 Y // Use a custom profile by filename (this must come _AFTER_ * * ./geomag70 IGRF12.COF 2020.00 D K2.180 28.7621661 -17.8920302 * -* MAGNET 30.576 23.571 // La Palma, 2020.0 -* ARRANG $(AZM_FROM_SOUTH) // La Palma, 2020.0 MAGNET 30.678 23.415 // La Palma, 01/12/2021 ARRANG -4.84 // La Palma, 01/12/2021 * @@ -158,13 +89,6 @@ ARRANG -4.84 // La Palma, 01/12/2021 * * [ Core range ] * -* #if defined(NSCAT) && defined(CSCAT) -* CSCAT $(NSCAT) $(CSCAT)E2 0. // Requires NSCAT (10/20) and CSCAT (in units of meters). -* #elif defined(DIFFUSE) -CSCAT 20 1900E2 0. // Use shower several times (protons+electrons+..., larger area for diffuse origin) -* #else -CSCAT 10 1400E2 0. // Use shower several times (gammas, point source only) -* #endif * * * [ Telescope positions, for IACT option ] @@ -182,14 +106,7 @@ CSCAT 10 1400E2 0. // Use shower several times (gammas, point source only) * * TELESCOPE 0.0E2 0.0E2 $(LSTHT) $(LSTRAD) # LST TELESCOPE 0.0E2 0.0E2 16.00E2 12.50E2 // LST -* CSCAT 1 12.3E2 0. # Just beyond dish -* CSCAT 1 11E2 0. # Close to edge of dish CSCAT 1 9.8E2 0. # To 80% of dish radius -* #ifdef NSCAT -* # ifdef CSCAT -* CSCAT $(NSCAT) $(CSCAT)E2 0. // Requires NSCAT (1/10/20) and CSCAT (in units of meters). -* # endif -* #endif * * [Interaction flags] * @@ -208,8 +125,6 @@ ECTMAP 1.E6 // Cut on gamma factor for printout * [ Cherenkov emission parameters ] * CERSIZ 5. // Not above 10 for super/ultra-bialkali QE; 7 is fairly OK; 5 should be safe. -* CERFIL F // No old-style Cherenkov output to extra file -* Fixing 'CERFIL F' for CORSIKA version 77100: CERFIL 0 CWAVLG 230. 900. // Cherenkov wavelength band * @@ -219,19 +134,8 @@ DEBUG F 6 F 1000000 // Debug flag and logical unit for output DATBAS yes // Write a file with parameters used DIRECT ./ PAROUT F F -* #ifdef WITHOUT_MULTIPIPE -* # ifdef STRICT_BASELINE -* TELFIL cta-prod5-lapalma-baseline.corsika.gz // If telescope simulation not done directly in pipe -TELFIL run${RUNNR}_${PRMNAME}_za${ZA}deg_azm${AZM}deg-lapalma-baseline${extra_suffix2}.corsika.zst -* # else -TELFIL run${RUNNR}_${PRMNAME}_za${ZA}deg_azm${AZM}deg-lapalma${extra_suffix2}.corsika.zst -* # endif -* #else -* TELFIL |${SIM_TELARRAY_PATH}/run_sim_cta // Telescope photon bunch output (eventio format) -* TELFIL @output_filename@_muon.corsika.gz // Telescope photon bunch output (eventio format) TELFIL run${RUNNR}_muon.corsika.gz IACT TELOPT -c cta-prod5-lapalma-baseline -* #endif IACT PRINT_EVENTS 100 100 1 * * [ IACT tuning parameters ]