|
36 | 36 | #include "TVector3.h" |
37 | 37 | #include "TPDGCode.h" |
38 | 38 |
|
39 | | - |
40 | 39 | using namespace o2; |
41 | 40 | using namespace o2::framework; |
42 | 41 | using namespace o2::framework::expressions; |
@@ -343,7 +342,7 @@ struct AngularCorrelationsInJets { |
343 | 342 |
|
344 | 343 | template <class T> |
345 | 344 | bool singleSpeciesTPCNSigma(T const& track, int species) // make cut configurable |
346 | | - { // reject any track that has nsigma < 3 for more than 1 species |
| 345 | + { // reject any track that has nsigma < 3 for more than 1 species |
347 | 346 | if (track.tpcNSigmaStoreEl() < 3.0 || track.tpcNSigmaStoreMu() < 3.0 || track.tpcNSigmaStorePi() < 3.0 || track.tpcNSigmaStoreKa() < 3.0 || track.tpcNSigmaStoreTr() < 3.0 || track.tpcNSigmaStoreAl() < 3.0) |
348 | 347 | return false; |
349 | 348 | switch (species) { |
@@ -376,12 +375,12 @@ struct AngularCorrelationsInJets { |
376 | 375 |
|
377 | 376 | // nsigma |
378 | 377 | double tofNsigma = track.hasTOF() ? track.tofNSigmaPr() : 999; |
379 | | - if ((track.pt() < fProtonTPCTOFpT && (TMath::Abs(track.tpcNSigmaPr()) > fProtonNsigma)) || (track.pt() > fProtonTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaPr()*track.tpcNSigmaPr() + tofNsigma*tofNsigma) > fProtonNsigma))) |
380 | | - if (TMath::Sqrt(track.tpcNSigmaPr()*track.tpcNSigmaPr() + tofNsigma*tofNsigma) > fProtonNsigma) |
381 | | - return false; |
| 378 | + if ((track.pt() < fProtonTPCTOFpT && (TMath::Abs(track.tpcNSigmaPr()) > fProtonNsigma)) || (track.pt() > fProtonTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaPr() * track.tpcNSigmaPr() + tofNsigma * tofNsigma) > fProtonNsigma))) |
| 379 | + if (TMath::Sqrt(track.tpcNSigmaPr() * track.tpcNSigmaPr() + tofNsigma * tofNsigma) > fProtonNsigma) |
| 380 | + return false; |
382 | 381 | if (!singleSpeciesTPCNSigma(track, 1)) |
383 | 382 | return false; |
384 | | - } else { // for yields |
| 383 | + } else { // for yields |
385 | 384 | // DCA |
386 | 385 | if (TMath::Abs(track.dcaXY()) > fProtonDCAxyYield) |
387 | 386 | return false; |
@@ -428,7 +427,7 @@ struct AngularCorrelationsInJets { |
428 | 427 |
|
429 | 428 | // nsigma |
430 | 429 | double tofNsigma = track.hasTOF() ? track.tofNSigmaPr() : 999; |
431 | | - if ((track.pt() < fAntiprotonTPCTOFpT && (TMath::Abs(track.tpcNSigmaPr()) > fAntiprotonNsigma)) || (track.pt() > fAntiprotonTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaPr()*track.tpcNSigmaPr() + tofNsigma*tofNsigma) > fAntiprotonNsigma))) |
| 430 | + if ((track.pt() < fAntiprotonTPCTOFpT && (TMath::Abs(track.tpcNSigmaPr()) > fAntiprotonNsigma)) || (track.pt() > fAntiprotonTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaPr() * track.tpcNSigmaPr() + tofNsigma * tofNsigma) > fAntiprotonNsigma))) |
432 | 431 | return false; |
433 | 432 | if (!singleSpeciesTPCNSigma(track, 1)) |
434 | 433 | return false; |
@@ -479,7 +478,7 @@ struct AngularCorrelationsInJets { |
479 | 478 |
|
480 | 479 | // nsigma |
481 | 480 | double tofNsigma = track.hasTOF() ? track.tofNSigmaDe() : 999; |
482 | | - if ((track.pt() < fNucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaDe()) > fNucleiNsigma)) || (track.pt() > fNucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaDe()*track.tpcNSigmaDe() + tofNsigma*tofNsigma) > fNucleiNsigma))) |
| 481 | + if ((track.pt() < fNucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaDe()) > fNucleiNsigma)) || (track.pt() > fNucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaDe() * track.tpcNSigmaDe() + tofNsigma * tofNsigma) > fNucleiNsigma))) |
483 | 482 | return false; |
484 | 483 | if (!singleSpeciesTPCNSigma(track, 2)) |
485 | 484 | return false; |
@@ -519,7 +518,7 @@ struct AngularCorrelationsInJets { |
519 | 518 |
|
520 | 519 | // nsigma |
521 | 520 | double tofNsigma = track.hasTOF() ? track.tofNSigmaHe() : 999; |
522 | | - if ((track.pt() < fNucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaHe()) > fNucleiNsigma)) || (track.pt() > fNucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaHe()*track.tpcNSigmaHe() + tofNsigma*tofNsigma) > fNucleiNsigma))) |
| 521 | + if ((track.pt() < fNucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaHe()) > fNucleiNsigma)) || (track.pt() > fNucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaHe() * track.tpcNSigmaHe() + tofNsigma * tofNsigma) > fNucleiNsigma))) |
523 | 522 | return false; |
524 | 523 | if (!singleSpeciesTPCNSigma(track, 3)) |
525 | 524 | return false; |
@@ -572,7 +571,7 @@ struct AngularCorrelationsInJets { |
572 | 571 |
|
573 | 572 | // nsigma |
574 | 573 | double tofNsigma = track.hasTOF() ? track.tofNSigmaDe() : 999; |
575 | | - if ((track.pt() < fAntinucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaDe()) > fAntinucleiNsigma)) || (track.pt() > fAntinucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaDe()*track.tpcNSigmaDe() + tofNsigma*tofNsigma) > fAntinucleiNsigma))) |
| 574 | + if ((track.pt() < fAntinucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaDe()) > fAntinucleiNsigma)) || (track.pt() > fAntinucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaDe() * track.tpcNSigmaDe() + tofNsigma * tofNsigma) > fAntinucleiNsigma))) |
576 | 575 | return false; |
577 | 576 | if (!singleSpeciesTPCNSigma(track, 2)) |
578 | 577 | return false; |
@@ -610,7 +609,7 @@ struct AngularCorrelationsInJets { |
610 | 609 |
|
611 | 610 | // nsigma |
612 | 611 | double tofNsigma = track.hasTOF() ? track.tofNSigmaHe() : 999; |
613 | | - if ((track.pt() < fAntinucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaHe()) > fAntinucleiNsigma)) || (track.pt() > fAntinucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaHe()*track.tpcNSigmaHe() + tofNsigma*tofNsigma) > fAntinucleiNsigma))) |
| 612 | + if ((track.pt() < fAntinucleiTPCTOFpT && (TMath::Abs(track.tpcNSigmaHe()) > fAntinucleiNsigma)) || (track.pt() > fAntinucleiTPCTOFpT && (TMath::Sqrt(track.tpcNSigmaHe() * track.tpcNSigmaHe() + tofNsigma * tofNsigma) > fAntinucleiNsigma))) |
614 | 613 | return false; |
615 | 614 | if (!singleSpeciesTPCNSigma(track, 3)) |
616 | 615 | return false; |
|
0 commit comments