The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada-French-Hawaii telescope (CFHT). We use a new set of data extracted from VESTIGE to study the impact of the hostile cluster environment on the star formation process down to the scale of HII regions (~ 50 pc). HII regions are identified and their parameters measured using the HIIphot code on a sample of 114 late-type galaxies spanning a wide range in morphological type (Sa-Sd, Im, BCD), stellar mass (10^6.5 <= M_star <= 10^11 Mo), and star formation activity (10^-3 <= SFR <= 10 Mo yr^-1). Owing to the exquisite average resolution of the VESTIGE data (0.65 arcsec), we detect 11302 HII regions with an Halpha luminosity L(Halpha) >= 10^37 erg s^-1. We show that the typical number of HII regions in gas-stripped objects is significantly lower than in healthy late-types of similar stellar mass. We also show that in these gas-stripped galaxies the number of HII regions significantly drops outside the effective radius, suggesting that the quenching process occurs outside-in, in agreement with other multifrequency observations. These new results consistently confirm that the main mechanism responsible for the decrease of the star formation activity observed in cluster galaxies is ram pressure, allowing us to discard other milder processes such as starvation or strangulation unable to reproduce the observed radially truncated profiles.

We present an optical to near-infrared selected astronomical catalogue covering 1270 deg.$^2$. This is the first attempt to systematically combine data from 23 of the premier extragalactic survey fields - the product of a vast investment of telescope time. The fields are those imaged by the Herschel Space Observatory which form the Herschel Extragalactic Legacy Project (HELP). Our catalogue of 170 million objects is constructed by a positional cross match of 51 public surveys. This high resolution optical, near-infrared, and mid-infrared catalogue is designed for photometric redshift estimation, extraction of fluxes in lower resolution far-infrared maps, and spectral energy distribution modelling. It collates, standardises, and provides value added derived quantities including corrected aperture magnitudes and astrometry correction over the Herschel extragalactic wide fields for the first time. $grizy$ fluxes are available on all fields with $g$ band data reaching $5\sigma$ point-source depths in a 2 arcsec aperture of 23.5, 24.4, and 24.6 (AB) mag at the 25th, 50th, and 75th percentiles, by area covered, across all HELP fields. It has $K$ or $K_s$ coverage over 1146 deg.$^2$ with depth percentiles of 20.2, 20.4, and 21.0 mag respectively. The IRAC Ch 1 band is available over 273 deg.$^2$ with depth percentiles of 17.7, 21.4, and 22.2 mag respectively. This paper defines the "masterlist" objects for the first data release (DR1) of HELP. This large sample of standardised total and corrected aperture fluxes, uniform quality flags, and completeness measures provides large well understood statistical samples over the full Herschel extragalactic sky.

The Herschel Extragalactic Legacy Project (HELP) focuses to publish an astronomical multiwavelength catalogue of millions of objects over 1300~deg$^2$ of the Herschel Space Observatory survey fields. Millions of galaxies with ultraviolet--far infrared photometry {make} HELP a perfect sample for testing spectral energy distribution fitting models, and to prepare tools for next-generation data. In the frame of HELP collaboration we estimated the main physical properties of all galaxies from the HELP database and we checked a new procedure to select peculiar galaxies from large galaxy sample and we investigated the influence of used modules for stellar mass estimation.

Context. Measuring how the physical properties of galaxies change across cosmic times is essential to understand galaxy formation and evolution. With the advent of numerous ground-based and space-borne instruments launched over the past few decades we now have exquisite multi-wavelength observations of galaxies from the FUV to the radio domain. To tap into this mine of data and obtain new insight into the formation and evolution of galaxies, it is essential that we are able to extract information from their SED. Aims. We present a completely new implementation of CIGALE. Written in python, its main aims are to easily and efficiently model the FUV to radio spectrum of galaxies and estimate their physical properties such as star formation rate, attenuation, dust luminosity, stellar mass, and many other physical quantities. Methods. To compute the spectral models, CIGALE builds composite stellar populations from simple stellar populations combined with highly flexible star formation histories, calculates the emission from gas ionised by massive stars, and attenuates both the stars and the ionised gas with a highly flexible attenuation curve. Based on an energy balance principle, the absorbed energy is then re-emitted by the dust in the mid- and far-infrared domains while thermal and non-thermal components are also included, extending the spectrum far into the radio range. A large grid of models is then fitted to the data and the physical properties are estimated through the analysis of the likelihood distribution. Results. CIGALE is a versatile and easy-to-use tool that makes full use of the architecture of multi-core computers, building grids of millions of models and analysing samples of thousands of galaxies, both at high speed. Beyond fitting the SEDs of galaxies and parameter estimations, it can also be used as a model-generation tool or serve as a library to build new applications.

