Hi, it seems to me that the PseudoRapidity in 2.23/04 has a bug. Is there a minus sign missing? See below the line I changed to get it working. Bye, Mike *** PHYSICS_Vector3.cxx.orig Wed Oct 20 14:32:22 1999 --- PHYSICS_Vector3.cxx Wed Oct 20 14:33:08 1999 *************** *** 274,280 **** //return 0.5*log( (m+fZ)/(m-fZ) ); // guard against Pt=0 double cosTheta = CosTheta(); ! if (cosTheta*cosTheta < 1) return 0.5* TMath::Log( (1.0-cosTheta)/(1.0+cosTh eta) ); Warning("PseudoRapidity","transvers momentum = 0! return +/- 10e10"); if (fZ > 0) return 10e10; else return -10e10; --- 274,280 ---- //return 0.5*log( (m+fZ)/(m-fZ) ); // guard against Pt=0 double cosTheta = CosTheta(); ! if (cosTheta*cosTheta < 1) return -0.5* TMath::Log( (1.0-cosTheta)/(1.0+cosT heta) ); Warning("PseudoRapidity","transvers momentum = 0! return +/- 10e10"); if (fZ > 0) return 10e10; else return -10e10; -------------------------------------------------------------------------- Michael Sievers Michael.Sievers@desy.de -------------------------------------------------------------------------- esa$ gcc -Wall -o ariane5 ariane4.c ariane4.c: 666: warning: long float implicitly truncated to unsigned type esa$ ariane5
This archive was generated by hypermail 2b29 : Tue Jan 04 2000 - 00:43:41 MET