Functions/Subroutines
program	democache

subroutine	getinvariants (MomInv6, ievent)

double complex function, dimension(0:n-2)	submasses (N, k, masses)

double complex function, dimension((n-1) *(n-2)/2)	submominv (N, k, MomInv)

Function/Subroutine Documentation

◆ democache()

program democache

Definition at line 10 of file democache.f90.

◆ getinvariants()

subroutine democache::getinvariants	(	double complex, dimension(1:15)	MomInv6,
		integer	ievent
	)

Definition at line 198 of file democache.f90.

  
     implicit none
     
     double precision p(1:6,0:3) 
     double complex MomInv6(1:15) 
     double precision :: s(0:5,0:5),fact=47d0/49d0
     integer :: N=6,k,i,newmom,ievent
     save s,fact,n
  
     newmom=0
  
     if (ievent.eq.1) then
       newmom=1
 ! use momentum set from Racoon
       p(1,0) =   -0.9500000000000000d+02
       p(1,1) =    0.0000000000000000d+00
       p(1,2) =    0.0000000000000000d+00
       p(1,3) =    0.9500000000000000d+02
       
       p(2,0) =   -0.9500000000000000d+02
       p(2,1) =    0.0000000000000000d+00
       p(2,2) =    0.0000000000000000d+00
       p(2,3) =   -0.9500000000000000d+02
       
       p(3,0) =    0.2046111001757171d+02
       p(3,1) =    0.1057734233089455d+02
       p(3,2) =   -0.2324961261504543d+01
       p(3,3) =    0.1736005205921753d+02
       
       p(4,0) =    0.3558305163378869d+01
       p(4,1) =    0.1436222934374051d+01
       p(4,2) =   -0.2174258125294355d+01
       p(4,3) =   -0.2423097382091398d+01
       
       p(5,0) =    0.8154540918019539d+02
       p(5,1) =   -0.5230395944682889d+02
       p(5,2) =    0.3083642435466509d+02
       p(5,3) =    0.5443403822581044d+02
       
       p(6,0) =    0.8443517563885433d+02
       p(6,1) =    0.4029039418156027d+02
       p(6,2) =   -0.2633720496786619d+02
       p(6,3) =   -0.6937099290293661d+02
       
  
     else if (ievent.eq.50) then
       newmom=1
 ! set with small momentum   1.9 versus 3.0 seconds
       p(1,0) =   -250.0000000000000d0       
       p(1,1) =    0.0000000000000000d+00
       p(1,2) =    0.0000000000000000d+00
       p(1,3) =   -250.0000000000000d0
       
       p(2,0) =   -250.0000000000000d0
       p(2,1) =    0.0000000000000000d+00
       p(2,2) =    0.0000000000000000d+00
       p(2,3) =    250.0000000000000d0
       
       p(3,0) =    239.91697591218110d0
       p(3,1) =   31.640107125154238d0
       p(3,2) =   -116.12928664566658d0
       p(3,3) =   -207.54047252312847d0
  
       p(4,0) =    10.276659578244022d0
       p(4,1) =    1.2584629023924521d0
       p(4,2) =    10.131487190720421d0
       p(4,3) =   -1.1742957526504370d0
  
       p(5,0) =    181.18800109198216d0
       p(5,1) =    -35.098898474333737d0
       p(5,2) =    109.26326872331617d0
       p(5,3) =    140.20947604742324d0
  
       p(6,0) =      68.618363417592718d0
       p(6,1) =    2.2003284467870472d0
       p(6,2) =    -3.2654692683700119d0
       p(6,3) =    68.505292228355671d0
  
     else if (ievent.eq.100) then
       newmom=1
         p(1,0) = -100.0000000000000d0
         p(1,1) = 0.000000000000000d0
         p(1,2) = 0.000000000000000d0
         p(1,3) = -100.0000000000000d0
  
         p(2,0) = -100.0000000000000d0
         p(2,1) = 0.000000000000000d0
         p(2,2) = 0.000000000000000d0
         p(2,3) = 100.0000000000000d0
  
         p(5,0) = 42.13045937773281d0
         p(5,1) = 38.71045260055064d0
         p(5,2) = 16.54764683624944d0
         p(5,3) = 1.628450497276325d0
  
         p(6,0) = 73.31699208541762d0
         p(6,1) = -66.17895647220249d0
         p(6,2) = 31.53024134467396d0
         p(6,3) = -1.253367244715983d0
  
         p(4,0) = 69.60509150033207d0
         p(4,1) = 34.42307953372518d0
         p(4,2) = -60.31848303205667d0
         p(4,3) = 4.647683605402023d0
  
         p(3,0) = 14.94745703651754d0
         p(3,1) = -6.954575662073326d0
         p(3,2) = 12.24059485113326d0
         p(3,3) = -5.022766857962365d0
     else if (ievent.eq.150) then
       newmom=1
       p(1,0) =   -100.0000000000000d0       
       p(1,1) =    0.0000000000000000d+00
       p(1,2) =    0.0000000000000000d+00
       p(1,3) =   -100.0000000000000d0
       
       p(2,0) =   -100.0000000000000d0
       p(2,1) =    0.0000000000000000d+00
       p(2,2) =    0.0000000000000000d+00
       p(2,3) =    100.0000000000000d0
       
       p(3,0) =    72.02905523601244d0
       p(3,1) =   32.33193950844802d0
       p(3,2) =   64.10893732125734d0
       p(3,3) =   5.733641195059672d0
  
       p(4,0) =    25.07519219892066d0
       p(4,1) =    2.906699172283076d0
       p(4,2) =    -24.65281899347465d0
       p(4,3) =    3.543286526606821d0
  
       p(5,0) =    60.19609558319355d0
       p(5,1) =    2.282219923475870d0
       p(5,2) =    -59.72688022443283d0
       p(5,3) =    -7.145710205299344d0
  
       p(6,0) =    42.69965698187335d0
       p(6,1) =    -37.52085860420697d0
       p(6,2) =    20.27076189665014d0
       p(6,3) =    -2.131217516367149d0
     else if (ievent.eq.200) then
       newmom=1
       p(1,0) =   -100.0000000000000d0       
       p(1,1) =    0.0000000000000000d+00
       p(1,2) =    0.0000000000000000d+00
       p(1,3) =   -100.0000000000000d0
       
       p(2,0) =   -100.0000000000000d0
       p(2,1) =    0.0000000000000000d+00
       p(2,2) =    0.0000000000000000d+00
       p(2,3) =    100.0000000000000d0
       
       p(3,0) =    69.63925440091566d0
       p(3,1) =   -68.33683461638704d0
       p(3,2) =   -3.803463998624776d0
       p(3,3) =   -12.85443307731553d0
  
       p(4,0) =    34.49634785978589d0
       p(4,1) =    21.54290146454439d0
       p(4,2) =    -15.60187514686951d0
       p(4,3) =    21.96549348532329d0
  
       p(5,0) =    66.20628864972379d0
       p(5,1) =    65.60907499667245d0
       p(5,2) =    -3.500460884404577d0
       p(5,3) =    -8.152834381263276d0
  
       p(6,0) =    29.65810908957467d0
       p(6,1) =    -18.81514184482980d0
       p(6,2) =    22.90580002989887d0
       p(6,3) =    -0.9582260267444891d0
  
     else if (ievent.eq.250) then
       newmom=1
         p( 1, 0) =      -159.2798637127597488d0
         p( 1, 1) =         0.0000000000000000d0
         p( 1, 2) =         0.0000000000000000d0
         p( 1, 3) =      -159.2798637127597488d0
         p( 2, 0) =      -237.6261878491357891d0
         p( 2, 1) =         0.0000000000000000d0
         p( 2, 2) =         0.0000000000000000d0
         p( 2, 3) =       237.6261878491357891d0
         p( 3, 0) =       133.6018149124976446d0
         p( 3, 1) =        60.5134748218625731d0
         p( 3, 2) =       112.6312302087531094d0
         p( 3, 3) =       -38.7526811273741174d0
         p( 4, 0) =        69.9256648527542950d0
         p( 4, 1) =        66.6463990867762135d0
         p( 4, 2) =         7.6538437283820429d0
         p( 4, 3) =       -19.7300473909055754d0
         p( 5, 0) =        85.7314709474216841d0
         p( 5, 1) =       -38.6625297933730749d0
         p( 5, 2) =       -63.4432116127613241d0
         p( 5, 3) =       -42.7791164126088006d0
         p( 6, 0) =       107.6471008492219283d0
         p( 6, 1) =       -88.4973441152657188d0
         p( 6, 2) =       -56.8418623243738281d0
         p( 6, 3) =        22.9155207945124531d0
  
  
     end if
  
     if (newmom.eq.1) then
       ! define invariants
       n=6 
       do k=0,n-1
         s(modulo(k+1,n),k) = p(k+1,0)**2
         do i=1,3
           s(modulo(k+1,n),k) = s(modulo(k+1,n),k) - p(k+1,i)**2
         end do
         if(abs(s(modulo(k+1,n),k)).lt.1d-14* abs(p(k+1,0))**2) s(modulo(k+1,n),k) = 0d0
       
         s(k,modulo(k+1,n)) = s(modulo(k+1,n),k)
       end do
  
       do k=0,n-1
         s(modulo(k+2,n),k) = (p(modulo(k+1,n)+1,0) + p(k+1,0))**2
         do i=1,3
           s(modulo(k+2,n),k) = s(modulo(k+2,n),k) &
               - (p(modulo(k+1,n)+1,i) + p(k+1,i))**2
         end do
         s(k,modulo(k+2,n)) = s(modulo(k+2,n),k)
       end do
       do k=0,(n-1)/2
         s(modulo(k+3,n),k) = (p(modulo(k+2,n)+1,0) + p(modulo(k+1,n)+1,0) + p(k+1,0))**2
         do i=1,3
           s(modulo(k+3,n),k) = s(modulo(k+3,n),k) &
               - (p(modulo(k+2,n)+1,i) + p(modulo(k+1,n)+1,i) + p(k+1,i))**2
         end do
         s(k,modulo(k+3,n)) = s(modulo(k+3,n),k)
       end do
       
       do k=1,n
         mominv6(k) = s(modulo(k,n),k-1)
       end do
       do k=1,n
         mominv6(k+n) = s(modulo(k+1,n),k-1)
       end do
       do k=1,n/2
         mominv6(k+2*n) = s(modulo(k+2,n),k-1)
       end do
       fact = 47d0/97d0  
  
     else
 !     rescale invariants
       fact = 4*fact*(1-fact)
  
       do k=1,n
         mominv6(k) = s(modulo(k,n),k-1)*fact
       end do
       do k=1,n
         mominv6(k+n) = s(modulo(k+1,n),k-1)*fact
       end do
       do k=1,n/2
         mominv6(k+2*n) = s(modulo(k+2,n),k-1)*fact
       end do
     end if
       

