Skip to content

Commit 87806fa

Browse files
author
Ilya Mandel
committed
Plot mod
1 parent b48aa8d commit 87806fa

File tree

1 file changed

+5
-5
lines changed

1 file changed

+5
-5
lines changed

compas_matlab_utils/ComparisonPlots.m

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -85,8 +85,8 @@ function ComparisonPlots(filename1, name1, filename2, name2)
8585
DWDplot(filename2, name2, 6, 'b', 20);
8686
end;
8787
figure(6); hold off; set(gca,'FontSize',20); title('Double White Dwarfs'); legend;
88-
xlabel('$M*a [M_\odot * R_\odot]$ @ ZAMS', 'Interpreter', 'latex');
89-
ylabel('$M*a [M_\odot * R_\odot]$ @ end', 'Interpreter', 'latex');
88+
xlabel('$log_{10}(M*a) [M_\odot * R_\odot]$ @ ZAMS', 'Interpreter', 'latex');
89+
ylabel('$log_{10}(M*a) [M_\odot * R_\odot]$ @ end', 'Interpreter', 'latex');
9090

9191

9292
[binariescount1, SNcount1, BHcompletecount1, SNbothcount1, SNonecount1, ...
@@ -189,7 +189,7 @@ function ComparisonPlots(filename1, name1, filename2, name2)
189189
'DisplayName', ['CE, ', name]); hold on;
190190
scatter(log10(P(BNS & ~isCE)), e(BNS & ~isCE), point, colour, 'DisplayName', ['Stable, ', name]);
191191

192-
%%%
192+
%%% TO CLEAN
193193
good = mergingBBH & isCE & OKCE;
194194
type1CE = h5read(file,'/BSE_Common_Envelopes/Stellar_Type(1)<CE');
195195
type2CE = h5read(file,'/BSE_Common_Envelopes/Stellar_Type(2)<CE');
@@ -330,8 +330,8 @@ function DWDplot(file, name, fignumberDWD, colour, point)
330330
SeedRLOF=h5read(file, '/BSE_RLOF/SEED');
331331
[hadRLOF,RLOFIndex]=ismember(ind,SeedRLOF);
332332
figure(fignumberDWD), hold on;
333-
scatter(MAZAMS(~hadRLOF), MAWD(~hadRLOF), point, 'filled', colour, 'DisplayName', ['DWDs without mass transfer, ', name]);
334-
scatter(MAZAMS(hadRLOF), MAWD(hadRLOF), point, colour, 'DisplayName', ['DWDs after mass transfer, ', name]);
333+
scatter(log10(MAZAMS(~hadRLOF)), log10(MAWD(~hadRLOF)), point, 'filled', colour, 'DisplayName', ['DWDs without mass transfer, ', name]);
334+
scatter(log10(MAZAMS(hadRLOF)), log10(MAWD(hadRLOF)), point, colour, 'DisplayName', ['DWDs after mass transfer, ', name]);
335335
end %end of DWDplot
336336

337337
%SN varieties -- just count

0 commit comments

Comments
 (0)