GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » EMC » emc cluster summary information (rev.2147)
emc cluster summary information (rev.2147) [message #5799] Mon, 04 February 2008 11:53 Go to next message
Dima Melnychuk is currently offline  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
Re: emc cluster summary information (rev.2147) [message #5910 is a reply to message #5799] Tue, 19 February 2008 14:22 Go to previous messageGo to next message
M.Babai is currently offline  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
Re: emc cluster summary information (rev.2147) [message #5915 is a reply to message #5910] Tue, 19 February 2008 16:43 Go to previous messageGo to next message
StefanoSpataro is currently offline  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.
Re: emc cluster summary information (rev.2147) [message #5916 is a reply to message #5915] Tue, 19 February 2008 17:00 Go to previous messageGo to next message
M.Babai is currently offline  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 .
Re: emc cluster summary information (rev.2147) [message #5917 is a reply to message #5916] Tue, 19 February 2008 17:18 Go to previous messageGo to next message
StefanoSpataro is currently offline  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.

Re: emc cluster summary information (rev.2147) [message #5918 is a reply to message #5917] Tue, 19 February 2008 17:27 Go to previous messageGo to next message
Dima Melnychuk is currently offline  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
Re: emc cluster summary information (rev.2147) [message #5931 is a reply to message #5918] Wed, 20 February 2008 11:37 Go to previous message
M.Babai is currently offline  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.
Previous Topic: EMC small updates
Next Topic: Bug correction in EmcMapper (rev. 2199)
Goto Forum:
  


Current Time: Mon May 20 21:32:42 CEST 2024

Total time taken to generate the page: 0.01068 seconds