TMBMCpart.hpp

Go to the documentation of this file.

#ifndef TMBMCpart_H
#define TMBMCpart_H

// file: TMBMCpart.hpp                                                  //
// TMBTree MCpart class 
//    has MCparticle information
//
// created: 21-March-2002  Serban Protopopescu
//                                                                      //

#include "tmb_tree/TMBLorentzVector.hpp"

#ifndef ROOT_TRef
#include "TRef.h"
#endif
#ifndef ROOT_TRefArray
#include "TRefArray.h"
#endif

#include <cmath>
class TMBMCvtx;

class TMBMCpart : public TMBLorentzVector {

private:
  Float_t        _q;            
  Int_t          _pdgid;        
  Int_t          _status;       
  TRef           _vtx;          
  TRef           _vtxend;       

public:
  TMBMCpart() {;}
  TMBMCpart(Float_t px, Float_t py, Float_t pz,Float_t E,
            Float_t q, Float_t pt, Float_t phi, Float_t eta, 
            Int_t pdgid, Int_t status, TRef vtx, TRef vtxend);

  Float_t charge()   const    {return _q;}

  inline Int_t pdgid() const {return _pdgid;}
  inline Int_t abspdgid() const {return pdgid() > 0 ? pdgid() : -pdgid();}
  const TMBMCvtx*       getPMCvtx() const {return _vtx.IsValid() ? (TMBMCvtx*)_vtx.GetObject() : 0;}
  const TMBMCvtx*       getDMCvtx() const {return _vtxend.IsValid() ? (TMBMCvtx*)_vtxend.GetObject() : 0;}
  bool isStable() const {bool stable=true; if(getDMCvtx()) stable=false; return stable;}
  bool isStableParticle() const {bool stable=true; if(getDMCvtx()) stable=false; if ( _status > 1 || _status < 0 ) stable=false; return stable;}
  bool isLepton()const{return abs(_pdgid)>10 && abs(_pdgid)<19;}

  Int_t status() const { return _status; }
  Float_t zFrag() const { return Float_t(_status/10)/1000.; }

  ~TMBMCpart() {;}

   ClassDef(TMBMCpart, 2)   //TMBTree MCpart class
};

#endif

---