@@ -84,6 +84,10 @@ struct phiflowder {
8484 histos.add (" hpQypvscent" , " hpQypvscent" , HistType::kTHnSparseF , {axisCent, axisV1}, true );
8585 histos.add (" hpQytvscent" , " hpQytvscent" , HistType::kTHnSparseF , {axisCent, axisV1}, true );
8686 histos.add (" hMixpairs" , " hMixpairs" , HistType::kTHnSparseF , {axisNPairs}, true );
87+ histos.add (" hSparseSameUx" , " Same-event #Sigma u_{x};M;pT;centrality;eta" , kTHnSparseF , {axisInvMass, axisPhiPt, axisCent, axisEta, axisV1}, true );
88+ histos.add (" hSparseSameUy" , " Same-event #Sigma u_{y};M;pT;centrality;eta" , kTHnSparseF , {axisInvMass, axisPhiPt, axisCent, axisEta, axisV1}, true );
89+ histos.add (" hSparseSameUQA" , " Same-event #Sigma u#upoint Q_{A};M;pT;centrality;eta" , kTHnSparseF , {axisInvMass, axisPhiPt, axisCent, axisEta, axisV1}, true );
90+ histos.add (" hSparseSameUQC" , " Same-event #Sigma u#upoint Q_{C};M;pT;centrality;eta" , kTHnSparseF , {axisInvMass, axisPhiPt, axisCent, axisEta, axisV1}, true );
8791 }
8892
8993 uint8_t getRequiredPidBit () const
@@ -291,12 +295,18 @@ struct phiflowder {
291295 auto uy = sinNPhi; // imaginary part of candidate q vector
292296 auto oddv1 = ux * (qxZDCA - qxZDCC) + uy * (qyZDCA - qyZDCC);
293297 // auto evenv1 = ux * (qxZDCA + qxZDCC) + uy * (qyZDCA + qyZDCC);
298+ auto uQA = ux * qxZDCA + uy * qyZDCA;
299+ auto uQC = ux * qxZDCC + uy * qyZDCC;
294300
295301 if (phiMass < massMin || phiMass > massMax) {
296302 continue ;
297303 }
298304
299305 histos.fill (HIST (" hSparseSame" ), phiMass, phiPt, centrality, etaCand, oddv1);
306+ histos.fill (HIST (" hSparseSameUx" ), phiMass, phiPt, centrality, etaCand, ux);
307+ histos.fill (HIST (" hSparseSameUy" ), phiMass, phiPt, centrality, etaCand, uy);
308+ histos.fill (HIST (" hSparseSameUQA" ), phiMass, phiPt, centrality, etaCand, uQA);
309+ histos.fill (HIST (" hSparseSameUQC" ), phiMass, phiPt, centrality, etaCand, uQC);
300310 // nPhiSame++;
301311 }
302312 }
0 commit comments