EcalPerEvtLaserAnalyzer Class Reference

#include <EcalPerEvtLaserAnalyzer.h>

Inheritance diagram for EcalPerEvtLaserAnalyzer:

Public Types
enum	VarCol { iBlue, iRed, nColor }
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	beginJob () override
	EcalPerEvtLaserAnalyzer (const edm::ParameterSet &iConfig)
void	endJob () override
	~EcalPerEvtLaserAnalyzer () override
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const noexcept
bool	wantsStreamRuns () const noexcept
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
const unsigned int	_channel
const std::string	_ecalPart
const int	_fedid
const unsigned int	_firstsample
const unsigned int	_firstsamplePN
const unsigned int	_lastsample
const unsigned int	_lastsamplePN
const unsigned int	_niter
const unsigned int	_nsamples
const unsigned int	_nsamplesPN
const unsigned int	_presample
const unsigned int	_presamplePN
const unsigned int	_timingcuthigh
const unsigned int	_timingcutlow
const unsigned int	_tower
double	adc [10]
std::string	ADCfile
TFile *	ADCFile
int	adcG [10]
TTree *	ADCtrees
double	apdAmpl
TFile *	APDFile
double	apdTime
TTree *	APDtrees
int	channelNumber
int	colMat
int	color
std::vector< int >	colors
int	dccID
const std::string	digiCollection_
const std::string	digiPNCollection_
const std::string	digiProducer_
int	doesRefFileExist
edm::EDGetTokenT< EBDigiCollection >	ebDigiToken_
edm::EDGetTokenT< EEDigiCollection >	eeDigiToken_
int	eta
int	event
const std::string	eventHeaderCollection_
const std::string	eventHeaderProducer_
int	evtMat
int	fedID
TTree *	header [2]
int	iEvent
int	IsFileCreated
int	IsHeaderFilled [nColor]
int	IsThereDataADC [nColor]
int	lightside
const edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	mappingToken_
TFile *	matacqFile
TTree *	matacqTree
unsigned int	nCrys
double	peak
double	peakMat
int	phi
double	pn [50]
double	pn0
double	pn1
double	pnAmpl
const edm::EDGetTokenT< EcalPnDiodeDigiCollection >	pnDiodeDigiToken_
const edm::EDGetTokenT< EcalRawDataCollection >	rawDataToken_
const std::string	refalphabeta_
const std::string	resdir_
std::string	resfile
int	runNum
int	runType
double	tt
double	ttMat
double	ttrig

Additional Inherited Members
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 43 of file EcalPerEvtLaserAnalyzer.h.

Member Enumeration Documentation

VarCol

enum EcalPerEvtLaserAnalyzer::VarCol

Enumerator
iBlue
iRed
nColor

Definition at line 52 of file EcalPerEvtLaserAnalyzer.h.

{ iBlue, iRed, nColor };

Constructor & Destructor Documentation

EcalPerEvtLaserAnalyzer()

EcalPerEvtLaserAnalyzer::EcalPerEvtLaserAnalyzer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 36 of file EcalPerEvtLaserAnalyzer.cc.

References _ecalPart, digiCollection_, digiProducer_, ebDigiToken_, eeDigiToken_, ProducerED_cfi::InputTag, nCrys, NCRYSEB, and NCRYSEE.

    : iEvent(0),
      eventHeaderCollection_(iConfig.getParameter<std::string>("eventHeaderCollection")),
      eventHeaderProducer_(iConfig.getParameter<std::string>("eventHeaderProducer")),
      digiCollection_(iConfig.getParameter<std::string>("digiCollection")),
      digiProducer_(iConfig.getParameter<std::string>("digiProducer")),
      digiPNCollection_(iConfig.getParameter<std::string>("digiPNCollection")),
      rawDataToken_(consumes<EcalRawDataCollection>(edm::InputTag(eventHeaderProducer_, eventHeaderCollection_))),
      pnDiodeDigiToken_(consumes<EcalPnDiodeDigiCollection>(edm::InputTag(digiProducer_, digiPNCollection_))),
      mappingToken_(esConsumes()),
      // framework parameters with default values
      _nsamples(iConfig.getUntrackedParameter<unsigned int>("nSamples", 10)),
      _presample(iConfig.getUntrackedParameter<unsigned int>("nPresamples", 3)),
      _firstsample(iConfig.getUntrackedParameter<unsigned int>("firstSample", 1)),
      _lastsample(iConfig.getUntrackedParameter<unsigned int>("lastSample", 2)),
      _nsamplesPN(iConfig.getUntrackedParameter<unsigned int>("nSamplesPN", 50)),
      _presamplePN(iConfig.getUntrackedParameter<unsigned int>("nPresamplesPN", 6)),
      _firstsamplePN(iConfig.getUntrackedParameter<unsigned int>("firstSamplePN", 7)),
      _lastsamplePN(iConfig.getUntrackedParameter<unsigned int>("lastSamplePN", 8)),
      _timingcutlow(iConfig.getUntrackedParameter<unsigned int>("timingCutLow", 3)),
      _timingcuthigh(iConfig.getUntrackedParameter<unsigned int>("timingCutHigh", 7)),
      _niter(iConfig.getUntrackedParameter<unsigned int>("nIter", 3)),
      _fedid(iConfig.getUntrackedParameter<unsigned int>("fedID", 0)),
      _tower(iConfig.getUntrackedParameter<unsigned int>("tower", 1)),
      _channel(iConfig.getUntrackedParameter<unsigned int>("channel", 1)),
      _ecalPart(iConfig.getUntrackedParameter<std::string>("ecalPart", "EB")),
      resdir_(iConfig.getUntrackedParameter<std::string>("resDir")),
      refalphabeta_(iConfig.getUntrackedParameter<std::string>("refAlphaBeta")),
      nCrys(NCRYSEB),
      IsFileCreated(0),
      runType(-1),
      runNum(0),
      dccID(-1),
      lightside(2),
      doesRefFileExist(0),
      ttMat(-1),
      peakMat(-1),
      peak(-1),
      evtMat(-1),
      colMat(-1)
