Skip to content

TM MC pipeline#452

Open
mghrear wants to merge 9 commits into
masterfrom
TM_for_EP
Open

TM MC pipeline#452
mghrear wants to merge 9 commits into
masterfrom
TM_for_EP

Conversation

@mghrear

@mghrear mghrear commented Jul 8, 2026

Copy link
Copy Markdown
Collaborator

Adding TM MC pipeline. This starts with lhe files from Emrys' generator. Displacement is given by 'ctau' which is included as a parameter in job.json and applied via hpsmc.tools.DisplaceTime. I did this a few months ago but don't want it to get buried as it may be useful in the future.

mghrear and others added 3 commits July 8, 2026 11:50
Co-authored-by: Emrys Peets <epeets@stanford.edu>
lcio_bin_jar = /home/groups/laurenat/majd/HPS/hps-lcio/target/lcio-2.7.4-SNAPSHOT-bin.jar

[EvioToLcio]
hps_java_bin_jar = /home/groups/laurenat/majd/HPS/hps-java/distribution/target/hps-distribution-5.2.2-SNAPSHOT-bin.jar

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Here and in many files below, there are a lot of user paths. Are these needing to be committed?

Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'll remove the file paths and leave a generic .hpsmc file.

/scratch/groups/laurenat/2021/v9_pass5/TM/slic/phiKK_3pt74_rot_97.slcio
/scratch/groups/laurenat/2021/v9_pass5/TM/slic/phiKK_3pt74_rot_98.slcio
/scratch/groups/laurenat/2021/v9_pass5/TM/slic/phiKK_3pt74_rot_99.slcio
/scratch/groups/laurenat/2021/v9_pass5/TM/slic/phiKK_3pt74_rot_9.slcio

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

These files needeD? Also under a TM directory but labelled as phiKK?

Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

They are only needed if someone wants to process my specific TM SLIC simulations. I can remove the signal_pulser_overlay_to recon subdirectory all together, this just applies the common pulser overlay job.

"detector": "{{ detector }}",
"output_files": {
"phi_rot.stdhep": "stdhep/phiKK_{{ run_params }}_rot_{{ job_id }}.stdhep",
"phi_rot.slcio": "slic/phiKK_{{ run_params }}_rot_{{ job_id }}.slcio"

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I guess that naming comes from this. Harmless, of course, but maybe worth updating?

Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yes, it harmless but should be updated. I'll updated it here and other places this shows up.

Comment thread tools/stdhep-tools/src/lhe_phi.cc Outdated
/**
* Takes input stdhep file, applies beam rotation and width,
* and writes to a new stdhep file
* /todo is this description correct?

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is it correct? I don't think so. Looks like it shuffles around PDG IDs?

Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Sorry, lhe_phi.cc is actually not relevant to this PR. It's a slightly modified version of lhe_tridents.cc that I used for my phi->KK simulations. I'll take it out of this PR.

Comment thread tools/stdhep-tools/src/lhe_phi.cc Outdated
vector<stdhep_entry> new_event;

int rseed = 0;
double decay_length = 0.0;

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is decay length ever used beyond a print out?

Comment thread tools/stdhep-tools/src/lhe_phi.cc Outdated
open_write(argv[optind+1], ostream, 1000);
nevhep = 1;

printf("Applying decay length of %f mm\n", decay_length > 0 ? decay_length : 0.0);

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Where does this get applied below?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants