<jupyter_start><jupyter_text>Compare Composite Lineshapes<jupyter_code>#Transitional Disks

#CI Tau data
citau_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/data/transitional_disks/citau_lineflux_data.p','rb'))
citau_lineshape_data=make_lineshape(citau_wave, citau_flux, citau_lineflux_data)

#LkHa 330 data
lkha330_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/data/transitional_disks/lkha330_lineflux_data.p','rb'))
lkha330_lineshape_data=make_lineshape(lkha330_wave, lkha330_flux, lkha330_lineflux_data)

#TW Hya data
twhya_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/data/transitional_disks/twhya_lineflux_data.p','rb'))
twhya_lineshape_data=make_lineshape(twhya_wave, twhya_flux, twhya_lineflux_data)

#DoAr 44 data
doar44_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/data/transitional_disks/doar44_lineflux_data.p','rb'))
doar44_lineshape_data=make_lineshape(doar44_wave, doar44_flux, doar44_lineflux_data)

#HD 135344 data
hd135344_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/data/transitional_disks/hd135344_lineflux_data.p','rb'))
hd135344_lineshape_data=make_lineshape(hd135344_wave, hd135344_flux, hd135344_lineflux_data)

#UX Tau data
uxtau_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/data/transitional_disks/uxtau_lineflux_data.p','rb'))
uxtau_lineshape_data=make_lineshape(uxtau_wave, uxtau_flux, uxtau_lineflux_data)

#Classical Disks

#DoAr 24 data
doar24_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/data/classical_disks/doar24_lineflux_data.p','rb'))
doar24_lineshape_data=make_lineshape(doar24_wave, doar24_flux, doar24_lineflux_data)

#DF Tau data
dftau_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/data/classical_disks/dftau_lineflux_data.p','rb'))
dftau_lineshape_data=make_lineshape(dftau_wave, dftau_flux, dftau_lineflux_data)

#DL Tau data
dltau_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/data/classical_disks/dltau_lineflux_data.p','rb'))
dltau_lineshape_data=make_lineshape(dltau_wave, dltau_flux, dltau_lineflux_data)


#Selina's Disks

# Classical Disks

#GK Tau data
gktau_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/classical/gktau_line.p','rb'))
gktau_lineshape_data=make_lineshape(gktau_wave, gktau_flux, gktau_lineflux_data)

#RNO 90 data
rno90_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/classical/rno90_line.p','rb'))
rno90_lineshape_data=make_lineshape(rno90_wave, rno90_flux, rno90_lineflux_data)


#Herbig Disks

#AB Aur data
abaur_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/herbig/abaur_line.p','rb'))
abaur_lineshape_data=make_lineshape(abaur_wave, abaur_flux, abaur_lineflux_data)

#HD 36112 data
hd36112_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/herbig/hd36112_line.p','rb'))
hd36112_lineshape_data=make_lineshape(hd36112_wave, hd36112_flux, hd36112_lineflux_data)

#HD 144432 data
hd144432_lineflux_data=pickle.load(open('/Users/erichegonzales/Desktop/eriche-thesis/herbig/hd144432_line.p','rb'))
hd144432_lineshape_data=make_lineshape(hd144432_wave, hd144432_flux, hd144432_lineflux_data)
#Plotting composite lineshapes 0
fig=plt.figure(figsize=(16,8))

ax1=fig.add_subplot(121)
ax1.plot(citau_lineshape_data[0], citau_lineshape_data[1], label='Ci Tau')
ax1.plot(lkha330_lineshape_data[0], lkha330_lineshape_data[1], label='DoAr 44')
ax1.plot(twhya_lineshape_data[0], twhya_lineshape_data[1], label='TW Hya')
ax1.plot(doar44_lineshape_data[0], doar44_lineshape_data[1], label='DoAr 44')
ax1.plot(hd135344_lineshape_data[0], hd135344_lineshape_data[1], label='HD 135344')
ax1.plot(uxtau_lineshape_data[0], uxtau_lineshape_data[1], label='UX Tau')