//========================================================================

{
  if (_ecalPart == "EB") {
    ebDigiToken_ = consumes<EBDigiCollection>(edm::InputTag(digiProducer_, digiCollection_));
  } else if (_ecalPart == "EE") {
    eeDigiToken_ = consumes<EEDigiCollection>(edm::InputTag(digiProducer_, digiCollection_));
  }

  // Define geometrical constants
  //
  if (_ecalPart == "EB") {
    nCrys = NCRYSEB;
  } else {
    nCrys = NCRYSEE;
  }
}

~EcalPerEvtLaserAnalyzer()

EcalPerEvtLaserAnalyzer::~EcalPerEvtLaserAnalyzer ( )

override

Definition at line 96 of file EcalPerEvtLaserAnalyzer.cc.

                                                  {
  //========================================================================

  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

analyze()

void EcalPerEvtLaserAnalyzer::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 149 of file EcalPerEvtLaserAnalyzer.cc.

References _channel, _ecalPart, _fedid, _firstsamplePN, _lastsamplePN, _nsamples, _nsamplesPN, _presample, _presamplePN, _timingcuthigh, _timingcutlow, _tower, adc, adcG, ADCtrees, cms::cuda::assert(), edm::SortedCollection< T, SORT >::begin(), DummyCfis::c, EcalElectronicsId::channelId(), channelNumber, colMat, color, MEEEGeom::crystal(), dccID, das-up-to-nevents::df, digiCollection_, digiPNCollection_, doesRefFileExist, TPNFit::doFit(), MillePedeFileConverter_cfg::e, CollectionTags_cfi::EBDigi, ebDigiToken_, CollectionTags_cfi::EEDigi, eeDigiToken_, MEEBGeom::electronic_channel(), edm::SortedCollection< T, SORT >::end(), eta, eventHeaderCollection_, evtMat, fedID, spr::find(), TPNFit::getAmpl(), mps_fire::i, iEvent, TPNFit::init(), IsFileCreated, edm::HandleBase::isValid(), dqmdumpme::k, EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_STD, lightside, MEEEGeom::lmmod(), MEEBGeom::lmmod(), MEEBGeom::localCoord(), mappingToken_, matacqFile, matacqTree, METSkim_cff::Max, nCrys, peak, peakMat, phi, MEEBGeom::pn(), MEEEGeom::pn(), pn, pn0, pn1, pnAmpl, pnDiodeDigiToken_, edm::Handle< T >::product(), rawDataToken_, resdir_, resfile, runNum, runType, MEEBGeom::side(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), nano_mu_digi_cff::strip, EcalElectronicsId::stripId(), ctpps_dqm_sourceclient-live_cfg::test, ecaldqm::towerID(), EcalElectronicsId::towerId(), ttMat, ttrig, EcalDCCHeaderBlock::EcalDCCEventSettings::wavelength, and EcalElectronicsId::xtalId().

                                                                               {
  //========================================================================

  ++iEvent;

  // retrieving DCC header
  edm::Handle<EcalRawDataCollection> pDCCHeader;
  const EcalRawDataCollection* DCCHeader = nullptr;
  e.getByToken(rawDataToken_, pDCCHeader);
  if (!pDCCHeader.isValid()) {
    edm::LogError("nodata") << "Error! can't get the product  retrieving DCC header" << eventHeaderCollection_.c_str();
  } else {
    DCCHeader = pDCCHeader.product();
  }

  // retrieving crystal data from Event
  edm::Handle<EBDigiCollection> pEBDigi;
  const EBDigiCollection* EBDigi = nullptr;
  edm::Handle<EEDigiCollection> pEEDigi;
  const EEDigiCollection* EEDigi = nullptr;
  if (_ecalPart == "EB") {
    e.getByToken(ebDigiToken_, pEBDigi);
    if (!pEBDigi.isValid()) {
      edm::LogError("nodata") << "Error! can't get the product retrieving EB crystal data " << digiCollection_.c_str();
    } else {
      EBDigi = pEBDigi.product();
    }
  } else {
    e.getByToken(eeDigiToken_, pEEDigi);
    if (!pEEDigi.isValid()) {
      edm::LogError("nodata") << "Error! can't get the product retrieving EE crystal data " << digiCollection_.c_str();
    } else {
      EEDigi = pEEDigi.product();
    }
  }

  // retrieving crystal PN diodes from Event
  edm::Handle<EcalPnDiodeDigiCollection> pPNDigi;
  const EcalPnDiodeDigiCollection* PNDigi = nullptr;
  e.getByToken(pnDiodeDigiToken_, pPNDigi);
  if (!pPNDigi.isValid()) {
    edm::LogError("nodata") << "Error! can't get the product " << digiPNCollection_.c_str();
  } else {
    PNDigi = pPNDigi.product();
  }

  // retrieving electronics mapping
  const auto& TheMapping = c.getData(mappingToken_);

  // ====================================
  // Decode Basic DCCHeader Information
  // ====================================

  for (EcalRawDataCollection::const_iterator headerItr = DCCHeader->begin(); headerItr != DCCHeader->end();
       ++headerItr) {
    // Get run type and run number

    int fed = headerItr->fedId();

    if (fed != _fedid && _fedid != -999)
      continue;

    runType = headerItr->getRunType();
    runNum = headerItr->getRunNumber();
    event = headerItr->getLV1();
    dccID = headerItr->getDccInTCCCommand();
    fedID = headerItr->fedId();

    // take event only if the fed corresponds to the DCC in TCC
    if (600 + dccID != fedID)
      continue;

    // Cut on runType
    if (runType != EcalDCCHeaderBlock::LASER_STD && runType != EcalDCCHeaderBlock::LASER_GAP)
      return;

    // Define output results files' names

    if (IsFileCreated == 0) {
      stringstream namefile2;
      namefile2 << resdir_ << "/APDAmpl_Run" << runNum << "_" << _fedid << "_" << _tower << "_" << _channel << ".root";
      resfile = namefile2.str();

      // Get Matacq ttrig

      stringstream namefile;
      namefile << resdir_ << "/Matacq-Run" << runNum << ".root";

      doesRefFileExist = 0;

      FILE* test;
      test = fopen(namefile.str().c_str(), "r");
      if (test)
        doesRefFileExist = 1;

      if (doesRefFileExist == 1) {
        matacqFile = new TFile((namefile.str().c_str()));
        matacqTree = (TTree*)matacqFile->Get("MatacqShape");

        matacqTree->SetBranchAddress("event", &evtMat);
        matacqTree->SetBranchAddress("color", &colMat);
        matacqTree->SetBranchAddress("peak", &peakMat);
        matacqTree->SetBranchAddress("ttrig", &ttMat);
      }

      IsFileCreated = 1;
    }

    // Retrieve laser color and event number

    EcalDCCHeaderBlock::EcalDCCEventSettings settings = headerItr->getEventSettings();
    int color = settings.wavelength;

    vector<int>::iterator iter = find(colors.begin(), colors.end(), color);
    if (iter == colors.end()) {
      colors.push_back(color);
      edm::LogVerbatim("EcalPerEvtLaserAnalyzer") << " new color found " << color << " " << colors.size();
    }
  }

  // cut on fedID

  if (fedID != _fedid && _fedid != -999)
    return;

  // ======================
  // Decode PN Information
  // ======================

  TPNFit* pnfit = new TPNFit();
  pnfit->init(_nsamplesPN, _firstsamplePN, _lastsamplePN);

  double chi2pn = 0;
  double ypnrange[50];
  double dsum = 0.;
  double dsum1 = 0.;
  double bl = 0.;
  double bl1 = 0.;
  double val_max = 0.;
  unsigned int samplemax = 0;
  unsigned int k;

  std::vector<double> allPNAmpl;

  for (EcalPnDiodeDigiCollection::const_iterator pnItr = PNDigi->begin(); pnItr != PNDigi->end();
       ++pnItr) {  // Loop on PNs

    for (int samId = 0; samId < (*pnItr).size(); samId++) {  // Loop on PN samples
      pn[samId] = (*pnItr).sample(samId).adc();
    }

    for (dsum = 0., k = 0; k < _presamplePN; k++) {
      dsum += pn[k];
    }
    bl = dsum / ((double)_presamplePN);

    for (val_max = 0., k = 0; k < _nsamplesPN; k++) {
      ypnrange[k] = pn[k] - bl;

      if (ypnrange[k] > val_max) {
        val_max = ypnrange[k];
        samplemax = k;
      }
    }

    chi2pn = pnfit->doFit(samplemax, &ypnrange[0]);

    if (chi2pn == 101 || chi2pn == 102 || chi2pn == 103)
      pnAmpl = 0.;
    else
      pnAmpl = pnfit->getAmpl();

    allPNAmpl.push_back(pnAmpl);
  }

  // ===========
  // Get Matacq
  // ===========

  ttrig = -1.;
  peak = -1.;

  if (doesRefFileExist == 1) {
    // FIXME
    if (color == 0)
      matacqTree->GetEntry(event - 1);
    else if (color == 3)
      matacqTree->GetEntry(matacqTree->GetEntries("color==0") + event - 1);
    ttrig = ttMat;
    peak = peakMat;
  }

  // ===========================
  // Decode EBDigis Information
  // ===========================

  double yrange[10];
  int adcGain = 0;
  int side = 0;

  if (EBDigi) {
    for (EBDigiCollection::const_iterator digiItr = EBDigi->begin(); digiItr != EBDigi->end();
         ++digiItr) {  // Loop on crystals

      EBDetId id_crystal(digiItr->id());
      EBDataFrame df(*digiItr);

      int etaG = id_crystal.ieta();  // global
      int phiG = id_crystal.iphi();  // global

      int etaL;  // local
      int phiL;  // local
      std::pair<int, int> LocalCoord = MEEBGeom::localCoord(etaG, phiG);

      etaL = LocalCoord.first;
      phiL = LocalCoord.second;

      eta = etaG;
      phi = phiG;

      side = MEEBGeom::side(etaG, phiG);

      EcalElectronicsId elecid_crystal = TheMapping.getElectronicsId(id_crystal);

      int towerID = elecid_crystal.towerId();
      // int channelID=elecid_crystal.channelId()-1;  // FIXME so far for endcap only
      int strip = elecid_crystal.stripId();
      int xtal = elecid_crystal.xtalId();
      int channelID = 5 * (strip - 1) + xtal - 1;  // FIXME

      int module = MEEBGeom::lmmod(etaG, phiG);

      std::pair<int, int> pnpair = MEEBGeom::pn(module);
      unsigned int MyPn0 = pnpair.first;
      unsigned int MyPn1 = pnpair.second;

      unsigned int channel = MEEBGeom::electronic_channel(etaL, phiL);
      assert(channel < nCrys);

      double adcmax = 0.0;

      if (towerID != int(_tower) || channelID != int(_channel) || dccID != int(_fedid - 600))
        continue;
      else
        channelNumber = channel;

      for (unsigned int i = 0; i < (*digiItr).size(); ++i) {  // Loop on adc samples

        EcalMGPASample samp_crystal(df.sample(i));
        adc[i] = samp_crystal.adc();
        adcG[i] = samp_crystal.gainId();

        if (i == 0)
          adcGain = adcG[i];
        if (i > 0)
          adcGain = TMath::Max(adcG[i], adcGain);

        if (adc[i] > adcmax) {
          adcmax = adc[i];
        }
      }

      for (dsum = 0., dsum1 = 0., k = 0; k < _presample; k++) {
        dsum += adc[k];
        if (k < _presample - 1)
          dsum1 += adc[k];
      }

      bl = dsum / ((double)_presample);
      bl1 = dsum1 / ((double)_presample - 1);

      for (val_max = 0., k = 0; k < _nsamples; k++) {
        yrange[k] = adc[k] - bl;
        if (yrange[k] > val_max) {
          val_max = yrange[k];
          samplemax = k;
        }
      }

      if (samplemax == 4 || samplemax == 5) {
        for (k = 0; k < _nsamples; k++) {
          yrange[k] = yrange[k] + bl - bl1;
        }
      }

      for (unsigned int k = 0; k < _nsamples; k++) {
        adc[k] = yrange[k];
      }

      pn0 = allPNAmpl[MyPn0];
      pn1 = allPNAmpl[MyPn1];

      if (samplemax >= _timingcutlow && samplemax <= _timingcuthigh && lightside == side)
        ADCtrees->Fill();
    }

  } else {
    for (EEDigiCollection::const_iterator digiItr = EEDigi->begin(); digiItr != EEDigi->end();
         ++digiItr) {  // Loop on crystals

      EEDetId id_crystal(digiItr->id());
      EEDataFrame df(*digiItr);

      phi = id_crystal.ix() - 1;
      eta = id_crystal.iy() - 1;
      side = 0;  // FIXME

      // Recover the TT id and the electronic crystal numbering from EcalElectronicsMapping

      EcalElectronicsId elecid_crystal = TheMapping.getElectronicsId(id_crystal);

      int towerID = elecid_crystal.towerId();
      int channelID = elecid_crystal.channelId() - 1;

      int module = MEEEGeom::lmmod(phi, eta);

      std::pair<int, int> pnpair = MEEEGeom::pn(module, _fedid);
      unsigned int MyPn0 = pnpair.first;
      unsigned int MyPn1 = pnpair.second;

      unsigned int channel = MEEEGeom::crystal(phi, eta);
      assert(channel < nCrys);

      double adcmax = 0.0;

      if (towerID != int(_tower) || channelID != int(_channel) || dccID != int(_fedid - 600))
        continue;
      else
        channelNumber = channel;

      for (unsigned int i = 0; i < (*digiItr).size(); ++i) {  // Loop on adc samples

        EcalMGPASample samp_crystal(df.sample(i));
        adc[i] = samp_crystal.adc();
        adcG[i] = samp_crystal.gainId();

        if (i == 0)
          adcGain = adcG[i];
        if (i > 0)
          adcGain = TMath::Max(adcG[i], adcGain);

        if (adc[i] > adcmax) {
          adcmax = adc[i];
        }
      }

      for (dsum = 0., dsum1 = 0., k = 0; k < _presample; k++) {
        dsum += adc[k];
        if (k < _presample - 1)
          dsum1 += adc[k];
      }

      bl = dsum / ((double)_presample);
      bl1 = dsum1 / ((double)_presample - 1);

      for (val_max = 0., k = 0; k < _nsamples; k++) {
        yrange[k] = adc[k] - bl;
        if (yrange[k] > val_max) {
          val_max = yrange[k];
          samplemax = k;
        }
      }

      if (samplemax == 4 || samplemax == 5) {
        for (k = 0; k < _nsamples; k++) {
          yrange[k] = yrange[k] + bl - bl1;
        }
      }

      for (unsigned int k = 0; k < _nsamples; k++) {
        adc[k] = yrange[k];
      }

      pn0 = allPNAmpl[MyPn0];
      pn1 = allPNAmpl[MyPn1];

      if (samplemax >= _timingcutlow && samplemax <= _timingcuthigh && lightside == side)
        ADCtrees->Fill();
    }
  }

}  // analyze

beginJob()

void EcalPerEvtLaserAnalyzer::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 104 of file EcalPerEvtLaserAnalyzer.cc.

References adc, ADCfile, ADCFile, ADCtrees, color, dccID, eta, IsFileCreated, Skims_PA_cff::name, peak, phi, pn0, pn1, resdir_, and ttrig.

                                       {
  //========================================================================

  // Define temporary files' names

  stringstream namefile1;
  namefile1 << resdir_ << "/ADCSamples.root";

  ADCfile = namefile1.str();

  // Create temporary file and trees to save adc samples

  ADCFile = new TFile(ADCfile.c_str(), "RECREATE");

  stringstream name;
  name << "ADCTree";

  ADCtrees = new TTree(name.str().c_str(), name.str().c_str());

  ADCtrees->Branch("iphi", &phi, "phi/I");
  ADCtrees->Branch("ieta", &eta, "eta/I");
  ADCtrees->Branch("dccID", &dccID, "dccID/I");
  ADCtrees->Branch("event", &event, "event/I");
  ADCtrees->Branch("color", &color, "color/I");
  ADCtrees->Branch("adc", &adc, "adc[10]/D");
  ADCtrees->Branch("ttrigMatacq", &ttrig, "ttrig/D");
  ADCtrees->Branch("peakMatacq", &peak, "peak/D");
  ADCtrees->Branch("pn0", &pn0, "pn0/D");
  ADCtrees->Branch("pn1", &pn1, "pn1/D");

  ADCtrees->SetBranchAddress("ieta", &eta);
  ADCtrees->SetBranchAddress("iphi", &phi);
  ADCtrees->SetBranchAddress("dccID", &dccID);
  ADCtrees->SetBranchAddress("event", &event);
  ADCtrees->SetBranchAddress("color", &color);
  ADCtrees->SetBranchAddress("adc", adc);
  ADCtrees->SetBranchAddress("ttrigMatacq", &ttrig);
  ADCtrees->SetBranchAddress("peakMatacq", &peak);
  ADCtrees->SetBranchAddress("pn0", &pn0);
  ADCtrees->SetBranchAddress("pn1", &pn1);

  IsFileCreated = 0;
}

endJob()

void EcalPerEvtLaserAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 533 of file EcalPerEvtLaserAnalyzer.cc.

References _channel, _fedid, _firstsample, _lastsample, _niter, _nsamples, _tower, adc, ADCfile, ADCFile, ADCtrees, simBeamSpotPI::alpha, apdAmpl, APDFile, apdTime, APDtrees, cms::cuda::assert(), HLT_2024v14_cff::beta, channelNumber, nano_mu_local_reco_cff::chi2, color, TkAlMuonSelectors_cfi::cut, dccID, PulseFitWithFunction::doFit(), RemoveAddSevLevel::flag, PulseFitWithFunction::getAmpl(), PulseFitWithFunction::getTime(), header, mps_fire::i, hcalRecHitTable_cff::ieta, PulseFitWithFunction::init(), hcalRecHitTable_cff::iphi, IsHeaderFilled, IsThereDataADC, hfnoseParametersInitialization_cfi::name2, nColor, peak, pn0, pn1, refalphabeta_, resfile, ecaldqm::towerID(), and ttrig.

Referenced by o2olib.O2ORunMgr::executeJob().

                                     {
  //========================================================================

  assert(colors.size() <= nColor);
  unsigned int nCol = colors.size();

  ADCtrees->Write();

  edm::LogVerbatim("EcalPerEvtLaserAnalyzer")
      << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
  edm::LogVerbatim("EcalPerEvtLaserAnalyzer")
      << "\t+=+       Analyzing laser data: getting per event               +=+";
  edm::LogVerbatim("EcalPerEvtLaserAnalyzer")
      << "\t+=+            APD Amplitudes and ADC samples                   +=+";
  edm::LogVerbatim("EcalPerEvtLaserAnalyzer")
      << "\t+=+    for fed:" << _fedid << ", tower:" << _tower << ", and channel:" << _channel;

  // Define temporary tree to save APD amplitudes

  APDFile = new TFile(resfile.c_str(), "RECREATE");

  int ieta, iphi, channelID, towerID, flag;
  double alpha, beta;

  colors.push_back(color);

  for (unsigned int i = 0; i < nCol; i++) {
    stringstream name1;
    name1 << "headerCol" << colors[i];

    header[i] = new TTree(name1.str().c_str(), name1.str().c_str());

    header[i]->Branch("alpha", &alpha, "alpha/D");
    header[i]->Branch("beta", &beta, "beta/D");
    header[i]->Branch("iphi", &iphi, "iphi/I");
    header[i]->Branch("ieta", &ieta, "ieta/I");
    header[i]->Branch("dccID", &dccID, "dccID/I");
    header[i]->Branch("towerID", &towerID, "towerID/I");
    header[i]->Branch("channelID", &channelID, "channelID/I");

    header[i]->SetBranchAddress("alpha", &alpha);
    header[i]->SetBranchAddress("beta", &beta);
    header[i]->SetBranchAddress("iphi", &iphi);
    header[i]->SetBranchAddress("ieta", &ieta);
    header[i]->SetBranchAddress("dccID", &dccID);
    header[i]->SetBranchAddress("towerID", &towerID);
    header[i]->SetBranchAddress("channelID", &channelID);
  }

  stringstream name2;
  name2 << "APDTree";
  APDtrees = new TTree(name2.str().c_str(), name2.str().c_str());

  //List of branches

  APDtrees->Branch("event", &event, "event/I");
  APDtrees->Branch("color", &color, "color/I");
  APDtrees->Branch("adc", &adc, "adc[10]/D");
  APDtrees->Branch("peakMatacq", &peak, "peak/D");
  APDtrees->Branch("ttrigMatacq", &ttrig, "ttrig/D");
  APDtrees->Branch("apdAmpl", &apdAmpl, "apdAmpl/D");
  APDtrees->Branch("apdTime", &apdTime, "apdTime/D");
  APDtrees->Branch("flag", &flag, "flag/I");
  APDtrees->Branch("pn0", &pn0, "pn0/D");
  APDtrees->Branch("pn1", &pn1, "pn1/D");

  APDtrees->SetBranchAddress("event", &event);
  APDtrees->SetBranchAddress("color", &color);
  APDtrees->SetBranchAddress("adc", adc);
  APDtrees->SetBranchAddress("peakMatacq", &peak);
  APDtrees->SetBranchAddress("ttrigMatacq", &ttrig);
  APDtrees->SetBranchAddress("apdAmpl", &apdAmpl);
  APDtrees->SetBranchAddress("apdTime", &apdTime);
  APDtrees->SetBranchAddress("flag", &flag);
  APDtrees->SetBranchAddress("pn0", &pn0);
  APDtrees->SetBranchAddress("pn1", &pn1);

  // Retrieve alpha and beta for APD amplitudes calculation

  TFile* alphaFile = new TFile(refalphabeta_.c_str());
  TTree* alphaTree[2];

  Double_t alphaRun, betaRun;
  int ietaRun, iphiRun, channelIDRun, towerIDRun, dccIDRun, flagRun;

  for (unsigned int i = 0; i < nCol; i++) {
    stringstream name3;
    name3 << "ABCol" << i;
    alphaTree[i] = (TTree*)alphaFile->Get(name3.str().c_str());
    alphaTree[i]->SetBranchStatus("*", false);
    alphaTree[i]->SetBranchStatus("alpha", true);
    alphaTree[i]->SetBranchStatus("beta", true);
    alphaTree[i]->SetBranchStatus("iphi", true);
    alphaTree[i]->SetBranchStatus("ieta", true);
    alphaTree[i]->SetBranchStatus("dccID", true);
    alphaTree[i]->SetBranchStatus("towerID", true);
    alphaTree[i]->SetBranchStatus("channelID", true);
    alphaTree[i]->SetBranchStatus("flag", true);

    alphaTree[i]->SetBranchAddress("alpha", &alphaRun);
    alphaTree[i]->SetBranchAddress("beta", &betaRun);
    alphaTree[i]->SetBranchAddress("iphi", &iphiRun);
    alphaTree[i]->SetBranchAddress("ieta", &ietaRun);
    alphaTree[i]->SetBranchAddress("dccID", &dccIDRun);
    alphaTree[i]->SetBranchAddress("towerID", &towerIDRun);
    alphaTree[i]->SetBranchAddress("channelID", &channelIDRun);
    alphaTree[i]->SetBranchAddress("flag", &flagRun);
  }

  PulseFitWithFunction* pslsfit = new PulseFitWithFunction();

  double chi2;

  for (unsigned int icol = 0; icol < nCol; icol++) {
    IsThereDataADC[icol] = 1;
    stringstream cut;
    cut << "color==" << colors.at(icol);
    if (ADCtrees->GetEntries(cut.str().c_str()) < 10)
      IsThereDataADC[icol] = 0;
    IsHeaderFilled[icol] = 0;
  }

  // Define submodule and channel number inside the submodule (as Patrice)

  for (Long64_t jentry = 0; jentry < ADCtrees->GetEntriesFast(); jentry++) {  // Loop on events
    ADCtrees->GetEntry(jentry);

    int iCry = channelNumber;

    // get back color

    unsigned int iCol = 0;
    for (unsigned int i = 0; i < nCol; i++) {
      if (color == colors[i]) {
        iCol = i;
        i = colors.size();
      }
    }

    alphaTree[iCol]->GetEntry(iCry);

    flag = flagRun;
    iphi = iphiRun;
    ieta = ietaRun;
    towerID = towerIDRun;
    channelID = channelIDRun;
    alpha = alphaRun;
    beta = betaRun;

    if (IsHeaderFilled[iCol] == 0) {
      header[iCol]->Fill();
      IsHeaderFilled[iCol] = 1;
    }
    // Amplitude calculation

    apdAmpl = 0;
    apdTime = 0;

    pslsfit->init(_nsamples, _firstsample, _lastsample, _niter, alphaRun, betaRun);
    chi2 = pslsfit->doFit(&adc[0]);

    if (chi2 < 0. || chi2 == 102 || chi2 == 101) {
      apdAmpl = 0;
      apdTime = 0;

    } else {
      apdAmpl = pslsfit->getAmpl();
      apdTime = pslsfit->getTime();
    }

    APDtrees->Fill();
  }

  alphaFile->Close();

  ADCFile->Close();

  APDFile->Write();
  APDFile->Close();

  // Remove unwanted files

  stringstream del;
  del << "rm " << ADCfile;
  system(del.str().c_str());

  edm::LogVerbatim("EcalPerEvtLaserAnalyzer")
      << "\t+=+    .................................................. done  +=+";
  edm::LogVerbatim("EcalPerEvtLaserAnalyzer")
      << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
}

Member Data Documentation

_channel

const unsigned int EcalPerEvtLaserAnalyzer::_channel

private

Definition at line 84 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

_ecalPart

const std::string EcalPerEvtLaserAnalyzer::_ecalPart

private

Definition at line 85 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

_fedid

const int EcalPerEvtLaserAnalyzer::_fedid

private

Definition at line 82 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

_firstsample

const unsigned int EcalPerEvtLaserAnalyzer::_firstsample

private

Definition at line 73 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

_firstsamplePN

const unsigned int EcalPerEvtLaserAnalyzer::_firstsamplePN

private

Definition at line 77 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

_lastsample

const unsigned int EcalPerEvtLaserAnalyzer::_lastsample

private

Definition at line 74 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

_lastsamplePN

const unsigned int EcalPerEvtLaserAnalyzer::_lastsamplePN

private

Definition at line 78 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

_niter

const unsigned int EcalPerEvtLaserAnalyzer::_niter

private

Definition at line 81 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

_nsamples

const unsigned int EcalPerEvtLaserAnalyzer::_nsamples

private

Definition at line 71 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

_nsamplesPN

const unsigned int EcalPerEvtLaserAnalyzer::_nsamplesPN

private

Definition at line 75 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

_presample

const unsigned int EcalPerEvtLaserAnalyzer::_presample

private

Definition at line 72 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

_presamplePN

const unsigned int EcalPerEvtLaserAnalyzer::_presamplePN

private

Definition at line 76 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

_timingcuthigh

const unsigned int EcalPerEvtLaserAnalyzer::_timingcuthigh

private

Definition at line 80 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

_timingcutlow

const unsigned int EcalPerEvtLaserAnalyzer::_timingcutlow

private

Definition at line 79 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

_tower

const unsigned int EcalPerEvtLaserAnalyzer::_tower

private

Definition at line 83 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

adc

double EcalPerEvtLaserAnalyzer::adc[10]

private

Definition at line 137 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

ADCfile

std::string EcalPerEvtLaserAnalyzer::ADCfile

private

Definition at line 92 of file EcalPerEvtLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

ADCFile

TFile* EcalPerEvtLaserAnalyzer::ADCFile

private

Definition at line 118 of file EcalPerEvtLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

adcG

int EcalPerEvtLaserAnalyzer::adcG[10]

private

Definition at line 138 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

ADCtrees

TTree* EcalPerEvtLaserAnalyzer::ADCtrees

private

Definition at line 119 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

apdAmpl

double EcalPerEvtLaserAnalyzer::apdAmpl

private

Definition at line 143 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

APDFile

TFile* EcalPerEvtLaserAnalyzer::APDFile

private

Definition at line 121 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

apdTime

double EcalPerEvtLaserAnalyzer::apdTime

private

Definition at line 144 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

APDtrees

TTree* EcalPerEvtLaserAnalyzer::APDtrees

private

Definition at line 123 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

channelNumber

int EcalPerEvtLaserAnalyzer::channelNumber

private

Definition at line 110 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

colMat

int EcalPerEvtLaserAnalyzer::colMat

private

Definition at line 130 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

color

int EcalPerEvtLaserAnalyzer::color

private

Definition at line 136 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), beginJob(), endJob(), and edmStreamStallGrapher.StreamInfoElement::unpack().

