Fork me on GitHub

Opened 5 years ago

Last modified 4 years ago

#1444 new How to

change default mass assumption for tracks, implement PID

Reported by: Leticia Cunqueiro Owned by:
Priority: major Milestone:
Component: Delphes code Version: Delphes 3
Keywords: PID, mass assumption for tracks Cc:

Description

Hello,
if I understand correctly the mass assumption for tracks is zero. Is that right? How can I change that to pion mass as default?
On top of that, I would like to implement some PID. If I have a proton that produces a track, I would like to assign the mass of the proton to the track with a given probability (and some missidentification probability too). How can this be done? Thanks a lot for your help,
Leticia

Change History (1)

comment:1 by Jaredb027, 4 years ago

Hello,

I am also interested in assigning the pion mass to the track as well. I have identified where the track mass is defined in the Track::P4() method inside DelphesClasses.cc (​https://github.com/delphes/delphes/blob/master/classes/DelphesClasses.cc#L102)

TLorentzVector Track::P4() const
{
  TLorentzVector vec;
  vec.SetPtEtaPhiM(PT, Eta, Phi, 0.0);
  return vec;
}

My thoughts at a solution would be to replace the 0.0 with the pion mass (139.57 GeV). I believe this is the correct way to change the pion mass.

However, I haven't been able to get it to work to know for certain. I am having issues with the event generation after implementing this, though I am not sure it is the result of this change.

Note: See TracTickets for help on using tickets.