We have observed the late-type peculiar galaxy NGC 4424 during VESTIGE, a blind narrow-band Halpha[NII] imaging survey of the Virgo cluster carried out with MegaCam at the CFHT. The presence of a 110 kpc long HI tail in the S direction indicates that this galaxy is undergoing a ram pressure stripping event. The deep narrow-band image revealed the presence of a low surface brightness ionised gas tail ~10 kpc long extending from the centre of the galaxy to the NW direction, in the direction opposite to the HI tail. Chandra and XMM X-rays data do not show any compact source in the nucleus nor the presence of an extended tail of hot gas, while MUSE spectroscopy indicates that the gas is photo-ionised in the inner regions and shock-ionised in the outer parts. IFU spectroscopy confirms that the ionised gas is kinematically decoupled from the stellar component and indicates the presence of two kinematically distinct structures in the stellar disc. The analysis of the SED of the galaxy indicates that the activity of star formation has been totally quenched in the outer disc ~ 250-280 Myr ago, while only reduced by ~80% in the central regions. All this observational evidence suggests that NGC 4424 is the remnant of an unequal-mass merger occurred <= 500 Myr ago, when the galaxy was already a member of the Virgo cluster, now undergoing a ram pressure stripping event which has removed the gas and quenched the activity of star formation in the outer disc. The tail of ionised gas probably results from the outflow produced by a central starburst fed by the collapse of gas induced by the merging episode. This outflow is sufficiently powerful to overcome the ram pressure induced by the intracluster medium on the disc of the galaxy crossing the cluster. This analysis thus suggests that feedback can participate in the quenching process of galaxies in high-density regions.

The HELP project focuses on the data from ESA's Herschel mission, which covered over 1300$deg^2$ and is preparing to publish a multi-wavelength catalogue of millions of objects. Our main goal is to find the best approach to simultaneously fitting SEDs of millions of galaxies across a wide redshift range to obtain homogeneous estimates of the main physical parameters of detected IR galaxies. We perform SED fitting on the UV/NIR to FIR emission of 42 047 galaxies from the pilot HELP field: ELAIS N1. We use the latest release of CIGALE, a galaxy SED fitting code relying on energy balance, to deliver the main physical parameters such as $M_{star}$, SFR, and $L_{dust}$. We implement additional quality criteria to the fits by calculating $\chi^2$ values for the stellar and dust part of the spectra independently. These criteria allow us to identify the best fits and to identify peculiar galaxies. We perform the SED fitting by assuming three different dust attenuation laws separately allowing us to test the impact of the assumed law on estimated physical parameters. We implemented two additional quality value checks for the SED fitting method based on $M_{star}$ estimation and energy budget. This method allows us to identify possible objects with incorrect matching in the catalogue and peculiar galaxies; we found 351 possible candidates of lensed galaxies using two complementary $\chi^2$s criteria (stellar and IR) and $z_{phot}$ calculated for the IR part of the spectrum only. We find that the attenuation law has an important impact on the $M_{star}$ estimate (on average leading to disparities of a factor of two). We derive the relation between $M_{star}$ estimates obtained by different attenuation laws and we find the best recipe for our sample. We also make independent estimates of the total $L_{dust}$ parameter from stellar emission by fitting the galaxies with and without IR data separately.

Dust attenuation shapes the spectral energy distribution of galaxies. It is particularly true for dusty galaxies in which stars experience a heavy attenuation. The combination of UV-to-IR photometry with the spectroscopic measurement of the H$\alpha$ recombination line helps to quantify dust attenuation of the whole stellar population and its wavelength dependence. We selected an IR complete sample of galaxies in the COSMOS 3D-HST CANDELS field detected with the Herschel satellite with a signal to noise ratio larger than five. Optical to NIR photometry is available as well as NIR spectroscopy for each source. We reduced the sample to the redshift range $0.6 < z < 1.6$ to include the H$\alpha$ line in the G141 grism spectra. We have used a new version of the CIGALE code to fit simultaneously the continuum and H$\alpha$ line emission of the 34 selected galaxies. Using flexible attenuation laws with free parameters, we are able to measure the shape of the attenuation curve for each galaxy as well as the amount of attenuation of each stellar population, the former being in general steeper than the starburst law in the UV-optical with a large variation of the slope among galaxies. The attenuation of young stars or nebular continuum is found on average about twice the attenuation affecting older stars, again with a large variation. Our model with power-laws, based on a modification of the Charlot and Fall recipe, gives results in better agreement with the radiative transfer models than the global modification of the slope of the Calzetti law.

