Skip to content

Commit

Permalink
Update config_LST1_muons.cfg
Browse files Browse the repository at this point in the history 
Cleaning up
  • Loading branch information
@Voutsi
Voutsi authored Mar 9, 2022
1 parent b35cab8 commit cfa787f
Showing 1 changed file with 6 additions and 102 deletions.
108 changes: 6 additions & 102 deletions config_LST1_muons.cfg
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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@
Expand Down Expand Up @@ -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
*
Expand All @@ -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 ]
Expand All @@ -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]
*
Expand All @@ -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
*
Expand All @@ -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 ]
Expand Down

0 comments on commit cfa787f

Please sign in to comment.