def __init__(self):
#Import Dazer_Files Class to import data from lines_logs
Txt_Files_Manager.__init__(self)
#Declare Hydrogen and Helium collisional lines for the analysis
#self.posHydrogen_Ions = ['Hdelta_6_2', 'Hgamma_5_2', 'Hbeta_4_2', 'Halpha_3_2']
self.posHydrogen_Lines = ['H1_4102A', 'H1_4340A', 'H1_4861A', 'H1_6563A']
self.Hydrogen_Wavelengths = [4101.742, 4340.471, 4862.683, 6562.819]
self.posHelium_Lines = ['He1_3889A', 'He1_4026A', 'He1_4387A', 'He1_4471A', 'He1_4686A', 'He1_4714A', 'He1_4922A', 'He1_5876A', 'He1_6678A', 'He1_7065A', 'He1_7281A', 'He1_10380A']
self.Helium_Wavelengths = [3889.0, 4026.0, 4387.0, 4471.0, 4686.0, 4714.0, 4922.0, 5876.0, 6678.0, 7065.0, 7281.0, 10830.0]
self.Cand_Hydrogen_Lines = []
self.Cand_Helium_Lines = []
self.nHydrogen = None
self.nHelium = None
#Define indexes and labels to speed up the code
self.H13889A_label = 'H1_3889A'
self.HBeta_label = 'H1_4861A'
self.He3889_label = 'He1_3889A'
self.He3889blended_label = 'He1_3889A_blended'
self.He3889_check = None
#Set up the right emissivities
pn.atomicData.setDataFile('he_i_rec_Pal12-Pal13.fits')
#Declare pyneb Hydrogen and Helium atoms to calculate emissivities
self.H1 = pn.RecAtom('H', 1)
self.He1 = pn.RecAtom('He', 1)
print 'I am using this atomic data'
print self.He1.printSources()
#Import collisional coefficients table
self.Coef_Kalpha_dict = self.Import_Coll_Coeff_Table(self.posHydrogen_Lines, None)
#Import Optical depth function
self.Coef_ftau_dict = self.Import_OpticalDepth_Coeff_Table(self.posHelium_Lines)
#Declare dictionaries to store the data
#WARNING: In this analysis flambda is mantained constant
self.Flux_dict = OrderedDict()
self.Error_dict = OrderedDict()
self.Wave_dict = OrderedDict()
self.PynebCode_dict = OrderedDict()
self.EqW_dict = OrderedDict()
self.EqWerror_dict = OrderedDict()
self.hlambda_dict = OrderedDict()
self.flambda_dict = self.get_flambda_dict(self.posHydrogen_Lines + self.posHelium_Lines, self.Hydrogen_Wavelengths + self.Helium_Wavelengths)
#Extra
self.EmptyRowFormat = 'nan'
def __init__(self):
self.Hydrogen_CollCoeff_TableAddress = '/home/vital/git/Dazer/Dazer/dazer/libraries/Astro_Libraries/Neutral_Hydrogen_Collisional_Correction_coef.txt'
self.Helium_CollCoeff_TableAddress = '/home/vital/git/Dazer/Dazer/dazer/libraries/Astro_Libraries/Neutral_Helium_Collisional_Correction_coef.txt'
self.Helium_OpticalDepth_TableAddress = '/home/vital/git/Dazer/Dazer/dazer/libraries/Astro_Libraries/Helium_OpticalDepthFunction_Coefficients.txt'
#Import Dazer_Files Class to import data from lines_logs
Txt_Files_Manager.__init__(self)
#Declare Hydrogen and Helium lines for the analysis
self.posHydrogen_Lines = ['H1_4102A', 'H1_4340A', 'H1_4861A', 'H1_6563A']
self.Hydrogen_Wavelengths = [4101.742, 4340.471, 4862.683, 6562.819]
self.posHelium_Lines = ['He1_3889A', 'He1_4026A', 'He1_4387A', 'He1_4471A', 'He1_4686A', 'He1_4714A', 'He1_4922A', 'He1_5876A', 'He1_6678A', 'He1_7065A', 'He1_7281A', 'He1_10830A']
self.Helium_Wavelengths = [3889.0, 4026.0, 4387.0, 4471.0, 4686.0, 4714.0, 4922.0, 5876.0, 6678.0, 7065.0, 7281.0, 10830.0]
self.Cand_Hydrogen_Lines = []
self.Cand_Helium_Lines = []
self.nHydrogen = None
self.nHelium = None
#Define indexes and labels to speed up the code
self.H13889A_label = 'H1_3889A'
self.HBeta_label = 'H1_4861A'
self.He3889_label = 'He1_3889A'
self.He3889_Check = None
#Set up the right emissivities
pn.atomicData.setDataFile('he_i_rec_Pal12-Pal13.fits')
#Declare pyneb Hydrogen and Helium atoms to calculate emissivities
self.H1 = pn.RecAtom('H', 1)
self.He1 = pn.RecAtom('He', 1)
print 'Helium emissivities:\t', self.He1.printSources()
#Import collisional coefficients table
self.Coef_Kalpha_dict = self.Import_Coll_Coeff_Table(self.posHydrogen_Lines, None)
#Import Optical depth function
self.Coef_ftau_dict = self.Import_OpticalDepth_Coeff_Table(self.posHelium_Lines)
#Declare dictionaries to store the data
self.Flux_dict = OrderedDict()
self.Error_dict = OrderedDict()
self.Wave_dict = OrderedDict()
self.PynebCode_dict = OrderedDict()
self.EqW_dict = OrderedDict()
self.EqWerror_dict = OrderedDict()
self.hlambda_dict = OrderedDict()
self.flambda_dict = self.get_flambda_dict(self.posHydrogen_Lines + self.posHelium_Lines, self.Hydrogen_Wavelengths + self.Helium_Wavelengths) #flambda does not form part of the inference predictions
