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