void JPetRawSignal::addPoint(const JPetSigCh& sigch) { if (sigch.getType() == JPetSigCh::Trailing) { fTrailingPoints.push_back(sigch); } else if (sigch.getType() == JPetSigCh::Leading) { fLeadingPoints.push_back(sigch); } else if (sigch.getType() == JPetSigCh::Charge) { fTOTPoint = sigch; } }
void TimeWindowCreator::exec() { //getting the data from event in apropriate format //const is commented because this class has inproper architecture: // all get-methods aren't tagged with const modifier if (auto evt = dynamic_cast </*const*/ EventIII * const > (getEvent())) { int ntdc = evt->GetTotalNTDCChannels(); getStatistics().getHisto1D("ChannelsPerEvt").Fill( ntdc ); JPetTimeWindow tslot; tslot.setIndex(fCurrEventNumber); auto tdcHits = evt->GetTDCChannelsArray(); for (int i = 0; i < ntdc; ++i) { //const is commented because this class has inproper architecture: // all get-methods aren't tagged with const modifier auto tdcChannel = dynamic_cast </*const*/ TDCChannel * const > (tdcHits->At(i)); auto tomb_number = tdcChannel->GetChannel(); if (tomb_number % 65 == 0) { // skip trigger signals from TRB continue; } if ( getParamBank().getTOMBChannels().count(tomb_number) == 0 ) { WARNING(Form("DAQ Channel %d appears in data but does not exist in the setup from DB.", tomb_number)); continue; } JPetTOMBChannel& tomb_channel = getParamBank().getTOMBChannel(tomb_number); // one TDC channel may record multiple signals in one TSlot // iterate over all signals from one TDC channel // analyze number of hits per channel getStatistics().getHisto1D("HitsPerEvtCh").Fill( tdcChannel->GetHitsNum() ); const int kNumHits = tdcChannel->GetHitsNum(); for (int j = 0; j < kNumHits; ++j) { // check for unreasonable times // the times should be negative (measured w.r.t end of time window) // and not smaller than -1*timeWindowWidth (which can vary for different) // data but shoudl not exceed 1 ms, i.e. 1.e6 ns) if ( tdcChannel->GetLeadTime(j) > fMaxTime || tdcChannel->GetLeadTime(j) < fMinTime )continue; if ( tdcChannel->GetTrailTime(j) > fMaxTime || tdcChannel->GetTrailTime(j) < fMinTime )continue; JPetSigCh sigChTmpLead = generateSigCh(tomb_channel, JPetSigCh::Leading); JPetSigCh sigChTmpTrail = generateSigCh(tomb_channel, JPetSigCh::Trailing); // finally, set the times in ps [raw times are in ns] sigChTmpLead.setValue(tdcChannel->GetLeadTime(j) * 1000.); sigChTmpTrail.setValue(tdcChannel->GetTrailTime(j) * 1000.); tslot.addCh(sigChTmpLead); tslot.addCh(sigChTmpTrail); } } saveTimeWindow(tslot); fCurrEventNumber++; } }
JPetSigCh TimeWindowCreator::generateSigCh(const JPetTOMBChannel& channel, JPetSigCh::EdgeType edge) const { JPetSigCh sigch; sigch.setDAQch(channel.getChannel()); sigch.setType(edge); sigch.setThresholdNumber(channel.getLocalChannelNumber()); sigch.setThreshold(channel.getThreshold()); sigch.setPM(channel.getPM()); sigch.setFEB(channel.getFEB()); sigch.setTRB(channel.getTRB()); sigch.setTOMBChannel(channel); return sigch; }
void JPetRawSignal::setTOTPoint(const JPetSigCh & totSigCh) { if (totSigCh.getType() == JPetSigCh::Charge) { fTOTPoint = totSigCh; } }
/** * Sets up Signal Channel fields */ JPetSigCh TimeWindowCreatorTools::generateSigCh( double tdcChannelTime, const JPetTOMBChannel& channel, map<unsigned int, vector<double>>& timeCalibrationMap, map<unsigned int, vector<double>>& thresholdsMap, JPetSigCh::EdgeType edge, bool setTHRValuesFromChannels ) { JPetSigCh sigCh; sigCh.setValue(1000.*(tdcChannelTime + UniversalFileLoader::getConfigurationParameter(timeCalibrationMap, channel.getChannel()) )); sigCh.setType(edge); sigCh.setTOMBChannel(channel); sigCh.setPM(channel.getPM()); sigCh.setFEB(channel.getFEB()); sigCh.setTRB(channel.getTRB()); sigCh.setDAQch(channel.getChannel()); sigCh.setThresholdNumber(channel.getLocalChannelNumber()); if(setTHRValuesFromChannels) { sigCh.setThreshold(channel.getThreshold()); } else { sigCh.setThreshold( UniversalFileLoader::getConfigurationParameter(thresholdsMap, channel.getChannel()) ); } return sigCh; }