plotting_graphMetrics.m

function plotting_graphMetrics(G, metrics, grainsVar, mergedGrains, connectedGrains)
    % Create comprehensive twin analysis plots
    % Inputs:
    %   G - graph object from twin network
    %   metrics - structure from calculateTwinMetrics
    %   grainsVar - grains or parentGrains variable

    imgSampPath = evalin('base', 'imgSampPath');
    samp_name = evalin('base', 'samp_name');
    methodName = evalin('base', 'methodName');

    phaseName = evalin('base', 'phaseName');
    if length(grainsVar.mineralList) > 1
        grainsVar = grainsVar(phaseName);
    end
    
    % Twin order distribution with centrality correlation
    plot_TwinOrder(G, metrics, imgSampPath, samp_name, methodName);
    
    % Max twin depth plots
    plot_MaxTwinDepth(metrics, grainsVar, mergedGrains, imgSampPath, samp_name, methodName);
    
    % Twin variant distribution
    plot_TwinVariants(metrics, imgSampPath, samp_name, methodName);
    
    % Twin degree connectivity
    plot_TwinDegree(metrics, grainsVar, mergedGrains, imgSampPath, samp_name, methodName);
    
    % Connected components
    plot_ConnComponents(metrics, grainsVar, connectedGrains, imgSampPath, samp_name, methodName);

end