Propagating Weight ID's from LHEF to HepMC

Hi,

This is technically a question for pythia8 authors. Since they wrote that interface.
But I do not think that they are good reason for not doing what you are doing.
The only important point is to keep the current convention for name associated to scale variation.

Now It would be nice if you can share what you did. I will then propose this for merging in the official version (or you can do the proposal yourself if you prefer).

Cheers,

Olivier

Sure, sharing is eazy. It's a bit faster for me if I just copy the relevant parts here as there aren't too many changes.

So one change here is to accomodate labels without spaces (and I also removed the transformation to uppercase)

// Function to convert aMC@NLO weight names to "HepMC convention".
string convertFromAMCATNLO( string input) {
  string output="";
  // Count number of blanks in weight name.
  int appearances = 0;
  for(int n = input.find(" ", 0); n != int(string::npos);
          n = input.find(" ", n)) {
    appearances++;
    n++;
  }
  // Cut string by position of blanks.
  vector <string> pieces;
  for(int i =0; i < appearances;++i) {
    int n = input.find(" ", 0);
    pieces.push_back(input.substr(0,n));
    input = input.substr(n+1,input.size());
  }
  // Now pieces contains the details of the weight name.
  for(int i = 4; i < int(pieces.size()); ++i) {
    // Make upper case.
    transform(pieces[i].begin(),pieces[i].end(),pieces[i].begin(), ::toupper);
    // Add piece to weight name, with "_" as separator.
    output+=pieces[i] + "_";
  }
  // S.D. For labels without spaces
  if(appearances==0){
   //transform(input.begin(), input.end(), input.begin(), ::toupper);
   output=input + "_";
  }

  return output;
}

The second change is within the #ifdef HEPMCHACK block

#ifdef HEPMC2HACK
           <....>

              // Get name.
              string name;
              for ( std::map<string,string>::const_iterator
                it_att = it_w->second.attributes.begin();
                it_att != it_w->second.attributes.end(); ++it_att)
                name += it_att->first + "=" + it_att->second + "_";

              // In aMC@NLO, the weight names should be extracted from the
              // weight tag contents, not its attributes.
              // S.D. id actually contains the weight name!
              if (it_w->second.attributes.size()==0)
                name = convertFromAMCATNLO(it_w->second.id);

             <....>
            } // Done looping through weights.
          } // Done looping through weight groups.
        } // Done looping through cuts.
      }

Above I simply have
name = convertFromAMCATNLO(it_w->second.id);
Instead of
name = convertFromAMCATNLO(it_w->second.contents);

--------------------------------------------------------------------------------------------------------

Thanks for mentioning the naming conventions, I had scale variations turned off for now.
These changes only work when use_syst = False in the run_card.dat.
When systematics are used the interface works as before. The conventions are in tact, but then the weight names are back to being cryptic again. So there needs to be some sort of separation between the two.

Regards,

Simonas

Hi again,

the interface was designed by Valentin & myself some years back. At that point, we agreed that it was reasonable to use the *content* of the weight tag, since for the conventional 1001-1009 scale variation weights, this content is more in line with the guidelines on how to name weights for HepMC (see https://arxiv.org/pdf/1405.1067.pdf, page 162). Maybe Valentin remembers better...

Now it would of course not be ideal to use "1001" as weight name for these conventional scale variations... but I do see the point that in some cases, using the weight id might be more human-readable. Just shows again that it would be very helpful to have a proper naming convention for weights in the HepMC file (something I've been pushing without success since 2013, it seems :D )

Note that in the last pythia 8.3 release, we opted for a scheme were we internally, in Pythia, translate 1001-1009 to the names that we should use, while keeping all the other weights with names determined from their id's, prefixed with a string "AUX_" (indicating that we did not know how to translate these ids internally). Also, all weights are automatically written by the HepMC-converter, so you'll never have to touch the weights explicitly. That might help for you as well.

So from my side:
- a complete change of the MG5_PY8 interface to use weight ids as name might be a bit dangerous
- maybe a translation could be added, and a flag that would allow switching from the current version to your update
- i would very much value feedback on Pythia 8.3003, which should have the behavior you expected - although the MG5_PY8 interface might have to be updated a little to comply with 8.3-syntax.

Cheers,
Stefan

Hi,

I ran the main44 example on my lhe file with Pythia-8303.
The weights are named just as they should be - the weights are named after their IDs in the lhe file and there is "AUX_" attached in front.
It would be great to have this or analogous functionality when running Pythia through MadGraph!

Regards,

Simonas