GSI Forum: EMC » emc cluster summary information (rev.2147)

Home » PANDA » PandaRoot » EMC » emc cluster summary information (rev.2147)

Show: Today's Messages :: Polls :: Message Navigator

emc cluster summary information (rev.2147) [message #5799]

Mon, 04 February 2008 11:53

Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...

first-grade participant

From: *fuw.edu.pl

Hi all,

I have ported from Babar framework several classes for emc cluster summary information.

PndEmcClusterEnergySums - gives information on energy of central crystal, central 9 and 25 crystals and ratios between them.

PndEmcClusterDistances - radial and angular distances from specific digi to cluster centroid.

PndEmcClusterMoments - information on different cluster moments.

PndEmcXClMoments - cluster Zernike moments.

Macro reco_analys3.C demonstrates some capabilities of PndEmcClusterEnergySums, i.e. plots energy of cluster, energy of central crystal, 9 and 25 central crystals.

Best regards,
Dima

Report message to a moderator

Re: emc cluster summary information (rev.2147) [message #5910 is a reply to message #5799]

Tue, 19 February 2008 14:22

M.Babai
Messages: 46
Registered: January 2008
Location: Netherlands

continuous participant

From: *KVI.nl

Hi all,

I'm trying to use these cluster properties to perform particle identification.

If I want to initialize a PndEmcXClMoments object, root starts to complain and produces's an error with the following message:

Wrong call for PndEmcStructure. It should be instantiated first with TGeoManager

And this is how I'm using this function:

//Select cluster 
PndEmcCluster* cl = (PndEmcCluster*) cluster_array->At(idWithHighestEnergy);

 PndEmcXClMoments zmom = cl->Xmoments();
 ..........

Any body, any idea, how to solve this?

Kind Wishes,
MB

Report message to a moderator

Re: emc cluster summary information (rev.2147) [message #5915 is a reply to message #5910]

Tue, 19 February 2008 16:43

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *physik.uni-giessen.de

I think you should open in your macro even the sim_emc.root file, the one with the geometry.

Report message to a moderator

Re: emc cluster summary information (rev.2147) [message #5916 is a reply to message #5915]

Tue, 19 February 2008 17:00

M.Babai
Messages: 46
Registered: January 2008
Location: Netherlands

continuous participant

From: *KVI.nl

I open the simulation file too, but the the error remains the same. Sad

Report message to a moderator

Re: emc cluster summary information (rev.2147) [message #5917 is a reply to message #5916]

Tue, 19 February 2008 17:18

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *physik.uni-giessen.de

Maybe the lines:

TFile* fsim = new TFile("sim_emc.root");
TTree *tsim=(TTree *) fsim->Get("cbmsim");
PndEmcMapper *emcMap=PndEmcMapper::Instance(1);

could help. I had a similar prolem, and I solved it in this way.
Try and let me know.

Report message to a moderator

Re: emc cluster summary information (rev.2147) [message #5918 is a reply to message #5917]

Tue, 19 February 2008 17:27

Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...

first-grade participant

From: *fuw.edu.pl

Hi Mohammad,

2 remarks:

1) update the code to the rev. 2198, I made the correction in PndEmcXClMoments.cxx

2) You should instantiate PndEmcStructure before accessing PndEmcXClMoments.

TFile* fsim = new TFile("sim_emc.root");
TTree *tsim=(TTree *) fsim->Get("cbmsim");
PndEmcMapper *emcMap=PndEmcMapper::Instance(1);

TGeoManager *geoMan = (TGeoManager*) gROOT->FindObject("CBMGeom");
PndEmcStructure::Instance(geoMan);

Dima

Report message to a moderator

Re: emc cluster summary information (rev.2147) [message #5931 is a reply to message #5918]

Wed, 20 February 2008 11:37

M.Babai
Messages: 46
Registered: January 2008
Location: Netherlands

continuous participant

From: *KVI.nl

Thank you for your replies.
this initialization seems to solve the problem

TGeoManager *geoMan = (TGeoManager*) gROOT->FindObject("CBMGeom");
PndEmcStructure::Instance(geoMan);

Kind Wishes,
MB.

Report message to a moderator

Previous Topic:	EMC small updates
Next Topic:	Bug correction in EmcMapper (rev. 2199)

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sat Feb 21 06:28:59 CET 2026

Total time taken to generate the page: 0.00589 seconds