Skip to content

Commit d87a34b

Browse files
feisenhufeisenhu
committed
checkMCPairTemplate.cxx: Add b2c2l_b2c2l_Hc_ctau histograms which select decays based on life time, fix clang
1 parent ae0edcd commit d87a34b

File tree

1 file changed

+112
-0
lines changed

1 file changed

+112
-0
lines changed

‎PWGEM/Dilepton/Tasks/checkMCPairTemplate.cxx‎

Lines changed: 112 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -497,6 +497,12 @@ struct checkMCPairTemplate {
497497
for (int ib = 0; ib < nb_b; ib++) {
498498
fRegistry.addClone("Generated/ccbar/c2l_c2l/", Form("Generated/bbbar/b2%s2l_b2%s2l/", sum_charmed_baryons[ib].data(), sum_charmed_baryons[ib].data()));
499499
}
500+
fRegistry.addClone("Generated/ccbar/c2l_c2l/", "Generated/bbbar/b2c2l_b2c2l_Hc_ctau50/");
501+
fRegistry.addClone("Generated/ccbar/c2l_c2l/", "Generated/bbbar/b2c2l_b2c2l_Hc_ctau130/");
502+
fRegistry.addClone("Generated/ccbar/c2l_c2l/", "Generated/bbbar/b2c2l_b2c2l_Hc_ctau300/");
503+
fRegistry.addClone("Generated/ccbar/c2l_c2l/", "Generated/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau130/");
504+
fRegistry.addClone("Generated/ccbar/c2l_c2l/", "Generated/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau300/");
505+
fRegistry.addClone("Generated/ccbar/c2l_c2l/", "Generated/bbbar/b2c2l_b2c2l_Hc_ctau130_ctau300/");
500506
}
501507

502508
// evaluate acceptance for polarization
@@ -611,6 +617,12 @@ struct checkMCPairTemplate {
611617
for (int ib = 0; ib < nb_b; ib++) {
612618
fRegistry.addClone("Pair/ccbar/c2l_c2l/", Form("Pair/bbbar/b2%s2l_b2%s2l/", sum_charmed_baryons[ib].data(), sum_charmed_baryons[ib].data()));
613619
}
620+
fRegistry.addClone("Pair/ccbar/c2l_c2l/", "Pair/bbbar/b2c2l_b2c2l_Hc_ctau50/");
621+
fRegistry.addClone("Pair/ccbar/c2l_c2l/", "Pair/bbbar/b2c2l_b2c2l_Hc_ctau130/");
622+
fRegistry.addClone("Pair/ccbar/c2l_c2l/", "Pair/bbbar/b2c2l_b2c2l_Hc_ctau300/");
623+
fRegistry.addClone("Pair/ccbar/c2l_c2l/", "Pair/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau130/");
624+
fRegistry.addClone("Pair/ccbar/c2l_c2l/", "Pair/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau300/");
625+
fRegistry.addClone("Pair/ccbar/c2l_c2l/", "Pair/bbbar/b2c2l_b2c2l_Hc_ctau130_ctau300/");
614626
}
615627

