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
18 changes: 9 additions & 9 deletions PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -336,7 +336,7 @@

template <class T>
bool singleSpeciesTPCNSigma(T const& track, int species) // make cut configurable
{ // reject any track that has nsigma < 3 for more than 1 species
{ // reject any track that has nsigma < 3 for more than 1 species
if (track.tpcNSigmaStoreEl() < nsigmaRejection || track.tpcNSigmaStoreMu() < nsigmaRejection || track.tpcNSigmaStorePi() < nsigmaRejection || track.tpcNSigmaStoreKa() < nsigmaRejection || track.tpcNSigmaStoreTr() < nsigmaRejection || track.tpcNSigmaStoreAl() < nsigmaRejection)
return false;
switch (species) {
Expand All @@ -357,9 +357,9 @@

if (tightCuts) { // for correlation function
// DCA
if (TMath::Abs(track.dcaXY()) > protonDCAxyCF)

Check warning on line 360 in PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root-entity]

Consider replacing ROOT entities with equivalents from standard C++ or from O2.
return false;
if (TMath::Abs(track.dcaZ()) > protonDCAzCF)

Check warning on line 362 in PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root-entity]

Consider replacing ROOT entities with equivalents from standard C++ or from O2.
return false;

registryData.fill(HIST("hTPCnsigmaProtonCF"), track.pt(), track.tpcNSigmaPr());
Expand All @@ -368,9 +368,9 @@

// nsigma
double tofNsigma = track.hasTOF() ? track.tofNSigmaPr() : 999;
if ((track.pt() < protonTPCTOFpT && (TMath::Abs(track.tpcNSigmaPr()) > protonNsigma)) || (track.pt() > protonTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaPr()*track.tpcNSigmaPr() + tofNsigma*tofNsigma) > protonNsigma)))
if (TMath::Sqrt(track.tpcNSigmaPr()*track.tpcNSigmaPr() + tofNsigma*tofNsigma) > protonNsigma)
return false;
if ((track.pt() < protonTPCTOFpT && (TMath::Abs(track.tpcNSigmaPr()) > protonNsigma)) || (track.pt() > protonTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaPr() * track.tpcNSigmaPr() + tofNsigma * tofNsigma) > protonNsigma)))

Check warning on line 371 in PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root-entity]

Consider replacing ROOT entities with equivalents from standard C++ or from O2.
if (TMath::Sqrt(track.tpcNSigmaPr() * track.tpcNSigmaPr() + tofNsigma * tofNsigma) > protonNsigma)

Check warning on line 372 in PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root-entity]

Consider replacing ROOT entities with equivalents from standard C++ or from O2.
return false;
if (!singleSpeciesTPCNSigma(track, 1))
return false;
} else { // for yields
Expand Down Expand Up @@ -418,7 +418,7 @@

// nsigma
double tofNsigma = track.hasTOF() ? track.tofNSigmaPr() : 999;
if ((track.pt() < antiprotonTPCTOFpT && (TMath::Abs(track.tpcNSigmaPr()) > antiprotonNsigma)) || (track.pt() > antiprotonTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaPr()*track.tpcNSigmaPr() + tofNsigma*tofNsigma) > antiprotonNsigma)))
if ((track.pt() < antiprotonTPCTOFpT && (TMath::Abs(track.tpcNSigmaPr()) > antiprotonNsigma)) || (track.pt() > antiprotonTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaPr() * track.tpcNSigmaPr() + tofNsigma * tofNsigma) > antiprotonNsigma)))
return false;
if (!singleSpeciesTPCNSigma(track, 1))
return false;
Expand Down Expand Up @@ -467,7 +467,7 @@

// nsigma
double tofNsigma = track.hasTOF() ? track.tofNSigmaDe() : 999;
if ((track.pt() < nucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaDe()) > nucleiNsigma)) || (track.pt() > nucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaDe()*track.tpcNSigmaDe() + tofNsigma*tofNsigma) > nucleiNsigma)))
if ((track.pt() < nucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaDe()) > nucleiNsigma)) || (track.pt() > nucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaDe() * track.tpcNSigmaDe() + tofNsigma * tofNsigma) > nucleiNsigma)))
return false;
if (!singleSpeciesTPCNSigma(track, 2))
return false;
Expand Down Expand Up @@ -507,7 +507,7 @@

