TRestG4Event::GetFirstPositionInVolume not working correctly?

lobis · July 12, 2019, 12:40pm

REST version : v2.2.11

Hello,

I am using the method TRestG4Event::GetFirstPositionInVolume(volume_id) to draw in 3D the points where the particle enters my volume, in my case a muon veto consisting of a rectangle on top of the detector.

I have found that sometimes I get back a position way inside the volume, even though when viewing individual events the step size seems to be very small inside the veto.

Also on some instances I get a position way outside the veto, which is the reason I am making this post. This can be seen in the images below. I have gotten all my 3D points as position = event->GetFirstPositionInVolume(3); where 3 is the id of the top veto.

The top rectangle is the volume corresponding to the veto.

jgalan · July 12, 2019, 2:55pm

Hi Luis, looking to the methods GetFirstPositionInVolume it seems strange that you get hits outside the volume.

It would be good if you could paste in this post the result of PrintEvent in one of those events that produce the observable outside the volumeID = 3. To determine manually the hit that should be returned.

The top rectangle, volId=3, it is a simple primitive in your GDML definition, or it is the result of more complex volume boolean operations?

lobis · July 14, 2019, 4:46pm

I recreated the simulation and know the positions I get are within bounds. There seemed to be a problem with that simulation in particular

index 230 event ID: 557222 Y pos in vol: 284.113 E dep in vol: 13804.4
index 232 event ID: 562909 Y pos in vol: 311.478 E dep in vol: 11182.2
index 233 event ID: 563658 Y pos in vol: 298.302 E dep in vol: 10847.8
index 240 event ID: 570929 Y pos in vol: 290.13 E dep in vol: 9695.85
index 241 event ID: 571037 Y pos in vol: 309.686 E dep in vol: 10210.2
index 242 event ID: 574147 Y pos in vol: 275 E dep in vol: 10642.3
index 243 event ID: 582437 Y pos in vol: 298.091 E dep in vol: 10438.4
index 244 event ID: 587826 Y pos in vol: 275 E dep in vol: 8815.01
index 245 event ID: 597219 Y pos in vol: 317.53 E dep in vol: 10213.8
index 247 event ID: 597973 Y pos in vol: 136.503 E dep in vol: 0.00212024
Error in <TBasket::Streamer>: The value of fNbytes is incorrect (-698) ; trying to recover by setting it to zero
Error in <TBranchElement::GetBasket>: File: /home/lobis/gitlab/IAXOD0-REST/simulations_data/cosmic_muons/analysis/RUN_ArIso_wall-4000-n1M_lobis_v2.11_0.root at byte:242690348, branch:fTrack.fTrack_ID, entry:248, badread=1, nerrors=1, basketnumber=84
index 248 event ID: 597973 Y pos in vol: 136.503 E dep in vol: 0.00212024
Error in <TBasket::Streamer>: The value of fNbytes is incorrect (-698) ; trying to recover by setting it to zero
Error in <TBranchElement::GetBasket>: File: /home/lobis/gitlab/IAXOD0-REST/simulations_data/cosmic_muons/analysis/RUN_ArIso_wall-4000-n1M_lobis_v2.11_0.root at byte:242690348, branch:fTrack.fTrack_ID, entry:249, badread=1, nerrors=2, basketnumber=84
index 249 event ID: 598317 Y pos in vol: 136.503 E dep in vol: 12646.1
Error in <TBasket::Streamer>: The value of fNbytes is incorrect (-698) ; trying to recover by setting it to zero
Error in <TBranchElement::GetBasket>: File: /home/lobis/gitlab/IAXOD0-REST/simulations_data/cosmic_muons/analysis/RUN_ArIso_wall-4000-n1M_lobis_v2.11_0.root at byte:242690348, branch:fTrack.fTrack_ID, entry:250, badread=1, nerrors=3, basketnumber=84
index 250 event ID: 599164 Y pos in vol: 136.503 E dep in vol: 10376.8
Error in <TBasket::Streamer>: The value of fNbytes is incorrect (-698) ; trying to recover by setting it to zero
Error in <TBranchElement::GetBasket>: File: /home/lobis/gitlab/IAXOD0-REST/simulations_data/cosmic_muons/analysis/RUN_ArIso_wall-4000-n1M_lobis_v2.11_0.root at byte:242690348, branch:fTrack.fTrack_ID, entry:251, badread=1, nerrors=4, basketnumber=84
index 251 event ID: 601480 Y pos in vol: 136.503 E dep in vol: 10079.5
Error in <TBasket::Streamer>: The value of fNbytes is incorrect (-698) ; trying to recover by setting it to zero
Error in <TBranchElement::GetBasket>: File: /home/lobis/gitlab/IAXOD0-REST/simulations_data/cosmic_muons/analysis/RUN_ArIso_wall-4000-n1M_lobis_v2.11_0.root at byte:242690348, branch:fTrack.fTrack_ID, entry:252, badread=1, nerrors=5, basketnumber=84
index 252 event ID: 606522 Y pos in vol: 136.503 E dep in vol: 9295.68
Error in <TBasket::Streamer>: The value of fNbytes is incorrect (-698) ; trying to recover by setting it to zero
Error in <TBranchElement::GetBasket>: File: /home/lobis/gitlab/IAXOD0-REST/simulations_data/cosmic_muons/analysis/RUN_ArIso_wall-4000-n1M_lobis_v2.11_0.root at byte:242690348, branch:fTrack.fTrack_ID, entry:253, badread=1, nerrors=6, basketnumber=84
index 253 event ID: 611994 Y pos in vol: 136.503 E dep in vol: 10376.9
Error in <TBasket::Streamer>: The value of fNbytes is incorrect (-698) ; trying to recover by setting it to zero
Error in <TBranchElement::GetBasket>: File: /home/lobis/gitlab/IAXOD0-REST/simulations_data/cosmic_muons/analysis/RUN_ArIso_wall-4000-n1M_lobis_v2.11_0.root at byte:242690348, branch:fTrack.fTrack_ID, entry:254, badread=1, nerrors=7, basketnumber=84
index 254 event ID: 615212 Y pos in vol: 136.503 E dep in vol: 11288.2
index 255 event ID: 616544 Y pos in vol: 275 E dep in vol: 9238.07
index 256 event ID: 624029 Y pos in vol: 313.25 E dep in vol: 29103
index 258 event ID: 625022 Y pos in vol: 322.503 E dep in vol: 13056
index 259 event ID: 625104 Y pos in vol: 315.214 E dep in vol: 18948.1
index 261 event ID: 628843 Y pos in vol: 275 E dep in vol: 9026.73
index 263 event ID: 633895 Y pos in vol: 276.991 E dep in vol: 19360.4
index 264 event ID: 634105 Y pos in vol: 309.346 E dep in vol: 20156
index 266 event ID: 638497 Y pos in vol: 324.144 E dep in vol: 12532.1
index 269 event ID: 646914 Y pos in vol: 320.688 E dep in vol: 11394.7
index 270 event ID: 653180 Y pos in vol: 315.507 E dep in vol: 12576.6
index 271 event ID: 670345 Y pos in vol: 279.878 E dep in vol: 10855.3
index 272 event ID: 670409 Y pos in vol: 316.106 E dep in vol: 11258.6
index 273 event ID: 674354 Y pos in vol: 276.775 E dep in vol: 11943.6

