g4Config.C

Go to the documentation of this file.

/********************************************************************************
 *    Copyright (C) 2014 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH    *
 *                                                                              *
 *              This software is distributed under the terms of the             *
 *              GNU Lesser General Public Licence (LGPL) version 3,             *
 *                  copied verbatim in the file "LICENSE"                       *
 ********************************************************************************/
// Configuration macro for Geant4 VirtualMC
void Config()
{

    Bool_t mtMode = FairRunSim::Instance()->IsMT();
    Bool_t specialStacking = false;
    FairFastSimRunConfiguration* runConfiguration = new FairFastSimRunConfiguration(
        "geomRoot", "QGSP_BERT_EMV", "stepLimiter+specialCuts+specialControls", specialStacking, mtMode);
    // TG4RunConfiguration* runConfiguration
    //     = new TG4RunConfiguration("geomRoot", "QGSP_BERT_EMV", "stepLimiter+specialCuts+specialControls",
    //                               specialStacking, mtMode);

    TGeant4* geant4 = new TGeant4("TGeant4", "The Geant4 Monte Carlo", runConfiguration);
    cout << "Geant4 has been created." << endl;

    FairStack* stack = new FairStack(1000);
    stack->StoreSecondaries(kTRUE);
    stack->SetMinPoints(0);
    geant4->SetStack(stack);

    if (FairRunSim::Instance()->IsExtDecayer()) {
        TVirtualMCDecayer* decayer = TVirtualMC::GetMC()->GetDecayer();
        //      TVirtualMCDecayer* decayer = TPythia6Decayer::Instance();
        geant4->SetExternalDecayer(decayer);
    }


    TString configm(gSystem->Getenv("VMCWORKDIR"));
    TString configm1 = configm + "/common/gconfig/g4config.in";
    cout << " -I g4Config() using g4conf  macro: " << configm1 << endl;

    // set geant4 specific stuff
    geant4->SetMaxNStep(10000);   // default is 30000
    geant4->ProcessGeantMacro(configm1.Data());
}