Skip to content

PWGUD: Analysis in photonuclear events #5095

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
rolavick merged 19 commits into from
Apr 25, 2024
Merged

PWGUD: Analysis in photonuclear events #5095

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
19 commits
Select commit Hold shift + click to select a range
6e3df8c
PWGUD: Analysis in photonuclear events
JosueMtzGar Mar 11, 2024
08afae7
formatting issues
JosueMtzGar Mar 11, 2024
7d4e9be
update with process switch
JosueMtzGar Mar 11, 2024
a5f22f0
;
JosueMtzGar Mar 12, 2024
447618f
Update
JosueMtzGar Mar 12, 2024
b2c43d4
format issue
JosueMtzGar Mar 13, 2024
7c03a63
Merge branch 'master' into localtest
JosueMtzGar Mar 19, 2024
ab7e666
new changes
JosueMtzGar Apr 14, 2024
880eec7
new changes
JosueMtzGar Apr 16, 2024
3584112
new changes
JosueMtzGar Apr 16, 2024
4981245
new changes
JosueMtzGar Apr 18, 2024
5203ef4
new changes
JosueMtzGar Apr 22, 2024
81b013c
new changes
JosueMtzGar Apr 22, 2024
f89e1d1
new changes
JosueMtzGar Apr 22, 2024
5e8054f
new changes
JosueMtzGar Apr 23, 2024
ea1bffc
new changes
JosueMtzGar Apr 23, 2024
30f32d1
new changes
JosueMtzGar Apr 23, 2024
5200557
new changes
JosueMtzGar Apr 23, 2024
552f07f
new changes
JosueMtzGar Apr 23, 2024
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
206 changes: 165 additions & 41 deletions PWGMM/Lumi/Tasks/lumiStability.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,14 +16,18 @@
#include "Framework/runDataProcessing.h"
#include "Framework/AnalysisTask.h"
#include "Framework/AnalysisDataModel.h"
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/TrackSelectionTables.h"
#include "Framework/ASoAHelpers.h"
#include "DataFormatsFDD/Digit.h"
#include "Framework/ASoA.h"

using namespace o2;
using namespace o2::framework;
using namespace o2::aod::evsel;

using BCsWithTimestamps = soa::Join<aod::BCs, aod::Timestamps>;
using BCsRun3 = soa::Join<aod::BCs, aod::Timestamps, aod::BcSels, aod::Run3MatchedToBCSparse>;

int nBCsPerOrbit = 3564;

