Skip to content
Merged
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
36 changes: 18 additions & 18 deletions PWGCF/TwoParticleCorrelations/TableProducer/identifiedBfFilter.cxx
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.

Check warning on line 1 in PWGCF/TwoParticleCorrelations/TableProducer/identifiedBfFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[doc/file]

Provide mandatory file documentation.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
Expand All @@ -8,7 +8,7 @@
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
#include "PWGCF/TwoParticleCorrelations/TableProducer/identifiedBfFilter.h"

Check warning on line 11 in PWGCF/TwoParticleCorrelations/TableProducer/identifiedBfFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[doc/file]

Documentation for \file is missing, incorrect or misplaced.

Check warning on line 11 in PWGCF/TwoParticleCorrelations/TableProducer/identifiedBfFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[doc/file]

Documentation for \brief is missing, incorrect or misplaced.

Check warning on line 11 in PWGCF/TwoParticleCorrelations/TableProducer/identifiedBfFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[doc/file]

Documentation for \author is missing, incorrect or misplaced.

#include <cmath>
#include <algorithm>
Expand Down Expand Up @@ -831,9 +831,9 @@
48, -6, 6,
ptbins, ptlow, ptup);
fhNSigmaTPCIdTrks[sp] = new TH2F(TString::Format("fhNSigmaTPC_IdTrks_%s", speciesName[sp]).Data(),
TString::Format("N Sigma from TPC vs P for Identified %s;N #sigma;p (GeV/c)", speciesTitle[sp]).Data(),
48, -6, 6,
ptbins, ptlow, ptup);
TString::Format("N Sigma from TPC vs P for Identified %s;N #sigma;p (GeV/c)", speciesTitle[sp]).Data(),
48, -6, 6,
ptbins, ptlow, ptup);
}
LOGF(info, "Making histos");

Expand Down Expand Up @@ -1379,15 +1379,15 @@
MatchRecoGenSpecies spMinNSigma = kWrongSpecies;
for (int sp = 0; sp < kIdBfNoOfSpecies; ++sp) {
if (std::fabs(nsigmas[sp]) < std::fabs(minNSigma)) { // Check if species nsigma is less than current nsigma
minNSigma = nsigmas[sp]; // If yes, set species nsigma to current nsigma
spMinNSigma = MatchRecoGenSpecies(sp); // set current species sp number to current sp
minNSigma = nsigmas[sp]; // If yes, set species nsigma to current nsigma
spMinNSigma = MatchRecoGenSpecies(sp); // set current species sp number to current sp
}
}
bool doublematch = false;
MatchRecoGenSpecies spDouble = kWrongSpecies;
if (minNSigma < maxPIDSigma && minNSigma > minPIDSigma) { // Check that current nsigma is in accpetance range
if (minNSigma < maxPIDSigma && minNSigma > minPIDSigma) { // Check that current nsigma is in accpetance range
for (int sp = 0; (sp < kIdBfNoOfSpecies) && !doublematch; ++sp) { // iterate over all species while there's no double match and we're in the list
if (sp != spMinNSigma) { // for species not current minimum nsigma species
if (sp != spMinNSigma) { // for species not current minimum nsigma species
if (nsigmas[sp] < maxRejectSigma && nsigmas[sp] > minRejectSigma) { // If secondary species is in rejection range
doublematch = true; // Set double match true
spDouble = MatchRecoGenSpecies(sp);
Expand Down Expand Up @@ -1494,7 +1494,7 @@
if (pdgpart != nullptr) {
charge = getCharge(pdgpart->Charge());
}
if ((particle.flags()&0x8)!=0x8){
if ((particle.flags() & 0x8) != 0x8) {
if (particle.isPhysicalPrimary() && std::fabs(charge) > 0.0) {
if ((particle.mcCollisionId() == 0) && traceCollId0) {
LOGF(info, "Particle %d passed isPhysicalPrimary", particle.globalIndex());
Expand Down Expand Up @@ -1534,10 +1534,10 @@
LOGF(info, "Particle %d NOT passed isPhysicalPrimary", particle.globalIndex());
}
}
}else{
if ((particle.mcCollisionId() == 0) && traceCollId0) {
LOGF(info, "Particle %d Out of Bunch Pileup", particle.globalIndex());
}
} else {
if ((particle.mcCollisionId() == 0) && traceCollId0) {
LOGF(info, "Particle %d Out of Bunch Pileup", particle.globalIndex());
}
}
return kWrongSpecies;
}
Expand Down Expand Up @@ -1675,15 +1675,15 @@
}

float dcaxy = std::sqrt((particle.vx() - collision.posX()) * (particle.vx() - collision.posX()) +
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY()));
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY()));
if (traceDCAOutliers.mDoIt && (traceDCAOutliers.mLowValue < dcaxy) && (dcaxy < traceDCAOutliers.mUpValue)) {
fhTrueDCAxyBid->Fill(TString::Format("%d", particle.pdgCode()).Data(), 1.0);
}

fhTrueDCAxyB->Fill(std::sqrt((particle.vx() - collision.posX()) * (particle.vx() - collision.posX()) +
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY())));
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY())));
fhTrueDCAxyzB->Fill(std::sqrt((particle.vx() - collision.posX()) * (particle.vx() - collision.posX()) +
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY())),
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY())),
(particle.vz() - collision.posZ()));
fhTrueDCAzB->Fill((particle.vz() - collision.posZ()));
}
Expand All @@ -1698,18 +1698,18 @@
fhTrueEtaA->Fill(particle.eta());
fhTruePhiA->Fill(particle.phi());
float dcaxy = std::sqrt((particle.vx() - collision.posX()) * (particle.vx() - collision.posX()) +
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY()));
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY()));
if (traceDCAOutliers.mDoIt && (traceDCAOutliers.mLowValue < dcaxy) && (dcaxy < traceDCAOutliers.mUpValue)) {
LOGF(info, "DCAxy outlier: Particle with index %d and pdg code %d assigned to MC collision %d, pT: %f, phi: %f, eta: %f",
particle.globalIndex(), particle.pdgCode(), particle.mcCollisionId(), particle.pt(), particle.phi(), particle.eta());
LOGF(info, " With status %d and flags %0X", particle.statusCode(), particle.flags());
}

fhTrueDCAxyA->Fill(std::sqrt((particle.vx() - collision.posX()) * (particle.vx() - collision.posX()) +
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY())));
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY())));
fhTrueDCAzA->Fill((particle.vz() - collision.posZ()));
fhTrueDCAxyzA->Fill(std::sqrt((particle.vx() - collision.posX()) * (particle.vx() - collision.posX()) +
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY())),
(particle.vy() - collision.posY()) * (particle.vy() - collision.posY())),
(particle.vz() - collision.posZ()));
}
fhTruePA[sp]->Fill(particle.p());
Expand Down
Loading