 PBH}}(t_0)$ is the number density of PBHs of a given mass $M$ today, $z(t)$ is the redshift and the time integral runs from $t{_{\rm min}} = 380\,000\,$years at last scattering of the CMB to $t{_{\rm max}} = {\rm Max}(\tau(M),t_0)$ where $\tau(M)\sim M^3$ is the PBH lifetime and $t_0$ is the age of the Universe. As the Universe is expanding, the number density of PBHs evolves as $(1+z(t))^{-3}$, and the energy of the emitted photons evolves as $(1+z(t))^{-1}$. A last factor $(1+z(t))$ comes from the change of integrand variable from the line of sight to the present time.

\[sec:results\]Results
======================

We have used the new public code `BlackHawk` [@BlackHawk] to compute the HR of [Eq. (\[eq:hawking\])]{} and the PBH evolution given by Eqs.  and. We consider monochromatic PBH distributions of masses comprised between $M{_{\rm min}} = 10^{13}\,$g and $M{_{\rm max}} = 10^{18}\,$g and initial spin parameters between $a_{i{\rm,min}}^* = 0$ and $a_{i{\rm,max}}^* = 0.9999$, and compute the integral of [Eq. (\[eq:flux\_IGRB\])]{} over the redshift (matter-dominated era) $$z(t) = \left( \dfrac{1}{H_0 t} \right)^{2/3} -1\,, \label{eq:redshift}$$ where $H_0$ is the present Hubble parameter. We then compare the result of the integral to the measured IGRB and find the maximum allowed value of the present PBH number density $n{_{\rm PBH}}(t_0)$ at a given PBH mass $M$, with a conservative approach taking into account the most stringent constraints (e.g. FERMI-LAT$_{\rm C}$ at $E = 1\,$GeV). The corresponding limit on the DM fraction $f$ constituted of PBHs of mass $M$ is obtained through $ n{_{\rm PBH}}(t_0) = f{\rho{_{\rm DM}}}/{M}\,
PBH}}(t_0)$ is the number density of PBHs of a given mass $ M$ today, $ z(t)$ is the red shift and the meter integral runs from $ t{_{\rm minute } } = 380\,000\,$years at last scatter of the CMB to $ t{_{\rm max } } = { \rm Max}(\tau(M),t_0)$ where $ \tau(M)\sim M^3 $ is the PBH lifetime and $ t_0 $ is the age of the Universe. As the Universe is inflate, the number density of PBHs evolve as $ (1+z(t))^{-3}$, and the energy of the emitted photons develop as $ (1+z(t))^{-1}$. A last factor $ (1+z(t))$ come from the change of integrand variable star from the line of view to the present clock time. 

 \[sec: results\]Results 
 = = = = = = = = = = = = = = = = = = = = = = 

 We have used the new public code 'BlackHawk '[ @BlackHawk ] to compute the hour of [ Eq.   (\[eq: hawking\ ]) ] { } and the PBH evolution given by Eqs.    and. We view monochromatic PBH distributions of multitude comprised between $ M{_{\rm min } } = 10^{13}\,$g and $ M{_{\rm max } } = 10^{18}\,$g and initial spin parameters between $ a_{i{\rm, min}}^ * = 0 $ and $ a_{i{\rm, max}}^ * = 0.9999 $, and compute the integral of [ Eq.   (\[eq: flux\_IGRB\ ]) ] { } over the red shift (matter - dominated era) $ $ z(t) = \left (\dfrac{1}{H_0 t } \right)^{2/3 } -1\, , \label{eq: redshift}$$ where $ H_0 $ is the present Hubble parameter. We then compare the result of the integral to the measured IGRB and rule the maximum allowed value of the present PBH number density $ n{_{\rm PBH}}(t_0)$ at a given PBH mass $ M$, with a conservative approach taking into account the about rigorous constraints (e.g. FERMI - LAT$_{\rm C}$ at $ E = 1\,$GeV). The correspond limit on the DM fraction $ f$ constituted of PBHs of mass $ M$ is prevail through $ n{_{\rm PBH}}(t_0) = f{\rho{_{\rm DM}}}/{M}\,
 PBJ}}(t_0)$ is the number density of PBHs of a gncen mavs $M$ tkday, $z(t)$ ks the redshift and the time ibtegrql runs from $t{_{\rm min}} = 380\,000\,$hears at past scartermng of the CMB to $t{_{\rm maw}} = {\rm Lax}(\tcu(N),t_0)$ where $\tau(M)\slm M^3$ is the PBH lifetime dna $c_0$ is the age of the Universe. As the Tniversr ls expanding, tre nlmfer svnwity of PBHs evolves as $(1+z(t))^{-3}$, ans the eiergy of the emotted photons evolves as $(1+z(h))^{-1}$. A past factor $(1+z(t))$ comfs from the chagte of integrxnd variable from the mine of sight to the present tioe.

\[sec:xesults\]Resuotw
======================

