Skip to content

Commit 4397f4c

Browse files
Archita-DashArchita-Dash
committed
added more event selection bins before and after VtxZ cut to estimate the vertex reconstruction bias and verter reconstruction efficiency
1 parent bebd0d8 commit 4397f4c

3 files changed

Lines changed: 97 additions & 13 deletions

File tree

‎PWGJE/DataModel/JetReducedData.h‎

Lines changed: 12 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -98,8 +98,14 @@ DECLARE_SOA_COLUMN(FullTriggerSel, fullTriggerSel, uint32_t);
9898
DECLARE_SOA_COLUMN(ChargedHFTriggerSel, chargedHFTriggerSel, uint8_t);
9999
DECLARE_SOA_COLUMN(ReadCounts, readCounts, std::vector<int>);
100100
DECLARE_SOA_COLUMN(ReadCountsWithTVX, readCountsWithTVX, std::vector<int>);
101+
DECLARE_SOA_COLUMN(ReadCountsWithTVXAndSel8, readCountsWithTVXAndSel8, std::vector<int>);
102+
DECLARE_SOA_COLUMN(ReadCountsWithTVXAndSel8AndIsGoodZvtxFT0vsPV, readCountsWithTVXAndSel8AndIsGoodZvtxFT0vsPV, std::vector<int>);
103+
DECLARE_SOA_COLUMN(ReadCountsWithTVXAndSel8Full, readCountsWithTVXAndSel8Full, std::vector<int>);
104+
DECLARE_SOA_COLUMN(ReadCountsWithTVXAndSel8FullAndIsGoodZvtxFT0vsPV, readCountsWithTVXAndSel8FullAndIsGoodZvtxFT0vsPV, std::vector<int>);
101105
DECLARE_SOA_COLUMN(ReadCountsWithTVXAndZVertexAndSel8, readCountsWithTVXAndZVertexAndSel8, std::vector<int>);
106+
DECLARE_SOA_COLUMN(ReadCountsWithTVXAndZVertexAndSel8AndIsGoodZvtxFT0vsPV, readCountsWithTVXAndZVertexAndSel8AndIsGoodZvtxFT0vsPV, std::vector<int>);
102107
DECLARE_SOA_COLUMN(ReadCountsWithTVXAndZVertexAndSel8Full, readCountsWithTVXAndZVertexAndSel8Full, std::vector<int>);
108+
DECLARE_SOA_COLUMN(ReadCountsWithTVXAndZVertexAndSel8FullAndIsGoodZvtxFT0vsPV, readCountsWithTVXAndZVertexAndSel8FullAndIsGoodZvtxFT0vsPV, std::vector<int>);
103109
DECLARE_SOA_COLUMN(ReadCountsWithTVXAndZVertexAndSel8FullPbPb, readCountsWithTVXAndZVertexAndSel8FullPbPb, std::vector<int>);
104110
DECLARE_SOA_COLUMN(ReadCountsWithTVXAndZVertexAndSelMC, readCountsWithTVXAndZVertexAndSelMC, std::vector<int>);
105111
DECLARE_SOA_COLUMN(ReadCountsWithTVXAndZVertexAndSelMCFull, readCountsWithTVXAndZVertexAndSelMCFull, std::vector<int>);
@@ -171,8 +177,14 @@ DECLARE_SOA_TABLE(JChHFTrigSels, "AOD", "JCHHFTRIGSEL",
171177
DECLARE_SOA_TABLE_STAGED(CollisionCounts, "COLLCOUNT",
172178
jcollision::ReadCounts,
173179
jcollision::ReadCountsWithTVX,
180+
jcollision::ReadCountsWithTVXAndSel8,
181+
jcollision::ReadCountsWithTVXAndSel8AndIsGoodZvtxFT0vsPV,
182+
jcollision::ReadCountsWithTVXAndSel8Full,
183+
jcollision::ReadCountsWithTVXAndSel8FullAndIsGoodZvtxFT0vsPV,
174184
jcollision::ReadCountsWithTVXAndZVertexAndSel8,
185+
jcollision::ReadCountsWithTVXAndZVertexAndSel8AndIsGoodZvtxFT0vsPV,
175186
jcollision::ReadCountsWithTVXAndZVertexAndSel8Full,
187+
jcollision::ReadCountsWithTVXAndZVertexAndSel8FullAndIsGoodZvtxFT0vsPV,
176188
jcollision::ReadCountsWithTVXAndZVertexAndSel8FullPbPb,
177189
jcollision::ReadCountsWithTVXAndZVertexAndSelMC,
178190
jcollision::ReadCountsWithTVXAndZVertexAndSelMCFull,

‎PWGJE/TableProducer/luminosityCalculator.cxx‎

Lines changed: 30 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -40,7 +40,7 @@ struct LuminosityCalculator {
4040
void init(InitContext&)
4141
{
4242

43-
std::vector<std::string> histLabels = {"BC", "BC+TVX", "BC+TVX+NoTFB", "BC+TVX+NoTFB+NoITSROFB", "Coll", "Coll+TVX", "Coll+TVX+VtxZ+Sel8", "Coll+TVX+VtxZ+Sel8Full", "Coll+TVX+VtxZ+Sel8FullPbPb", "Coll+TVX+VtxZ+SelMC", "Coll+TVX+VtxZ+SelMCFull", "Coll+TVX+VtxZ+SelMCFullPbPb", "Coll+TVX+VtxZ+SelUnanchoredMC", "Coll+TVX+VtxZ+SelTVX", "Coll+TVX+VtxZ+Sel7", "Coll+TVX+VtxZ+Sel7KINT7", "custom"};
43+
std::vector<std::string> histLabels = {"BC", "BC+TVX", "BC+TVX+NoTFB", "BC+TVX+NoTFB+NoITSROFB", "Coll", "Coll+TVX", "Coll+TVX+Sel8", "Coll+TVX+Sel8+IsGoodZvtx", "Coll+TVX+Sel8Full", "Coll+TVX+VtxZ+Sel8", "Coll+TVX+VtxZ+Sel8+IsGoodZvtx", "Coll+TVX+VtxZ+Sel8Full", "Coll+TVX+VtxZ+Sel8Full+IsGoodZvtx", "Coll+TVX+VtxZ+Sel8FullPbPb", "Coll+TVX+VtxZ+SelMC", "Coll+TVX+VtxZ+SelMCFull", "Coll+TVX+VtxZ+SelMCFullPbPb", "Coll+TVX+VtxZ+SelUnanchoredMC", "Coll+TVX+VtxZ+SelTVX", "Coll+TVX+VtxZ+Sel7", "Coll+TVX+VtxZ+Sel7KINT7", "custom"};
4444
registry.add("counter", "BCs and Collisions", HistType::kTH1D, {{static_cast<int>(histLabels.size()), -0.5, static_cast<double>(histLabels.size()) - 0.5}});
4545
auto counter = registry.get<TH1>(HIST("counter"));
4646
for (std::vector<std::string>::size_type iCounter = 0; iCounter < histLabels.size(); iCounter++) {
@@ -64,8 +64,14 @@ struct LuminosityCalculator {
6464

6565
int readCollision = 0;
6666
int readCollisionWithTVXCounter = 0;
67+
int readCollisionWithTVXAndSel8Counter = 0;
68+
int readCollisionWithTVXAndSel8AndIsGoodZvtxFT0vsPVCounter = 0;
69+
int readCollisionWithTVXAndSel8FullCounter = 0;
70+
int readCollisionWithTVXAndSel8FullAndIsGoodZvtxFT0vsPVCounter = 0;
6771
int readCollisionWithTVXAndZVertexAndSel8Counter = 0;
72+
int readCollisionWithTVXAndZVertexAndSel8AndIsGoodZvtxFT0vsPVCounter = 0;
6873
int readCollisionWithTVXAndZVertexAndSel8FullCounter = 0;
74+
int readCollisionWithTVXAndZVertexAndSel8FullAndIsGoodZvtxFT0vsPVCounter = 0;
6975
int readCollisionWithTVXAndZVertexAndSel8FullPbPbCounter = 0;
7076
int readCollisionWithTVXAndZVertexAndSelMCCounter = 0;
7177
int readCollisionWithTVXAndZVertexAndSelMCFullCounter = 0;
@@ -79,8 +85,14 @@ struct LuminosityCalculator {
7985
for (const auto& collisionCount : collisionCounts) {
8086
readCollision += collisionCount.readCounts().front();
8187
readCollisionWithTVXCounter += collisionCount.readCountsWithTVX().front();
88+
readCollisionWithTVXAndSel8Counter += collisionCount.readCountsWithTVXAndSel8().front();
89+
readCollisionWithTVXAndSel8AndIsGoodZvtxFT0vsPVCounter += collisionCount.readCountsWithTVXAndSel8AndIsGoodZvtxFT0vsPV().front();
90+
readCollisionWithTVXAndSel8FullCounter += collisionCount.readCountsWithTVXAndSel8Full().front();
91+
readCollisionWithTVXAndSel8FullAndIsGoodZvtxFT0vsPVCounter += collisionCount.readCountsWithTVXAndSel8FullAndIsGoodZvtxFT0vsPV().front();
8292
readCollisionWithTVXAndZVertexAndSel8Counter += collisionCount.readCountsWithTVXAndZVertexAndSel8().front();
93+
readCollisionWithTVXAndZVertexAndSel8AndIsGoodZvtxFT0vsPVCounter += collisionCount.readCountsWithTVXAndZVertexAndSel8AndIsGoodZvtxFT0vsPV().front();
8394
readCollisionWithTVXAndZVertexAndSel8FullCounter += collisionCount.readCountsWithTVXAndZVertexAndSel8Full().front();
95+
readCollisionWithTVXAndZVertexAndSel8FullAndIsGoodZvtxFT0vsPVCounter += collisionCount.readCountsWithTVXAndZVertexAndSel8FullAndIsGoodZvtxFT0vsPV().front();
8496
readCollisionWithTVXAndZVertexAndSel8FullPbPbCounter += collisionCount.readCountsWithTVXAndZVertexAndSel8FullPbPb().front();
8597
readCollisionWithTVXAndZVertexAndSelMCCounter += collisionCount.readCountsWithTVXAndZVertexAndSelMC().front();
8698
readCollisionWithTVXAndZVertexAndSelMCFullCounter += collisionCount.readCountsWithTVXAndZVertexAndSelMCFull().front();
@@ -98,17 +110,23 @@ struct LuminosityCalculator {
98110
registry.get<TH1>(HIST("counter"))->SetBinContent(4, registry.get<TH1>(HIST("counter"))->GetBinContent(4) + readBCWithTVXAndNoTFBAndNoITSROFB);
99111
registry.get<TH1>(HIST("counter"))->SetBinContent(5, registry.get<TH1>(HIST("counter"))->GetBinContent(5) + readCollision);
100112
registry.get<TH1>(HIST("counter"))->SetBinContent(6, registry.get<TH1>(HIST("counter"))->GetBinContent(6) + readCollisionWithTVXCounter);
101-
registry.get<TH1>(HIST("counter"))->SetBinContent(7, registry.get<TH1>(HIST("counter"))->GetBinContent(7) + readCollisionWithTVXAndZVertexAndSel8Counter);
102-
registry.get<TH1>(HIST("counter"))->SetBinContent(8, registry.get<TH1>(HIST("counter"))->GetBinContent(8) + readCollisionWithTVXAndZVertexAndSel8FullCounter);
103-
registry.get<TH1>(HIST("counter"))->SetBinContent(9, registry.get<TH1>(HIST("counter"))->GetBinContent(9) + readCollisionWithTVXAndZVertexAndSel8FullPbPbCounter);
104-
registry.get<TH1>(HIST("counter"))->SetBinContent(10, registry.get<TH1>(HIST("counter"))->GetBinContent(10) + readCollisionWithTVXAndZVertexAndSelMCCounter);
105-
registry.get<TH1>(HIST("counter"))->SetBinContent(11, registry.get<TH1>(HIST("counter"))->GetBinContent(11) + readCollisionWithTVXAndZVertexAndSelMCFullCounter);
106-
registry.get<TH1>(HIST("counter"))->SetBinContent(12, registry.get<TH1>(HIST("counter"))->GetBinContent(12) + readCollisionWithTVXAndZVertexAndSelMCFullPbPbCounter);
107-
registry.get<TH1>(HIST("counter"))->SetBinContent(13, registry.get<TH1>(HIST("counter"))->GetBinContent(13) + readCollisionWithTVXAndZVertexAndSelUnanchoredMCCounter);
108-
registry.get<TH1>(HIST("counter"))->SetBinContent(14, registry.get<TH1>(HIST("counter"))->GetBinContent(14) + readCollisionWithTVXAndZVertexAndSelTVXCounter);
109-
registry.get<TH1>(HIST("counter"))->SetBinContent(15, registry.get<TH1>(HIST("counter"))->GetBinContent(15) + readCollisionWithTVXAndZVertexAndSel7Counter);
110-
registry.get<TH1>(HIST("counter"))->SetBinContent(16, registry.get<TH1>(HIST("counter"))->GetBinContent(16) + readCollisionWithTVXAndZVertexAndSel7KINT7Counter);
111-
registry.get<TH1>(HIST("counter"))->SetBinContent(16, registry.get<TH1>(HIST("counter"))->GetBinContent(17) + readCollisionWithCustomCounter);
113+
registry.get<TH1>(HIST("counter"))->SetBinContent(7, registry.get<TH1>(HIST("counter"))->GetBinContent(7) + readCollisionWithTVXAndSel8Counter);
114+
registry.get<TH1>(HIST("counter"))->SetBinContent(8, registry.get<TH1>(HIST("counter"))->GetBinContent(8) + readCollisionWithTVXAndSel8AndIsGoodZvtxFT0vsPVCounter);
115+
registry.get<TH1>(HIST("counter"))->SetBinContent(9, registry.get<TH1>(HIST("counter"))->GetBinContent(9) + readCollisionWithTVXAndSel8FullCounter);
116+
registry.get<TH1>(HIST("counter"))->SetBinContent(10, registry.get<TH1>(HIST("counter"))->GetBinContent(10) + readCollisionWithTVXAndSel8FullAndIsGoodZvtxFT0vsPVCounter);
117+
registry.get<TH1>(HIST("counter"))->SetBinContent(11, registry.get<TH1>(HIST("counter"))->GetBinContent(11) + readCollisionWithTVXAndZVertexAndSel8Counter);
118+
registry.get<TH1>(HIST("counter"))->SetBinContent(12, registry.get<TH1>(HIST("counter"))->GetBinContent(12) + readCollisionWithTVXAndZVertexAndSel8AndIsGoodZvtxFT0vsPVCounter);
119+
registry.get<TH1>(HIST("counter"))->SetBinContent(13, registry.get<TH1>(HIST("counter"))->GetBinContent(13) + readCollisionWithTVXAndZVertexAndSel8FullCounter);
120+
registry.get<TH1>(HIST("counter"))->SetBinContent(14, registry.get<TH1>(HIST("counter"))->GetBinContent(14) + readCollisionWithTVXAndZVertexAndSel8FullAndIsGoodZvtxFT0vsPVCounter);
121+
registry.get<TH1>(HIST("counter"))->SetBinContent(15, registry.get<TH1>(HIST("counter"))->GetBinContent(15) + readCollisionWithTVXAndZVertexAndSel8FullPbPbCounter);
122+
registry.get<TH1>(HIST("counter"))->SetBinContent(16, registry.get<TH1>(HIST("counter"))->GetBinContent(16) + readCollisionWithTVXAndZVertexAndSelMCCounter);
123+
registry.get<TH1>(HIST("counter"))->SetBinContent(17, registry.get<TH1>(HIST("counter"))->GetBinContent(17) + readCollisionWithTVXAndZVertexAndSelMCFullCounter);
124+
registry.get<TH1>(HIST("counter"))->SetBinContent(18, registry.get<TH1>(HIST("counter"))->GetBinContent(18) + readCollisionWithTVXAndZVertexAndSelMCFullPbPbCounter);
125+
registry.get<TH1>(HIST("counter"))->SetBinContent(19, registry.get<TH1>(HIST("counter"))->GetBinContent(19) + readCollisionWithTVXAndZVertexAndSelUnanchoredMCCounter);
126+
registry.get<TH1>(HIST("counter"))->SetBinContent(20, registry.get<TH1>(HIST("counter"))->GetBinContent(20) + readCollisionWithTVXAndZVertexAndSelTVXCounter);
127+
registry.get<TH1>(HIST("counter"))->SetBinContent(21, registry.get<TH1>(HIST("counter"))->GetBinContent(21) + readCollisionWithTVXAndZVertexAndSel7Counter);
128+
registry.get<TH1>(HIST("counter"))->SetBinContent(22, registry.get<TH1>(HIST("counter"))->GetBinContent(22) + readCollisionWithTVXAndZVertexAndSel7KINT7Counter);
129+
registry.get<TH1>(HIST("counter"))->SetBinContent(23, registry.get<TH1>(HIST("counter"))->GetBinContent(23) + readCollisionWithCustomCounter);
112130
}
113131
PROCESS_SWITCH(LuminosityCalculator, processCalculateLuminosity, "calculate ingredients for luminosity and fill a histogram", true);
114132
};

0 commit comments

Comments
 (0)