mc_truth/interactions

[YifanC]

Fill mc_truth/interactions with vertices when there is no mc_hdr available.

[cuddandr]

In raw_event_generator.py a couple warning messages are added, but one of them uses print() and the other uses warning.warn(). Use print() or warn() for both.

Also, what is the purpose of switching between file_vertex_id and vertex_id? It looks like it's there as a fallback if the file_ version is not present in the dataset.

[YifanC]

@cuddandr addressed the warning comment.
For the "file_vertex_id" one, please also check this PR: DUNE/larnd-sim#234
An interaction is considered as a unit in raw_event_generator. It is written with the assumption that vertex_id is unique. If we pass non-unique vertex here, the event=interaction selection would be wrong. (

ndlar_flow/src/proto_nd_flow/reco/charge/raw_event_generator.py

Line 275 in 9800a96

intr_evid = self.mc_events[self.vertex_id_name][:]

) (It is not so critical to have unique traj_id I think..)
In particle bomb adaption in edepsim, the vertex_id restart from 0 every event and the traj_id restart from 0 every interaction. In the beam simulation vertex_id is unique. I have not looked at the cosmic simulation yet.