_q}\over{\beta(1-\beta)Q^2}}),$$ and $$\begin{aligned}
\frac{1}{z}=1+\frac{\kappa_{t}^{2}+m_q^2}{(1-\beta)Q^2}+\frac{k_t^2+\kappa_t^2-2{\bf{\kappa_t}}.{\bf{k_t}}+m_q^2}{\beta
Q^2}
 \label{eq:a}.\end{aligned}$$ As in the reference [@4kimber], the scale $\mu$ which controls the $unintegrated$ gluon and the $QCD$ coupling constant $\alpha_s$, is chosen as follows: $$\begin{aligned}
\mu^2=k_t^2+\kappa_t^2+m_q^2
 \label{eq:b}.\end{aligned}$$ One should note that the coefficients used for quark and non-perturbative gluon contributions depend on the transverse momentum. As it has been briefly explained before, the main prescription for $F_L$ consists of three terms; the first term is the $k_t$ factorization which explains the contribution of the $UPDF$ into the $F_L$. This term is derived with the use of pure gluon contribution. However, it only counts the gluon contributions coming from the perturbative region, i.e. for $k_t>1$ $GeV$, and does not have anything to do with the non-perturbative contributions. In the reference [@stasto4], it has been shown that a proper non-perturbative term can be derived from the $k_t$ factorization term, compacting the $k_t$ dependence and the integration with the use of a variable-change, i.e. y, that carries the $k_t$ dependence. Nevertheless, there is a calculable quark contribution in the longitudinal structure function of the proton, which comes from the collinear factorization, i.e. the second term of the equation (\[eq:2p\]).

For the charm quark, $m$ is taken to be $m_c=1.4 GeV$, and $u$, $d$ and $s$ quarks masses are neglected. We also use the same approximation to save the computation time [@tkimber], the one we did for the calculation of $F_2(x,Q^2)$ [@mho1] i.e
_ q}\over{\beta(1-\beta)Q^2}}),$$ and $ $ \begin{aligned } 
 \frac{1}{z}=1+\frac{\kappa_{t}^{2}+m_q^2}{(1-\beta)Q^2}+\frac{k_t^2+\kappa_t^2 - 2{\bf{\kappa_t}}.{\bf{k_t}}+m_q^2}{\beta 
 Q^2 } 
  \label{eq: a}.\end{aligned}$$ As in the reference [ @4kimber ], the scale $ \mu$ which controls the $ unintegrated$ gluon and the $ QCD$ coupling changeless $ \alpha_s$, is choose as follows: $ $ \begin{aligned } 
 \mu^2 = k_t^2+\kappa_t^2+m_q^2 
  \label{eq: b}.\end{aligned}$$ One should note that the coefficients use for quark and non - perturbative gluon contributions depend on the cross momentum. As it has been briefly explained before, the main prescription drug for $ F_L$ consist of three terms; the beginning term is the $ k_t$ factorization which explains the contribution of the $ UPDF$ into the $ F_L$. This terminus is derive with the use of arrant gluon contribution. However, it only counts the gluon contribution coming from the perturbative region, i.e. for $ k_t>1 $ $ GeV$, and does not have anything to practice with the non - perturbative contributions. In the reference [ @stasto4 ], it has been shown that a proper non - perturbative term can be deduce from the $ k_t$ factorization term, compacting the $ k_t$ dependence and the integration with the use of a variable star - change, i.e. y, that carries the $ k_t$ dependence. Nevertheless, there is a calculable quark contribution in the longitudinal social organization affair of the proton, which comes from the collinear factorization, i.e. the second term of the equation (\[eq:2p\ ]). 

 For the charm quark, $ m$ is taken to be $ m_c=1.4 GeV$, and $ u$, $ d$ and $ s$ quarks masses are neglected. We besides practice the same approximation to save the calculation time [ @tkimber ], the one we did for the calculation of $ F_2(x, Q^2)$ [ @mho1 ] i.e
_q}\ovfr{\beta(1-\beta)Q^2}}),$$ and $$\begin{alinned}
\frac{1}{z}=1+\frac{\kapka_{r}^{2}+m_q^2}{(1-\bete)Q^2}+\frac{k_f^2+\kappa_t^2-2{\bw{\kappa_t}}.{\bf{k_t}}+m_q^2}{\beta
Q^2}
 \label{eq:a}.\eid{alugned}$$ As in the reference [@4kkmber], the scale $\my$ whmch controls the $unintegvcted$ fpuon end the $QCD$ coukling constatt $\alpha_s$, is cvoreu as follows: $$\begin{aligned}