◆ submasses()

double complex function, dimension(0:n-2) democache::submasses	(	integer	N,
		integer	k,
		double complex, dimension(0:n-1)	masses
	)

Definition at line 468 of file democache.f90.

  
     integer :: N,k,i
     double complex :: masses(0:N-1), SubMasses(0:N-2)
  
     do i=0,n-2
       if (i.lt.k) then
         submasses(i) = masses(i)
       else
         submasses(i) = masses(i+1)
       end if
     end do    
  

◆ submominv()

double complex function, dimension((n-1)*(n-2)/2) democache::submominv	(	integer, intent(in)	N,
		integer, intent(in)	k,
		double complex, dimension(n*(n-1)/2)	MomInv
	)

Definition at line 495 of file democache.f90.

  
     integer, intent(in) :: N,k
     integer :: i,cnt,moms,limit
     double complex :: MomInv(N*(N-1)/2), SubMomInv((N-1)*(N-2)/2)
  
     if (n.eq.2) return
  
     cnt = 1
     do moms=1,(n-3)/2
       do i=1,n
         if (i.ne.k+1) then
           if ((k.ge.i).and.(i.ge.k-moms+1).or.(n+k.ge.i).and.(i.ge.n+k-moms+1)) then
             submominv(cnt) = mominv(n*moms+i)
           else
             submominv(cnt) = mominv(n*(moms-1)+i)
           end if
  
           cnt = cnt+1
         end if
       end do
     end do
  
     if (mod(n,2).eq.1) then
       moms = (n-1)/2
       if(k.le.(n-1)/2) then 
         limit = (n+1)/2
       else
         limit = (n-1)/2
       end if
       do i=1,limit
         if (i.ne.k+1) then
           if ((k.ge.i).and.(i.ge.k-moms+1).or.(n+k.ge.i).and.(i.ge.n+k-moms+1)) then
             submominv(cnt) = mominv(n*moms+i-(n-1)/2)
           else
             submominv(cnt) = mominv(n*(moms-1)+i)
           end if
  
           cnt = cnt+1
         end if
       end do
  
     else
       do i=1,n
         if (i.ne.k+1) then
           moms = (n/2-1)
           if ((k.ge.i).and.(i.ge.k-moms+1).or.(n+k.ge.i).and.(i.ge.n+k-moms+1)) then
             if (i.gt.n/2) then
               submominv(cnt) = mominv(n*moms-n/2+i)
             else
               submominv(cnt) = mominv(n*moms+i)
             endif
           else
             submominv(cnt) = mominv(n*(moms-1)+i)
           end if
  
           cnt = cnt+1
         end if
       end do
  
     end if  
  
  

Functions/Subroutines

Function/Subroutine Documentation

◆ democache()

◆ getinvariants()

◆ submasses()

◆ submominv()