直線で描く場合。 import matplotlib.pyplot as plt fig, ax = plt.subplots() ax.plot([0,1],[0,1]) plt.show() dashes 変数を使うと、点線の色あり、色なしの部分の幅を指定できる。以下はその例。 import matplotlib.pyplot as plt fig, ax = plt.subplot…

Google Apps Script(GAS)でTOTPを生成するには、SHA1を実装している外部ライブラリを使うのが良いだろう。 ウェブアプリとしてJSON形式で表示させるところまで実装したコード function myFunction() { const secret = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"; const c…

Google Authenticator（Google認証システム、認証アプリ）にインポートした 時間ベースのワンタイムパスワード TOTP のシークレットキーを確認したいことがあるのでその方法。 Google Authenticatorから、アカウントを移行→アカウントのエクスポートで、必要…

あらゆるものをWindows上で動かさないと気がすまない筆者なので、頑張ってExperimental Physics and Industrial Control System (EPICS) R3.15をWindows上で動かすことにする。7.0用の記事もあるが、Windowsで動かすには色々バグがあるかもしれない。 ビルド…

あらゆるものをWindows上で動かさないと気がすまない筆者なので、頑張ってExperimental Physics and Industrial Control System (EPICS) 7.0をWindows上で動かすことにする。 必要なアイテムは4つ Windows用のmake (PATHを通す必要がある) Windows用のperl V…

インストール済み UWP (ユニバーサル Windows プラットフォーム) アプリ表示するPowerShell コマンド Get-AppxPackage | Sort-Object Name | Select Name,PackageFullName Microsoft StoreからはアンインストールできないUWPアプリをアンインストールするコ…

const periodicTable = [ "H", "He", "Li", "Be", "B", "C", "N", "O", "F", "Ne", "Na", "Mg", "Al", "Si", "P", "S", "Cl", "Ar", "K", "Ca", "Sc", "Ti", "V", "Cr", "Mn", "Fe", "Co", "Ni", "Cu", "Zn", "Ga", "Ge", "As", "Se", "Br", "Kr", "Rb", "Sr…

Experimental Physics and Industrial Control System (EPICS) とは何か? Wikipediaには EPICS（Experimental Physics and Industrial Control System）は、加速器、望遠鏡、その他の大規模は実験用機器を運用する分散制御システムを開発・実装するのに使わ…

一時的にファイルの変更を無視するためのコマンドは git update-index --skip-worktree [ファイル名] 似たコマンドで --assume-unchanged がある。 qiita.com --skip-worktree は変更が 「ある」 ものとして無視するが、--assume-unchanged は変更が 「ない…

def get_nuclei_yield(filename): lines = open(filename).readlines() flag=False nuclei_yield={} for i, line in enumerate(lines): if line=="[Calculations]\n": flag=True continue if flag: str1=line.split()[0] str2=line.split("=")[1].split(",")…

下記のコードではJupyter Labを使っています。 matplotlibの他に、ipymplが必要なので、pip install ipymplでインストールしてください。 コードを完全に理解していないため、解説は他の記事を参照してください。 マウスでクリックした時の、回数、座標、そ…

root.cern.ch Better palettes like kBird are recommended. と書かれているので使ってみよう。 import matplotlib red_bird = [ 0.2082, 0.0592, 0.0780, 0.0232, 0.1802, 0.5301, 0.8186, 0.9956, 0.9764] green_bird = [ 0.1664, 0.3599, 0.5041, 0.6419,…

同じデータを使って、CERN ROOT、Pythonのscipy、lmfit を使って正規分布でフィッティングするコードと結果を比較する。 まずは、データ生成部分 C++用 std::vector<double> vx{-3.65,-3.45,-3.35,-3.25,-3.15,-3.05,-2.95,-2.85,-2.75,-2.65,-2.55,-2.45,-2.35,-2.2</double>…

グリッド間隔、線の種類、線の幅、線の色、用紙サイズは自由に設定できるようにした。 要点は次の3つ plt.subplots_adjust()で枠線の外側をなくしてfigureを紙面全体にする plt.gca().axis("off") で枠線、目盛り、目盛りラベルを消す plt.plot(transform=pl…

matplotlib公式のlegendのチュートリアルの一部をこの記事で紹介する。 matplotlib.org matplotlibにおける凡例はplt.legend() とする。 label で指定したグラフの凡例が、グラフの隙間に自動的に描画される。 import numpy as np import scipy import matpl…

普通に描画する fig, ax = plt.subplots() ax.plot([0,1],[0,1]) plt.show() 軸を全て消す。 fig, ax = plt.subplots() ax.plot([0,1],[0,1]) ax.axis("off") plt.show() subplotsを使わない場合はaxのかわりにplt.gca()とする。 plt.plot([0,1],[0,1]) plt.…

ポッドキャストも聞けるSpotify(スポティファイ) は、開発者向けのAPIを公開しています。 SpotifyのAPIでPodcastの情報を取得するためのサンプルコードを公開します。 Client IDとClient secret は Spotify for Developers で取得します。 pip install spoti…

ChatGPTに書かせた def split_numbers(total, ratios): # 各割合に対応する数値を計算 numbers = [total * ratio for ratio in ratios] # 小数点以下を四捨五入 numbers = [round(number) for number in numbers] # 残りの誤差を調整 total_numbers = sum(nu…

ヒストグラムのビン(bin)の値を直接操作した後で、ヒストグラムっぽいグラフを作りたいことがあるだろう。 以下では、plt.hist を使った場合と、np.histgramで計算させた、binの位置 X とCounts Y を使って plt.hist 以外で描画する方法を紹介する。 データ…

matplotlib.patheffectsを使いこなすと良い。 黒抜き文字の例 import matplotlib.patheffects as patheffects import matplotlib.pyplot as plt plt.plot([0,1],[0,1],lw=50) plt.text(0.5,0.5,"testテスト",ha="center",va="center",color="black",fontsize…

改行を含む日本語を書きたい場合は\begin{document}の前に \usepackage{CJKutf8} \long\def\Ja#1{\begin{CJK}{UTF8}{ipxm}#1\end{CJK}} または \usepackage{CJKutf8} \newcommand{\Ja}[1]{\begin{CJK}{UTF8}{ipxm}#1\end{CJK}} と書く。使い方は \Ja{ ここに…

Yahama ルーター 192.168.0.200 から Windowsサーバー 192.168.0.100 へsyslogを送ることにする。 Yamahaルーターの設定。syslog debugは膨大なのでoffが望ましいだろう。 syslog host 192.168.0.100 syslog facility local1 syslog notice on syslog info o…

gigazine.net セットアップ時に、Shift+F10 でコマンドプロンプトを立ち上げて cd oobe で C:\Windows\System32\oobe ディレクトリに移動し BypassNRO.cmd を実行すればOK。 自動的に再起動されるので、「ネットワークに接続しましょう」をスキップして、ロ…

Office 365 (Microsoft 365) や、 OutlookのWebメールで、常に特定のメールアドレスにBCCしたい (いわゆる Always BCC) 需要があった。 Chromeの拡張機能を使うと簡単に実装できそうなのでやってみた。Chromeの拡張機能の作り方は以下のサイトを参考にしてほ…

放射性崩壊などの時間発展をPythonで計算してみよう。 寿命 life=10 として、微分方程式は dx/dt=-x/life なので、 import numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt def system(x, t): life = 10 dxdt = -x/life re…

まず、Wi-Fi(無線LAN)で記録されたプロファイル(SSID)を全て表示する。必要なのはプロファイル名。 >netsh wlan show profiles インターフェイス Wi-Fi のプロファイル: グループ ポリシー プロファイル (読み取り専用) -----------------------------------…

CRediT (Contributor Roles Taxonomy) つまり、貢献者の役割の分類法は、Elsevierによると、 個々の著者の貢献を認識し、著者の論争を減らし、共同作業を促進する目的 the intention of recognizing individual author contributions, reducing authorship d…

多面体クラスである polyhedron を使えば良いが、一見どう使えばいいのか分からない。 例えば、y=0の平面上の三角形から高さ100の三角柱を描く場合は以下のPythonコードのようにする。 from solid import * from solid.utils import * x0 = 0 x1 = 20 x2 = -…

from PyPDF2 import PdfWriter, PdfReader def overlay(): file1 = r"file1.pdf" file2 = r"file2.pdf" file3 = r"overlay.pdf" back_ground = PdfReader(open(file1, "rb")) page1 = back_ground.pages[0] fore_ground = PdfReader(open(file2, "rb")) page…

ROOTのTH2DやTH2Fなどのクラスでは、2次元ヒストグラムを、legoっぽく図示できる。 ご存知の通り、Python + matplotlib の hist2dには、このような描画をするコマンドはない(はずである)。 というわけで、作ってみた。3次元バーを描画する関数であるbar3dを…