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 final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const 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
bool	wantsStreamRuns () const
	~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
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > 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
bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices 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 ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
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 () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
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

enum EcalPerEvtLaserAnalyzer::VarCol

Enumerator
iBlue
iRed
nColor

Definition at line 52 of file EcalPerEvtLaserAnalyzer.h.

{ iBlue, iRed, nColor };

Constructor & Destructor Documentation

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

explicit

Definition at line 36 of file EcalPerEvtLaserAnalyzer.cc.

References _ecalPart, digiCollection_, digiProducer_, ebDigiToken_, eeDigiToken_, HLT_FULL_cff::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 ( )

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

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(), edm::DataFrameContainer::begin(), EcalElectronicsId::channelId(), channelNumber, colMat, color, colors, MEEEGeom::crystal(), dccID, isotrackApplyRegressor::df, digiCollection_, digiPNCollection_, doesRefFileExist, TPNFit::doFit(), ebDigiToken_, eeDigiToken_, MEEBGeom::electronic_channel(), edm::SortedCollection< T, SORT >::end(), edm::DataFrameContainer::end(), eta, event, eventHeaderCollection_, evtMat, fedID, spr::find(), TPNFit::getAmpl(), edm::Event::getByToken(), edm::EventSetup::getData(), mps_fire::i, iEvent, TPNFit::init(), IsFileCreated, edm::HandleBase::isValid(), isotrackApplyRegressor::k, EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_STD, lightside, MEEEGeom::lmmod(), MEEBGeom::lmmod(), MEEBGeom::localCoord(), mappingToken_, matacqFile, matacqTree, Max(), callgraph::module, nCrys, peak, peakMat, phi, MEEBGeom::pn(), MEEEGeom::pn(), pn, pn0, pn1, pnAmpl, pnDiodeDigiToken_, edm::Handle< T >::product(), rawDataToken_, resdir_, resfile, runNum, runType, MEEBGeom::side(), digitizers_cfi::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) {
        val_max = val_max + bl - bl1;
        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) {
        val_max = val_max + bl - bl1;
        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

void EcalPerEvtLaserAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 104 of file EcalPerEvtLaserAnalyzer.cc.

References adc, ADCfile, ADCFile, ADCtrees, color, dccID, eta, event, IsFileCreated, mergeVDriftHistosByStation::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;
}

void EcalPerEvtLaserAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 535 of file EcalPerEvtLaserAnalyzer.cc.

References _channel, _fedid, _firstsample, _lastsample, _niter, _nsamples, _tower, adc, ADCfile, ADCFile, ADCtrees, alpha, apdAmpl, APDFile, apdTime, APDtrees, cms::cuda::assert(), HLT_FULL_cff::beta, channelNumber, HLT_FULL_cff::chi2, color, colors, GOODCOLL_filter_cfg::cut, dccID, PulseFitWithFunction::doFit(), event, PulseFitWithFunction::getAmpl(), PulseFitWithFunction::getTime(), header, mps_fire::i, PulseFitWithFunction::init(), IsHeaderFilled, IsThereDataADC, 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)

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

    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

const unsigned int EcalPerEvtLaserAnalyzer::_channel

private

Definition at line 84 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

const std::string EcalPerEvtLaserAnalyzer::_ecalPart

private

Definition at line 85 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

const int EcalPerEvtLaserAnalyzer::_fedid

private

Definition at line 82 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

const unsigned int EcalPerEvtLaserAnalyzer::_firstsample

private

Definition at line 73 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

const unsigned int EcalPerEvtLaserAnalyzer::_firstsamplePN

private

Definition at line 77 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

const unsigned int EcalPerEvtLaserAnalyzer::_lastsample

private

Definition at line 74 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

const unsigned int EcalPerEvtLaserAnalyzer::_lastsamplePN

private

Definition at line 78 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

const unsigned int EcalPerEvtLaserAnalyzer::_niter

private

Definition at line 81 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

const unsigned int EcalPerEvtLaserAnalyzer::_nsamples

private

Definition at line 71 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

const unsigned int EcalPerEvtLaserAnalyzer::_nsamplesPN

private

Definition at line 75 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

const unsigned int EcalPerEvtLaserAnalyzer::_presample

private

Definition at line 72 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

const unsigned int EcalPerEvtLaserAnalyzer::_presamplePN

private

Definition at line 76 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

const unsigned int EcalPerEvtLaserAnalyzer::_timingcuthigh

private

Definition at line 80 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

const unsigned int EcalPerEvtLaserAnalyzer::_timingcutlow

private

Definition at line 79 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

const unsigned int EcalPerEvtLaserAnalyzer::_tower

private

Definition at line 83 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

double EcalPerEvtLaserAnalyzer::adc[10]

private

Definition at line 137 of file EcalPerEvtLaserAnalyzer.h.

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

std::string EcalPerEvtLaserAnalyzer::ADCfile

private

Definition at line 92 of file EcalPerEvtLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

TFile* EcalPerEvtLaserAnalyzer::ADCFile

private

Definition at line 118 of file EcalPerEvtLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

int EcalPerEvtLaserAnalyzer::adcG[10]

private

Definition at line 138 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

TTree* EcalPerEvtLaserAnalyzer::ADCtrees

private

Definition at line 119 of file EcalPerEvtLaserAnalyzer.h.

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