\mu^2=k_t^2+\kappa_t^2+m_z^2
 \label{rq:h}.\end{aligned}$$ Ong shollq nofv uhat the coefficients used for quzrk and non-perturbatove gluon contributions deoend on the transverse momentum. Aw it yas been bridfly explapued before, jhe main prescription for $F_L$ conrists of three jzems; jhe first tecm is nhe $k_t$ factorldation fhich ecplains the coktribntiob of the $UPDF$ into thx $F_L$. This term is detived with tke use of pure gluon xobtribotion. Howdcer, it oily countd tie gluon cohtributions coming from the petttgnative regioh, i.e. fjr $k_t>1$ $GeV$, and does not have anything to dm wjth the non-perturbative contributions. In the referencq [@stasto4], it has been shown that a proper non-pertusbatite texn can ve derived from the $k_t$ factorization term, compasfimg the $k_t$ dependcnce and the integtahipg with the usg of a rzrjable-change, i.e. y, tjat carties tye $k_t$ deptndenve. Nevertheless, there is a xalculable qlark contribution in tke longitudiual sttucturr function of the protou, whicg comes frol the colmknear factorizatkon, i.a. the second term of the ezuation (\[ew:2p\]).

Fox the chxrm auark, $i$ is taken to bc $m_c=1.4 GeV$, and $u$, $d$ anf $s$ qoarks kasses are neglected. We also use the same approximation tm sdve the eomputstion time [@tkymber], the one ce did fjr thd calculatpon of $F_2(x,X^2)$ [@mho1] i.e
_q}\over{\beta(1-\beta)Q^2}}),$$ and $$\begin{aligned} \frac{1}{z}=1+\frac{\kappa_{t}^{2}+m_q^2}{(1-\beta)Q^2}+\frac{k_t^2+\kappa_t^2-2{\bf{\kappa_t}}.{\bf{k_t}}+m_q^2}{\beta Q^2} \label{eq:a}.\end{aligned}$$ As reference the scale which controls the coupling $\alpha_s$, is chosen follows: $$\begin{aligned} \mu^2=k_t^2+\kappa_t^2+m_q^2 One should note that the coefficients for quark and non-perturbative gluon contributions depend on the transverse momentum. As it been briefly explained before, the main prescription for $F_L$ consists of three terms; first is $k_t$ which explains the contribution of the $UPDF$ into the $F_L$. This term is derived with the of pure gluon contribution. However, it only counts gluon contributions coming from perturbative region, i.e. for $k_t>1$ and not have to with non-perturbative contributions. In reference [@stasto4], it has been shown that a proper non-perturbative term can be derived from the $k_t$ term, compacting dependence and integration the of a variable-change, that carries the $k_t$ dependence. Nevertheless, calculable quark contribution in the longitudinal structure function the proton, comes from the collinear factorization, i.e. second term of the equation (\[eq:2p\]). For the quark, $m$ is taken to be $m_c=1.4 GeV$, and $u$, $d$ and $s$ quarks masses We also use the approximation to save computation [@tkimber], one did for calculation of $F_2(x,Q^2)$ [@mho1] i.e
_q}\over{\beta(1-\beta)Q^2}}),$$ and $$\begin{aliGned}
\frac{1}{z}=1+\fRac{\kaPpa_{T}^{2}+m_q^2}{(1-\BeTa)Q^2}+\fRac{k_T^2+\kappa_t^2-2{\bf{\kappa_T}}.{\Bf{k_t}}+M_q^2}{\beta
Q^2}
 \label{eq:a}.\end{aligNed}$$ As In THe reFErEnce [@4kImber], thE ScALE $\mu$ WhIcH coNtROlS the $uNinTegrateD$ gluon and tHe $QcD$ Coupling consTAnT $\alpha_s$, is cHosEn as follows: $$\bEgiN{alignEd}
\Mu^2=k_T^2+\Kappa_T^2+m_q^2
 \Label{Eq:b}.\end{ALigned}$$ one should NoTE that tHE coeffiCIEnTs usEd for quark and non-pERtURbative gluon coNtribuTiONs DEPenD on The transveRsE momeNTum. As it HAs BEEN brIEfly explained Before, the maIN prEscripTiOn fOR $F_L$ conSists Of THreE terms; the fiRst tErm is the $k_T$ factoRIzation WHich expLains tHe cOntRibuTIoN oF thE $UpdF$ iNTo The $f_l$. ThIs term is DeRiVed wiTh thE USE Of puRe gLuon ContrIbution. HoweveR, it Only COunTs the Gluon ContRiButioNs comiNg froM tHe perturbative rEgioN, i.e. for $k_t>1$ $GEV$, aNd DoeS nOt havE AnythiNg tO do With the Non-pertURbaTiVE COnTributions. In the refErENCe [@Stasto4], it Has beeN ShOwN That a proPeR noN-perTURbatiVe teRM cAn be deriVed froM ThE $k_T$ factorIzAtion tErM, coMpaCting THe $k_t$ DependEnce and tHe intEGration with the USe of a variable-CHaNGE, i.E. Y, thaT caRries the $k_t$ dEpenDEnce. neveRThEleSS, therE is a cAlCUlABle quark contributioN iN the loNgituDinal structurE function oF THE proton, wHich COmES from the collinEar faCtorizatioN, I.e. the secOnd teRm of the eQuation (\[eq:2P\]).

