si_syst::ParticleSIManager Class Reference

Manager. More...

#include <SecondaryInteractionSystematic.hxx>

Inheritance diagram for si_syst::ParticleSIManager:

Public Member Functions
	ParticleSIManager ()
	ctor
virtual	~ParticleSIManager ()
	dtor
virtual void	SetPropagator (ParticleSIPropagator *propagator)
	Set the propagator.
si_syst::ParticleSIPropagator *	GetPropagator () const
	Get propagator.
si_syst::SISystInput *	CollectParticleInteractions (AnaTrueParticleB **, int, const SubDetId::SubDetEnum &) const
	Calculates the information needed to compute an event weight.

Protected Member Functions
virtual si_syst::SISystInput *	GetRelevantParticles (AnaTrueParticleB **, int) const

Protected Attributes
ParticleSIPropagator *	_propagator

Detailed Description

Manager.

Definition at line 433 of file SecondaryInteractionSystematic.hxx.

Member Function Documentation

GetRelevantParticles()

si_syst::SISystInput * si_syst::ParticleSIManager::GetRelevantParticles ( AnaTrueParticleB **  allTrajs, int  nTrajs  ) const

protectedvirtual

Get relevant particles Given array of input true tracks, its size The propagator provides the information on relevant tracks, interactions and VOI, so one has to set them first Default method just splits between inelastic and no-interaction in VOI

Fill the interaction types

Definition at line 560 of file SecondaryInteractionSystematic.cxx.

                                                                                                               {
  //********************************************************************

#ifdef DEBUG
  std::cout << " ======= " << std::endl;
  std::cout << " si_syst::ParticleSIManager::GetRelevantParticles " << std::endl;
#endif 

  assert(GetPropagator());
 
  assert(GetPropagator()->IsInitialized());
  
  
  si_syst::SISystInput* result = new si_syst::SISystInput();

  // We want to count all particles of interest (PDG) that originate in particular volume 
  // and group them (if exist) with the daughters of their probable interactions
  si_syst::ParentsAndDaughters ParticlesAndInteractions;

  for (Int_t it = 0; it < nTrajs; it++){

    AnaTrueParticleB* track = allTrajs[it];

    // Should be inside PDGs of interest
    // not using set::find here since then need to cast int to PdgEnum,  not sure it will always succeed
    bool ok = false;
    std::set<ParticleId::PdgEnum>::const_iterator it_set = GetPropagator()->GetParticlePDGs().begin();

    for (; it_set != GetPropagator()->GetParticlePDGs().end(); it_set++){
      if (*it_set == track->PDG){
        ok = true;
        break;
      } 
    }

    if (!ok)
      continue;

    // Should start in the VOI
    if (!GetPropagator()->InVOI(track->Position)) continue;

    ParticlesAndInteractions.push_back(std::make_pair(track, std::vector<AnaTrueParticleB*>()));

    // Now loop over trajectories and find those that mention this particle as parent 
    // those daughtes should start in the volume of interest as well
    for (Int_t jt = 0; jt < nTrajs; jt++){
      AnaTrueParticleB* track2 = allTrajs[jt];

      if (track->ID == track2->ParentID && 
          GetPropagator()->InVOI(track->Position)){

        // Add this trajectory to interaction list
        ParticlesAndInteractions.back().second.push_back(track2);

      }

    } //second loop over true tracks
  }//first loop over true tracks

  //now we have all primary particles grouped together with their possible daughters
  //we can so arrange them into interactions

  result->nParticles = ParticlesAndInteractions.size();

  result->Particles  = new si_syst::ParticleHistory[result->nParticles];


  /// Fill the interaction types
  for(Int_t i = 0; i < result->nParticles; i++){

    result->Particles[i] = si_syst::ParticleHistory(si_syst::kInelastic, *(ParticlesAndInteractions[i].first));

    // No interactions in the volume
    if (ParticlesAndInteractions[i].second.size() == 0)
      result->Particles[i].intType = si_syst::kNoInteraction;  

  }

#ifdef DEBUG
  std::cout << " nParticles: " << result->nParticles << std::endl;
#endif 

#ifdef DEBUG
  std::cout << " ======= " << std::endl;
#endif 
  return result;
}

---

The documentation for this class was generated from the following files:

• /hep/T2K/nd280rep/ANALYSIS/HIGHLAND2/ANAPARTICLE/psyche/psycheND280Utils/v3r12/src/SecondaryInteractionSystematic.hxx
• /hep/T2K/nd280rep/ANALYSIS/HIGHLAND2/ANAPARTICLE/psyche/psycheND280Utils/v3r12/src/SecondaryInteractionSystematic.cxx