#Setting labels, limits, legend
ax1.set_xlabel('Velocity [km/s]')
ax1.set_ylabel('Arbitrary flux')
ax1.set_title('Composite Lineshape of Transitional Disks')
ax1.legend(loc="upper center", bbox_to_anchor=(0.5, 1.2), ncol=2)
ax1.set_ylim(0, 4)

#Plotting composite lineshapes
ax2=fig.add_subplot(122)
ax2.plot(doar24_lineshape_data[0], doar24_lineshape_data[1], label='DoAr 24')
ax2.plot(dftau_lineshape_data[0], dftau_lineshape_data[1], label='DF Tau')
ax2.plot(dltau_lineshape_data[0], dltau_lineshape_data[1], label='DL Tau')
ax2.plot(gktau_lineshape_data[0], gktau_lineshape_data[1], label='GK Tau')
ax2.plot(rno90_lineshape_data[0], rno90_lineshape_data[1], label='RNO 90')

#Setting labels, limits, legend
ax2.set_xlabel('Velocity [km/s]')
ax2.set_ylabel('Arbitrary flux')
ax2.set_title('Composite Lineshape of Classical Disks')
ax2.legend(loc="upper center", bbox_to_anchor=(0.5, 1.2), ncol=2)
ax2.set_ylim(0, 4)
#Plotting composite lineshapes 1
fig=plt.figure(figsize=(16,8))

ax1=fig.add_subplot(121)
ax1.plot(citau_lineshape_data[0], citau_lineshape_data[1]-1.55, label='Ci Tau')
ax1.plot(lkha330_lineshape_data[0], lkha330_lineshape_data[1]-1.03, label='DoAr 44')
ax1.plot(twhya_lineshape_data[0], twhya_lineshape_data[1]-0.18, label='TW Hya')
ax1.plot(doar44_lineshape_data[0], doar44_lineshape_data[1]-0.43, label='DoAr 44')
ax1.plot(hd135344_lineshape_data[0], hd135344_lineshape_data[1]-3.18, label='HD 135344')
ax1.plot(uxtau_lineshape_data[0], uxtau_lineshape_data[1]-0.33, label='UX Tau')

#Setting labels, limits, legend
ax1.set_xlabel('Velocity [km/s]')
ax1.set_ylabel('Arbitrary flux')
ax1.set_title('Composite Lineshape of Transitional Disks')
ax1.legend(loc="upper center", bbox_to_anchor=(0.5, 1.2), ncol=2)
ax1.set_ylim(0, 2)
ax1.set_xlim(-150, 150)
#Plotting composite lineshapes 2
fig=plt.figure(figsize=(16,8))

ax1=fig.add_subplot(121)
ax1.plot(citau_lineshape_data[0], (citau_lineshape_data[1]-1.61)*2.17, label='Ci Tau')
ax1.plot(lkha330_lineshape_data[0], (lkha330_lineshape_data[1]-1.1)*5.6, label='LkHa 330')
ax1.plot(twhya_lineshape_data[0], (twhya_lineshape_data[1]-0.242)*6.5, label='TW Hya')
ax1.plot(doar44_lineshape_data[0], (doar44_lineshape_data[1]-0.505)*7.8, label='DoAr 44')
ax1.plot(hd135344_lineshape_data[0], (hd135344_lineshape_data[1]-3.23)*3.15, label='HD 135344')
ax1.plot(uxtau_lineshape_data[0], (uxtau_lineshape_data[1]-0.416)*54, label='UX Tau')

#Setting labels, limits, legend
ax1.set_xlabel('Velocity [km/s]')
ax1.set_ylabel('Arbitrary flux')
ax1.set_title('Normalized Composite Lineshape of Transitional Disks')
ax1.legend(loc="upper center", bbox_to_anchor=(0.5, 1.2), ncol=2)
ax1.set_ylim(0, 1)
ax1.set_xlim(-150, 150)
#Plotting composite lineshapes 3
fig=plt.figure(figsize=(16,8))