colors

std::vector<int> EcalPerEvtLaserAnalyzer::colors

private

Definition at line 112 of file EcalPerEvtLaserAnalyzer.h.

dccID

int EcalPerEvtLaserAnalyzer::dccID

private

Definition at line 106 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

digiCollection_

const std::string EcalPerEvtLaserAnalyzer::digiCollection_

private

Definition at line 59 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

digiPNCollection_

const std::string EcalPerEvtLaserAnalyzer::digiPNCollection_

private

Definition at line 61 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

digiProducer_

const std::string EcalPerEvtLaserAnalyzer::digiProducer_

private

Definition at line 60 of file EcalPerEvtLaserAnalyzer.h.

Referenced by EcalPerEvtLaserAnalyzer().

doesRefFileExist

int EcalPerEvtLaserAnalyzer::doesRefFileExist

private

Definition at line 125 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

ebDigiToken_

edm::EDGetTokenT<EBDigiCollection> EcalPerEvtLaserAnalyzer::ebDigiToken_

private

Definition at line 64 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

eeDigiToken_

edm::EDGetTokenT<EEDigiCollection> EcalPerEvtLaserAnalyzer::eeDigiToken_

private

Definition at line 65 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

eta

int EcalPerEvtLaserAnalyzer::eta

private

Definition at line 134 of file EcalPerEvtLaserAnalyzer.h.

