|
[Rivet] Bug/change in RIVET IdentifiedFinalState()?Marco Delmastro Marco.Delmastro at cern.chWed Jul 8 10:05:03 BST 2015
Hello, You'll find the analysis attached here. As I mentioned, I'm using rivet 2.2.0 via the ATLAS Athena framework, M. -------------- next part -------------- A non-text attachment was scrubbed... Name: ATLAS_2012_DIPHOTONS_bug.cc Type: application/octet-stream Size: 12163 bytes Desc: not available URL: <https://www.hepforge.org/lists-archive/rivet/attachments/20150708/ef2b0777/attachment.obj> -------------- next part -------------- On 08 Jul 2015, at 11:02, David Grellscheid <david.grellscheid at durham.ac.uk> wrote: > Hi Marco, > > can you please send a minimal analysis and a short HepMC file that shows > the problem? Also, which Rivet version are you using? > > Thanks, > > David > > > On 08/07/15 09:27, Marco Delmastro wrote: >> Hello, >> >> Andy Buckley suggested that I report this bug (?) I encountered when using Rivet within the ATLAS Athena framework. >> >> I'm writing the RIVET routines for the ATLAS 8 TeV inclusive diphoton and diphoton + jet analyses. In the first tests I made, running on a Pythia diphoton samples, I noticed something strange: only photons with eta > 1 were used, and the initial pT selection was not applied. I tracked down the problem to this method: >> >> IdentifiedFinalState photonfs(Cuts::abseta < 2.37 && Cuts::pT > 30.*GeV); // 30 GeV cut on subleading photon >> >> Apparently the 30 GeV cut on the photon pT is not applied, and or some reason only objects with eta > 1 are kept. From some conversation on the ATLAS MC Generator mailing list, it seems that the problem might be related to the handling of the && operator. If I change the line to: >> >> IdentifiedFinalState photonfs(-2.37,2.37,30.*GeV); // etamin, etamax, ptmin >> >> everything seems to work fine. I'm using rivet 2.2.0 from Athena 19.2.1: >> >> [mdelmast at lxplus0157 Rivet]# which rivet >> /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.1/AtlasSimulation/19.2.1/InstallArea/share/bin/rivet >> >> [mdelmast at lxplus0157 Rivet]# rivet >> Rivet 2.2.0 running on machine lxplus0157.cern.ch (x86_64) >> >> Thanks for your help! >> >> Ciao, M. >> >> >> >> _______________________________________________ >> Rivet mailing list >> Rivet at projects.hepforge.org >> https://www.hepforge.org/lists/listinfo/rivet >> -- ------------------------------------------------- Marco DELMASTRO | phone @ LAPP: +33 4 50 09 17 85 CNRS/LAPP ATLAS | phone @ CERN: +41 76 4 87 86 90 ------------------------------------------------- -------------- next part -------------- A non-text attachment was scrubbed... Name: smime.p7s Type: application/pkcs7-signature Size: 3034 bytes Desc: not available URL: <https://www.hepforge.org/lists-archive/rivet/attachments/20150708/ef2b0777/attachment.bin>
More information about the Rivet mailing list |