Merge pull request #626 from galacticusorg/featJHLFAnalysis

Add an `outputAnalysis` class to compute high-z H and J-band luminosity functions

doc/Galacticus.bib

@@ -7100,6 +7100,25 @@ @article{srisawat_sussing_2013
 	file = {Sussing Merger Trees: The Merger Trees Comparison Project:/home/abenson/.mozilla/firefox/f54gqgdx.default/zotero/storage/WGFQ8UIZ/stt1545.html:text/html}
 }
 
+@article{stefanon_evolution_2013,
+	title = {The evolution of the rest-frame {J}- and {H}-band luminosity function of galaxies to z=3.5},
+	volume = {429},
+	issn = {0035-8711},
+	url = {https://ui.adsabs.harvard.edu/abs/2013MNRAS.429..881S},
+	doi = {10.1093/mnras/sts388},
+	abstract = {We present the rest-frame J- and H-band luminosity function (LF) of field galaxies, based on a deep multiwavelength composite sample from the Multiwavelength Survey by Yale-Chile, Faint Infrared Extragalactic Survey and GOODS Chandra Deep Field-South survey public catalogues, covering a total area of 450 arcmin2. The availability of flux measurements in the Spitzer Infrared Array Camera 3.6, 4.5, 5.8 and 8 μm channels allows us to compute absolute magnitudes in the rest-frame J and H bands up to z = 3.5 minimizing the dependence on the stellar evolution models. We compute the LF in the four redshift bins: 1.5 {\textless} z {\textless} 2.0; 2.0 {\textless} z {\textless} 2.5; 2.5 {\textless} z {\textless} 3.0 and 3.0 {\textless} z {\textless} 3.5. Combining our results with those already available at lower redshifts, we find that (1) the faint-end slope is consistent with being constant up to z = 3.5, with α = -1.05 ± 0.03 for the rest-frame J band and α = -1.15 ± 0.02 for the rest-frame H band; (2) the normalization φ* decreases by a factor of 4-6 between z = 0 and z ≃ 1.75 and by a factor of 2-3 between z ≃ 1.75 and z = 3.25; (3) the characteristic magnitude M* shows a brightening from z = 0 to z ≃ 2 followed by a slow dimming to z = 3.25. We finally compute the luminosity density (LD) in both rest-frame J and H bands. The analysis of our results together with those available in the literature shows that the LD is approximately constant up to z ≈ 1, and it then decreases by a factor of 6 up to z = 3.5.},
+	urldate = {2024-05-23},
+	journal = {Monthly Notices of the Royal Astronomical Society},
+	author = {Stefanon, Mauro and Marchesini, Danilo},
+	month = feb,
+	year = {2013},
+	note = {Publisher: OUP
+ADS Bibcode: 2013MNRAS.429..881S},
+	keywords = {Astrophysics - Cosmology and Nongalactic Astrophysics, galaxies: evolution, galaxies: high-redshift, galaxies: luminosity function, galaxies: statistics, infrared: galaxies, mass function},
+	pages = {881--894},
+	file = {Full Text PDF:/home/abensonca/.mozilla/firefox/f54gqgdx.default/zotero/storage/XZLM9DV5/Stefanon and Marchesini - 2013 - The evolution of the rest-frame J- and H-band lumi.pdf:application/pdf},
+}
+
 @article{story_measurement_2012,
        author = {{Story}, K.~T. and {Reichardt}, C.~L. and {Hou}, Z. and {Keisler}, R. and {Aird}, K.~A. and {Benson}, B.~A. and {Bleem}, L.~E. and {Carlstrom}, J.~E. and {Chang}, C.~L. and {Cho}, H. -M. and {Crawford}, T.~M. and {Crites}, A.~T. and {de Haan}, T. and {Dobbs}, M.~A. and {Dudley}, J. and {Follin}, B. and {George}, E.~M. and {Halverson}, N.~W. and {Holder}, G.~P. and {Holzapfel}, W.~L. and {Hoover}, S. and {Hrubes}, J.~D. and {Joy}, M. and {Knox}, L. and {Lee}, A.~T. and {Leitch}, E.~M. and {Lueker}, M. and {Luong-Van}, D. and {McMahon}, J.~J. and {Mehl}, J. and {Meyer}, S.~S. and {Millea}, M. and {Mohr}, J.~J. and {Montroy}, T.~E. and {Padin}, S. and {Plagge}, T. and {Pryke}, C. and {Ruhl}, J.~E. and {Sayre}, J.~T. and {Schaffer}, K.~K. and {Shaw}, L. and {Shirokoff}, E. and {Spieler}, H.~G. and {Staniszewski}, Z. and {Stark}, A.~A. and {van Engelen}, A. and {Vanderlinde}, K. and {Vieira}, J.~D. and {Williamson}, R. and {Zahn}, O.},
         title = "{A Measurement of the Cosmic Microwave Background Damping Tail from the 2500-Square-Degree SPT-SZ Survey}",

source/objects.stellar_luminosities.F90

@@ -1160,6 +1160,9 @@ integer function Stellar_Luminosities_Index_From_Properties(filterName,filterTyp
        write (label,'(f7.4)') redshiftBand
        message=message//' : '//label
     end if
+    if (present(postprocessChain)) then
+       message=message//' : '//postprocessChain
+    end if
     message=message//']'
     do i=1,luminosityCount
        message=message//char(10)//i//" of "//luminosityCount