When iterating over the events I got this warning which I assumed was normal but I guess its a sign that the file got corrupted somehow. The correct simulation was produced with a newer version of REST than the corrupted one so I assume if there was some problem in the code it got fixed by some commit, or maybe its just bad luck.

Anyway after repeating the simulation I got all my first positions in volume within correct values.

However there is something that still puzzles me:

The veto is a rectangular prism situated at y=300 with a width of 50 so compatible y values would be 275-325. I have checked and I got all the y component of the first position in volume within this value but I have seen that instead of getting y’s very close to 325 (muons enter from the top of the veto) I get a more or less uniform distribution of values within 275 and 325, I get many values at exactly 275.

I have checked manually one of this events and I can visually see that the first position where energy was deposited in the veto is close to 325, so I have no idea why I get back the value 275.

jgalan · July 15, 2019, 5:55am

The method TRestG4Event::GetFirstPositionInVolume returns the first registered track that deposits energy on that volume. It might happen that there are several tracks, due to secondaries, and the first track in the event is around 275? Other tracks are then ignored.

I guess to understand it we need to see the details using PrintEvent on that particular event.

lobis · July 15, 2019, 7:34am

*******************************************************
 EVENT ID : 3582633
 TIME : Thu, 01 Jan 1970 00:00:00 +0000 (GMT) +        0 nsec
 SUB-EVENT ID : 0
 Status : OK
