TMCParticleStatus.h

Go to the documentation of this file.

// -----------------------------------------------------------------------
// Copyright (C) 2019 CERN and copyright holders of VMC Project.
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "LICENSE".
//
// See https://github.com/vmc-project/vmc for full licensing information.
// -----------------------------------------------------------------------
 
// Authors: Benedikt Volkel 07/03/2019
 
/*************************************************************************
 * Copyright (C) 2019, Rene Brun and Fons Rademakers.                    *
 * Copyright (C) 2019, ALICE Experiment at CERN.                         *
 * All rights reserved.                                                  *
 *                                                                       *
 * For the licensing terms see $ROOTSYS/LICENSE.                         *
 * For the list of contributors see $ROOTSYS/README/CREDITS.             *
 *************************************************************************/
 
#ifndef ROOT_TMCParticleStatus
#define ROOT_TMCParticleStatus
 
// Class TMCParticleStatus
// ---------------------
// additional information on the current status of a TParticle
//
 
#include <iostream>
 
#include "TVector3.h"
#include "TLorentzVector.h"
#include "TParticle.h"
#include "TError.h"
 
struct TMCParticleStatus {

   TMCParticleStatus() = default;

   void InitFromParticle(const TParticle *particle)
   {
      particle->ProductionVertex(fPosition);
      particle->Momentum(fMomentum);
      particle->GetPolarisation(fPolarization);
      fWeight = particle->GetWeight();
   }
 
   virtual ~TMCParticleStatus() = default;
 
   //
   // verbosity
   //

   void Print() const
   {
      ::Info("Print", "Status of track");
      std::cout << "\t"
                << "ID: " << fId << "\n"
                << "\t"
                << "parentID: " << fParentId << "\n"
                << "\t"
                << "weight: " << fWeight << "\n"
                << "\t"
                << "geo state index: " << fGeoStateIndex << "\n"
                << "\t"
                << "step number: " << fStepNumber << "\n"
                << "\t"
                << "track length: " << fTrackLength << "\n"
                << "\t"
                << "position" << std::endl;
      fPosition.Print();
      std::cout << "\t"
                << "momentum" << std::endl;
      fMomentum.Print();
      std::cout << "\t"
                << "polarization" << std::endl;
      fPolarization.Print();
   }

   Int_t fStepNumber = 0;
   Double_t fTrackLength = 0.;
   TLorentzVector fPosition;
   TLorentzVector fMomentum;
   TVector3 fPolarization;
   Double_t fWeight = 1.;
   UInt_t fGeoStateIndex = 0;
   Int_t fId = -1;
   Int_t fParentId = -1;
   Bool_t fIsOutside;
 
private:
   TMCParticleStatus(const TMCParticleStatus &);
   TMCParticleStatus &operator=(const TMCParticleStatus &);
 
   ClassDef(TMCParticleStatus, 1)
};
 
#endif /* ROOT_TMCParticleStatus */