c Subroutine Convpagp(nrclus,Ierr) c c Author G.Cases 14/12/93 c Implicit None C #include "partap.inc" #include "clucom.inc" #include "genjet.inc" #include "morsub.inc" #include "conecut.inc" C integer i,j,k,nrclus,ierr,nfinal real theta,phi,E,ET,eta,pt real deg,atg2 real x,y,z,px,py,pz,rcl deg=180./3.14159 nfinal=0 njet1=0 Do i=Npart-NrClus+1,Npart theta=2*atan(exp(-PC(i,1))) phi=PC(i,2) E=PC(i,4) ET=PC(i,3) c---- p(1,903)= ET*cos(phi) p(2,903)= ET*sin(phi) p(3,903) = E*cos(theta) p(4,903) = E p(5,903) = 0. c PT = SQRT(P(1,903)**2 + P(2,903)**2) theta = acos(P(3,903) / sqrt(PT**2+P(3,903)**2)) ETA = -alog(tan(theta / 2.)) c if(ET.gt.Conejet_ptmin) then if (eta.gt.Conejet_etamin.and.eta.lt.Conejet_etamax) then nfinal=nfinal+1 do j=1,5 PGJET1(nfinal,j)=p(j,903) enddo PGJET1(nfinal,6)=PC(i,1) PGJET1(nfinal,7)=PT PGJET1(nfinal,8)= + atan2(PGJET1(nfinal,2),PGJET1(nfinal,1)) PGJET1(nfinal,9)=PC(i,1) PGJET1(nfinal,10)=ET PGJET1(nfinal,11)=phi endif endif c enddo njet1=nfinal ierr=0 return end