MCAppParticle.cxx

Go to the documentation of this file.

//
// MCAppParticle
//
// MCAppParticle is a base class to store data for a particle primary or
// secondary produced during particle transport.
//
// Author:  S. Kasahara 10/2006
//
// Notes: The principle requirement of this class is that it HasA TParticle
//        with particle data information, since TParticle is the format
//        expected by the VMC.
//

#include <TParticlePDG.h>

#include "MessageService/MsgService.h"
#include "MCApplication/MCAppParticle.h"

ClassImp(MCAppParticle)

CVSID("$Id: MCAppParticle.cxx,v 1.11 2009/06/22 04:11:21 kasahara Exp $");

std::ostream& operator << (std::ostream& os, const MCAppParticle& mp)
                          { return mp.Print(os); }

//_____________________________________________________________________________
MCAppParticle::MCAppParticle(Int_t id,Int_t pdg, Int_t statuscode, 
                    Double_t px, Double_t py, Double_t pz, Double_t e, 
                    Double_t vx, Double_t vy, Double_t vz, Double_t tof,
                    TMCProcess process) 
  : fID(id), fTParticle(pdg,statuscode,-1,-1,-1,-1,
    px,py,pz,e,vx,vy,vz,tof), fProcess(process), 
    fPrintPreface(""),  fPrintOption("")  {
  // Normal constructor, Recommended

}

//_____________________________________________________________________________
MCAppParticle::MCAppParticle(Int_t id,TParticle* particle, TMCProcess process) 
  : fID(id), fTParticle(), fProcess(process), 
    fPrintPreface(""),fPrintOption("") {
  // Normal constructor
  // The input particle is copied, but not adopted

  if ( particle ) fTParticle = *particle;  // copy the input particle

}

//_____________________________________________________________________________
MCAppParticle::MCAppParticle() : fID(-1), fTParticle(), fProcess(kPNoProcess),
                                 fPrintPreface(""), fPrintOption("") {
  // Default constructor

}


//_____________________________________________________________________________
MCAppParticle::~MCAppParticle() {
  // Destructor.  Delete all owned sub-objects

}

//_____________________________________________________________________________
std::ostream& MCAppParticle::Print(std::ostream& os) const {
  // Print particle data on ostream.
  // GetPrintOption() == "full" will print child particles
  // GetPrintOption() == "" prints this particle only
  //
  
  os << GetPrintPreface().c_str() << GetID() << ")";
  os << GetPdgCode() << "/" << fTParticle.GetName()
     << " s: " << GetStatusCode() << "/" 
     << GetStatusCodeName().c_str() 
     << " v: " << fTParticle.Vx() 
     << " "    << fTParticle.Vy() 
     << " "    << fTParticle.Vz() 
     << " t: " << fTParticle.T()
     << " p: " << fTParticle.Px() 
     << " "    << fTParticle.Py()
     << " "    << fTParticle.Pz()
     << " P:"  << GetParentID(0) 
     << ","    << GetParentID(1) 
     << " C:"  << GetChildID(0) 
     << " to " << GetChildID(GetNChildren()-1)
     << "."    << endl;

  if ( !strcmp(GetPrintOption().c_str(),"full") ) {
    Int_t nd = GetNChildren();
    if ( nd > 0 ) {
      os << GetPrintPreface().c_str() << "  " << nd;
      if ( nd == 1 ) os << " Child:" << endl;
      else os  << " Children:" << endl;
      // Recursive loop over children, children's children,...
      for ( Int_t id = 0; id < nd; id++ ) {
        MCAppParticle* cdpart = const_cast<MCAppParticle*>(GetChild(id));
        std::string saveoption  = cdpart->GetPrintOption();
        std::string savepreface = cdpart->GetPrintPreface();
        cdpart -> SetPrintOption(GetPrintOption());
        cdpart -> SetPrintPreface(GetPrintPreface() + "  ");
        os << *(cdpart);
        cdpart -> SetPrintOption(saveoption);
        cdpart -> SetPrintPreface(savepreface);
      }
    }
  }


  return os;

}

