# plotting carrier (electron) lifetime

v_index=1; 
x =  -0.025e-6;   
z = linspace(getnamed('silicon','z min'),getnamed('cathode','z min'),401);
recomb = getresult('CHARGE','recombination');
vtx = [recomb.x, recomb.z];   
taun = -log10(pinch(recomb.taun,3,1));  
taun_src = interptri(recomb.elements,vtx,taun,x,z,0); 
plot(z*1e6,taun_src,'depth (um)','-log10(taun) (sec)');


# plotting Trap assisted recombination rate

z = linspace(-0.28e-6,-0.1e-6,101);
Rsrh_3 = abs(pinch(recomb.Rsrh,3,3));
Rsrh_3_src = interptri(recomb.elements,vtx,Rsrh_3,x,z,0);
Rsrh_4 = abs(pinch(recomb.Rsrh,3,4));
Rsrh_4_src = interptri(recomb.elements,vtx,Rsrh_4,x,z,0);
Rsrh_5 = abs(pinch(recomb.Rsrh,3,5));
Rsrh_5_src = interptri(recomb.elements,vtx,Rsrh_5,x,z,0);
Rsrh_6 = abs(pinch(recomb.Rsrh,3,6));
Rsrh_6_src = interptri(recomb.elements,vtx,Rsrh_6,x,z,0);
plot(z,Rsrh_3_src,Rsrh_4_src,Rsrh_5_src,Rsrh_6_src,'z (um)','log(abs(Rsrh))','','logplot'); 
legend('0.8 V','1.2 V', '1.6 V', '2 V');