Abstract:
The invention concerns a novel method adapted to experimental conditions, for determining the function defining the actual detecting sensitivity of a solid-state detector of nuclear traces and which represents the detecting performances for each type of nucleus emitting energy alpha particles Ei, located in each position inside a sample to be analyzed, and hence precision dosimetry of the radioactive nuclei in a sample using nuclear trace solid-state detectors. The method provides an improvement of the quality of sample analyses, a significant reduction of analysis costs and an important widening of the scope of application of this type of detectors. Thus, the method will bring about the development of a new generation of detectors having a better price-quality ratio which can be found in the field of radioactive nuclei analyses in different types of samples.

Description:
RELATED U.S. APPLICATIONS  
       [0001]     Not applicable.  
       STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT  
       [0002]     Not applicable.  
       REFERENCE TO MICROFICHE APPENDIX  
       [0003]     Not applicable.  
       FIELD OF THE INVENTION  
       [0004]     The present invention relates to a method for direct dosimetry of radioactive nuclei in different types of samples using solid-state detectors of nuclear traces.  
       BACKGROUND OF THE INVENTION  
       [0005]     The history of solid state nuclear track detectors began in 1958 after initial observations by D. A. Young of fission fragment tracks of lithium fluoride (LiF) irradiated by thermal neutrons. In 1959, Silk and Barnes made these same observations on mica using an electron microscope. Numerous observations have been made since then on a large number of materials. The detectors have the ability to record the passage of particles such as alpha particles, protons, and fission fragments.  
         [0006]     In the years that have passed, work by the international scientific community on solid state nuclear track detectors made it possible to develop several dosimetry methods for radioactive nuclei with the use of these detectors. Though these methods use basic knowledge of calculation of the probability that an alpha particle emitted from a radioactive nucleus of a sample can reach and be recorded on a solid state nuclear track detector, they are based on estimation of the global critical angle depending on the sample, the characteristics of the particles to be detected, and the detector. The estimation thus induces an error in the calculation and makes the solid state nuclear track detectors (SSNTDs) imprecise.  
         [0007]     Other methods of detection that do not use solid state nuclear track detectors already exist. They involve, for example, geochemical analyses and instrumental techniques such as analysis by neutron activation, analysis by mass spectometry or analysis by gamma spectrometry and liquid scintillation analysis. All of these techniques, often destructive to the samples analyzed, use equipment that is heavy, sophisticated, and especially costly.  
       BRIEF SUMMARY OF THE INVENTION  
       [0008]     In order to correct these disadvantages, a new method has been developed. It is based on the implementation of a new technique specified for calculating the sensitivity of a solid state nuclear track detector for each type of nucleus emitting alpha particles having energy E 1  located in each position in the sample to be analyzed. This method does not allow the estimation but instead the calculation without errors of the critical angle of incidence depending on the sample, the characteristics of the alpha particles to be detected, and the detector. Thus, this method is based on the actual function of the critical angle of incidence which makes it possible to solve the problem of the sample-detector geometry, and as a result, determination of the function that defines the real sensitivity of detection of a solid state nuclear track detector (SSTND) and which represents the efficiency of detection for each type of nucleus emitting alpha particles having energy E 1  located in each postion in the middle of the sample to be analyzed.  
         [0009]     Applications of this method allow an improvement in the quality of sample analyses, a considerable reduction of costs for analysis and a great enlargement in the field of application for this type of detector. Thus, a new generation of detectors will be born having a better price-to-quality ratio than has existed until now in the domain of analyses of radioactive nuclei in different types of samples. 
     
    
     BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS  
       [0010]    
       FIG. 1 
     
         [0011]    
       FIG. 2 
     
         [0012]    
       FIG. 3 
     
         [0013]    
       FIG. 4 
     
         [0014]    
       FIG. 5 
     
     
    
     DETAILED DESCRIPTION OF THE INVENTION  
       [0015]     II. Methodology and Procedures for Evaluation of Different Parameters  
         [0016]     Let (dv) be the volume element of a cylinder having a radius between r and r+dr and a height between h and h+dh inside the sample analyzed, see FIG. II- 1  in the drawings section where: (1: SSNTD, 2: plastic capsule, 3: sample analyzed). 
 
 dv= 2π rdrdh  
 
         [0017]     Let dN i  be the number of the alpha particle having energy E i  emitted by the radioactive element having index i contained in the volume dv and for the exposure time te: 
 
 dN   i   =A   ci   te   ·dv  
 
         [0018]     Where A ci   te  is the alpha activity of the radioactive nuclei of the emitter element i of the alpha particle having energy E i  per unit volume and for the exposure time te.  
         [0019]     It is thus necessary to multiply the number dN i  of alpha particles emitted in the volume dv by a probability of detection P i (r,h), [P i (r,h) is the probability that an alpha particle of energy E i  emitted by the nucleus of index i contained in the volume dv can be detected by the solid state nuclear track detector] and which also represents the efficiency of detection of alpha particles of energy E i  emitted by the nucleus of index i contained in the volume dv.  
         [0020]     As a result, for radioactive nuclei of element i contained in the volume dv, the number of tracks developed will be: 
 
 dn   i   =A   ci   te   ·P   i ( r,h )· dv  
 
         [0021]     Where: dn i  is the number of alpha particles having energy E i  emitted by the nuclei of index i contained in the volume dv, for the time of exposure (te) and which are detected in our solid state nuclear track detector.  
         [0022]     Actually, the total number of tracks observed due to the radioactive nuclei of index i contained in the total volume is: 
 
 n   it =2π· A   ci   ie   ·∫∫P   i ( r,h )· rdrdh  
 
 II-1 Method for Calculating the Probability P i (r,h). 
 
         [0023]     Let a nucleus of index i denoted K i , an emitter of alpha particles of energy E i , be located in a position of coordinates (r,h) (see FIG. II- 5  in the drawings section) such that these emitted particles give rise to observable tracks on the solid state nuclear track detector. The particles must satisfy two conditions: 
        1—they are emitted at a well-defined angle of incidence by means of the function of the critical angle of incidence (θ ic ).     2—they reach the detector. 
 
 II-1-1 Identification and Evaluation of the Function (θ ic ). 
       
 
         [0026]     A) Identification of the Function (θ ic ).  
         [0027]     FIG. II- 2  in the drawings section shows the experimental and intrinsic function of a SSNTD describing the variation of the critical angle (θ c ) resolving the tracks as a function of the residual energy for conditions enabling chemical resolution and reading of well-defined tracks. This function is from a solid state nuclear track detector that is currently the most used (CR-39; LR-115). Three possible cases can be distinguished: 
 