*******************************************************
Total energy : 5.611e+05 keV
Sensitive volume energy : 38.49 keV
Source origin : (35.9,355,-12.29) mm
Source 0 Particle name : mu-
Source 0 direction : (-0.138,-0.9902,0.02115)
Source 0 energy : 1.377e+06 keV
Number of active volumes : 10
Active volume 0 has been stored.
Total energy deposit in volume 0 : 38.49 keV
Active volume 1 has been stored.
Total energy deposit in volume 1 : 0 keV
Active volume 2 has been stored.
Total energy deposit in volume 2 : 0 keV
Active volume 3 has been stored.
Total energy deposit in volume 3 : 9077 keV
Active volume 4 has been stored.
Total energy deposit in volume 4 : 1.149e+04 keV
Active volume 5 has been stored.
Total energy deposit in volume 5 : 0 keV
Active volume 6 has been stored.
Total energy deposit in volume 6 : 0 keV
Active volume 7 has been stored.
Total energy deposit in volume 7 : 5.366e+05 keV
Active volume 8 has been stored.
Total energy deposit in volume 8 : 0 keV
Active volume 9 has been stored.
Total energy deposit in volume 9 : 3827 keV
---------------------------------------------------------------------------
Total number of tracks : 16109
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 SubEvent ID : 0 Global timestamp : 0 seconds
 Track ID : 1 Parent ID : 0 Particle : mu- Time track length : 0.0029145 us
 Origin : X = 35.9mm Y = 355mm Z = -12.286mm
 Ekin : 1.3768e+06 keV
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Hit 0 process : TRansportation volume : 3 X : 24.749 Y : 275 Z : -10.588 mm Edep : 9050.4 keV
Hit 1 process : muIoni volume : 7 X : 24.54 Y : 273.49 Z : -10.55 mm Edep : 1714.4 keV
Hit 2 process : muIoni volume : 7 X : 23.864 Y : 268.84 Z : -10.429 mm Edep : 5443.4 keV
Hit 3 process : muIoni volume : 7 X : 22.738 Y : 261.15 Z : -10.247 mm Edep : 9971.5 keV
Hit 4 process : muIoni volume : 7 X : 17.228 Y : 221.48 Z : -8.8873 mm Edep : 49270 keV
Hit 5 process : muIoni volume : 7 X : 12.299 Y : 196.23 Z : -7.9806 mm Edep : 28883 keV
Hit 6 process : muIoni volume : 7 X : 8.8778 Y : 176.76 Z : -7.0234 mm Edep : 22021 keV
Hit 7 process : muIoni volume : 7 X : -2.0795 Y : 130.44 Z : -5.3524 mm Edep : 56776 keV
Hit 8 process : muIoni volume : 7 X : -3.8015 Y : 122 Z : -5.1383 mm Edep : 10196 keV
Hit 9 process : TRansportation volume : 7 X : -12.815 Y : 75 Z : -3.9657 mm Edep : 57794 keV
Hit 10 process : mu-Step volume : 0 X : -17.834 Y : 44.471 Z : -4.8954 mm Edep : 0 keV
Hit 11 process : mu-Step volume : 0 X : -17.849 Y : 44.372 Z : -4.8983 mm Edep : 0 keV
Hit 12 process : mu-Step volume : 0 X : -17.864 Y : 44.273 Z : -4.9013 mm Edep : 0 keV
Hit 13 process : mu-Step volume : 0 X : -17.878 Y : 44.175 Z : -4.9042 mm Edep : 0 keV
Hit 14 process : mu-Step volume : 0 X : -17.893 Y : 44.076 Z : -4.9071 mm Edep : 0.034578 keV
Hit 15 process : mu-Step volume : 0 X : -17.908 Y : 43.977 Z : -4.9101 mm Edep : 0 keV
Hit 16 process : mu-Step volume : 0 X : -17.923 Y : 43.878 Z : -4.913 mm Edep : 0 keV
Hit 17 process : mu-Step volume : 0 X : -17.938 Y : 43.779 Z : -4.9159 mm Edep : 0 keV
Hit 18 process : mu-Step volume : 0 X : -17.953 Y : 43.68 Z : -4.9188 mm Edep : 0 keV
Hit 19 process : mu-Step volume : 0 X : -17.968 Y : 43.581 Z : -4.9218 mm Edep : 0 keV
Hit 20 process : mu-Step volume : 0 X : -17.982 Y : 43.483 Z : -4.9247 mm Edep : 0 keV
Hit 21 process : mu-Step volume : 0 X : -17.997 Y : 43.384 Z : -4.9276 mm Edep : 0 keV
Hit 22 process : mu-Step volume : 0 X : -18.012 Y : 43.285 Z : -4.9306 mm Edep : 0 keV
Hit 23 process : mu-Step volume : 0 X : -18.027 Y : 43.186 Z : -4.9335 mm Edep : 0 keV
Hit 24 process : mu-Step volume : 0 X : -18.042 Y : 43.087 Z : -4.9364 mm Edep : 0 keV
Hit 25 process : mu-Step volume : 0 X : -18.057 Y : 42.988 Z : -4.9394 mm Edep : 0 keV
Hit 26 process : mu-Step volume : 0 X : -18.071 Y : 42.89 Z : -4.9423 mm Edep : 0 keV
Hit 27 process : mu-Step volume : 0 X : -18.086 Y : 42.791 Z : -4.9452 mm Edep : 0.21324 keV
Hit 28 process : mu-Step volume : 0 X : -18.101 Y : 42.692 Z : -4.9481 mm Edep : 0 keV
Hit 29 process : mu-Step volume : 0 X : -18.116 Y : 42.593 Z : -4.9511 mm Edep : 0 keV
Hit 30 process : mu-Step volume : 0 X : -18.131 Y : 42.494 Z : -4.954 mm Edep : 0 keV
Hit 31 process : mu-Step volume : 0 X : -18.146 Y : 42.395 Z : -4.9569 mm Edep : 0 keV
Hit 32 process : mu-Step volume : 0 X : -18.161 Y : 42.296 Z : -4.9599 mm Edep : 0 keV
Hit 33 process : mu-Step volume : 0 X : -18.175 Y : 42.198 Z : -4.9628 mm Edep : 0 keV
Hit 34 process : mu-Step volume : 0 X : -18.19 Y : 42.099 Z : -4.9657 mm Edep : 0.064655 keV
Hit 35 process : mu-Step volume : 0 X : -18.205 Y : 42 Z : -4.9686 mm Edep : 0.079525 keV
Hit 36 process : mu-Step volume : 0 X : -18.22 Y : 41.901 Z : -4.9716 mm Edep : 0 keV
Hit 37 process : mu-Step volume : 0 X : -18.235 Y : 41.802 Z : -4.9745 mm Edep : 0 keV
Hit 38 process : mu-Step volume : 0 X : -18.25 Y : 41.703 Z : -4.9774 mm Edep : 0.031791 keV
Hit 39 process : mu-Step volume : 0 X : -18.264 Y : 41.604 Z : -4.9804 mm Edep : 0 keV
Hit 40 process : mu-Step volume : 0 X : -18.279 Y : 41.506 Z : -4.9833 mm Edep : 0 keV
Hit 41 process : mu-Step volume : 0 X : -18.294 Y : 41.407 Z : -4.9862 mm Edep : 0 keV
Hit 42 process : mu-Step volume : 0 X : -18.309 Y : 41.308 Z : -4.9892 mm Edep : 0 keV
Hit 43 process : mu-Step volume : 0 X : -18.324 Y : 41.209 Z : -4.9921 mm Edep : 0 keV
Hit 44 process : mu-Step volume : 0 X : -18.339 Y : 41.11 Z : -4.995 mm Edep : 0 keV
Hit 45 process : mu-Step volume : 0 X : -18.354 Y : 41.011 Z : -4.9979 mm Edep : 0 keV
Hit 46 process : mu-Step volume : 0 X : -18.368 Y : 40.913 Z : -5.0009 mm Edep : 0 keV
Hit 47 process : mu-Step volume : 0 X : -18.383 Y : 40.814 Z : -5.0038 mm Edep : 0 keV
Hit 48 process : mu-Step volume : 0 X : -18.398 Y : 40.715 Z : -5.0067 mm Edep : 0 keV
Hit 49 process : mu-Step volume : 0 X : -18.413 Y : 40.616 Z : -5.0097 mm Edep : 0 keV
Hit 50 process : mu-Step volume : 0 X : -18.428 Y : 40.517 Z : -5.0126 mm Edep : 0 keV
Hit 51 process : mu-Step volume : 0 X : -18.443 Y : 40.418 Z : -5.0155 mm Edep : 0 keV
Hit 52 process : mu-Step volume : 0 X : -18.457 Y : 40.319 Z : -5.0185 mm Edep : 0.21032 keV
Hit 53 process : mu-Step volume : 0 X : -18.472 Y : 40.221 Z : -5.0214 mm Edep : 0.27068 keV
Hit 54 process : mu-Step volume : 0 X : -18.487 Y : 40.122 Z : -5.0243 mm Edep : 0.015967 keV
Hit 55 process : mu-Step volume : 0 X : -18.502 Y : 40.023 Z : -5.0273 mm Edep : 0 keV
Hit 56 process : mu-Step volume : 0 X : -18.517 Y : 39.924 Z : -5.0302 mm Edep : 0 keV
Hit 57 process : mu-Step volume : 0 X : -18.532 Y : 39.825 Z : -5.0331 mm Edep : 0 keV
Hit 58 process : mu-Step volume : 0 X : -18.547 Y : 39.726 Z : -5.036 mm Edep : 0 keV
Hit 59 process : mu-Step volume : 0 X : -18.561 Y : 39.628 Z : -5.039 mm Edep : 0 keV
Hit 60 process : mu-Step volume : 0 X : -18.576 Y : 39.529 Z : -5.0419 mm Edep : 0 keV
Hit 61 process : mu-Step volume : 0 X : -18.591 Y : 39.43 Z : -5.0448 mm Edep : 0.055844 keV
Hit 62 process : mu-Step volume : 0 X : -18.606 Y : 39.331 Z : -5.0478 mm Edep : 0 keV
Hit 63 process : mu-Step volume : 0 X : -18.621 Y : 39.232 Z : -5.0507 mm Edep : 0 keV
Hit 64 process : mu-Step volume : 0 X : -18.636 Y : 39.133 Z : -5.0536 mm Edep : 0.030316 keV
Hit 65 process : mu-Step volume : 0 X : -18.65 Y : 39.034 Z : -5.0566 mm Edep : 0 keV
Hit 66 process : mu-Step volume : 0 X : -18.665 Y : 38.936 Z : -5.0595 mm Edep : 0.026102 keV
Hit 67 process : mu-Step volume : 0 X : -18.68 Y : 38.837 Z : -5.0624 mm Edep : 0 keV
Hit 68 process : mu-Step volume : 0 X : -18.695 Y : 38.738 Z : -5.0654 mm Edep : 0 keV
Hit 69 process : mu-Step volume : 0 X : -18.71 Y : 38.639 Z : -5.0683 mm Edep : 0 keV
Hit 70 process : mu-Step volume : 0 X : -18.725 Y : 38.54 Z : -5.0712 mm Edep : 0 keV
Hit 71 process : mu-Step volume : 0 X : -18.74 Y : 38.441 Z : -5.0741 mm Edep : 0.01576 keV
Hit 72 process : mu-Step volume : 0 X : -18.754 Y : 38.342 Z : -5.0771 mm Edep : 0 keV
Hit 73 process : mu-Step volume : 0 X : -18.769 Y : 38.244 Z : -5.08 mm Edep : 0 keV
Hit 74 process : mu-Step volume : 0 X : -18.784 Y : 38.145 Z : -5.0829 mm Edep : 0 keV
Hit 75 process : mu-Step volume : 0 X : -18.799 Y : 38.046 Z : -5.0859 mm Edep : 0.088732 keV
Hit 76 process : mu-Step volume : 0 X : -18.814 Y : 37.947 Z : -5.0888 mm Edep : 0 keV
Hit 77 process : mu-Step volume : 0 X : -18.829 Y : 37.848 Z : -5.0917 mm Edep : 0 keV
Hit 78 process : mu-Step volume : 0 X : -18.844 Y : 37.749 Z : -5.0947 mm Edep : 0 keV
Hit 79 process : mu-Step volume : 0 X : -18.858 Y : 37.651 Z : -5.0976 mm Edep : 0.79189 keV
Hit 80 process : mu-Step volume : 0 X : -18.873 Y : 37.552 Z : -5.1005 mm Edep : 0 keV
Hit 81 process : mu-Step volume : 0 X : -18.888 Y : 37.453 Z : -5.1035 mm Edep : 0.078219 keV
Hit 82 process : mu-Step volume : 0 X : -18.903 Y : 37.354 Z : -5.1064 mm Edep : 0.021643 keV
Hit 83 process : mu-Step volume : 0 X : -18.918 Y : 37.255 Z : -5.1093 mm Edep : 0 keV
Hit 84 process : mu-Step volume : 0 X : -18.933 Y : 37.156 Z : -5.1123 mm Edep : 0.075441 keV
Hit 85 process : mu-Step volume : 0 X : -18.947 Y : 37.057 Z : -5.1152 mm Edep : 0 keV
Hit 86 process : mu-Step volume : 0 X : -18.962 Y : 36.959 Z : -5.1181 mm Edep : 0 keV
Hit 87 process : mu-Step volume : 0 X : -18.977 Y : 36.86 Z : -5.1211 mm Edep : 0 keV
Hit 88 process : mu-Step volume : 0 X : -18.992 Y : 36.761 Z : -5.124 mm Edep : 0 keV
Hit 89 process : mu-Step volume : 0 X : -19.007 Y : 36.662 Z : -5.1269 mm Edep : 0 keV
Hit 90 process : mu-Step volume : 0 X : -19.022 Y : 36.563 Z : -5.1299 mm Edep : 0.044557 keV
Hit 91 process : mu-Step volume : 0 X : -19.037 Y : 36.464 Z : -5.1328 mm Edep : 0 keV
Hit 92 process : mu-Step volume : 0 X : -19.051 Y : 36.366 Z : -5.1357 mm Edep : 0.12433 keV
Hit 93 process : mu-Step volume : 0 X : -19.066 Y : 36.267 Z : -5.1387 mm Edep : 0 keV
Hit 94 process : mu-Step volume : 0 X : -19.081 Y : 36.168 Z : -5.1416 mm Edep : 0.026623 keV
Hit 95 process : mu-Step volume : 0 X : -19.096 Y : 36.069 Z : -5.1445 mm Edep : 0 keV
Hit 96 process : mu-Step volume : 0 X : -19.111 Y : 35.97 Z : -5.1475 mm Edep : 0 keV
Hit 97 process : mu-Step volume : 0 X : -19.126 Y : 35.871 Z : -5.1504 mm Edep : 0 keV
Hit 98 process : mu-Step volume : 0 X : -19.14 Y : 35.772 Z : -5.1533 mm Edep : 0.064288 keV
Hit 99 process : mu-Step volume : 0 X : -19.155 Y : 35.674 Z : -5.1563 mm Edep : 0 keV
Hit 100 process : mu-Step volume : 0 X : -19.17 Y : 35.575 Z : -5.1592 mm Edep : 0 keV
Hit 101 process : mu-Step volume : 0 X : -19.185 Y : 35.476 Z : -5.1621 mm Edep : 0 keV
Hit 102 process : mu-Step volume : 0 X : -19.2 Y : 35.377 Z : -5.1651 mm Edep : 0 keV
Hit 103 process : mu-Step volume : 0 X : -19.215 Y : 35.278 Z : -5.168 mm Edep : 0 keV
Hit 104 process : mu-Step volume : 0 X : -19.23 Y : 35.179 Z : -5.1709 mm Edep : 0 keV
Hit 105 process : mu-Step volume : 0 X : -19.244 Y : 35.081 Z : -5.1739 mm Edep : 0 keV
Hit 106 process : mu-Step volume : 0 X : -19.259 Y : 34.982 Z : -5.1768 mm Edep : 0 keV
Hit 107 process : mu-Step volume : 0 X : -19.274 Y : 34.883 Z : -5.1797 mm Edep : 0 keV
Hit 108 process : mu-Step volume : 0 X : -19.289 Y : 34.784 Z : -5.1827 mm Edep : 0 keV
Hit 109 process : mu-Step volume : 0 X : -19.304 Y : 34.685 Z : -5.1856 mm Edep : 0 keV
Hit 110 process : mu-Step volume : 0 X : -19.319 Y : 34.586 Z : -5.1886 mm Edep : 0 keV
Hit 111 process : mu-Step volume : 0 X : -19.333 Y : 34.487 Z : -5.1915 mm Edep : 0.083928 keV
Hit 112 process : mu-Step volume : 0 X : -19.348 Y : 34.389 Z : -5.1944 mm Edep : 0 keV
Hit 113 process : mu-Step volume : 0 X : -19.363 Y : 34.29 Z : -5.1974 mm Edep : 0 keV
Hit 114 process : mu-Step volume : 0 X : -19.378 Y : 34.191 Z : -5.2003 mm Edep : 0 keV
Hit 115 process : mu-Step volume : 0 X : -19.393 Y : 34.092 Z : -5.2032 mm Edep : 0.044144 keV
Hit 116 process : mu-Step volume : 0 X : -19.408 Y : 33.993 Z : -5.2062 mm Edep : 0 keV
Hit 117 process : mu-Step volume : 0 X : -19.423 Y : 33.894 Z : -5.2091 mm Edep : 0 keV
Hit 118 process : mu-Step volume : 0 X : -19.437 Y : 33.795 Z : -5.212 mm Edep : 0 keV
Hit 119 process : mu-Step volume : 0 X : -19.452 Y : 33.697 Z : -5.215 mm Edep : 0 keV
Hit 120 process : mu-Step volume : 0 X : -19.467 Y : 33.598 Z : -5.2179 mm Edep : 0 keV
Hit 121 process : mu-Step volume : 0 X : -19.482 Y : 33.499 Z : -5.2208 mm Edep : 0 keV
Hit 122 process : mu-Step volume : 0 X : -19.497 Y : 33.4 Z : -5.2238 mm Edep : 0 keV
Hit 123 process : mu-Step volume : 0 X : -19.512 Y : 33.301 Z : -5.2267 mm Edep : 0 keV
Hit 124 process : mu-Step volume : 0 X : -19.527 Y : 33.202 Z : -5.2296 mm Edep : 0 keV
Hit 125 process : mu-Step volume : 0 X : -19.541 Y : 33.104 Z : -5.2326 mm Edep : 0.018657 keV
Hit 126 process : mu-Step volume : 0 X : -19.556 Y : 33.005 Z : -5.2355 mm Edep : 0.151 keV
Hit 127 process : mu-Step volume : 0 X : -19.571 Y : 32.906 Z : -5.2384 mm Edep : 0 keV
Hit 128 process : mu-Step volume : 0 X : -19.586 Y : 32.807 Z : -5.2414 mm Edep : 0 keV
Hit 129 process : mu-Step volume : 0 X : -19.601 Y : 32.708 Z : -5.2443 mm Edep : 0 keV
Hit 130 process : mu-Step volume : 0 X : -19.616 Y : 32.609 Z : -5.2473 mm Edep : 0 keV
Hit 131 process : mu-Step volume : 0 X : -19.63 Y : 32.51 Z : -5.2502 mm Edep : 0 keV
Hit 132 process : mu-Step volume : 0 X : -19.645 Y : 32.412 Z : -5.2531 mm Edep : 0.25534 keV
Hit 133 process : mu-Step volume : 0 X : -19.66 Y : 32.313 Z : -5.2561 mm Edep : 0 keV
Hit 134 process : mu-Step volume : 0 X : -19.675 Y : 32.214 Z : -5.259 mm Edep : 0 keV
Hit 135 process : mu-Step volume : 0 X : -19.69 Y : 32.115 Z : -5.2619 mm Edep : 0 keV
Hit 136 process : mu-Step volume : 0 X : -19.705 Y : 32.016 Z : -5.2649 mm Edep : 0.15654 keV
Hit 137 process : mu-Step volume : 0 X : -19.72 Y : 31.917 Z : -5.2678 mm Edep : 0 keV
Hit 138 process : mu-Step volume : 0 X : -19.734 Y : 31.819 Z : -5.2708 mm Edep : 0 keV
Hit 139 process : mu-Step volume : 0 X : -19.749 Y : 31.72 Z : -5.2737 mm Edep : 0 keV
Hit 140 process : mu-Step volume : 0 X : -19.764 Y : 31.621 Z : -5.2766 mm Edep : 0 keV
Hit 141 process : mu-Step volume : 0 X : -19.779 Y : 31.522 Z : -5.2796 mm Edep : 0 keV
Hit 142 process : mu-Step volume : 0 X : -19.794 Y : 31.423 Z : -5.2825 mm Edep : 0 keV
Hit 143 process : mu-Step volume : 0 X : -19.809 Y : 31.324 Z : -5.2854 mm Edep : 0 keV
Hit 144 process : mu-Step volume : 0 X : -19.823 Y : 31.225 Z : -5.2884 mm Edep : 0.024715 keV
Hit 145 process : mu-Step volume : 0 X : -19.838 Y : 31.127 Z : -5.2913 mm Edep : 0 keV
Hit 146 process : mu-Step volume : 0 X : -19.853 Y : 31.028 Z : -5.2942 mm Edep : 0 keV
Hit 147 process : mu-Step volume : 0 X : -19.868 Y : 30.929 Z : -5.2972 mm Edep : 0.028604 keV
Hit 148 process : mu-Step volume : 0 X : -19.883 Y : 30.83 Z : -5.3001 mm Edep : 0 keV
Hit 149 process : mu-Step volume : 0 X : -19.898 Y : 30.731 Z : -5.303 mm Edep : 0 keV
Hit 150 process : mu-Step volume : 0 X : -19.913 Y : 30.632 Z : -5.306 mm Edep : 0 keV
Hit 151 process : mu-Step volume : 0 X : -19.927 Y : 30.533 Z : -5.3089 mm Edep : 0 keV
Hit 152 process : mu-Step volume : 0 X : -19.942 Y : 30.435 Z : -5.3118 mm Edep : 0 keV
Hit 153 process : mu-Step volume : 0 X : -19.957 Y : 30.336 Z : -5.3148 mm Edep : 0 keV
Hit 154 process : mu-Step volume : 0 X : -19.972 Y : 30.237 Z : -5.3177 mm Edep : 0 keV
Hit 155 process : mu-Step volume : 0 X : -19.987 Y : 30.138 Z : -5.3206 mm Edep : 0 keV
Hit 156 process : mu-Step volume : 0 X : -20.002 Y : 30.039 Z : -5.3236 mm Edep : 0 keV
Hit 157 process : mu-Step volume : 0 X : -20.016 Y : 29.94 Z : -5.3265 mm Edep : 0.15779 keV
Hit 158 process : mu-Step volume : 0 X : -20.031 Y : 29.842 Z : -5.3294 mm Edep : 0 keV
Hit 159 process : mu-Step volume : 0 X : -20.046 Y : 29.743 Z : -5.3324 mm Edep : 0 keV
Hit 160 process : mu-Step volume : 0 X : -20.061 Y : 29.644 Z : -5.3353 mm Edep : 0 keV
Hit 161 process : mu-Step volume : 0 X : -20.076 Y : 29.545 Z : -5.3382 mm Edep : 0 keV
Hit 162 process : mu-Step volume : 0 X : -20.091 Y : 29.446 Z : -5.3412 mm Edep : 0 keV
Hit 163 process : mu-Step volume : 0 X : -20.106 Y : 29.347 Z : -5.3441 mm Edep : 0 keV
Hit 164 process : mu-Step volume : 0 X : -20.12 Y : 29.248 Z : -5.347 mm Edep : 0.053267 keV
Hit 165 process : mu-Step volume : 0 X : -20.135 Y : 29.15 Z : -5.35 mm Edep : 0 keV
Hit 166 process : mu-Step volume : 0 X : -20.15 Y : 29.051 Z : -5.3529 mm Edep : 0.025637 keV
Hit 167 process : mu-Step volume : 0 X : -20.165 Y : 28.952 Z : -5.3558 mm Edep : 0 keV
Hit 168 process : mu-Step volume : 0 X : -20.18 Y : 28.853 Z : -5.3587 mm Edep : 0 keV
Hit 169 process : mu-Step volume : 0 X : -20.195 Y : 28.754 Z : -5.3617 mm Edep : 0 keV
Hit 170 process : mu-Step volume : 0 X : -20.209 Y : 28.655 Z : -5.3646 mm Edep : 0.016285 keV
Hit 171 process : mu-Step volume : 0 X : -20.224 Y : 28.557 Z : -5.3675 mm Edep : 0 keV
Hit 172 process : mu-Step volume : 0 X : -20.239 Y : 28.458 Z : -5.3705 mm Edep : 0 keV
Hit 173 process : mu-Step volume : 0 X : -20.254 Y : 28.359 Z : -5.3734 mm Edep : 0 keV
Hit 174 process : mu-Step volume : 0 X : -20.269 Y : 28.26 Z : -5.3763 mm Edep : 0 keV
Hit 175 process : mu-Step volume : 0 X : -20.284 Y : 28.161 Z : -5.3793 mm Edep : 0 keV
Hit 176 process : mu-Step volume : 0 X : -20.299 Y : 28.062 Z : -5.3822 mm Edep : 0 keV
Hit 177 process : mu-Step volume : 0 X : -20.313 Y : 27.963 Z : -5.3851 mm Edep : 0 keV
Hit 178 process : mu-Step volume : 0 X : -20.328 Y : 27.865 Z : -5.388 mm Edep : 0 keV
Hit 179 process : mu-Step volume : 0 X : -20.343 Y : 27.766 Z : -5.391 mm Edep : 0 keV
Hit 180 process : mu-Step volume : 0 X : -20.358 Y : 27.667 Z : -5.3939 mm Edep : 0 keV
Hit 181 process : mu-Step volume : 0 X : -20.373 Y : 27.568 Z : -5.3968 mm Edep : 0 keV
Hit 182 process : mu-Step volume : 0 X : -20.388 Y : 27.469 Z : -5.3998 mm Edep : 0 keV
Hit 183 process : mu-Step volume : 0 X : -20.402 Y : 27.37 Z : -5.4027 mm Edep : 0 keV
Hit 184 process : mu-Step volume : 0 X : -20.417 Y : 27.272 Z : -5.4056 mm Edep : 0 keV