//_____________________________________________________________________________
void MCAppParticle::Print(Option_t* /* option */) const {
  // Print particle in form supported by TObject::Print
  // The user can also control the print by specifying a preface using
  // the SetPrintPreface method, e.g. to set an indentation.

  Print(std::cout);
  return;

}

//_____________________________________________________________________________
const MCAppParticle* MCAppParticle::GetChild(UInt_t id) const {
  // Return pointer to id-th child or null if out-of-range

  if ( id >= fChildren.size() ) return 0;
  
  const MCAppParticle* child = fChildren[id];

  return child;
  
}


//_____________________________________________________________________________
const MCAppParticle* MCAppParticle::GetParent(UInt_t id) const {
  // Return pointer to id-th parent (0 or 1) or null if out-of-range

  if ( id >= fParents.size() ) return 0;
  
  const MCAppParticle* parent = fParents[id];

  return parent;

}

//_____________________________________________________________________________
const MCAppParticle* MCAppParticle::GetSibling(UInt_t id) const {
  // Return pointer to id-th sibling or null if out-of-range
  // Siblings are particles all generated in the same production
  // mechanism event.  For example, decay products from a single decay
  // are siblings.

  if ( id >= fSiblings.size() ) return 0;
  
  const MCAppParticle* sibling = fSiblings[id];

  return sibling;

}

//_____________________________________________________________________________
void MCAppParticle::AddChild(const MCAppParticle* child) {
  // Add child to list of indices.

  if ( !child ) {
    MSG("MCApp",Msg::kWarning) 
      << "MCAppParticle::AddChild received null child ptr" << endl;
    return;
  }

  fChildren.push_back(child);
  
}

//_____________________________________________________________________________
void MCAppParticle::AddParent(const MCAppParticle* parent) {
  // Add parent to list of indices.

  if ( !parent ) {
    MSG("MCApp",Msg::kWarning) 
      << "MCAppParticle::AddParent received null parent ptr" << endl;
    return;
  }
  
  fParents.push_back(parent);
  
}

//_____________________________________________________________________________
void MCAppParticle::AddSibling(const MCAppParticle* sibling) {
  // Add sibling to list of indices.
  // Siblings are particles all generated in the same production
  // mechanism event.  For example, decay products from a single decay
  // are siblings.

  if ( !sibling ) {
    MSG("MCApp",Msg::kWarning) 
      << "MCAppParticle::AddSibling received null sibling ptr" << endl;
    return;
  }
  
  fSiblings.push_back(sibling);
  
}

//_____________________________________________________________________________
Double_t MCAppParticle::GetCharge() { 
  return ( fTParticle.GetPDG() ) ?  fTParticle.GetPDG() -> Charge() : 0; 
}

//_____________________________________________________________________________
void MCAppParticle::ProductionVertex(TLorentzVector& lv) const {
  lv.SetXYZT(0,0,0,0);
  fTParticle.ProductionVertex(lv); 
}

//_____________________________________________________________________________
void MCAppParticle::ProductionMomentum(TLorentzVector& lv) const {
  lv.SetPxPyPzE(0,0,0,0);
  fTParticle.Momentum(lv); 
}

//_____________________________________________________________________________
void MCAppParticle::GetPolarisation(TVector3& vec) {
  vec.SetXYZ(0,0,0);
  fTParticle.GetPolarisation(vec); 
}

//_____________________________________________________________________________
Int_t MCAppParticle::GetParentID(UInt_t i) const {
  const MCAppParticle* parent = GetParent(i);
  if ( parent ) return parent->GetID();
  else return -1;
}

//_____________________________________________________________________________
Int_t MCAppParticle::GetChildID(UInt_t i) const {
  const MCAppParticle* child = GetChild(i);
  if ( child ) return child->GetID();
  else return -1;
}

//_____________________________________________________________________________
Int_t MCAppParticle::GetSiblingID(UInt_t i) const {
  const MCAppParticle* sibling = GetSibling(i);
  if ( sibling ) return sibling->GetID();
  else return -1;
}

---