Skip to content

Commit c83005d

Browse files
dsekihatdsekihat
authored
[PWGEM/Dilepton] upadte dalitzEEQC (#10264)
1 parent 0fcf71d commit c83005d

File tree

2 files changed

+8
-36
lines changed

2 files changed

+8
-36
lines changed

‎PWGEM/PhotonMeson/Tasks/dalitzEEQC.cxx‎

Lines changed: 4 additions & 18 deletions
Original file line numberDiff line numberDiff line change
@@ -82,7 +82,6 @@ struct DalitzEEQC {
8282
Configurable<float> cfg_min_mass{"cfg_min_mass", 0.0, "min mass"};
8383
Configurable<float> cfg_max_mass{"cfg_max_mass", 0.5, "max mass"};
8484
Configurable<bool> cfg_apply_phiv{"cfg_apply_phiv", true, "flag to apply phiv cut"};
85-
Configurable<bool> cfg_apply_pf{"cfg_apply_pf", false, "flag to apply phiv prefilter"};
8685
Configurable<bool> cfg_require_itsib_any{"cfg_require_itsib_any", true, "flag to require ITS ib any hits"};
8786
Configurable<bool> cfg_require_itsib_1st{"cfg_require_itsib_1st", false, "flag to require ITS ib 1st hit"};
8887
Configurable<float> cfg_phiv_slope{"cfg_phiv_slope", 0.0185, "slope for m vs. phiv"};
@@ -175,21 +174,8 @@ struct DalitzEEQC {
175174
o2::aod::pwgem::photonmeson::utils::eventhistogram::addEventHistograms(&fRegistry);
176175

177176
// pair info
178-
std::vector<double> ptbins;
179-
std::vector<double> massbins;
180-
181-
for (int i = 0; i < 51; i++) {
182-
massbins.emplace_back(0.01 * (i - 0) + 0.0); // every 0.01 GeV/c2 from 0.0 to 0.5 GeV/c2
183-
}
184-
const AxisSpec axis_mass{massbins, "m_{ee} (GeV/c^{2})"};
185-
186-
for (int i = 0; i < 50; i++) {
187-
ptbins.emplace_back(0.1 * (i - 0) + 0.0); // every 0.1 GeV/c from 0.0 to 5.0 GeV/c
188-
}
189-
for (int i = 50; i < 61; i++) {
190-
ptbins.emplace_back(0.5 * (i - 50) + 5.0); // every 0.5 GeV/c from 5.0 to 10 GeV/c
191-
}
192-
const AxisSpec axis_pt{ptbins, "p_{T,ee} (GeV/c)"};
177+
const AxisSpec axis_mass{200, 0, 0.2, "m_{ee} (GeV/c^{2})"};
178+
const AxisSpec axis_pt{200, 0, 2, "p_{T,ee} (GeV/c)"};
193179

194180
fRegistry.add("Pair/same/hMvsPt", "m_{ee} vs. p_{T,ee};m_{ee} (GeV/c^{2});p_{T,ee} (GeV/c)", kTH2F, {axis_mass, axis_pt}, true);
195181
fRegistry.add("Pair/same/hMvsPhiV", "m_{ee} vs. #varphi_{V};#varphi (rad.);m_{ee} (GeV/c^{2})", kTH2F, {{90, 0, M_PI}, {100, 0.0f, 0.1f}}, true);
@@ -331,8 +317,8 @@ struct DalitzEEQC {
331317

332318
SliceCache cache;
333319
Preslice<MyTracks> perCollision_track = aod::emprimaryelectron::emeventId;
334-
Filter trackFilter = static_cast<float>(dileptoncuts.cfg_min_pt_track) < o2::aod::track::pt && nabs(o2::aod::track::eta) < static_cast<float>(dileptoncuts.cfg_max_eta_track) && o2::aod::track::tpcChi2NCl < static_cast<float>(dileptoncuts.cfg_max_chi2tpc) && o2::aod::track::itsChi2NCl < static_cast<float>(dileptoncuts.cfg_max_chi2its) && nabs(o2::aod::track::dcaXY) < static_cast<float>(dileptoncuts.cfg_max_dcaxy) && nabs(o2::aod::track::dcaZ) < static_cast<float>(dileptoncuts.cfg_max_dcaz);
335-
Filter pidFilter = (static_cast<float>(dileptoncuts.cfg_min_TPCNsigmaEl) < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < static_cast<float>(dileptoncuts.cfg_max_TPCNsigmaEl)) && (o2::aod::pidtpc::tpcNSigmaPi < static_cast<float>(dileptoncuts.cfg_min_TPCNsigmaPi) || static_cast<float>(dileptoncuts.cfg_max_TPCNsigmaPi) < o2::aod::pidtpc::tpcNSigmaPi);
320+
Filter trackFilter = dileptoncuts.cfg_min_pt_track < o2::aod::track::pt && nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track && o2::aod::track::tpcChi2NCl < dileptoncuts.cfg_max_chi2tpc && o2::aod::track::itsChi2NCl < dileptoncuts.cfg_max_chi2its && nabs(o2::aod::track::dcaXY) < dileptoncuts.cfg_max_dcaxy && nabs(o2::aod::track::dcaZ) < dileptoncuts.cfg_max_dcaz;
321+
Filter pidFilter = dileptoncuts.cfg_min_TPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < dileptoncuts.cfg_max_TPCNsigmaEl && (o2::aod::pidtpc::tpcNSigmaPi < dileptoncuts.cfg_min_TPCNsigmaPi || dileptoncuts.cfg_max_TPCNsigmaPi < o2::aod::pidtpc::tpcNSigmaPi);
336322
using FilteredMyTracks = soa::Filtered<MyTracks>;
337323
Partition<FilteredMyTracks> posTracks = o2::aod::emprimaryelectron::sign > int8_t(0);
338324
Partition<FilteredMyTracks> negTracks = o2::aod::emprimaryelectron::sign < int8_t(0);

‎PWGEM/PhotonMeson/Tasks/dalitzEEQCMC.cxx‎

Lines changed: 4 additions & 18 deletions
Original file line numberDiff line numberDiff line change
@@ -86,7 +86,6 @@ struct DalitzEEQCMC {
8686
Configurable<float> cfg_min_mass{"cfg_min_mass", 0.0, "min mass"};
8787
Configurable<float> cfg_max_mass{"cfg_max_mass", 1e+10, "max mass"};
8888
Configurable<bool> cfg_apply_phiv{"cfg_apply_phiv", true, "flag to apply phiv cut"};
89-
Configurable<bool> cfg_apply_pf{"cfg_apply_pf", false, "flag to apply phiv prefilter"};
9089
Configurable<bool> cfg_require_itsib_any{"cfg_require_itsib_any", true, "flag to require ITS ib any hits"};
9190
Configurable<bool> cfg_require_itsib_1st{"cfg_require_itsib_1st", false, "flag to require ITS ib 1st hit"};
9291
Configurable<float> cfg_phiv_slope{"cfg_phiv_slope", 0.0185, "slope for m vs. phiv"};
@@ -131,21 +130,8 @@ struct DalitzEEQCMC {
131130
// event info
132131
o2::aod::pwgem::photonmeson::utils::eventhistogram::addEventHistograms(&fRegistry);
133132

134-
std::vector<double> ptbins;
135-
std::vector<double> massbins;
136-
137-
for (int i = 0; i < 51; i++) {
138-
massbins.emplace_back(0.01 * (i - 0) + 0.0); // every 0.01 GeV/c2 from 0.0 to 0.5 GeV/c2
139-
}
140-
const AxisSpec axis_mass{massbins, "m_{ee} (GeV/c^{2})"};
141-
142-
for (int i = 0; i < 50; i++) {
143-
ptbins.emplace_back(0.1 * (i - 0) + 0.0); // every 0.1 GeV/c from 0.0 to 5.0 GeV/c
144-
}
145-
for (int i = 50; i < 61; i++) {
146-
ptbins.emplace_back(0.5 * (i - 50) + 5.0); // every 0.5 GeV/c from 5.0 to 10 GeV/c
147-
}
148-
const AxisSpec axis_pt{ptbins, "p_{T,ee} (GeV/c)"};
133+
const AxisSpec axis_mass{200, 0, 0.2, "m_{ee} (GeV/c^{2})"};
134+
const AxisSpec axis_pt{200, 0, 2, "p_{T,ee} (GeV/c)"};
149135

150136
// generated info
151137
fRegistry.add("Generated/sm/Pi0/hMvsPt", "m_{ee} vs. p_{T,ee} ULS", kTH2F, {axis_mass, axis_pt}, true);
@@ -431,8 +417,8 @@ struct DalitzEEQCMC {
431417
std::vector<int> used_trackIds;
432418
SliceCache cache;
433419
Preslice<MyMCTracks> perCollision_track = aod::emprimaryelectron::emeventId;
434-
Filter trackFilter = static_cast<float>(dileptoncuts.cfg_min_pt_track) < o2::aod::track::pt && nabs(o2::aod::track::eta) < static_cast<float>(dileptoncuts.cfg_max_eta_track) && o2::aod::track::tpcChi2NCl < static_cast<float>(dileptoncuts.cfg_max_chi2tpc) && o2::aod::track::itsChi2NCl < static_cast<float>(dileptoncuts.cfg_max_chi2its) && nabs(o2::aod::track::dcaXY) < static_cast<float>(dileptoncuts.cfg_max_dcaxy) && nabs(o2::aod::track::dcaZ) < static_cast<float>(dileptoncuts.cfg_max_dcaz);
435-
Filter pidFilter = (static_cast<float>(dileptoncuts.cfg_min_TPCNsigmaEl) < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < static_cast<float>(dileptoncuts.cfg_max_TPCNsigmaEl)) && (o2::aod::pidtpc::tpcNSigmaPi < static_cast<float>(dileptoncuts.cfg_min_TPCNsigmaPi) || static_cast<float>(dileptoncuts.cfg_max_TPCNsigmaPi) < o2::aod::pidtpc::tpcNSigmaPi);
420+
Filter trackFilter = dileptoncuts.cfg_min_pt_track < o2::aod::track::pt && nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track && o2::aod::track::tpcChi2NCl < dileptoncuts.cfg_max_chi2tpc && o2::aod::track::itsChi2NCl < dileptoncuts.cfg_max_chi2its && nabs(o2::aod::track::dcaXY) < dileptoncuts.cfg_max_dcaxy && nabs(o2::aod::track::dcaZ) < dileptoncuts.cfg_max_dcaz;
421+
Filter pidFilter = dileptoncuts.cfg_min_TPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < dileptoncuts.cfg_max_TPCNsigmaEl && (o2::aod::pidtpc::tpcNSigmaPi < dileptoncuts.cfg_min_TPCNsigmaPi || dileptoncuts.cfg_max_TPCNsigmaPi < o2::aod::pidtpc::tpcNSigmaPi);
436422
using FilteredMyMCTracks = soa::Filtered<MyMCTracks>;
437423
Partition<FilteredMyMCTracks> posTracks = o2::aod::emprimaryelectron::sign > int8_t(0);
438424
Partition<FilteredMyMCTracks> negTracks = o2::aod::emprimaryelectron::sign < int8_t(0);

0 commit comments

Comments
 (0)