// nsigma
double tofNsigma = track.hasTOF() ? track.tofNSigmaHe() : 999;
if ((track.pt() < nucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaHe()) > nucleiNsigma)) || (track.pt() > nucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaHe()*track.tpcNSigmaHe() + tofNsigma*tofNsigma) > nucleiNsigma)))
if ((track.pt() < nucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaHe()) > nucleiNsigma)) || (track.pt() > nucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaHe() * track.tpcNSigmaHe() + tofNsigma * tofNsigma) > nucleiNsigma)))
return false;
if (!singleSpeciesTPCNSigma(track, 3))
return false;
Expand Down Expand Up @@ -558,7 +558,7 @@

// nsigma
double tofNsigma = track.hasTOF() ? track.tofNSigmaDe() : 999;
if ((track.pt() < antinucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaDe()) > antinucleiNsigma)) || (track.pt() > antinucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaDe()*track.tpcNSigmaDe() + tofNsigma*tofNsigma) > antinucleiNsigma)))
if ((track.pt() < antinucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaDe()) > antinucleiNsigma)) || (track.pt() > antinucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaDe() * track.tpcNSigmaDe() + tofNsigma * tofNsigma) > antinucleiNsigma)))
return false;
if (!singleSpeciesTPCNSigma(track, 2))
return false;
Expand Down Expand Up @@ -598,7 +598,7 @@

// nsigma
double tofNsigma = track.hasTOF() ? track.tofNSigmaHe() : 999;
if ((track.pt() < antinucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaHe()) > antinucleiNsigma)) || (track.pt() > antinucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaHe()*track.tpcNSigmaHe() + tofNsigma*tofNsigma) > antinucleiNsigma)))
if ((track.pt() < antinucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaHe()) > antinucleiNsigma)) || (track.pt() > antinucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaHe() * track.tpcNSigmaHe() + tofNsigma * tofNsigma) > antinucleiNsigma)))
return false;
if (!singleSpeciesTPCNSigma(track, 3))
return false;
Expand Down Expand Up @@ -781,7 +781,7 @@
if (px == 0 && py != 0) {

uy = -(pz * pz) / py;
ux = sign * sqrt(py * py - (pz * pz * pz * pz) / (py * py));

Check warning on line 784 in PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
uz = pz;
u.SetXYZ(ux, uy, uz);
return;
Expand All @@ -791,7 +791,7 @@
if (py == 0 && px != 0) {

ux = -(pz * pz) / px;
uy = sign * sqrt(px * px - (pz * pz * pz * pz) / (px * px));

Check warning on line 794 in PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
uz = pz;
u.SetXYZ(ux, uy, uz);
return;
Expand All @@ -809,7 +809,7 @@
}

// Solutions
ux = (-b + sign * sqrt(delta)) / (2.0 * a);

Check warning on line 812 in PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
uy = (-pz * pz - px * ux) / py;
uz = pz;
u.SetXYZ(ux, uy, uz);
Expand Down Expand Up @@ -886,13 +886,13 @@
TVector3 particleDir(track.px(), track.py(), track.pz());
double deltaEtaJet = particleDir.Eta() - pJet.Eta();
double deltaPhiJet = getDeltaPhi(particleDir.Phi(), pJet.Phi());
double deltaRJet = sqrt(deltaEtaJet * deltaEtaJet + deltaPhiJet * deltaPhiJet);

Check warning on line 889 in PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
double deltaEtaUE1 = particleDir.Eta() - UEAxis1.Eta();
double deltaPhiUE1 = getDeltaPhi(particleDir.Phi(), UEAxis1.Phi());
double deltaRUE1 = sqrt(deltaEtaUE1 * deltaEtaUE1 + deltaPhiUE1 * deltaPhiUE1);

Check warning on line 892 in PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
double deltaEtaUE2 = particleDir.Eta() - UEAxis2.Eta();
double deltaPhiUE2 = getDeltaPhi(particleDir.Phi(), UEAxis2.Phi());
double deltaRUE2 = sqrt(deltaEtaUE2 * deltaEtaUE2 + deltaPhiUE2 * deltaPhiUE2);

Check warning on line 895 in PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.

if (deltaRJet < Rmax) {
if (deltaPhiJet != -999)
Expand Down
Loading