PRO dustem_run_plugins, p_dim ,$
param_descs,$
param_values,$
param_func,$
scopes,$
st=st,$
avoid=avoid,$
dustem_run=dustem_run,$
help=help

;+
; NAME:
;    dustem_run_plugins
; PURPOSE:
;    runs the plugins according to their scopes
;    a N_ELEMENTS(PLUGINS)X N_ELEMENTS(SCOPES) matrix is created to achieve this 
; CATEGORY:
;    DUSTEM Wrapper, Distributed, Mid-Level, Plugins
; CALLING SEQUENCE:
;    dustem_run_plugins,p_dim,param_descs,param_values,scopes,[,st][,avoid][,dustem_run][,/help] 
; INPUTS:
;    p_dim  = parameter values 
;    param_descs = parameter description vector
;    param_values = current parameter values
;    param-func   = plugin indices array 
;    scopes       = scopes of plugins 
; OPTIONAL INPUT PARAMETERS:
;    avoid      = scopes should be avoided
;    dustem_run = Dustem should be run
; OUTPUTS:
;    None
; OPTIONAL OUTPUT PARAMETERS:
;    st         = Dustem output structure
; ACCEPTED KEY-WORDS:
;    help                  = if set, print this help
; COMMON BLOCKS:
;    None
; SIDE EFFECTS:
;    None
; RESTRICTIONS:
;    The DustEM fortran code must be installed
;    The DustEMWrap idl code must be installed
; PROCEDURES AND SUBROUTINES USED:
; EXAMPLES
;    dustem_activate_plugins, p_min, pram_descs, param_values, param_func, scopes, st=st, avoid=avoid, dustem_run=dustem_run 
; MODIFICATION HISTORY:
;    Written by IC 
;    Evolution details on the DustEMWrap gitlab.
;    See http://dustemwrap.irap.omp.eu/ for FAQ and help.  
;-



;Because I don't want to run this block each time the procedure is run I will test on the presence of all the scope tags in !dustem_plugin
;if all of them are present then this block has been run. 

IF (tag_names(*!dustem_plugin))(0) EQ 'NONE' THEN goto, the_end 

;Just to be sure 
if not keyword_set(avoid) then avoid=0
if not keyword_set(dustem_run) then dustem_run=0


test_all_scopes = 1. 
FOR i=0L,n_tags(*!dustem_plugin)-1 DO BEGIN

	IF isa(((*!dustem_plugin).(i).scope)) then test=1 else test=0
	test_all_scopes*=test 

ENDFOR

IF ~test_all_scopes THEN BEGIN
	;FIRST LOOP TO RETRIEVE THE SCOPES OF ALL THE PLUGINS
	f=1 ; Initializing the index that is associated to each plugin 
	FOR i=0L,n_elements(param_descs)-1 DO BEGIN
	   
	   parameter_type=dustem_parameter_description2type(param_descs[i],string_name=string_name) ; Looping over the parameter description vector (here saved in a system variable) to sort out its different elements with repect to their types
	   
	   IF parameter_type EQ 'PLUGIN' THEN BEGIN   ; Selecting the plugins
	        
	        ftn = strmid(param_descs(i),0) ; String containing the name of the plugin and the keyword used (ie: dustem_create_continuum_2)    
	        ii = strsplit(ftn,'_',count=countx) & ii = ii(countx-1)-1 ; Locating the last underscore to automate the extraction of the plugin's keyword
	        ftn = strmid(ftn,0,ii) ; String containing the name of the plugin without the associated keyword
	        
	        k=where(strmid(tag_names(*!dustem_plugin),0,8) eq strmid(strupcase(strmid(ftn,14)),0,8),counte) ; Selecting a plugin through matching the string name of the plugin from the scope system variable with the one read from the parameter description vector
	        
	        tmp = where(param_func eq f, count) ; Array of indices helping with the sorting out of the different plugins (as a plugin can be called several times with different keywords so as to fit several quantities) 
	        PDO_tmp = param_descs(tmp) ; Array containing the different calls of each plugin at a time
	         
	        p_dim_tmp = param_values(tmp) ; Array containing the values of the fitted parameters
	        
	        vari = strmid(PDO_tmp(0),ii+1) ; Test variable to be used to alter the initialization of the index and value arrays below as a keyword can be a scalar or a string
	        
	        ; Initialization of the index and value arrays. These arrays will be used (below) when calling the plugins (with the default key=key and val=val keywords)
	        
	        if strlen(strtrim(string(vari),2)) gt 1 then begin
	            index = strarr(count) & value = fltarr(count)  
	        endif else  begin
	            index = fltarr(count) & value = fltarr(count)
	        endelse
	        
	        ; Filling the index and value arrays for one plugin at a time  
	        FOR l=0, count-1 DO BEGIN
	           index[l] = strmid(PDO_tmp[l],ii+1) 
	           value[l] = p_dim_tmp[l]                      
	        ENDFOR    
	        
	        ;==============Dry run of the plugins to obtain their scopes and run them accordingly==============
	        
	        str='toto='+ftn+'(scope=scope)' & str=str(0)        
	        toto=execute(str)            
	        ((*!dustem_plugin).(k).scope)=ptr_new(scope)
	           
	        str='toto='+ftn+'(paramtag=paramtag)' & str=str(0)        
            toto=execute(str)            
            ((*!dustem_plugin).(k).paramtag)=ptr_new(paramtag)    
	            
	            
	        f=f+1 & i=i+count-1 ; Incrementing the parameter and same-type plugin indices  
	  ENDIF
	ENDFOR
