Abstract. An earthquake is an event that releases energy in the earth's crustal rocks, which can trigger aftershocks. The economic losses caused by earthquakes can be mitigated through earthquake insurance. However, conventional premium calculations only take into account the main earthquake (mainshock), while aftershocks are excluded through a declustering process. This article presents the development of Clustering Probabilistic Seismic Hazard Analysis (CPSHA) to evaluate earthquake risk by incorporating aftershock effects using a temporal Epidemic-Type Aftershock Sequence (ETAS) model. The approach is an improvement on traditional earthquake insurance premium calculations, which generally only consider the main earthquake. The ETAS model, with parameters θ=(μ,K,c,α,p), is estimated through Maximum Likelihood Estimation (MLE). This analysis combines seismic hazard ETAS with structural damage data to determine the pure premium of earthquake insurance with aftershock. The study uses West Java earthquake data from 21 November 2022 to 18 February 2023, sourced from the BMKG catalogue. The results indicate that customers in the Cianjur area must pay a total insurance premium (TP) of Rp 28.263.945,00 per year for clusters. This approach enhances the comprehension of earthquake risk and can aid in the calculation of more precise insurance premiums. Abstrak. Gempa bumi merupakan peristiwa pelepasan sejumlah energi pada batuan kerak bumi yang dapat memicu gempa susulan (aftershock). Dampak dari gempa bumi dapat menyebabkan kerugian ekonomi. Salah satu cara untuk menanggulangi hal tersebut dengan asuransi gempa bumi. Perhitungan premi konvensional hanya memperhitungkan gempa utama (mainshock) sementara aftershock disisihkan melalui proses declustering. Artikel ini menyajikan pengembangan Clustering Probabilistic Seismic Hazard Analysis (CPSHA) untuk mengevaluasi risiko gempa bumi dengan memasukkan efek aftershock menggunakan model Epidemic-Type Aftershock Sequence (ETAS) secara temporal. Perhitungan premi asuransi gempa bumi umumnya hanya mempertimbangkan gempa utama, namun pendekatan ini menyertakan aftershock dalam estimasi premi. Model ETAS dengan parameter θ=(μ,K,c,α,p) diestimasi melalui Maximum Likelihood Estimation (MLE). Analisis ini menggabungkan seismic hazard ETAS dengan data kerusakan struktur untuk menentukan premi murni asuransi gempa bumi dengan aftershock. Data yang digunakan adalah data gempa bumi Jawa Barat 21 November 2022 sampai 18 Februari 2023 melalui katalog BMKG. Hasil yang diperoleh menunjukkan total premi (TP) asuransi yang harus dibayarkan oleh nasabah untuk cluster di wilayah Cianjur sebesar Rp 28.263.946,00 per tahun. Pendekatan ini memberikan kontribusi pada pemahaman risiko gempa bumi dan dapat digunakan dalam perhitungan premi asuransi yang lebih akurat.