Referenced by Particle.Particle::__str__(), analyze(), beginJob(), datamodel.Object::DeltaR(), Jet.Jet::jetID(), datamodel.Object::p4(), and Jet.Jet::puJetId().

event

int EcalPerEvtLaserAnalyzer::event

private

Definition at line 135 of file EcalPerEvtLaserAnalyzer.h.

Referenced by Types.EventID::cppID(), looper.Looper::process(), and core.AutoHandle.AutoHandle::product().

eventHeaderCollection_

const std::string EcalPerEvtLaserAnalyzer::eventHeaderCollection_

private

Definition at line 57 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

eventHeaderProducer_

const std::string EcalPerEvtLaserAnalyzer::eventHeaderProducer_

private

Definition at line 58 of file EcalPerEvtLaserAnalyzer.h.

evtMat

int EcalPerEvtLaserAnalyzer::evtMat

private

Definition at line 130 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

fedID

int EcalPerEvtLaserAnalyzer::fedID

private

Definition at line 107 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

header

TTree* EcalPerEvtLaserAnalyzer::header[2]

private

Definition at line 122 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

iEvent

int EcalPerEvtLaserAnalyzer::iEvent

private

Definition at line 55 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

IsFileCreated

int EcalPerEvtLaserAnalyzer::IsFileCreated