double EcalPerEvtLaserAnalyzer::apdAmpl

private

Definition at line 143 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

TFile* EcalPerEvtLaserAnalyzer::APDFile

private

Definition at line 121 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

double EcalPerEvtLaserAnalyzer::apdTime

private

Definition at line 144 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

TTree* EcalPerEvtLaserAnalyzer::APDtrees

private

Definition at line 123 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

int EcalPerEvtLaserAnalyzer::channelNumber

private

Definition at line 110 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

int EcalPerEvtLaserAnalyzer::colMat

private

Definition at line 130 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

int EcalPerEvtLaserAnalyzer::color

private

Definition at line 136 of file EcalPerEvtLaserAnalyzer.h.

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

std::vector<int> EcalPerEvtLaserAnalyzer::colors

private

Definition at line 112 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

int EcalPerEvtLaserAnalyzer::dccID

private

Definition at line 106 of file EcalPerEvtLaserAnalyzer.h.

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

const std::string EcalPerEvtLaserAnalyzer::digiCollection_

private

Definition at line 59 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

const std::string EcalPerEvtLaserAnalyzer::digiPNCollection_

private

Definition at line 61 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

const std::string EcalPerEvtLaserAnalyzer::digiProducer_

private

Definition at line 60 of file EcalPerEvtLaserAnalyzer.h.

Referenced by EcalPerEvtLaserAnalyzer().

int EcalPerEvtLaserAnalyzer::doesRefFileExist

private

Definition at line 125 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<EBDigiCollection> EcalPerEvtLaserAnalyzer::ebDigiToken_

private

Definition at line 64 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

edm::EDGetTokenT<EEDigiCollection> EcalPerEvtLaserAnalyzer::eeDigiToken_

private

Definition at line 65 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

int EcalPerEvtLaserAnalyzer::eta

private

Definition at line 134 of file EcalPerEvtLaserAnalyzer.h.

Referenced by Particle.Particle::__str__(), analyze(), beginJob(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

int EcalPerEvtLaserAnalyzer::event

private

Definition at line 135 of file EcalPerEvtLaserAnalyzer.h.

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

const std::string EcalPerEvtLaserAnalyzer::eventHeaderCollection_

private

Definition at line 57 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

const std::string EcalPerEvtLaserAnalyzer::eventHeaderProducer_

private

Definition at line 58 of file EcalPerEvtLaserAnalyzer.h.

int EcalPerEvtLaserAnalyzer::evtMat

private

Definition at line 130 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

int EcalPerEvtLaserAnalyzer::fedID

private

Definition at line 107 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

TTree* EcalPerEvtLaserAnalyzer::header[2]

private

Definition at line 122 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

int EcalPerEvtLaserAnalyzer::iEvent

private

Definition at line 55 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

int EcalPerEvtLaserAnalyzer::IsFileCreated

private

Definition at line 100 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and beginJob().

int EcalPerEvtLaserAnalyzer::IsHeaderFilled[nColor]

private

Definition at line 127 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

int EcalPerEvtLaserAnalyzer::IsThereDataADC[nColor]

private

Definition at line 126 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

int EcalPerEvtLaserAnalyzer::lightside

private

Definition at line 108 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 67 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

TFile* EcalPerEvtLaserAnalyzer::matacqFile

private

Definition at line 115 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

TTree* EcalPerEvtLaserAnalyzer::matacqTree

private

Definition at line 116 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

unsigned int EcalPerEvtLaserAnalyzer::nCrys

private

Definition at line 98 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

double EcalPerEvtLaserAnalyzer::peak

private

Definition at line 129 of file EcalPerEvtLaserAnalyzer.h.

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

double EcalPerEvtLaserAnalyzer::peakMat

private

Definition at line 129 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

int EcalPerEvtLaserAnalyzer::phi

private

Definition at line 134 of file EcalPerEvtLaserAnalyzer.h.

Referenced by Particle.Particle::__str__(), analyze(), beginJob(), and ntupleDataFormat.Track::phiPull().

double EcalPerEvtLaserAnalyzer::pn[50]

private

Definition at line 142 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

double EcalPerEvtLaserAnalyzer::pn0

private

Definition at line 141 of file EcalPerEvtLaserAnalyzer.h.

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

double EcalPerEvtLaserAnalyzer::pn1

private

Definition at line 141 of file EcalPerEvtLaserAnalyzer.h.

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

double EcalPerEvtLaserAnalyzer::pnAmpl

private

Definition at line 145 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 66 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 63 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

const std::string EcalPerEvtLaserAnalyzer::refalphabeta_

private

Definition at line 88 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

const std::string EcalPerEvtLaserAnalyzer::resdir_

private

Definition at line 87 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and beginJob().

std::string EcalPerEvtLaserAnalyzer::resfile

private

Definition at line 93 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

int EcalPerEvtLaserAnalyzer::runNum

private

Definition at line 105 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

int EcalPerEvtLaserAnalyzer::runType

private

Definition at line 104 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

double EcalPerEvtLaserAnalyzer::tt

private

Definition at line 139 of file EcalPerEvtLaserAnalyzer.h.

double EcalPerEvtLaserAnalyzer::ttMat

private

Definition at line 129 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

double EcalPerEvtLaserAnalyzer::ttrig

private

Definition at line 140 of file EcalPerEvtLaserAnalyzer.h.

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