Skip to content
Open
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
18 changes: 9 additions & 9 deletions PWGJE/Tasks/jetCorrelationD0.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -212,9 +212,9 @@ struct JetCorrelationD0 {
if (jetBase.pt() > pTHatMaxMCD * pTHat || pTHat < pTHatAbsoluteMin) { // reject events that are too hard / soft
return;
}
if (jetBase.has_matchedJetGeo()) { // geometric matching
if (jetBase.has_matchedJetGeo()) { // geometric matching
for (auto& jetTag : jetBase.template matchedJetGeo_as<std::decay_t<U>>()) {
if (jetTag.pt() > pTHatMaxMCP * pTHat) { // cuts overly hard jets from jettag (mcp) sample
if (jetTag.pt() > pTHatMaxMCP * pTHat) { // cuts overly hard jets from jettag (mcp) sample
continue;
}
registry.fill(HIST("hPtMatched"), jetBase.pt() - (rho * jetBase.area()), jetTag.pt(), weight);
Expand Down Expand Up @@ -278,7 +278,7 @@ struct JetCorrelationD0 {

for (const auto& d0Candidate : d0Candidates) {
const auto scores = d0Candidate.mlScores();
if(d0Candidate.pt() < d0PtCutMin) {
if (d0Candidate.pt() < d0PtCutMin) {
return;
}
fillD0Histograms(d0Candidate, scores);
Expand All @@ -291,11 +291,11 @@ struct JetCorrelationD0 {
d0Candidate.eta(),
d0Candidate.phi());
for (const auto& jet : jets) {
if(jet.pt() < jetPtCutMin) {
if (jet.pt() < jetPtCutMin) {
return;
}
float dphi = RecoDecay::constrainAngle(jet.phi() - d0Candidate.phi());
if (abs(dphi - M_PI) > (M_PI /2)) {
if (abs(dphi - M_PI) > (M_PI / 2)) {
return;
}
fillJetHistograms(jet, dphi);
Expand All @@ -318,7 +318,7 @@ struct JetCorrelationD0 {
tableCollision(collision.posZ());

for (const auto& d0MCPCandidate : d0MCPCandidates) {
if(d0MCPCandidate.pt() < d0PtCutMin) {
if (d0MCPCandidate.pt() < d0PtCutMin) {
return;
}
tableD0MCParticle(tableCollision.lastIndex(),
Expand All @@ -328,11 +328,11 @@ struct JetCorrelationD0 {
d0MCPCandidate.phi());

for (const auto& jet : jets) {
if(jet.pt() < jetPtCutMin) {
if (jet.pt() < jetPtCutMin) {
return;
}
float dphi = RecoDecay::constrainAngle(jet.phi() - d0MCPCandidate.phi());
if (abs(dphi - M_PI) > (M_PI /2)) {
if (abs(dphi - M_PI) > (M_PI / 2)) {
return;
}
fillJetHistograms(jet, dphi);
Expand Down Expand Up @@ -406,7 +406,7 @@ struct JetCorrelationD0 {
LOGF(info, "Collision ID %i, D0 pt %.2f, D0 eta %.2f, D0 phi %.2f, MCP origin %hhd, Reflection %i", CollIdx, particleMother.pt(), particleMother.eta(), particleMother.phi(), particleMother.originMcGen(), reflection);

// Jet matching
fillMatchedHistograms(mcdJets, mcpJets); // Do I need to include pthat cuts in loop rather than this function to actually do jet matching?
fillMatchedHistograms(mcdJets, mcpJets); // Do I need to include pthat cuts in loop rather than this function to actually do jet matching?

for (const auto& mcpJet : mcpJets) {
float dphi = RecoDecay::constrainAngle(mcpJet.phi() - d0MCDCandidate.phi());
Expand Down