EventSummaryFill.cpp

Go to the documentation of this file.

#include "Containers/EventSummary.h"

namespace NAIA {
bool EventSummaryData::Fill(AMSEventR *evPtr) {
  NAntiCluster = evPtr->NAntiCluster();
  NTofCluster = evPtr->NTofCluster();
  NTrRecHit = evPtr->NTrRecHit();
  NTrTrack = evPtr->NTrTrack();
  NRichRing = evPtr->NRichRing();
  NEcalShower = evPtr->NEcalShower();
  NParticle = evPtr->NParticle();

  HighestEcalEnergyParticle = static_cast<short>(std::distance(
      begin(evPtr->Particle()),
      std::max_element(begin(evPtr->Particle()), end(evPtr->Particle()), [](ParticleR &p1, ParticleR &p2) {
        if (p2.pEcalShower() == p1.pEcalShower())
          return false;

        double e1 = p1.pEcalShower() ? p1.pEcalShower()->EnergyD : 0;
        double e2 = p2.pEcalShower() ? p2.pEcalShower()->EnergyD : 0;
        return e1 < e2;
      })));

  HighestBetaHChargeParticle = static_cast<short>(std::distance(
      begin(evPtr->Particle()),
      std::max_element(begin(evPtr->Particle()), end(evPtr->Particle()), [](ParticleR &p1, ParticleR &p2) {
        if (p2.pBetaH() == p1.pBetaH())
          return false;

        int nlayers = 0;
        float qrms = 0;

        double q1 = p1.pBetaH() ? p1.pBetaH()->GetQ(nlayers, qrms, 2, TofClusterHR::DefaultQOptIon) : 0;
        double q2 = p2.pBetaH() ? p2.pBetaH()->GetQ(nlayers, qrms, 2, TofClusterHR::DefaultQOptIon) : 0;
        return q1 < q2;
      })));

  HighestInnerChargeParticle = static_cast<short>(std::distance(
      begin(evPtr->Particle()),
      std::max_element(begin(evPtr->Particle()), end(evPtr->Particle()), [](ParticleR &p1, ParticleR &p2) {
        if (p2.pTrTrack() == p1.pTrTrack())
          return false;

        double q1 = p1.pTrTrack() ? p1.pTrTrack()->GetInnerQ_all().Mean : 0;
        double q2 = p2.pTrTrack() ? p2.pTrTrack()->GetInnerQ_all().Mean : 0;
        return q1 < q2;
      })));

  if (evPtr->pLevel1(0)) {
    JMembPatt = evPtr->pLevel1(0)->JMembPatt;
    if (_isMC) {
      // how not to design a datamodel:
      // PhysBPatt only uses the lower 8 bits to store information
      // however the API requires an int and only an int
      // so if you want to save space you have to cast the result
      // Another excellent demonstration of how using parameters
      // as return values in c++ is inherently stupid
      int dummy = 0;
      int otherDummy = 0;

      evPtr->pLevel1(0)->RebuildTrigPatt(dummy, otherDummy);
      PhysBPatt = static_cast<unsigned short>(otherDummy);
    } else {
      // evPtr->pLevel1(0)->RestorePhysBPat();
      PhysBPatt = evPtr->pLevel1(0)->PhysBPatt;
    }

    // neat trick to count how many bits are set in AntiPatt
    NAcc = std::bitset<8>(evPtr->pLevel1(0)->AntiPatt).count();
  }

  return true;
}
} // namespace NAIA