ToPrintXS.f

Go to the documentation of this file.

!!!!!!!
!     include this inside Gencol.f, if you want to check Sibyll
!     xs. 
!     The place is after the next comment line in Gencol.f

!     !!!!    you may include  ToPrintXS.f here to print out E vs XS for Sibyl
      
!     normal Gencol output will not be obtained. After printing the table
!     execution stops immediately.
!     In the  output table, you will find lines starting from "xsxs"
!                      
!             Ep           Log10(Ep)  pO xs    pAir xs 
!     xsxs     80.00         1.903   299.902   272.651
!     ... 
!     pO:   p may be pi or K depending on input collision spec.
!           O may be N depending on ..
!     pAir: p may be pi or K depending on input collision spec.
!     Air target case is normmaly ~ 0.78 N target + 0.22 O target
!      
      
         if(activemdl .eq. 'sibyll' ) then
            pj%fm%p(4)=80.
            if( pj%code == 6 ) then
               l = 1
            elseif(pj%code == 5 ) then
               l = 3
            elseif(pj%code == 4 ) then
               l = 2
            endif
               
           do while(   pj%fm%p(4)< 1e11)
              roots = sqrt(2*masp*pj%fm%p(4) + masp**2 + pj%mass**2)
              call sib_sigma_hair (l, roots, xsair, sigbdif)


              pj%fm%p(3)=sqrt(pj%fm%p(4)**2-pj%mass**2)
              call csibyllxs(pj, tg, xs)
              if( xs /= xs ) then
                 write(*,*) ' nan by if '
              endif
              if( isnan(xs) ) then
                 write(*,*) ' nan by isnan '
              endif
              write(*,'(a, 1p, g13.4, 0p, 3f10.3)')
     *              "xsxs ", pj%fm%p(4),log10(pj%fm%p(4)), xs, xsair
              pj%fm%p(4) =  pj%fm%p(4)*10.0d0**0.2
           enddo
           stop
         endif
!!!!!!!!!!!!