#Plotting composite lineshapes
ax2=fig.add_subplot(122)
ax2.plot(doar24_lineshape_data[0], doar24_lineshape_data[1]-2.02, label='DoAr 24')
ax2.plot(dftau_lineshape_data[0], dftau_lineshape_data[1]-1.4, label='DF Tau')
ax2.plot(dltau_lineshape_data[0], dltau_lineshape_data[1]-0.85, label='DL Tau')
ax2.plot(gktau_lineshape_data[0], gktau_lineshape_data[1]-0.58, label='GK Tau')
ax2.plot(rno90_lineshape_data[0], rno90_lineshape_data[1]-2.32, label='RNO 90')

#Setting labels, limits, legend
ax2.set_xlabel('Velocity [km/s]')
ax2.set_ylabel('Arbitrary flux')
ax2.set_title('Composite Lineshape of Classical Disks')
ax2.legend(loc="upper center", bbox_to_anchor=(0.5, 1.2), ncol=2)
ax2.set_ylim(0, 2)
ax2.set_xlim(-150, 150)
#Plotting composite lineshapes 4
fig=plt.figure(figsize=(16,8))

#Plotting composite lineshapes
ax2=fig.add_subplot(122)
ax2.plot(doar24_lineshape_data[0], (doar24_lineshape_data[1]-2.02)*1.05, label='DoAr 24')
ax2.plot(dftau_lineshape_data[0], (dftau_lineshape_data[1]-1.43)*1.66, label='DF Tau')
ax2.plot(dltau_lineshape_data[0], (dltau_lineshape_data[1]-0.895)*4.75, label='DL Tau')
ax2.plot(gktau_lineshape_data[0], (gktau_lineshape_data[1]-0.643)*6.8, label='GK Tau')
ax2.plot(rno90_lineshape_data[0], (rno90_lineshape_data[1]-2.32)*0.95, label='RNO 90')

#Setting labels, limits, legend
ax2.set_xlabel('Velocity [km/s]')
ax2.set_ylabel('Arbitrary flux')
ax2.set_title('Normalized Composite Lineshape of Classical Disks')
ax2.legend(loc="upper center", bbox_to_anchor=(0.5, 1.2), ncol=2)
ax2.set_ylim(0, 1)
ax2.set_xlim(-150, 150)
#Plotting composite lineshapes 5
fig=plt.figure(figsize=(16,8))

#Plotting composite lineshapes
ax3=fig.add_subplot(121)
ax3.plot(abaur_lineshape_data[0], abaur_lineshape_data[1]-9.38, label='AB Aur')
ax3.plot(hd36112_lineshape_data[0], hd36112_lineshape_data[1]-2.85, label='HD 36112')
ax3.plot(hd144432_lineshape_data[0], hd144432_lineshape_data[1]-2.05, label='HD 144432')

#Setting labels, limits, legend
ax3.set_xlabel('Velocity [km/s]')
ax3.set_ylabel('Arbitrary flux')
ax3.set_title('Composite Lineshape of Herbig Ae/Be Disks')
ax3.legend(loc="upper center", bbox_to_anchor=(0.5, 1.2), ncol=1)
ax3.set_ylim(0, 2)
ax3.set_xlim(-150, 150)
#Plotting composite lineshapes 6
fig=plt.figure(figsize=(16,8))

#Plotting composite lineshapes
ax3=fig.add_subplot(121)
ax3.plot(abaur_lineshape_data[0], (abaur_lineshape_data[1]-9.36)*0.55, label='AB Aur')
ax3.plot(hd36112_lineshape_data[0], (hd36112_lineshape_data[1]-2.91)*1.85, label='HD 36112')
ax3.plot(hd144432_lineshape_data[0], (hd144432_lineshape_data[1]-2.11)*5.4, label='HD 144432')

#Setting labels, limits, legend
ax3.set_xlabel('Velocity [km/s]')
ax3.set_ylabel('Arbitrary flux')
ax3.set_title('Normalized Composite Lineshape of Herbig Ae/Be Disks')
ax3.legend(loc="upper center", bbox_to_anchor=(0.5, 1.2), ncol=1)
ax3.set_ylim(0, 1)
ax3.set_xlim(-150, 150)
#Plotting composite lineshapes 7
fig=plt.figure(figsize=(16,8))

