@@ -420,7 +420,10 @@ struct DptDptCorrelationsTask {
420420 fhSum2DptDpt_vsDEtaDPhi[track1.trackacceptedid ()][track2.trackacceptedid ()]->AddBinContent (globalbin, dptdptw);
421421 fhSum2PtPt_vsDEtaDPhi[track1.trackacceptedid ()][track2.trackacceptedid ()]->AddBinContent (globalbin, track1.pt () * track2.pt () * corr);
422422 if constexpr (doinvmass) {
423- fhInvMass[track1.trackacceptedid ()][track2.trackacceptedid ()]->Fill (std::sqrt (getInvMassSquared (track1, poimass[track1.trackacceptedid ()], track2, poimass[track2.trackacceptedid ()])) * 1000 .0f );
423+ if ((track1.trackacceptedid () % 2 ) != (track2.trackacceptedid () % 2 )) {
424+ /* only opposite charge invariant mass*/
425+ fhInvMass[track1.trackacceptedid ()][track2.trackacceptedid ()]->Fill (std::sqrt (getInvMassSquared (track1, poimass[static_cast <int >(track1.trackacceptedid () / 2 )], track2, poimass[static_cast <int >(track2.trackacceptedid () / 2 )])) * 1000 .0f );
426+ }
424427 }
425428 }
426429 fhN2_vsPtPt[track1.trackacceptedid ()][track2.trackacceptedid ()]->Fill (track1.pt (), track2.pt (), corr);
@@ -865,7 +868,7 @@ struct DptDptCorrelationsTask {
865868 tnames.push_back (std::string (TString::Format (" %sP" , pidselector.getSpeciesFName (ix)).Data ()));
866869 tnames.push_back (std::string (TString::Format (" %sM" , pidselector.getSpeciesFName (ix)).Data ()));
867870 poimass.push_back (pidselector.getSpeciesMass (ix));
868- LOGF (info, " Incorporated species name %s to the analysis" , poinames[ix].c_str ());
871+ LOGF (info, " Incorporated species name %s with mass %f to the analysis" , poinames[ix].c_str (), poimass[ix] );
869872 }
870873 }
871874 uint ntracknames = tnames.size ();
0 commit comments