void EUTelTrackAnalysis::plotIncidenceAngles(EUTelTrack track) {
streamlog_out(DEBUG2) << " EUTelTrackAnalysis::plotIncidenceAngles------------------------------BEGIN"<< std::endl;
std::vector<EUTelState> states = track.getStates();
for(size_t i=0; i<states.size(); ++i) {
EUTelState state = states.at(i);
state.print();
TVectorD stateVec = state.getStateVec();
float incidenceXZ = stateVec[1];
typedef std::map<int , AIDA::IHistogram1D * >::iterator it_type;
for(it_type iterator =_mapFromSensorIDToIncidenceXZ.begin(); iterator != _mapFromSensorIDToIncidenceXZ.end(); iterator++) {
if(iterator->first == state.getLocation()) {
streamlog_out(DEBUG2) << "Add incidence XZ : " << incidenceXZ << std::endl;
_mapFromSensorIDToIncidenceXZ[ state.getLocation() ] ->fill(incidenceXZ);
break;
}
}
}
for(size_t i=0; i<states.size(); ++i) {
EUTelState state = states.at(i);
state.print();
TVectorD stateVec = state.getStateVec();
float incidenceYZ = stateVec[2];
typedef std::map<int , AIDA::IHistogram1D * >::iterator it_type;
for(it_type iterator =_mapFromSensorIDToIncidenceYZ.begin(); iterator != _mapFromSensorIDToIncidenceYZ.end(); iterator++) {
if(iterator->first == state.getLocation()) {
streamlog_out(DEBUG2) << "Add incidence YZ : " << incidenceYZ << std::endl;
_mapFromSensorIDToIncidenceYZ[ state.getLocation() ] ->fill(incidenceYZ);
break;
}
}
}
streamlog_out(DEBUG2) << " EUTelTrackAnalysis::plotIncidenceAngles------------------------------END"<< std::endl;
}
bool EUTelState::operator!=(const EUTelState compareState ) const {
if(getLocation() == compareState.getLocation() and getPosition()[0] == compareState.getPosition()[0] and getPosition()[1] == compareState.getPosition()[1] and getPosition()[2] == compareState.getPosition()[2]){
return false;
}else{
return true;
}
}
std::vector<int> EUTelTrack::getPlaIDs() const {
//Using an iterator did not seem to work?
std::vector<int> planes;
for(size_t i = 0; i < _states.size() ; ++i){
EUTelState state = _states.at(i);
planes.push_back(state.getLocation());
}
return planes;
}
void EUTelTrackAnalysis::plotResidualVsPosition(EUTelTrack track) {
streamlog_out(DEBUG2) << " EUTelTrackAnalysis::plotResidualVsPosition------------------------------BEGIN"<< std::endl;
std::vector<EUTelState> states = track.getStates();
for(size_t i=0; i<states.size(); ++i) {
EUTelState state = states.at(i);
state.print();
if(!state.getIsThereAHit()) {
continue;
}
EVENT::TrackerHit* hit = state.getHit();
float* statePosition = state.getPosition();
const double* hitPosition = hit->getPosition();
float residual[2];
streamlog_out(DEBUG2) << "State position: " << statePosition[0]<<","<<statePosition[1]<<","<<statePosition[2]<< std::endl;
streamlog_out(DEBUG2) << "Hit position: " << hitPosition[0]<<","<<hitPosition[1]<<","<<hitPosition[2]<< std::endl;
typedef std::map<int ,AIDA::IProfile2D* >::iterator it_type;
for(it_type iterator = _mapFromSensorIDToHistogramX.begin(); iterator != _mapFromSensorIDToHistogramX.end(); iterator++) {
if(iterator->first == state.getLocation()) {
residual[0]=statePosition[0]-hitPosition[0];
streamlog_out(DEBUG2) << "Add residual X : " << residual[0]<< std::endl;
_mapFromSensorIDToHistogramX[ state.getLocation() ] -> fill( hitPosition[0], hitPosition[1], residual[0], 1 );
break;
}
}
for(it_type iterator = _mapFromSensorIDToHistogramY.begin(); iterator != _mapFromSensorIDToHistogramY.end(); iterator++) {
if(iterator->first == state.getLocation()) {
residual[1]=statePosition[1]-hitPosition[1];
streamlog_out(DEBUG2) << "Add residual Y : " << residual[1]<< std::endl;
_mapFromSensorIDToHistogramY[ state.getLocation() ] -> fill( hitPosition[0], hitPosition[1], residual[1], 1 );
break;
}
}
}
streamlog_out(DEBUG2) << " EUTelTrackAnalysis::plotResidualVsPosition------------------------------END"<< std::endl;
}
//Overload operators.
bool EUTelState::operator<(const EUTelState compareState ) const {
return getPosition()[2]<compareState.getPosition()[2];
}
