Commit 68846b07 authored by Marcin's avatar Marcin
Browse files

CoreE0 Energies finally working: bug in the size of the Corre0 pid gated array...

CoreE0 Energies finally working: bug in the size of the Corre0 pid gated array corrected, cleaning/changing names of CoreE0 histograms
parent e52000a2
......@@ -442,25 +442,25 @@ void PlotResultsSelector::SlaveBegin(TTree * /*tree*/)
#ifdef USE_CORE_ENERGIES
// CoreE0 histograms
AGATA_CoreE0_h = BuildHist1D("AGATA_CoreE0",EGamma2,fAgataHistList);
CoreE0_h = BuildHist1D("CoreE0",EGamma2,fAgataHistList);
l = new TList;l->SetName("CoreE0_ngated");fAgataHistList->Add(l);
l = new TList;l->SetName("ngated");fAgataHistList->Add(l);
for (Int_t i=0; i<MaxNeutrons; i++)
{
TString hist_name = Form("CoreE0_ngated%i",i);
AGATA_CoreE0_ngated_h[i] = BuildHist1D(hist_name,EGamma2,l);
TString hist_name = Form("ngated%i",i);
CoreE0_ngated_h[i] = BuildHist1D(hist_name,EGamma2,l);
}
l = new TList;l->SetName("CoreE0_agated");fAgataHistList->Add(l);
l = new TList;l->SetName("agated");fAgataHistList->Add(l);
for (Int_t i=0; i<MaxAlphas; i++)
{
TString hist_name = Form("CoreE0_agated%i",i);
AGATA_CoreE0_agated_h[i] = BuildHist1D(hist_name,EGamma2,l);
TString hist_name = Form("agated%i",i);
CoreE0_agated_h[i] = BuildHist1D(hist_name,EGamma2,l);
}
l = new TList;l->SetName("CoreE0_pdx");fAgataHistList->Add(l);
l = new TList;l->SetName("pdx");fAgataHistList->Add(l);
for (Int_t i=0; i<MaxParticleidx; i++)
{
TString hist_name = Form("CoreE0_pdx%i",i);
AGATA_CoreE0_pdx_h[i] = BuildHist1D(hist_name,EGamma2,l);
TString hist_name = Form("pdx%03i",i);
CoreE0_pdx_h[i] = BuildHist1D(hist_name,EGamma2,l);
}
#endif
......@@ -949,12 +949,12 @@ void PlotResultsSelector::FillHistos()
}
}
#ifdef USE_CORE_ENERGIES //this piece of code is not finished yet
#ifdef USE_CORE_ENERGIES
for ( int ig1 =0 ; ig1<nbCores; ig1++) {
AGATA_CoreE0_h->Fill (CoreE0[ig1]);
AGATA_CoreE0_ngated_h[in]->Fill (CoreE0[ig1]);
AGATA_CoreE0_agated_h[ia]->Fill (CoreE0[ig1]);
AGATA_CoreE0_pdx_h[particle_idx]->Fill (CoreE0[ig1]);
CoreE0_h->Fill (CoreE0[ig1]);
CoreE0_ngated_h[in]->Fill (CoreE0[ig1]);
CoreE0_agated_h[ia]->Fill (CoreE0[ig1]);
CoreE0_pdx_h[particle_idx]->Fill (CoreE0[ig1]);
}
#endif
......
......@@ -33,11 +33,11 @@
#include "TEnv.h"
// #define USE_NEDA_NZCO
// #define USE_CORE_ENERGIES
#define USE_CORE_ENERGIES
#define USE_NEDA_CCTOF_CUTS
#define USE_NEDA_ENERGYCC_CUTS
#define USE_NEDA_TIME_LIMS
//#define REMOVE_EMPTY_HISTOGRAMS
//#define REMOVE_EMPgTY_HISTOGRAMS
using namespace std;
......@@ -291,10 +291,10 @@ public :
TH1F *AGATA_TrackE_particleidx_h[MaxParticleidx];
#ifdef USE_CORE_ENERGIES
TH1F *AGATA_CoreE0_h;
TH1F *AGATA_CoreE0_ngated_h[MaxNeutrons];
TH1F *AGATA_CoreE0_agated_h[MaxAlphas];
TH1F *AGATA_CoreE0_pdx_h[MaxNeutrons];
TH1F *CoreE0_h;
TH1F *CoreE0_ngated_h[MaxNeutrons];
TH1F *CoreE0_agated_h[MaxAlphas];
TH1F *CoreE0_pdx_h[MaxParticleidx];
#endif
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment