Skip to content
Closed
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
35 changes: 18 additions & 17 deletions PWGUD/AQC/udQC.cxx
Original file line number Diff line number Diff line change
@@ -1,3 +1,3 @@
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
Expand All @@ -12,7 +12,7 @@
/// \brief A task for Asynchronus Quality Control for Ultra-perimpheral and Diffraction (AQC-UD)
/// \author Anisa Khatun, anisa.khatun@cern.ch
/// \author Paul Buehler, paul.buehler@oeaw.ac.at
/// \since 04.05.2023

Check warning on line 15 in PWGUD/AQC/udQC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[doc/file]

Documentation for \file is missing, incorrect or misplaced.

#include "Framework/runDataProcessing.h"
#include "Framework/AnalysisTask.h"
Expand Down Expand Up @@ -45,7 +45,7 @@

// get a DGCutparHolder
DGCutparHolder diffCuts = DGCutparHolder();
Configurable<DGCutparHolder> DGCuts{"DGCuts", {}, "DG event cuts"};

Check warning on line 48 in PWGUD/AQC/udQC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
Configurable<bool> withAmbTrackAnalysis{"ambiguousTracks", false, "with ambiguous tracks analysis"};
Configurable<bool> withAmbFwdTrackAnalysis{"ambiguousFwdTracks", false, "with ambiguous forward tracks analysis"};
Configurable<bool> doCleanFITBC{"doCleanFITBC", false, "Require cleanFIT in compatible BCs"};
Expand Down Expand Up @@ -189,15 +189,16 @@
TCs const& tracks, FWs const& fwdtracks, ATs const& /*ambtracks*/, AFTs const& /*ambfwdtracks*/,
aod::FT0s const& /*ft0s*/, aod::FV0As const& /*fv0as*/, aod::FDDs const& /*fdds*/,
aod::Zdcs& zdcs, aod::Calos& calos,
aod::V0s const& v0s, aod::Cascades const& cascades)

Check warning on line 192 in PWGUD/AQC/udQC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-process]

Argument aod::Calos& calos is not const&.

Check warning on line 192 in PWGUD/AQC/udQC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-process]

Argument aod::Zdcs& zdcs is not const&.
{
//LOGF(debug, "<UDQC. Collision %d", collision.globalIndex());
//LOGF(debug, "<UDQC> Start %i", abcrs.size());
// LOGF(debug, "<UDQC. Collision %d", collision.globalIndex());
// LOGF(debug, "<UDQC> Start %i", abcrs.size());

if (!tracks.size())
return;
if (collision.numContrib() < 1)
return;


if (!tracks.size()) return;
if (collision.numContrib() < 1) return;

bool isDGcandidate = true;
registry.get<TH1>(HIST("collisions/Stat"))->Fill(0., isDGcandidate * 1.);

Expand Down Expand Up @@ -252,15 +253,15 @@

if (track.tpcSignal() > maxdEdxTPC) {
maxdEdxTPC = track.tpcSignal();
// LOGF(debug, "<UDQC> New maxdEdx TPC %f", maxdEdxTPC);
// LOGF(debug, "<UDQC> New maxdEdx TPC %f", maxdEdxTPC);
}

// TOF hit?
if (track.hasTOF()) {
registry.get<TH2>(HIST("tracks/dEdxTOF"))->Fill(track.p() / track.sign(), track.beta());
if (track.tofSignal() > maxdEdxTOF) {
maxdEdxTOF = track.tofSignal();
//LOGF(debug, "<UDQC> New maxdEdx TOF %f", maxdEdxTOF);
// LOGF(debug, "<UDQC> New maxdEdx TOF %f", maxdEdxTOF);
}

// No vertex track with TOF hit?
Expand All @@ -285,7 +286,7 @@
if (collision.numContrib() > 0) {
rgtrwTOF /= collision.numContrib();
}
//LOGF(debug, "<UDQC> PV tracks with TOF: %f [1]", rgtrwTOF);
// LOGF(debug, "<UDQC> PV tracks with TOF: %f [1]", rgtrwTOF);
registry.get<TH2>(HIST("collisions/tResvsrTOFTracks"))->Fill(collision.collisionTimeRes(), rgtrwTOF);
registry.get<TH2>(HIST("collisions/tResvsTOFTrkNoPV"))->Fill(collision.collisionTimeRes(), norgtrwTOF);

Expand Down Expand Up @@ -433,7 +434,7 @@
}
// define Lorentz vector to create invariant mass
lvtmp.SetPtEtaPhiM(track.pt(), track.eta(), track.phi(), mass2Use);
//LOGF(debug, "mass %f track pt %f/%f eta %f/%f", mass2Use, track.pt(), lvtmp.Perp(), track.eta(), lvtmp.Eta());
// LOGF(debug, "mass %f track pt %f/%f eta %f/%f", mass2Use, track.pt(), lvtmp.Perp(), track.eta(), lvtmp.Eta());
if (track.pt() <= diffCuts.minPt() || track.pt() >= diffCuts.maxPt()) {
goodpts = false;
}
Expand Down Expand Up @@ -525,7 +526,7 @@
registry.get<TH2>(HIST("DG/etaminus"))->Fill(track.eta(), track.phi());
}

