GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » EMC » Access to clusters properties
Access to clusters properties [message #11757] Tue, 10 May 2011 08:56 Go to previous message
Ronald Kunne is currently offline  Ronald Kunne
Messages: 32
Registered: October 2009
continuous participant
From: *in2p3.fr
Good morning!

I have two questions on access to clusters properties.
The PandaRoot version I use currently is 11787, dated 5/5/11.

Thanks in advance,
Ronald Kunne

1) What is the simplest way to access the contents of the individual crystals, that make up a cluster?

My goal is to distinguish between clusters consisting of one or two gammas, i.e. pizeros.

I would like to do something like this:
  TString inPidTPCFile  = "pid_tpccombi.root";

  TFile *inFile = TFile::Open(inPidTPCFile,"READ");
  TTree *lhe=(TTree *) inFile->Get("cbmsim") ;

...

  TClonesArray* cluster_array=new TClonesArray("PndEmcCluster");
  lhe->SetBranchAddress("EmcCluster",&cluster_array);

...

// loop over  events
    NTevents=lhe->GetEntriesFast();
    for (Int_t j=0; j< NTevents ; j++){
       lhe->GetEntry(j);   // kinematics

// loop over clusters
       Int_t nclusters= cluster_array->GetEntriesFast();
       for (Int_t nc = 0; nc < nclusters; nc++) {
           PndEmcCluster *cl = (PndEmcCluster *)cluster_array->At(nc);
// access to energy and moments of cluster
           Double_t Energy = cl->energy(); 
           Double_t Z20 = cl->Z20(); 
           Double_t Z53 = cl->Z53(); 
           Double_t LatMom = cl->LatMom(); 

// (hypothetical) loop over crystals in this cluster
           Int_t ncrystals= crystal_array->GetEntriesFast();
           for (Int_t nx = 0; nx < ncrystals; nx++) {

              ...etc...
           }
        }
     }


2) The above method gives me access to the moments Z20, Z53 and LatMom.
How do I get access to *all* the moments available in

/emc/EmcData/PndEmcXClMoments.h ?

 
Read Message
Read Message
Read Message
Previous Topic: Uniform particle generation in solid angle
Next Topic: New FwEndCap geometry implemented
Goto Forum:
  


Current Time: Thu Oct 03 14:55:49 CEST 2024

Total time taken to generate the page: 0.00758 seconds