During pilot observations of the Virgo Environmental Survey Tracing Galaxy Evolution (VESTIGE), a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the CFHT, we have observed the spiral galaxy NGC 4254 (M99). Deep Halpha+[NII] narrow-band and GALEX UV images revealed the presence of 60 compact (70-500 pc radius) star forming regions up to ~ 20 kpc outside the optical disc of the galaxy. These regions are located along a tail of HI gas stripped from the disc of the galaxy after a rapid gravitational encounter with another Virgo cluster member that simulations indicate occurred 280-750 Myr ago. We have combined the VESTIGE data with multifrequency data from the UV to the far-infrared to characterise the stellar populations of these regions and study the star formation process in an extreme environment such as the tails of stripped gas embedded in the hot intracluster medium. The colour, spectral energy distribution (SED), and linear size consistently indicate that these regions are coeval and have been formed after a single burst of star formation that occurred ~< 100 Myr ago. These regions might become free floating objects within the cluster potential well, and be the local analogues of compact sources produced after the interaction of gas-rich systems that occurred during the early formation of clusters.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey carried out with MegaCam at the Canada-France-Hawaii Telescope. The survey covers the whole Virgo cluster region from its core to one virial radius (104 deg^2). The sensitivity of the survey is of f(Halpha) ~ 4 x 10^-17 erg sec-1 cm^-2 (5 sigma detection limit) for point sources and Sigma (Halpha) ~ 2 x 10^-18 erg sec^-1 cm^-2 arcsec^-2 (1 sigma detection limit at 3 arcsec resolution) for extended sources, making VESTIGE the deepest and largest blind narrow-band survey of a nearby cluster. This paper presents the survey in all its technical aspects, including the survey design, the observing strategy, the achieved sensitivity in both the narrow-band Halpha+[NII] and in the broad-band r filter used for the stellar continuum subtraction, the data reduction, calibration, and products, as well as its status after the first observing semester. We briefly describe the Halpha properties of galaxies located in a 4x1 deg^2 strip in the core of the cluster north of M87, where several extended tails of ionised gas are detected. This paper also lists the main scientific motivations of VESTIGE, which include the study of the effects of the environment on galaxy evolution, the fate of the stripped gas in cluster objects, the star formation process in nearby galaxies of different type and stellar mass, the determination of the Halpha luminosity function and of the Halpha scaling relations down to ~ 10^6 Mo stellar mass objects, and the reconstruction of the dynamical structure of the Virgo cluster. This unique set of data will also be used to study the HII luminosity function in hundreds of galaxies, the diffuse Halpha+[NII] emission of the Milky Way at high Galactic latitude, and the properties of emission line galaxies at high redshift.

In this work we investigate the far-UV to NIR shape of the dust attenuation curve of a sample of IR selected dust obscured (U)LIRGs at z$\sim$2. The spectral energy distributions (SEDs) are fitted with CIGALE, a physically-motivated spectral synthesis model based on energy balance. Its flexibility allows us to test a wide range of different analytical prescriptions for the dust attenuation curve, including the well-known Calzetti and Charlot & Fall curves, and modified versions of them. The attenuation curves computed under the assumption of our reference double power-law model are in very good agreement with those derived, in previous works, with radiative transfer (RT) SED fitting. We investigate the position of our galaxies in the IRX-$\beta$ diagram and find this to be consistent with grayer slopes, on average, in the UV. We also find evidence for a flattening of the attenuation curve in the NIR with respect to more classical Calzetti-like recipes. This larger NIR attenuation yields larger derived stellar masses from SED fitting, by a median factor of $\sim$ 1.4 and up to a factor $\sim$10 for the most extreme cases. The star formation rate appears instead to be more dependent on the total amount of attenuation in the galaxy. Our analysis highlights the need for a flexible attenuation curve when reproducing the physical properties of a large variety of objects.

