10 #include "otsdaq-core/MonicelliInterface/Event.h"
17 using namespace monicelli;
35 void Event::setRawData (
int trig,
const plaqMapDef& theRawData)
37 theRawData_ = theRawData;
41 void Event::clear (
void )
46 clustersHits_.clear() ;
47 trackCandidates_.clear() ;
48 fittedTracksChi2_.clear() ;
49 fittedTrackResiduals_.clear() ;
50 fittedTracks_.clear() ;
51 fittedTracksCovariance_.clear() ;
54 Event::residualsMapDef Event::makeFittedTrackDeviations (
int trackNumber)
56 std::stringstream ss_;
57 fittedTrackDeviations_.clear();
58 if ( trackCandidates_.size() != fittedTracks_.size() )
61 ss_ <<
"track candidates size: " << trackCandidates_.size()
62 <<
" fitted tracks size: " << fittedTracks_.size() ;
68 for(
unsigned int i=0; i < trackCandidates_.size(); i++ )
70 this->makeSingleTrackDeviations(i);
75 this->makeSingleTrackDeviations(trackNumber);
78 return fittedTrackDeviations_;
82 void Event::makeSingleTrackDeviations(
int trackNum )
84 for(Event::alignedHitsCandidateMapDef::iterator pit = trackCandidates_[trackNum].begin(); pit!= trackCandidates_[trackNum].end() ; pit++)
86 double resX = ( (*pit).second[
"x"] - fittedTracks_[trackNum][0]*(*pit).second[
"z"] - fittedTracks_[trackNum][1] ) ;
87 double resY = ( (*pit).second[
"y"] - fittedTracks_[trackNum][2]*(*pit).second[
"z"] - fittedTracks_[trackNum][3] ) ;
88 fittedTrackDeviations_[trackNum][(*pit).first] = std::make_pair(resX, resY);
93 void Event::setClustersHits(
const clustersHitsMapDef& clustersHits)
95 clustersHits_ = clustersHits ;
99 void Event::setClusters(
const clustersMapDef& clusters)
101 clusters_ = clusters;
105 void Event::addUnconstrainedFittedTrack(
unsigned int trackN, std::string detector,
const vectorDef& fittedTrack,
const matrixDef& covarianceMatrix,
double fittedTrackChi2)
107 while(unconstrainedFittedTracks_.size() <= trackN)
109 unconstrainedFittedTracks_ .push_back(std::map<std::string,vectorDef>());
110 unconstrainedFittedTracksCovariance_.push_back(std::map<std::string,matrixDef>());
111 unconstrainedFittedTracksChi2_ .push_back(std::map<std::string,double>());
113 unconstrainedFittedTracks_ [trackN][detector] = fittedTrack;
114 unconstrainedFittedTracksCovariance_[trackN][detector] = covarianceMatrix;
115 unconstrainedFittedTracksChi2_ [trackN][detector] = fittedTrackChi2;