....
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

Here is the top part of the log (its very long, hits character limit). Only in the first line there is a “hit” regarding volume 3 (the veto of interest).

Hit 0 process : TRansportation volume : 3 X : 24.749 Y : 275 Z : -10.588 mm Edep : 9050.4 keV

I will take a look at the code. I guess when it says “hit 0 process” its really meaning track 0? edit: nevermind, its printing info about all hits on that track (which is really really long).

lobis · July 15, 2019, 7:58am

I don’t really understand it. I have taken a look at track 0, where almost all of the energy is deposited in the veto (9050.4/9077 keV). There is only one hit (hit 0) that leaves energy in the track but on the picture its clear to me there are multiple hits in that volume (unless we don’t agree on what a “hit” is). Maybe it has to do with the type of process “TRansportation”?

lobis · July 15, 2019, 8:10am

Okay, I think now I understand what happens.

In the image all those pink crosses are electrons, which are tracks on their own, stemming directly from the main muon track. ID of subtracks (if this is the proper term) are added backwards, so the first subtrack (the first electron in the muon veto) ID is in this case the highest one (16109). This track has the correct position (around y=325) but it leaves no energy in the volume/ no hits.

edit: however there are clearly energy deposits near y=325 so maybe this is what we should look for.

This is the information I am looking for, I guess there is no straight forward way to access it and it has to be implemented.

edit: Using GetLastPositionFromVolume gives me a more appropiate result, this has to do with the reverse ordering of subtracks, is this intentional?

At the moment (v2.2.11_dev@b5a56035) if I plot the points using “GetLastPositionInVolume” I get what I would expect so maybe we should either reverse the order in which tracks are stored or swap GetLastPosition with GetFirstPosition methods @jgalan?

jgalan · July 15, 2019, 5:12pm

I think here there is some confusion, because in fact, in the drawing, the pink markers are the position of each new G4Track, and not a hit (or energy deposit). Hits are drawn in yellow.

These methods were introduced by @luzon, perhaps we should identify the tracks by the lowest track ID? For the first hit? Or the hit belonging to the highest energy track? main track? I don’t know which is the most generic way to do this.

In this particular case, it is clear that if I identify the muon I will get a good answer by getting the first hit?

Then, it will be better to attach an output file.