21cm_gal_uncertainties/do_calculations_new_aa4 at main · annehutter/21cm_gal_uncertainties · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
#! /bin/bash
if [ $# -lt 0 ]; then
    echo "Usage: ./do_calculations"
else

    redshift=$1
    frequency=$2
    crosscorrXHI=$3
    crosscorrDeltaHI=$4
    crosscorrRegionMean=$5
    crosscorrRegionStd=$6
    galMinLum=$7
    galSigmaRedshift=$8
    surveyFieldOfView=$9
    simulationBaseDir=${10}
    simulationLAEDir=${11}
    simulationName=${12}
    snapshot=${13}
    outputDir=${14}

    snapshotIon=$((snapshot - 1))

    echo $simulationName

#   ----------------------------------------------------------
#   Create directory
#   ----------------------------------------------------------

    thisDirName=$simulationName'_galMinLum'$galMinLum'_LAEsigmaz'$galSigmaRedshift'_FoV'$surveyFieldOfView

    if [ ! -d $outputDir ]; then
      mkdir $outputDir
    fi
    cd $outputDir
    if [ ! -d $thisDirName ]; then
      mkdir $thisDirName
    fi
    cd $thisDirName
    cp ../../sample_iniFile/iniFile.ini iniFile.ini

#   ----------------------------------------------------------
#   Adjust iniFile.ini
#   ----------------------------------------------------------

    oldStringRedshift='redshift = 7.0'
    newStringRedshift='redshift = '$redshift
    sed -i "s+${oldStringRedshift}+${newStringRedshift}+g" iniFile.ini

    #--------------------------------------
    # For estimates from analytic functions (units are Mpc/h)
    #--------------------------------------

    oldStringCrosscorrXHI='crosscorrXHI = 0.49'
    newStringCrosscorrXHI='crosscorrXHI = '$crosscorrXHI
    sed -i "s+${oldStringCrosscorrXHI}+${newStringCrosscorrXHI}+g" iniFile.ini

    oldStringCrosscorrDeltaHI='crosscorrDeltaHI = 0.7'
    newStringCrosscorrDeltaHI='crosscorrDeltaHI = '$crosscorrDeltaHI
    sed -i "s+${oldStringCrosscorrDeltaHI}+${newStringCrosscorrDeltaHI}+g" iniFile.ini

    oldStringCrosscorrRegionMean='crosscorrRegionMean = 5.48'
    newStringCrosscorrRegionMean='crosscorrRegionMean = '$crosscorrRegionMean
    sed -i "s+${oldStringCrosscorrRegionMean}+${newStringCrosscorrRegionMean}+g" iniFile.ini

    oldStringCrosscorrRegionStd='crosscorrRegionStd = 0.80'
    newStringCrosscorrRegionStd='crosscorrRegionStd = '$crosscorrRegionStd
    sed -i "s+${oldStringCrosscorrRegionStd}+${newStringCrosscorrRegionStd}+g" iniFile.ini

    oldStringGalSigmaRedshift='galSigmaRedshift = 0.01'
    newStringGalSigmaRedshift='galSigmaRedshift = '$galSigmaRedshift
    sed -i "s+${oldStringGalSigmaRedshift}+${newStringGalSigmaRedshift}+g" iniFile.ini

    #--------------------------------------
    # Survey specifics
    #--------------------------------------

    oldStringSurveyFieldOfView='surveyFieldOfView = 100.'
    newStringSurveyFieldOfView='surveyFieldOfView = '$surveyFieldOfView
    sed -i "s+${oldStringSurveyFieldOfView}+${newStringSurveyFieldOfView}+g" iniFile.ini

    oldStringSurveyFile='surveyFile = input/SKA1_Low.drift_blmin0_blmax5003_0.177GHz_arrayfile.npz'
    newStringSurveyFile='surveyFile = input/SKA1_Lowaa4.track_6.0hr_blmin0_blmax1701_'$frequency'GHz_arrayfile.npz'
    sed -i "s+${oldStringSurveyFile}+${newStringSurveyFile}+g" iniFile.ini

    #--------------------------------------
    # For estimates from simulations (default units are in Mpc/h)
    #--------------------------------------

    oldStringTbFile='TbFile = /home/hutter/PostdocDAWN/Projects/Astraeus_LAEs_2023/MHDEC/analysis/galaxy_analysis/LAEs/LyaProfileGAUSS_LyaLum1.e42_dustSIMPLE_faTofc1.0/crosscorrfunc_21cm_LAE/crosscorrfunc_21cm_LAE_z7.00.dat_21'
    newStringTbFile='TbFile = '$simulationBaseDir'/'$simulationLAEDir'/crosscorrfunc_21cm_LAE/crosscorrfunc_21cm_LAE_z'$redshift'.dat_21'
    sed -i "s+${oldStringTbFile}+${newStringTbFile}+g" iniFile.ini

    oldStringIonFile='ionFile = /home/hutter/PostdocDAWN/Projects/Astraeus_LAEs_2023/MHDEC/XHII.out_55'
    newStringIonFile='ionFile = '$simulationBaseDir'/XHII.out_'$snapshotIon
    sed -i "s+${oldStringIonFile}+${newStringIonFile}+g" iniFile.ini

    oldStringDensFile='densFile = /home/hutter/PostdocDAWN/Projects/MD_simulations/VSMD_160mpc_3840/fields/512/dmdens_cic_512.dat_056'
    newStringDensFile='densFile = /home/hutter/PostdocDAWN/Projects/MD_simulations/VSMD_160mpc_3840/fields/512/dmdens_cic_512.dat_0'$snapshot
    sed -i "s+${oldStringDensFile}+${newStringDensFile}+g" iniFile.ini

    oldStringGalFile='galFile = /home/hutter/PostdocDAWN/Projects/Astraeus_LAEs_2023/MHDEC/analysis/galaxy_analysis/LAEs/LyaProfileGAUSS_LyaLum1.e42_dustSIMPLE_faTofc1.0/Tables'
    newStringGalFile='galFile = '$simulationBaseDir'/'$simulationLAEDir'/Tables'
    sed -i "s+${oldStringGalFile}+${newStringGalFile}+g" iniFile.ini

    oldStringGalMinLum='galMinLum = 6.e42'
    newStringGalMinLum='galMinLum = '$galMinLum
    sed -i "s+${oldStringGalMinLum}+${newStringGalMinLum}+g" iniFile.ini

    #--------------------------------------
    # Specifics for statistics
    #--------------------------------------

    #--------------------------------------
    # Output directory and base name
    #--------------------------------------

    oldStringOutfile='output = output/z7.0'
    newStringOutfile='output = '$outputDir'/'$thisDirName'/z'$redshift
    sed -i "s+${oldStringOutfile}+${newStringOutfile}+g" iniFile.ini

    oldStringOutfileStat='outputStat = output/z7.0'
    newStringOutfileStat='outputStat = '$outputDir'/'$thisDirName'/z'$redshift
    sed -i "s+${oldStringOutfileStat}+${newStringOutfileStat}+g" iniFile.ini

#   ----------------------------------------------------------
#   Move iniFile.ini and execute
#   ----------------------------------------------------------

    mv iniFile.ini iniFile_z$redshift.ini
    cd ../..

    python src/main.py $outputDir/$thisDirName/iniFile_z$redshift.ini

#   ----------------------------------------------------------
#   Make plots
#   ----------------------------------------------------------

    python plots/crosspowerspec.py $outputDir/$thisDirName $outputDir/$thisDirName/crossps_21cm-LAE.png
    python plots/covariances.py $outputDir/$thisDirName $outputDir/$thisDirName/covariance_21cm-LAE.png
#     python plots/crosscorrfunc.py $outputDir/$thisDirName $outputDir/$thisDirName/crosscorrfunc_21cm-LAE.png

fi