#### kaggle-dark-worlds /HaloSize/halosize.py

 #carbonads{margin: 0 auto;padding: 1em;max-width: 26em !important;border: solid 2px #428bca;border-radius: 3px;font-size: .9em;font-family: Verdana, sans-serif;line-height: 1.5;} #carbonads * {display: block;overflow: hidden;} .carbon-img{float: left; margin-right: 1em !important; margin-bottom: .5em;} .carbon-text{float: left;max-width: 12em; text-align: left;} .carbon-poweredby{margin-top: -1em;text-align: right;font-size: .9em !important;} Language Python Lines 34 MD5 Hash bc8ae6e9b678956c05fd4163de1d60d4 Estimated Cost \$562 (why?) Repository https://bitbucket.org/paheliwala/kaggle-dark-worlds.git View Raw File View Project SPDX Find Similar Files View File Tree
 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31``` ```import pandas import cmath import math import numpy halos = pandas.read_csv( "../Raw/Training_halos.csv" ) def distance( ): pass def gal_heading( v ): a = math.atan2( v['e2'],v['e1'] )/2 + math.pi/2; return a def estimate_size( sky_name ): filename = "../Raw/Train_Skies/Training_%s.csv" % sky_name h1 = halos[halos.SkyId == sky_name ] halo_z = float( h1.halo_x1 ) + 1j*float( h1.halo_y1 ) sky = pandas.read_csv( filename ) sky['z1'] = sky.x+1j*sky.y - halo_z sky['l1'] = sky['z1'].map(lambda x: abs(x)) sky['a1'] = sky['z1'].map(lambda x: cmath.phase(x)) # Force angle sky_sub = sky[sky.l1 < 700] sky_sub['ga'] = sky_sub.apply(gal_heading,1) sky_sub['sa'] = numpy.abs( numpy.sin( sky_sub['ga'] - sky_sub['a1'] ) ) return ( sky_sub['sa'].mean() ) #print(halos) for i in xrange(1,301) : sky_name = "Sky%d" % i print i, estimate_size( sky_name ) ```