Skip to content

Commit

Permalink
cfnetcharfluct
Browse files Browse the repository at this point in the history
  • Loading branch information
@nidamalikk
nidamalikk committed Jan 25, 2025
1 parent 1939ff9 commit e15b3fa
Showing 1 changed file with 107 additions and 102 deletions.
209 changes: 107 additions & 102 deletions PWGCF/EbyEFluctuations/Tasks/NetchargeFluctuations.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -13,8 +13,6 @@
// Department of Physics, Aligarh Muslim University, India
// to study the net charge fluctuations by observable, #nu_dyn



#include "Framework/AnalysisTask.h"
#include "Framework/runDataProcessing.h"
#include "Common/DataModel/EventSelection.h"
Expand All @@ -24,146 +22,153 @@
#include "Common/Core/TrackSelection.h"
#include "Common/DataModel/TrackSelectionTables.h"


using namespace o2;
using namespace o2::framework;
using namespace o2::framework::expressions;
using namespace std;


namespace o2::aod
{
<<<<<<< HEAD
using MyCollisions = soa::Join<aod::Collisions,aod::EvSels,aod::CentRun2V0Ms,aod::Mults>;
using MyCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::CentRun2V0Ms, aod::Mults>;
=======
using MyCollisions = soa::Join<aod::Collisions,aod::EvSels,aod::CentRun2V0Ms,aod::Mults>;
using MyCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::CentRun2V0Ms, aod::Mults>;
>>>>>>> 03c805b2635be6256d2a2a4ad4e9fa12fbdd067d
using MyCollision = MyCollisions::iterator;
using MyCollision = MyCollisions::iterator;

using MyTracks = soa::Join<aod::Tracks,aod::TracksExtra,aod::TracksDCA,aod::TrackSelection>;
using MyTrack = MyTracks::iterator;
using MyTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection>;
using MyTrack = MyTracks::iterator;
<<<<<<< HEAD
}
=======
}
>>>>>>> 03c805b2635be6256d2a2a4ad4e9fa12fbdd067d