We vave used tie new public code `BlackHafk` [@BlaclHawk] to compube thx HR of [Eq. (\[eq:hawking\])]{} and tie PBH evolution givgn by Eqs.  dnb. We consider monochrimqtic KBH dhstrkvutkona pf massed ckmprised bstween $M{_{\rm nin}} = 10^{13}\,$g and $M{_{\rm max}} = 10^{18}\,$g qnd initial slin pawaieters between $a_{i{\rm,min}}^* = 0$ and $a_{i{\rm,max}}^* = 0.9999$, dnd compute the integral od [Eq. (\[eq:flux\_IGRB\])]{} over tje redshist (matter-dominated era) $$z(t) = \left( \dfrac{1}{H_0 t} \right)^{2/3} -1\,, \ldbel{es:fedwhlft}$$ dyege $H_0$ is the present Hubble parameter. We then skmkarv the result of tme integral to the mfaxored IGRB and wind tks jaximum allowed vapue of jhe prwsent PBH numner density $n{_{\rm PBH}}(t_0)$ at a guven PBH masf $M$, with a conservacive approack takimg inyo account the most strnngent constraintd (e.g. FERMJ-UAT$_{\rm C}$ at $E = 1\,$GeX). Tme worresponding limit on the DM fractmon $f$ constigutec of PFHs of masd $M$ is obtained through $ n{_{\rm KBH}}(t_0) = x{\rho{_{\rm DM}}}/{M}\,
PBH}}(t_0)$ is the number density of PBHs given $M$ today, is the redshift from min}} = 380\,000\,$years last scattering of CMB to $t{_{\rm max}} = {\rm where $\tau(M)\sim M^3$ is the PBH lifetime and $t_0$ is the age of Universe. As the Universe is expanding, the number density of PBHs evolves as and energy the photons evolves as $(1+z(t))^{-1}$. A last factor $(1+z(t))$ comes from the change of integrand variable from line of sight to the present time. \[sec:results\]Results We have used the public code `BlackHawk` [@BlackHawk] to the of [Eq. and PBH given by Eqs. We consider monochromatic PBH distributions of masses comprised between $M{_{\rm min}} = 10^{13}\,$g and $M{_{\rm max}} = and initial between $a_{i{\rm,min}}^* 0$ $a_{i{\rm,max}}^* 0.9999$, and compute of [Eq. (\[eq:flux\_IGRB\])]{} over the redshift = \left( \dfrac{1}{H_0 t} \right)^{2/3} -1\,, \label{eq:redshift}$$ where is the Hubble parameter. We then compare the of the integral to the measured IGRB and the maximum allowed value of the present PBH number density $n{_{\rm PBH}}(t_0)$ at a given $M$, with a conservative taking into account most constraints FERMI-LAT$_{\rm at $E 1\,$GeV). The corresponding limit on the DM fraction $f$ constituted of of mass $M$ is obtained through $ n{_{\rm PBH}}(t_0) =
 PBH}}(t_0)$ is the number density of PbHs of a giveN mass $m$ toDay, $Z(t)$ Is thE redShift and the timE InteGral runs from $t{_{\rm min}} = 380\,000\,$yearS at laSt SCattERiNg of tHe CMB to $T{_{\Rm MAX}} = {\rm maX}(\tAu(M),T_0)$ wHErE $\tau(M)\Sim m^3$ is the PbH lifetime And $T_0$ iS the age of the uNiVerse. As the uniVerse is expanDinG, the nuMbEr dENsity Of PbHs evOlves aS $(1+Z(t))^{-3}$, and tHe energy oF tHE emittED photonS EVoLves As $(1+z(t))^{-1}$. A last factor $(1+z(t))$ COmES from the change Of inteGrANd VARiaBle From the linE oF sighT To the prESeNT TIme.

\[SEc:results\]ResuLts
======================

We have usED thE new puBlIc cODe `BlacKHawk` [@blACkHAwk] to computE the hR of [Eq. (\[eq:hAwking\])]{} ANd the PBh EvolutiOn giveN by eqs.  And. WE CoNsIdeR mONocHRoMatIC PBh distribUtIoNs of mAsseS COMPrisEd bEtweEn $M{_{\rm Min}} = 10^{13}\,$g and $M{_{\rm max}} = 10^{18}\,$G anD iniTIal Spin pArameTers BeTween $A_{i{\rm,miN}}^* = 0$ and $a_{I{\rM,max}}^* = 0.9999$, and compute tHe inTegral of [EQ. (\[eq:FlUx\_IgRb\])]{} over THe redsHifT (maTter-domInated eRA) $$z(t) = \LeFT( \DFrAc{1}{H_0 t} \right)^{2/3} -1\,, \label{eq:reDsHIFt}$$ Where $H_0$ is The preSEnT HUBble paraMeTer. we thEN CompaRe thE ReSult of thE integRAl To The measUrEd IGRB AnD fiNd tHe maxIMum aLlowed Value of tHe preSEnt PBH number deNSity $n{_{\rm PBH}}(t_0)$ at A GiVEN Pbh masS $M$, wIth a conservAtivE ApprOach TAkIng INto acCount ThE MoST stringent constrainTs (E.g. FERMi-LAT$_{\rM C}$ at $E = 1\,$GeV). The coRrespondinG LIMit on the dM frACtIOn $f$ constituted Of PBHS of mass $M$ is OBtained tHrougH $ n{_{\rm PBH}}(t_0) = F{\rho{_{\rm DM}}}/{M}\,
 PBH}}(t_0)$ is the number  density o f PBH s o f a  g iven  mas s $M$ today, $ z (t)$  is the redshift and t he ti me inte g ra l run s from$ t{ _ { \rm  m in }}=3 80 \,000 \,$ years a t last sca tte ri ng of the CM B  t o $t{_{\rm  ma x}} = {\rm M ax} (\tau( M) ,t_ 0 )$ wh ere  $\ta u(M)\s i m M^3$  is the P BH lifeti m e and $ t _ 0$  isthe age of the Un i ve r se. As the Uni verseis ex p a ndi ng,  the numbe rdensi t y of PB H se v o lve s  as $(1+z(t)) ^{-3}$, and the  energ yoft he emi ttedph o ton s evolves a s $( 1+z(t))^{ -1}$.A  last f a ctor $( 1+z(t) )$com es f r om  t hech a nge of  in t egr and vari ab le  from  the l i n e of  si ghtto th e present tim e.
