HcalTestAnalysis.h

Go to the documentation of this file.

#ifndef SimG4CMS_HcalTestAnalysis_H
#define SimG4CMS_HcalTestAnalysis_H
// File: HcalTestAnalysis.h
// Analysis of simhits inside the OSCAR framework

#include "SimG4Core/Notification/interface/Observer.h"
#include "SimG4Core/Watcher/interface/SimWatcher.h"

#include "SimDataFormats/CaloHit/interface/CaloHit.h"
#include "SimDataFormats/CaloTest/interface/HcalTestHistoClass.h"
#include "SimG4CMS/Calo/interface/HcalQie.h"
#include "SimG4CMS/Calo/interface/HcalTestHistoManager.h"
#include "SimG4CMS/Calo/interface/HcalTestNumberingScheme.h"
#include "Geometry/HcalCommonData/interface/HcalNumberingFromDDD.h"
#include "Geometry/HcalCommonData/interface/HcalDDDSimConstants.h"

#include <iostream>
#include <memory>
#include <vector>
#include <string>

class G4Step;
class BeginOfJob;
class BeginOfRun;
class BeginOfEvent;
class EndOfEvent;

namespace CLHEP {
  class HepRandomEngine;
}

class HcalTestAnalysis : public SimWatcher,
                         public Observer<const BeginOfJob *>,
                         public Observer<const BeginOfRun *>,
                         public Observer<const BeginOfEvent *>,
                         public Observer<const EndOfEvent *>,
                         public Observer<const G4Step *> {
public:
  HcalTestAnalysis(const edm::ParameterSet &p);
  ~HcalTestAnalysis() override;

private:
  // observer classes
  void update(const BeginOfJob *run) override;
  void update(const BeginOfRun *run) override;
  void update(const BeginOfEvent *evt) override;
  void update(const EndOfEvent *evt) override;
  void update(const G4Step *step) override;

  // analysis-related stuff
  std::vector<int> layerGrouping(int);
  std::vector<int> towersToAdd(int centre, int nadd);
  void fill(const EndOfEvent *ev);
  void qieAnalysis(CLHEP::HepRandomEngine *);
  void layerAnalysis();
  double timeOfFlight(int det, int layer, double eta);

private:
  //Keep parameters to instantiate HcalTestHistoClass later
  std::string fileName_;

  // Qie Analysis
  std::unique_ptr<HcalQie> myqie_;
  int addTower_;

  // Private Tuples
  std::unique_ptr<HcalTestHistoManager> tuplesManager_;
  HcalTestHistoClass *tuples_;

  // Numbering scheme
  std::unique_ptr<HcalNumberingFromDDD> numberingFromDDD_;
  const HcalDDDSimConstants *hcons_;
  HcalTestNumberingScheme *org_;

  // Hits for qie analysis
  std::vector<CaloHit> caloHitCache_;
  std::vector<int> group_, tower_;
  int nGroup_, nTower_;

  // to read from ParameterSet
  std::vector<std::string> names_;
  double eta0_, phi0_;
  int centralTower_;

  // some private members for ananlysis
  unsigned int count_;
  double edepEB_, edepEE_, edepHB_, edepHE_;
  double edepHO_, edepl_[20];
  double mudist_[20];  // Distance of muon from central part
};

#endif