Skip to content

Commit eae3f16

Browse files
dsekihatdsekihat
committed
[PWGEM/Dilepton] update task v1
1 parent 520c3f2 commit eae3f16

1 file changed

Lines changed: 23 additions & 7 deletions

File tree

‎PWGEM/Dilepton/Tasks/dimuonV1.cxx‎

Lines changed: 23 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -88,7 +88,7 @@ struct dimuonV1 {
8888

8989
o2::framework::ConfigurableAxis ConfPtllBins{"ConfPtllBins", {10, 0, 10}, "pTll bins for output histograms"};
9090
o2::framework::ConfigurableAxis ConfYllBins{"ConfYllBins", {3, -4.0, -2.5}, "yll bins for output histograms"};
91-
o2::framework::ConfigurableAxis ConfUQBins{"ConfUQBins", {200, -1, 1}, "uQ bins for output histograms"};
91+
o2::framework::ConfigurableAxis ConfUQBins{"ConfUQBins", {400, -1, 1}, "uQ bins for output histograms"};
9292
o2::framework::Configurable<int> cfgNrotation{"cfgNrotation", 1, "number of rotation bkg"};
9393
o2::framework::Configurable<int> cfgRandomSeed{"cfgRandomSeed", 1, "randam seed for rotation bkg"};
9494
o2::framework::Configurable<float> cfgRotationMin{"cfgRotationMin", -M_PI / 4, "min. rotation angle for rotation bkg"};
@@ -147,7 +147,7 @@ struct dimuonV1 {
147147
o2::framework::Configurable<float> cfg_min_deta{"cfg_min_deta", 0.02, "min deta between 2 muons (elliptic cut)"};
148148
o2::framework::Configurable<float> cfg_min_dphi{"cfg_min_dphi", 0.02, "min dphi between 2 muons (elliptic cut)"};
149149

150-
o2::framework::Configurable<uint8_t> cfg_track_type{"cfg_track_type", 3, "muon track type [0: MFT-MCH-MID, 3: MCH-MID]"};
150+
// o2::framework::Configurable<uint8_t> cfg_track_type{"cfg_track_type", 3, "muon track type [0: MFT-MCH-MID, 3: MCH-MID]"};
151151
o2::framework::Configurable<float> cfg_min_pt_track{"cfg_min_pt_track", 0.8, "min pT for single track"};
152152
o2::framework::Configurable<float> cfg_max_pt_track{"cfg_max_pt_track", 1e+10, "max pT for single track"};
153153
o2::framework::Configurable<float> cfg_min_eta_track{"cfg_min_eta_track", -4.0, "min eta for single track"};
@@ -222,10 +222,15 @@ struct dimuonV1 {
222222
// event info
223223
o2::aod::pwgem::dilepton::utils::eventhistogram::addEventHistograms<-1>(&fRegistry);
224224

225-
fRegistry.add("Event/before/ZDC/hQxtQxp", "Q_{x}^{t} #upoint Q_{x}^{p} vs. centrality;centrality FT0C (%);Q_{x}^{t} #upoint Q_{x}^{p}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {2000, -1, +1}}, false);
226-
fRegistry.add("Event/before/ZDC/hQytQyp", "Q_{y}^{t} #upoint Q_{y}^{p} vs. centrality;centrality FT0C (%);Q_{y}^{t} #upoint Q_{y}^{p}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {2000, -1, +1}}, false);
227-
fRegistry.add("Event/before/ZDC/hQxtQyp", "Q_{x}^{t} #upoint Q_{y}^{p} vs. centrality;centrality FT0C (%);Q_{x}^{t} #upoint Q_{y}^{p}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {2000, -1, +1}}, false);
228-
fRegistry.add("Event/before/ZDC/hQytQxp", "Q_{y}^{t} #upoint Q_{x}^{p} vs. centrality;centrality FT0C (%);Q_{y}^{t} #upoint Q_{x}^{p}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {2000, -1, +1}}, false);
225+
fRegistry.add("Event/before/ZDC/hQxt", "Q_{x}^{t} vs. centrality;centrality FT0C (%);Q_{x}^{t}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {200, -1, +1}}, false);
226+
fRegistry.add("Event/before/ZDC/hQyt", "Q_{y}^{t} vs. centrality;centrality FT0C (%);Q_{y}^{t}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {200, -1, +1}}, false);
227+
fRegistry.add("Event/before/ZDC/hQxp", "Q_{x}^{p} vs. centrality;centrality FT0C (%);Q_{x}^{p}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {200, -1, +1}}, false);
228+
fRegistry.add("Event/before/ZDC/hQyp", "Q_{y}^{p} vs. centrality;centrality FT0C (%);Q_{y}^{p}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {200, -1, +1}}, false);
229+
230+
fRegistry.add("Event/before/ZDC/hQxtQxp", "Q_{x}^{t} #upoint Q_{x}^{p} vs. centrality;centrality FT0C (%);Q_{x}^{t} #upoint Q_{x}^{p}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {1000, -1, +1}}, false);
231+
fRegistry.add("Event/before/ZDC/hQytQyp", "Q_{y}^{t} #upoint Q_{y}^{p} vs. centrality;centrality FT0C (%);Q_{y}^{t} #upoint Q_{y}^{p}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {1000, -1, +1}}, false);
232+
fRegistry.add("Event/before/ZDC/hQxtQyp", "Q_{x}^{t} #upoint Q_{y}^{p} vs. centrality;centrality FT0C (%);Q_{x}^{t} #upoint Q_{y}^{p}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {1000, -1, +1}}, false);
233+
fRegistry.add("Event/before/ZDC/hQytQxp", "Q_{y}^{t} #upoint Q_{x}^{p} vs. centrality;centrality FT0C (%);Q_{y}^{t} #upoint Q_{x}^{p}", o2::framework::HistType::kTH2D, {{110, 0, 110}, {1000, -1, +1}}, false);
229234
fRegistry.addClone("Event/before/ZDC/", "Event/after/ZDC/");
230235

231236
// pair info
@@ -279,7 +284,8 @@ struct dimuonV1 {
279284
fDimuonCut.SetMindEtadPhi(dimuoncuts.cfg_apply_detadphi, dimuoncuts.cfg_min_deta, dimuoncuts.cfg_min_dphi);
280285

281286
// for track
282-
fDimuonCut.SetTrackType(dimuoncuts.cfg_track_type);
287+
// fDimuonCut.SetTrackType(dimuoncuts.cfg_track_type);
288+
fDimuonCut.SetTrackType(3);
283289
fDimuonCut.SetTrackPtRange(dimuoncuts.cfg_min_pt_track, dimuoncuts.cfg_max_pt_track);
284290
fDimuonCut.SetTrackEtaRange(dimuoncuts.cfg_min_eta_track, dimuoncuts.cfg_max_eta_track);
285291
fDimuonCut.SetTrackPhiRange(dimuoncuts.cfg_min_phi_track, dimuoncuts.cfg_max_phi_track);
@@ -362,6 +368,11 @@ struct dimuonV1 {
362368
float QytQxp = collision.qxZDCA() * collision.qyZDCC();
363369
float QxtQyp = collision.qxZDCC() * collision.qyZDCA();
364370

371+
fRegistry.fill(HIST("Event/before/ZDC/hQxt"), centrality, collision.qxZDCC());
372+
fRegistry.fill(HIST("Event/before/ZDC/hQyt"), centrality, collision.qyZDCC());
373+
fRegistry.fill(HIST("Event/before/ZDC/hQxp"), centrality, collision.qxZDCA());
374+
fRegistry.fill(HIST("Event/before/ZDC/hQyp"), centrality, collision.qyZDCA());
375+
365376
fRegistry.fill(HIST("Event/before/ZDC/hQxtQxp"), centrality, QxtQxp);
366377
fRegistry.fill(HIST("Event/before/ZDC/hQytQyp"), centrality, QytQyp);
367378
fRegistry.fill(HIST("Event/before/ZDC/hQxtQyp"), centrality, QxtQyp);
@@ -379,6 +390,11 @@ struct dimuonV1 {
379390
fRegistry.fill(HIST("Event/before/hCollisionCounter"), o2::aod::pwgem::dilepton::utils::eventhistogram::nbin_ev); // accepted
380391
fRegistry.fill(HIST("Event/after/hCollisionCounter"), o2::aod::pwgem::dilepton::utils::eventhistogram::nbin_ev); // accepted
381392

393+
fRegistry.fill(HIST("Event/after/ZDC/hQxt"), centrality, collision.qxZDCC());
394+
fRegistry.fill(HIST("Event/after/ZDC/hQyt"), centrality, collision.qyZDCC());
395+
fRegistry.fill(HIST("Event/after/ZDC/hQxp"), centrality, collision.qxZDCA());
396+
fRegistry.fill(HIST("Event/after/ZDC/hQyp"), centrality, collision.qyZDCA());
397+
382398
fRegistry.fill(HIST("Event/after/ZDC/hQxtQxp"), centrality, QxtQxp);
383399
fRegistry.fill(HIST("Event/after/ZDC/hQytQyp"), centrality, QytQyp);
384400
fRegistry.fill(HIST("Event/after/ZDC/hQxtQyp"), centrality, QxtQyp);

0 commit comments

Comments
 (0)