*E ind &lt;E min , *E ind ≦E min ≦E min , *E ind &gt;E max  
 
 where: 
        E ind  is the residual energy at normal incidence on the detector of the alpha particle having energy E i  emitted from a nucleus located at a height h from the plane of the detector.     E min  is the minimum energy detectable at normal incidence.     E max  is the maximum energy detectable at normal incidence.         Comment: The identification of (θ ic )=f(h), (θ ic1 )=f a (h), and (θ ic2 )=f b (h) can also be made by means of the inverse function of the function (θ c )=f(E res ):E res =f(θ c ) which shows the minimum and maximum energies detectable for each incidence (θ i ). 
 
 A-1—Case of Nuclei Emitters of Alpha Particles Having Energy E i  and with E ind &lt;E min , 
   
 
         [0032]     These are nuclei located in the volume for contribution to the tracks defined by h∈[R i −R min ,R i ] for the case without an absorbent medium, given R i  and the total trajectory in the sample of the alpha particle having energy E i  in the sample, and given R min  and the total trajectory in the sample of an alpha particle having energy E=E min . All alpha particles having energy E i  emitted from these nuclei will reach the detector with an energy lower than E min  and as a result will not be detected.  
         [0000]     A-2—Case of Nuclei Emitters of Alpha Particles Having Energy E i  and E min ≦E ind ≦E max .  
         [0033]     These are nuclei located in the volume for contribution to the tracks defined by h∈[R i −R max , R i −R min ] for the case without an absorbent medium, given R max  and the total trajectory in the sample of an alpha particle having energy E=E max .  
         [0034]     For a better explanation, according to FIG. II. 2 , there are two cases: E ind  less than the intermediate energy denoted (E m ) and E ind  greater than the intermediate energy E m .  
         [0035]     1—case where E ind  is less than E m :  
         [0036]     FIG. II- 3  in the drawings section shows the variation in the sample space, from (θ i ) and (θ c ) as a function of the residual energy and trajectory (X), of an alpha particle of energy E i  and E ind ∈[E min , E m ] emitted of a radioactive core located at the height h from the surface of the detector. It is observed that: 
        when the residual energy decreases, (θ c ) decreases and (θ i ) increases.     for (θ i )=0 (normal incidence for each position h), the corresponding (θ c ) has the larger value for each position h, and when (θ i ) increases, (θ c ) decreases. Thus, (θ ic ) for each position h is none other than the intersection between the angle of incidence of the particle to be detected (θ i ) and the resolution angle of the detector (θ c ), because from this value, (θ i ) is greater than (θ c ).        
 
         [0039]     2—case where E ind  is greater than E m ;  
         [0040]     At first, in the interval of residual energy [E m , E ind ], it is observed that when the residual energy decreases, (θ i ) increases and (θ c ) increases.  
         [0041]     Secondly, in the interval of residual energy [θ, E m ], it is observed that when the residual energy decreases, (θ i ) increases and (θ c ) decreases.  
         [0042]     For (θ i )=0 (normal incidence for each position h), the corresponding angle (θc) is also greater than zero for each position h and when (θ i ) increases, (θ c ) increases and then decreases, thus when the function appears as (θ c )=F(E res ) {strong phase of decrease of (θ c ) for E res ∈[E m , E i ]}, (θ i ) will have one and only one intersection with (θ c ) for each position h, i.e. in the first interval or in the second interval of energy, from this value of intersection, (θ i ) is greater than (θ c ).  
         [0043]     Thus, in a general manner, when the function appears as (θ c )=F(E res ), for each position h where E min ≦E ind ≦E max , (θ ic ) is none other than the intersection between the angle of incidence of the particle to be detected (θ i ) and the resolution angle of the detector (θ c ).  
         [0000]     A-3—Case of Nuclei Emitters of Alpha Particles Having Energy E i  and for which E ind &gt;E max ,  
         [0044]     These are nuclei located in the volume for contribution to the tracks defined by h∈[0, R i −R max ], for the case without an absorbent medium, where R max  and the total trajectory in the sample of an alpha particle having energy E=E max .  
         [0045]     In the same manner as above, when the function has the appearance (θ c )=F(E res ), there will be either zero intersections, or 2 intersections between (θ i ) and (θ c ) for each (h) which will be denoted (θ ic1 ) and (θ ic2 ) which correspond respectively to the maximum energy and to the one that is the minimum detectable at oblique incidence for each (h).  
         [0046]     Thus, for E ind &gt;E max , only the particles emitted at angles between (θ ic1 ) and (θ ic2 ) will be detected if they reach the detector.  
         [0047]     B) Evaluation of the Function (θ ic ).  
         [0048]     In order to determine (θ ic ), (θ ic1 ), and (θ ic2 ) as a function of h, for the case without an absorbent medium, proceed as follows: 
 
