Skip to Content.
Sympa Menu

sphenix-software-l - Re: [Sphenix-software-l] [Sphenix-jet-structure-l] fast jet -> pseudo jet area?

sphenix-software-l AT lists.bnl.gov

Subject: sPHENIX discussion of software

List archive

Chronological Thread  
  • From: Chris Pinkenburg <pinkenburg AT bnl.gov>
  • To: sphenix-software-l AT lists.bnl.gov
  • Subject: Re: [Sphenix-software-l] [Sphenix-jet-structure-l] fast jet -> pseudo jet area?
  • Date: Wed, 17 May 2017 17:22:22 -0400

if it helps - the fastjet sources for our version is under

/phenix/u/phnxbld/workarea/sPHENIX/fastjet-3.2.1

the message comes from:
/phenix/u/phnxbld/workarea/sPHENIX/fastjet-3.2.1/src/PseudoJet.cc

//----------------------------------------------------------------------
// the following ones require a computation of the area in the
// associated ClusterSequence (See ClusterSequenceAreaBase for details)
//----------------------------------------------------------------------

#ifndef __FJCORE__

//----------------------------------------------------------------------
// if possible, return a valid ClusterSequenceAreaBase pointer; otherwise
// throw an error
const ClusterSequenceAreaBase * PseudoJet::validated_csab() const {
const ClusterSequenceAreaBase *csab = dynamic_cast<const ClusterSequenceAreaBa
se*>(validated_cs());
if (csab == NULL) throw Error("you requested jet-area related information, but the PseudoJet does not have associated area information.");
return csab;
}

If I interpret this correctly we have to tell fastjet at some point explicitly to calculate the area (or do it yourself somewhere)

Chris

On 05/17/2017 05:13 PM, Perepelitsa, Dennis wrote:
Hi Rosi,

To understand what you’re doing, is there a git commit (or a text file with
diffs, or a pointer to a working version) we could look at? I think that
would help me, at least, better understand the problem...

Dennis

On May 17, 2017, at 2:30 PM, Rosi Reed <rjr215 AT lehigh.edu> wrote:

Hello All,

I am implementing the STAR/ALICE background subtraction algorithm in
sPHENIX and ran into an issue that I am sure is rather simple, but I
am not sure what to do about it.

For this algorithm, I need to get the area of the jet. Normally I
access this from a vector of pseudo jets via:
jets[i].area()

I have added all the necessary machinery (and it compiles, thanks to
Chris!) but when I run the fun4all macro, I get the error:
fastjet::Error: you requested jet-area related information, but the
PseudoJet does not have associated area information.

Any thoughts? (Yes, I changed FastJetAlgo.C)

Cheers,
Rosi



--
Rosi Reed
Assistant Professor, Physics Department
Lehigh University
(610)758-3907
16 Memorial Drive East Office 406
Bethlehem, PA 18015
_______________________________________________
sPHENIX-jet-structure-l mailing list
sPHENIX-jet-structure-l AT lists.bnl.gov
https://lists.bnl.gov/mailman/listinfo/sphenix-jet-structure-l
Dennis V. Perepelitsa
Assistant Professor, Physics Department
University of Colorado Boulder



_______________________________________________
sPHENIX-software-l mailing list
sPHENIX-software-l AT lists.bnl.gov
https://lists.bnl.gov/mailman/listinfo/sphenix-software-l


--

*********************************************************

Christopher H. Pinkenburg ; pinkenburg AT bnl.gov

Brookhaven National Laboratory ; phone: (631) 344-5692
Physics Department Bldg 510 C ; fax: (631) 344-3253
Upton, NY 11973-5000

http://www.phenix.bnl.gov/~pinkenbu/

*********************************************************





Archive powered by MHonArc 2.6.24.

Top of Page