ENDIF


;LOCATING THE PLUGINS THAT SHOULD BE RUN

indices = fltarr(n_tags(*!dustem_plugin),n_elements(scopes)) ;


FOR i=0L,n_elements(scopes)-1 DO BEGIN

	;test if there is a '*' string. This means that all scopes containing that string will be either avoided or chosen.
	;this means that we need the index of the scope corresponding to this string.
		
	;first we need to extract the string between the two '*' and then match it with the existing plugin 
	; this enables us to get the indices of the plugins to run...

	IF strmid(scopes(i),0,1) EQ '*' THEN BEGIN

		;Now I need to locate the second '*': 
		ind_ast = STRPOS(scopes(i), '*', /REVERSE_SEARCH)
		scope_to_match = strmid(scopes(i),1,ind_ast-1) 
		;
	ENDIF ELSE BEGIN ;assuming the string is identical to the plugin's scope
	
    	scope_to_match = strupcase(scopes(i))  
	
	ENDELSE	
	;Now locating the indices of this scope in !dustem_plugin. 
    
	FOR j=0L,n_tags(*!dustem_plugin)-1 DO BEGIN
        
        test = strcmp((*((*!dustem_plugin).(j).scope)),scope_to_match,strlen(scope_to_match),/FOLD_CASE) EQ 1
        
        IF test and keyword_set(avoid) THEN BEGIN
            if (size(indices))[0] GT 1 then indices[j,i] *= 0 ELSE indices[j] *= 0
 		ENDIF ;ELSE BEGIN
;     		if (size(indices))[0] GT 1 then indices[j,i] = 1 ELSE indices[j] = 1
; 		ENDELSE 
; 		
		IF test and not keyword_set(avoid)then  BEGIN
    		if (size(indices))[0] GT 1 then indices[j,i] += 1 ELSE indices[j] += 1
    		
		ENDIF
		
		if ~test and keyword_set(avoid) then begin
            if (size(indices))[0] GT 1 then indices[j,i] += 1 ELSE indices[j] += 1

        endif
        		
	ENDFOR		

ENDFOR


;RUNNING THE PLUGINS ACCORDING TO THEIR SCOPES

f=1 ; Initializing the index that is associated to each plugin 
FOR i=0L,n_elements(param_descs)-1 DO BEGIN
	   
	 parameter_type=dustem_parameter_description2type(param_descs[i],string_name=string_name) ; Looping over the parameter description vector (here saved in a system variable) to sort out its different elements with repect to their types
	   
	  IF parameter_type EQ 'PLUGIN' THEN BEGIN   ; Selecting the plugins
	        
	       ftn = strmid(param_descs(i),0) ; String containing the name of the plugin and the keyword used (ie: dustem_create_continuum_2)    
	       ii = strsplit(ftn,'_',count=countx) & ii = ii(countx-1)-1 ; Locating the last underscore to automate the extraction of the plugin's keyword
	       ftn = strmid(ftn,0,ii) ; String containing the name of the plugin without the associated keyword
	        
	       k=where(strmid(tag_names(*!dustem_plugin),0,8) eq strmid(strupcase(strmid(ftn,14)),0,8),counte) ; Selecting a plugin through matching the string name of the plugin from the scope system variable with the one read from the parameter description vector
	        
	       tmp = where(param_func eq f, count) ; Array of indices helping with the sorting out of the different plugins (as a plugin can be called several times with different keywords so as to fit several quantities) 
	       PDO_tmp = param_descs(tmp) ; Array containing the different calls of each plugin at a time
	         
	       p_dim_tmp = param_values(tmp) ; Array containing the values of the parameters used for the run
	        
	       vari = strmid(PDO_tmp(0),ii+1) ; Test variable to be used to alter the initialization of the index and value arrays below for a keyword can be a scalar or a string
	        
	        ; Initialization of the index and value arrays. These arrays will be used (below) when calling the plugins (with the default key=key and val=val keywords)
	        
	       if strlen(strtrim(string(vari),2)) gt 1 then begin
	           index = strarr(count) & value = fltarr(count)  
	       endif else  begin
	           index = fltarr(count) & value = fltarr(count)
	       endelse
	        
	        ; Filling the index and value arrays for one plugin at a time  
	       FOR l=0, count-1 DO BEGIN
	          index[l] = strmid(PDO_tmp[l],ii+1) 
	          value[l] = p_dim_tmp[l]                      
	       ENDFOR    
	        
	       ;==============Runnig of the plugins according to their scopes==============
	   
	       
	       if keyword_set(avoid) then idd=product(indices[k,*]) else idd=total(indices[k,*])
            	   
	       IF  idd GE 1 THEN BEGIN
	           
               str='((*!dustem_plugin).('+strtrim(k,2)+').spec)=ptr_new('+ftn+'(key=index,val=value)'+')' & str=str(0)
               toto=execute(str) & IF !dustem_verbose NE 0 THEN message,strupcase(strmid(ftn,7)),/info
               
	      ENDIF
    	   
    	     
	          
	       f=f+1 & i=i+count-1 ; Incrementing the parameter and same-type plugin indices  
    	  
	  ENDIF
	ENDFOR
	
if keyword_set(dustem_run) then begin
    
    st = dustem_run(p_dim) 
    !dustem_current = ptr_new(st)
    
endif 
	
the_end:
END