Add files via upload

Author: Tszon

InitialMassFunction/HR_Diagram_AutomatedPlans.png

@@ -0,0 +1,68 @@
+# Cluster Match
+
+(Pre-requisites: Python >=3.7, Astropy, Numpy, Matplotlib)
+
+To install pre-requisities run: 
+`pip install astropy numpy matplotlib`
+
+For LRC computers you will need to specify `C:\Python39\python` instead of just `python` as the path to run python 3 (by default `python` uses python 2.7).
+
+To install astropy on LRC computers run `C:\Python39\Scripts\pip install astropy`
+
+If using Spyder you can add parameters by clicking Run>Configuration per file, then entering the parameters (not including the script name itself) into the "Command line options" box in the "General Settings" section (be careful NOT to use the "Command line options" box in the "Console" section higher up).
+
+Images will need to be plate-solved for these scripts (and APT) to run successfully.
+
+
+**Step 1:**
+
+Download the UCAC4 stars covering the area around an input file:
+If your files have spaces, don't forget to enclose the entire path in double quotes.
+
+Use:
+
+  `download.py <input file>`
+  
+e.g.
+
+  `>> python download.py "NGC 457_B_15s_B2_T3_72939.fit"`
+
+**Step 2:**
+
+For each filter (V and B), match stars in the image to the UCAC4 stars, and estimate the zero point:
+
+Use:
+
+  `match-zp.py <catalog file> <APT source list> <output file name> <filter (V or B)>`
+  
+e.g.
+
+ ` >> python match-zp.py "ngc 457_V_15s_B2_T3_72941-UCAC4.tbl" "NGC457-V-inst-mags.csv" "NGC457-V-cal-mags.csv" V`
+  
+  `>> python match-zp.py "ngc 457_V_15s_B2_T3_72941-UCAC4.tbl" "NGC457-B-inst-mags.csv" "NGC457-B-cal-mags.csv" B`
+  
+**Step 3:**
+
+Match stars from the V and B files, and calculate B-V colour:
+  
+Use:
+
+  `match-bv.py <B mag file> <V mag file> <output file>`
+  
+e.g.
+
+  `>> python match-bv.py "NGC457-B-cal-mags.csv" "NGC457-V-cal-mags.csv" "NGC457-BV.csv"`
+  
+**Step 4:**
+
+The CSV file can be opened in Excel or other software to plot the results.
+
+You can make a basic plot of B-V colours against apparent V magnitudes to create an HR diagram:
+  
+Use:
+
+  `plot.py <BV file>`
+    
+e.g.
+
+   `>> python plot.py "NGC457-BV.csv"`

InitialMassFunction/IMFplot_histogram.png

@@ -0,0 +1,55 @@
+from urllib import request, parse
+from astropy.io import fits
+import os
+import sys
+
+#read the filename from the arguments
+filename = sys.argv[1:][0]
+
+#Open the fits file header data unit
+hdul = fits.open(filename)
+
+try:
+    #Extract Right ascension of centre of the image
+    RA= float(hdul[0].header['CRVAL1'])	
+    print("RA: "+str(RA))
+
+    #Extract declination of centre of the image
+    DEC= float(hdul[0].header['CRVAL2'])	
+    print("DEC: "+str(DEC))
+    
+except:
+    print("Error: Could not find CRVAL1 or CRVAL2 in FITS header. This likely means the image has not been plate-solved.")
+    print("See https://observatory.herts.ac.uk/wiki/Plate_Solving for help with plate-solving images.")
+    quit()
+
+
+postData = {  '-source' : 'I/322A/out', '-out.max' : '10000',  '-out.form' : '| -Separated-Values',  '-oc.form' : 'dec',  '-c.eq' : 'J2000', '-c.r' : '++2', '-c.u' : 'arcmin', '-c.geom' : 'r', '-out.src' : 'I/322A/out', '-out.orig' : 'standard',  '-out' : 'RAJ2000, DEJ2000, Vmag, Bmag', 'Bmag' : '>0',  'Vmag' : '>0', 'RAJ2000' : str(RA)+'+/-0.5',  'DEJ2000' : str(DEC)+'+/-0.5'}
+
+#URL encode data
+
+data = parse.urlencode(postData).encode()
+
+
+#Create web request
+req =  request.Request('http://vizier.u-strasbg.fr/viz-bin/asu-tsv', data=data)
+
+#print("Downloading UCAC4 data")
+
+#Request web query
+response = request.urlopen(req)
+
+#Get return data
+return_data = response.read()
+
+#print(response)
+
+#print(return_data)
+
+print("Writing to file: "+os.path.splitext(filename)[0]+"-UCAC4.tbl")
+
+#write the return data to file
+with open(os.path.splitext(filename)[0]+"-UCAC4.tbl", 'wb') as s:
+    s.write(return_data)
+	
+print("Done")

