pgs muon cleaning script (customize)

Asked by Pablo

Hi,

I would like to run pgs without the muon cleaning script already included in madgraph
I guess it should be possible changing something in
   /Template/bin/internal/run_pgs

where it says

   $pgsdir/pgs --stdhep pythia_events.hep --nev 100000 --L0 --detector ../Cards/pgs_card.dat pgs_uncleaned_events.lhco \
    && $pgsdir/clean_output -muon pgs_uncleaned_events.lhco pgs_events.lhco && \
   touch pgs.done

I tried to comment the seond line, but then lhco file has no events..

I guess that I should be able also to change this in the pgs folder, where the cleaning file is, however I would like to be able to run pgs with no muon cleaning and then apply my own cuts in the ExRoot Analysis

Is it possible to do this just changin the run_pgs file?

Thanks in advance

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Johan Alwall
Solved:
Last query:
Last reply:
Revision history for this message
Johan Alwall (johan-alwall) said :
#1

Hello Pablo,

I guess what you want in the run_pgs file is

   $pgsdir/pgs --stdhep pythia_events.hep --nev 100000 --L0 --detector ../Cards/pgs_card.dat pgs_events.lhco \
   touch pgs.done

(note the change in the output file name). Was this what you did?

All the best,
Johan

Revision history for this message
Pablo (pablosp1988) said :
#2

That was my guess ... however as you mentioned I didn't remove the "uncleaned" part in the original

Original:
    $pgsdir/pgs --stdhep pythia_events.hep --nev 100000 --L0 --detector ../Cards/pgs_card.dat pgs_uncleaned_events.lhco \
     && $pgsdir/clean_output -muon pgs_uncleaned_events.lhco pgs_events.lhco && \
     touch pgs.done

To: $pgsdir/pgs --stdhep pythia_events.hep --nev 100000 --L0 --detector ../Cards/pgs_card.dat pgs_events.lhco \
   touch pgs.done

However now, wether I used uncleaned or not I get the error Fail to create LHCO events, here is part of the output:

    Running hep2lhe
    pythia_events.lhe file generated, with simple jet clustering
    Warning! Never use this file for detector studies!
    Creating Pythia LHE Root File
    ** Calculating number of events to process. Please wait...
    ** Input file contains 1000 events
    [#########################] (100.00%) : 1000/1000 entries processed
    ** Exiting...
    finish
    prepare PGS run
    Running PGS
    Fail to create LHCO events

Thanks for your help,
Pablo

Revision history for this message
Best Johan Alwall (johan-alwall) said :
#3

I'm sorry, you need
$pgsdir/pgs --stdhep pythia_events.hep --nev 100000 --L0 --detector ../Cards/pgs_card.dat pgs_events.lhco \
 && touch pgs.done
(note the "&&" before "touch pgs.done").

All the best,
Johan

Revision history for this message
Pablo (pablosp1988) said :
#4

That codes works in the sense it gives no error.
However, when I get into the *lhco pgs file , the lhco events are empty. I mena,it got the headres and so, but the place where the events shoudl appear, it's empty. Of course, the root file it's also empty.

Thanks for the support,
Pablo

Revision history for this message
Johan Alwall (johan-alwall) said :
#5

Do you get events when you run with the unmodified version?

Johan

Revision history for this message
Pablo (pablosp1988) said :
#6

Thanks Johan Alwall, that solved my question.

Revision history for this message
Pablo (pablosp1988) said :
#7

Sorry Johan, I forgot again to replace the uncleaned so your code

$pgsdir/pgs --stdhep pythia_events.hep --nev 100000 --L0 --detector ../Cards/pgs_card.dat pgs_events.lhco \
 && touch pgs.done

is working perfectly

Thanks a lot, the best for you,
Pablo