HcalPedestalAnalysis.h

Go to the documentation of this file.

#ifndef HcalPedestalAnalysis_H
#define HcalPedestalAnalysis_H


#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"

#include "CondFormats/HcalObjects/interface/AllObjects.h"

#include "TH1F.h"
#include "TF1.h"

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

// User switches for HcalPedestalAnalysis to set in cfg are:
//   nevtsample - number of events per sample in which data will be divided
//                for stability checks (default: 0 = do not use),
//   hiSaveflag - flag to save histos of charge per cap-id (default: 0),
//   pedValflag - pedestal validation flag:
//                0 - write out current raw constants (default)
//                1 - write out validated constants

class HcalDbService;
class TFile;

class HcalPedestalAnalysis{
  
public:
  
  HcalPedestalAnalysis(const edm::ParameterSet& ps);  
  ~HcalPedestalAnalysis();

  void setup(const std::string& m_outputFileROOT);

  void SampleAnalysis();

  void setTopology(const HcalTopology* htopo) {fTopology = htopo;}

  int done(const HcalPedestals* fInputPedestals, 
        const HcalPedestalWidths* fInputWidths,
        HcalPedestals* fOutputPedestals, 
        HcalPedestalWidths* fOutputWidths);

  void processEvent(const HBHEDigiCollection& hbhe,
            const HODigiCollection& ho,
            const HFDigiCollection& hf,
            const HcalDbService& cond);

// pedestal validation: HcalPedVal=-1 means not validated,
//                                  0 everything OK,
//                                  N>0 : mod(N,100000) drifts + width changes
//                                        int(N/100000) missing channels
  static int HcalPedVal(int nstat[4], const HcalPedestals* fRefPedestals,
            const HcalPedestalWidths* fRefPedestalWidths,
            HcalPedestals* fRawPedestals,
            HcalPedestalWidths* fRawPedestalWidths,
            HcalPedestals* fValPedestals,
            HcalPedestalWidths* fValPedestalWidths);

protected:
  
  
private:
  //###
  //#  PEDBUNCH is used in map<HcalDetId,map<int, PEDBUNCH > > PEDTRENDS;
  //#  For each HcalDetId (channel) a map<int, PEDBUNCH> is associated;
  //#  int is cap-id (1-4);
  //#  PEDBUNCH is a pair - first element is the main 
  //#  histo with the pedestal distribution and second one is another pair;
  //#  this pair contains map<int, std::vector<double> > as a first element;
  //#  int is cap-id, and vector contains some useful variables;
  //#  the second element is a vector of histos (pointers);
  //#  for the "trend" analysis the main histo (with pedestals) is reset every 
  //#  nevt_ped events and info is put in the other part of the PEDBUNCH;
  //#  so at the end we have the trends for the variables in concern
  //#  which are written in THE vector<TH1F*>; 
  //###  
  typedef std::pair<TH1F*,std::pair<std::map<int, std::vector<double> >,std::vector<TH1F*> > > PEDBUNCH;

  void per2CapsHists(int flag, int id, const HcalDetId detid, const HcalQIESample& qie1, const HcalQIESample& qie2, std::map<HcalDetId, std::map<int,PEDBUNCH> > &toolT,const HcalDbService& cond);

  void GetPedConst(std::map<HcalDetId,std::map<int, PEDBUNCH > > &toolT, TH1F* PedMeans, TH1F* PedWidths);

  void Trendings(std::map<HcalDetId,std::map<int, PEDBUNCH > > &toolT, TH1F* Chi2, TH1F* CapidAverage, TH1F* CapidChi2);

  void AllChanHists(const HcalDetId detid, const HcalQIESample& qie0, const HcalQIESample& qie1, const HcalQIESample& qie2, const HcalQIESample& qie3, const HcalQIESample& qie4, const HcalQIESample& qie5, std::map<HcalDetId, std::map<int,PEDBUNCH> > &toolT);

  TFile* m_file;

  std::string m_outputFileROOT;
  std::string m_outputFileMean;
  std::string m_outputFileWidth;
  std::ofstream m_logFile;
  int m_startTS;
  int m_endTS;
  int m_nevtsample;
  int m_pedsinADC;
  int m_hiSaveflag;
  int m_pedValflag;
  int m_AllPedsOK;
  
  const HcalQIEShape* m_shape;
  const HcalQIECoder* m_coder;
  struct{
    std::map<HcalDetId,std::map<int, PEDBUNCH > > PEDTRENDS;
    TH1F* ALLPEDS;
    TH1F* PEDRMS;
    TH1F* PEDMEAN;
    TH1F* CHI2;
    TH1F* CAPID_AVERAGE;
    TH1F* CAPID_CHI2;
  } hbHists, hfHists, hoHists;
  std::map<HcalDetId,std::map<int, PEDBUNCH > >::iterator _meot;
  const HcalPedestals* fRefPedestals;
  const HcalPedestalWidths* fRefPedestalWidths;
  HcalPedestals* fRawPedestals;
  HcalPedestalWidths* fRawPedestalWidths;
  HcalPedestals* fValPedestals;
  HcalPedestalWidths* fValPedestalWidths;
  int evt;
  int sample;
  int evt_curr;
  float m_stat[4];
  std::vector<bool> state;
  const HcalTopology* fTopology;
// flag to make gaussian fits to charge dists
  static const int fitflag=0;
};

#endif