ax1=fig.add_subplot(121)
ax1.plot(citau_lineshape_data[0], citau_lineshape_data[1]-1.55, label='Transitional Disks', color='red')
ax1.plot(lkha330_lineshape_data[0], lkha330_lineshape_data[1]-1.03, color='red')
ax1.plot(twhya_lineshape_data[0], twhya_lineshape_data[1]-0.18, color='red')
ax1.plot(doar44_lineshape_data[0], doar44_lineshape_data[1]-0.43, color='red')
ax1.plot(hd135344_lineshape_data[0], hd135344_lineshape_data[1]-3.18, color='red')
ax1.plot(uxtau_lineshape_data[0], uxtau_lineshape_data[1]-0.33, color='red')


#Plotting composite lineshapes
ax1.plot(doar24_lineshape_data[0], doar24_lineshape_data[1]-2.02, label='Classical Disks', color='green')
ax1.plot(dftau_lineshape_data[0], dftau_lineshape_data[1]-1.4, color='green')
ax1.plot(dltau_lineshape_data[0], dltau_lineshape_data[1]-0.85, color='green')
ax1.plot(gktau_lineshape_data[0], gktau_lineshape_data[1]-0.58, color='green')
ax1.plot(rno90_lineshape_data[0], rno90_lineshape_data[1]-2.32, color='green')


#Plotting composite lineshapes
ax1.plot(abaur_lineshape_data[0], abaur_lineshape_data[1]-9.38, label='Herbig Ae/Be Disks', color='blue')
ax1.plot(hd36112_lineshape_data[0], hd36112_lineshape_data[1]-2.85, color='blue')
ax1.plot(hd144432_lineshape_data[0], hd144432_lineshape_data[1]-2.05, color='blue')

#Setting labels, limits, legend
ax1.set_xlabel('Velocity [km/s]')
ax1.set_ylabel('Arbitrary flux')
ax1.set_title('Composite Lineshapes of Protoplanetary Disks')
ax1.legend(loc="upper center", bbox_to_anchor=(0.5, 1.2), ncol=1)
ax1.set_ylim(0, 2)
ax1.set_xlim(-150, 150)
#Plotting composite lineshapes 8
fig=plt.figure(figsize=(16,8))

ax1=fig.add_subplot(121)

#Plotting composite lineshapes - Transitional
ax1=fig.add_subplot(121)
ax1.plot(citau_lineshape_data[0], (citau_lineshape_data[1]-1.61)*2.17, label='Transitional Disks', color='red')
ax1.plot(lkha330_lineshape_data[0], (lkha330_lineshape_data[1]-1.1)*5.6, color='red')
ax1.plot(twhya_lineshape_data[0], (twhya_lineshape_data[1]-0.242)*6.5, color='red')
ax1.plot(doar44_lineshape_data[0], (doar44_lineshape_data[1]-0.505)*7.8, color='red')
ax1.plot(hd135344_lineshape_data[0], (hd135344_lineshape_data[1]-3.23)*3.15, color='red')
ax1.plot(uxtau_lineshape_data[0], (uxtau_lineshape_data[1]-0.416)*54, color='red')

#Plotting composite lineshapes - Classical
ax1.plot(doar24_lineshape_data[0], (doar24_lineshape_data[1]-2.02)*1.05, label='Classical Disks', color='green')
ax1.plot(dftau_lineshape_data[0], (dftau_lineshape_data[1]-1.43)*1.66, color='green')
ax1.plot(dltau_lineshape_data[0], (dltau_lineshape_data[1]-0.895)*4.75, color='green')
ax1.plot(gktau_lineshape_data[0], (gktau_lineshape_data[1]-0.643)*6.8, color='green')
ax1.plot(rno90_lineshape_data[0], (rno90_lineshape_data[1]-2.32)*0.95, color='green')


