Skip to content

Commit 449fa6f

Browse files
dsekihatdsekihat
authored
PWGEM/PhotonMeson: update for TPConly tracks (#4097)
1 parent 937f1c7 commit 449fa6f

9 files changed

Lines changed: 132 additions & 68 deletions

File tree

‎PWGEM/PhotonMeson/Core/CutsLibrary.cxx‎

Lines changed: 11 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -48,8 +48,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
4848
cut->SetV0PtRange(0.1f, 1e10f);
4949
cut->SetV0EtaRange(-0.9, +0.9);
5050
cut->SetMinCosPA(0.99);
51-
cut->SetMaxPCA(0.5);
52-
cut->SetRxyRange(1, 180);
51+
cut->SetMaxPCA(1.5);
52+
cut->SetRxyRange(1, 90);
5353
return cut;
5454
}
5555
if (!nameStr.compare("qc_ITSTPC")) {
@@ -67,8 +67,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
6767
cut->SetV0EtaRange(-0.9, +0.9);
6868
cut->SetMinCosPA(0.99);
6969
cut->SetMaxPCA(0.5);
70-
cut->SetRxyRange(1, 180);
71-
cut->SetAPRange(0.95, 0.02);
70+
cut->SetRxyRange(1, 90);
71+
cut->SetAPRange(0.95, 0.01);
7272
return cut;
7373
}
7474
if (!nameStr.compare("qc_ITSonly")) {
@@ -82,8 +82,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
8282
cut->SetV0EtaRange(-0.9, +0.9);
8383
cut->SetMinCosPA(0.99);
8484
cut->SetMaxPCA(0.5);
85-
cut->SetRxyRange(1, 180);
86-
cut->SetAPRange(0.95, 0.02);
85+
cut->SetRxyRange(1, 90);
86+
cut->SetAPRange(0.95, 0.01);
8787
return cut;
8888
}
8989
if (!nameStr.compare("qc_TPConly")) {
@@ -101,7 +101,8 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
101101
cut->SetV0EtaRange(-0.9, +0.9);
102102
cut->SetMinCosPA(0.99);
103103
cut->SetMaxPCA(1.5);
104-
cut->SetRxyRange(1, 180);
104+
cut->SetRxyRange(1, 90);
105+
cut->SetAPRange(0.95, 0.01);
105106
return cut;
106107
}
107108
if (!nameStr.compare("wwire")) { // conversion only on tungstate wire
@@ -121,7 +122,7 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
121122
cut->SetMaxPCA(0.5);
122123
cut->SetOnWwireIB(true);
123124
cut->SetOnWwireOB(true);
124-
cut->SetAPRange(0.95, 0.02);
125+
cut->SetAPRange(0.95, 0.01);
125126
return cut;
126127
}
127128
if (!nameStr.compare("wwire_ib")) { // conversion only on tungstate wire outside of ITSib
@@ -141,7 +142,7 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
141142
cut->SetMaxPCA(0.5);
142143
cut->SetOnWwireIB(true);
143144
cut->SetOnWwireOB(false);
144-
cut->SetAPRange(0.95, 0.02);
145+
cut->SetAPRange(0.95, 0.01);
145146
return cut;
146147
}
147148
if (!nameStr.compare("wwire_ob")) { // conversion only on tungstate wire outside of ITSob (middle layer)
@@ -161,7 +162,7 @@ V0PhotonCut* o2::aod::pcmcuts::GetCut(const char* cutName)
161162
cut->SetMaxPCA(0.5);
162163
cut->SetOnWwireIB(false);
163164
cut->SetOnWwireOB(true);
164-
cut->SetAPRange(0.95, 0.02);
165+
cut->SetAPRange(0.95, 0.01);
165166
return cut;
166167
}
167168
if (!nameStr.compare("nopid")) {

‎PWGEM/PhotonMeson/Core/HistogramsLibrary.cxx‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -64,7 +64,7 @@ void o2::aod::emphotonhistograms::DefineHistograms(THashList* list, const char*
6464
list->Add(new TH2F("hZX", "Z vs. X;Z (cm);X (cm)", 200, -100, 100, 100, 0, 100));
6565
list->Add(new TH2F("hZY", "Z vs. Y;Z (cm);Y (cm)", 200, -100, 100, 100, -50, 50));
6666
list->Add(new TH2F("hDCAxyZ", "DCAxy vs. Z;Z (cm);DCA_{xy} (cm)", 200, -100, +100, 100, -50, 50));
67-
list->Add(new TH2F("hXZ_tgl", "correlation of tgl vs. z/x;tgl;z/x", 300, -1.5, 1.5, 300, -1.5, 1.5));
67+
list->Add(new TH2F("hXZ_tgl", "correlation of tgl vs. z/x;tgl;z/x - tgl", 300, -1.5, 1.5, 300, -1.5, 1.5));
6868
if (TString(subGroup) == "mc") {
6969
list->Add(new TH2F("hPtGen_DeltaPtOverPtGen", "electron p_{T} resolution;p_{T}^{gen} (GeV/c);(p_{T}^{rec} - p_{T}^{gen})/p_{T}^{gen}", 1000, 0, 10, 400, -1.0f, 1.0f));
7070
list->Add(new TH2F("hPtGen_DeltaEta", "electron #eta resolution;p_{T}^{gen} (GeV/c);#eta^{rec} - #eta^{gen}", 1000, 0, 10, 400, -1.0f, 1.0f));

‎PWGEM/PhotonMeson/Core/V0PhotonCut.h‎

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -386,17 +386,17 @@ class V0PhotonCut : public TNamed
386386
return mMinChi2PerClusterITS < track.itsChi2NCl() && track.itsChi2NCl() < mMaxChi2PerClusterITS;
387387

388388
case V0PhotonCuts::kIsWithinBeamPipe: {
389-
// return track.isWithinBeamPipe();
390-
if (abs(track.y()) > abs(track.x() * TMath::Tan(10.f * TMath::DegToRad())) + 15.f) {
391-
return false;
392-
}
389+
// if (abs(track.y()) > abs(track.x() * TMath::Tan(10.f * TMath::DegToRad())) + 15.f) {
390+
// return false;
391+
// }
393392

394-
if (track.x() > 82.9 && abs(track.y()) > abs(track.x() * TMath::Tan(10.f * TMath::DegToRad())) + 5.f) {
395-
return false;
396-
}
393+
// if (track.x() > 82.9 && abs(track.y()) > abs(track.x() * TMath::Tan(10.f * TMath::DegToRad())) + 5.f) {
394+
// return false;
395+
// }
397396

398397
const float slope = TMath::Tan(2 * TMath::ATan(TMath::Exp(-0.5)));
399-
return !(track.x() > 82.9 && abs(track.y()) < 40.f && abs(abs(track.z()) - track.x() / slope) < 3.5f && 15.f < abs(track.dcaXY()));
398+
return !(track.x() > 82.9 && abs(track.y()) < 15.f && abs(abs(track.z()) - track.x() / slope) < 7.f && 15.f < abs(track.dcaXY()));
399+
// return !(track.x() > 82.9 && abs(track.y()) < 40.f && abs(abs(track.z()) - 47.f) < 3.f && 15.f < abs(track.dcaXY()));
400400
}
401401
case V0PhotonCuts::kRequireITSTPC:
402402
return isITSTPCTrack(track);

‎PWGEM/PhotonMeson/DataModel/gammaTables.h‎

Lines changed: 10 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -269,14 +269,14 @@ DECLARE_SOA_COLUMN(Px, px, float); //! px f
269269
DECLARE_SOA_COLUMN(Py, py, float); //! py for photon kf
270270
DECLARE_SOA_COLUMN(Pz, pz, float); //! pz for photon kf
271271
DECLARE_SOA_COLUMN(MGamma, mGamma, float); //! invariant mass of dielectron at SV
272-
DECLARE_SOA_COLUMN(MGammaPV, mGammaPV, float); //! invariant mass of dielectron at PV
273-
DECLARE_SOA_COLUMN(DCAxyV0ToPV, dcaXYv0topv, float); //! DCAxy of V0 to PV
274-
DECLARE_SOA_COLUMN(DCAzV0ToPV, dcaZv0topv, float); //! DCAz of V0 to PV
275-
DECLARE_SOA_COLUMN(CosPA, cospa, float); //!
276-
DECLARE_SOA_COLUMN(PCA, pca, float); //!
277-
DECLARE_SOA_COLUMN(Alpha, alpha, float); //!
278-
DECLARE_SOA_COLUMN(QtArm, qtarm, float); //!
279-
DECLARE_SOA_COLUMN(ChiSquareNDF, chiSquareNDF, float); // Chi2 / NDF of the reconstructed V0
272+
// DECLARE_SOA_COLUMN(MGammaPV, mGammaPV, float); //! invariant mass of dielectron at PV
273+
DECLARE_SOA_COLUMN(DCAxyV0ToPV, dcaXYv0topv, float); //! DCAxy of V0 to PV
274+
DECLARE_SOA_COLUMN(DCAzV0ToPV, dcaZv0topv, float); //! DCAz of V0 to PV
275+
DECLARE_SOA_COLUMN(CosPA, cospa, float); //!
276+
DECLARE_SOA_COLUMN(PCA, pca, float); //!
277+
DECLARE_SOA_COLUMN(Alpha, alpha, float); //!
278+
DECLARE_SOA_COLUMN(QtArm, qtarm, float); //!
279+
DECLARE_SOA_COLUMN(ChiSquareNDF, chiSquareNDF, float); // Chi2 / NDF of the reconstructed V0
280280

281281
DECLARE_SOA_DYNAMIC_COLUMN(E, e, [](float px, float py, float pz, float m = 0) -> float { return RecoDecay::sqrtSumOfSquares(px, py, pz, m); }); //! energy of v0 photn, mass to be given as argument when getter is called!
282282
DECLARE_SOA_DYNAMIC_COLUMN(Pt, pt, [](float px, float py) -> float { return RecoDecay::sqrtSumOfSquares(px, py); });
@@ -289,7 +289,8 @@ DECLARE_SOA_TABLE(V0PhotonsKF, "AOD", "V0PHOTONKF", //!
289289
o2::soa::Index<>, v0photonkf::CollisionId, v0photonkf::PosTrackId, v0photonkf::NegTrackId,
290290
v0photonkf::Vx, v0photonkf::Vy, v0photonkf::Vz,
291291
v0photonkf::Px, v0photonkf::Py, v0photonkf::Pz,
292-
v0photonkf::MGamma, v0photonkf::MGammaPV,
292+
v0photonkf::MGamma,
293+
// v0photonkf::MGammaPV,
293294
v0photonkf::DCAxyV0ToPV, v0photonkf::DCAzV0ToPV,
294295
v0photonkf::CosPA, v0photonkf::PCA,
295296
v0photonkf::Alpha, v0photonkf::QtArm,

0 commit comments

Comments
 (0)