InitialMassFunction/PolynomialFit.png

@@ -0,0 +1,120 @@
+import os
+import math
+import sys
+
+BMags= sys.argv[1:][0] # "NGC457-B-cal-mags.csv"
+
+VMags= sys.argv[1:][1] #"NGC457-V-cal-mags.csv"
+
+outFile= sys.argv[1:][2] #"NGC457-BV.csv"
+
+print("B magnitudes: "+BMags)
+
+print("V magnitudes: "+VMags)
+
+print("Output file: "+outFile)
+
+if os.path.exists(outFile):
+	print("Output file already exists")
+	sys.exit()	
+
+#Maximum distance from expected coordinates for a source to be considered (arcseconds)
+maxDist = 2
+
+#Open the corresponding source list for reading
+VMagsFile = open(VMags, "r")
+#skip headers
+Vfirst_line = VMagsFile.readline()
+VLines = VMagsFile.readlines()
+
+BMagsFile = open(BMags, "r")
+Bfirst_line = BMagsFile.readline()
+BLines = BMagsFile.readlines()
+
+
+
+
+lineNum=0
+			
+outlist = []
+
+x=0
+
+matches = 0
+
+print("---------------------")
+print("Matching stars...")
+
+#Loop through all lines in the file
+for Vline in VLines:
+
+
+	partsV = Vline.split(",")
+	
+	#Extract useful values from the source list
+	raV = float(partsV[0]) 	#Right ascensions (degrees)
+	decV = float(partsV[1]) 	#Declination (degrees)
+	Vmag = float(partsV[2].strip())			#Source instrument magnitude
+
+	firstLine = 0
+	
+	closest = 360
+	
+	bestB =0
+
+	for Bline in BLines:
+	
+	
+		partsB = Bline.split(",")
+
+		raB = float(partsB[0]) 	#Right ascensions (degrees)
+		decB = float(partsB[1]) 	#Declination (degrees)
+		Bmag= float(partsB[2].strip())			#Source instrument magnitude
+		
+		ra1 = math.radians(raB)
+		ra2 = math.radians(raV)
+		d1 = math.radians(decB)
+		d2 = math.radians(decV)			
+		
+		#Calculate the distance from this source to the target
+				
+		if abs(ra1-ra2)==0 and abs(d1-d2)==0:
+			angSep = 0
+		else:		
+			angSep = 3600*math.degrees(math.acos(math.sin(d1)*math.sin(d2)+math.cos(d1)*math.cos(d2)*math.cos(ra1-ra2)))
+		
+		if angSep < maxDist and angSep<closest:
+				
+			closest = angSep
+			
+			bestB = Bmag
+			
+	Bmag = bestB 
+
+	if closest<maxDist:
+	
+		string = str(raV)+","+str(decV)+","+str(round(Vmag,4))+","+str(round(Bmag,4))+","+str(round(Bmag-Vmag,4))
+		outlist.append(string)
+		matches = matches + 1
+
+			
+			
+print(str(matches) + " matches found")
+
+print("---------------------")
+print("Writing to file: " +outFile)	
+
+f= open(outFile,"a")	
+
+#write header
+f.write("RA,Dec,Vmag,Bmag,B-V\n")
+
+lineNum=0	
+	
+for line in outlist:
+
+	f.write(line+"\n")
+
+f.close()
+		
+print("Done")