[Rivet] Rivet Query

[Debarati Roy]

Hi Andy,
           No problem. Thanks.

Regards,
Debarati
________________________________________
From: Andy Buckley
Sent: 03 September 2014 01:34
To: Debarati Roy; rivet at projects.hepforge.org
Cc: Albert Knutsson; Maxime Gouzevitch; Lars Sonnenschein
Subject: Re: [Rivet] Rivet Query

Hi Debarati,

Sorry, I don't have time to debug this. The merge isn't just a copy &
paste -- it will require fixing up the other changes to #includes etc.
which you can work out by looking at the diffs in the web interface.

If that's too much then you could either try to check out and build the
trunk version (we don't offer support for that, though!) or to wait for
the 2.2.0 release... which I hope to have out within two weeks.

Andy


On 02/09/14 19:24, Debarati Roy wrote:
> Hi Andy,
>             Thanks for the files.
> But it seems  that the copy & pasting are not working.
> I set up the stand alone Rivet version 2.1.2 in lxplus SLC6 by following the commands from
> https://rivet.hepforge.org/trac/wiki/GettingStarted
> Now I just tried to do the "make" here /Rivet-2.1.2/src/Projections with the two new files(keeping Hemispheres.cc in the same path and Hemispheres.hh in path Rivet-2.1.2/include/Rivet/Projections).
> The code just crashed. The corresponding error file I have kept here, /afs/cern.ch/user/d/deroy/public/RivetErr.txt
>
> Can you please take a look at the txt file and suggest where I am wrong?
>
> Also if I do the "make" with the old files(old hemispheres.cc and corresponding .hh) it compiled at once.
>
> Regards,
> Debarati
>
>
> ________________________________________
> From: Andy Buckley
> Sent: 26 August 2014 20:32
> To: Debarati Roy; rivet at projects.hepforge.org
> Cc: Albert Knutsson; Maxime Gouzevitch; Lars Sonnenschein
> Subject: Re: [Rivet] Rivet Query
>
> 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
>


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