EcalPerEvtLaserAnalyzer.h

Go to the documentation of this file.

// $Id: EcalPerEvtLaserAnalyzer.h

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

class TTree;
class TFile;

// Define geometrical constants
// NOT the same for "EB" and "EE"
//
//     "EB"       "EE"
//
//      0          0
//   1     2    1     2
//   3     4
//   5     6
//   7     8
//
//
#define NSIDES 2  // Number of sides (0, 1)

// "EB" geometry
#define NCRYSEB 1700  // Number of crystals per EB supermodule
#define NTTEB 68      // Number of EB Trigger Towers
#define NPNEB 10      // Number of PN per EB supermodule

// "EE" geometry
#define NCRYSEE 825  // Number of crystals per EE supermodule
#define NTTEE 33     // Number of EE Trigger Towers
#define NPNEE 4      // Number of PN per EE supermodule

class EcalPerEvtLaserAnalyzer : public edm::EDAnalyzer {
public:
  explicit EcalPerEvtLaserAnalyzer(const edm::ParameterSet &iConfig);
  ~EcalPerEvtLaserAnalyzer() 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 _samplemin;
  unsigned int _samplemax;
  unsigned int _nsamplesPN;
  unsigned int _presamplePN;
  unsigned int _firstsamplePN;
  unsigned int _lastsamplePN;
  unsigned int _timingcutlow;
  unsigned int _timingcuthigh;
  unsigned int _niter;
  int _fedid;
  unsigned int _tower;
  unsigned int _channel;
  std::string _ecalPart;

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

  // Output file names

  std::string alphafile;
  std::string ADCfile;
  std::string APDfile;
  std::string resfile;

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

  unsigned int nCrys;
  unsigned int nTT;
  unsigned int nSides;

  int IsFileCreated;

  // Define number of channels (68 or 1700) for alpha and beta calculation

  unsigned int nCh;

  // Identify run type

  int runType;
  int runNum;
  int dccID;
  int fedID;
  int lightside;

  int channelNumber;

  std::vector<int> colors;

  // Temporary root files and trees
  TFile *matacqFile;
  TTree *matacqTree;

  TFile *ADCFile;
  TTree *ADCtrees;

  TFile *APDFile;
  TTree *header[2];
  TTree *APDtrees;

  int doesRefFileExist;
  int IsThereDataADC[nColor];
  int IsHeaderFilled[nColor];

  double ttMat, peakMat, peak;
  int evtMat, colMat;

  // Declaration of leaves types for temporary trees

  int phi, eta;
  int event;
  int color;
  double adc[10];
  int adcG[10];
  double tt;
  double ttrig;
  double pn0, pn1;
  double pn[50];
  double apdAmpl;
  double apdTime;
  double pnAmpl;
};