QualityFlagAnalyzer.h

Go to the documentation of this file.

 
#ifndef TRIGSCINT_QUALITYFLAGANALYZER_H
#define TRIGSCINT_QUALITYFLAGANALYZER_H
 
// LDMX Framework
#include "Framework/Configure/Parameters.h"
#include "Framework/EventProcessor.h"  //Needed to declare processor
#include "TH1.h"
#include "TH2.h"
#include "TrigScint/Event/EventReadout.h"
#include "TrigScint/Event/TestBeamHit.h"
 
namespace trigscint {
 
class QualityFlagAnalyzer : public framework::Analyzer {
 public:
  QualityFlagAnalyzer(
      const std::string& name,
      framework::Process& process);  // : framework::Analyzer(name, process) {}
  virtual ~QualityFlagAnalyzer() = default;
  void configure(framework::config::Parameters& parameters) override;
 
  void analyze(const framework::Event& event) override;
 
  void onProcessStart() override;
 
  void onProcessEnd() override;
 
 private:
  std::vector<std::vector<TH1F*> > vChargeVsTime;
 
  // configurable parameters
  std::string inputEventCol_;  // full event stream input
  std::string inputEventPassName_{""};
  std::string inputHitCol_;  // hit collection
  std::string inputHitPassName_{""};
  std::vector<double> peds_;
  std::vector<double> gain_;
  int startSample_{0};
 
  // plotting stuff
  const int nEv{200};
  static constexpr int nChannels{16};
  const int nFlags{6};
  int peFillNb{0};
 
  // make sure to match constants above
  const int flags[6] = {16, 8, 4,
                        2,  1, 0};  // this order just makes looping easier
  int nEvDrawn[6] = {0};  // keep a counter for each flag type to get good stats
  TH1F* hOut[200][16];
  TH1F* hOutPE[200][16];
  TH1F* hOutFlag[6][200][16];  // for 4 quality flags and 0 (no flag)
  TH1F* hPE[16];
 
  TH2F* hTDCfireChanvsEvent;
};
}  // namespace trigscint
 
#endif /* TRIGSCINT_QUALITYFLAGANALYZER_H */