9#include <tracking/trackFindingCDC/findlets/complete/AsicBackgroundLibraryCreator.h>
10#include <tracking/trackFindingCDC/eventdata/hits/CDCWireHit.h>
11#include <tracking/trackFindingCDC/eventdata/tracks/CDCTrack.h>
12#include <tracking/trackFindingCDC/eventdata/trajectories/CDCTrajectory2D.h>
13#include <tracking/trackFindingCDC/eventdata/trajectories/CDCTrajectorySZ.h>
14#include <tracking/trackFindingCDC/topology/CDCWire.h>
15#include <tracking/trackFindingCDC/utilities/StringManipulation.h>
16#include <framework/core/ModuleParamList.templateDetails.h>
17#include <cdc/dataobjects/CDCHit.h>
18#include <framework/logging/Logger.h>
28using namespace TrackFindingCDC;
44 if (wireHit->isAxial()) {
45 recoPos2D = wireHit->reconstruct2D(trajectory2D);
50 const double z = trajectorySZ.
mapSToZ(arcLength);
53 const double driftLength = wireHit->getRefDriftLength();
56 recoPos2D = wirePos2DAtZ + disp2D;
58 const float distanceToHit = trajectory2D.
getDist2D(recoPos2D);
59 return std::abs(distanceToHit);
68 if ((*m_channelMapFromDB).isValid()) {
69 B2DEBUG(29,
"CDC Channel map is valid");
71 B2FATAL(
"CDC Channel map is not valid");
76 auto leavesCreator = [
this](TTree & tree) {
79 tree.Branch(
"ADC", &
m_adc_sig,
"ADC_Sig/S:ADC_bg/S");
82 tree.Branch(
"Board", &
m_board,
"board/b");
83 tree.Branch(
"Channel", &
m_channel,
"channel/b");
84 tree.Branch(
"Nhit", &
m_n_hit,
"nhit/b");
86 "TDC0/S:ADC0/S:TOT0/S:TDC1/S:ADC1/S:TOT1/S:TDC2/S:ADC2/S:TOT2/S:TDC3/S:ADC3/S:TOT3/S:TDC4/S:ADC4/S:TOT4/S:TDC5/S:ADC5/S:TOT5/S:TDC6/S:ADC6/S:TOT6/S:TDC7/S:ADC7/S:TOT7/S");
101 const int isl = cm.getISuperLayer();
102 const int il = cm.getILayer();
103 const int iw = cm.getIWire();
104 const int iBoard = cm.getBoardID();
105 const int iCh = cm.getBoardChannel();
106 const WireID wireId(isl, il, iw);
113 return "Finds suitable ASICs with a single hit attached to a track and uses them to create the library";
119 map< pair<int, int>, vector<const CDCWireHit*>> groupedByAsic;
122 B2ASSERT(
"Channel map NOT found for the channel",
m_map.count(eWire) > 0);
123 auto board =
m_map[eWire].first;
124 auto channel =
m_map[eWire].second;
125 auto asicID = pair<int, int>(board, channel / 8);
126 groupedByAsic[asicID].push_back(&wireHit);
128 for (
auto& asicList : groupedByAsic) {
138 moduleParamList->
addParameter(TrackFindingCDC::prefixed(prefix,
"minimalHitNumberASIC"),
140 "Required number of hits per ASIC for library creation",
143 moduleParamList->
addParameter(TrackFindingCDC::prefixed(prefix,
"AsicLibraryFileName"),
145 "ASIC library file name",
148 moduleParamList->
addParameter(TrackFindingCDC::prefixed(prefix,
"maximalDistanceSignal"),
150 "maximal distance in cm from track to signal hit",
153 moduleParamList->
addParameter(TrackFindingCDC::prefixed(prefix,
"minimalDistanceBackground"),
155 "minimal distance in cm from track to background hit",
158 moduleParamList->
addParameter(TrackFindingCDC::prefixed(prefix,
"useAxialHitsOnly"),
160 "use axial layers only",
163 moduleParamList->
addParameter(TrackFindingCDC::prefixed(prefix,
"writeExtraVars"),
165 "Write extra variables to the library",
168 moduleParamList->
addParameter(TrackFindingCDC::prefixed(prefix,
"minimalHitsOnTrack"),
170 "Required number of hits on track for library creation",
174 moduleParamList->
getParameter<std::string>(
"inputTracks").setDefaultValue(
"CDCTrackVector");
175 moduleParamList->
getParameter<std::string>(
"inputWireHits").setDefaultValue(
"CDCWireHitVector");
191 if (wireHits.size() > 8) {
192 B2ERROR(
"Number of hits per asic should not exceed 8, observe too many hits." <<
LogVar(
"nHits", wireHits.size()));
201 for (
auto& hit : wireHits) {
202 if (!(*hit)->hasBackgroundFlag() && (*hit)->hasTakenFlag()) {
221 const CDCTrack* signalTrack =
nullptr;
222 for (
auto& track : tracks) {
223 for (
auto& hit : track) {
224 if (&hit.getWireHit() == signal) {
225 signalTrack = &track;
229 if (signalTrack !=
nullptr) {
234 if (signalTrack ==
nullptr) {
235 B2DEBUG(29,
"AsicBackgroundLibraryCreator::No track found for the signal hit");
244 const auto& trajectory = signalTrack->getStartTrajectory3D();
248 for (
auto& hit : wireHits) {
250 const float dist = getDist2D(trajectory, hit);
277 for (
auto& hit : wireHits) {
278 auto eWire = hit->getWireID().getEWire();
279 auto channel =
m_map[eWire].second;
280 auto asicCH = channel % 8;
281 m_asic_info[asicCH].ADC = hit->getHit()->getADCCount();
282 m_asic_info[asicCH].TDC = hit->getHit()->getTDCCount();
290 auto eWire = signal->getWireID().getEWire();
293 m_adc_sig = signal->getHit()->getADCCount();
The Module parameter list class.
double m_distance_signal_max
maximal distance from track to signal hit
Short_t m_adc_sig
ADC of the signal.
std::unique_ptr< DBArray< CDCChannelMap > > m_channelMapFromDB
Channel map retrieved from DB.
asicChannel m_asic_info[8]
all 8 channels
UShort_t m_n_hit_track
Number of hits on the track.
void initialize() final
Access database here, open library for writing:
void beginRun() final
Reload channel map if needed.
bool m_use_axial_hits_only
use axial layers only
size_t m_minimal_hits_on_track
min. number of hits on the track
float m_dist_bg
min. distance to non-linked hits
UChar_t m_board
signal boardID
float m_dist_signal
Distance to signal hit.
std::string m_library_name
output library name
std::map< int, std::pair< int, int > > m_map
map from ewire to board/channel ID
std::string getDescription() final
Short description of the findlet.
void selectAsic(const std::vector< const CDCWireHit * > &wireHits, const std::vector< CDCTrack > &tracks)
Algorithm to select suitable ASIC for library creation.
void exposeParameters(ModuleParamList *moduleParamList, const std::string &prefix) final
Expose the parameters to a module.
UChar_t m_n_hit
For debugging, store also number of channels with hits.
UChar_t m_channel
signal channelID
bool m_write_extra_vars
extra vars to the library
std::unique_ptr< Recorder > m_recorder
Recorder for the root output.
double m_distance_background_min
minimal distance from track to background hit
size_t m_minimal_hit_number
min. number of hits in ASIC for background check
Short_t m_adc_max_bg
Max. ADC background.
void terminate() final
write out the library
void apply(const std::vector< CDCWireHit > &wireHits, const std::vector< CDCTrack > &tracks) final
Main algorithm marking hit as background.
Class representing a sequence of three dimensional reconstructed hits.
Particle trajectory as it is seen in xy projection represented as a circle.
double calcArcLength2D(const Vector2D &point) const
Calculate the travel distance from the start position of the trajectory.
Vector2D getClosest(const Vector2D &point) const
Calculates the closest approach on the trajectory to the given point.
double getDist2D(const Vector2D &point) const
Calculates the distance from the point to the trajectory as seen from the xy projection.
Particle full three dimensional trajectory.
CDCTrajectory2D getTrajectory2D() const
Getter for the two dimensional trajectory.
CDCTrajectorySZ getTrajectorySZ() const
Getter for the sz trajectory.
Linear trajectory in sz space.
double mapSToZ(const double s=0) const
Translates the travel distance to the z coordinate.
Class representing a hit wire in the central drift chamber.
const WireID & getWireID() const
Getter for the WireID of the wire the hit is located on.
Class representing a sense wire in the central drift chamber.
Vector2D getWirePos2DAtZ(const double z) const
Gives the xy projected position of the wire at the given z coordinate.
void initialize() override
void terminate() override
virtual void exposeParameters(ModuleParamList *moduleParamList, const std::string &prefix)
Class to fill a tree from a set of variables.
A two dimensional vector which is equipped with functions for correct handling of orientation relate...
double normalizeTo(const double toLength)
Normalizes the vector to the given length.
Class to identify a wire inside the CDC.
unsigned short getEWire() const
Getter for encoded wire number.
Class to store variables with their name which were sent to the logging service.
ModuleParam< T > & getParameter(const std::string &name) const
Returns a reference to a parameter.
void addParameter(const std::string &name, T ¶mVariable, const std::string &description, const T &defaultValue)
Adds a new parameter to the module list.
Abstract base class for different kinds of events.