616628
// for correlated bkg due to mis-identified hadrons, and true combinatorial bkg
@@ -1328,6 +1340,56 @@ struct checkMCPairTemplate {
13281340
fRegistry.fill(HIST("Generated/bbbar/b2Lcpm2l_b2Lcpm2l/lsmm/hs"), mass, pt, weight);
13291341
}
13301342
}
1343+
1344+
if ((std::abs(pdgMotherC1) == kLambdaCPlus || std::abs(pdgMotherC1) == kXiC0) && (std::abs(pdgMotherC2) == kLambdaCPlus || std::abs(pdgMotherC2) == kXiC0)) {
1345+
if (sign1 * sign2 < 0) { // ULS
1346+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau50/uls/hs"), mass, pt, weight);
1347+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1348+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau50/lspp/hs"), mass, pt, weight);
1349+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1350+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau50/lsmm/hs"), mass, pt, weight);
1351+
}
1352+
} else if ((std::abs(pdgMotherC1) == kD0 || std::abs(pdgMotherC1) == kDS || std::abs(pdgMotherC1) == kXiCPlus) && (std::abs(pdgMotherC2) == kD0 || std::abs(pdgMotherC2) == kDS || std::abs(pdgMotherC2) == kXiCPlus)) {
1353+
if (sign1 * sign2 < 0) { // ULS
1354+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau130/uls/hs"), mass, pt, weight);
1355+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1356+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau130/lspp/hs"), mass, pt, weight);
1357+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1358+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau130/lsmm/hs"), mass, pt, weight);
1359+
}
1360+
} else if ((std::abs(pdgMotherC1) == kDPlus) && (std::abs(pdgMotherC2) == kDPlus)) {
1361+
if (sign1 * sign2 < 0) { // ULS
1362+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau300/uls/hs"), mass, pt, weight);
1363+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1364+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau300/lspp/hs"), mass, pt, weight);
1365+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1366+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau300/lsmm/hs"), mass, pt, weight);
1367+
}
1368+
} else if (((std::abs(pdgMotherC1) == kLambdaCPlus || std::abs(pdgMotherC1) == kXiC0) && (std::abs(pdgMotherC2) == kD0 || std::abs(pdgMotherC2) == kDS || std::abs(pdgMotherC2) == kXiCPlus)) || ((std::abs(pdgMotherC2) == kLambdaCPlus || std::abs(pdgMotherC2) == kXiC0) && (std::abs(pdgMotherC1) == kD0 || std::abs(pdgMotherC1) == kDS || std::abs(pdgMotherC1) == kXiCPlus))) {
1369+
if (sign1 * sign2 < 0) { // ULS
1370+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau130/uls/hs"), mass, pt, weight);
1371+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1372+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau130/lspp/hs"), mass, pt, weight);
1373+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1374+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau130/lsmm/hs"), mass, pt, weight);
1375+
}
1376+
} else if (((std::abs(pdgMotherC1) == kLambdaCPlus || std::abs(pdgMotherC1) == kXiC0) && (std::abs(pdgMotherC2) == kDPlus)) || ((std::abs(pdgMotherC2) == kLambdaCPlus || std::abs(pdgMotherC2) == kXiC0) && (std::abs(pdgMotherC1) == kDPlus))) {
1377+
if (sign1 * sign2 < 0) { // ULS
1378+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau300/uls/hs"), mass, pt, weight);
1379+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1380+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau300/lspp/hs"), mass, pt, weight);
1381+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1382+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau300/lsmm/hs"), mass, pt, weight);
1383+
}
1384+
} else if (((std::abs(pdgMotherC1) == kD0 || std::abs(pdgMotherC1) == kDS || std::abs(pdgMotherC1) == kXiCPlus) && (std::abs(pdgMotherC2) == kDPlus)) || ((std::abs(pdgMotherC2) == kD0 || std::abs(pdgMotherC2) == kDS || std::abs(pdgMotherC2) == kXiCPlus) && (std::abs(pdgMotherC1) == kDPlus))) {
1385+
if (sign1 * sign2 < 0) { // ULS
1386+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau130_ctau300/uls/hs"), mass, pt, weight);
1387+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1388+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau130_ctau300/lspp/hs"), mass, pt, weight);
1389+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1390+
fRegistry.fill(HIST("Generated/bbbar/b2c2l_b2c2l_Hc_ctau130_ctau300/lsmm/hs"), mass, pt, weight);
1391+
}
1392+
}
13311393
}
13321394
}
13331395

