{"cells":[{"metadata":{"_uuid":"8f2839f25d086af736a60e9eeb907d3b93b6e0e5","_cell_guid":"b1076dfc-b9ad-4769-8c92-a6c4dae69d19","trusted":true,"collapsed":true},"cell_type":"code","source":"import numpy as np\nimport pandas as pd\nimport matplotlib.pyplot as plt\nfrom mpl_toolkits.mplot3d import Axes3D\nfrom trackml.dataset import load_event\nimport os\n\npath_to_train = \"../input/train_1\"\nevent_prefix = \"event000001000\"\n\nhits, cells, particles, truth = load_event(os.path.join(path_to_train, event_prefix))\n\ndef plot_track(axis,particle_id):\n    p_traj = (truth[truth.particle_id == particle_id][['tx', 'ty', 'tz']]).sort_values(by='tz')\n    axis.plot(\n        xs=p_traj.tx,\n        ys=p_traj.ty,\n        zs=p_traj.tz)\n\ndef show_3Dplot(ax):\n    ax.set_xlabel('X (mm)')\n    ax.set_ylabel('Y (mm)')\n    ax.set_zlabel('Z  (mm) -- Detection layers')\n    plt.show()\n\ndef plot_sample_track(N):\n    fig, ax = plt.subplots()\n    ax = Axes3D(fig)\n    sample_particles = particles.sample(N)\n    if N > 1:\n        for particle_id in sample_particles.particle_id.values:\n            plot_track(ax,particle_id)\n    else:\n        plot_track(ax,sample_particles.particle_id.values[0])\n    show_3Dplot(ax)\n\nplot_sample_track(10)","execution_count":null,"outputs":[]},{"metadata":{"_cell_guid":"79c7e3d0-c299-4dcb-8224-4455121ee9b0","collapsed":true,"_uuid":"d629ff2d2480ee46fbb7e2d37f6b5fab8052498a","trusted":false},"cell_type":"code","source":"","execution_count":null,"outputs":[]}],"metadata":{"kernelspec":{"display_name":"Python 3","language":"python","name":"python3"},"language_info":{"name":"python","version":"3.6.5","mimetype":"text/x-python","codemirror_mode":{"name":"ipython","version":3},"pygments_lexer":"ipython3","nbconvert_exporter":"python","file_extension":".py"}},"nbformat":4,"nbformat_minor":1}