We present XID+ a new generation of software for prior-based photometry extraction in the Herschel SPIRE maps. Based on a Bayesian framework, XID+ allows the inclusion of prior information and gives access to the full posterior probability distribution of fluxes. XID+ is developed within the Herschel Extragalactic Legacy Project (HELP) and is available at https://github.com/H-E-L-P/XID_plus.

Low surface brightness galaxies (LSBGs) represent a significant percentage of local galaxies but their formation and evolution remain elusive. They may hold crucial information for our understanding of many key issues (i.e., census of baryonic and dark matter, star formation in the low density regime, mass function). The most massive examples - the so called giant LSBGs - can be as massive as the Milky Way, but with this mass being distributed in a much larger disk. Malin 1 is an iconic giant LSBG, perhaps the largest disk galaxy known. We attempt to bring new insights on its structure and evolution on the basis of new images covering a wide range in wavelength. We have computed surface brightness profiles (and average surface brightnesses in 16 regions of interest), in six photometric bands (FUV, NUV, u, g, i, z). We compared these data to various models, testing a variety of assumptions concerning the formation and evolution of Malin 1. We find that the surface brightness and color profiles can be reproduced by a long and quiet star-formation history due to the low surface density; no significant event, such as a collision, is necessary. Such quiet star formation across the giant disk is obtained in a disk model calibrated for the Milky Way, but with an angular momentum approximately 20 times larger. Signs of small variations of the star-formation history are indicated by the diversity of ages found when different regions within the galaxy are intercompared.For the first time, panchromatic images of Malin 1 are used to constrain the stellar populations and the history of this iconic example among giant LSBGs. Based on our model, the extreme disk of Malin 1 is found to have a long history of relatively low star formation (about 2 Msun/yr). Our model allows us to make predictions on its stellar mass and metallicity.

We study the star formation quenching mechanism in cluster galaxies by fitting the SED of the Herschel Reference Survey, a complete volume-limited K-band-selected sample of nearby galaxies including objects in different density regions, from the core of the Virgo cluster to the general field. The SED are fitted using the CIGALE SED modelling code. The truncated activity of cluster galaxies is parametrised using a specific SFH with 2 free parameters, the quenching age QA and the quenching factor QF. These 2 parameters are crucial for the identification of the quenching mechanism which acts on long timescales if starvation while rapid and efficient if ram pressure. To be sensitive to an abrupt and recent variation of the star formation activity, we combine in a new way 20 UV to FIR photometric bands with 3 age-sensitive Balmer line absorption indices extracted from available medium-resolution integrated spectroscopy and with Halpha narrow band imaging data. The use of a truncated SFH significantly increases the quality of the fit in those objects whose atomic gas content has been removed during the interaction with the hostile cluster environment. The typical QA of the perturbed late-type galaxies is QA < 300 Myr whenever the activity of star formation is reduced by 50% < QF <= 80% and QA < 500 Myr for QF > 80%, while that of the quiescent early-types is QA ~ 1-3 Gyr. The fraction of late-types with a star formation activity reduced by QF > 80% and with an HI-deficiency parameter HI-def > 0.4 drops by a factor of ~ 5 from the inner half virial radius of the Virgo cluster, where the hot diffuse X-ray emitting gas of the cluster is located, to the outer regions. The efficient quenching of the star formation activity observed in Virgo suggests that the dominant stripping process is ram pressure. We discuss the implication of this result in the cosmological context of galaxy evolution.

Gas flows play a fundamental role in galaxy formation and evolution, providing the fuel for the star formation process. These mechanisms leave an imprint in the amount of heavy elements. Thus, the analysis of this metallicity signature provides additional constraint on the galaxy formation scenario. We aim to discriminate between four different galaxy formation models based on two accretion scenarios and two different star formation recipes. We address the impact of a bimodal accretion scenario and a strongly regulated star formation recipe. We present a new extension of the eGalICS model, which allows us to track the metal enrichment process. Our new chemodynamical model is applicable for situations ranging from metal-free primordial accretion to very enriched interstellar gas contents. We use this new tool to predict the metallicity evolution of both the stellar populations and gas phase. We also address the evolution of the gas metallicity with the star formation rate (SFR). We then focus on a sub-sample of Milky Way-like galaxies. We compare both the cosmic stellar mass assembly and the metal enrichment process of such galaxies with observations and detailed chemical evolution models. Our models, based on a strong star formation regulation, allow us to reproduce well the stellar mass to gas-phase metallicity relation observed in the local universe. However, we observe a systematic shift towards high masses. Our $Mstar-Zg-SFR relation is in good agreement with recent measurements: our best model predicts a clear dependence with the SFR. Both SFR and metal enrichment histories of our Milky Way-like galaxies are consistent with observational measurements and detailed chemical evolution models. We finally show that Milky Way progenitors start their evolution below the observed main sequence and progressively reach this observed relation at z = 0.

