Using ROOT input with multithreading

Particles, Track, Event, Run and Biasing
1

Hi all,

I am trying to run a simulation where the properties of the primary particles are read from a tree in a ROOT file. This is done through a dedicated FileReader class, of which the PrimaryGeneratorAction class has a static member. Since the file reading is done only as the PrimaryGeneratorAction is initialised (a PrimaryGeneratorMessenger macro command is used to supply the file name), it is done locally on a worker thread. Using a mutex, this should not be an issue as only one thread is needed to fill the FileReader.

// PrimaryGeneratorAction.hh
class PrimaryGeneratorAction : public G4VUserPrimaryGeneratorAction {
    // Other stuff...
    public:
        void ImportBeamFromFile(G4String&);

    private:
        G4ParticleGun* fParticleGun;
        static FileReader* fFileReader;
};

//PrimaryGeneratorAction.cc
namespace { G4Mutex rootPrimGenMutex = G4MUTEX_INITIALIZER; }
FileReader* PrimaryGeneratorAction::fFileReader = 0;

// Other stuff...

void PrimaryGeneratorAction::ImportBeamFromFile(G4String& fname) {
    G4AutoLock lock(&rootPrimGenMutex);
    if(!fFileReader) fFileReader->TakeBeamFromFile(fname);
}

// FileReader.cc
void FileReader::TakeBeamFromFile(G4String& fname) {
    auto analysisReader = G4RootAnalysisReader::Instance();
    analysisReader->SetFileName(fname);
    
    // Reading of file happens here...
}

However, since the FileReader uses the G4RootAnalysisReader class, which derives from G4VBaseFileManager, the file name supplied is appended with “_{$THREAD_ID}” (member function SetFullFileName) and so an error is thrown that the file cannot be found.

Short of manually creating all of the variants of the ROOT file with each thread id, is there any way to stop the file name being changed? Any thoughts or ideas on this would be greatly appreciated!

Thanks,
KF