(θ ic )= f   1 ( E   res )= f   2 ( x ),  E   res   ∈E[ 0 , E   i   ], x∈E[ 0 ,R   i ]
        R i  is the total trajectory of the particle of energy E i  in the medium analyzed.     x: is the trajectory of the particle of energy E i  in the medium analyzed, this particle reaches the detector with an energy E res , x will be determined from the relation between energy and particle in the medium analyzed [2]. 
 
(θ i )= arc cos( h/x )= g   2 ( h,x )= g   1 ( h,E   res ),  x∈[h,R   i ]
    h is the position of the radioactive nucleus emitter of the particle of energy E i  relative to the surface of the detector.     x: is the trajectory of the particle of energy E i  in the medium analyzed and varies from h to R i  for each position h.        
 
         [0053]     In fact, (θ ic ) will be determined by solving the equation: g 2 (h,x)=f 2 (x) for the heights h i  where it is confirmed that E min ≦E ind ≦E max , (θ ic1 ) and (θ ic2 ) will be determined by solving the equation: g 2 (h,x)=f 2 (x) for the heights h i  where it is confirmed that E ind &gt;E max . Thus, a determination is done of (θ ic )=f(h) for the values of h for which E min ≦E ind ≦E max , (θ ic1 )=f a (h) and (θ ic2 )=f n (h) for the values of h for which E ind &gt;E max .  
         [0054]     It is possible to convert from a function f 1 (E res ) to its equivalent f 2 (x) for each type of alpha particle having energy (E i ) and total trajectory in the sample analyzed (R i ), emitted from a nucleus located in a position (h) within the sample to be analyzed, either by the slope of the function of Linear Energy Transfer (L.E.T.=dE/dx) of the alpha particles in a medium, or by the Energy-Trajectory relation [E=f(R)] of the alpha particles in a medium: 
        For the case without an absorbing medium:     Using the function of Linear Energy Transfer (dE/dx) in the medium to be analyzed:  
         E   res     =         ∫   0       R   i     -   x       ⁢           ⁢         ⅆ   E       ⅆ   x       ·     ⅆ   x         =       E   i     -       ∫       R   i     -   x       R   i       ⁢         ⅆ   E       ⅆ   x       ·           ⁢     ⅆ   x                 
 
 thus an E res  corresponds for each x and, as a result, the conversion of a function f 1 (E res ) to its equivalent f 2 (x) and vice-versa. 
    Using the Energy-Trajectory relation, [E=f(R)] in the medium to be analyzed:        
 
         [0058]     The function E=f(R) is used and the change in the variable (R) is made by (R i −x), this function becomes: E res = 31   i −x), thus for each x, the corresponding E res  will be determined and as a result, the conversion from one function f 1 (E res ) to its equivalent f 2 (x) and vice-versa. 
        For the case with an absorbing medium:     Using the Energy-Trajectory relation [E=f(R)]:        
 
         [0061]     According to FIG. II- 4  shown in the drawings section, where [(M1): medium to be analyzed. (M2): absorbent medium. (h a ): thickness of the absorbent medium. x 1  and x 2  are, respectively, the trajectory of the alpha particle in the analyzed medium and the one in the absorbent medium.]:  
         [0062]     Given:  
           cos   ⁡     (     θ   i     )       =         h   -     h   a         x   1       =         h   a       x   2       =     h   x           ,     x   =       x   1     +     x   2             
 
 i.e.: 
        E=f(R), the Energy-Trajectory function of alpha particles in the medium to be analyzed.     R=g(E), the Energy-Trajectory function of alpha particles in the absorbent medium.     E res1  and E res  are, respectively, the residual energy of the incident alpha particle at the input of the absorbant medium and that at the input of the SSNTD.        
 
         [0066]     Given E res1 =f(R i −x i ), thus R lab =g(E res1 )=g∘f(R i −x 1 ). Where R lab  is the total trajectory of the alpha particle having energy E res1  in the absorbent medium. 
 
 As a result, [R lab −x 2 ]=[g∘f(R i −x 1 )−x 2 ]=g(E res ), or  
           x   2     =         x   1     ·       h   a       h   -     h   a           =     x   ·       h   a     h           ,       
 
 In fact,  
           g   ·     f   ⁡     (       R   i     -     x   ·       h   a       h   -     h   a             )         -     x   ·       h   a     h         =     g   ⁡     (     E   res     )           
 
         [0067]     As a result, for each (h), each x corresponds to a residual energy E res  at the input of the SSNTD. Thus, for each position (h i ), the function g 2 (h i ,x) will be transformed into its equivalent function g 1 (h i ,E res ) and vice-versa. As a result, the resolution for each (h i ) of the equation g 1 (h i ,E res )=f 1 (E res ) will make it possible to determine θ ic (h i ) for each 
 