@@ -1677,6 +1739,56 @@ struct checkMCPairTemplate {
16771739
fRegistry.fill(HIST("Pair/bbbar/b2Lcpm2l_b2Lcpm2l/lsmm/hs"), mass, pt, pair_dca, weight);
16781740
}
16791741
}
1742+
1743+
if ((std::abs(pdgMotherC1) == kLambdaCPlus || std::abs(pdgMotherC1) == kXiC0) && (std::abs(pdgMotherC2) == kLambdaCPlus || std::abs(pdgMotherC2) == kXiC0)) {
1744+
if (sign1 * sign2 < 0) { // ULS
1745+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau50/uls/hs"), mass, pt, pair_dca, weight);
1746+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1747+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau50/lspp/hs"), mass, pt, pair_dca, weight);
1748+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1749+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau50/lsmm/hs"), mass, pt, pair_dca, weight);
1750+
}
1751+
} else if ((std::abs(pdgMotherC1) == kD0 || std::abs(pdgMotherC1) == kDS || std::abs(pdgMotherC1) == kXiCPlus) && (std::abs(pdgMotherC2) == kD0 || std::abs(pdgMotherC2) == kDS || std::abs(pdgMotherC2) == kXiCPlus)) {
1752+
if (sign1 * sign2 < 0) { // ULS
1753+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau130/uls/hs"), mass, pt, pair_dca, weight);
1754+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1755+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau130/lspp/hs"), mass, pt, pair_dca, weight);
1756+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1757+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau130/lsmm/hs"), mass, pt, pair_dca, weight);
1758+
}
1759+
} else if (std::abs(pdgMotherC1) == kDPlus && std::abs(pdgMotherC2) == kDPlus) {
1760+
if (sign1 * sign2 < 0) { // ULS
1761+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau300/uls/hs"), mass, pt, pair_dca, weight);
1762+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1763+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau300/lspp/hs"), mass, pt, pair_dca, weight);
1764+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1765+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau300/lsmm/hs"), mass, pt, pair_dca, weight);
1766+
}
1767+
} else if (((std::abs(pdgMotherC1) == kLambdaCPlus || std::abs(pdgMotherC1) == kXiC0) && (std::abs(pdgMotherC2) == kD0 || std::abs(pdgMotherC2) == kDS || std::abs(pdgMotherC2) == kXiCPlus)) || ((std::abs(pdgMotherC2) == kLambdaCPlus || std::abs(pdgMotherC2) == kXiC0) && (std::abs(pdgMotherC1) == kD0 || std::abs(pdgMotherC1) == kDS || std::abs(pdgMotherC1) == kXiCPlus))) {
1768+
if (sign1 * sign2 < 0) { // ULS
1769+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau130/uls/hs"), mass, pt, pair_dca, weight);
1770+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1771+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau130/lspp/hs"), mass, pt, pair_dca, weight);
1772+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1773+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau130/lsmm/hs"), mass, pt, pair_dca, weight);
1774+
}
1775+
} else if (((std::abs(pdgMotherC1) == kLambdaCPlus || std::abs(pdgMotherC1) == kXiC0) && (std::abs(pdgMotherC2) == kDPlus)) || ((std::abs(pdgMotherC2) == kLambdaCPlus || std::abs(pdgMotherC2) == kXiC0) && (std::abs(pdgMotherC1) == kDPlus))) {
1776+
if (sign1 * sign2 < 0) { // ULS
1777+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau300/uls/hs"), mass, pt, pair_dca, weight);
1778+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1779+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau300/lspp/hs"), mass, pt, pair_dca, weight);
1780+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1781+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau50_ctau300/lsmm/hs"), mass, pt, pair_dca, weight);
1782+
}
1783+
} else if (((std::abs(pdgMotherC1) == kD0 || std::abs(pdgMotherC1) == kDS || std::abs(pdgMotherC1) == kXiCPlus) && (std::abs(pdgMotherC2) == kDPlus)) || ((std::abs(pdgMotherC2) == kD0 || std::abs(pdgMotherC2) == kDS || std::abs(pdgMotherC2) == kXiCPlus) && (std::abs(pdgMotherC1) == kDPlus))) {
1784+
if (sign1 * sign2 < 0) { // ULS
1785+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau130_ctau300/uls/hs"), mass, pt, pair_dca, weight);
1786+
} else if (sign1 > 0 && sign2 > 0) { // LS++
1787+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau130_ctau300/lspp/hs"), mass, pt, pair_dca, weight);
1788+
} else if (sign1 < 0 && sign2 < 0) { // LS--
1789+
fRegistry.fill(HIST("Pair/bbbar/b2c2l_b2c2l_Hc_ctau130_ctau300/lsmm/hs"), mass, pt, pair_dca, weight);
1790+
}
1791+
}
16801792
}
16811793
}
16821794

0 commit comments

Comments
 (0)