3737namespace dimu
3838{
3939// dimuon
40+ DECLARE_SOA_COLUMN (RunNumber, runNumber, int );
4041DECLARE_SOA_COLUMN (M, m, float );
4142DECLARE_SOA_COLUMN (Energy, energy, float );
4243DECLARE_SOA_COLUMN (Px, px, float );
@@ -73,6 +74,7 @@ DECLARE_SOA_COLUMN(Nclass, nclass, int);
7374namespace o2 ::aod
7475{
7576DECLARE_SOA_TABLE (DiMu, " AOD" , " DIMU" ,
77+ dimu::RunNumber,
7678 dimu::M, dimu::Energy, dimu::Px, dimu::Py, dimu::Pz, dimu::Pt, dimu::Rap, dimu::Phi,
7779 dimu::PhiAv, dimu::PhiCh,
7880 dimu::EnergyP, dimu::Pxp, dimu::Pyp, dimu::Pzp, dimu::Ptp, dimu::Etap, dimu::Phip,
@@ -112,6 +114,7 @@ DECLARE_SOA_TABLE(GenDimu, "AOD", "GENDIMU",
112114namespace recodimu
113115{
114116// dimuon
117+ DECLARE_SOA_COLUMN (RunNumber, runNumber, int );
115118DECLARE_SOA_COLUMN (M, m, float );
116119DECLARE_SOA_COLUMN (Pt, pt, float );
117120DECLARE_SOA_COLUMN (Rap, rap, float );
@@ -141,6 +144,7 @@ DECLARE_SOA_COLUMN(GenPhin, genPhin, float);
141144namespace o2 ::aod
142145{
143146DECLARE_SOA_TABLE (RecoDimu, " AOD" , " RECODIMU" ,
147+ recodimu::RunNumber,
144148 recodimu::M, recodimu::Pt, recodimu::Rap, recodimu::Phi,
145149 recodimu::PhiAv, recodimu::PhiCh,
146150 recodimu::Ptp, recodimu::Etap, recodimu::Phip,
@@ -371,8 +375,8 @@ struct FwdMuonsUPC {
371375 }
372376 }
373377
374- // template function that fills a map with the collision id of each udmccollision as key
375- // and a vector with the tracks
378+ // template function that fills a map with the collision id of each udcollision as key
379+ // and a vector with the tracks and corresponding geneated particles
376380 // map == (key, element) == (udCollisionId, vector(track1, mcPart1, track2, mcPart2))
377381 template <typename TTracks>
378382 void collectRecoCandID (std::unordered_map<int32_t , std::vector<int32_t >>& tracksPerCand, TTracks& tracks)
@@ -621,13 +625,15 @@ struct FwdMuonsUPC {
621625
622626 // store the event to save it into a tree
623627 if (tr1.sign () > 0 ) {
624- dimuSel (p.M (), p.E (), p.Px (), p.Py (), p.Pz (), p.Pt (), p.Rapidity (), p.Phi (),
628+ dimuSel (cand.runNumber (),
629+ p.M (), p.E (), p.Px (), p.Py (), p.Pz (), p.Pt (), p.Rapidity (), p.Phi (),
625630 phiAverage, phiCharge,
626631 p1.E (), p1.Px (), p1.Py (), p1.Pz (), p1.Pt (), p1.PseudoRapidity (), p1.Phi (),
627632 p2.E (), p2.Px (), p2.Py (), p2.Pz (), p2.Pt (), p2.PseudoRapidity (), p2.Phi (),
628633 zdc.timeA , zdc.enA , zdc.timeC , zdc.enC , znClass);
629634 } else {
630- dimuSel (p.M (), p.E (), p.Px (), p.Py (), p.Pz (), p.Pt (), p.Rapidity (), p.Phi (),
635+ dimuSel (cand.runNumber (),
636+ p.M (), p.E (), p.Px (), p.Py (), p.Pz (), p.Pt (), p.Rapidity (), p.Phi (),
631637 phiAverage, phiCharge,
632638 p2.E (), p2.Px (), p2.Py (), p2.Pz (), p2.Pt (), p2.PseudoRapidity (), p2.Phi (),
633639 p1.E (), p1.Px (), p1.Py (), p1.Pz (), p1.Pt (), p1.PseudoRapidity (), p1.Phi (),
@@ -836,7 +842,8 @@ struct FwdMuonsUPC {
836842
837843 // store the event to save it into a tree
838844 if (tr1.sign () > 0 ) {
839- dimuReco (p.M (), p.Pt (), p.Rapidity (), p.Phi (),
845+ dimuReco (cand.runNumber (),
846+ p.M (), p.Pt (), p.Rapidity (), p.Phi (),
840847 phiAverage, phiCharge,
841848 p1.Pt (), p1.PseudoRapidity (), p1.Phi (),
842849 p2.Pt (), p2.PseudoRapidity (), p2.Phi (),
@@ -845,7 +852,8 @@ struct FwdMuonsUPC {
845852 p1Mc.Pt (), p1Mc.PseudoRapidity (), p1Mc.Phi (),
846853 p2Mc.Pt (), p2Mc.PseudoRapidity (), p2Mc.Phi ());
847854 } else {
848- dimuReco (p.M (), p.Pt (), p.Rapidity (), p.Phi (),
855+ dimuReco (cand.runNumber (),
856+ p.M (), p.Pt (), p.Rapidity (), p.Phi (),
849857 phiAverage, phiCharge,
850858 p2.Pt (), p2.PseudoRapidity (), p2.Phi (),
851859 p1.Pt (), p1.PseudoRapidity (), p1.Phi (),
0 commit comments