CastorPedestalAnalysis.h

Go to the documentation of this file.

#ifndef CastorPedestalAnalysis_H
#define CastorPedestalAnalysis_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/HcalDetId/interface/HcalCastorDetId.h"
#include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"


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

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

// User switches for CastorPedestalAnalysis 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 CastorPedestals;
class CastorPedestalWidths;
class CastorDbService;
class CastorQIEShape;
class CastorQIECoder;
class TFile;
class CastorPedestalAnalysis{
  
public:
  
  CastorPedestalAnalysis(const edm::ParameterSet& ps);  
  ~CastorPedestalAnalysis();

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

  void SampleAnalysis();

  int done(const CastorPedestals* fInputPedestals, 
        const CastorPedestalWidths* fInputWidths,
        CastorPedestals* fOutputPedestals, 
        CastorPedestalWidths* fOutputWidths);

  void processEvent(const CastorDigiCollection& castor,
            const CastorDbService& cond);

// pedestal validation: CastorPedVal=-1 means not validated,
//                                  0 everything OK,
//                                  N>0 : mod(N,100000) drifts + width changes
//                                        int(N/100000) missing channels
  static int CastorPedVal(int nstat[4], const CastorPedestals* fRefPedestals,
            const CastorPedestalWidths* fRefPedestalWidths,
            CastorPedestals* fRawPedestals,
            CastorPedestalWidths* fRawPedestalWidths,
            CastorPedestals* fValPedestals,
            CastorPedestalWidths* 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 CastorDbService& 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 CastorQIEShape* m_shape;
  const CastorQIECoder* 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;
  } castorHists;
  std::map<HcalDetId,std::map<int, PEDBUNCH > >::iterator _meot;
  const CastorPedestals* fRefPedestals;
  const CastorPedestalWidths* fRefPedestalWidths;
  CastorPedestals* fRawPedestals;
  CastorPedestalWidths* fRawPedestalWidths;
  CastorPedestals* fValPedestals;
  CastorPedestalWidths* fValPedestalWidths;
  int evt;
  int sample;
  int evt_curr;
  float m_stat[4];
  std::vector<bool> state;

// flag to make gaussian fits to charge dists
  static const int fitflag=0;
};

#endif