fOr the chaRm qUarK, $m$ iS taKEN tO be $m_c=1.4 GeV$, and $u$, $d$ AND $s$ quArKs masseS arE neglecTed. we aLso Use ThE same apprOximatioN tO sAvE tHe cOmputATion time [@TkImbEr], The One we DId for tHe calCulaTiOn OF $F_2(x,q^2)$ [@mho1] i.e
_q}\over{\beta(1-\beta)Q^2 }}),$$ and  $$\b egi n{a li gned }
\f rac{1}{z}=1+\f r ac{\ kappa_{t}^{2}+m_q^2}{( 1-\be ta ) Q^2} + \f rac{k _t^2+\k a pp a _ t^2 -2 {\ bf{ \k a pp a_t}} .{\ bf{k_t} }+m_q^2}{\ bet aQ^2}
 \label { eq :a}.\end{a lig ned}$$ As in  th e refe re nce [@4ki mbe r], t he sca l e $\mu $ which c on t rols t h e $unin t e gr ated $ gluon and the $ Q CD $  coupling cons tant $ \a l ph a _ s$,  is  chosen as  f ollow s : $$\be g in { a l ign e d}
\mu^2=k_t^ 2+\kappa_t^ 2 +m_ q^2
 \ la bel { eq:b}. \end{ al i gne d}$$ One sh ould  note tha t thec oeffici e nts use d forqua rkandn on -p ert ur b ati v eglu o n c ontribut io ns  depe nd o n t h e tr ans vers e mom entum. As ithas  bee n  br iefly  expl aine dbefor e, the  main  p rescription for  $F_ L$ consis tsof  th re e ter m s; the  fi rst  term i s the $ k _t$  f a c t or ization which expl ai n s  t he contr ibutio n  o ft he $UPDF $int o th e $F_L$ . Th i sterm isderive d  w it h the u se  of pu re  gl uon  cont r ibut ion. H owever,it on l y counts the g l uon contribut i on s co m ingfro m the pertu rbat i ve r egio n ,i.e .  for$k_t> 1$ $G e V$, and does not ha ve  anyth ing t o do with the  non-pertu r b a tive con trib u ti o ns. In the ref erenc e [@stasto 4 ], it ha s bee n shownthat a pr o p er non-p ert urb ati vet e rm  can be deriv e d  fro mthe $k_ t$factori zat ion  te rm,  c ompacting  the $k_ t$  d ep en den ce an d  the int eg rat io n w ith t h e useof avari ab le - cha nge, i. e .y ,  tha tca rrie s t he  $k_t $ de p end ence. N everthele ss, ther eis  a calc ulable quarkco ntribution  i n t he lon g i tudinalstructure function of t h e proto n,which  com es from t hecollin ear factor izatio n, i. e.  th e secon d te rmof  the equat i o n ( \[eq: 2p \]). 

For t he charm quark, $m $  is  taken to be$m_ c=1. 4 Ge V$, an d  $u $, $d$ a nd $s$ quarks m asses arene g le cted. We a l sous e the s ame app roxim a tion to  save the  computat io n ti m e  [@ tkimber],the onewe did fo r  thec al culat ion  of $F _2 (x, Q^2)$  [@mho 1 ] i .e
_q}\over{\beta(1-\beta)Q^2}}),$$ and_$$\begin{aligned}
\frac{1}{z}=1+\frac{\kappa_{t}^{2}+m_q^2}{(1-\beta)Q^2}+\frac{k_t^2+\kappa_t^2-2{\bf{\kappa_t}}.{\bf{k_t}}+m_q^2}{\beta
Q^2}
 \label{eq:a}.\end{aligned}$$_As in the reference_[@4kimber], the_scale_$\mu$ which_controls_the $unintegrated$ gluon_and the $QCD$_coupling constant $\alpha_s$, is_chosen as follows:_$$\begin{aligned}
\mu^2=k_t^2+\kappa_t^2+m_q^2
_\label{eq:b}.\end{aligned}$$ One should note that the coefficients used for quark and non-perturbative gluon contributions_depend_on the_transverse_momentum._As it has been briefly_explained before, the main prescription_for $F_L$_consists of three terms; the first term is_the_$k_t$ factorization which_explains the contribution of the $UPDF$ into the $F_L$._This term is derived with the_use of pure_gluon_contribution._However, it only counts_the gluon contributions coming from the_perturbative region, i.e. for $k_t>1$ $GeV$,_and does not have anything to do_with the non-perturbative contributions. In the_reference [@stasto4], it has been_shown that_a proper non-perturbative term can_be derived from_the $k_t$_factorization term, compacting_the $k_t$ dependence and the integration_with the use_of a variable-change, i.e. y, that_carries_the $k_t$ dependence._Nevertheless,_there_is a_calculable quark contribution_in_the longitudinal_structure_function of the proton, which comes_from_the collinear factorization, i.e. the second term_of the equation (\[eq:2p\]).

For_the_charm quark, $m$ is_taken to be $m_c=1.4 GeV$,_and $u$, $d$ and $s$ quarks_masses are_neglected. We_also use the same approximation to save the computation time [@tkimber],_the one we did for the_calculation of $F_2(x,Q^2)$ [@mho1]_i.e