Project

General

Profile

Simulation » History » Version 3

Alexandre Camsonne, 11/01/2021 02:02 AM

1 1 Alexandre Camsonne
h1. Simulation
2
3
root -l s3https://dtn01.sdcc.bnl.gov:9000/eictest/ATHENA/FULL/master/SR/SR.10GeV_5kVthreshold_hepmc/25098.0001.root -e 'events->Draw("@VertexBarrelHits.size()")'
4
5
https://dtn01.sdcc.bnl.gov:9000/minio/eictest/ATHENA/FULL/master/SR/SR.10GeV_5kVthreshold_hepmc/
6
7
eicS3read for both u/p
8
9
https://indico.bnl.gov/event/12817/
10
11
12
h1. Running 
13
14
Wouter Deconinck (he/him):juggling: Sep 2nd at 10:18 AM
15
Here are instructions (thread)
16
3 replies
17
18
Wouter Deconinck (he/him):juggling:  2 months ago
19
Mirror SR files over to a local directory (254 MB):
20
mc mirror S3/eictest/ATHENA/EVGEN/SR/ SR/
21
22
Wouter Deconinck (he/him):juggling:  2 months ago
23
Install eic-shell to current directory:
24
curl -L get.athena-eic.org | bash
25
26
Wouter Deconinck (he/him):juggling:  2 months ago
27
Run the simulation on one of the files inside container:
28
cat << EOF | ./eic-shell
29
source /opt/detector/setup.sh
30
npsim \
31
  --runType batch \
32
  --printLevel WARNING \
33
  --skipNEvents 0 \
34
  --numberOfEvents 1000 \
35
  --part.minimalKineticEnergy 1*TeV \
36
  --hepmc3.useHepMC3 true \
37
  --compactFile /opt/detector/share/athena/athena.xml \
38
  --inputFiles SR/SR.10GeV_5kVthreshold_hepmc/25098.hepmc \
39
  --outputFile 25098.root
40 2 Alexandre Camsonne
41
42
h1. https://indico.bnl.gov/event/12809/
43
44
Input SR files are at https://dtn01.sdcc.bnl.gov:9000/minio/eictest/ATHENA/EVGEN/SR/SR.10GeV_5kVthreshold_hepmc/, which can be synced with mc mirror S3/eictest/ATHENA/EVGEN/SR/SR.10GeV_5kVthreshold_hepmc/ SR/ (it's not that large). More on setting up mc to access S3 at https://doc.athena-eic.org/en/latest/howto/s3_file_storage.html. mc is installed inside the container below, but you will still need to setup the config with username eicS3read and password eicS3read.
45
After syncing the hepmc files to the SR/ directory, you'll need to install the ATHENA eic-shell container without cvmfs (you'll need to modify it, which isn't easily done using cvmfs containers). So, instead of the simplest version of the following, use curl -L get.athena-eic.org | bash -s -- --no-cvmfs (with additional no-cvmfs option). This will download a container image to your current directory (~2.5 GB). Next you should edit the script eic-shell and add two flags to the last line: singularity exec --writable --fakeroot (unfortunately a rather slow step).
46
When you start the container with ./eic-shell, first source /opt/detector/setup.sh and then run the full simulation with the files in SR/ for example as
47
npsim --runType batch --inputFiles SR/25098.hepmc --outputFile SR/25098.root --compactFile /opt/detector/share/athena/athena.xml --numberOfEvents 10
48
You probably want to use --filter.tracker edep0 since default is edep1kev though I have not seen it have any impact (might have been due to the gold layer).
49
All this up till now is without making any changes. To make changes, in particular gold layer, change gold_thickness in /opt/detector/share/athena/ip6/beampipe.xml (that's the only reason why this is more complicated than usual...).
50
For analysis (also inside container), you can start with something like
51
root -l SR/25098.root -e 'events->Scan("VertexBarrelHits.energyDeposit")'
52
But with edep0 filters, and with 0 um gold layer, there is still nothing in at least the first 10 events I tested.
53
Next steps:
54
Maybe run npsim with the --part.enableDetailedHitsAndParticleInfo flag to see what is actually happening (for few events),
55
Create verbose.mac with
56
/tracking/verbose 2
57
/run/beamOn 10
58
and run this with --runType run --macroFile verbose.mac instead. That should definitely show where the SR photons are going. If you're familiar with geant4 detailed output...
59
I see some tracks like this:
60
*********************************************************************************************************
61
* G4Track Information:   Particle = gamma,   Track ID = 1,   Parent ID = 0
62
*********************************************************************************************************
63
64
Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
65
    0     29.5      0.5 4.31e+03   0.00752        0        0         0 AV_934!BeamPipe_assembly_17#17!v_downstream_vacuum_12#12 initStep
66
    1       30    0.503 4.47e+03   0.00752        0      161       161 AV_934!BeamPipe_assembly_17#17!v_downstream_matter_11#11 Transportation
67
    2       30    0.503 4.47e+03         0  0.00155   0.0221       161 AV_934!BeamPipe_assembly_17#17!v_downstream_matter_11#11 phot
68
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  1 ---------------
69
    :        30     0.503  4.47e+03   0.00597                 e-
70
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
71
The beampipe volumes are v_downstream_vacuum for the vacuum, and v_downstream_matter for the aluminum matter. So at least this photon is being simulated and it hit the beam pipe at z = 4.47m.
72 3 Alexandre Camsonne
73
74
75
spack:
76
  specs:
77
  - acts +dd4hep +digitization +examples +fatras +identification +json +tgeo +ipo
78
  - cmake
79
  - clhep cxxstd=17
80
  - dd4hep +geant4 +assimp +hepmc3 +lcio +ipo
81
  - eigen
82
  - gaudi@36.0 +ipo
83
  - geant4 +ipo +qt +opengl -python +vecgeom cxxstd=17
84
  - genfit +ipo
85
  - hepmc3 +interfaces +python +rootio
86
  - lcio +ipo
87
  - mesa -llvm swr=none
88
  - opencascade +ipo +tbb -vtk
89
  - podio +ipo
90
  - pythia8 +fastjet
91
  - root cxxstd=17 +fftw +fortran +gdml +mlp +pythia8 +root7 +tmva +vc +xrootd +ssl
92
  - xrootd cxxstd=17 +python
93
  concretization: together
94
  config:
95
    install_missing_compilers: true
96
    install_tree:
97
      root: /opt/software
98
  view: /opt/local