HepSim
Repository with Monte Carlo simulations for particle physics
- Nov 1, 2022, Higgs, WW, ZZ production at 240, 250 GeV (e+e-,mu+mu-)
- Oct 27, 2022, Higgs, WW, ZZ production at 380 GeV
- Oct 17, 2022: hs-tools was corrected for https
- Sep 27, 2022: Interface improvements
- Feb 01, 2022: Migration to ALCF's Eagle filesystem
- Feb 14, 2021: ROOT6 support for Jas4pp program
- Nov 25, 2020: Moving to a new web server.
- Oct 22, 2020: Added OSG Connect storage (Snowmass21)
- Aug 26, 2020: Radion model Wkk+gg.
- June 5, 2020: Moving to a new (larger) web storage
- May 21, 2020: Filtered 2-leptons in multijet QCD
- Apr, 24, 2020: MG5/Pythia8 samples for KK+radion model
- Feb 20, 2020: Several samples for dark QCD processes
- Nov 28, 2019: DOI (Digital Object Identifier) by OSTI DOE (see osti.gov)
- Apr 15, 2019: Moving to Globus (Petrel)
- Sep.10 2018: Zprime/DM event samples
- Mar.15 2018: Charged Higgs event samples
- Sep,22 2017: Z+Higgs → nunu+XX event samples
- Sep,15 2017: Higgs → mu+mu- event samples
- Sep,10 2017: rfull059 tag with improved tracking strategy
- Aug.25, 2017: rfull300 for the DPF (ATLAS/CMS -like) detector
- Aug.7, 2017: J/Psi and Upsilon(S1) event samples for ep 45 GeV
- Jun.29, 2017: rfull058 tag with improved tracking strategy from D.Blyth
- Jun.20, 2017: rfull057 tag with alternative tracking strategy from D.Blyth
- Jun.9, 2017: Reprocessing rfull009 - rfull015 tags after correcting timing problem in SLIC. Using modifications for low-memory footprint
- Jun.9, 2017: Reprocessing rfull009 - rfull015 tags after correcting timing problem in SLIC. Using modifications for low-memory footprint
- May.16, 2017: Production of rfull056 using SiEIC(v5) detector for EIC
- Apr.30, 2017: CLIC 3 TeV e+e- samples using Pythia8
- Apr.20, 2017: Started production of ep 35 GeV samples
- Apr.13, 2017: New rfull015 with Geant10.3p1 using SiFCC(v7)
- Apr.3, 2017: rfull054 and rfull055 using SiEIC(v4) detector
- Feb.3, 2017: CLIC event samples for e+e- at 380 GeV. Link to events
- Feb.1, 2017: Updated rfull053 using SiEIC(v3) detector including track timing
- Nov.12, 2016: Production of rfull053 using SiEIC(v3) detector for EIC
- Nov.2, 2016: Production of rfull101 using SiCEPC(v2) detector for CEPC
- Nov.2, 2016: Production of Higgs+V for different CM energies (8-100 TeV pp)
- Oct.31,2016: Production of rfull052 using SiEIC(v2) detector for EIC
- Oct.14, 2016: Production of rfast005 for FCC-hh (pp 100 TeV) using Delphes-3.3.3
- Sep.23, 2016: Production of rfull051 using SiEIC(v1) detector for EIC
- Sep.15, 2016: Z'(5 TeV) to different channels using several SiFCC(v7) geometries
- Aug.28, 2016: rfull010, rfull011, rfull012and rfull013 for SiFCC(v7) using HCAL cells from 1 to 20 cm
- Aug.11, 2016: Production of rfull009 for 100 TeV (pp) with SiFCC-hh (v7) detector using new Pandora
- Aug.9, 2016: OSG grid pack with new (fast) PandoraPFA from J.Marshall
- Jul.27, 2016: Simulation of SiFCC-hh (v7) detector for 100 TeV (pp) (rfull008)
- Jul.24, 2016: Inclusive jets (100 TeV pp) tev100_qcd_pythia8_ptall
- Jul.13, 2016: Increase in statistics for ttbar+b (13 TeV pp) to 2.1 ab-1 tev13_mg5_ttbar_bjet
- Jun.20, 2016: Samples with single and double K-long for calorimeter studies. See KL samples
- May 19, 2016: Creating rfull007 for the SiFCC-hh (v5) detector with coarse HCAL granularity
- May 19, 2016: Re-processing rfull006 for SiFCC-hh (v4) after fixing endcap.
- Apr 8, 2016: H+ttbar (MG5) for 13 TeV (pp) (link)
- Apr 3, 2016: A new tag for fast simulation of 14 TeV (pp) (rfast004)
- Mar.29, 2016: Simulation of SiFCC-hh (v4) detector for 100 TeV (pp) (rfull006)
- Mar.26, 2016. All data sources were redirected to OSG due to a problem at ANL
- Mar.9, 2016: Fast simulation of ttbar+N jet process (pp, 14 TeV, MG5) (link)
- Mar.4, 2016: Full simulation of SiFCC-hh (v3) detector for 100 TeV (pp) (rfull005)
- Feb.5, 2016: Single particles for ITK studies (ATLAS phase II upgrade) (link)
- Feb.1, 2016: Z' with M=10,20,40 TeV decaying to qqbar, ttbar, WW for full simulations (link)
- Jan.19, 2016: 10 TeV Z' using a full simulation with 40 and 64 HCAL layers (link)
- Jan.14, 2016: TTbar+N jet process (pp, 14 TeV, MG5) (link)
- Jan.06, 2016: Heavy Higgs simulation (mu+mu-, 5 TeV) (link)
- Dec.17, 2015: Full SiD detector simulation of Zprime (10 TeV) to WW (link)
- Dec.17, 2015: Heavy higgs simulation for pp at 100 TeV (link)
- Dec.07, 2015: Full SiD detector simulation of Zprime to tautau (link)
- Nov.25, 2015: Particle gun samples for detector performance studies (pgun)
- Nov.18, 2015: Simulation of ttbar+bjet at 13,14,100 TeV (mg5_ttbar_bjet)
- Nov.9, 2015: Full simulation for e+e- (250 GeV) for SiD-CC (rfull002)
- Nov.6, 2015: Fast simulation of DIS events for EIC (141gev%rfast001)
- Oct.22, 2015: DIS events at the EIC collider (141 GeV)
- Oct.16, 2015: Delphes 3.3 fast simulation for ATLAS-like (13tev%rfast002) and CMS-like (13tev%rfast003) detectors. Same for 14 TeV.
- Oct.16, 2015: b-tagging was corrected for the tag rfast002
- Oct.15, 2015: Please update hs-toolkit.tgz
- Oct.9, 2015: Delphes 3.3 simulation of pp events (100 TeV) using the FCC detector (rfast002)
- Oct.6, 2015: Full simulation. e+e- events (250 GeV) for the SiD detector (rfull001)
- Sep.27, 2015: Fast simulation. e+e- events (250 GeV) for the ILD detector (rfast001)
truth_q2_dislowq2_lcio.py (raw text file)# Validation script: Read LCIO files
# Use Jas4pp as: fpad script.py
# Task: Extract truth-level Q2 variable and compare with Q2 from PFA
# S.Chekanov (ANL)
from org.lcsim.lcio import LCIOReader
from hep.io.sio import SIOReader
from hep.lcio.implementation.sio import SIOLCReader
from hep.lcio.implementation.io import LCFactory
from hep.lcio.event import *
from hep.lcio.io import *
from java.awt import Color,Font
from jhplot import * # import graphics
from jhplot.io import HBook
from hephysics.particle import LParticle
from java.lang import System
from java.io import *
from java.net import URL
from hepsim import HepSim
import sys,math
from java.awt import Color
from org.apache.commons.io import FileUtils
NFilesMax=10 # max number of files
dataset="gev35ep_lepto6ard_dislowq2"
tag="rfull056"
www=""
if len(sys.argv)>1:
flist=FileList.get(sys.argv[1],"slcio")
if (len(sys.argv)>2): FilesMax=int(sys.argv[2])
else:
sites=HepSim.urlRedirector(dataset)
www=sites[0]+"/"+tag
flist=HepSim.getList(www)
if len(flist)==0: print "Error: No input file!"; sys.exit(0)
else: print "Reading "+str(len(flist))+" files. Nr of files= ",NFilesMax
h1=H1D("Q2(ele true)",100,0,100)
h2=H1D("Q2(ele reco)",100,0,100)
factory = LCFactory.getInstance()
nEvent=0
for n in range(len(flist)):
if (n+1>NFilesMax): break
print "Open file=", flist[n]
url=URL(www+"/"+flist[n])
xfile=File("/tmp/tmp.slcio");
FileUtils.copyURLToFile(url, xfile)
reader = factory.createLCReader()
reader.open(xfile.getAbsolutePath())
while(1):
evt=reader.readNextEvent()
if (evt == None): break
nEvent=nEvent+1
if (nEvent%100==0): print "# Event: ",nEvent
params=evt.getParameters()
Q2true=params.getFloatVal("EVGEN:DIS:Q2")
h1.fill(Q2true);
# get electron beam energy
col = evt.getCollection("MCParticle")
nMc=col.getNumberOfElements()
e_beam=0;
for i in range(nMc): # loop over all particles
par=col.getElementAt(i)
momentum = par.getMomentum()
ee=par.getEnergy()
if (abs(par.getPDG())==11):
e_beam=ee;
break
# get reconstructed electrons using PFA
col = evt.getCollection("PandoraPFOCollection")
nPFA = col.getNumberOfElements()
particles=[]
for i in range(nPFA):
# http://www.lcsim.org/sites/lcsim/apidocs/org/lcsim/lcio/SIOReconstructedParticle.html
# http://www.lcsim.org/sites/lcsim/apidocs/org/lcsim/event/ReconstructedParticle.html
pa=col.getElementAt(i)
charge=pa.getCharge()
p4=pa.getMomentum()
ee=pa.getEnergy()
p=LParticle(p4[0],p4[1],p4[2],ee)
p.setCharge(charge)
if (abs(par.getPDG())==11):
particles.append(p)
break
if (len(particles)>0):
ele=particles[0]
mag=ele.mag()
theta=math.atan(ele.perp()/p.pz())
if(theta<0): theta+=3.14159
Q2ele=2.*e_beam*ele.e()*(1+math.cos(theta))
# print Q2true,Q2ele,e_beam
h2.fill(Q2ele);
del col
del evt
reader.close() # close the file
del reader
xfile.delete()
System.gc() # if you low on memory, manual JVM cleanup for every 5th file
c1=HPlot() # plotting part
c1.visible()
c1.setAutoRange()
c1.setMarginLeft(90)
c1.setNameX("Q^{2} [GeV^{2}]")
c1.setNameY("Events")
c1.setRangeX(0,10)
c1.draw(h1)
p2=P1D(h2)
p2.setColor(Color.blue)
c1.draw(p2) # show reco Q2
l2=HLabel("=HepSim=",0.9,0.86, "NDC")
l2.setColor(Color.gray)
l2.setFont(Font("Helvetica", Font.PLAIN, 14))
c1.add(l2)
name="output"
if len(sys.argv[0])>0: name=sys.argv[0].replace(".py","")
file=HBook(name+".jdat","w"); print name+".jdat created"
file.write(h1)
file.close()
c1.export(name+".svg"); print name+".svg created"
# xsec.toFile(name+".txt"); print name+".txt created"
# sys.exit(0)