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
12 changes: 2 additions & 10 deletions PWGLF/Tasks/Nuspex/multiplicityPt.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -138,7 +138,7 @@ struct MultiplicityPt {

// Custom track cuts matching spectraTOF
TrackSelection customTrackCuts;

// TF1 pointers for phi cuts
TF1* fphiCutLow = nullptr;
TF1* fphiCutHigh = nullptr;
Expand Down Expand Up @@ -178,7 +178,6 @@ struct MultiplicityPt {
static constexpr int PDGKaon = kKPlus;
static constexpr int PDGProton = kProton;


// Get magnetic field from CCDB
int getMagneticField(uint64_t timestamp)
{
Expand Down Expand Up @@ -335,7 +334,7 @@ struct MultiplicityPt {

if (!passedNClTPCPIDCut(track))
return false;

// Add phi cut with magnetic field
if (!passedPhiCut(track, magField))
return false;
Expand Down Expand Up @@ -411,7 +410,6 @@ struct MultiplicityPt {
return true;
}


void processData(CollisionTableData::iterator const& collision,
TrackTableData const& tracks,
BCsRun3 const& bcs);
Expand All @@ -426,7 +424,6 @@ struct MultiplicityPt {
BCsRun3 const& bcs);
PROCESS_SWITCH(MultiplicityPt, processMC, "process MC", true);


void init(InitContext const&);

void endOfStream(EndOfStreamContext& /*eos*/)
Expand All @@ -441,13 +438,11 @@ struct MultiplicityPt {
}
};


WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
{
return WorkflowSpec{adaptAnalysisTask<MultiplicityPt>(cfgc)};
}


void MultiplicityPt::init(InitContext const&)
{
LOG(info) << "==================================================";
Expand Down Expand Up @@ -759,7 +754,6 @@ void MultiplicityPt::processMC(TrackTableMC const& tracks,
LOG(info) << "Total collision labels: " << labels.size();
LOG(info) << "Total centrality entries: " << centTable.size();


LOG(info) << "\n=== CENTRALITY DEBUG - RAW DATA ===";
LOG(info) << "First 20 centrality values from centTable:";
int debugCount = 0;
Expand Down Expand Up @@ -790,7 +784,6 @@ void MultiplicityPt::processMC(TrackTableMC const& tracks,
LOG(info) << "Checking if centrality might be inverted...";
LOG(info) << "Will check correlation with multiplicity in the next step.";


std::map<int64_t, int> mcCollisionToNch;
std::map<int64_t, float> mcCollisionVz;
std::set<int64_t> physicsSelectedMCCollisions;
Expand Down Expand Up @@ -908,7 +901,6 @@ void MultiplicityPt::processMC(TrackTableMC const& tracks,
LOG(info) << "recoToMcMap size: " << recoToMcMap.size();
LOG(info) << "recoToCentMap size: " << recoToCentMap.size();


LOG(info) << "\n=== CENTRALITY VS MULTIPLICITY DEBUG ===";

// Create temporary vectors to check correlation
Expand Down