(h i )∈[h o ,R i +h o −f −1 ∘g −1 (h o )]. 
        Using the function for Linear Energy Transfer (dE/dx);        
 
         [0069]     The calculation methodology is the same as the one using the Energy-Trajectory relation, only instead of the functions Energy-Trajectory and Trajectory-Energy in the analyzed medium and the absorbent medium, they are replaced with the concept of the integral of the function for Linear Energy Transfer as described in the case without absorbent medium.  
         [0000]     II-1-2—Determination of the Probability P i (r, h).  
         [0000]     II-1-2-a—Case of Nucleus Emitting Alpha Particles of Energy E i  and for which E min ≦E ind ≦E max .  
         [0070]     According to FIG. (II- 5 ) shown in the drawings section, the proportion of particles emitted from a nucleus located in a position of coordinates (r, h) and at angles of incidence less than the critical angle of incidence (θ ic ) which reach the detector is the ratio of the surface of intersection {between the surface described by the angle (θ ic ) on the plane of the detector denoted s(θ ic ) and that of the detector denoted S(R)} and the surface described by the angle (θ ic ) on the plane of the detector denoted s(θ ic ). This proportion is described by the following equation:  
       {         S   ⁡     (     θ   ic     )       ⋂     S   ⁡     (   R   )           S   ⁡     (     θ   ic     )         }       
 
         [0071]     Thus, the probability that a particle of energy E i  emitted by the nucleus having index i contained in the volume dv and confirmed by E min ≦E ind ≦E max  can be detected by the solid state nuclear track detector is:  
         Pi   ⁡     (     r   ,   h     )       =       {       1   -     cos   ⁢           ⁢     θ   ic         2     }     ⁢     {         S   ⁡     (     θ   ic     )       ⋂     S   ⁡     (   R   )           S   ⁡     (     θ   ic     )         }           
        Evaluation of  
         {         S   ⁡     (     θ   ic     )       ⋂     S   ⁡     (   R   )           S   ⁡     (     θ   ic     )         }     ,       
 
 with S(θ ic )=πR 2   c , where R c =h·tg(θ ic ). (See FIG. II- 5 ) 
 
 a) for (θ ic )≧(θ 2 ), S(θ ic )∩S(R)=πR 2 , where R: radius of SSNTD 
 
 b) for (θ ic )≦(θ 2 ): 2 cases are present, (θ ic )≧(θ 1 ) and (θ ic )≦(θ 1 ): 
    if (θ ic )≧(θ 1 )  
           S   ⁡     (     θ   ic     )       ⋂     S   ⁡     (   R   )         =     2   ⁡     [         ∫     r   -   Rc     I     ⁢           R   c   2     -       (     x   -   r     )     2         .           ⁢     ⅆ   x         +       ∫   I   R     ⁢           R   2     -     x   2         .           ⁢     ⅆ   x           ]           
     where     
       I   =         R   2     +     r   2     -     R   c   2         2   ⁢   r           
    if (θ ic )≦(θ 1 ): S(θ ic )∩S(R)=S(θ ic ) 
 
 II-1-2-b—Case of Nuclei Emitting Alpha Particles of Energy E i  and for which E ind &gt;E max  
       
 
         [0075]     According to FIG. (II- 5 ) shown in the drawings section, the proportion of particles emitted from a nucleus located in a position of coordinates (r, h) and at angles of incidence between (θ ic1 ) and (θ ic2 ) which reach the detector is the ratio of the surface of intersection {between the surface of the ring described by the angles (θ ic1 ) and (θ ic2 ) on the plane of the detector denoted s(θ ic1 , θ ic2 ) and that of the detector denoted S(R)} and the surface of the ring described by the angles (θ ic1 ) and (θ ic2 ) on the plane of the detector denoted s(θ ic1 , θ ic2 ). This proportion is described by the following equation:  
       {         S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )       ⋂     S   ⁡     (   R   )           S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )         }       
 
         [0076]     Thus, the probability that a particle of energy E i  emitted by the nucleus having index i contained in the volume dv and confirmed by E ind &gt;E max  can be detected by the solid state nuclear track detector is:  
         Pi   ⁡     (     r   ,   h     )       =       {         cos   ⁢           ⁢     θ     ic   ⁢           ⁢   1         -     cos   ⁢           ⁢     θ     ic   ⁢           ⁢   2           2     }     ⁢     {         S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )       ⋂     S   ⁡     (   R   )           S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )         }           
       ◆   ⁢           ⁢   Evaluation   ⁢           ⁢   of   ⁢           ⁢     {         S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )       ⋂     S   ⁡     (   R   )           S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )         }         
 
 with: −S(θ ic1 ,θ ic2 )=π(R c2   2 −R c2   2 −R c1   2 ) where R c1 =h·tg(θ ic1 ), R c2 =h·tg(θ ic2 ) 
        −S(θ ic1 ,θ ic2 )∩S(R)=S(θ ic2 )∩S(R)−S(θ ic1 )∩S(R) 
 
 a) for (θ ic1 )≧(θ 2 ),S(θ ic1 ,θ ic2 )∩S(R)=πR 2 −πR 2 =0 
    where R: radius of the SSNTD 
 
 b) for (θ 1 )≦(θ ic1 )≦(θ 2 ): there are 2 cases, (θ ic2 )≧(θ 2 ) and (θ ic2 )≦(θ 2 ) 
    if (θ ic2 )≧(θ 2 ):  
           S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )       ⋂     S   ⁡     (   R   )         =       π   ⁢           ⁢     R   2       -     2   ⁡     [         ∫     r   -     R     c   ⁢           ⁢   1         I     ⁢           R     c   ⁢           ⁢   1     2     -       (     x   -   r     )     2         ·           ⁢     ⅆ   x         +       ∫     I   1     R     ⁢           R   2     -     x   2         ·           ⁢     ⅆ   x           ]             
         where   ⁢           ⁢     I   1       =         R   2     +     r   2     -     R     c   ⁢           ⁢   1     2         2   ⁢   r           
           ◆if   ⁢           ⁢     (     θ     ic   ⁢           ⁢   2       )       ≤         (     θ   2     )     ⁢     :     ⁢     S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )         ⋂     S   ⁡     (   R   )           =       S   ⁡     (     θ     ic   ⁢           ⁢   2       )       ⋂     S   ⁡     (   R   )       -     S   ⁡     (     θ     ic   ⁢           ⁢   1       )       ⋂     S   ⁡     (   R   )             
       where   ⁢     :         
           S   ⁡     (     θ     ic   ⁢           ⁢   2       )       ⋂     S   ⁡     (   R   )         =     2   ⁡     [         ∫     r   -     R     c   ⁢           ⁢   2           I   2       ⁢           R     c   ⁢           ⁢   2     2     -       (     x   -   r     )     2         ·           ⁢     ⅆ   x         +       ∫     I   2     R     ⁢           R   2     -     x   2         ·           ⁢     ⅆ   x           ]           
         I   2     =       R   2     +     r   2     -       R     c   ⁢           ⁢   2     2     ⁢           [   sic   ]           
     and     
           S   ⁡     (     θ     ic   ⁢           ⁢   1       )       ⋂     S   ⁡     (   R   )         =     2   ⁡     [         ∫     r   -     R     c   ⁢           ⁢   1           I   1       ⁢           R     c   ⁢           ⁢   1     2     -       (     x   -   r     )     2         ·           ⁢     ⅆ   x         +       ∫     I   1     R     ⁢           R   2     -     x   2         ·           ⁢     ⅆ   x           ]           
         I   1     =         R   2     +     r   2     -     R     c   ⁢           ⁢   1     2         2   ⁢   r           
 
 c) for (θ ic1 )≦(θ i ): there are 2 cases, (θ ic2 )≧(θ 1 ) and (θ ic2 )≦(θ 1 ): 
 
 1—for (θ ic2 )≧(θ 1 ): there are 2 cases, (θ ic2 )≧(θ 2 ) and (θ ic2 )≦(θ 2 ): 
    if (θ ic2 )≧(θ 2 ): S(θ ic1 , θ ic2 )∩S(R)=πR 2 −πR c1   2   
           ◆if   ⁢           ⁢     (     θ     ic   ⁢           ⁢   2       )       ≤         (     θ   2     )     ⁢     :     ⁢     S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )         ⋂     S   ⁡     (   R   )           =       2   ⁡     [         ∫     r   -     R     c   ⁢           ⁢   2           I   2       ⁢           R     c   ⁢           ⁢   2     2     -       (     x   -   r     )     2         ·           ⁢     ⅆ   x         +       ∫     I   2     R     ⁢           R   2     -     x   2         ·           ⁢     ⅆ   x           ]       -     πR     c   ⁢           ⁢   1     2           
             2   ⁢     -     ⁢           ⁢   for   ⁢           ⁢     (     θ     ic   ⁢           ⁢   2       )       ≤         (     θ   1     )     ⁢     :     ⁢     S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )         ⋂     S   ⁡     (   R   )           =       S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )       =       π   ⁢           ⁢     R     c   ⁢           ⁢   2     2       -     π   ⁢           ⁢     R     c   ⁢           ⁢   1     2             ⁢               
 
 II-2—Generalization for the Total Contribution Volume for Resolution of Tracks: 
       
 
         [0081]     According to the above, the total number of tracks observed and which are due to the radioactive nuclei of index i contained in the total volume: 
 
 n   it =2π· A   ci   te   ·∫∫P   i ( r,h )· rdrdh  
 
         [0082]     For the case E min ≦E ind ≦E max , a comparison of (θ ic ) with (θ 1 ) and (θ 2 ) allows to determine respectively r 1 (h) and r 2 (h) where these latter are characteristics of respectively the radii defining the interval of the volume where (θ ic )≦(θ 1 ) from the one where (θ ic )≧(θ 1 ) and the ones defining the interval of the volume where (θ ic )≦(θ 2 ) from the one where (θ ic )≧(θ 2 ).  
         [0083]     For the case E ind &gt;E max , a comparison of (θ ic1 ) with (θ 1 ) and (θ 2 ) allows to determine respectively r 1,1 (h) and r 2,1 (h) where these latter are characteristics of respectively the radii defining the interval of the volume where (θ ic )≦(θ 1 ) [illegible] (θ ic )≧(θ 1 ) and the ones defining the interval of the volume where (θ ic1 )≦(θ 2 ) [illegible] (θ ic1 )≧(θ 2 ). And a comparison of (θ ic2 ) with (θ 1 ) and (θ 2 ) allows us to determine respectively r 1,2 (h) and r 2,2 (h) where these latter are characteristics of respectively the radii defining the interval of the volume where (θ ic2 )≦(θ 1 ) from the one where (θ ic2 )≧(θ 1 ) and the ones defining the interval of the volume where (θ ic2 )≦(θ 2 ) from the one where (θ ic2 )≧(θ 2 ).  
         [0084]     The next part of the procedure is an integration by parts of the interval over r and h.  
         [0085]     To determine r 1 (h), r 2 (h), r 1,1 (h), r 2,1 (h), r 1,2 (h), and r 2,2 (h), it is simply a matter of solving the following equations: 
 
 It has: −(θ ic )=f(h); (θ ic1 )=f a (h); (θ ic2 )=f b (h)  
           -     (     θ   1     )       =     arc   ⁢           ⁢     cos   ⁢           [     h         h   2     +       (     R   -   r     )     2           ]         ;       (     θ   2     )     =     arc   ⁢           ⁢     cos   [     h         h   2     +       (     R   +   r     )     2           ]             
 
 thus: 
        *(θ ic )=(θ 1 )         r 1 (h); *(θ ic )=(θ 2 )         r 2 (h)     *(θ ic1 )=(θ 1 )         r 1,1 (h); *(θ ic1 )=(θ 2 )         r 2,1 (h)     *(θ ic2 )=(θ 1 )         r 1,2 (h); *(θ ic2 )=(θ 2 )         r 2,2 (h)        
 
         [0089]     Actually, the total volume of the sample contributing to the tracks (V t ) will be the total of a series of the volume V f  (V t =ΣV j ), where each V j  itself corresponds to a probability P i,j (r,h).  
         [0090]     In conclusion and in a general manner, the total number of tracks observed and which are due to the radioactive nuclei of index i contained in the total volume:  
         n   it     =       A   ci   te     ·       ∑   j     ⁢     ∫     ∫     ∫         P   ij     ⁡     (     r   ,   h     )       ·     ⅆ     v   j                       
 
 where: Pij(r,h) and vj will be defined as a function of r 1 (h), r 2 (h), r 1,1 (h), r 2,1 (h), r 1,2 (h), and r 2,2 (h) 
 
 with:  
           Pi   ⁡     (     r   ,   h     )       =       (       1   -     cos   ⁢           ⁢     θ   ic         2     )     ⁢     {         S   ⁡     (     θ   ic     )       ⋂     S   ⁡     (   R   )           S   ⁡     (     θ   ic     )         }         ;       for   ⁢           ⁢   E   ⁢           ⁢   min     ≤     E   ⁢           ⁢   ind     ≤     E   ⁢           ⁢   max           
           Pi   ⁡     (     r   ,   h     )       =       (         cos   ⁢           ⁢     θ     ic   ⁢           ⁢   1         -     cos   ⁢           ⁢     θ     ic   ⁢           ⁢   2           2     )     ⁢     {         S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )       ⋂     S   ⁡     (   R   )           S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )         }         ;       for   ⁢           ⁢   E   ⁢           ⁢   ind     &gt;     E   ⁢           ⁢   max           
 
 where:  
         {         S   ⁡     (     θ   ic     )       ⋂     S   ⁡     (   R   )           S   ⁡     (     θ   ic     )         }     ⁢           ⁢   and   ⁢           ⁢     {         S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )       ⋂     S   ⁡     (   R   )           S   ⁡     (       θ     ic   ⁢           ⁢   1       ,     θ     ic   ⁢           ⁢   2         )         }         
 
 are well defined for c [sic]
 
 III—Application Example: 
 
         [0091]     Taking the simplest case, involving the media where the total trajectory of the particle α i  of energy E i  is much less than the radius of the detector (R i &lt;&lt;R) {case most reported due to the fact that this condition is verifiable for alpha particles in almost all liquid and solid samples} and a detector having a response function for conditions of developing and reading the well-defined tracks such as E min =0 and E max &gt;E j .  
         [0092]     Following the above gives:  
         cos   ⁢           ⁢     (     θ   2     )       =     h         h   2     +       (     R   +   r     )     2               
 
 Thus, (θ 2min ) for (Ri&lt;&lt;R) approaches  
         π   2     ,       
 
 and as a result, (θ ic ) is always less than (θ 2 ) 
 
         [0093]     Actually, a comparison of (θ ic ) with (θ 1 ) allows us to determine r 1 (h) and, as a result:  
         n   it     =     2   ⁢           ⁢     π   ·     A   ci   te     ·     {         ∫   0     R       1   ′     ⁢   i         ⁢       ∫   0     (   h   )       ⁢         1   -     cos   ⁢           ⁢     θ   ic         2     ⁢   r   ⁢           ⁢     ⅆ   r     ⁢           ⁢     ⅆ   h           +       ∫   0     R   1       ⁢       ∫       r   1     ⁡     (   h   )       R     ⁢         1   -     cos   ⁢           ⁢     θ   ic         2     ⁢     (         S   ⁡     (     θ   ic     )       ⋂     S   ⁡     (   R   )           S   ⁡     (     θ   ic     )         )     ⁢   r   ⁢           ⁢     ⅆ   r     ⁢           ⁢     ⅆ   h             }             
       where   ⁢     :         
           _S   ⁢     (     θ   ic     )       =     π   ⁢           ⁢     R   c   2         ,       R   c     =     h   ·     tg   ⁡     (     θ   ic     )               
           _S   ⁢     (     θ   ic     )       ⋂     S   ⁡     (   R   )         =     2   ⁡     [         ∫     r   -     R     c   ⁢                 I     ⁢           R     c   ⁢             2     -       (     x   -   r     )     2         ·           ⁢     ⅆ   x         +       ∫   I   R     ⁢           R   2     -     x   2         ·           ⁢     ⅆ   x           ]           
       I   =         R   2     +     r   2     -     R   c   2         2   ⁢   r           
 
 IV—Solution to the Problem Involving Media Having Several Types of Nuclei Emitting Particles with Different Energies and Activities. 
 
         [0094]     A technique for extracting energy is used whereby an absorbent medium must be installed in contact with the solid state nuclear track detector in a manner such that only those particles having the largest energy can reach the detector and the other energies follow in sequence.  
         [0095]     And if it appears that for a given medium, the ability to separate the energies by the SSNTD using an absorbent medium is not feasible, then an absorbent medium of thickness (h a ) is installed in contact with the detector with different thicknesses (h a ) of the absorbent medium, in a manner such that the number of cases of thicknesses (h a ) of the absorbent medium coincides with the total number of different alpha emitters of the medium to be analyzed. As a result, one has the number of equations linking the total number of tracks observed for each (h a ) to the activities of the different alpha emitters, [illegible] of the alpha activities of the different emitters (the unknowns). Thus, [illegible] alpha activities of the different emitters.  
         [0096]     For the determination of (θ ic (h)), (θ ic1 (h)), (θ ic2 (h)) of each type of nucleus in the presence of an absorbent medium, the method is to use two equations (θ c =f(E res )) and (θ i =f(E res )) (see part B of paragraph II-1-1) and the calculation methodology is the same while taking into account the presence of the absorbing medium.  
         [0000]     V—Conclusion:  
         [0097]     This direct technique allowing the dosimetry of the radioactive nucleii in different types of samples has been developed at home and in an individual manner. It is based on the use of a novel method adapted to the experimental conditions using the solid state nuclear track detectors.  
         [0098]     It consists in calculating the probability P i (r,h) that an alpha particle emitted by a radioactive nucleus contained in the volume dv of a sample can be detected by a solid state nuclear track detector, where (dv) is the volume element of a cylinder having a radius between r and r+dr and a height between h and h+dh inside the sample to be analyzed.  
         [0099]     This method is easy to implement even though its formalism appears quite complex. It is precise, absolute, and not costly.  
         [0100]     Solely for reasons of simplification of the calculations the radius is set of the detector equal to the radius of the irradiation cell, the methodology of calculation is the same for a radius of the detector less than the radius of the irradiation cell with small modifications at the step of evaluation of the surfaces of intersection {S(θ ic )∩S(R)} and {S(θ ic1 ,θ ic2 )∩S(R)}, this time, the calculation must be made of the surfaces of intersection for r (characterizing the position of the emitting nucleus) less than the radius of the detector and for r greater than the radius of the detector.  
         [0101]     In a general manner, regardless of the appearance of the experimental function of the intrinsic critical angle of a detector θ c =F(E res ) and for the position (h i ), more than 2 intersections with θ i (h i ) are present, the efficiency of the detection for these positions (h i ) will be a series of combinations between the two general formulas for efficiencies of detection described in this paper, the one corresponding to a single intersection and the other corresponding to two intersections, where the nature of this combination depends on the number of intersections.  
       REFERENCES  
       [0000]    
       
          [1] G. Somogyi, Nucl. Instr. Meth., 173, 21 (1980).  
          [2] J. F. Ziegler, Programme TRIM.