lfhcal_studiesProcessor.h

Go to the documentation of this file.

// Copyright 2023, Friederike Bock
// Subject to the terms in the LICENSE file found in the top-level directory.
//
//  Sections Copyright (C) 2023 Friederike Bock
//  under SPDX-License-Identifier: LGPL-3.0-or-later
 
#pragma once
 
#include <DDSegmentation/BitFieldCoder.h>
#include <JANA/JEventProcessor.h>
#include <JANA/Utils/JTypeInfo.h>
#include <TDirectory.h>
#include <TH2.h>
#include <TH3.h>
#include <TTree.h>
#include <spdlog/logger.h>
#include <memory>
#include <string>
 
class lfhcal_studiesProcessor : public JEventProcessor {
public:
  lfhcal_studiesProcessor() { SetTypeName(NAME_OF_THIS); }
 
  void Init() override;
  //     void InitWithGlobalRootLock() override;
  //     void ProcessSequential(const std::shared_ptr<const JEvent>& event) override;
  void Process(const std::shared_ptr<const JEvent>& event) override;
  void Finish() override;
  //     void FinishWithGlobalRootLock() override;
  TDirectory* m_dir_main;
  // Declare histogram and tree pointers here. e.g.
  TH2D* hMCEnergyVsEta;
  TH3D* hClusterEcalib_E_eta;
  TH3D* hClusterNCells_E_eta;
  TH3D* hClusterEcalib_E_phi;
  TH2D* hPosCaloHitsXY;
  TH2D* hPosCaloHitsZX;
  TH2D* hPosCaloHitsZY;
  TH3D* hClusterESimcalib_E_eta;
  TH3D* hClusterSimNCells_E_eta;
  TH3D* hClusterESimcalib_E_phi;
  TH2D* hCellESim_layerZ;
  TH2D* hCellESim_layerX;
  TH2D* hCellESim_layerY;
  TH2D* hCellTSim_layerZ;
  TH2D* hPosCaloSimHitsXY;
  TH2D* hPosCaloSimHitsZX;
  TH2D* hPosCaloSimHitsZY;
  TH3D* hRecClusterEcalib_E_eta;
  TH3D* hRecNClusters_E_eta;
  TH3D* hRecClusterEcalib_Ehigh_eta;
  TH3D* hRecClusterNCells_Ehigh_eta;
  TH3D* hRecFClusterEcalib_E_eta;
  TH3D* hRecFNClusters_E_eta;
  TH3D* hRecFClusterEcalib_Ehigh_eta;
  TH3D* hRecFClusterNCells_Ehigh_eta;
  TH3D* hRecFEmClusterEcalib_E_eta;
  TH3D* hRecFEmNClusters_E_eta;
  TH3D* hRecFEmClusterEcalib_Ehigh_eta;
  TH2D* hSamplingFractionEta;
 
  bool enableTree = true;
  TTree* event_tree;
  const int maxNTowers = 65000;
  int t_lFHCal_towers_N;
  short* t_lFHCal_towers_cellIDx;
  short* t_lFHCal_towers_cellIDy;
  short* t_lFHCal_towers_cellIDz;
  short* t_lFHCal_towers_clusterIDA;
  short* t_lFHCal_towers_clusterIDB;
  float* t_lFHCal_towers_cellE;
  float* t_lFHCal_towers_cellT;
  int* t_lFHCal_towers_cellTrueID;
 
  bool enableTreeCluster = true;
  bool enableECalCluster = true;
  TTree* cluster_tree;
  const int maxNCluster = 50;
  const int maxNMC      = 50;
  int t_mc_N;
  float* t_mc_E;
  float* t_mc_Phi;
  float* t_mc_Eta;
  int t_lFHCal_clusters_N;
  float* t_lFHCal_cluster_E;
  int* t_lFHCal_cluster_NCells;
  float* t_lFHCal_cluster_Phi;
  float* t_lFHCal_cluster_Eta;
  int t_fEMC_clusters_N;
  float* t_fEMC_cluster_E;
  int* t_fEMC_cluster_NCells;
  float* t_fEMC_cluster_Phi;
  float* t_fEMC_cluster_Eta;
 
  int nEventsWithCaloHits = 0;
  bool isLFHCal           = true;
  std::shared_ptr<spdlog::logger> m_log;
  dd4hep::DDSegmentation::BitFieldCoder* m_decoder;
  std::string nameSimHits       = "LFHCALHits";
  std::string nameRecHits       = "LFHCALRecHits";
  std::string nameClusters      = "LFHCALClusters";
  std::string nameProtoClusters = "LFHCALIslandProtoClusters";
  short iPassive;
  short iLx;
  short iLy;
  short iLz;
};