PParticle.h

Go to the documentation of this file.

/* Copyright (C) 1998-2019 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Marios A. Kagarlis, Florian Uhlig [committer] */
 
// Author: Marios A. Kagarlis
// Written: 15.11.98
// Revised: 20.07.99
// Revised: 30.08.00  R. Holzmann
// Revised: 07.09.00  J. Ritman (added parent-,siblings-,daughter- index
//                    and decayTime)
// Revised: 15.12.00  R. Holzmann  (creation time added)
// Revised: 22.03.05  R. Holzmann  (get/setParent() added)
// Revised: 23.07.07  IF (new framework)
 
// PParticle Class Header
 
#ifndef _PPARTICLE_H_
#define _PPARTICLE_H_
 
#define MAX_DAUGHTERS 7
 
#include "TLorentzVector.h"
#include "TMath.h"
#include "TString.h"
 
#include "PValues.h"
 
class PParticle : public TLorentzVector {
 
public:
  PParticle(Int_t id = 0, Double_t T = 0., Double_t w = 1.);
  // id, lab kinetic energy (GeV), weight
 
  /*
  PParticle(const char *, Double_t T=0., Double_t w=1.);
  // name, lab kinetic energy (GeV), weight
 
  PParticle(Int_t, Double_t, Double_t, Double_t, Double_t m=0., Double_t w=1.);
  // id, Px, Py, Pz (GeV/c), mass (GeV/c**2) overrides default, weight
 
  PParticle(const char *, Double_t, Double_t, Double_t,Double_t m=0., Double_t w=1.);
  // name, Px, Py, Pz (GeV/c), mass (GeV/c**2) overrides default, weight
 
  PParticle(Int_t, const TVector3 &, Double_t m=0., Double_t w=1.);
  // id, 3-momentum vector (GeV/c), mass (GeV/c**2) overrides default, weight
 
  PParticle(Int_t, Double_t *, Double_t w=1.);
  // id, pointer to 4-dim array (Px, Py, Pz, E) (GeV/c, GeV), weight
  
  PParticle(Int_t, float *, Double_t w=1.);
  // id, pointer to 4-dim array (Px, Py, Pz, E) (GeV/c, GeV), weight
 
  PParticle(const PParticle &);
  // copy constructor
 
  PParticle(const PParticle *);
  // copy constructor
*/
 
  virtual ~PParticle()
  {
    if (qParticle1) delete qParticle1;
    if (qParticle2) delete qParticle2;
    if (values) delete values;
  }
 
  Int_t ID() const { return pid; }
 
  TLorentzVector Vect4() const { return TLorentzVector(Vect(), E()); }
 
  void SetVertex(Double_t x, Double_t y, Double_t z, Double_t t)
  {
    fV.SetXYZ(x, y, z);
    fT = t;
  }
 
  inline Int_t GetParentIndex() const { return parentIndex; }
  inline Int_t GetDaughterIndex() const { return daughterIndex; }
 
  inline TVector3& GetVertex() { return fV; }
  inline TVector3& getVertex() { return fV; }
 
  inline void SetParentId(Int_t pId) { parentId = pId; }
  inline void SetSourceId(Int_t sId) { sourceId = sId; }
  inline void SetDaughterIndex(Int_t dInd) { daughterIndex = dInd; }
  inline void SetSiblingIndex(Int_t sInd) { siblingIndex = sInd; }
  inline void SetParentIndex(Int_t pInd) { parentIndex = pInd; }
  inline void SetSibling(PParticle* p) { sParticle = p; }
 
  inline void ResetDaughters(void)
  {
    for (int i = 0; i < (MAX_DAUGHTERS + 1); i++)
      daughters[i] = nullptr;
  };
 
protected:
  void defaults(void);
 
  Int_t pid;                    //  particle code  (partially conforms with Geant 3)
  Int_t sourceId;               //  Source ID
  Int_t parentId;               //  parent ID
  Int_t parentIndex;            //  parent index in particle array
  Int_t decayModeIndex;         //  decay Mode index (for decayAll option)
  Int_t destroyDecayModeIndex;  // save only for data file
  Int_t daughterIndex;          //  daughter index
  Int_t siblingIndex;           //  sibling index
  Int_t spectator;              
 
  Double_t decayTime;  //  proper time until decay  (in mm/c)
  Double_t wt;         //  weight
  Double_t genwt;      
  Double_t invgenwt;   
  Double_t mult;       
 
  TVector3 fV;                              //  creation vertex (in mm)
  Double_t fT;                              //  creation time (in mm/c)
  Bool_t active;                            
  Int_t index;                              
  PParticle* pParticle;                     
  PParticle* qParticle1;                    
  PParticle* qParticle2;                    
  PParticle* sParticle;                     
  PParticle* daughters[MAX_DAUGHTERS + 1];  
  TString debug;                            
 
  PValues* values;  
  Int_t status;     
 
  Bool_t make_new_qParticle;  
 
  ClassDef(PParticle, 4)  // Pluto Particle Class
};
#endif  // _PPARTICLE_H_