MCTruthUtil Class Reference

#include <MCTruthUtil.h>

Static Public Member Functions
static bool	isInBTL (const G4Track *)
static void	primary (G4Track *aPrimary)
static void	secondary (G4Track *aSecondary, const G4Track &mother, int)

Detailed Description

Definition at line 13 of file MCTruthUtil.h.

Member Function Documentation

isInBTL()

bool MCTruthUtil::isInBTL ( const G4Track *  aTrack )

static

Definition at line 79 of file MCTruthUtil.cc.

References MillePedeFileConverter_cfg::out, and simpleEdmComparison::tName.

Referenced by secondary().

                                               {
  bool out = false;
  G4String tName(aTrack->GetVolume()->GetLogicalVolume()->GetRegion()->GetName());
  if (tName == "FastTimerRegionBTL" || tName == "FastTimerRegionSensBTL") {
    out = true;
  }
  return out;
}

primary()

void MCTruthUtil::primary ( G4Track *  aPrimary )

static

Definition at line 7 of file MCTruthUtil.cc.

References TrackInformation::setGenParticleP(), TrackInformation::setGenParticlePID(), TrackInformation::setMCTruthID(), TrackInformation::setPrimary(), and TrackInformation::setStoreTrack().

Referenced by KillSecondariesStackingAction::ClassifyNewTrack(), and StackingAction::ClassifyNewTrack().

                                         {
  TrackInformation *trkInfo = new TrackInformation();
  trkInfo->setPrimary(true);
  trkInfo->setStoreTrack();
  trkInfo->setGenParticlePID(aTrack->GetDefinition()->GetPDGEncoding());
  trkInfo->setGenParticleP(aTrack->GetMomentum().mag());
  trkInfo->setMCTruthID(aTrack->GetTrackID());
  aTrack->SetUserInformation(trkInfo);
}

secondary()

void MCTruthUtil::secondary ( G4Track *  aSecondary, const G4Track &  mother, int  flag  )

static

Definition at line 17 of file MCTruthUtil.cc.

References RemoveAddSevLevel::flag, isInBTL(), and LogTrace.

Referenced by KillSecondariesStackingAction::ClassifyNewTrack(), and StackingAction::ClassifyNewTrack().

                                                                            {
  auto motherInfo = static_cast<const TrackInformation *>(mother.GetUserInformation());
  auto trkInfo = new TrackInformation();

  // Take care of cascade decays
  if (flag == 1) {
    trkInfo->setPrimary(true);
    trkInfo->setStoreTrack();
    trkInfo->setGenParticlePID(aTrack->GetDefinition()->GetPDGEncoding());
    trkInfo->setGenParticleP(aTrack->GetMomentum().mag());
    trkInfo->setMCTruthID(aTrack->GetTrackID());
  } else {
    // secondary
    trkInfo->setGenParticlePID(motherInfo->genParticlePID());
    trkInfo->setGenParticleP(motherInfo->genParticleP());
    trkInfo->setMCTruthID(motherInfo->mcTruthID());
  }

  // Store if decay or conversion
  if (flag > 0) {
    trkInfo->setStoreTrack();
    trkInfo->setIDonCaloSurface(aTrack->GetTrackID(),
                                motherInfo->getIDCaloVolume(),
                                motherInfo->getIDLastVolume(),
                                aTrack->GetDefinition()->GetPDGEncoding(),
                                aTrack->GetMomentum().mag());
  } else {
    // transfer calo ID from mother (to be checked in TrackingAction)
    trkInfo->setIDonCaloSurface(motherInfo->getIDonCaloSurface(),
                                motherInfo->getIDCaloVolume(),
                                motherInfo->getIDLastVolume(),
                                motherInfo->caloSurfaceParticlePID(),
                                motherInfo->caloSurfaceParticleP());
  }

  // for Run1 and Run2
  if (motherInfo->hasCastorHit()) {
    trkInfo->setCastorHitPID(motherInfo->getCastorHitPID());
  }

  // for MTD
  if (!trkInfo->isPrimary() && !isInBTL(aTrack)) {
    trkInfo->setExtSecondary();
  }
  if (motherInfo->isExtSecondary()) {
    trkInfo->setExtSecondary();
  }
  if (motherInfo->isBTLlooper()) {
    trkInfo->setBTLlooper();
  }
  if (motherInfo->isInTrkFromBackscattering()) {
    trkInfo->setInTrkFromBackscattering();
  }

  aTrack->SetUserInformation(trkInfo);
#ifdef EDM_ML_DEBUG
  LogTrace("SimG4CoreApplication") << "MCTruthUtil called for " << aTrack->GetTrackID() << " mother "
                                   << motherInfo->isPrimary() << " flag " << flag;
  trkInfo->Print();
#endif
}