private

Definition at line 100 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and beginJob().

IsHeaderFilled

int EcalPerEvtLaserAnalyzer::IsHeaderFilled[nColor]

private

Definition at line 127 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

IsThereDataADC

int EcalPerEvtLaserAnalyzer::IsThereDataADC[nColor]

private

Definition at line 126 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

lightside

int EcalPerEvtLaserAnalyzer::lightside

private

Definition at line 108 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

mappingToken_

const edm::ESGetToken<EcalElectronicsMapping, EcalMappingRcd> EcalPerEvtLaserAnalyzer::mappingToken_

private

Definition at line 67 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

matacqFile

TFile* EcalPerEvtLaserAnalyzer::matacqFile

private

Definition at line 115 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

matacqTree

TTree* EcalPerEvtLaserAnalyzer::matacqTree

private

Definition at line 116 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

nCrys

unsigned int EcalPerEvtLaserAnalyzer::nCrys

private

Definition at line 98 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

peak

double EcalPerEvtLaserAnalyzer::peak

private

Definition at line 129 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

peakMat

double EcalPerEvtLaserAnalyzer::peakMat

private

Definition at line 129 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

phi

int EcalPerEvtLaserAnalyzer::phi

private

Definition at line 134 of file EcalPerEvtLaserAnalyzer.h.

