genie::utils::fragmrec Namespace Reference

Simple utilities for the Fragmentation Event Record. More...

Functions
int	NParticles (int pdg_code, const TClonesArray *const particle_list)
int	NParticles (int pdg_code, int status, const TClonesArray *const particle_list)
int	NPositives (const TClonesArray *const particle_list)
int	NNegatives (const TClonesArray *const particle_list)
void	Print (const TClonesArray *const part_list)

Detailed Description

Simple utilities for the Fragmentation Event Record.

The Fragmentation event record is a TClonesArray of TMCParticles - equivalent to PYTHIA's PYJETS.

Author

Costas Andreopoulos <costas.andreopoulos stfc.ac.uk> University of Liverpool & STFC Rutherford Appleton Lab

November 26, 2004

Copyright (c) 2003-2019, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org or see $GENIE/LICENSE

Function Documentation

int genie::utils::fragmrec::NNegatives

(

const TClonesArray *const

particle_list

)

Definition at line 77 of file FragmRecUtils.cxx.

References genie::PDGLibrary::Find(), and genie::PDGLibrary::Instance().

{
// Find out the number of negative particles in the particle container

  TIter piter(part_list);

  TMCParticle * p = 0;
  int nneg = 0;

  while( (p = (TMCParticle *) piter.Next()) )
         if( PDGLibrary::Instance()->Find(p->GetKF())->Charge() < 0 ) nneg++;

  return nneg;
}

int genie::utils::fragmrec::NParticles	(	int	pdg_code,
		const TClonesArray *const	particle_list
	)

Definition at line 32 of file FragmRecUtils.cxx.

Referenced by beamlinerecoana::WCTrackRecoAnalysis::analyze(), and main().

{
  int nparticles = 0;
  TMCParticle* p = 0;

  TObjArrayIter particle_iter(particle_list);

  while( (p = (TMCParticle *) particle_iter.Next()) ) {
    if(p->GetKF() == pdg_code)  {
      if(p->GetKS()<10) nparticles++;
    }
  }
  return nparticles;
}

int genie::utils::fragmrec::NParticles	(	int	pdg_code,
		int	status,
		const TClonesArray *const	particle_list
	)

Definition at line 48 of file FragmRecUtils.cxx.

References fabricate::status.

{
  int nparticles = 0;
  TMCParticle* p = 0;

  TObjArrayIter particle_iter(particle_list);

  while( (p = (TMCParticle *) particle_iter.Next()) )
              if(p->GetKF() == pdg_code && p->GetKS() == status) nparticles++;

  return nparticles;
}

int genie::utils::fragmrec::NPositives

(

const TClonesArray *const

particle_list

)

Definition at line 62 of file FragmRecUtils.cxx.

References genie::PDGLibrary::Find(), and genie::PDGLibrary::Instance().

{
// Find out the number of negative particles in the particle container

  TIter piter(part_list);

  TMCParticle * p = 0;
  int npos = 0;

  while( (p = (TMCParticle *) piter.Next()) )
         if( PDGLibrary::Instance()->Find(p->GetKF())->Charge() > 0 ) npos++;

  return npos;
}

void genie::utils::fragmrec::Print

(

const TClonesArray *const

part_list

)

Definition at line 92 of file FragmRecUtils.cxx.

References MECModelEnuComparisons::i, pINFO, and SLOG.

Referenced by genie::PythiaHadronization::Hadronize(), and main().

{
  TIter piter(part_list);

  unsigned int i=0;
  TMCParticle * particle = 0;

  double sum_px = 0, sum_py = 0, sum_pz = 0, sum_E = 0;

  
  while( (particle = (TMCParticle *) piter.Next()) ) {

    sum_E  += (particle->GetEnergy());
    sum_px += (particle->GetPx());
    sum_py += (particle->GetPy());
    sum_pz += (particle->GetPz());

    SLOG("FragmRecUtils", pINFO)
        << "-> " << i++ << " " << particle->GetName()
        << " KF = " << particle->GetKF()
        << " KS = " << particle->GetKS()
        << " mom = " << particle->GetParent()
        << " kids = {" 
        << particle->GetFirstChild() << ", " << particle->GetLastChild() 
        << "}(E = "  << particle->GetEnergy()
        << ",Px = "  << particle->GetPx()
        << ",Py = "  << particle->GetPy()
        << ",Pz = "  << particle->GetPz() << ")";
  }

  SLOG("FragmRecUtils", pINFO)
       << "SUMS: E = " << sum_E
       << ", px = " << sum_px << ", py = " << sum_py << ", pz = " << sum_pz;

}