SUBROUTINE delphi(dphi) $$INCLUDE 'C$INC:BCS.INC' $$INCLUDE 'YBOS$LIBRARY:ERRCOD.INC' REAL dphi ! difference in phi's REAL phi1 ! phi of muon REAL phi2 ! phi of electron C................................. dphi = -9999. status = blocat(iw,'MUON',1,ind,indmuo) IF(status .NE. yesucc) RETURN status = blocat(iw,'ELEC',1,ind,indele) IF(status .NE. yesucc) RETURN phi1 = rw(indmuo+15) phi2 = rw(indele+5) dphi = ABS(phi1-phi2) RETURN END
This routine will probably crash.