We use the GALEX data of the GUViCS survey to construct the NUV luminosity function of the Virgo cluster over ~ 300 deg.2, an area covering the cluster and its surrounding regions up to ~ 1.8 virial radii. The NUV luminosity function is also determined for galaxies of different morphological type and NUV-i colour, and for the different substructures within the cluster. These luminosity functions are robust vs. statistical corrections since based on a sample of 833 galaxies mainly identified as cluster members with spectroscopic redshift (808) or high-quality optical scaling relations (10). We fit these luminosity functions with a Schechter function, and compare the fitted parameters with those determined for other nearby clusters and for the field. The faint end slope of the Virgo NUV luminosity function (alpha = -1.19), here sampled down to ~ NUV = -11.5 mag, is significantly flatter than the one measured in other nearby clusters and similar to the field one. Similarly M* = -17.56 is one-to-two magnitudes fainter than measured in Coma, A1367, the Shapley supercluster, and the field. These differences seem due to the quite uncertain statistical corrections and the small range in absolute magnitude sampled in these clusters. We do not observe strong systematic differences in the overall NUV luminosity function of the core of the cluster with respect to that of its periphery. We notice, however, that the relative contribution of red-to-blue galaxies at the faint end is inverted, with red quiescent objects dominating the core of the cluster and star forming galaxies dominating beyond one virial radius. These observational evidences are discussed in the framework of galaxy evolution in dense environments.

We study the role of the environment on galaxy evolution using a sample of 868 galaxies in the Virgo cluster and in its surrounding regions selected from the GUViCS Survey with the purpose of understanding the origin of the red sequence in dense environments. We collected multifrequency data covering the whole electromagnetic spectrum for most of the galaxies. We identify the different dynamical substructures composing the Virgo cluster and we calculate the local density of galaxies using different methods. We then study the distribution of galaxies belonging to the red sequence, the green valley, and the blue cloud within the different cluster substructures. Our analysis indicates that all the most massive galaxies are slow rotators and are the dominant galaxies of the different cluster substructures generally associated with a diffuse X-ray emission. They are probably the result of major merging events that occurred at early epochs. Slow rotators of lower stellar mass are also preferentially located within the different high-density substructures of the cluster. They are virialised within the cluster, thus Virgo members since its formation. They have been shaped by gravitational perturbations occurring within the infalling groups that later formed the cluster. On the contrary, low-mass star-forming systems are extremely rare in the inner regions of the Virgo cluster A, where the density of the intergalactic medium is at its maximum. Our ram pressure stripping models consistently indicate that these star-forming systems can be rapidly deprived of their interstellar medium during their interaction with the intergalactic medium. The lack of gas quenches their star formation activity transforming them into quiescent dwarf ellipticals. This mild transformation does not perturb the kinematic properties of these galaxies which still have rotation curves typical of star-forming systems.

