GeneratorTGenerator.cxx

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
#include "SimulationDataFormat/MCGenProperties.h"
#include "SimulationDataFormat/ParticleStatus.h"
#include "Generators/GeneratorTGenerator.h"
#include "Generators/GeneratorExternalParam.h"
#include <fairlogger/Logger.h>
#include "FairPrimaryGenerator.h"
#include "TGenerator.h"
#include "TClonesArray.h"
#include "TParticle.h"
#include "TMCProcess.h" // for VMC Particle Production Process
 
namespace o2
{
namespace eventgen
{
 
/*****************************************************************/
/*****************************************************************/
 
GeneratorTGenerator::GeneratorTGenerator() : Generator("ALICEo2", "ALICEo2 TGenerator Generator"),
                                             mTGenerator(nullptr),
                                             mCloneParticles(nullptr)
{
  mInterface = reinterpret_cast<void*>(mTGenerator);
  mInterfaceName = "tgenerator";
  mCloneParticles = new TClonesArray("TParticle");
  mCloneParticles->SetOwner(kTRUE);
}
 
/*****************************************************************/
 
GeneratorTGenerator::GeneratorTGenerator(const Char_t* name, const Char_t* title) : Generator(name, title),
                                                                                    mTGenerator(nullptr),
                                                                                    mCloneParticles(nullptr)
{
  mInterface = reinterpret_cast<void*>(mTGenerator);
  mInterfaceName = "tgenerator";
  mCloneParticles = new TClonesArray("TParticle");
  mCloneParticles->SetOwner(kTRUE);
}
 
/*****************************************************************/
 
GeneratorTGenerator::~GeneratorTGenerator()
{
  delete mCloneParticles;
}
 
/*****************************************************************/
 
Bool_t
  GeneratorTGenerator::generateEvent()
{
  mTGenerator->GenerateEvent();
 
  return kTRUE;
}
 
/*****************************************************************/
 
Bool_t
  GeneratorTGenerator::importParticles()
{
  mTGenerator->ImportParticles(mCloneParticles, "All");
  auto nparticles = mCloneParticles->GetEntries();
  for (Int_t iparticle = 0; iparticle < nparticles; iparticle++) {
    auto particle = (TParticle*)mCloneParticles->At(iparticle);
    // Some particle comes in from somewhere.
    // To support the user and in case the status code is found to be not encoded, do it here,
    // assuming that the current code is the HepMC status (like it used to be before).
    auto statusCode = particle->GetStatusCode();
    if (!mcgenstatus::isEncoded(statusCode)) {
      particle->SetStatusCode(mcgenstatus::MCGenStatusEncoding(statusCode, 0).fullEncoding);
    }
    mParticles.push_back(*particle);
    // Set the transport bit according to the HepMC status code as it always used to be
    mParticles.back().SetBit(ParticleStatus::kToBeDone, mcgenstatus::getHepMCStatusCode(particle->GetStatusCode()) == 1);
 
    if (GeneratorExternalParam::Instance().markAllAsPrimary) {
      // The production process will be set to kPPrimary which is the expected default for generator-level particles in O2
      // not doing this might confuse transport and selection of particles to keep in AOD
      mParticles.back().SetUniqueID(kPPrimary);
    }
  }
 
  return kTRUE;
}
 
/*****************************************************************/
 
Bool_t
  GeneratorTGenerator::Init()
{
  Generator::Init();
 
  if (!mTGenerator) {
    LOG(fatal) << "No TGenerator inteface assigned" << std::endl;
    return kFALSE;
  }
 
  return kTRUE;
}
 
/*****************************************************************/
/*****************************************************************/
 
} /* namespace eventgen */
} /* namespace o2 */