[Rivet] Bug/change in RIVET IdentifiedFinalState()?

Marco Delmastro Marco.Delmastro at cern.ch
Wed 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