The Herschel Multi-tiered Extragalactic Survey (HerMES) is the largest Guaranteed Time Key Programme on the Herschel Space Observatory. With a wedding cake survey strategy, it consists of nested fields with varying depth and area totalling ~380 deg^2. In this paper, we present deep point source catalogues extracted from Herschel-SPIRE observations of all HerMES fields, except for the later addition of the 270 deg^2 HeLMS field. These catalogues constitute the second Data Release (DR2) made in October 2013. A subset of these catalogues, which consists of bright sources extracted from Herschel-SPIRE observations completed by May 1, 2010 (covering ~ 74 deg^2) were released earlier in the first extensive Data Release (DR1) in March 2012. Two different methods are used to generate the point source catalogues, the SUSSEXtractor (SXT) point source extractor used in two earlier data releases (EDR and EDR2) and a new source detection and photometry method. The latter combines an iterative source detection algorithm, StarFinder (SF), and a De-blended SPIRE Photometry (DESPHOT) algorithm. We use end-to-end Herschel-SPIRE simulations with realistic number counts and clustering properties to characterise basic properties of the point source catalogues, such as the completeness, reliability, photometric and positional accuracy. Over 500, 000 catalogue entries in HerMES fields (except HeLMS) are released to the public through the HeDAM website (http://hedam.oamp.fr/herMES).

We study the link between observed ultraviolet luminosity, stellar mass, and dust attenuation within rest-frame UV-selected samples at z~ 4, 3, and 1.5. We measure by stacking at 250, 350, and 500 um in the Herschel/SPIRE images from the HerMES program the average infrared luminosity as a function of stellar mass and UV luminosity. We find that dust attenuation is mostly correlated with stellar mass. There is also a secondary dependence with UV luminosity: at a given UV luminosity, dust attenuation increases with stellar mass, while at a given stellar mass it decreases with UV luminosity. We provide new empirical recipes to correct for dust attenuation given the observed UV luminosity and the stellar mass. Our results also enable us to put new constraints on the average relation between star formation rate and stellar mass at z~ 4, 3, and 1.5. The star formation rate-stellar mass relations are well described by power laws (SFR~ M^0.7), with the amplitudes being similar at z~4 and z~3, and decreasing by a factor of 4 at z~1.5 at a given stellar mass. We further investigate the evolution with redshift of the specific star formation rate. Our results are in the upper range of previous measurements, in particular at z~3, and are consistent with a plateau at 3<z<4. Current model predictions (either analytic, semi-analytic or hydrodynamic) are inconsistent with these values, as they yield lower predictions than the observations in the redshift range we explore. We use these results to discuss the star formation histories of galaxies in the framework of the Main Sequence of star-forming galaxies. Our results suggest that galaxies at high redshift (2.5<z<4) stay around 1 Gyr on the Main Sequence. With decreasing redshift, this time increases such that z=1 Main Sequence galaxies with 10^8<M_*/Msun<10^10 stay on the Main Sequence until z=0.

We present pcigale, the port to Python of CIGALE (Code Investigating Galaxy Emission) a Fortran spectral energy distribution (SED) fitting code developed at the Laboratoire d'Astrophysique de Marseille. After recalling the specifics of the SED fitting method, we show the gains in modularity and versatility offered by Python, as well as the drawbacks compared to the compiled code.

We study the far-infrared (IR) and sub-millimeter properties of a sample of ultraviolet (UV) selected galaxies at z\sim1.5. Using stacking at 250, 350 and 500 um from Herschel Space Observatory SPIRE imaging of the COSMOS field obtained within the HerMES key program, we derive the mean IR luminosity as a function of both UV luminosity and slope of the UV continuum beta. The IR to UV luminosity ratio is roughly constant over most of the UV luminosity range we explore. We also find that the IR to UV luminosity ratio is correlated with beta. We observe a correlation that underestimates the correlation derived from low-redshift starburst galaxies, but is in good agreement with the correlation derived from local normal star-forming galaxies. Using these results we reconstruct the IR luminosity function of our UV-selected sample. This luminosity function recovers the IR luminosity functions measured from IR selected samples at the faintest luminosities (Lir ~ 10^{11} L_sun), but might underestimate them at the bright-end (Lir > 5.10^{11} L_sun). For galaxies with 10^{11}<Lir/L_sun<10^{13}, the IR luminosity function of a UV selection recovers (given the differences in IR-based estimates) 52-65 to 89-112 per cent of the star-formation rate density derived from an IR selection. The cosmic star-formation rate density derived from this IR luminosity function is 61-76 to 100-133 per cent of the density derived from IR selections at the same epoch. Assuming the latest Herschel results and conservative stacking measurements, we use a toy model to fully reproduce the far IR luminosity function from our UV selection at z\sim 1.5. This suggests that a sample around 4 magnitudes deeper (i.e. reaching u \sim 30 mag) and a large dispersion of the IR to UV luminosity ratio are required.

We present CIGALE (Burgarella et al. 2005; Noll et al. 2009), a software developed at the Laboratoire d'Astrophysique de Marseille to fit galaxy spectral energy distributions from the rest-frame far-UV to far-IR wavelength range, and to derive some of their physical parameters. We also give some examples of scientific results obtained with CIGALE.