@@ -24,7 +24,7 @@ struct HFMuonFwdTracks {
2424 if (!muonTracksOut.is_open ()) {
2525 LOGF (fatal, " Failed to open muonTracks.csv for writing" );
2626 }
27- muonTracksOut << " trackID,trackType,phi,tgl,signed1Pt,nClusters,pDCA,rAtAbsorberEnd,sign,chi2,chi2MatchMCHMID,chi2MatchMCHMFT,trackTime,eta,pt,p" << std::endl;
27+ muonTracksOut << " trackID,trackType,firstMotherPDG, phi,tgl,signed1Pt,nClusters,pDCA,rAtAbsorberEnd,sign,chi2,chi2MatchMCHMID,chi2MatchMCHMFT,trackTime,eta,pt,p" << std::endl;
2828
2929 // define axes you want to use
3030 const AxisSpec axisCounter{1 , 0 , +1 , " " };
@@ -124,10 +124,10 @@ struct HFMuonFwdTracks {
124124 if (muTrackType == 3 ) histos.fill (HIST (" muPtHistReco" ), muPt); // look at pT for standalone tracks
125125
126126 // save all muon tracks to the output file
127- muonTracksOut << recoTrackID << " ," << muTrackType << " ," << track. phi () << " ," << track.tgl () << " ," << track.signed1Pt () << " ,"
128- << static_cast < int64_t >( track.nClusters ()) << " ," << muDca << " ," << track. rAtAbsorberEnd () << " ,"
129- << static_cast < int64_t >( track.sign ()) << " ," << muChi2 << " ," << muChi2MatchMCHMID << " ," << muChi2MatchMCHMFT << " ,"
130- << track.trackTime () << " ," << muEta << " ," << muPt << " ," << track.p () << std::endl;
127+ muonTracksOut << recoTrackID << " ," << muTrackType << " ," << muMotherPDG << " ," << track.phi () << " ," << track.tgl () << " ,"
128+ << track.signed1Pt () << " ," << static_cast < int64_t >(track. nClusters ()) << " ," << muDca << " ,"
129+ << track.rAtAbsorberEnd () << " ," << static_cast < int64_t >(track. sign ()) << " ," << muChi2 << " ," << muChi2MatchMCHMID << " ,"
130+ << muChi2MatchMCHMFT << " , " << track.trackTime () << " ," << muEta << " ," << muPt << " ," << track.p () << std::endl;
131131 };
132132
133133 // check whether is prompt HF hadron
0 commit comments