Skip to content
Merged
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
11 changes: 5 additions & 6 deletions EventFiltering/PWGJE/jetFilter.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -46,7 +46,7 @@ static const std::vector<std::string> highPtObjectsNames = {"JetChLowPt", "JetCh
struct jetFilter {
enum { kJetChLowPt = 0,
kJetChHighPt = 1,
kHighPtObjects = 2,
kHighPtObjects = 2,
kAllMBEvents = 3,
kAllObjects };

Expand Down Expand Up @@ -127,7 +127,6 @@ struct jetFilter {
spectra.add("hRho", "Underlying event density #rho", HistType::kTH1F,
{{200, 0., +20., "#rho (GeV/#it{c})"}});


auto scalers{std::get<std::shared_ptr<TH1>>(spectra.add(
"fProcessedEvents", ";;Number of filtered events", HistType::kTH1F,
{{kAllObjects, -0.5, kAllObjects - 0.5}}))};
Expand All @@ -148,20 +147,20 @@ struct jetFilter {
// collision process loop
bool keepEvent[kHighPtObjects]{false};
spectra.fill(HIST("fCollZpos"), collision.posZ());
spectra.fill(HIST("fProcessedEvents"), kAllMBEvents); //all minimum bias events
spectra.fill(HIST("fProcessedEvents"), kAllMBEvents); // all minimum bias events

// FILL SPECTRA OF INCLUSIVE JETS IN FIDUCIAL VOLUME
if (TMath::Abs(collision.posZ()) < cfgZvtx) {
spectra.fill(HIST("fProcessedEvents"), kHighPtObjects); //minimum bias events |z_vtx|<10 cm
spectra.fill(HIST("hRho"), collision.rho());
spectra.fill(HIST("fProcessedEvents"), kHighPtObjects); // minimum bias events |z_vtx|<10 cm
spectra.fill(HIST("hRho"), collision.rho());

for (const auto& jet : jets) { // jets are ordered by pT
for (unsigned int ir = 0; ir < jetIntR.size(); ir++) {
if (jet.r() == jetIntR[ir]) {
if (TMath::Abs(jet.eta()) < jetRFidVolume[ir]) {
float jetr = (jet.r() / 100. + 1e-5);
spectra.fill(HIST("hPtAKTJetsInclusive"), jetr, jet.pt());
spectra.fill(HIST("hPtAKTJetsInclusiveBgSubtr"), jetr, jet.pt() - (collision.rho() * jet.area()));
spectra.fill(HIST("hPtAKTJetsInclusiveBgSubtr"), jetr, jet.pt() - (collision.rho() * jet.area()));

if (jet.pt() > 10.) {
spectra.fill(HIST("hEtaAKTJetsInclusive"), jetr, jet.eta());
Expand Down