sepd/sepd_eventplanecalib/EventPlaneData.cc

0001 #include "EventPlaneData.h"
0002
0003 #include <cmath>
0004
0005 EventPlaneData::EventPlaneData()
0006 {
0007   sepd_charge.fill(0);
0008 }
0009
0010 void EventPlaneData::Reset()
0011 {
0012   event_id = 0;
0013   event_zvertex = std::numeric_limits<double>::quiet_NaN();
0014   event_centrality = std::numeric_limits<double>::quiet_NaN();
0015   sepd_totalcharge = std::numeric_limits<double>::quiet_NaN();
0016   sepd_charge.fill(0);
0017 }
0018
0019 void EventPlaneData::identify(std::ostream& os) const
0020 {
0021   os << "--- EventPlaneData Identify ---" << std::endl;
0022   os << "Event ID: " << event_id << std::endl;
0023   os << "Z-Vertex: " << event_zvertex << std::endl;
0024   os << "Centrality: " << event_centrality << std::endl;
0025   os << "sEPD Total Charge: " << sepd_totalcharge << std::endl;
0026   os << "-------------------------------" << std::endl;
0027 }
0028
0029 int EventPlaneData::isValid() const
0030 {
0031   // An object is considered invalid if the Z-vertex is still NaN
0032   // or if the event ID hasn't been set (remains 0).
0033   if (std::isnan(event_zvertex))
0034   {
0035     return 0;
0036   }
0037
0038   if (event_id == 0)
0039   {
0040     return 0;
0041   }
0042
0043   return 1;
0044 }