EcalMatacqAnalyzer.h

Go to the documentation of this file.

// $Id: EcalMatacqAnalyzer.h

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

class TFile;
class TTree;
class TMTQ;

#define N_samples 2560
#define N_channels 1
#define NSIDES     2     // Number of sides

class EcalMatacqAnalyzer: public edm::EDAnalyzer{  

 public:
  
  explicit EcalMatacqAnalyzer(const edm::ParameterSet& iConfig);  
  ~EcalMatacqAnalyzer() override;
  
  
  void analyze( const edm::Event & e, const  edm::EventSetup& c) override;
  void beginJob() override;
  void endJob() override;
  
  enum VarCol  { iBlue, iRed, nColor }; 
  enum VarSide { iSide0, iSide1, nSide }; 
  
 private:
    
  int iEvent;

//
// Framework parameters
//

  double       _presample;
  unsigned int _nsamplesaftmax;
  unsigned int _noiseCut;
  unsigned int _parabnbefmax;
  unsigned int _parabnaftmax;
  unsigned int _thres;
  unsigned int _lowlev;
  unsigned int _highlev;
  unsigned int _nevlasers;
  unsigned int _timebefmax;
  unsigned int _timeaftmax;
  double       _cutwindow;
  unsigned int _nsamplesshape;
  unsigned int _presampleshape;
  unsigned int _slide;
  int          _fedid;
  int          _debug;

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

  std::string outfile;
  std::string sampfile;


  // Identify run type

  unsigned int nSides;
  int lightside;
  int runType;
  int runNum;
  int  dccID;
  int  fedID;

  // Count Laser Events
  int laserEvents;
  bool isThereMatacq;

  //Declaration of leaves types
  
  int   event ;
  int   color ;
  double  matacq[N_samples]; 
  int  maxsamp; 
  int  nsamples; 
  double tt;

  TFile *sampFile;
  TTree *tree;

  TMTQ *MTQ[nColor][nSide];
  TTree *meanTree[nColor];

  std::vector<int> colors;

  TFile *outFile;
  int status;
  double peak, sigma, fit, ampl, trise, fwhm, fw20, fw80, ped, pedsig, ttrig, sliding;
  TTree* mtqShape;


};