New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Vertex generation used with BHgen #192
Comments
There is a vertex volume option in MCwrapper…..may be able to leverage it:
VERTEX_VOLUME=ccdb:29.5 #defines the VERTEX card in geant. The first value should be a comma seperated string x,y,z,var_xx,var_xy,var_yy,dxdz,dydz and the second part is the length (only the length need be given)
Thomas Britton
Staff Scientist, Scientific Computing
Jefferson Lab
From: ***@***.***>
Sent: Wednesday, July 7, 2021 2:47 PM
To: ***@***.***>
Cc: ***@***.***>
Subject: [EXTERNAL] [JeffersonLab/HDGeant4] Vertex generation used with BHgen (#192)
We have realized that BHgen currently generates pairs only in a defined hydrogen target. This needs to be expanded for CPP, since this experiment does not have a hydrogen target. The CPP target is a solid 208Pb target placed at z=1cm. (The entire target and start counter are removed and reconfigured to hold the solid target —still to be designed). Therefore, no pairs are currently generated if we have the CPP geometry defined via the environment variable JANA_GEOMETRY_URL. The next step is to make the system more flexible, particularly for the CPP experiment.
Possible solutions: There are at least two things to consider:
1. Vertex generation. One could use the command cards in the control.in file that are used for the particle gun. The SCAP and TGTWIDTH parameters seem to be appropriate and are already operational. Defining these separately from the geometry file is useful, since for testing and for systematic studies one may wish to generate events that are not uniformly distributed along the target. For example, one might wish to study events generated from windows along the beamline.
2. We also talked about updating the generator to use nuclear targets (i.e. 208Pb specifically). One possibility would be to dig into the defined geometry and use materials specified there. However, that might limit flexibility and need to make additional assumptions. Another alternative would be to add a new command card, or add parameters to one of the above commands. In general I favor having the event generators to be as self-contained as possible.
—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub<https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_JeffersonLab_HDGeant4_issues_192&d=DwMFaQ&c=CJqEzB1piLOyyvZjb8YUQw&r=ccFffx721N71hPpKcJGvJIqY9RM4gBTuzp9ir7rze5Q&m=ZwWBwa1Rx0FAWMsdbV3hm5vGRA75beCujh2_dJ-Slds&s=IZTELzPk9zkzfMcUy_HutLPZgKG_Pxe5gxAwFNQRJQo&e=>, or unsubscribe<https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_AFRO2BGBOCU3YAV4GPUMOOTTWSOMZANCNFSM477G4XLQ&d=DwMFaQ&c=CJqEzB1piLOyyvZjb8YUQw&r=ccFffx721N71hPpKcJGvJIqY9RM4gBTuzp9ir7rze5Q&m=ZwWBwa1Rx0FAWMsdbV3hm5vGRA75beCujh2_dJ-Slds&s=pbkZuyA14k6nxUkx-K1LVEwTg4z4IZooQMDLNSICSjI&e=>.
|
Hello Elton,
If you have a set of predefined volume names for your targets, with
assigned materials, I can add these to the internal pairs generator and
select the appropriate atomic/nuclear form factors to go with that target.
It would be helpful if you could provide me with ready-made c++ methods for
these two form factors so I do not have to develop my own and test them.
…-Richard Jones
On Wed, Jul 7, 2021 at 2:47 PM eltonssmith ***@***.***> wrote:
We have realized that BHgen currently generates pairs only in a defined
hydrogen target. This needs to be expanded for CPP, since this experiment
does not have a hydrogen target. The CPP target is a solid 208Pb target
placed at z=1cm. (The entire target and start counter are removed and
reconfigured to hold the solid target —still to be designed). Therefore, no
pairs are currently generated if we have the CPP geometry defined via the
environment variable JANA_GEOMETRY_URL. The next step is to make the system
more flexible, particularly for the CPP experiment.
Possible solutions: There are at least two things to consider:
1. Vertex generation. One could use the command cards in the control.in
file that are used for the particle gun. The SCAP and TGTWIDTH parameters
seem to be appropriate and are already operational. Defining these
separately from the geometry file is useful, since for testing and for
systematic studies one may wish to generate events that are not uniformly
distributed along the target. For example, one might wish to study events
generated from windows along the beamline.
2. We also talked about updating the generator to use nuclear targets
(i.e. 208Pb specifically). One possibility would be to dig into the defined
geometry and use materials specified there. However, that might limit
flexibility and need to make additional assumptions. Another alternative
would be to add a new command card, or add parameters to one of the above
commands. In general I favor having the event generators to be as
self-contained as possible.
—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub
<#192>, or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AB3YKWHXNLFX5MHNQBX6ID3TWSOMPANCNFSM477G4XLQ>
.
|
Richard, I assume the information in the ccdb xml file contains all the information you need on the CPP target: /GEOMETRY> cat TargetCPP_HDDS.xml | |
Accidentally closed issue, so I reopened it. |
|
Elton, perfect, thanks!
-Richard
…On Wed, Jul 14, 2021 at 10:51 AM eltonssmith ***@***.***> wrote:
Richard, I assume the information in the ccdb xml file contains all the
information you need on the CPP target:
/GEOMETRY> cat TargetCPP_HDDS.xml
+------------+
| xml |
| string |
+------------+
|
|
+------------+
/GEOMETRY>
—
You are receiving this because you commented.
Reply to this email directly, view it on GitHub
<#192 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AB3YKWCLYOPPPGBN3VMMOK3TXWP7TANCNFSM477G4XLQ>
.
|
The latest version of hdgeant4 (2.26.0) can now generate electron and muon pairs from the CPP vertex. Thanks to Richard! Dump for four events using GENBEAM 'BHgen' 0.5 event= 1 beam type= Photon E= 9.89018821716 px= 0.000307428475935 py= -8.97805512068e-05 pz= 9.89018821716 event= 2 beam type= Photon E= 8.39611148834 px= 0.000193996616872 py= 2.03653107747e-05 pz= 8.39611148834 event= 3 beam type= Photon E= 6.4943614006 px= 0.000148683568113 py= 0.000158074893989 pz= 6.4943614006 event= 4 beam type= Photon E= 9.87286281586 px= -0.000124752055854 py= -0.000148085571709 pz= 9.87286281586 |
We have realized that BHgen currently generates pairs only in a defined hydrogen target. This needs to be expanded for CPP, since this experiment does not have a hydrogen target. The CPP target is a solid 208Pb target placed at z=1cm. (The entire target and start counter are removed and reconfigured to hold the solid target —still to be designed). Therefore, no pairs are currently generated if we have the CPP geometry defined via the environment variable JANA_GEOMETRY_URL. The next step is to make the system more flexible, particularly for the CPP experiment.
Possible solutions: There are at least two things to consider:
The text was updated successfully, but these errors were encountered: