ó
ô’§Qc           @   sÚ   d  Z  d GHd d l Z d e j d <e j d d d ƒe j d	 d
 e ƒd d l Td d l Td d l Td d l m	 Z	 d d l
 Z
 d d l m Z d d l Z d GHd „  Z d „  Z d „  Z e d k rÖ e ƒ  n  d S(   s.    CREATED: 7/2012
    AUTHOR: MICHAL ODSTRÄŒIL
s   importing modules iÿÿÿÿNt   Aggt   backendt   fontt   sizet   10t   textt   usetex(   t   *(   t   Shot(   t
   mquantiless   importing modules done c         C   s  t  ƒ  } | j } t d d^ d |  ƒ } | d } | rN | d } | d } n d } d	 } t j ƒ  } d } d
 }	 | j d ƒ rÁ | d \ }
 } | |	 9} t | |
 | k |
 | k  @d ƒ } n  d } xa d d d g D]P } y@ | | \ } } t | t | | | k | | k  @d ƒ ƒ } Wq× q× Xq× Wd t | d | d ƒ g } d t | ƒ } t d d d d d d  g ƒt d d d d d d  g d t
 ƒ| r¸t d d d d d d t
 ƒn t d  d! d" d d d t
 ƒt d d# d$ d | ƒt d d% d$ d | ƒg t d d& d$ d | ƒg | j d ƒ r\| d' d( k r\t d d) d$ d | d |	 d t
 ƒn d  g | r°| j d* ƒ r°| d+ d, k r°| d- d. k  r°t d* d/ d0 d d1 ƒn d  g g } d2 GHt | | d3 |  d4 GHd5 d6 l m } m } | d7 | ƒ  ƒ GHt d d d ƒ t d d d d t
 ƒt d  d8 d" d d d t
 ƒt d d# d$ ƒ t d d9 d$ ƒ g g } t | | d: |  t d d d d; d d  g ƒt d d d d; d d  g d t
 ƒt d d d< d; d d  g d d d t
 ƒt d= d! d< d; d d  g d d d t
 ƒt d  d8 d< d; d d  g d d d t
 ƒg t d d# d$ d; d d  g ƒt d d9 d$ d; d d  g ƒg g } t | | d> d? |  t d d d d; d d  g ƒt d d d d; d d  g d t
 ƒ| rt d  d8 d@ d; d d  g d d d t
 ƒt d d d d; d d  g d d d t
 ƒg n' t d  d8 d" d; d d  g d d d t
 ƒt dA d# dB d; d d  g ƒt dC d9 dB d; d d  g ƒg g } t | | dD dE |  t dF d dB d; d d  g ƒt dG dH dI d; d d  g ƒt dJ dK dL d; d d  g ƒt dM dN dO d; d d  g ƒg t dA d# dB d; d d  g ƒt dC d9 dB d; d d  g ƒg g } t | | dP dQ |  t d dR dR dS dR ƒ| rªt d  dR dR dS dR ƒg t d dR dR dS dR ƒg n t d  dR dR dS dR ƒg g } t | dR dT d_ dW ƒ t d dR d d; d d  g d d d  g dS dR ƒ| r‰t d  dR d@ d; d d  g d d d  g dS dR d d d t
 ƒt d dR d d; d d  g d d d  g dS dR d d d t
 ƒg n9 t d  dR d" d; d d  g d d d  g dS dR d d d t
 ƒg } t | dR dX d` dZ d[ ƒ d\ t j ƒ  | GHt d] d ƒ d  S(a   Nt   figsizei	   i   t	   file_typet   plasmat   plasma_startt
   plasma_endi    g{®Gáz¤?gš™™™™™™?t
   hxr_smoothg®Gáz®ï?g333333Ó?t
   photodiodet   photodiode_alphat   photodiode_otherg333333ó?s   Golem shot No:t   loop_voltages   Loop voltages   U [V]t   ylimt   toroidal_fields   Toroidal mag. fields	   B$_t$ [T]t	   reductiont   plasma_currents   Plasma currents   I$_{p}$ [kA]t   data_rescalegü©ñÒMbP?t   rogowski_currents   Chamber currents   I$_{ch}$ [kA]t   Visibles   Intensity [a.u.]s
   H$_\alpha$t   Othert   hxr_meang{®Gáz„?t   HXRt   electron_densityt   electron_density_meang €à7yÃACs   electron_density:reliabilityg      Ð?s   Electron density ($n_e$)s    n$_{e}$ [10$^{19}\cdot$m$^{-3}$]g¬Ò¶OÉƒý;s   start plottingt	   graphpress   end plottingiÿÿÿÿ(   t   gmtimet   strftimes   %Y-%m-%d %H:%M:%Ss   Total currents	   H$_\alphat
   graphprintt   xlims   I [kA]t   chamber_currents    - final data t   graphpresfulls   I$_p$+I$_{ch}$ [kA]t   photos   [DAS V]t   haphotos    - integrated datat
   graphpresit   uloopt   btors   Derivative  of mag. fields   dB$_t$/dt [DAS V]t   irogs   Raw Rogowski signals   dI$_{p+ch}$/dt [DAS V]t   PlasmaDetects   Derivative of I$_{p}$s   dI$_{p}$/dt [a.u.]s    - raw datat
   graphpresbt    t   xlabelt   iconi   i   i(   t   graphicg      ø?id   t
   horizontals	    Time: %gt   status(   i	   i   (   i   i   (   i	   g      ø?(   R   t   shot_numt   dictt   timet   existR	   t   maxt   strt   get_datat   Nonet   Truet	   multiplotR!   R"   t   paralel_multiplott	   saveconst(   R   t   St   shott   plot_paramsR   t   startt   endt   t0t   max_HXRt	   HXR_constt   tvec_HXRR   t	   max_photot   sigt   tvecR   t   ph_ranget
   shot_titlet   dataR!   R"   (    (    s8   ../includes/analysis/Basics/0411ShotHomepage.ONN/main.pyt	   plot_data   s”    		


&0!$FJY!''*$!Z'$$:-½c    3   	   B   sÚ  d d l  Te j ƒ  }  e d ƒ } e ƒ  e ƒ  e ƒ  \ } } } } e ƒ  \ } } } }	 e |	 ƒ \ }
 } } } e | |
 | ƒ \ } } e	 | | ƒ \ } } } } | rÌ e | |
 | | | ƒ \ } } n  xˆ e
 d d d g d d d	 g ƒ D]e \ } } y, e | | | | ƒ } e | | | | ƒ Wqî e k
 rR} d
 | d Ge | ƒ GHqî Xqî We | | | ƒ } e | | | ƒ } e | | | ƒ } e | | | ƒ } e | | ƒ } e | | | | | | | | ƒ } | r;e | | ƒ } e | | ƒ } e | | | | | ƒ \ }  }! }" }# e | ƒ }$ e | | |$ | ƒ }% e | | | | ƒ \ }& }' }( e | ƒ }) e | | | | ƒ } e d d d d d d g ƒ d }* e |* ƒ }+ e |+ ƒ }, xT e  |+ ƒ D]F }- e! | | | | |* |- ƒ \ }. |, |- <d |* |- |. |, |- f GHq¾We" |, e# |, ƒ <e! | | | | |* e$ |, ƒ ƒ n  e% | ƒ }/ e& | | | ƒ \ } }0 e' | ƒ }1 e( | | | ƒ }2 e) | |
 | | | ƒ y e* ƒ  Wn" e k
 rÂ} d Ge | ƒ GHn Xd Ge j ƒ  |  GHd  S(   Niÿÿÿÿ(   R   t   Aktual_PfeifferMerkaVakuat   photovrbR(   R'   t
   PhotodVrbat   PhotodHalphat   Photods   Photodiode s    failed gš™™™™™©?gš™™™™™¹?g333333Ã?gš™™™™™É?g333333Ó?gš™™™™™Ù?gü©ñÒMbP?s   win %g rate %g err %gs   BreakdownProba faileds   time of basic diagn generation(+   t   basic_diagnR7   t	   loadconstt   getDatet   save_configt   getBtoroidalt   getUloopt   getIrogowskit
   getIplasmaR-   t   zipt	   getPhotodt   getMeanPhotodt	   ExceptionR:   t	   getMeanBtt   getTotalCharget   getMeanUloopt   getMeanCurrentt   getOhmicHeatingPowert   Failurest   getQedget   getMeanElectronTemperaturet   getBreakDownVoltaget   getStateEqElectronDensityt%   getElectronConfinementTimeFirstApproxt   getElectronTemperaturet   getGreenwaldDensityt   getQedgeTimet   arrayt   lent   zerost   ranget   getBreakDownRatet   inft   isnant   argmint   getChamberResistancet   getOhmicHeatingPowerTimet   getMagneticFluxt   getTransformatorSaturationt   EnergyBalancet   BreakdownProba(3   t   tRQ   t   Btort   dBtort   BtMaxt   BtMeant   Uloopt   UloopMaxt	   UloopMeant
   ReversedCDt   Irogt   dIdt_rogMaxt   IrogMaxt   I_startt   Iplat   Icht   Plasmat   PlasmaStartt	   PlasmaEndt   PlasmaTimeLengtht   namet   fnameR   t   et   MeanBtt   TotalCharget	   MeanUloopt   MeanIplat   OhmicHeatingPowert   PlasmaStatust   Qedget   ElectronTemperaturet   Umaxt   Btimet   Bbreakt   Biplat   StateEqElectronDensityt"   ElectronConfinementTimeFirstApproxt   ElectronTempTimet   medElectronTempt   maxTt   GreenwaldDensityt   trangest   ntt   errt   it   ratet   ChamberResistancet   OhmicHeatingChambert   MagneticFluxt
   Saturation(    (    s8   ../includes/analysis/Basics/0411ShotHomepage.ONN/main.pyt   prepare_data¨   sf    
!.!$"&#c           C   s…   t  j d d k r t ƒ  nd t  j d d k rJ t d ƒ t d d ƒ n7 t  j d d k r t d ƒ t d ƒ t d d ƒ n  d  S(	   Ni   t   acquisitiont   plotst   pngR4   i    t   postanalysist   svgz(   t   syst   argvR¯   RP   R@   (    (    (    s8   ../includes/analysis/Basics/0411ShotHomepage.ONN/main.pyt   mainó   s    



t   __main__(   t   __doc__t
   matplotlibt   rcParamst   rcR=   t   numpyt   pygolem_lite.configt   pygolem_lite.modulest   pygolem_liteR   Rµ   t   scipy.stats.mstatsR	   R7   RP   R¯   R·   t   __name__(    (    (    s8   ../includes/analysis/Basics/0411ShotHomepage.ONN/main.pyt   <module>   s$   


	Œ	K	