34 B2INFO(
"Creating efficiencies for every layer");
36 unsigned short wireID;
37 unsigned short layerID;
42 efftree->SetBranchAddress(
"wireID", &wireID);
43 efftree->SetBranchAddress(
"layerID", &layerID);
44 efftree->SetBranchAddress(
"z", &z);
45 efftree->SetBranchAddress(
"isFound", &isFound);
48 B2INFO(
"Building empty efficiency objects");
50 B2INFO(
"Wire Topology found");
52 B2INFO(
"CDC Geometry found");
55 unsigned short layerNo = wireLayer.getICLayer();
61 unsigned short nzbins = 30 - layerNo / 7;
63 double widbins[2] = { -0.5, cdcgeo.
nWiresInLayer(layerNo) - 0.5};
64 double zbins[2] = {wireLayer.getBackwardZ(), wireLayer.getForwardZ()};
66 TEfficiency* effInLayer =
new TEfficiency(Form(
"effLayer_%d", layerNo), Form(
"Hit efficiency pf L%d ; z (cm) ; IWire ", layerNo),
67 nzbins, zbins[0], zbins[1], nwidbins, widbins[0], widbins[1]);
69 B2INFO(
"Teff for layer " << layerNo <<
" was successfully listed.");
71 TFile* outputCollection =
new TFile(
"LayerEfficiencies.root",
"RECREATE");
74 B2INFO(
"Filling the efficiencies");
75 const Long64_t nEntries = efftree->GetEntries();
76 B2INFO(
"Number of entries in tree: " << nEntries);
77 for (Long64_t i = 0; i < nEntries; i++) {
79 TEfficiency* efficiencyInLayer = (TEfficiency*)
m_efficiencyList->At(layerID);
80 efficiencyInLayer->Fill(isFound, z, wireID);
83 B2INFO(
"TEfficiencies successfully filled.");
84 outputCollection->Close();
85 B2INFO(
"TEfficiencies successfully saved");
87 if (nEntries > 1000)
return true;
93 B2INFO(
"Beginning detection of bad wires");
95 B2INFO(
"Wire Topology found");
98 unsigned short layerNo = wireLayer.getICLayer();
99 B2INFO(
"Checking layer " << layerNo);
103 auto passed = (TH2F*)efficiencyInLayer->GetPassedHistogram();
104 auto total = (TH2F*)efficiencyInLayer->GetTotalHistogram();
107 if (!total->GetEntries())
continue;
110 double minFitRange = wireLayer.getBackwardZ() + 30;
111 double maxFitRange = wireLayer.getForwardZ() - 30;
112 TF1* constantFunction =
new TF1(
"constantFunction",
"pol0", minFitRange, maxFitRange);
115 auto passedProjectedX = passed->ProjectionX(
"projectionx1");
116 auto totalProjectedX = total->ProjectionX(
"projectionx2");
117 TEfficiency* efficiencyProjectedX =
new TEfficiency(*passedProjectedX, *totalProjectedX);
119 unsigned short minFitBin = passedProjectedX->FindBin(minFitRange) + 1;
120 unsigned short maxFitBin = passedProjectedX->FindBin(maxFitRange) - 1;
123 auto passedProjectedY = passed->ProjectionY(
"projectiony1", minFitRange, maxFitRange);
124 auto totalProjectedY = total->ProjectionY(
"projectiony2", minFitRange, maxFitRange);
125 TEfficiency* efficiencyProjectedY =
new TEfficiency(*passedProjectedY, *totalProjectedY);
128 float totalAverage = 0;
129 int nonZeroWires = 0;
130 for (
int i = 0; i <= passedProjectedY->GetNbinsX(); ++i) {
131 float efficiencyAtBin = efficiencyProjectedY->GetEfficiency(i);
132 if (efficiencyAtBin > 0.2) {
133 totalAverage += efficiencyAtBin;
138 nonZeroWires > 0 ? totalAverage /= nonZeroWires : totalAverage == 0;
140 TGraphAsymmErrors* graphEfficiencyProjected = efficiencyProjectedX->CreateGraph();
142 B2INFO(
"Successfully determined average properties of layer " << layerNo);
143 B2INFO(
"Looping over wires");
146 for (
int i = 0; i <= passed->GetNbinsY(); ++i) {
149 auto singleWirePassed = passed->ProjectionX(
"single wire projection passed", i, i);
150 auto singleWireTotal = total->ProjectionX(
"single wire projection total", i, i);
153 if (!singleWirePassed->Integral(minFitBin, maxFitBin)) {
154 double wireID = passed->GetYaxis()->GetBinCenter(i);
156 if (wireID < 0 || wireID >= maxWires) {
157 B2ERROR(
"Invalid wireID: " << wireID <<
" for LayerID: " << layerNo
158 <<
". Max wires: " << maxWires);
165 TEfficiency* singleWireEfficiency =
new TEfficiency(*singleWirePassed, *singleWireTotal);
166 TGraphAsymmErrors* graphSingleWireEfficiency = singleWireEfficiency->CreateGraph();
169 TFitResultPtr singleWireFitResults = graphSingleWireEfficiency->Fit(constantFunction,
"SQR");
170 double singleWireEfficiencyFromFit = (singleWireFitResults.Get())->Value(0);
171 double singleWireUpErrorFromFit = (singleWireFitResults.Get())->UpperError(0);
174 float p_value =
chiTest(graphSingleWireEfficiency, graphEfficiencyProjected, minFitRange, maxFitRange);
176 bool averageCondition = (0.95 * totalAverage) > (singleWireEfficiencyFromFit + 3 * singleWireUpErrorFromFit);
177 bool pvalueCondition = p_value < 0.01;
178 bool generalCondition = singleWireEfficiencyFromFit < 0.4;
179 if (generalCondition || (averageCondition && pvalueCondition)) {
180 double wireID = passed->GetYaxis()->GetBinCenter(i);
181 m_badWireList->setWire(layerNo, round(wireID), singleWireEfficiencyFromFit);
184 B2INFO(
"Bad wires for " << layerNo <<
" recorded");
188 B2INFO(
"Bad wire list successfully saved.");
196 unsigned short ndof = 0;
199 int numOfEntries1 = graph1->GetN();
200 int numOfEntries2 = graph2->GetN();
204 for (
int index1 = 0; index1 < numOfEntries1; ++index1) {
206 if (graph1->GetX()[index1] < minValue)
continue;
207 if (graph1->GetX()[index1] > maxValue)
continue;
208 for (
int index2 = 0; index2 < numOfEntries2; ++index2) {
209 if (graph1->GetX()[index1] == graph2->GetX()[index2]) {
211 double chiNumerator = pow(graph1->GetY()[index1] - graph2->GetY()[index2], 2);
212 double chiDenominator = pow(graph1->GetErrorYhigh(index1), 2) + pow(graph1->GetErrorYlow(index1), 2);
213 chi += chiNumerator / chiDenominator;
219 return TMath::Prob(chi, ndof);