Problem in calculation of y corrected.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 19 Jul 2002 08:11:54 +0000 (08:11 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 19 Jul 2002 08:11:54 +0000 (08:11 +0000)
EVGEN/TestPrimaries.C

index a369b695f4d42589fcd793280f93246f7142602a..f0f02f1102a20d22d84d65544876abc1fce5809f 100644 (file)
@@ -71,7 +71,7 @@ void TestPrimaries(Int_t evNumber1=0, Int_t evNumber2=0)
            Float_t theta = MPart->Theta();
            Float_t phi   = MPart->Phi()-TMath::Pi();
            Float_t eta   = -TMath::Log(TMath::Tan(theta/2.));
            Float_t theta = MPart->Theta();
            Float_t phi   = MPart->Phi()-TMath::Pi();
            Float_t eta   = -TMath::Log(TMath::Tan(theta/2.));
-           Float_t y     = 0.5*TMath::Log((E+Pz)/(E-Pz+1.e-13));
+           Float_t y     = 0.5*TMath::Log((E+Pz+1.e-13)/(E-Pz+1.e-13));
            
            if (child1 >= 0) continue;
            if (mpart == kPi0 || mpart == kGamma) continue;
            
            if (child1 >= 0) continue;
            if (mpart == kPi0 || mpart == kGamma) continue;