MaterialBudgetHcalHistos.h

Go to the documentation of this file.

#ifndef Validation_Geometry_MaterialBudgetHcalHistos_h
#define Validation_Geometry_MaterialBudgetHcalHistos_h 1

#include "DetectorDescription/Core/interface/DDFilteredView.h"
#include "DetectorDescription/Core/interface/DDsvalues.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "G4Step.hh"
#include "G4Track.hh"

#include <TH1F.h>
#include <TH2F.h>
#include <TProfile.h>
#include <TProfile2D.h>

#include <string>
#include <vector>

class MaterialBudgetHcalHistos {

public:
  
  MaterialBudgetHcalHistos(const edm::ParameterSet &p);
  virtual ~MaterialBudgetHcalHistos() { hend(); }
  
  void fillBeginJob(const DDCompactView &);
  void fillStartTrack(const G4Track*);
  void fillPerStep(const G4Step *);
  void fillEndTrack();
  
private:
  
  void book(); 
  void fillHisto(int ii);
  void fillLayer();
  void hend();
  std::vector<std::string> getNames(DDFilteredView& fv);
  std::vector<double>      getDDDArray(const std::string & str,
                       const DDsvalues_type & sv);
  bool isSensitive(std::string);
  bool isItHF(const G4VTouchable*);
  bool isItEC(std::string);
  
private:

  static const int         maxSet = 25, maxSet2 = 9;
  std::vector<std::string> sensitives, hfNames, sensitiveEC;
  std::vector<int>         hfLevels;
  bool                     fillHistos, printSum;
  int                      binEta, binPhi;
  double                   maxEta, etaLow, etaHigh;
  std::vector<std::string> matList;
  std::vector<double>      stepLength, radLength, intLength;
  TH1F                     *me400[maxSet], *me800[maxSet], *me1300[maxSet2];
  TH2F                     *me1200[maxSet],*me1400[maxSet2];
  TProfile                 *me100[maxSet], *me200[maxSet], *me300[maxSet];
  TProfile                 *me500[maxSet], *me600[maxSet], *me700[maxSet];
  TProfile                 *me1500[maxSet2];
  TProfile2D               *me900[maxSet], *me1000[maxSet],*me1100[maxSet];
  int                      id, layer, steps;
  double                   radLen, intLen, stepLen;
  double                   eta, phi;
  int                      nlayHB, nlayHE, nlayHO, nlayHF;
};


#endif