Expand All @@ -33,39 +37,42 @@ struct lumiStabilityTask {

void init(InitContext const&)
{
const AxisSpec axisFDDTriggger{nBCsPerOrbit, -0.5f, nBCsPerOrbit - 0.5f};
const AxisSpec axisFT0Triggger{nBCsPerOrbit, -0.5f, nBCsPerOrbit - 0.5f};
const AxisSpec axisFV0Triggger{nBCsPerOrbit, -0.5f, nBCsPerOrbit - 0.5f};
const AxisSpec axisTriggger{nBCsPerOrbit, -0.5f, nBCsPerOrbit - 0.5f};

// histo about triggers
histos.add("FDD/bcVertexTrigger", "vertex trigger per BC (FDD);BC in FDD; counts", kTH1F, {axisFDDTriggger});
histos.add("FDD/bcVertexTriggerCoincidence", "vertex trigger per BC (FDD) with coincidences;BC in FDD; counts", kTH1F, {axisFDDTriggger});
histos.add("FDD/bcSCentralTrigger", "scentral trigger per BC (FDD);BC in FDD; counts", kTH1F, {axisFDDTriggger});
histos.add("FDD/bcSCentralTriggerCoincidence", "scentral trigger per BC (FDD) with coincidences;BC in FDD; counts", kTH1F, {axisFDDTriggger});
histos.add("FDD/bcVSCTrigger", "vertex and scentral trigger per BC (FDD);BC in FDD; counts", kTH1F, {axisFDDTriggger});
histos.add("FDD/bcVSCTriggerCoincidence", "vertex and scentral trigger per BC (FDD) with coincidences;BC in FDD; counts", kTH1F, {axisFDDTriggger});
histos.add("FDD/bcCentralTrigger", "central trigger per BC (FDD);BC in FDD; counts", kTH1F, {axisFDDTriggger});
histos.add("FDD/bcCentralTriggerCoincidence", "central trigger per BC (FDD) with coincidences;BC in FDD; counts", kTH1F, {axisFDDTriggger});
histos.add("FDD/bcVCTrigger", "vertex and central trigger per BC (FDD);BC in FDD; counts", kTH1F, {axisFDDTriggger});
histos.add("FDD/bcVCTriggerCoincidence", "vertex and central trigger per BC (FDD) with coincidences;BC in FDD; counts", kTH1F, {axisFDDTriggger});

histos.add("FT0/bcVertexTrigger", "vertex trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisFT0Triggger});
histos.add("FT0/bcSCentralTrigger", "Scentral trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisFT0Triggger});
histos.add("FT0/bcVSCTrigger", "vertex and Scentral trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisFT0Triggger});
histos.add("FT0/bcCentralTrigger", "central trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisFT0Triggger});
histos.add("FT0/bcVCTrigger", "vertex and central trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisFT0Triggger});
histos.add("FT0/bcSCentralCentralTrigger", "Scentral and central trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisFT0Triggger});

histos.add("FV0/bcOutTrigger", "Out trigger per BC (FV0);BC in V0; counts", kTH1F, {axisFV0Triggger});
histos.add("FV0/bcInTrigger", "In trigger per BC (FV0);BC in V0; counts", kTH1F, {axisFV0Triggger});
histos.add("FV0/bcSCenTrigger", "SCen trigger per BC (FV0);BC in V0; counts", kTH1F, {axisFV0Triggger});
histos.add("FV0/bcCenTrigger", "Out trigger per BC (FV0);BC in V0; counts", kTH1F, {axisFV0Triggger});
histos.add("FV0/bcSCenCenTrigger", "SCen and Cen trigger per BC (FV0);BC in V0; counts", kTH1F, {axisFV0Triggger});
histos.add("FDD/bcVertexTrigger", "vertex trigger per BC (FDD);BC in FDD; counts", kTH1F, {axisTriggger});
histos.add("FDD/bcVertexTriggerCoincidence", "vertex trigger per BC (FDD) with coincidences;BC in FDD; counts", kTH1F, {axisTriggger});
histos.add("FDD/bcSCentralTrigger", "scentral trigger per BC (FDD);BC in FDD; counts", kTH1F, {axisTriggger});
histos.add("FDD/bcSCentralTriggerCoincidence", "scentral trigger per BC (FDD) with coincidences;BC in FDD; counts", kTH1F, {axisTriggger});
histos.add("FDD/bcVSCTrigger", "vertex and scentral trigger per BC (FDD);BC in FDD; counts", kTH1F, {axisTriggger});
histos.add("FDD/bcVSCTriggerCoincidence", "vertex and scentral trigger per BC (FDD) with coincidences;BC in FDD; counts", kTH1F, {axisTriggger});
histos.add("FDD/bcCentralTrigger", "central trigger per BC (FDD);BC in FDD; counts", kTH1F, {axisTriggger});
histos.add("FDD/bcCentralTriggerCoincidence", "central trigger per BC (FDD) with coincidences;BC in FDD; counts", kTH1F, {axisTriggger});
histos.add("FDD/bcVCTrigger", "vertex and central trigger per BC (FDD);BC in FDD; counts", kTH1F, {axisTriggger});
histos.add("FDD/bcVCTriggerCoincidence", "vertex and central trigger per BC (FDD) with coincidences;BC in FDD; counts", kTH1F, {axisTriggger});

histos.add("FT0/bcVertexTrigger", "vertex trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisTriggger});
histos.add("FT0/bcSCentralTrigger", "Scentral trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisTriggger});
histos.add("FT0/bcVSCTrigger", "vertex and Scentral trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisTriggger});
histos.add("FT0/bcCentralTrigger", "central trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisTriggger});
histos.add("FT0/bcVCTrigger", "vertex and central trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisTriggger});
histos.add("FT0/bcSCentralCentralTrigger", "Scentral and central trigger per BC (FT0);BC in FT0; counts", kTH1F, {axisTriggger});

histos.add("FV0/bcOutTrigger", "Out trigger per BC (FV0);BC in V0; counts", kTH1F, {axisTriggger});
histos.add("FV0/bcInTrigger", "In trigger per BC (FV0);BC in V0; counts", kTH1F, {axisTriggger});
histos.add("FV0/bcSCenTrigger", "SCen trigger per BC (FV0);BC in V0; counts", kTH1F, {axisTriggger});
histos.add("FV0/bcCenTrigger", "Out trigger per BC (FV0);BC in V0; counts", kTH1F, {axisTriggger});
histos.add("FV0/bcSCenCenTrigger", "SCen and Cen trigger per BC (FV0);BC in V0; counts", kTH1F, {axisTriggger});
}

bool checkAnyCoincidence(const std::vector<int>& channels)
{
std::map<int, int> channelPairs = {{0, 4}, {1, 5}, {2, 6}, {4, 7}};
constexpr std::pair<int, int> pair0 = {0, 4};
constexpr std::pair<int, int> pair1 = {1, 5};
constexpr std::pair<int, int> pair2 = {2, 6};
constexpr std::pair<int, int> pair3 = {3, 7};
constexpr std::array<std::pair<int, int>, 4> channelPairs = {pair0, pair1, pair2, pair3};
// std::map<int, int> channelPairs = {{0, 4}, {1, 5}, {2, 6}, {3, 7}};
for (const auto& pair : channelPairs) {
if (std::find(channels.begin(), channels.end(), pair.first) != channels.end() &&
std::find(channels.begin(), channels.end(), pair.second) != channels.end()) {
Expand All @@ -75,17 +82,59 @@ struct lumiStabilityTask {
return false;
}

void processFDDFT0(aod::FT0s const& ft0s, aod::FDDs const& fdds, aod::BCsWithTimestamps const&)
void processMain(aod::FT0s const& ft0s, aod::FDDs const& fdds, aod::FV0As const& fv0s, BCsRun3 const& bcs)
{
int CountNormal(0), CountPastProtec(0);
for (auto const& fdd : fdds) {
auto bc = fdd.bc_as<BCsWithTimestamps>();
if (bc.timestamp() == false) {
auto bc = fdd.bc_as<BCsRun3>();
if (bc.timestamp() == 0) {
continue;
}

Long64_t globalBC = bc.globalBC();
int localBC = globalBC % nBCsPerOrbit;

int deltaIndex = 0; // backward move counts
int deltaBC = 0; // current difference wrt globalBC
int maxDeltaBC = 5; // maximum difference
bool pastActivityFDD = false;
while (deltaBC < maxDeltaBC) {
deltaIndex++;
if (bc.globalIndex() - deltaIndex < 0) {
break;
}
const auto& bc_past = bcs.iteratorAt(bc.globalIndex() - deltaIndex);
deltaBC = globalBC - bc_past.globalBC();
if (deltaBC < maxDeltaBC) {
pastActivityFDD |= bc_past.has_fdd();
}
}
deltaIndex = 0;
deltaBC = 0;

bool futureActivityFDD = false;
while (deltaBC < maxDeltaBC) {
deltaIndex++;
if (bc.globalIndex() + deltaIndex >= bcs.size()) {
break;
}
const auto& bc_future = bcs.iteratorAt(bc.globalIndex() + deltaIndex);
deltaBC = bc_future.globalBC() - globalBC;
if (deltaBC < maxDeltaBC) {
futureActivityFDD |= bc_future.has_fdd();
}
}

CountNormal++;
if (pastActivityFDD == true || futureActivityFDD == true) {
CountPastProtec++;
continue;
}
/*if (pastActivityFDD == true) {
CountPastProtec++;
continue;
}*/

std::bitset<8> fddTriggers = fdd.triggerMask();
bool vertex = fddTriggers[o2::fdd::Triggers::bitVertex];
bool scentral = fddTriggers[o2::fdd::Triggers::bitSCen];
Expand Down Expand Up @@ -144,14 +193,52 @@ struct lumiStabilityTask {
} // loop over FDD events

for (auto const& ft0 : ft0s) {
auto bc = ft0.bc_as<BCsWithTimestamps>();
if (bc.timestamp() == false) {
auto bc = ft0.bc_as<BCsRun3>();
if (bc.timestamp() == 0) {
continue;
}

Long64_t globalBC = bc.globalBC();
int localBC = globalBC % nBCsPerOrbit;

int deltaIndex = 0; // backward move counts
int deltaBC = 0; // current difference wrt globalBC
int maxDeltaBC = 5; // maximum difference
bool pastActivityFT0 = false;
while (deltaBC < maxDeltaBC) {
deltaIndex++;
if (bc.globalIndex() - deltaIndex < 0) {
break;
}
const auto& bc_past = bcs.iteratorAt(bc.globalIndex() - deltaIndex);
deltaBC = globalBC - bc_past.globalBC();
if (deltaBC < maxDeltaBC) {
pastActivityFT0 |= bc_past.has_ft0();
}
}
deltaIndex = 0;
deltaBC = 0;

bool futureActivityFT0 = false;
while (deltaBC < maxDeltaBC) {
deltaIndex++;
if (bc.globalIndex() + deltaIndex >= bcs.size()) {
break;
}
const auto& bc_future = bcs.iteratorAt(bc.globalIndex() + deltaIndex);
deltaBC = bc_future.globalBC() - globalBC;
if (deltaBC < maxDeltaBC) {
futureActivityFT0 |= bc_future.has_ft0();
}
}

if (pastActivityFT0 == true || futureActivityFT0 == true) {
continue;
}
/*if (pastActivityFT0 == true) {
continue;
}*/

std::bitset<8> fT0Triggers = ft0.triggerMask();
bool vertex = fT0Triggers[o2::fdd::Triggers::bitVertex];

Expand Down Expand Up @@ -179,21 +266,55 @@ struct lumiStabilityTask {
}
}
} // loop over FT0 events
} // end processFDDFT0

PROCESS_SWITCH(lumiStabilityTask, processFDDFT0, "Process FDD and FT0 to lumi stability analysis", true);

void processV0(aod::FV0As const& fv0s, aod::BCsWithTimestamps const&)
{
for (auto const& fv0 : fv0s) {
auto bc = fv0.bc_as<BCsWithTimestamps>();
if (bc.timestamp() == false) {
auto bc = fv0.bc_as<BCsRun3>();
if (bc.timestamp() == 0) {
continue;
}

Long64_t globalBC = bc.globalBC();
int localBC = globalBC % nBCsPerOrbit;

int deltaIndex = 0; // backward move counts
int deltaBC = 0; // current difference wrt globalBC
int maxDeltaBC = 5; // maximum difference
bool pastActivityV0A = false;
while (deltaBC < maxDeltaBC) {
deltaIndex++;
if (bc.globalIndex() - deltaIndex < 0) {
break;
}
const auto& bc_past = bcs.iteratorAt(bc.globalIndex() - deltaIndex);
deltaBC = globalBC - bc_past.globalBC();
if (deltaBC < maxDeltaBC) {
pastActivityV0A |= bc_past.has_fv0a();
}
}
deltaIndex = 0;
deltaBC = 0;

bool futureActivityV0A = false;
while (deltaBC < maxDeltaBC) {
deltaIndex++;
if (bc.globalIndex() + deltaIndex >= bcs.size()) {
break;
}
const auto& bc_future = bcs.iteratorAt(bc.globalIndex() + deltaIndex);
deltaBC = bc_future.globalBC() - globalBC;
if (deltaBC < maxDeltaBC) {
futureActivityV0A |= bc_future.has_fv0a();
}
}

if (pastActivityV0A == true || futureActivityV0A == true) {
continue;
}

/*if (pastActivityV0A == true) {
continue;
}*/

std::bitset<8> fv0Triggers = fv0.triggerMask();
bool aOut = fv0Triggers[o2::fdd::Triggers::bitAOut];
bool aIn = fv0Triggers[o2::fdd::Triggers::bitAIn];
Expand All @@ -219,9 +340,12 @@ struct lumiStabilityTask {
}
}
} // loop over V0 events
} // end processV0
std::cout << "************ >>>>>>>>>>>>>> "
<< "Whithout Past Protection: " << CountNormal << " "
<< "Avoided Whith Past Protection: " << CountPastProtec << "<<<<<<<<<<<< ********************" << std::endl;
} // end processMain

PROCESS_SWITCH(lumiStabilityTask, processV0, "Process V0 to lumi stability analysis", true);
PROCESS_SWITCH(lumiStabilityTask, processMain, "Process FDD and FT0 to lumi stability analysis", true);
};

WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
Expand Down
4 changes: 3 additions & 1 deletion PWGMM/Mult/DataModel/ReducedTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,6 +11,8 @@

#ifndef PWGMM_MULT_DATAMODEL_REDUCEDTABLES_H_
#define PWGMM_MULT_DATAMODEL_REDUCEDTABLES_H_
#include <vector>

#include "Framework/AnalysisDataModel.h"
#include "Common/DataModel/Multiplicity.h"
#include "Common/DataModel/Centrality.h"
Expand All @@ -34,7 +36,7 @@ namespace rcol
{
DECLARE_SOA_INDEX_COLUMN(RBC, rbc);
DECLARE_SOA_COLUMN(MapEtaPhi, mapetaphi, std::vector<int>);
}
} // namespace rcol

#define Ccols o2::soa::Index<>, \
rcol::RBCId, \
Expand Down
5 changes: 5 additions & 0 deletions PWGUD/Tasks/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -92,3 +92,8 @@ o2physics_add_dpl_workflow(polarisation-rho
SOURCES PolarisationRho.cxx
PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2Physics::DGPIDSelector
COMPONENT_NAME Analysis)

o2physics_add_dpl_workflow(upc-photonuclear-jmg
SOURCES upcPhotonuclearAnalysisJMG.cxx
PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2::ReconstructionDataFormats O2::DetectorsBase O2::DetectorsCommonDataFormats
COMPONENT_NAME Analysis)
Loading