Bash parallelization of MicrOmegas

[xiekeping]

Hi Florian,

I have question about the possibility of parallelization of MicrOmegas with simple BASH script. We need to run the import different spectrum and rotation matrix files generated by SPheno into MicrOmegas. As default, the MicrOmegas code generated by SARAH can only read the file name as SPheno.spc.MODEL. How can I modify the code, use an argument to specify the spectrum file when run MicrOmegas, such as CalcOmega_with_DDetection_MOv4.3?

Did you think about this possibility before? With Bash, I have already parallelized the SPheno running very easily.

Best regards,
Keping

[FStaub]

Hi,

the problem is that MO can not really parallelised because it works with absolute paths at some places. You need to create several MO projects for your model in different sub-directories.

Cheers,
Florian

[xiekeping]

Dear Florian,

Thank your for your reply. Just one follow up question. Did you consider the latest MadDM for DM study? I generate the UFO files for MadGraph, which is workable. But when I switch to the MadDM, I met some incompatibilities between SPheno output as param_cart.dat and MadDM reader as follows. I attached my UFO files, which includes the param_card.dat as well.

Traceback (most recent call last):
File "./maddm.py", line 23, in <module>
dm.FindCoannParticles()
File "/home/kepingx/MG5_aMC/MG5_aMC_v2_6_1/maddm/darkmatter.py", line 443, in FindCoannParticles
self.ChangeParamCard(prompts)
File "/home/kepingx/MG5_aMC/MG5_aMC_v2_6_1/maddm/darkmatter.py", line 535, in ChangeParamCard
self._fullmodel.set_parameters_and_couplings(param_card=self._paramcard)
File "/home/kepingx/MG5_aMC/MG5_aMC_v2_6_1/models/model_reader.py", line 162, in set_parameters_and_couplings
raise MadGraph5Error, msg
madgraph.MadGraph5Error: Invalid restriction card (not same block)
set(['ye', 'umix', 'tu', 'snumix', 'uulmix', 'sphenolowenergy', 'fwcoef', 'msu2', 'loopmsoft', 'gauge', 'nmix', 'lsp', 'msd2', 'chargemix', 'phases', 'sminputs', 'msl2', 'decay', 'higgslhc14', 'minpar', 'selmix', 'flavorkitqfv', 'higgslhc13', 'doublechargemix', 'higgslhc8', 'effhiggscouplings', 'td', 'te', 'higgslhc7', 'usqmix', 'imfwcoef', 'uelmix', 'spheno', 'imnmix', 'msoft', 'msq2', 'higgsfcc100', 'dsqmix', 'yd', 'modsel', 'treemsoft', 'yu', 'mse2', 'pseudoscalarmix', 'none', 'gaugegut', 'uurmix', 'udlmix', 'extpar', 'vmix', 'udrmix', 'imphases', 'mass', 'uermix', 'flavorkitlfv', 'hmix', 'scalarmix']) != set(['imyu', 'umix', 'imudlmix', 'snumix', 'imuelmix', 'imuurmix', 'msoft', 'doublechargemix', 'imvmix', 'imyd', 'imte', 'imdoublechargemix', 'uulmix', 'effhiggscouplings', 'chargemix', 'imnmix', 'phases', 'udrmix', 'decay', 'tu', 'selmix', 'imsnumix', 'imuulmix', 'immsoft', 'td', 'te', 'imtd', 'usqmix', 'imselmix', 'uelmix', 'dsqmix', 'imudrmix', 'ye', 'yd', 'imhmix', 'imuermix', 'imtu', 'sminputs', 'yu', 'pseudoscalarmix', 'imye', 'none', 'uurmix', 'udlmix', 'imumix', 'imusqmix', 'vmix', 'nmix', 'imphases', 'mass', 'imdsqmix', 'uermix', 'hmix', 'scalarmix']).
Missing block: imyu,imye,imtd,imte,imselmix,imuelmix,imuurmix,imudrmix,imudlmix,imsnumix,imumix,imhmix,imusqmix,imvmix,imyd,imtu,imdoublechargemix,imdsqmix,imuulmix,immsoft,imuermix
Unknown block : sphenolowenergy,fwcoef,msu2,loopmsoft,gauge,lsp,msd2,msl2,higgslhc14,minpar,flavorkitqfv,higgslhc13,higgslhc8,higgslhc7,imfwcoef,spheno,higgsfcc100,modsel,treemsoft,mse2,gaugegut,extpar,msq2,flavorkitlfv

Best regards,
Keping

[FStaub]

Hi,

I never tried MadDM. However, since it is based on MG, you need most likely to set flag

Code: Select all

 78 1               # Output for MadGraph (writes also vanishing blocks)   
 

in the SPheno input to get all kind of useless information in the output.

Cheers,
Florian