@@ -236,15 +236,15 @@ class FemtoDreamCollisionSelection
236236 }
237237
238238 mHistogramQn = registry;
239- mHistogramQn ->add (" Event/hN2allQn" , " ;centrality; #sum Re(Q_{2,A} Q_{2,B}^{*})" , kTH1F , {{centBins, 0 , 100 }});
240- mHistogramQn ->add (" Event/hD2allQn" , " ;centrality; #sum (W_{A} W_{B})" , kTH1F , {{centBins, 0 , 100 }});
239+ mHistogramQn ->add (" Event/hN2allQn" , " ;centrality; #sum Re(Q_{2,A} Q_{2,B}^{*})" , kTH1F , {{centBins, 0 , 100 }});
240+ mHistogramQn ->add (" Event/hD2allQn" , " ;centrality; #sum (W_{A} W_{B})" , kTH1F , {{centBins, 0 , 100 }});
241241 mHistogramQn ->get <TH1>(HIST (" Event/hN2allQn" ))->Sumw2 ();
242242 mHistogramQn ->get <TH1>(HIST (" Event/hD2allQn" ))->Sumw2 ();
243243
244244 if (doQnSeparation) {
245245 for (int iqn (0 ); iqn < mumQnBins; ++iqn) {
246246 hN2.push_back (mHistogramQn ->add ((" Qn/hN2_" + std::to_string (iqn)).c_str (), " ;centrality; #sum Re(Q_{2,A} Q_{2,B}^{*})" , kTH1F , {{centBins, 0 , 100 }}));
247- hD2.push_back (mHistogramQn ->add ((" Qn/hD2_" + std::to_string (iqn)).c_str (), " ;centrality; #sum (W_{A} W_{B})" , kTH1F , {{centBins, 0 , 100 }}));
247+ hD2.push_back (mHistogramQn ->add ((" Qn/hD2_" + std::to_string (iqn)).c_str (), " ;centrality; #sum (W_{A} W_{B})" , kTH1F , {{centBins, 0 , 100 }}));
248248 }
249249 for (int iqn (0 ); iqn < mumQnBins; ++iqn) {
250250 std::get<std::shared_ptr<TH1>>(hN2[iqn])->Sumw2 ();
@@ -441,29 +441,29 @@ class FemtoDreamCollisionSelection
441441 int nA = 0 , nB = 0 ;
442442
443443 for (auto const & trk : tracks) {
444- const double pt = trk.pt ();
444+ const double pt = trk.pt ();
445445 const double eta = trk.eta ();
446446 if (pt < ptMin || pt > ptMax) {
447447 continue ;
448448 }
449-
450- const double w = 1.0 ; // TODO: NUA/NUE weight if you want
449+
450+ const double w = 1.0 ; // TODO: NUA/NUE weight if you want
451451 const double phi = trk.phi ();
452- const double c = w * TMath::Cos (harmonic * phi);
453- const double s = w * TMath::Sin (harmonic * phi);
452+ const double c = w * TMath::Cos (harmonic * phi);
453+ const double s = w * TMath::Sin (harmonic * phi);
454454
455455 if (eta > fEtaGap ) {
456- Q2A_re += c;
457- Q2A_im += s;
458- WA += w;
456+ Q2A_re += c;
457+ Q2A_im += s;
458+ WA += w;
459459 nA++;
460- } else if (eta < -1 * fEtaGap ) {
461- Q2B_re += c;
462- Q2B_im += s;
463- WB += w;
460+ } else if (eta < -1 * fEtaGap ) {
461+ Q2B_re += c;
462+ Q2B_im += s;
463+ WB += w;
464464 nB++;
465465 }
466- }
466+ }
467467
468468 // need at least 1 track on each side to form pairs; for stability, require >=2
469469 if (nA < 2 || nB < 2 ) {
0 commit comments