Add new component MDIReader

[aciarma]

Added MDIReader, to use as input ascii files containing lists of particles.
Formats supported at the moment:

• GuineaPig++ pairs production
• XTrack beam losses

[aciarma]

Not ready yet, when running a test job for MDIReader I get the following error:

WARNING: cannot load libk4Gen.so for factory MDIReader
WARNING: libHepPDT.so: cannot open shared object file: No such file or directory
AlgorithmManager    ERROR Algorithm of type MDIReader is unknown (No factory available).

[vvolkl]

Use gaudi logging instead of std::cout

Suggested change

	std::cout<<"Fine di MDIreader"<<std::endl;
	debug() << "Finalization of MDIREADER" << endmsg;

[aciarma]

Thanks, I changed also the other std::cout to debug()

[vvolkl]

You should use HepMC3 here (see https://github.com/HEP-FCC/k4Gen/blob/main/k4Gen/src/components/HepMCFileReader.h)

[vvolkl]

Please also add an example options file and test, like here: https://github.com/HEP-FCC/k4Gen/blob/main/k4Gen/CMakeLists.txt#L49

[aciarma]

I think it should be ok now. Please let me know if this request can be accepted or if I have to do other modifications.

[vvolkl]

Yes, looks good. I'd like to split the CMake changes from this PR and properly add the RPATH changes to this repo, I'll take care of that and merge this afterwards.

[vvolkl]

Hi @aciarma,

I took the liberty of rebasing this onto the RPATH changes in master and added a quick fix. But I have a couple of questions more:

• Can you remove the geant algorithms from the test? This will complicate the dependencies, here it is enough just to run MDIReader
• As far as I see, the MDIReader can only ever process one event (one event per file, right?) Maybe allow reading a list of events so more than one can be processed during the same job?
• I'm surprised you can set mdi_converter.LongitudinalCut = 0 in the job options without declaring this and other floats as Gaudi::Property, can you check this really works?

[aciarma]

Removed the Geant4 call from the example and verified that the job options work correctly (even if not declared as Gaudi::Property).

[vvolkl]

Looks good, thank you very much!