struct NetchargeFluctuations
{
struct NetchargeFluctuations {
HistogramRegistry histos{"Histos", {}, OutputObjHandlingPolicy::AnalysisObject};
void init(o2::framework::InitContext&)
{
void init(o2::framework::InitContext&)
{
<<<<<<< HEAD
AxisSpec vtxZAxis = {100, -20, 20, "Z (cm)"};
AxisSpec dcaAxis = {1000, -100, 100, "DCA_{xy} (cm)"};
AxisSpec dcazAxis = {1000, -100, 100, "DCA_{z} (cm)"};
AxisSpec ptAxis = {40, 0.0, 4.0, "#it{p}_{T} (GeV/#it{c})"};
AxisSpec etaAxis = {30, -1.5, 1.5, "#eta"};
AxisSpec centAxis = {100, 0., 100., "centrality"};
AxisSpec multAxis = {2000, 0., 2000., "multiplicity"};

histos.add("hVertexZ_bef", "", kTH1F, {vtxZAxis});
histos.add("hVertexZ_aft", "", kTH1F, {vtxZAxis});
histos.add("hVertexZ_aft_sel", "", kTH1D, {vtxZAxis});
histos.add("hDCAxy_bef", "", kTH1D, {dcaAxis});
histos.add("hDCAxy_aft", "", kTH1D, {dcaAxis});
histos.add("hDCAz_bef", "", kTH1D, {dcazAxis});
histos.add("hDCAz_aft", "", kTH1D, {dcazAxis});
histos.add("hCentrality", "", kTH1D, {centAxis});
histos.add("hMultiplicity", "", kTH1D, {multAxis});
histos.add("hEta", "", kTH1F, {etaAxis});
histos.add("hPt", "", kTH1F, {ptAxis});
AxisSpec vtxZAxis = {100, -20, 20, "Z (cm)"};
AxisSpec dcaAxis = {1000, -100, 100, "DCA_{xy} (cm)"};
AxisSpec dcazAxis = {1000, -100, 100, "DCA_{z} (cm)"};
AxisSpec ptAxis = {40, 0.0, 4.0, "#it{p}_{T} (GeV/#it{c})"};
AxisSpec etaAxis = {30, -1.5, 1.5, "#eta"};
AxisSpec centAxis = {100, 0., 100., "centrality"};
AxisSpec multAxis = {2000, 0., 2000., "multiplicity"};

histos.add("hVertexZ_bef", "", kTH1F, {vtxZAxis});
histos.add("hVertexZ_aft", "", kTH1F, {vtxZAxis});
histos.add("hVertexZ_aft_sel", "", kTH1D, {vtxZAxis});
histos.add("hDCAxy_bef", "", kTH1D, {dcaAxis});
histos.add("hDCAxy_aft", "", kTH1D, {dcaAxis});
histos.add("hDCAz_bef", "", kTH1D, {dcazAxis});
histos.add("hDCAz_aft", "", kTH1D, {dcazAxis});
histos.add("hCentrality", "", kTH1D, {centAxis});
histos.add("hMultiplicity", "", kTH1D, {multAxis});
histos.add("hEta", "", kTH1F, {etaAxis});
histos.add("hPt", "", kTH1F, {ptAxis});
=======
AxisSpec vtxZAxis = {100, -20, 20, "Z (cm)"};
AxisSpec dcaAxis = {1000, -100, 100, "DCA_{xy} (cm)"};
AxisSpec dcazAxis = {1000, -100, 100, "DCA_{z} (cm)"};
AxisSpec ptAxis = {40, 0.0, 4.0, "#it{p}_{T} (GeV/#it{c})"};
AxisSpec etaAxis = {30, -1.5, 1.5, "#eta"};
AxisSpec centAxis = {100, 0., 100., "centrality"};
AxisSpec multAxis = {2000, 0., 2000., "multiplicity"};

histos.add("hVertexZ_bef", "", kTH1F, {vtxZAxis});
histos.add("hVertexZ_aft", "", kTH1F, {vtxZAxis});
histos.add("hVertexZ_aft_sel", "", kTH1D, {vtxZAxis});
histos.add("hDCAxy_bef", "", kTH1D, {dcaAxis});
histos.add("hDCAxy_aft", "", kTH1D, {dcaAxis});
histos.add("hDCAz_bef", "", kTH1D, {dcazAxis});
histos.add("hDCAz_aft", "", kTH1D, {dcazAxis});
histos.add("hCentrality", "", kTH1D, {centAxis});
histos.add("hMultiplicity", "", kTH1D, {multAxis});
histos.add("hEta", "", kTH1F, {etaAxis});
histos.add("hPt", "", kTH1F, {ptAxis});
AxisSpec vtxZAxis = {100, -20, 20, "Z (cm)"};
AxisSpec dcaAxis = {1000, -100, 100, "DCA_{xy} (cm)"};
AxisSpec dcazAxis = {1000, -100, 100, "DCA_{z} (cm)"};
AxisSpec ptAxis = {40, 0.0, 4.0, "#it{p}_{T} (GeV/#it{c})"};
AxisSpec etaAxis = {30, -1.5, 1.5, "#eta"};
AxisSpec centAxis = {100, 0., 100., "centrality"};
AxisSpec multAxis = {2000, 0., 2000., "multiplicity"};

histos.add("hVertexZ_bef", "", kTH1F, {vtxZAxis});
histos.add("hVertexZ_aft", "", kTH1F, {vtxZAxis});
histos.add("hVertexZ_aft_sel", "", kTH1D, {vtxZAxis});
histos.add("hDCAxy_bef", "", kTH1D, {dcaAxis});
histos.add("hDCAxy_aft", "", kTH1D, {dcaAxis});
histos.add("hDCAz_bef", "", kTH1D, {dcazAxis});
histos.add("hDCAz_aft", "", kTH1D, {dcazAxis});
histos.add("hCentrality", "", kTH1D, {centAxis});
histos.add("hMultiplicity", "", kTH1D, {multAxis});
histos.add("hEta", "", kTH1F, {etaAxis});
histos.add("hPt", "", kTH1F, {ptAxis});
>>>>>>> 03c805b2635be6256d2a2a4ad4e9fa12fbdd067d
}

}

void process(aod::MyCollision const& coll, aod::MyTracks const& inputTracks)
void process(aod::MyCollision const& coll, aod::MyTracks const& inputTracks)
{
<<<<<<< HEAD
histos.fill(HIST("hVertexZ_bef"), coll.posZ());
histos.fill(HIST("hVertexZ_bef"), coll.posZ());

if (fabs(coll.posZ()) > 10.f) {return;}
histos.fill(HIST("hVertexZ_aft"), coll.posZ());
if (fabs(coll.posZ()) > 10.f) {
return;
}
histos.fill(HIST("hVertexZ_aft"), coll.posZ());

if (!coll.sel7()) { return;}
if (!coll.sel7()) {
return;
}
histos.fill(HIST("hVertexZ_aft_sel"), coll.posZ());
histos.fill(HIST("hCentrality"), coll.centRun2V0M());
histos.fill(HIST("hMultiplicity"), coll.multFV0M());


for (auto track : inputTracks)
{
if (fabs(track.eta()) > 0.8) continue;
if (!(track.pt() > 0.2 && track.pt() < 2.)) continue;

histos.fill(HIST("hDCAxy_bef"), track.dcaXY());
histos.fill(HIST("hDCAz_bef"), track.dcaZ());

if (!track.isGlobalTrack()) continue;

histos.fill(HIST("hDCAxy_aft"), track.dcaXY());
histos.fill(HIST("hDCAz_aft"), track.dcaZ());
histos.fill(HIST("hPt"), track.pt());
histos.fill(HIST("hEta"), track.eta());
}
histos.fill(HIST("hCentrality"), coll.centRun2V0M());
histos.fill(HIST("hMultiplicity"), coll.multFV0M());

for (auto track : inputTracks) {
if (fabs(track.eta()) > 0.8)
continue;
if (!(track.pt() > 0.2 && track.pt() < 2.))
continue;

histos.fill(HIST("hDCAxy_bef"), track.dcaXY());
histos.fill(HIST("hDCAz_bef"), track.dcaZ());

if (!track.isGlobalTrack())
continue;

histos.fill(HIST("hDCAxy_aft"), track.dcaXY());
histos.fill(HIST("hDCAz_aft"), track.dcaZ());
histos.fill(HIST("hPt"), track.pt());
histos.fill(HIST("hEta"), track.eta());
}
}
};
WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
{
WorkflowSpec workflow{adaptAnalysisTask<NetchargeFluctuations>(cfgc)};
return workflow;
return workflow;
=======
histos.fill(HIST("hVertexZ_bef"), coll.posZ());
histos.fill(HIST("hVertexZ_bef"), coll.posZ());

if (fabs(coll.posZ()) > 10.f) {return;}
histos.fill(HIST("hVertexZ_aft"), coll.posZ());
if (fabs(coll.posZ()) > 10.f) {
return;
}
histos.fill(HIST("hVertexZ_aft"), coll.posZ());

if (!coll.sel7()) { return;}
if (!coll.sel7()) {
return;
}
histos.fill(HIST("hVertexZ_aft_sel"), coll.posZ());
histos.fill(HIST("hCentrality"), coll.centRun2V0M());
histos.fill(HIST("hMultiplicity"), coll.multFV0M());


for (auto track : inputTracks)
{
if (fabs(track.eta()) > 0.8) continue;
if (!(track.pt() > 0.2 && track.pt() < 2.)) continue;

histos.fill(HIST("hDCAxy_bef"), track.dcaXY());
histos.fill(HIST("hDCAz_bef"), track.dcaZ());

if (!track.isGlobalTrack()) continue;

histos.fill(HIST("hDCAxy_aft"), track.dcaXY());
histos.fill(HIST("hDCAz_aft"), track.dcaZ());
histos.fill(HIST("hPt"), track.pt());
histos.fill(HIST("hEta"), track.eta());
}
histos.fill(HIST("hCentrality"), coll.centRun2V0M());
histos.fill(HIST("hMultiplicity"), coll.multFV0M());

for (auto track : inputTracks) {
if (fabs(track.eta()) > 0.8)
continue;
if (!(track.pt() > 0.2 && track.pt() < 2.))
continue;

histos.fill(HIST("hDCAxy_bef"), track.dcaXY());
histos.fill(HIST("hDCAz_bef"), track.dcaZ());

if (!track.isGlobalTrack())
continue;

histos.fill(HIST("hDCAxy_aft"), track.dcaXY());
histos.fill(HIST("hDCAz_aft"), track.dcaZ());
histos.fill(HIST("hPt"), track.pt());
histos.fill(HIST("hEta"), track.eta());
}
}
};
WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
{
WorkflowSpec workflow{adaptAnalysisTask<NetchargeFluctuations>(cfgc)};
return workflow;
return workflow;
>>>>>>> 03c805b2635be6256d2a2a4ad4e9fa12fbdd067d
}

0 comments on commit e15b3fa

Please sign in to comment.