Homec4science
Diffusion ARRAKIHS (master) Browse

arrakihs/Workspace/CumNumOfSat/Readme
d192e801c42bmaster

Last Change
Raw File

Readme
View Options

# note:
Next to classy, you need to pull and pip install this:
https://bitbucket.org/aurelschneider/genmassfct/src/master/
# Note : predictions of the haloes in WDM/CDM
# https://arxiv.org/abs/2304.09810
# test the halo function distribution (pour Courbin)
./satComputeNvsM.py
# test the Mstar vs Mhalo relation
./satComputeMstarvsMhalo.py --LvvsMh_model rj2018
./satComputeMstarvsMhalo.py --LvvsMh_model model1 -p LvvsMhm1.png
./satComputeMstarvsMhalo.py --LvvsMh_model model2 -p LvvsMhm2.png
./satComputeMstarvsMhalo.py --LvvsMh_model model3 -p LvvsMhm3.png
./satComputeMstarvsMhalo.py --LvvsMh_model model4 -p LvvsMhm4.png
# do Ngal realisation and plot the cumulative numbers of satellites
./satComputeNCum.py --Ngal 80
./satComputeNCumWithErrors.py --Nrealisations=100 --Ngal 100 -p NhCum100.png
./satComputeNCumWithErrors.py --Nrealisations=100 --Ngal 75 -p NhCum075.png
./satComputeNCumWithErrors.py --Nrealisations=100 --Ngal 50 -p NhCum050.png
./satComputeNCumWithErrors.py --Nrealisations=100 --Ngal 25 -p NhCum025.png
# do Ngal realisation and plot the cumulative numbers of satellites vs Luminosity
./satComputeNCumvsLv.py --LvvsMh_file MsvsMh_2018.pkl --Ngal 80
# same but now we run several realisations
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model1 --data --Nrealisations=100 --Ngal 100 -p LCum100m1.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model1 --data --Nrealisations=100 --Ngal 75 -p LCum075m1.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model1 --data --Nrealisations=100 --Ngal 50 -p LCum050m1.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model1 --data --Nrealisations=100 --Ngal 25 -p LCum025m1.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model2 --data --Nrealisations=100 --Ngal 100 -p LCum100m2.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model2 --data --Nrealisations=100 --Ngal 75 -p LCum075m2.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model2 --data --Nrealisations=100 --Ngal 50 -p LCum050m2.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model2 --data --Nrealisations=100 --Ngal 25 -p LCum025m2.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 100 -p LCum100m3.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 -p LCum075m3.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 50 -p LCum050m3.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 25 -p LCum025m3.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model4 --data --Nrealisations=100 --Ngal 100 -p LCum100m4.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model4 --data --Nrealisations=100 --Ngal 75 -p LCum075m4.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model4 --data --Nrealisations=100 --Ngal 50 -p LCum050m4.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model4 --data --Nrealisations=100 --Ngal 25 -p LCum025m4.png
# variation of the DMmass
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 1.1 -p LCum075m3DM1.1.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 2.0 -p LCum075m3DM2.0.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 3.0 -p LCum075m3DM3.0.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 4.0 -p LCum075m3DM4.0.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 5.0 -p LCum075m3DM5.0.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 6.0 -p LCum075m3DM6.0.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 7.0 -p LCum075m3DM7.0.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 8.0 -p LCum075m3DM8.0.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 9.0 -p LCum075m3DM9.0.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 10.0 -p LCum075m3DM10.0.png
# distribution of Main halo mass
./satComputeNCumvsLvWithErrorsWithM0varied.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 2.0 --M0_var=0 -p LCum075m3DM2.0.d0.0.png
./satComputeNCumvsLvWithErrorsWithM0varied.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 2.0 --M0_var=0.1 -p LCum075m3DM2.0.d0.1.png
./satComputeNCumvsLvWithErrorsWithM0varied.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 2.0 --M0_var=0.5 -p LCum075m3DM2.0.d0.5.png
./satComputeNCumvsLvWithErrorsWithM0varied.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --DMmass 2.0 --M0_var=1.0 -p LCum075m3DM2.0.d1.0.png
# in magnitude
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --filter VISeuclid --DMmass 1.1 -p LCum_VIS_075m3DM1.1.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --filter VISeuclid --DMmass 2.0 -p LCum_VIS_075m3DM2.0.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --filter VISeuclid --DMmass 3.0 -p LCum_VIS_075m3DM3.0.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 75 --filter VISeuclid --DMmass 4.0 -p LCum_VIS_075m3DM4.0.png
./satComputeNCumvsLvWithErrors.py --LvvsMh_model model3 --data --Nrealisations=100 --Ngal 30 --filter VISeuclid --DMmass 4.0 -p LCum_VIS_075m3DM4.0.png
# comments
- no variation of the MW mass, in all realisation, we assume the same total mass and take
values from Sawala 2017
- the Lv vs Mh relation is obtained from CDM simulations (no WDM)
- the Lv vs Mh relation is quite uncertain
-> f : Subhalo mass function (SHMF) supression
see :
Nadler 2019 -> eq 7 - from Schneider 2012
Lovell et al. (2014),
dmo resolution (Forouhar): 4 × 10^5
###########################
# plots for Oliver
###########################
/home/revaz/ownCloud/papers/2023M83
###########################
# tests of statistics
###########################
#################
# 75 galaxies
#################
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e5 --DMmass -1 -o DM-1.00_75_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e5 --DMmass 5.0 -o DM05.00_75_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e5 --DMmass 7.0 -o DM07.00_75_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e5 --DMmass 10.0 -o DM10.00_75_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e6 --DMmass -1 -o DM-1.00_75_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e6 --DMmass 5.0 -o DM05.00_75_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e6 --DMmass 7.0 -o DM07.00_75_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e6 --DMmass 10.0 -o DM10.00_75_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e7 --DMmass -1 -o DM-1.00_75_1e7.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e7 --DMmass 5.0 -o DM05.00_75_1e7.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e7 --DMmass 7.0 -o DM07.00_75_1e7.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 75 --Lmin 1e7 --DMmass 10.0 -o DM10.00_75_1e7.pkl
./satPlotStatistics.py DM*_75_1e5.pkl
./satPlotStatistics.py DM*_75_1e6.pkl
./satPlotStatistics.py DM*_75_1e7.pkl
#################
# 150 galaxies
#################
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e5 --DMmass -1 -o DM-1.00_150_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e5 --DMmass 5.0 -o DM05.00_150_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e5 --DMmass 7.0 -o DM07.00_150_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e5 --DMmass 10.0 -o DM10.00_150_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e6 --DMmass -1 -o DM-1.00_150_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e6 --DMmass 5.0 -o DM05.00_150_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e6 --DMmass 7.0 -o DM07.00_150_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e6 --DMmass 10.0 -o DM10.00_150_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e7 --DMmass -1 -o DM-1.00_150_1e7.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e7 --DMmass 5.0 -o DM05.00_150_1e7.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e7 --DMmass 7.0 -o DM07.00_150_1e7.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 150 --Lmin 1e7 --DMmass 10.0 -o DM10.00_150_1e7.pkl
./satPlotStatistics.py DM*_150_1e5.pkl
./satPlotStatistics.py DM*_150_1e6.pkl
./satPlotStatistics.py DM*_150_1e7.pkl
#################
# 150 galaxies
#################
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e5 --DMmass -1 -o DM-1.00_15_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e5 --DMmass 5.0 -o DM05.00_15_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e5 --DMmass 7.0 -o DM07.00_15_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e5 --DMmass 10.0 -o DM10.00_15_1e5.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e6 --DMmass -1 -o DM-1.00_15_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e6 --DMmass 5.0 -o DM05.00_15_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e6 --DMmass 7.0 -o DM07.00_15_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e6 --DMmass 10.0 -o DM10.00_15_1e6.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e7 --DMmass -1 -o DM-1.00_15_1e7.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e7 --DMmass 5.0 -o DM05.00_15_1e7.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e7 --DMmass 7.0 -o DM07.00_15_1e7.pkl
./satComputeStatistics.py --LvvsMh_model model3 --Nrealisations=1000 --Ngal 15 --Lmin 1e7 --DMmass 10.0 -o DM10.00_15_1e7.pkl
./satPlotStatistics.py DM*_15_1e5.pkl
./satPlotStatistics.py DM*_15_1e6.pkl
./satPlotStatistics.py DM*_15_1e7.pkl
# other model
./satComputeStatistics.py --LvvsMh_model model1 --Nrealisations=1000 --Ngal 75 --Lmin 1e6 --DMmass -1 -o DM-1.00_75_1e6_M1.pkl
./satComputeStatistics.py --LvvsMh_model model1 --Nrealisations=1000 --Ngal 75 --Lmin 1e6 --DMmass 5.0 -o DM05.00_75_1e6_M1.pkl
c4science · Help