import sys

# Read the name of the file containing the TransmissionSpectrum object
filename1 = sys.argv[1]

# Read the name of the file containing the InelasticTransmissionSpectrum object
filename2 = sys.argv[2]

# Read the TransmissionSpectrum object
transmission = nlread(filename1,TransmissionSpectrum)[-1]

# Read InelasticTransmissionSpectrum object
transmission_xloe = nlread(filename2,InelasticTransmissionSpectrum)[-1]

# Retrieve the bias voltage at which the calculation has been performed
bias = transmission.bias()

# Calculate the elastic contribution to the current
current_elastic  = transmission_xloe.elasticCurrent(bias=bias)

# Calculate the inelastic contribution to the current
current_inelastic  = transmission_xloe.inelasticCurrent(bias=bias)

# Calculate the current without interactions
current_wointeractions = current_elastic

# Calculate the current with interactions
current_winteractions = current_elastic + current_inelastic

# Print the results
print('--------------------------------------------------------------------------------')
print(' TransmissionSpectrum read from file :')
print(' %s' %(filename1))
print(' InelasticTransmissionSpectrum read from file :')
print(' %s' %(filename2))
print('--------------------------------------------------------------------------------')
print(' Current (w/o interactions)    = %.2e nA' %(current_wointeractions*10**9))
print(' Current (w interactions)      = %.2e nA' %(current_winteractions*10**9))
print('--------------------------------------------------------------------------------')