Changeset 18275

Timestamp:

07/26/15 22:44:30 (9 years ago)

Author:

Daniela Dorner

Message:

fixed bug in calculation of fConcCore

File:

• trunk/Mars/mimage/MNewImagePar.cc (modified) (1 diff)

trunk/Mars/mimage/MNewImagePar.cc

@@ -220 +220 @@
         const Double_t dzx   =  hillas.GetCosDelta()*dx + hillas.GetSinDelta()*dy; // [mm]
         const Double_t dzy   = -hillas.GetSinDelta()*dx + hillas.GetCosDelta()*dy; // [mm]
-        const Double_t dz    =  gpix.GetT()*gpix.GetT()/4;
-        const Double_t tana  =  dzy*dzy/(dzx*dzx);
-        const Double_t distr =  (1+tana)/(rl + tana*rw);
-        if (distr>dist0-dz || dzx==0)
+        const Double_t dz    =  gpix.GetT()/2;
+        const Double_t distr =  (dzy*dzy+dzx*dzx)/(dzx*dzx*rl + dzy*dzy*rw);
+        if (sqrt(distr)>sqrt(dist0)-dz || (dzx==0 && dzy==0))
             fConcCore += nphot;