std::vector<JPetHit> TaskD::createHits(std::vector<JPetPhysSignal>& signals)
{
std::vector<JPetHit> hits;
for (auto i = signals.begin(); i != signals.end(); ++i) {
for (auto j = i; ++j != signals.end(); /**/) {
if (i->getPM().getScin() == j->getPM().getScin()) {
// found 2 signals from the same scintillator
JPetHit hit;
if (i->getPM().getSide() == JPetPM::SideA
&& j->getPM().getSide() == JPetPM::SideB) {
hit.setSignalA(*i);
hit.setSignalB(*j);
} else if (j->getPM().getSide() == JPetPM::SideA
&& i->getPM().getSide() == JPetPM::SideB) {
hit.setSignalA(*j);
hit.setSignalB(*i);
} else {
WARNING("TWO hits on the same scintillator side we ignore it");
// if two hits on the same side, ignore
continue;
}
hit.setScintillator(i->getPM().getScin());
getStatistics().getHisto1D("Scins multiplicity").Fill((*i).getPM().getScin().getID());
double dt = hit.getSignalA().getTime() - hit.getSignalB().getTime();
hit.setTimeDiff(dt);
double t = 0.5
* (hit.getSignalA().getTime() + hit.getSignalB().getTime());
hit.setTime(t);
// fill the appropriate dt histogram for this scintillator
getStatistics().getHisto1D(
Form("dt for scin %d", hit.getScintillator().getID())
).Fill(hit.getTimeDiff());
hits.push_back(hit);
// increment the counter of found hits
getStatistics().getCounter("No. found hits")++;
}
}
}
return hits;
}
vector<JPetHit> TaskC::createHits(const vector<JPetRawSignal>& signals)
{
vector<JPetHit> hits;
for (auto i = signals.begin(); i != signals.end(); ++i) {
for (auto j = i; ++j != signals.end();) {
if (i->getPM().getScin() == j->getPM().getScin()) {
// found 2 signals from the same scintillator
// wrap the RawSignal objects into RecoSignal and PhysSignal
// for now this is just wrapping opne object into another
// in the future analyses it will involve more logic like
// reconstructing the signal's shape, charge, amplitude etc.
JPetRecoSignal recoSignalA;
JPetRecoSignal recoSignalB;
JPetPhysSignal physSignalA;
JPetPhysSignal physSignalB;
// assign sides A and B properly
if (
(i->getPM().getSide() == JPetPM::SideA)
&& (j->getPM().getSide() == JPetPM::SideB)
) {
recoSignalA.setRawSignal(*i);
recoSignalB.setRawSignal(*j);
} else if (
(j->getPM().getSide() == JPetPM::SideA)
&& (i->getPM().getSide() == JPetPM::SideB)
) {
recoSignalA.setRawSignal(*j);
recoSignalB.setRawSignal(*i);
} else {
// if two hits on the same side, ignore
WARNING("TWO hits on the same scintillator side we ignore it");
continue;
}
if ( recoSignalA.getRawSignal().getNumberOfPoints(JPetSigCh::Leading) < 4 ) continue;
if ( recoSignalB.getRawSignal().getNumberOfPoints(JPetSigCh::Leading) < 4 ) continue;
bool thresholds_ok = true;
for (int i = 1; i <= 4; ++i) {
if ( recoSignalA.getRawSignal().getTimesVsThresholdNumber(JPetSigCh::Leading).count(i) < 1 ) {
thresholds_ok = false;
}
if ( recoSignalB.getRawSignal().getTimesVsThresholdNumber(JPetSigCh::Leading).count(i) < 1 ) {
thresholds_ok = false;
}
}
if (thresholds_ok == false) {
continue;
}
physSignalA.setRecoSignal(recoSignalA);
physSignalB.setRecoSignal(recoSignalB);
auto leading_points_a = physSignalA.getRecoSignal().getRawSignal().getTimesVsThresholdNumber(JPetSigCh::Leading);
auto leading_points_b = physSignalB.getRecoSignal().getRawSignal().getTimesVsThresholdNumber(JPetSigCh::Leading);
//skip signals with no information on 1st threshold
// if(leading_points_a.count(1) == 0) continue;
// if(leading_points_b.count(1) == 0) continue;
physSignalA.setTime(leading_points_a.at(1));
physSignalB.setTime(leading_points_b.at(1));
JPetHit hit;
hit.setSignalA(physSignalA);
hit.setSignalB(physSignalB);
hit.setScintillator(i->getPM().getScin());
hit.setBarrelSlot(i->getPM().getScin().getBarrelSlot());
physSignalA.setTime(physSignalA.getRecoSignal().getRawSignal().getTimesVsThresholdNumber(JPetSigCh::Leading).at(1));
physSignalB.setTime(physSignalB.getRecoSignal().getRawSignal().getTimesVsThresholdNumber(JPetSigCh::Leading).at(1));
hit.setTime( 0.5 * ( hit.getSignalA().getTime() + hit.getSignalB().getTime()) );
hits.push_back(hit);
getStatistics().getCounter("No. found hits")++;
}
}
}
return hits;
}
