[Rivet] Rivet Query

[Andy Buckley]

Hi Debarati,

That is perhaps the easiest way. Apart from a few changes in e.g. the
headers which it includes, I think the new Hemispheres projection code
should be fully backward compatible with the latest stable release.

Here are the new files:

https://rivet.hepforge.org/trac/browser/include/Rivet/Projections/Hemispheres.hh
https://rivet.hepforge.org/trac/browser/src/Projections/Hemispheres.cc

The alternative would be to build the full development version but I
think that is more likely to cause problems. Hopefully the copy &
pasting will work well enough for you!

Andy


On 26/08/14 14:57, Debarati Roy wrote:
> Hi Andy,
>              Many thanks. 
> For the  test I need a feedback. 
> Should I try the new projection with the existing Rivet version by modifying the existing Rivet projection according to  the new one or I need to do something else?
> 
> Regards,
> Debarati
>     
> ________________________________________
> From: Andy Buckley
> Sent: 26 August 2014 13:17
> To: Debarati Roy; rivet at projects.hepforge.org
> Cc: Albert Knutsson; Maxime Gouzevitch; Lars Sonnenschein
> Subject: Re: [Rivet] Rivet Query
> 
> Hi Debarati,
> 
> I have added several calc() methods including one that accepts a Jets
> list as input, and a M2sum() method even though this is trivial. I will
> push it to our development branch later today if you want to test it or
> just take a look at the new interface. Hope it works for you. It will
> appear in the next Rivet release which should not be far away.
> 
> Andy
> 
> 
> On 18/08/14 05:37, Debarati Roy wrote:
>> Hi Andy,
>>             That is nice.
>>
>>
>> Thanks,
>> Debarati
>> ________________________________________
>> From: Andy Buckley
>> Sent: 18 August 2014 00:23
>> To: Debarati Roy; rivet at projects.hepforge.org
>> Cc: Albert Knutsson; Maxime Gouzevitch; Lars Sonnenschein
>> Subject: Re: [Rivet] Rivet Query
>>
>> Hi Debarati,
>>
>> I will see if I can provide a calc() method on Hemispheres which will
>> accept a vector of FourMomentum. You could that way assemble whatever
>> you wanted as the input, with z components as original or zeroed. I'll
>> let you know how that goes.
>>
>> I'll also look into providing an M2sum method.
>>
>> Andy
>>
>>
>> On 15/08/14 11:01, Debarati Roy wrote:
>>> Hi Andy, Thanks for your feedback. First of all I would like to give
>>> you some more information regarding these hadronic event shape
>>> calculations applied in this analysis.
>>>
>>> In this analysis for the calculation of the five hadronic event
>>> shapes always the input variables we used are jet Px, Py, Pz and E.
>>> In some of them we do not require that Pz component like the thrust,
>>> total jet broadening.
>>>
>>> Now regarding your comments, For the total jet broadening I think
>>> your calc method can be useful since we have used only transverse
>>> component of momentum along with jet eta and phi. But for the total
>>> jet mass the concept is very similar as  used in Hemispheres.cc for
>>> mass2With and mass2Against. The difference exists in the input(always
>>> jet's four momentum) and then total jet mass is the sum of  mass2With
>>> and mass2Against(this one I called as M2sum() in my previous mail).
>>> For the total transverse jet mass I think the same method can work
>>> only again in transverse plane(Z component set to 0).  So over all I
>>> would say the jet hemisphere projection will cover all:
>>> broadening(setZ=0), total jet mass, total transverse jet mass(setZ=0)
>>> with jet's four momentum as inputs. If the jet hemisphere  you do not
>>> think should be added as a projection then I can try with the
>>> analysis code.
>>>
>>> Please let us know your view.
>>>
>>> Regards, Debarati ________________________________________ From: Andy
>>> Buckley Sent: 15 August 2014 01:09 To: Debarati Roy;
>>> rivet at projects.hepforge.org Cc: Albert Knutsson; Maxime Gouzevitch;
>>> Lars Sonnenschein Subject: Re: [Rivet] Rivet Query
>>>
>>> Hi Debarati,
>>>
>>> Those hemispheres are an e+e- observable, so I'm a bit surprised that
>>> a non-transverse decomposition work at a hadron collider... but then
>>> I'm not a hadronic event shapes expert. Are you using the jets to
>>> build a "3D" hemispheres observable, or do you want to set the jet
>>> z-components to zero and then run the hemispheres algorithm?
>>>
>>> If the latter, I guess we could expose a "calc()" method on the
>>> projection so it can be used directly rather than via applyProjection
>>> -- this is sometimes done with the FastJets projection, for example.
>>>
>>> If the former, then maybe we should supply a JetHemispheres
>>> projection instead. It might be easier to just implement in the
>>> analysis code, though, since projections are only really necessary
>>> for code that's likely to appear in several analyses and which would
>>> benefit from caching.
>>>
>>> I'd be happy to add a M2sum() method if you can define what it should
>>> do.
>>>
>>> Cheers, Andy
>>>
>>>
>>> On 14/08/14 16:10, Debarati Roy wrote:
>>>> Hi Experts, I need some feedback regarding an analysis which we
>>>> are planning to introduce in Rivet. It is a CMS analysis consists
>>>> of the study of hadronic event shapes (arxiv link :
>>>> http://arxiv.org/abs/arXiv:1407.2856). In this analysis five event
>>>> shapes are studied, thrust, total jet broadening, total jet mass,
>>>> total transverse jet mass, jet resolution parameter where each of
>>>> the event shape variable is calculated by giving the input of final
>>>> state jets momenta . I found that Rivet already has the thrust
>>>> projection where jets momenta of an event can be given as input to
>>>> calculate thrust.  But for other three variables (total jet
>>>> broadening, total jet mass, total transverse jet mass) I observed
>>>> that you have a similar but not identical projection named,
>>>> https://rivet.hepforge.org/trac/browser/src/Projections/Hemispheres.cc
>>>>
>>>>
>> In this projection you are taking the final state particles to calculate
>>>> broadening and mass if I am not wrong whereas in this analysis the
>>>> final state jets momenta are used. So is it possible to modify the
>>>> existing projection to be more flexible so that we can run it over
>>>> the final state jets too? Or is it possible that a different
>>>> projection with a different name can be added? Also here
>>>> https://rivet.hepforge.org/trac/browser/include/Rivet/Projections/Hemispheres.hh
>>>>
>>>>
>> I could see that there exists different methods  M2high(), M2low(),
>>>> M2diff(). Is it possible to introduce M2sum() as per this analysis
>>>> requirement?
>>>>
>>>> Please let us know.
>>>>
>>>> Thanks, Debarati
>>>>
>>>>
>>>> _______________________________________________ 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
>>>
>>
>>
>> --
>> Dr Andy Buckley, Royal Society University Research Fellow
>> Particle Physics Expt Group, University of Glasgow / PH Dept, CERN
>>
>> --
>> Dr Andy Buckley, Royal Society University Research Fellow
>> Particle Physics Expt Group, University of Glasgow / PH Dept, CERN
>>
> 
> 
> --
> Dr Andy Buckley, Royal Society University Research Fellow
> Particle Physics Expt Group, University of Glasgow / PH Dept, CERN
> 


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