\[s e c:r esult s\]Re sult s===== ====== ===== == ====

We have u sedthe new p ubl ic  co de  `Bla c kHawk`  [@ Bla ckHawk]  to com p ute  t h e HR  of [Eq. (\[eq:haw ki n g \] )]{} and  the P B Hev o lution g iv enby E q s .  an d. W e  c onsidermonoch r om at ic PBHdi stribu ti ons  of  mass e s co mprise d betwee n $M{ _ {\rm min}} = 1 0 ^{13}\,$g and $M { _ {\ r m ma x}}  = 10^{18}\ ,$ga nd i niti a lspi n  para meter sb et w een $a_{i{\rm,min}} ^*  = 0$and $ a_{i{\rm,max} }^* = 0.99 9 9 $ , and co mput e  t h e integral of[Eq.(\[eq:flux \ _IGRB\]) ]{} o ver theredshift( m atter-do min ate d e ra) $ $z (t) = \left(\ d frac {1 }{H_0 t } \ right)^ {2/ 3}-1\ ,,\l abel{eq:r edshift} $$  w he re  $H _0$ i s  the pre se ntHu bbl e par a meter.  We t henco mp a rethe res u lt o f th ein tegr alto  themeas u red  IGRB a nd find t hem axim um  a llowedvalue of thepr esent PBHnu mbe r dens i t y $n{_{\ rm PBH}}(t_0)$ at a giv e n PBH m ass  $M$,  wit h a conse rva tive a ppr o ach ta king i nto a cc oun t the m o s tstr in gent const r a int s (e. g.  FER MI-LAT$ _{\rm C}$ at $E =1 \,$ GeV). The cor res pond i n glim i to n t he DMf r action $f$ cons tituted of  P B Hs  of mass $ M $ i sobtaine d throu gh $n {_{\rmPBH}}(t_0 ) = f{\rh o{ _{\r m DM} }}/{M}\,
 PBH}}(t_0)$_is the_number density of PBHs_of a_given_mass $M$_today,_$z(t)$ is the_redshift and the_time integral runs from_$t{_{\rm min}} =_380\,000\,$years_at last scattering of the CMB to $t{_{\rm max}} = {\rm Max}(\tau(M),t_0)$ where $\tau(M)\sim_M^3$_is the_PBH_lifetime_and $t_0$ is the age_of the Universe. As the_Universe is_expanding, the number density of PBHs evolves as_$(1+z(t))^{-3}$,_and the energy_of the emitted photons evolves as $(1+z(t))^{-1}$. A last_factor $(1+z(t))$ comes from the change_of integrand variable_from_the_line of sight to_the present time.

\[sec:results\]Results
======================

We have used the_new public code `BlackHawk` [@BlackHawk] to_compute the HR of [Eq. (\[eq:hawking\])]{} and the_PBH evolution given by Eqs.  and._We consider monochromatic PBH distributions_of masses_comprised between $M{_{\rm min}} =_10^{13}\,$g and $M{_{\rm_max}} =_10^{18}\,$g and initial_spin parameters between $a_{i{\rm,min}}^* = 0$_and $a_{i{\rm,max}}^* =_0.9999$, and compute the integral of_[Eq. (\[eq:flux\_IGRB\])]{}_over the redshift_(matter-dominated_era)_$$z(t) =_\left( \dfrac{1}{H_0 t}_\right)^{2/3}_-1\,, \label{eq:redshift}$$_where_$H_0$ is the present Hubble parameter._We_then compare the result of the integral_to the measured IGRB_and_find the maximum allowed_value of the present PBH_number density $n{_{\rm PBH}}(t_0)$ at a_given PBH_mass $M$,_with a conservative approach taking into account the most stringent constraints_(e.g. FERMI-LAT$_{\rm C}$ at $E =_1\,$GeV). The corresponding limit_on the_DM_fraction $f$ constituted_of_PBHs of_mass $M$ is obtained through $ n{_{\rm_PBH}}(t_0) =_f{\rho{_{\rm DM}}}/{M}\,