@@ -322,12 +322,12 @@ struct lmeelfcocktail {
322322 addHistogram1D (" OpAng" , opAng_axis, i);
323323 addHistogram1D (" Mee" , mass_axis, i);
324324 addHistogram1D (" Ptee" , ptee_axis, i);
325- // addHistogram1D_stage("Dca", dca_axis, i, "rec/");
326- // addHistogram1D_stage("Dcaee", dcaee_axis, i, "rec/");
325+ // addHistogram1D_stage("Dca", dca_axis, i, "rec/");
326+ // addHistogram1D_stage("Dcaee", dcaee_axis, i, "rec/");
327327 i = -1 ;
328- // addHistogram2D_stage("DcaVsPt", dca_axis, pt_axis, i, "rec/");
329- // addHistogram2D_stage("DcaeeVsPtee", dcaee_axis, ptee_axis, i, "rec/");
330- // addHistogram2D_stage("DcaeeVsMee", dcaee_axis, mass_axis, i, "rec/");
328+ // addHistogram2D_stage("DcaVsPt", dca_axis, pt_axis, i, "rec/");
329+ // addHistogram2D_stage("DcaeeVsPtee", dcaee_axis, ptee_axis, i, "rec/");
330+ // addHistogram2D_stage("DcaeeVsMee", dcaee_axis, mass_axis, i, "rec/");
331331 i = -1 ;
332332 addHistogramND (" MeeVsPteeVsCos2DPhiRP" , std::vector<AxisSpec>{mass_axis, ptee_axis, cos2dphi_axis}, i);
333333 }
@@ -354,7 +354,7 @@ struct lmeelfcocktail {
354354 int nPos = 0 ;
355355
356356 ROOT::Math::PtEtaPhiMVector pEleGen, pPosGen, pEleRec, pPosRec;
357- float weight (1 .), effEle (1 .), effPos (1 .);// , dcaEle(0.), dcaPos(0.);
357+ float weight (1 .), effEle (1 .), effPos (1 .); // , dcaEle(0.), dcaPos(0.);
358358 for (const auto & daughter : particle.daughters_as <McParticlesSmeared>()) {
359359 int temp_pdg = daughter.pdgCode ();
360360 if (temp_pdg == 11 ) {
@@ -438,12 +438,12 @@ struct lmeelfcocktail {
438438
439439 if (s == kRec ) { // dca only at rec. level
440440 if (acceptedEle) {
441- // fillHistogram1D("Dca", s, pdg, other_daughter_pdg, dcaEle, weightEle);
442- // fillHistogram2D("DcaVsPt", s, pdg, other_daughter_pdg, dcaEle, pEle.Pt(), weightEle);
441+ // fillHistogram1D("Dca", s, pdg, other_daughter_pdg, dcaEle, weightEle);
442+ // fillHistogram2D("DcaVsPt", s, pdg, other_daughter_pdg, dcaEle, pEle.Pt(), weightEle);
443443 }
444444 if (acceptedPos) {
445- // fillHistogram1D("Dca", s, pdg, other_daughter_pdg, dcaPos, weightPos);
446- // fillHistogram2D("DcaVsPt", s, pdg, other_daughter_pdg, dcaPos, pPos.Pt(), weightPos);
445+ // fillHistogram1D("Dca", s, pdg, other_daughter_pdg, dcaPos, weightPos);
446+ // fillHistogram2D("DcaVsPt", s, pdg, other_daughter_pdg, dcaPos, pPos.Pt(), weightPos);
447447 }
448448 }
449449
@@ -461,7 +461,7 @@ struct lmeelfcocktail {
461461 float ptee = p12.Pt ();
462462 float opAng = o2::aod::pwgem::dilepton::utils::pairutil::getOpeningAngle (pPos.Px (), pPos.Py (), pPos.Pz (), pEle.Px (), pEle.Py (), pEle.Pz ());
463463 float phiV = o2::aod::pwgem::dilepton::utils::pairutil::getPhivPair (pPos.Px (), pPos.Py (), pPos.Pz (), pEle.Px (), pEle.Py (), pEle.Pz (), 1 , -1 , 1 );
464- // float dcaee = sqrt((pow(dcaEle, 2) + pow(dcaPos, 2)) / 2);
464+ // float dcaee = sqrt((pow(dcaEle, 2) + pow(dcaPos, 2)) / 2);
465465 float cos2dphi = std::cos (2 .f * p12.Phi ()); // PsiRP = 0 rad.
466466 double values[3 ] = {mee, ptee, cos2dphi};
467467 fillHistogramND (" MeeVsPteeVsCos2DPhiRP" , s, pdg, other_daughter_pdg, values, pairWeight);
@@ -470,9 +470,9 @@ struct lmeelfcocktail {
470470 fillHistogram1D (" PhiV" , s, pdg, other_daughter_pdg, phiV, pairWeight);
471471 fillHistogram1D (" OpAng" , s, pdg, other_daughter_pdg, opAng, pairWeight);
472472 if (s == kRec ) { // dca only at rec. level
473- // fillHistogram1D("Dcaee", s, pdg, other_daughter_pdg, dcaee, pairWeight);
474- // fillHistogram2D("DcaeeVsPtee", s, pdg, other_daughter_pdg, dcaee, ptee, pairWeight);
475- // fillHistogram2D("DcaeeVsMee", s, pdg, other_daughter_pdg, dcaee, mee, pairWeight);
473+ // fillHistogram1D("Dcaee", s, pdg, other_daughter_pdg, dcaee, pairWeight);
474+ // fillHistogram2D("DcaeeVsPtee", s, pdg, other_daughter_pdg, dcaee, ptee, pairWeight);
475+ // fillHistogram2D("DcaeeVsMee", s, pdg, other_daughter_pdg, dcaee, mee, pairWeight);
476476 }
477477 }
478478 }
0 commit comments