[Rivet] ATLAS prompt photon analysis

[Mike Hance]

Hi Frank, All,

>> I think you can add the analysis to Rivet's SVN now - I assume (since
>> it's SVN) that if we find some problems later that we can change the
>> code there. But things look good from my side, and it seems as though
>> they're good on your side too.
> 
> Yes, of course it's still modifiable.

Giovanni found a silly mistake on my part when extracting the uncertainties associated with our measurement (thanks 
Giovanni!) and I found another problem with my flat->aida conversion (I was listing the positive uncertainty first in 
the .dat file, but to make the .aida file correct the negative uncertainty should be first).  Those problems have been 
fixed in the attached tarball (same contents as before, only the .raw, .dat, and .aida files have changed).

So I think that with these final changes the data should be ready to insert into the HEPDATA database.  How can I/we go 
about doing that?  I think all the necessary information should be in the tarball, but if it's not, I can provide 
whatever is required.  If I should bother one person in particular to get the data inserted, please just let me know who 
that is (and I'll stop bothering all of you).

Thanks again for the help!

-Mike



> I have one more question: Is there any reason why you number the third 
> histogram as d01-x01-y04 instead of d01-x01-y03? I suppose on HepData 
> it's going to appear as d01-x01-y03 and since Rivet is supposed to be 
> compatible with the files exported from HepData that's also what it 
> would have to be in Rivet. If there is no reason for this, could you 
> make this small adaptation?
> I guess it also affects the following line in the code?
>   if(fabs(_eta_bins[i] - 1.37) < .0001) continue;
> 
> Cheers,
> Frank
> 
>>
>>
>>
>> Frank Siegert wrote:
>>> Hi Mike,
>>>
>>> Apart from the nb/pb issue the analysis seems to work well. Take a
>>> look at
>>> http://fsiegert.web.cern.ch/fsiegert/atlas/yjets/
>>> where I'm comparing it to the Sherpa 1.2.3 gamma+jets validation
>>> sample that we talked about in the MC validation meeting today:
>>> http://indico.cern.ch/conferenceDisplay.py?confId=127609
>>> Ignore the first one or two bins, since the sample has a pT_gamma>15
>>> GeV cut (on the ME level) which affects those. But apart from that it
>>> looks pretty good.
>>>
>>> If you want, I could add the analysis to Rivet's SVN now, or did you
>>> want to do any more checks?
>>>
>>> Cheers,
>>> Frank
>>>
>>> On 17/02/11 21:22, Mike Hance wrote:
>>>> Hi Frank,
>>>>
>>>> You're right, we should just put everything in pb. I took the numbers
>>>> from the tables in the appendix of the paper, which were in nb, and
>>>> mentally normalized away the extra 10^3 terms I saw in the rivet 
>>>> output.
>>>> Thanks for spotting that!
>>>>
>>>> I've attached a new tarball with that fixed, and with the
>>>> files/functions/plots/etc renamed to have the proper SPIRES ID (thanks
>>>> Jon!). I'm sure there are other mistakes though, so feel free to point
>>>> them out. In the mean time, I hope Giovanni (and others) will be 
>>>> able to
>>>> quickly spot-check the numbers in these files to make sure we've given
>>>> you the right numbers.
>>>>
>>>> Cheers,
>>>>
>>>> -Mike
>>>>
>>>>
>>>>
>>>> Frank Siegert wrote:
>>>>> Hi Mike,
>>>>>
>>>>> Thanks a lot and don't worry about it taking longer, it's a
>>>>> complicated analysis after all!
>>>>>
>>>>> One issue that I noticed without having done detailed testing yet: The
>>>>> .aida reference data seems to be in nb while Rivet's crossSection() is
>>>>> by default in pb and thus that's what the MC histos are being filled
>>>>> the way the code is now. Since the plot in the paper is also in pb I
>>>>> would suggest adapting the reference .aida file, or what do you think?
>>>>>
>>>>> Cheers,
>>>>> Frank
>>>>>
>>>>> On 17/02/11 17:05, Mike Hance wrote:
>>>>>> Hi All,
>>>>>>
>>>>>> I've attached a tarball that (I think) contains everything needed
>>>>>> to put
>>>>>> the ATLAS prompt photon results into HEPDATA, as well as a
>>>>>> corresponding
>>>>>> RIVET routine that seems like it might be doing the right thing.
>>>>>>
>>>>>> One thing I should note up front is that, as far as I can tell, the
>>>>>> paper doesn't have a proper SPIRES ID. I suspect this will change 
>>>>>> when
>>>>>> PRD finally publishes the paper (it's been accepted for
>>>>>> publication, but
>>>>>> I don't think it's actually made it into an issue of the journal 
>>>>>> yet),
>>>>>> so in the mean time, I've called the analysis
>>>>>> "ATLAS_2010_S1234567". My
>>>>>> attempts at calling the analysis something like
>>>>>> "ATLAS_2010_PROMPTPHOTONS" were rebuffed by RIVET, so I'm not sure of
>>>>>> the right thing to do in this case. If you want me to relabel it as
>>>>>> something else, just let me know; otherwise, feel free to call it
>>>>>> whatever you like.
>>>>>>
>>>>>> The contents of the tarball are:
>>>>>>
>>>>>> ATLAS_2010_S1234567.info - The .info file that describes the 
>>>>>> analysis.
>>>>>>
>>>>>> ATLAS_2010_S1234567.raw - A flat file that breaks down the different
>>>>>> (E_T, eta) bins and the corresponding measurements, uncertainties, 
>>>>>> and
>>>>>> so forth.
>>>>>>
>>>>>> ATLAS_2010_S1234567.dat - This is the 'final' flat .dat file
>>>>>> corresponding to the ATLAS measurements, for comparison with the 
>>>>>> RIVET
>>>>>> output.
>>>>>>
>>>>>> ATLAS_2010_S1234567.aida - The .aida file produced from the .dat file
>>>>>> above.
>>>>>>
>>>>>> ATLAS_2010_S1234567.plot - The .plot file for the three histograms
>>>>>> associated with this analysis.
>>>>>>
>>>>>> ATLAS_2010_S1234567.cc - The RIVET routine.
>>>>>>
>>>>>>
>>>>>> I've tested this running PYTHIA standalone, but the cross sections 
>>>>>> for
>>>>>> prompt photon production are low enough that the histogram never
>>>>>> really
>>>>>> gets populated. I also ran the routine using the ATHENA RIVET
>>>>>> interface
>>>>>> against some ATLAS photon+jet (EVGEN) samples, and the results seemed
>>>>>> reasonable, but I wouldn't be surprised if there are still some bugs
>>>>>> lurking in there. I'll try to test this stuff against some filtered
>>>>>> dijet samples next week.
>>>>>>
>>>>>> In any event, please feel free to let me know if there's anything
>>>>>> unclear/wrong with what I've done, and I'll try to fix it ASAP.
>>>>>>
>>>>>> Finally, thanks to everyone who helped (and bugged) us to get this
>>>>>> done
>>>>>> - I apologize it took us so long. Next time should be much faster, 
>>>>>> now
>>>>>> that we know what we're doing.
>>>>>>
>>>>>> Thanks!
>>>>>>
>>>>>> -Mike
>>>>>>
>>>>>>
>>>>>>
>>>>>> _______________________________________________
>>>>>> Rivet mailing list
>>>>>> Rivet at projects.hepforge.org
>>>>>> http://www.hepforge.org/lists/listinfo/rivet
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ATLAS_prompt_photon_analysis.tar
Type: application/x-tar
Size: 30720 bytes
Desc: not available
URL: <http://www.hepforge.org/lists-archive/rivet/attachments/20110223/5f1df0ea/attachment.tar>