note:

Next to classy, you need to pull and pip install this: https://bitbucket.org/aurelschneider/genmassfct/src/master/

1. Note : predictions of the haloes in WDM/CDM
2. 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

1. 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

2. plots for Oliver #######

/home/revaz/ownCloud/papers/2023M83