EcalABAnalyzer.h

Go to the documentation of this file.

#ifndef EcalABAnalyzer_h_
#define EcalABAnalyzer_h_

// $Id: EcalABAnalyzer.h

#include <memory>
#include <vector>
#include <map>

#include <FWCore/Framework/interface/EDAnalyzer.h>

class TShapeAnalysis;
class TAPDPulse;
class TMom;

// Define geometrical constants
// NOT the same for "EB" and "EE"
//
//     "EB"       "EE"
//
//      0          0
//   1     2    1     2
//   3     4
//   5     6
//   7     8
//
//

// "EB" geometry
#define NCRYSEB 1700  // Number of crystals per EB supermodule

// "EE" geometry
#define NCRYSEE 830  // Number of crystals per EE supermodule

class EcalABAnalyzer : public edm::EDAnalyzer {
public:
  explicit EcalABAnalyzer(const edm::ParameterSet &iConfig);
  ~EcalABAnalyzer() override;

  void analyze(const edm::Event &e, const edm::EventSetup &c) override;
  void beginJob() override;
  void endJob() override;

  enum VarCol { iBlue, iRed, nColor };

private:
  int iEvent;

  // Framework parameters

  unsigned int _nsamples;
  unsigned int _presample;
  unsigned int _firstsample;
  unsigned int _lastsample;
  unsigned int _timingcutlow;
  unsigned int _timingcuthigh;
  unsigned int _timingquallow;
  unsigned int _timingqualhigh;
  double _ratiomincutlow;
  double _ratiomincuthigh;
  double _ratiomaxcutlow;
  double _presamplecut;
  unsigned int _niter;
  double _alpha;
  double _beta;
  unsigned int _nevtmax;
  double _noise;
  double _chi2cut;
  std::string _ecalPart;
  int _fedid;
  double _qualpercent;
  int _debug;

  TAPDPulse *APDPulse;
  TMom *Delta01;
  TMom *Delta12;

  std::string resdir_;
  std::string digiCollection_;
  std::string digiProducer_;
  std::string eventHeaderCollection_;
  std::string eventHeaderProducer_;

  // Output file names

  std::string alphafile;
  std::string alphainitfile;

  TShapeAnalysis *shapana;
  unsigned int nevtAB[NCRYSEB];

  //  Define geometrical constants
  //  Default values correspond to "EB" geometry (1700 crystals)

  unsigned int nCrys;
  bool doesABTreeExist;

  bool _fitab;
  // Identify run type

  int runType;
  int runNum;
  int fedID;
  int dccID;
  int side;
  int lightside;
  int iZ;

  // Temporary root files and trees

  std::vector<int> colors;
  std::map<int, int> channelMapEE;
  std::vector<int> dccMEM;
  std::vector<int> modules;

  // Declaration of leaves types for temporary trees

  int phi, eta;
  int event;
  int color;
  double adc[10];
  int adcG[10];
  int channelIteratorEE;

  int iEta[NCRYSEB], iPhi[NCRYSEB];
  int iTowerID[NCRYSEB], iChannelID[NCRYSEB], idccID[NCRYSEB], iside[NCRYSEB];

  // Quality Checks variables and flags

  int nEvtBadGain[NCRYSEB];
  int nEvtBadTiming[NCRYSEB];
  int nEvtTot[NCRYSEB];

  bool wasGainOK[NCRYSEB];
  bool wasTimingOK[NCRYSEB];

  bool isGainOK;
  bool isTimingOK;
};

#endif