Changeset 184 in JEM-EUSO for esaf_lal/branches/camille/packages/simulation/detector/G4Detector/G4fresnellens/src/FresnelSurface.cc

Timestamp:

Jul 15, 2013, 2:08:33 PM (11 years ago)

Author:

moretto

Message:

changes from biktem r:3032

File:

• esaf_lal/branches/camille/packages/simulation/detector/G4Detector/G4fresnellens/src/FresnelSurface.cc (modified) (3 diffs)

esaf_lal/branches/camille/packages/simulation/detector/G4Detector/G4fresnellens/src/FresnelSurface.cc

@@ -46 +46 @@
     }
     fSurfacesIns.clear();
+}
+
+//__________________________________________________________________________________
+int FresnelSurface::FindTooth( double rho ) { 
+    return BinarySearch(theToothEdges, rho, 0, theToothN); 
+} 
+
+//__________________________________________________________________________________
+double FresnelSurface::PointOnSurf(double rho) { 
+    int n_tooth = FindTooth( rho );
+    if( n_tooth<0 ) return kInfinity;
+    return fSurfacesIns[ n_tooth*2 ]->PointOnSurf( rho );   
 }
 
@@ -180 +192 @@
       }
     #endif
-    int n_tooth = BinarySearch(theToothEdges,rho,0,theToothN);
+    int n_tooth = FindTooth( rho );
     if( n_tooth<0 ) {
         printf("Warning. Tooth not found. Shouldn't happen (rho=%g of [%g, %g]).",
@@ -255 +267 @@
     }// approach to the limiting surface
 
-    int theCurrentToothIndex=BinarySearch(theToothEdges, perp, 0, theToothN);
+    int theCurrentToothIndex=FindTooth( perp );
     if ( theCurrentToothIndex<0 ) return kInfinity;
     vector<SSurface*>::const_iterator it=fSurfaces.begin()+theCurrentToothIndex*2;