Referenced by Particle.Particle::__str__(), analyze(), beginJob(), datamodel.Object::DeltaR(), datamodel.Object::p4(), and ntupleDataFormat.Track::phiPull().

pn

double EcalPerEvtLaserAnalyzer::pn[50]

private

Definition at line 142 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

pn0

double EcalPerEvtLaserAnalyzer::pn0

private

Definition at line 141 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

pn1

double EcalPerEvtLaserAnalyzer::pn1

private

Definition at line 141 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

pnAmpl

double EcalPerEvtLaserAnalyzer::pnAmpl

private

Definition at line 145 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

pnDiodeDigiToken_

const edm::EDGetTokenT<EcalPnDiodeDigiCollection> EcalPerEvtLaserAnalyzer::pnDiodeDigiToken_

private

Definition at line 66 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

rawDataToken_

const edm::EDGetTokenT<EcalRawDataCollection> EcalPerEvtLaserAnalyzer::rawDataToken_

private

Definition at line 63 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

refalphabeta_

const std::string EcalPerEvtLaserAnalyzer::refalphabeta_

private

Definition at line 88 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

resdir_

const std::string EcalPerEvtLaserAnalyzer::resdir_

private

Definition at line 87 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and beginJob().

resfile

std::string EcalPerEvtLaserAnalyzer::resfile

private

Definition at line 93 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

runNum

int EcalPerEvtLaserAnalyzer::runNum

private

Definition at line 105 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

runType

int EcalPerEvtLaserAnalyzer::runType

private

Definition at line 104 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

tt

double EcalPerEvtLaserAnalyzer::tt

private

Definition at line 139 of file EcalPerEvtLaserAnalyzer.h.

ttMat

double EcalPerEvtLaserAnalyzer::ttMat

private

Definition at line 129 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

ttrig

double EcalPerEvtLaserAnalyzer::ttrig

private

Definition at line 140 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().