#Plotting composite lineshapes - Herbig
ax1.plot(abaur_lineshape_data[0], (abaur_lineshape_data[1]-9.36)*0.55, label='Herbig Ae/Be Disks', color='blue')
ax1.plot(hd36112_lineshape_data[0], (hd36112_lineshape_data[1]-2.91)*1.85, color='blue')
ax1.plot(hd144432_lineshape_data[0], (hd144432_lineshape_data[1]-2.11)*5.4, color='blue')

#Setting labels, limits, legend
ax1.set_xlabel('Velocity [km/s]')
ax1.set_ylabel('Arbitrary flux')
ax1.set_title('Normalized Composite Lineshapes of Protoplanetary Disks')
ax1.legend(loc="upper center", bbox_to_anchor=(0.5, 1.2), ncol=1)
ax1.set_ylim(0, 1)
ax1.set_xlim(-150, 150)
#Plotting composite lineshapes 9
fig=plt.figure(figsize=(16,8))

ax1=fig.add_subplot(121)

#Plotting composite lineshapes - Transitional
ax1=fig.add_subplot(121)
ax1.plot(citau_lineshape_data[0], (citau_lineshape_data[1]-1.61)*2.17, label='Transitional Disks', color='red')
ax1.plot(lkha330_lineshape_data[0], (lkha330_lineshape_data[1]-1.1)*5.6, color='red')
ax1.plot(twhya_lineshape_data[0], (twhya_lineshape_data[1]-0.242)*6.5, color='red')
ax1.plot(doar44_lineshape_data[0], (doar44_lineshape_data[1]-0.505)*7.8, color='red')
ax1.plot(hd135344_lineshape_data[0], (hd135344_lineshape_data[1]-3.23)*3.15, color='red')
ax1.plot(uxtau_lineshape_data[0], (uxtau_lineshape_data[1]-0.416)*54, color='red')

#Plotting composite lineshapes - Herbig
ax1.plot(abaur_lineshape_data[0], (abaur_lineshape_data[1]-9.36)*0.55, label='Herbig Ae/Be Disks', color='blue')
ax1.plot(hd36112_lineshape_data[0], (hd36112_lineshape_data[1]-2.91)*1.85, color='blue')
ax1.plot(hd144432_lineshape_data[0], (hd144432_lineshape_data[1]-2.11)*5.4, color='blue')

#Setting labels, limits, legend
ax1.set_xlabel('Velocity [km/s]')
ax1.set_ylabel('Arbitrary flux')
ax1.set_title('Normalized Composite Lineshapes of Transitional and Herbig Ae/Be Disks')
ax1.legend(loc="upper center", bbox_to_anchor=(0.5, 1.2), ncol=1)
ax1.set_ylim(0, 1)
ax1.set_xlim(-150, 150)<jupyter_output><ipython-input-179-39c61d727549>:7: MatplotlibDeprecationWarning: Adding an axes using the same arguments as a previous axes currently reuses the earlier instance.  In a future version, a new instance will always be created and returned.  Meanwhile, this warning can be suppressed, and the future behavior ensured, by passing a unique label to each axes instance.
  ax1=fig.add_subplot(121)<jupyter_text>Compare Star Properties and Parameters<jupyter_code>data = pd.read_csv("/Users/erichegonzales/Desktop/eriche-thesis/data/star_data.csv")
print(data)
from matplotlib import image

img = image.imread("/Users/erichegonzales/Desktop/ASTR300/doar_p35.png")
plt.imshow(img)

plt.show()<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.1532064928153273, 'fraction_numerical': 0.09073975518893028, 'mean_word_length': 9.594080338266386, 'pattern_counts': {'":': 0, '<': 9, '<?xml version=': 0, '>': 9, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '2', 'dir': 'jupyter-structured-clean-dedup', 'id': '479349', 'license': 'MIT', 'n_tokens_mistral': 7396, 'n_tokens_neox': 6574, 'n_words': 813, 'path': 'data/compare_data/compare_linewidths.ipynb', 'repo_name': 'erichegonzales/eriche-thesis', 'size': '15032'}