[Rivet] Bug in CMS_2013_I1209721

[Andy Buckley]

Thanks Giulio,

I've merged that into our release branch. We should release a Rivet
2.1.2 fairly soon to address this and other issues.

Cheers,
Andy


On 02/05/14 17:22, Piergiulio Lenzi wrote:
> Dear David,
>     thanks for your prompt reply.
> Looking more carefully into CMS_2013_I1209721 I found two other problems.
> The first is the histogram booking, which should be:
> 
> // Book histograms from data
>       for (size_t i = 0; i < 2; ++i) {
>         _histDeltaPhiZJ1_1[i]  = bookHisto1D(1+i*9, 1, 1);
>         _histDeltaPhiZJ1_2[i]  = bookHisto1D(2+i*9, 1, 1);
>         _histDeltaPhiZJ1_3[i]  = bookHisto1D(4+i*9, 1, 1);
>         _histDeltaPhiZJ2_3[i]  = bookHisto1D(5+i*9, 1, 1);
>         _histDeltaPhiZJ3_3[i]  = bookHisto1D(3+i*9, 1, 1);
>         _histDeltaPhiJ1J2_3[i] = bookHisto1D(6+i*9, 1, 1);
>         _histDeltaPhiJ1J3_3[i] = bookHisto1D(7+i*9, 1, 1);
>         _histDeltaPhiJ2J3_3[i] = bookHisto1D(8+i*9, 1, 1);
>         _histTransvThrust[i]   = bookHisto1D(9+i*9, 1, 1);
>       }
> 
> currently the "+i*9" is missing and the "boosted regime" histograms are
> never booked (or booked incorrectly with the non-boosted regime binning)
> 
> Then the part in which the momenta for the transverse thrust calculation
> are prepared is missing a piece
> 
> // Now compute the thrust
>       // Collect Z and jets transverse momenta to calculate transverse
> thrust
>       vector<Vector3> momenta;
>       momenta.clear();
>       Vector3 mom = z[0].momentum().p();
>       mom.setZ(0);
>       momenta.push_back(mom);
> 
>       for (size_t i = 0; i < cleanedJets.size(); ++i) {
>         Vector3 mj = cleanedJets[i]->momentum().vector3();
>         mj.setZ(0);
>         momenta.push_back(mj);
>       }
> // THE PART BELOW IS MISSING
>       if (momenta.size() <= 2){
>         // We need to use a ghost so that Thrust.calc() doesn't return 1.
>         momenta.push_back(Vector3(0.0000001,0.0000001,0.));
>       }
> 
> This part is needed for a correct thust calculation in Z+1 jet events
> (momenta.size()=2)
> 
> I'm also attaching the fixed analysis.
> 
> Cheers
> Giulio
> 
> On 01/05/2014 17:54, David Grellscheid wrote:
>> Hi Giulio,
>>
>> thanks for spotting this! It's fixed in the repository now. The patch
>> will be released in the next version and is also available at
>>
>>   https://rivet.hepforge.org/hg/rivet/raw-rev/94600cbc5503
>>
>> See you,
>>
>>    David
>>
>>
>> On 01/05/14 16:03, Piergiulio Lenzi wrote:
>>> Dear Rivet developers,
>>>       I noticed a but in CMS_2013_I1209721,
>>> the Z projections read:
>>>         // Z finders for electrons and muons
>>>         const ZFinder zfe(fs, -2.4, 2.4, 20*GeV, PID::ELECTRON, 71*GeV,
>>> 111*GeV);
>>>         const ZFinder zfm(fs, -2.4, 2.4, 20*GeV, PID::ELECTRON, 71*GeV,
>>> 111*GeV);
>>> while it should be
>>>
>>>         // Z finders for electrons and muons
>>>         const ZFinder zfe(fs, -2.4, 2.4, 20*GeV, PID::ELECTRON, 71*GeV,
>>> 111*GeV);
>>>         const ZFinder zfm(fs, -2.4, 2.4, 20*GeV, PID::MUON, 71*GeV,
>>> 111*GeV);
>>>
>>> Cheers
>>> Giulio
>>> _______________________________________________
>>> Rivet mailing list
>>> Rivet at projects.hepforge.org
>>> https://www.hepforge.org/lists/listinfo/rivet
>>>
> 
> 
> 
> _______________________________________________
> Rivet mailing list
> Rivet at projects.hepforge.org
> https://www.hepforge.org/lists/listinfo/rivet
> 


-- 
Dr Andy Buckley, Royal Society University Research Fellow
Particle Physics Expt Group, University of Glasgow / PH Dept, CERN