// LOGF(debug, "dEdx TPC %f TOF %i %f", track.tpcSignal(), track.hasTOF(), track.hasTOF() ? track.tofSignal() : 0.);
// LOGF(debug, "dEdx TPC %f TOF %i %f", track.tpcSignal(), track.hasTOF(), track.hasTOF() ? track.tofSignal() : 0.);
if (collision.numContrib() == 2) {
registry.get<TH2>(HIST("DG/dEdxTPC"))->Fill(track.tpcInnerParam() / track.sign(), track.tpcSignal());
registry.get<TH1>(HIST("DG/trkDCAxy"))->Fill(track.dcaXY());
Expand Down Expand Up @@ -566,11 +567,11 @@
void processCleanFIT1(CC const& collision, BCs const& bct0s,
aod::FT0s const& /*ft0s*/, aod::FV0As const& /*fv0as*/, aod::FDDs const& /*fdds*/)
{
//LOGF(debug, "<CleanFit. Collision %d", collision.globalIndex());
// LOGF(debug, "<CleanFit. Collision %d", collision.globalIndex());

// test influence of BCrange width using a series of NDtcoll
float ampFV0A, ampFT0A, ampFT0C, ampFDDA, ampFDDC;
auto FITlims = std::vector<float>(5, 1000000.);

Check warning on line 574 in PWGUD/AQC/udQC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
bool isDGcandidate = true;
for (int NDtcoll = 0; NDtcoll < 20; NDtcoll++) {
auto bcSlice = udhelpers::compatibleBCs(collision, NDtcoll, bct0s, 0);
Expand Down Expand Up @@ -622,7 +623,7 @@
void processCleanFIT2(CC const& collision, BCs const& bct0s,
aod::FT0s const& /*ft0s*/, aod::FV0As const& /*fv0as*/, aod::FDDs const& /*fdds*/)
{
//LOGF(debug, "<CleanFit. Collision %d", collision.globalIndex());
// LOGF(debug, "<CleanFit. Collision %d", collision.globalIndex());
uint64_t bcnum = 0;
if (collision.has_foundBC()) {
auto collbc = collision.foundBC_as<BCs>();
Expand Down Expand Up @@ -681,12 +682,12 @@
// ...............................................................................................................
void processFV0(aod::FV0As const& fv0s, BCs const&)
{
//LOGF(info, "<FV0Signals> %d", fv0s.size());
// LOGF(info, "<FV0Signals> %d", fv0s.size());
if (fv0s.size() <= 0) {
return;
}

for (auto fv0 : fv0s) {

Check warning on line 690 in PWGUD/AQC/udQC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
registry.get<TH1>(HIST("FV0/hV0A"))->Fill(fv0.time());
// side A
for (size_t ind = 0; ind < fv0.channel().size(); ind++) {
Expand All @@ -699,7 +700,7 @@
// ...............................................................................................................
void processFT0(aod::FT0s const& ft0s, aod::FT0sCorrected const& ft0scorr, BCs const&)
{
//LOGF(debug, "<processFT0> %d", ft0s.size());
// LOGF(debug, "<processFT0> %d", ft0s.size());
for (auto const& collision : ft0scorr) {

if (collision.t0ACorrectedValid()) {
Expand All @@ -713,7 +714,7 @@
registry.get<TH1>(HIST("FT0/hT0AC"))->Fill(collision.t0AC());
}
}
for (auto ft0 : ft0s) {

Check warning on line 717 in PWGUD/AQC/udQC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
registry.get<TH1>(HIST("FT0/hT0A"))->Fill(ft0.timeA());
registry.get<TH1>(HIST("FT0/hT0C"))->Fill(ft0.timeC());

Expand All @@ -733,9 +734,9 @@
// ...............................................................................................................
void processFDD(aod::FDDs const& fdds, BCs const&)
{
//LOGF(debug, "<FDDSignals> %d", fdds.size());
// LOGF(debug, "<FDDSignals> %d", fdds.size());

for (auto fdd : fdds) {

Check warning on line 739 in PWGUD/AQC/udQC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.

registry.get<TH1>(HIST("FDD/hFDDA"))->Fill(fdd.timeA());
registry.get<TH1>(HIST("FDD/hFDDC"))->Fill(fdd.timeC());
Expand All @@ -755,7 +756,7 @@
// ...............................................................................................................
void processZDC(aod::Zdc const& zdc)
{
//LOGF(debug, "<ZDCSignals> %d", zdc.size());
// LOGF(debug, "<ZDCSignals> %d", zdc.size());

// Zdc energies
registry.get<TH2>(HIST("ZdcEnergies"))->Fill(0., zdc.energyZEM1());
Expand Down
Loading