EcalLaserAnalyzer Class Reference

#include <EcalLaserAnalyzer.h>

Inheritance diagram for EcalLaserAnalyzer:

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
	EcalLaserAnalyzer (const edm::ParameterSet &iConfig)
void	endJob () override
void	setGeomEB (int etaG, int phiG, int module, int tower, int strip, int xtal, int apdRefTT, int channel, int lmr)
void	setGeomEE (int etaG, int phiG, int iX, int iY, int iZ, int module, int tower, int ch, int apdRefTT, int channel, int lmr)
	~EcalLaserAnalyzer () 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
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 double	_alpha
const double	_beta
const double	_chi2cut
const int	_debug
const bool	_docorpn
const std::string	_ecalPart
const int	_fedid
const unsigned int	_firstsample
const unsigned int	_firstsamplePN
bool	_fitab
const unsigned int	_lastsample
const unsigned int	_lastsamplePN
const unsigned int	_nevtmax
const unsigned int	_niter
const double	_noise
const unsigned int	_nsamples
const unsigned int	_nsamplesPN
unsigned int	_presample
const double	_presamplecut
const unsigned int	_presamplePN
const double	_qualpercent
const double	_ratiomaxcutlow
const double	_ratiomincuthigh
const double	_ratiomincutlow
const bool	_saveallevents
const unsigned int	_timingcuthigh
const unsigned int	_timingcutlow
const unsigned int	_timingqualhigh
const unsigned int	_timingquallow
double	adc [10]
std::string	ADCfile
TFile *	ADCFile
int	adcG [10]
double *	adcNoPed
TTree *	ADCtrees [1700]
std::string	alphafile
std::string	alphainitfile
double	APD [6]
double	apdAmpl
double	apdAmplA
double	apdAmplB
TAPD *	APDAnal [1700][nColor]
std::string	APDfile
TFile *	APDFile
TAPD *	APDFirstAnal [1700][nColor]
double	APDoAPDA [6]
double	APDoAPDB [6]
double	APDoPN [6]
double	APDoPNA [6]
double	APDoPNB [6]
TAPDPulse *	APDPulse
std::map< int, unsigned int >	apdRefMap [2]
double	apdTime
TTree *	APDtrees [1700]
int	channelID
int	channelIteratorEE
std::map< unsigned int, unsigned int >	channelMapEE
int	color
int	colorref
std::vector< int >	colors
int	dccID
TMom *	Delta01
TMom *	Delta12
const std::string	digiCollection_
const std::string	digiPNCollection_
const std::string	digiProducer_
bool	doesABTreeExist
edm::EDGetTokenT< EBDigiCollection >	ebDigiToken_
edm::EDGetTokenT< EEDigiCollection >	eeDigiToken_
int	eta
int	event
const std::string	eventHeaderCollection_
const std::string	eventHeaderProducer_
int	eventref
int	fedID
unsigned int	firstChanMod [22]
int	flag
int	flagAB
int	iChannelID [1700]
int	idccID [1700]
int	iEta [1700]
int	iEvent
unsigned int	iModule [1700]
int	iPhi [1700]
unsigned int	isFirstChanModFilled [22]
bool	isGainOK
int	iside [1700]
int	IsThereDataADC [1700][nColor]
bool	isTimingOK
int	iTowerID [1700]
int	iZ
int	laserEvents
int	lightside
const edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	mappingToken_
TMem *	Mem
int	moduleID
std::vector< int >	modules
unsigned int	nCrys
unsigned int	nevtAB [1700]
int	nEvtBadGain [1700]
int	nEvtBadTiming [1700]
int	nEvtTot [1700]
unsigned int	nMod
unsigned int	nPNPerMod
unsigned int	nRefChan
unsigned int	nRefTrees
unsigned int	nSides
int	phi
double	pn [50]
double	PN [6]
double	pn0
double	pn1
double	pnAmpl
TPN *	PNAnal [22][2][nColor]
const std::string	pncorfile_
TPNCor *	pnCorrector
const edm::EDGetTokenT< EcalPnDiodeDigiCollection >	pnDiodeDigiToken_
TPN *	PNFirstAnal [22][2][nColor]
int	pnG [50]
int	pnID
double *	pnNoPed
double	PNoPN [6]
double	PNoPNA [6]
double	PNoPNB [6]
TPNPulse *	PNPulse
const edm::EDGetTokenT< EcalRawDataCollection >	rawDataToken_
TTree *	RefAPDtrees [2][22]
const std::string	resdir_
std::string	resfile
TFile *	resFile
TTree *	respntrees [nColor]
TTree *	restrees [nColor]
int	runNum
int	runType
TShapeAnalysis *	shapana
int	side
double	Time [6]
int	towerID
bool	wasABCalcOK [1700]
bool	wasGainOK [1700]
bool	wasTimingOK [1700]

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 54 of file EcalLaserAnalyzer.h.

Member Enumeration Documentation

VarCol

enum EcalLaserAnalyzer::VarCol

Enumerator
iBlue
iRed
nColor

Definition at line 67 of file EcalLaserAnalyzer.h.

{ iBlue, iRed, nColor };

Constructor & Destructor Documentation

EcalLaserAnalyzer()

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

explicit

Definition at line 46 of file EcalLaserAnalyzer.cc.

References _ecalPart, _fedid, _firstsample, _lastsample, _nsamples, _nsamplesPN, _presample, _presamplePN, _ratiomaxcutlow, _ratiomincuthigh, _ratiomincutlow, _timingcuthigh, _timingcutlow, _timingqualhigh, _timingquallow, APDPulse, Delta01, Delta12, digiCollection_, digiProducer_, ebDigiToken_, eeDigiToken_, firstChanMod, iChannelID, idccID, iEta, cuy::ii, iModule, ProducerED_cfi::InputTag, iPhi, isFirstChanModFilled, isGainOK, iside, isTimingOK, iTowerID, iZ, dqmiolumiharvest::j, ME::lmmodFromDcc(), Mem, nCrys, NCRYSEB, NCRYSEE, nMod, NREFCHAN, nRefChan, pncorfile_, pnCorrector, PNPulse, wasGainOK, and wasTimingOK.

    : 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_))),
      mappingToken_(esConsumes()),
      // Framework parameters with default values
      _nsamples(iConfig.getUntrackedParameter<unsigned int>("nSamples", 10)),
      _presample(iConfig.getUntrackedParameter<unsigned int>("nPresamples", 2)),
      _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", 2)),
      _timingcuthigh(iConfig.getUntrackedParameter<unsigned int>("timingCutHigh", 9)),
      _timingquallow(iConfig.getUntrackedParameter<unsigned int>("timingQualLow", 3)),
      _timingqualhigh(iConfig.getUntrackedParameter<unsigned int>("timingQualHigh", 8)),
      _ratiomincutlow(iConfig.getUntrackedParameter<double>("ratioMinCutLow", 0.4)),
      _ratiomincuthigh(iConfig.getUntrackedParameter<double>("ratioMinCutHigh", 0.95)),
      _ratiomaxcutlow(iConfig.getUntrackedParameter<double>("ratioMaxCutLow", 0.8)),
      _presamplecut(iConfig.getUntrackedParameter<double>("presampleCut", 5.0)),
      _niter(iConfig.getUntrackedParameter<unsigned int>("nIter", 3)),
      _fitab(iConfig.getUntrackedParameter<bool>("fitAB", false)),
      _alpha(iConfig.getUntrackedParameter<double>("alpha", 1.5076494)),
      _beta(iConfig.getUntrackedParameter<double>("beta", 1.5136036)),
      _nevtmax(iConfig.getUntrackedParameter<unsigned int>("nEvtMax", 200)),
      _noise(iConfig.getUntrackedParameter<double>("noise", 2.0)),
      _chi2cut(iConfig.getUntrackedParameter<double>("chi2cut", 10.0)),
      _ecalPart(iConfig.getUntrackedParameter<std::string>("ecalPart", "EB")),
      _docorpn(iConfig.getUntrackedParameter<bool>("doCorPN", false)),
      _fedid(iConfig.getUntrackedParameter<int>("fedID", -999)),
      _saveallevents(iConfig.getUntrackedParameter<bool>("saveAllEvents", false)),
      _qualpercent(iConfig.getUntrackedParameter<double>("qualPercent", 0.2)),
      _debug(iConfig.getUntrackedParameter<int>("debug", 0)),
      resdir_(iConfig.getUntrackedParameter<std::string>("resDir")),
      pncorfile_(iConfig.getUntrackedParameter<std::string>("pnCorFile")),
      nCrys(NCRYSEB),
      nPNPerMod(NPNPERMOD),
      nMod(NMODEE),
      nSides(NSIDES),
      runType(-1),
      runNum(0),
      fedID(-1),
      dccID(-1),
      side(2),
      lightside(2),
      iZ(1),
      phi(-1),
      eta(-1),
      event(0),
      color(-1),
      pn0(0),
      pn1(0),
      apdAmpl(0),
      apdAmplA(0),
      apdAmplB(0),
      apdTime(0),
      pnAmpl(0),
      pnID(-1),
      moduleID(-1),
      channelIteratorEE(0)

//========================================================================

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

  // Geometrical constants initialization

  if (_ecalPart == "EB") {
    nCrys = NCRYSEB;
  } else {
    nCrys = NCRYSEE;
  }
  iZ = 1;
  if (_fedid <= 609)
    iZ = -1;
  modules = ME::lmmodFromDcc(_fedid);
  nMod = modules.size();
  nRefChan = NREFCHAN;

  for (unsigned int j = 0; j < nCrys; j++) {
    iEta[j] = -1;
    iPhi[j] = -1;
    iModule[j] = 10;
    iTowerID[j] = -1;
    iChannelID[j] = -1;
    idccID[j] = -1;
    iside[j] = -1;
    wasTimingOK[j] = true;
    wasGainOK[j] = true;
  }

  for (unsigned int j = 0; j < nMod; j++) {
    int ii = modules[j];
    firstChanMod[ii - 1] = 0;
    isFirstChanModFilled[ii - 1] = 0;
  }

  // Quality check flags

  isGainOK = true;
  isTimingOK = true;

  // PN linearity corrector

  pnCorrector = new TPNCor(pncorfile_);

  // Objects dealing with pulses

  APDPulse = new TAPDPulse(_nsamples,
                           _presample,
                           _firstsample,
                           _lastsample,
                           _timingcutlow,
                           _timingcuthigh,
                           _timingquallow,
                           _timingqualhigh,
                           _ratiomincutlow,
                           _ratiomincuthigh,
                           _ratiomaxcutlow);
  PNPulse = new TPNPulse(_nsamplesPN, _presamplePN);

  // Object dealing with MEM numbering

  Mem = new TMem(_fedid);

  // Objects needed for npresample calculation

  Delta01 = new TMom();
  Delta12 = new TMom();
}

~EcalLaserAnalyzer()

EcalLaserAnalyzer::~EcalLaserAnalyzer ( )

override

Definition at line 191 of file EcalLaserAnalyzer.cc.

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

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

Member Function Documentation

analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 307 of file EcalLaserAnalyzer.cc.

References _debug, _docorpn, _ecalPart, _fedid, _firstsamplePN, _fitab, _lastsamplePN, _nevtmax, _nsamplesPN, adc, adcG, ADCtrees, TMom::addEntry(), APDPulse, MEEBGeom::apdRefTower(), MEEEGeom::apdRefTower(), cms::cuda::assert(), edm::SortedCollection< T, SORT >::begin(), HltBtagPostValidation_cff::c, EcalElectronicsId::channelId(), channelID, channelIteratorEE, channelMapEE, color, alignCSCRings::corr, dccID, MEEEGeom::dee(), Delta01, Delta12, hgcalPerformanceValidation::df, digiCollection_, digiPNCollection_, doesABTreeExist, TPNFit::doFit(), MillePedeFileConverter_cfg::e, CollectionTags_cfi::EBDigi, ebDigiToken_, CollectionTags_cfi::EEDigi, eeDigiToken_, MEEBGeom::electronic_channel(), edm::SortedCollection< T, SORT >::end(), eta, eventHeaderCollection_, eventHeaderProducer_, fedID, spr::find(), TPNPulse::getAdcWithoutPedestal(), TPNFit::getAmpl(), TAPDPulse::getDelta(), TPNPulse::getMaxSample(), TPNCor::getPNCorrectionFactor(), ecalpyutils::hashedIndex(), mps_fire::i, EcalPnDiodeDetId::iDCCId(), iEvent, TPNFit::init(), createfilelist::int, EcalPnDiodeDetId::iPnId(), TMem::isMemRelevant(), TAPDPulse::isPulseOK(), TAPDPulse::isTimingQualOK(), edm::HandleBase::isValid(), iZ, EcalDCCHeaderBlock::LASER_DELAY_SCAN, EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_POWER_SCAN, EcalDCCHeaderBlock::LASER_STD, laserEvents, lightside, MEEEGeom::lmmod(), MEEBGeom::lmmod(), MEEBGeom::lmr(), MEEEGeom::lmr(), MEEBGeom::localCoord(), mappingToken_, METSkim_cff::Max, TMem::Mem(), Mem, nCrys, nevtAB, nEvtBadGain, nEvtBadTiming, nEvtTot, phi, MEEBGeom::pn(), MEEEGeom::pn(), pn, pn0, pn1, pnAmpl, pnCorrector, pnDiodeDigiToken_, pnG, pnNoPed, PNPulse, edm::Handle< T >::product(), TShapeAnalysis::putAllVals(), rawDataToken_, runNum, runType, setGeomEB(), setGeomEE(), TPNPulse::setPulse(), TAPDPulse::setPulse(), shapana, side, findQualityFiles::size, nano_mu_digi_cff::strip, EcalElectronicsId::stripId(), l1tHGCalTowerProducer_cfi::tower, EcalElectronicsId::towerId(), towerID, 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()
                            << " " << eventHeaderProducer_.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 if (_ecalPart == "EE") {
    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();
    }
  } else {
    edm::LogError("cfg_error") << " Wrong ecalPart in cfg file ";
    return;
  }

  // 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 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();
    lightside = headerItr->getRtHalf();

    // Check fed corresponds to the DCC in TCC

    if (600 + dccID != fedID)
      continue;

    // Cut on runType

    if (runType != EcalDCCHeaderBlock::LASER_STD && runType != EcalDCCHeaderBlock::LASER_GAP &&
        runType != EcalDCCHeaderBlock::LASER_POWER_SCAN && runType != EcalDCCHeaderBlock::LASER_DELAY_SCAN)
      return;

    // Retrieve laser color and event number

    EcalDCCHeaderBlock::EcalDCCEventSettings settings = headerItr->getEventSettings();
    color = settings.wavelength;
    if (color < 0)
      return;

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

  // Cut on fedID

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

  // Count laser events

  laserEvents++;

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

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

  double chi2pn = 0;
  unsigned int samplemax = 0;
  int pnGain = 0;

  std::map<int, std::vector<double> > allPNAmpl;
  std::map<int, std::vector<double> > allPNGain;

  // Loop on PNs digis

  for (EcalPnDiodeDigiCollection::const_iterator pnItr = PNDigi->begin(); pnItr != PNDigi->end(); ++pnItr) {
    EcalPnDiodeDetId pnDetId = EcalPnDiodeDetId((*pnItr).id());

    if (_debug == 1)
      edm::LogVerbatim("EcalLaserAnalyzer")
          << "-- debug -- Inside PNDigi - pnID=" << pnDetId.iPnId() << ", dccID=" << pnDetId.iDCCId();

    // Skip MEM DCC without relevant data

    bool isMemRelevant = Mem->isMemRelevant(pnDetId.iDCCId());
    if (!isMemRelevant)
      continue;

    // Loop on PN samples

    for (int samId = 0; samId < (*pnItr).size(); samId++) {
      pn[samId] = (*pnItr).sample(samId).adc();
      pnG[samId] = (*pnItr).sample(samId).gainId();
      if (samId == 0)
        pnGain = pnG[samId];
      if (samId > 0)
        pnGain = int(TMath::Max(pnG[samId], pnGain));
    }

    if (pnGain != 1)
      edm::LogVerbatim("EcalLaserAnalyzer") << "PN gain different from 1";

    // Calculate amplitude from pulse

    PNPulse->setPulse(pn);
    pnNoPed = PNPulse->getAdcWithoutPedestal();
    samplemax = PNPulse->getMaxSample();
    chi2pn = pnfit->doFit(samplemax, &pnNoPed[0]);
    if (chi2pn == 101 || chi2pn == 102 || chi2pn == 103)
      pnAmpl = 0.;
    else
      pnAmpl = pnfit->getAmpl();

    // Apply linearity correction

    double corr = 1.0;
    if (_docorpn)
      corr = pnCorrector->getPNCorrectionFactor(pnAmpl, pnGain);
    pnAmpl *= corr;

    // Fill PN ampl vector

    allPNAmpl[pnDetId.iDCCId()].push_back(pnAmpl);

    if (_debug == 1)
      edm::LogVerbatim("EcalLaserAnalyzer") << "-- debug -- Inside PNDigi - PNampl=" << pnAmpl << ", PNgain=" << pnGain;
  }

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

  int adcGain = 0;

  if (EBDigi) {
    // Loop on crystals
    //===================

    for (EBDigiCollection::const_iterator digiItr = EBDigi->begin(); digiItr != EBDigi->end(); ++digiItr) {
      // Retrieve geometry
      //===================

      EBDetId id_crystal(digiItr->id());
      EBDataFrame df(*digiItr);
      EcalElectronicsId elecid_crystal = TheMapping.getElectronicsId(id_crystal);

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

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

      int etaL = LocalCoord.first;   // local
      int phiL = LocalCoord.second;  // local

      int strip = elecid_crystal.stripId();
      int xtal = elecid_crystal.xtalId();

      int module = MEEBGeom::lmmod(etaG, phiG);
      int tower = elecid_crystal.towerId();

      int apdRefTT = MEEBGeom::apdRefTower(module);

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

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

      setGeomEB(etaG, phiG, module, tower, strip, xtal, apdRefTT, channel, lmr);

      if (_debug == 1)
        edm::LogVerbatim("EcalLaserAnalyzer")
            << "-- debug -- Inside EBDigi - towerID:" << towerID << " channelID:" << channelID << " module:" << module
            << " modules:" << modules.size();

      // APD Pulse
      //===========

      // Loop on adc samples

      for (unsigned int i = 0; i < (*digiItr).size(); ++i) {
        EcalMGPASample samp_crystal(df.sample(i));
        adc[i] = samp_crystal.adc();
        adcG[i] = samp_crystal.gainId();
        adc[i] *= adcG[i];
        if (i == 0)
          adcGain = adcG[i];
        if (i > 0)
          adcGain = TMath::Max(adcG[i], adcGain);
      }

      APDPulse->setPulse(adc);

      // Quality checks
      //================

      if (adcGain != 1)
        nEvtBadGain[channel]++;
      if (!APDPulse->isTimingQualOK())
        nEvtBadTiming[channel]++;
      nEvtTot[channel]++;

      // Associate PN ampl
      //===================

      int mem0 = Mem->Mem(lmr, 0);
      int mem1 = Mem->Mem(lmr, 1);

      if (allPNAmpl[mem0].size() > MyPn0)
        pn0 = allPNAmpl[mem0][MyPn0];
      else
        pn0 = 0;
      if (allPNAmpl[mem1].size() > MyPn1)
        pn1 = allPNAmpl[mem1][MyPn1];
      else
        pn1 = 0;

      // Fill if Pulse is fine
      //=======================

      if (APDPulse->isPulseOK() && lightside == side) {
        ADCtrees[channel]->Fill();

        Delta01->addEntry(APDPulse->getDelta(0, 1));
        Delta12->addEntry(APDPulse->getDelta(1, 2));

        if (nevtAB[channel] < _nevtmax && _fitab) {
          if (doesABTreeExist)
            shapana->putAllVals(channel, adc, eta, phi);
          else
            shapana->putAllVals(channel, adc, eta, phi, dccID, side, towerID, channelID);
          nevtAB[channel]++;
        }
      }
    }

  } else if (EEDigi) {
    // Loop on crystals
    //===================

    for (EEDigiCollection::const_iterator digiItr = EEDigi->begin(); digiItr != EEDigi->end(); ++digiItr) {
      // Retrieve geometry
      //===================

      EEDetId id_crystal(digiItr->id());
      EEDataFrame df(*digiItr);
      EcalElectronicsId elecid_crystal = TheMapping.getElectronicsId(id_crystal);

      int etaG = id_crystal.iy();
      int phiG = id_crystal.ix();

      int iX = (phiG - 1) / 5 + 1;
      int iY = (etaG - 1) / 5 + 1;

      int tower = elecid_crystal.towerId();
      int ch = elecid_crystal.channelId() - 1;

      int module = MEEEGeom::lmmod(iX, iY);
      if (module >= 18 && side == 1)
        module += 2;
      int lmr = MEEEGeom::lmr(iX, iY, iZ);
      int dee = MEEEGeom::dee(lmr);
      int apdRefTT = MEEEGeom::apdRefTower(lmr, module);

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

      int hashedIndex = 100000 * eta + phi;
      if (channelMapEE.count(hashedIndex) == 0) {
        channelMapEE[hashedIndex] = channelIteratorEE;
        channelIteratorEE++;
      }
      unsigned int channel = channelMapEE[hashedIndex];
      assert(channel < nCrys);

      setGeomEE(etaG, phiG, iX, iY, iZ, module, tower, ch, apdRefTT, channel, lmr);

      if (_debug == 1)
        edm::LogVerbatim("EcalLaserAnalyzer")
            << "-- debug -- Inside EEDigi - towerID:" << towerID << " channelID:" << channelID << " module:" << module
            << " modules:" << modules.size();

      // APD Pulse
      //===========

      if ((*digiItr).size() > 10)
        edm::LogVerbatim("EcalLaserAnalyzer") << "SAMPLES SIZE > 10!" << (*digiItr).size();

      // Loop on adc samples

      for (unsigned int i = 0; i < (*digiItr).size(); ++i) {
        EcalMGPASample samp_crystal(df.sample(i));
        adc[i] = samp_crystal.adc();
        adcG[i] = samp_crystal.gainId();
        adc[i] *= adcG[i];

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

      APDPulse->setPulse(adc);

      // Quality checks
      //================

      if (adcGain != 1)
        nEvtBadGain[channel]++;
      if (!APDPulse->isTimingQualOK())
        nEvtBadTiming[channel]++;
      nEvtTot[channel]++;

      // Associate PN ampl
      //===================

      int mem0 = Mem->Mem(lmr, 0);
      int mem1 = Mem->Mem(lmr, 1);

      if (allPNAmpl[mem0].size() > MyPn0)
        pn0 = allPNAmpl[mem0][MyPn0];
      else
        pn0 = 0;
      if (allPNAmpl[mem1].size() > MyPn1)
        pn1 = allPNAmpl[mem1][MyPn1];
      else
        pn1 = 0;

      // Fill if Pulse is fine
      //=======================

      if (APDPulse->isPulseOK() && lightside == side) {
        ADCtrees[channel]->Fill();

        Delta01->addEntry(APDPulse->getDelta(0, 1));
        Delta12->addEntry(APDPulse->getDelta(1, 2));

        if (nevtAB[channel] < _nevtmax && _fitab) {
          if (doesABTreeExist)
            shapana->putAllVals(channel, adc, eta, phi);
          else
            shapana->putAllVals(channel, adc, eta, phi, dccID, side, towerID, channelID);
          nevtAB[channel]++;
        }
      }
    }
  }
}

beginJob()

void EcalLaserAnalyzer::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 199 of file EcalLaserAnalyzer.cc.

References _alpha, _beta, _chi2cut, _firstsample, _fitab, _lastsample, _nevtmax, _noise, _nsamples, _presample, adc, ADCfile, ADCFile, ADCtrees, alphafile, alphainitfile, APDfile, channelID, color, dccID, doesABTreeExist, eta, mps_fire::i, laserEvents, Skims_PA_cff::name, nCrys, nevtAB, phi, pn0, pn1, resdir_, resfile, TShapeAnalysis::set_const(), shapana, side, ctpps_dqm_sourceclient-live_cfg::test, and towerID.

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

  // Create temporary files and trees to save adc samples
  //======================================================

  ADCfile = resdir_;
  ADCfile += "/APDSamplesLaser.root";

  APDfile = resdir_;
  APDfile += "/APDPNLaserAllEvents.root";

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

  for (unsigned int i = 0; i < nCrys; i++) {
    std::stringstream name;
    name << "ADCTree" << i + 1;
    ADCtrees[i] = new TTree(name.str().c_str(), name.str().c_str());

    ADCtrees[i]->Branch("ieta", &eta, "eta/I");
    ADCtrees[i]->Branch("iphi", &phi, "phi/I");
    ADCtrees[i]->Branch("side", &side, "side/I");
    ADCtrees[i]->Branch("dccID", &dccID, "dccID/I");
    ADCtrees[i]->Branch("towerID", &towerID, "towerID/I");
    ADCtrees[i]->Branch("channelID", &channelID, "channelID/I");
    ADCtrees[i]->Branch("event", &event, "event/I");
    ADCtrees[i]->Branch("color", &color, "color/I");
    ADCtrees[i]->Branch("adc", &adc, "adc[10]/D");
    ADCtrees[i]->Branch("pn0", &pn0, "pn0/D");
    ADCtrees[i]->Branch("pn1", &pn1, "pn1/D");

    ADCtrees[i]->SetBranchAddress("ieta", &eta);
    ADCtrees[i]->SetBranchAddress("iphi", &phi);
    ADCtrees[i]->SetBranchAddress("side", &side);
    ADCtrees[i]->SetBranchAddress("dccID", &dccID);
    ADCtrees[i]->SetBranchAddress("towerID", &towerID);
    ADCtrees[i]->SetBranchAddress("channelID", &channelID);
    ADCtrees[i]->SetBranchAddress("event", &event);
    ADCtrees[i]->SetBranchAddress("color", &color);
    ADCtrees[i]->SetBranchAddress("adc", adc);
    ADCtrees[i]->SetBranchAddress("pn0", &pn0);
    ADCtrees[i]->SetBranchAddress("pn1", &pn1);

    nevtAB[i] = 0;
  }

  // Define output results filenames and shape analyzer object (alpha,beta)
  //=====================================================================

  // 1) AlphaBeta files

  doesABTreeExist = true;

  std::stringstream nameabinitfile;
  nameabinitfile << resdir_ << "/ABInit.root";
  alphainitfile = nameabinitfile.str();

  std::stringstream nameabfile;
  nameabfile << resdir_ << "/AB.root";
  alphafile = nameabfile.str();

  FILE* test;
  if (_fitab)
    test = fopen(alphainitfile.c_str(), "r");
  else
    test = fopen(alphafile.c_str(), "r");
  if (test == nullptr) {
    doesABTreeExist = false;
    _fitab = true;
  };
  delete test;

  TFile* fAB = nullptr;
  TTree* ABInit = nullptr;
  if (doesABTreeExist) {
    fAB = new TFile(nameabinitfile.str().c_str());
  }
  if (doesABTreeExist && fAB) {
    ABInit = (TTree*)fAB->Get("ABCol0");
  }

  // 2) Shape analyzer

  if (doesABTreeExist && fAB && ABInit && ABInit->GetEntries() != 0) {
    shapana = new TShapeAnalysis(ABInit, _alpha, _beta, 5.5, 1.0);
    doesABTreeExist = true;
  } else {
    shapana = new TShapeAnalysis(_alpha, _beta, 5.5, 1.0);
    doesABTreeExist = false;
    _fitab = true;
  }
  shapana->set_const(_nsamples, _firstsample, _lastsample, _presample, _nevtmax, _noise, _chi2cut);

  if (doesABTreeExist && fAB)
    fAB->Close();

  //  2) APD file

  std::stringstream nameapdfile;
  nameapdfile << resdir_ << "/APDPN_LASER.root";
  resfile = nameapdfile.str();

  // Laser events counter

  laserEvents = 0;
}

endJob()

void EcalLaserAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 709 of file EcalLaserAnalyzer.cc.

References _debug, _ecalPart, _firstsample, _fitab, _lastsample, _niter, _nsamples, _presample, _presamplecut, _qualpercent, _saveallevents, adc, ADCfile, ADCFile, adcNoPed, ADCtrees, TAPD::addEntry(), TPN::addEntry(), simBeamSpotPI::alpha, alphafile, alphainitfile, APD, apdAmpl, apdAmplA, apdAmplB, APDAnal, APDfile, APDFile, APDFirstAnal, APDoAPDA, APDoAPDB, APDoPN, APDoPNA, APDoPNB, APDPulse, apdRefMap, apdTime, APDtrees, cms::cuda::assert(), HLT_2023v12_cff::beta, channelID, channelMapEE, nano_mu_local_reco_cff::chi2, color, colorref, TShapeAnalysis::computeShape(), DMR_cfg::cut, dccID, Delta01, Delta12, doesABTreeExist, PulseFitWithFunction::doFit(), event, eventref, firstChanMod, flag, flagAB, TAPDPulse::getAdcWithoutPedestal(), PulseFitWithFunction::getAmpl(), TAPD::getAPD(), TAPD::getAPDoAPD0(), TAPD::getAPDoAPD1(), TAPD::getAPDoPN(), TAPD::getAPDoPN0(), TAPD::getAPDoPN1(), TMom::getMean(), TPN::getPN(), TPN::getPNoPN(), TPN::getPNoPN0(), TPN::getPNoPN1(), PulseFitWithFunction::getTime(), TAPD::getTime(), edmStreamStallGrapher::getTime(), TShapeAnalysis::getVals(), mps_fire::i, iChannelID, idccID, hcalRecHitTable_cff::ieta, iEta, iModule, PulseFitWithFunction::init(), createfilelist::int, hcalRecHitTable_cff::iphi, iPhi, isGainOK, iside, TAPDPulse::isPulseOK(), IsThereDataADC, isTimingOK, iTowerID, laserEvents, moduleID, Skims_PA_cff::name, nColor, nCrys, nEvtBadGain, nEvtBadTiming, nEvtTot, nMod, nPNPerMod, nRefChan, PN, pn0, pn1, PNAnal, PNFirstAnal, pnID, PNoPN, PNoPNA, PNoPNB, RefAPDtrees, resfile, resFile, respntrees, restrees, TAPD::set2DAPDoAPD0Cut(), TAPD::set2DAPDoAPD1Cut(), TShapeAnalysis::set_nch(), TShapeAnalysis::set_presample(), TAPD::setAPDCut(), TAPD::setAPDoPN0Cut(), TAPD::setAPDoPN1Cut(), TAPD::setAPDoPNCut(), TPN::setPNCut(), TAPDPulse::setPresamples(), TAPDPulse::setPulse(), TAPD::setTimeCut(), shapana, side, Time, towerID, wasGainOK, and wasTimingOK.

Referenced by o2olib.O2ORunMgr::executeJob().

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

  // Adjust channel numbers for EE
  //===============================

  if (_ecalPart == "EE") {
    nCrys = channelMapEE.size();
    shapana->set_nch(nCrys);
  }

  // Set presamples number
  //======================

  double delta01 = Delta01->getMean();
  double delta12 = Delta12->getMean();
  if (delta12 > _presamplecut) {
    _presample = 2;
    if (delta01 > _presamplecut)
      _presample = 1;
  }

  APDPulse->setPresamples(_presample);
  shapana->set_presample(_presample);

  //  Get alpha and beta
  //======================

  if (_fitab) {
    edm::LogVerbatim("EcalLaserAnalyzer") << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
    edm::LogVerbatim("EcalLaserAnalyzer") << "\t+=+     Analyzing data: getting (alpha, beta)     +=+";
    TFile* fAB = nullptr;
    TTree* ABInit = nullptr;
    if (doesABTreeExist) {
      fAB = new TFile(alphainitfile.c_str());
    }
    if (doesABTreeExist && fAB) {
      ABInit = (TTree*)fAB->Get("ABCol0");
    }
    shapana->computeShape(alphafile, ABInit);
    edm::LogVerbatim("EcalLaserAnalyzer") << "\t+=+    .................................... done  +=+";
    edm::LogVerbatim("EcalLaserAnalyzer") << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
  }

  // Don't do anything if there is no events
  //=========================================

  if (laserEvents == 0) {
    ADCFile->Close();
    std::stringstream del;
    del << "rm " << ADCfile;
    system(del.str().c_str());
    edm::LogVerbatim("EcalLaserAnalyzer") << " No Laser Events ";
    return;
  }

  // Set quality flags for gains and timing
  //=========================================

  double BadGainEvtPercentage = 0.0;
  double BadTimingEvtPercentage = 0.0;

  int nChanBadGain = 0;
  int nChanBadTiming = 0;

  for (unsigned int i = 0; i < nCrys; i++) {
    if (nEvtTot[i] != 0) {
      BadGainEvtPercentage = double(nEvtBadGain[i]) / double(nEvtTot[i]);
      BadTimingEvtPercentage = double(nEvtBadTiming[i]) / double(nEvtTot[i]);
    }
    if (BadGainEvtPercentage > _qualpercent) {
      wasGainOK[i] = false;
      nChanBadGain++;
    }
    if (BadTimingEvtPercentage > _qualpercent) {
      wasTimingOK[i] = false;
      nChanBadTiming++;
    }
  }

  double BadGainChanPercentage = double(nChanBadGain) / double(nCrys);
  double BadTimingChanPercentage = double(nChanBadTiming) / double(nCrys);

  if (BadGainChanPercentage > _qualpercent)
    isGainOK = false;
  if (BadTimingChanPercentage > _qualpercent)
    isTimingOK = false;

  // Analyze adc samples to get amplitudes
  //=======================================

  edm::LogVerbatim("EcalLaserAnalyzer") << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
  edm::LogVerbatim("EcalLaserAnalyzer") << "\t+=+     Analyzing laser data: getting APD, PN, APD/PN, PN/PN    +=+";

  if (!isGainOK)
    edm::LogVerbatim("EcalLaserAnalyzer") << "\t+=+ ............................ WARNING! APD GAIN WAS NOT 1    +=+";
  if (!isTimingOK)
    edm::LogVerbatim("EcalLaserAnalyzer") << "\t+=+ ............................ WARNING! TIMING WAS BAD        +=+";

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

  int ieta, iphi, channelID, towerID, flag;

  for (unsigned int i = 0; i < nCrys; i++) {
    std::stringstream name;
    name << "APDTree" << i + 1;

    APDtrees[i] = new TTree(name.str().c_str(), name.str().c_str());

    //List of branches

    APDtrees[i]->Branch("event", &event, "event/I");
    APDtrees[i]->Branch("color", &color, "color/I");
    APDtrees[i]->Branch("iphi", &iphi, "iphi/I");
    APDtrees[i]->Branch("ieta", &ieta, "ieta/I");
    APDtrees[i]->Branch("side", &side, "side/I");
    APDtrees[i]->Branch("dccID", &dccID, "dccID/I");
    APDtrees[i]->Branch("towerID", &towerID, "towerID/I");
    APDtrees[i]->Branch("channelID", &channelID, "channelID/I");
    APDtrees[i]->Branch("apdAmpl", &apdAmpl, "apdAmpl/D");
    APDtrees[i]->Branch("apdTime", &apdTime, "apdTime/D");
    if (_saveallevents)
      APDtrees[i]->Branch("adc", &adc, "adc[10]/D");
    APDtrees[i]->Branch("flag", &flag, "flag/I");
    APDtrees[i]->Branch("flagAB", &flagAB, "flagAB/I");
    APDtrees[i]->Branch("pn0", &pn0, "pn0/D");
    APDtrees[i]->Branch("pn1", &pn1, "pn1/D");

    APDtrees[i]->SetBranchAddress("event", &event);
    APDtrees[i]->SetBranchAddress("color", &color);
    APDtrees[i]->SetBranchAddress("iphi", &iphi);
    APDtrees[i]->SetBranchAddress("ieta", &ieta);
    APDtrees[i]->SetBranchAddress("side", &side);
    APDtrees[i]->SetBranchAddress("dccID", &dccID);
    APDtrees[i]->SetBranchAddress("towerID", &towerID);
    APDtrees[i]->SetBranchAddress("channelID", &channelID);
    APDtrees[i]->SetBranchAddress("apdAmpl", &apdAmpl);
    APDtrees[i]->SetBranchAddress("apdTime", &apdTime);
    if (_saveallevents)
      APDtrees[i]->SetBranchAddress("adc", adc);
    APDtrees[i]->SetBranchAddress("flag", &flag);
    APDtrees[i]->SetBranchAddress("flagAB", &flagAB);
    APDtrees[i]->SetBranchAddress("pn0", &pn0);
    APDtrees[i]->SetBranchAddress("pn1", &pn1);
  }

  for (unsigned int iref = 0; iref < nRefChan; iref++) {
    for (unsigned int imod = 0; imod < nMod; imod++) {
      int jmod = modules[imod];

      std::stringstream nameref;
      nameref << "refAPDTree" << imod << "_" << iref;

      RefAPDtrees[iref][jmod] = new TTree(nameref.str().c_str(), nameref.str().c_str());

      RefAPDtrees[iref][jmod]->Branch("eventref", &eventref, "eventref/I");
      RefAPDtrees[iref][jmod]->Branch("colorref", &colorref, "colorref/I");
      if (iref == 0)
        RefAPDtrees[iref][jmod]->Branch("apdAmplA", &apdAmplA, "apdAmplA/D");
      if (iref == 1)
        RefAPDtrees[iref][jmod]->Branch("apdAmplB", &apdAmplB, "apdAmplB/D");

      RefAPDtrees[iref][jmod]->SetBranchAddress("eventref", &eventref);
      RefAPDtrees[iref][jmod]->SetBranchAddress("colorref", &colorref);
      if (iref == 0)
        RefAPDtrees[iref][jmod]->SetBranchAddress("apdAmplA", &apdAmplA);
      if (iref == 1)
        RefAPDtrees[iref][jmod]->SetBranchAddress("apdAmplB", &apdAmplB);
    }
  }

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

  // Declare PN stuff
  //===================

  for (unsigned int iM = 0; iM < nMod; iM++) {
    unsigned int iMod = modules[iM] - 1;

    for (unsigned int ich = 0; ich < nPNPerMod; ich++) {
      for (unsigned int icol = 0; icol < nCol; icol++) {
        PNFirstAnal[iMod][ich][icol] = new TPN(ich);
        PNAnal[iMod][ich][icol] = new TPN(ich);
      }
    }
  }

  // Declare function for APD ampl fit
  //===================================

  PulseFitWithFunction* pslsfit = new PulseFitWithFunction();
  double chi2;

  for (unsigned int iCry = 0; iCry < nCrys; iCry++) {
    for (unsigned int icol = 0; icol < nCol; icol++) {
      // Declare APD stuff
      //===================

      APDFirstAnal[iCry][icol] = new TAPD();
      IsThereDataADC[iCry][icol] = 1;
      std::stringstream cut;
      cut << "color==" << colors[icol];
      if (ADCtrees[iCry]->GetEntries(cut.str().c_str()) < 10)
        IsThereDataADC[iCry][icol] = 0;
    }

    unsigned int iMod = iModule[iCry] - 1;
    double alpha, beta;

    // Loop on events
    //================

    for (Long64_t jentry = 0; jentry < ADCtrees[iCry]->GetEntriesFast(); jentry++) {
      ADCtrees[iCry]->GetEntry(jentry);

      // Get back color

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

      // Retreive alpha and beta

      std::vector<double> abvals = shapana->getVals(iCry);
      alpha = abvals[0];
      beta = abvals[1];
      flagAB = int(abvals[4]);
      iphi = iPhi[iCry];
      ieta = iEta[iCry];
      towerID = iTowerID[iCry];
      channelID = iChannelID[iCry];

      // Amplitude calculation

      APDPulse->setPulse(adc);
      adcNoPed = APDPulse->getAdcWithoutPedestal();

      apdAmpl = 0;
      apdAmplA = 0;
      apdAmplB = 0;
      apdTime = 0;

      if (APDPulse->isPulseOK()) {
        pslsfit->init(_nsamples, _firstsample, _lastsample, _niter, alpha, beta);
        chi2 = pslsfit->doFit(&adcNoPed[0]);

        if (chi2 < 0. || chi2 == 102 || chi2 == 101) {
          apdAmpl = 0;
          apdTime = 0;
          flag = 0;
        } else {
          apdAmpl = pslsfit->getAmpl();
          apdTime = pslsfit->getTime();
          flag = 1;
        }
      } else {
        apdAmpl = 0;
        apdTime = 0;
        flag = 0;
      }

      if (_debug == 1)
        edm::LogVerbatim("EcalLaserAnalyzer") << "-- debug test -- apdAmpl=" << apdAmpl << ", apdTime=" << apdTime;
      double pnmean;
      if (pn0 < 10 && pn1 > 10) {
        pnmean = pn1;
      } else if (pn1 < 10 && pn0 > 10) {
        pnmean = pn0;
      } else
        pnmean = 0.5 * (pn0 + pn1);

      if (_debug == 1)
        edm::LogVerbatim("EcalLaserAnalyzer") << "-- debug test -- pn0=" << pn0 << ", pn1=" << pn1;

      // Fill PN stuff
      //===============

      if (firstChanMod[iMod] == iCry && IsThereDataADC[iCry][iCol] == 1) {
        for (unsigned int ichan = 0; ichan < nPNPerMod; ichan++) {
          PNFirstAnal[iMod][ichan][iCol]->addEntry(pnmean, pn0, pn1);
        }
      }

      // Fill APD stuff
      //================

      if (APDPulse->isPulseOK()) {
        APDFirstAnal[iCry][iCol]->addEntry(apdAmpl, pnmean, pn0, pn1, apdTime);
        APDtrees[iCry]->Fill();

        // Fill reference trees
        //=====================

        if (apdRefMap[0][iMod + 1] == iCry || apdRefMap[1][iMod + 1] == iCry) {
          apdAmplA = 0.0;
          apdAmplB = 0.0;
          eventref = event;
          colorref = color;

          for (unsigned int ir = 0; ir < nRefChan; ir++) {
            if (apdRefMap[ir][iMod + 1] == iCry) {
              if (ir == 0)
                apdAmplA = apdAmpl;
              else if (ir == 1)
                apdAmplB = apdAmpl;
              RefAPDtrees[ir][iMod + 1]->Fill();
            }
          }
        }
      }
    }
  }

  delete pslsfit;

  ADCFile->Close();

  // Remove temporary file
  //=======================
  std::stringstream del;
  del << "rm " << ADCfile;
  system(del.str().c_str());

  // Create output trees
  //=====================

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

  for (unsigned int iColor = 0; iColor < nCol; iColor++) {
    std::stringstream nametree;
    nametree << "APDCol" << colors[iColor];
    std::stringstream nametree2;
    nametree2 << "PNCol" << colors[iColor];

    restrees[iColor] = new TTree(nametree.str().c_str(), nametree.str().c_str());
    respntrees[iColor] = new TTree(nametree2.str().c_str(), nametree2.str().c_str());

    restrees[iColor]->Branch("iphi", &iphi, "iphi/I");
    restrees[iColor]->Branch("ieta", &ieta, "ieta/I");
    restrees[iColor]->Branch("side", &side, "side/I");
    restrees[iColor]->Branch("dccID", &dccID, "dccID/I");
    restrees[iColor]->Branch("moduleID", &moduleID, "moduleID/I");
    restrees[iColor]->Branch("towerID", &towerID, "towerID/I");
    restrees[iColor]->Branch("channelID", &channelID, "channelID/I");
    restrees[iColor]->Branch("APD", &APD, "APD[6]/D");
    restrees[iColor]->Branch("Time", &Time, "Time[6]/D");
    restrees[iColor]->Branch("APDoPN", &APDoPN, "APDoPN[6]/D");
    restrees[iColor]->Branch("APDoPNA", &APDoPNA, "APDoPNA[6]/D");
    restrees[iColor]->Branch("APDoPNB", &APDoPNB, "APDoPNB[6]/D");
    restrees[iColor]->Branch("APDoAPDA", &APDoAPDA, "APDoAPDA[6]/D");
    restrees[iColor]->Branch("APDoAPDB", &APDoAPDB, "APDoAPDB[6]/D");
    restrees[iColor]->Branch("flag", &flag, "flag/I");

    respntrees[iColor]->Branch("side", &side, "side/I");
    respntrees[iColor]->Branch("moduleID", &moduleID, "moduleID/I");
    respntrees[iColor]->Branch("pnID", &pnID, "pnID/I");
    respntrees[iColor]->Branch("PN", &PN, "PN[6]/D");
    respntrees[iColor]->Branch("PNoPN", &PNoPN, "PNoPN[6]/D");
    respntrees[iColor]->Branch("PNoPNA", &PNoPNA, "PNoPNA[6]/D");
    respntrees[iColor]->Branch("PNoPNB", &PNoPNB, "PNoPNB[6]/D");

    restrees[iColor]->SetBranchAddress("iphi", &iphi);
    restrees[iColor]->SetBranchAddress("ieta", &ieta);
    restrees[iColor]->SetBranchAddress("side", &side);
    restrees[iColor]->SetBranchAddress("dccID", &dccID);
    restrees[iColor]->SetBranchAddress("moduleID", &moduleID);
    restrees[iColor]->SetBranchAddress("towerID", &towerID);
    restrees[iColor]->SetBranchAddress("channelID", &channelID);
    restrees[iColor]->SetBranchAddress("APD", APD);
    restrees[iColor]->SetBranchAddress("Time", Time);
    restrees[iColor]->SetBranchAddress("APDoPN", APDoPN);
    restrees[iColor]->SetBranchAddress("APDoPNA", APDoPNA);
    restrees[iColor]->SetBranchAddress("APDoPNB", APDoPNB);
    restrees[iColor]->SetBranchAddress("APDoAPDA", APDoAPDA);
    restrees[iColor]->SetBranchAddress("APDoAPDB", APDoAPDB);
    restrees[iColor]->SetBranchAddress("flag", &flag);

    respntrees[iColor]->SetBranchAddress("side", &side);
    respntrees[iColor]->SetBranchAddress("moduleID", &moduleID);
    respntrees[iColor]->SetBranchAddress("pnID", &pnID);
    respntrees[iColor]->SetBranchAddress("PN", PN);
    respntrees[iColor]->SetBranchAddress("PNoPN", PNoPN);
    respntrees[iColor]->SetBranchAddress("PNoPNA", PNoPNA);
    respntrees[iColor]->SetBranchAddress("PNoPNB", PNoPNB);
  }

  // Set Cuts for PN stuff
  //=======================

  for (unsigned int iM = 0; iM < nMod; iM++) {
    unsigned int iMod = modules[iM] - 1;

    for (unsigned int ich = 0; ich < nPNPerMod; ich++) {
      for (unsigned int icol = 0; icol < nCol; icol++) {
        PNAnal[iMod][ich][icol]->setPNCut(PNFirstAnal[iMod][ich][icol]->getPN().at(0),
                                          PNFirstAnal[iMod][ich][icol]->getPN().at(1));
      }
    }
  }

  // Build ref trees indexes
  //========================
  for (unsigned int imod = 0; imod < nMod; imod++) {
    int jmod = modules[imod];
    if (RefAPDtrees[0][jmod]->GetEntries() != 0 && RefAPDtrees[1][jmod]->GetEntries() != 0) {
      RefAPDtrees[0][jmod]->BuildIndex("eventref");
      RefAPDtrees[1][jmod]->BuildIndex("eventref");
    }
  }

  // Final loop on crystals
  //=======================

  for (unsigned int iCry = 0; iCry < nCrys; iCry++) {
    unsigned int iMod = iModule[iCry] - 1;

    // Set cuts on APD stuff
    //=======================

    for (unsigned int iCol = 0; iCol < nCol; iCol++) {
      std::vector<double> lowcut;
      std::vector<double> highcut;
      double cutMin;
      double cutMax;

      cutMin = APDFirstAnal[iCry][iCol]->getAPD().at(0) - 2.0 * APDFirstAnal[iCry][iCol]->getAPD().at(1);
      if (cutMin < 0)
        cutMin = 0;
      cutMax = APDFirstAnal[iCry][iCol]->getAPD().at(0) + 2.0 * APDFirstAnal[iCry][iCol]->getAPD().at(1);

      lowcut.push_back(cutMin);
      highcut.push_back(cutMax);

      cutMin = APDFirstAnal[iCry][iCol]->getTime().at(0) - 2.0 * APDFirstAnal[iCry][iCol]->getTime().at(1);
      cutMax = APDFirstAnal[iCry][iCol]->getTime().at(0) + 2.0 * APDFirstAnal[iCry][iCol]->getTime().at(1);
      lowcut.push_back(cutMin);
      highcut.push_back(cutMax);

      APDAnal[iCry][iCol] = new TAPD();
      APDAnal[iCry][iCol]->setAPDCut(APDFirstAnal[iCry][iCol]->getAPD().at(0),
                                     APDFirstAnal[iCry][iCol]->getAPD().at(1));
      APDAnal[iCry][iCol]->setAPDoPNCut(APDFirstAnal[iCry][iCol]->getAPDoPN().at(0),
                                        APDFirstAnal[iCry][iCol]->getAPDoPN().at(1));
      APDAnal[iCry][iCol]->setAPDoPN0Cut(APDFirstAnal[iCry][iCol]->getAPDoPN0().at(0),
                                         APDFirstAnal[iCry][iCol]->getAPDoPN0().at(1));
      APDAnal[iCry][iCol]->setAPDoPN1Cut(APDFirstAnal[iCry][iCol]->getAPDoPN1().at(0),
                                         APDFirstAnal[iCry][iCol]->getAPDoPN1().at(1));
      APDAnal[iCry][iCol]->setTimeCut(APDFirstAnal[iCry][iCol]->getTime().at(0),
                                      APDFirstAnal[iCry][iCol]->getTime().at(1));
      APDAnal[iCry][iCol]->set2DAPDoAPD0Cut(lowcut, highcut);
      APDAnal[iCry][iCol]->set2DAPDoAPD1Cut(lowcut, highcut);
    }

    // Final loop on events
    //=======================

    for (Long64_t jentry = 0; jentry < APDtrees[iCry]->GetEntriesFast(); jentry++) {
      APDtrees[iCry]->GetEntry(jentry);

      double pnmean;
      if (pn0 < 10 && pn1 > 10) {
        pnmean = pn1;
      } else if (pn1 < 10 && pn0 > 10) {
        pnmean = pn0;
      } else
        pnmean = 0.5 * (pn0 + pn1);

      // Get back color
      //================

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

      // Fill PN stuff
      //===============

      if (firstChanMod[iMod] == iCry && IsThereDataADC[iCry][iCol] == 1) {
        for (unsigned int ichan = 0; ichan < nPNPerMod; ichan++) {
          PNAnal[iMod][ichan][iCol]->addEntry(pnmean, pn0, pn1);
        }
      }

      // Get ref amplitudes
      //===================

      if (_debug == 1)
        edm::LogVerbatim("EcalLaserAnalyzer") << "-- debug test -- Last Loop event:" << event << " apdAmpl:" << apdAmpl;
      apdAmplA = 0.0;
      apdAmplB = 0.0;

      for (unsigned int iRef = 0; iRef < nRefChan; iRef++) {
        RefAPDtrees[iRef][iMod + 1]->GetEntryWithIndex(event);
      }

      if (_debug == 1)
        edm::LogVerbatim("EcalLaserAnalyzer")
            << "-- debug test -- Last Loop apdAmplA:" << apdAmplA << " apdAmplB:" << apdAmplB << ", event:" << event
            << ", eventref:" << eventref;

      // Fill APD stuff
      //===============

      APDAnal[iCry][iCol]->addEntry(apdAmpl, pnmean, pn0, pn1, apdTime, apdAmplA, apdAmplB);
    }

    moduleID = iMod + 1;

    if (moduleID >= 20)
      moduleID -= 2;  // Trick to fix endcap specificity

    // Get final results for APD
    //===========================

    for (unsigned int iColor = 0; iColor < nCol; iColor++) {
      std::vector<double> apdvec = APDAnal[iCry][iColor]->getAPD();
      std::vector<double> apdpnvec = APDAnal[iCry][iColor]->getAPDoPN();
      std::vector<double> apdpn0vec = APDAnal[iCry][iColor]->getAPDoPN0();
      std::vector<double> apdpn1vec = APDAnal[iCry][iColor]->getAPDoPN1();
      std::vector<double> timevec = APDAnal[iCry][iColor]->getTime();
      std::vector<double> apdapd0vec = APDAnal[iCry][iColor]->getAPDoAPD0();
      std::vector<double> apdapd1vec = APDAnal[iCry][iColor]->getAPDoAPD1();

      for (unsigned int i = 0; i < apdvec.size(); i++) {
        APD[i] = apdvec.at(i);
        APDoPN[i] = apdpnvec.at(i);
        APDoPNA[i] = apdpn0vec.at(i);
        APDoPNB[i] = apdpn1vec.at(i);
        APDoAPDA[i] = apdapd0vec.at(i);
        APDoAPDB[i] = apdapd1vec.at(i);
        Time[i] = timevec.at(i);
      }

      // Fill APD results trees
      //========================

      iphi = iPhi[iCry];
      ieta = iEta[iCry];
      dccID = idccID[iCry];
      side = iside[iCry];
      towerID = iTowerID[iCry];
      channelID = iChannelID[iCry];

      if (!wasGainOK[iCry] || !wasTimingOK[iCry] || IsThereDataADC[iCry][iColor] == 0) {
        flag = 0;
      } else
        flag = 1;

      restrees[iColor]->Fill();
    }
  }

  // Get final results for PN
  //==========================

  for (unsigned int iM = 0; iM < nMod; iM++) {
    unsigned int iMod = modules[iM] - 1;

    side = iside[firstChanMod[iMod]];

    for (unsigned int ch = 0; ch < nPNPerMod; ch++) {
      pnID = ch;
      moduleID = iMod + 1;

      if (moduleID >= 20)
        moduleID -= 2;  // Trick to fix endcap specificity

      for (unsigned int iColor = 0; iColor < nCol; iColor++) {
        std::vector<double> pnvec = PNAnal[iMod][ch][iColor]->getPN();
        std::vector<double> pnopnvec = PNAnal[iMod][ch][iColor]->getPNoPN();
        std::vector<double> pnopn0vec = PNAnal[iMod][ch][iColor]->getPNoPN0();
        std::vector<double> pnopn1vec = PNAnal[iMod][ch][iColor]->getPNoPN1();

        for (unsigned int i = 0; i < pnvec.size(); i++) {
          PN[i] = pnvec.at(i);
          PNoPN[i] = pnopnvec.at(i);
          PNoPNA[i] = pnopn0vec.at(i);
          PNoPNB[i] = pnopn1vec.at(i);
        }

        // Fill PN results trees
        //========================

        respntrees[iColor]->Fill();
      }
    }
  }

  // Remove temporary files
  //========================
  if (!_saveallevents) {
    APDFile->Close();
    std::stringstream del2;
    del2 << "rm " << APDfile;
    system(del2.str().c_str());

  } else {
    APDFile->cd();
    APDtrees[0]->Write();

    APDFile->Close();
    resFile->cd();
  }

  // Save results
  //===============

  for (unsigned int i = 0; i < nCol; i++) {
    restrees[i]->Write();
    respntrees[i]->Write();
  }

  resFile->Close();

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

setGeomEB()

void EcalLaserAnalyzer::setGeomEB	(	int	etaG,
		int	phiG,
		int	module,
		int	tower,
		int	strip,
		int	xtal,
		int	apdRefTT,
		int	channel,
		int	lmr
	)

Definition at line 1336 of file EcalLaserAnalyzer.cc.

References ME::apdRefChannels(), apdRefMap, cms::cuda::assert(), channelID, submitPVResolutionJobs::count, dccID, eta, firstChanMod, iChannelID, idccID, iEta, iModule, iPhi, isFirstChanModFilled, iside, iTowerID, callgraph::module, nRefChan, phi, MEEBGeom::side(), side, nano_mu_digi_cff::strip, l1tHGCalTowerProducer_cfi::tower, and towerID.

Referenced by analyze().

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

  assert(module >= *min_element(modules.begin(), modules.end()) &&
         module <= *max_element(modules.begin(), modules.end()));

  eta = etaG;
  phi = phiG;
  channelID = 5 * (strip - 1) + xtal - 1;
  towerID = tower;

  std::vector<int> apdRefChan = ME::apdRefChannels(module, lmr);
  for (unsigned int iref = 0; iref < nRefChan; iref++) {
    if (channelID == apdRefChan[iref] && towerID == apdRefTT && apdRefMap[iref].count(module) == 0) {
      apdRefMap[iref][module] = channel;
    }
  }

  if (isFirstChanModFilled[module - 1] == 0) {
    firstChanMod[module - 1] = channel;
    isFirstChanModFilled[module - 1] = 1;
  }

  iEta[channel] = eta;
  iPhi[channel] = phi;
  iModule[channel] = module;
  iTowerID[channel] = towerID;
  iChannelID[channel] = channelID;
  idccID[channel] = dccID;
  iside[channel] = side;
}

setGeomEE()

void EcalLaserAnalyzer::setGeomEE	(	int	etaG,
		int	phiG,
		int	iX,
		int	iY,
		int	iZ,
		int	module,
		int	tower,
		int	ch,
		int	apdRefTT,
		int	channel,
		int	lmr
	)

Definition at line 1369 of file EcalLaserAnalyzer.cc.

References ME::apdRefChannels(), apdRefMap, cms::cuda::assert(), channelID, submitPVResolutionJobs::count, dccID, eta, firstChanMod, iChannelID, idccID, iEta, iModule, iPhi, isFirstChanModFilled, iside, iTowerID, iZ, callgraph::module, nRefChan, phi, MEEEGeom::side(), side, l1tHGCalTowerProducer_cfi::tower, and towerID.

Referenced by analyze().

                                                                                                                   {
  side = MEEEGeom::side(iX, iY, iZ);

  assert(module >= *min_element(modules.begin(), modules.end()) &&
         module <= *max_element(modules.begin(), modules.end()));

  eta = etaG;
  phi = phiG;
  channelID = ch;
  towerID = tower;

  std::vector<int> apdRefChan = ME::apdRefChannels(module, lmr);
  for (unsigned int iref = 0; iref < nRefChan; iref++) {
    if (channelID == apdRefChan[iref] && towerID == apdRefTT && apdRefMap[iref].count(module) == 0) {
      apdRefMap[iref][module] = channel;
    }
  }

  if (isFirstChanModFilled[module - 1] == 0) {
    firstChanMod[module - 1] = channel;
    isFirstChanModFilled[module - 1] = 1;
  }

  iEta[channel] = eta;
  iPhi[channel] = phi;
  iModule[channel] = module;
  iTowerID[channel] = towerID;
  iChannelID[channel] = channelID;
  idccID[channel] = dccID;
  iside[channel] = side;
}

Member Data Documentation

_alpha

const double EcalLaserAnalyzer::_alpha

private

Definition at line 104 of file EcalLaserAnalyzer.h.

Referenced by beginJob().

_beta

const double EcalLaserAnalyzer::_beta

private

Definition at line 105 of file EcalLaserAnalyzer.h.

Referenced by beginJob().

_chi2cut

const double EcalLaserAnalyzer::_chi2cut

private

Definition at line 108 of file EcalLaserAnalyzer.h.

Referenced by beginJob().

_debug

const int EcalLaserAnalyzer::_debug

private

Definition at line 114 of file EcalLaserAnalyzer.h.

Referenced by FrontierCondition_GT_autoExpress_cfi.Tier0Handler::_queryTier0DataSvc(), tier0.Tier0Handler::_queryTier0DataSvc(), analyze(), uploads.uploader::check_response_for_error_key(), endJob(), FrontierCondition_GT_autoExpress_cfi.Tier0Handler::setDebug(), tier0.Tier0Handler::setDebug(), FrontierCondition_GT_autoExpress_cfi.Tier0Handler::unsetDebug(), and tier0.Tier0Handler::unsetDebug().

_docorpn

const bool EcalLaserAnalyzer::_docorpn

private

Definition at line 110 of file EcalLaserAnalyzer.h.

Referenced by analyze().

_ecalPart

const std::string EcalLaserAnalyzer::_ecalPart

private

Definition at line 109 of file EcalLaserAnalyzer.h.

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

_fedid

const int EcalLaserAnalyzer::_fedid

private

Definition at line 111 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and EcalLaserAnalyzer().

_firstsample

const unsigned int EcalLaserAnalyzer::_firstsample

private

Definition at line 88 of file EcalLaserAnalyzer.h.

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

_firstsamplePN

const unsigned int EcalLaserAnalyzer::_firstsamplePN

private

Definition at line 92 of file EcalLaserAnalyzer.h.

Referenced by analyze().

_fitab

bool EcalLaserAnalyzer::_fitab

private

Definition at line 103 of file EcalLaserAnalyzer.h.

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

_lastsample

const unsigned int EcalLaserAnalyzer::_lastsample

private

Definition at line 89 of file EcalLaserAnalyzer.h.

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

_lastsamplePN

const unsigned int EcalLaserAnalyzer::_lastsamplePN

private

Definition at line 93 of file EcalLaserAnalyzer.h.

Referenced by analyze().

_nevtmax

const unsigned int EcalLaserAnalyzer::_nevtmax

private

Definition at line 106 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and beginJob().

_niter

const unsigned int EcalLaserAnalyzer::_niter

private

Definition at line 102 of file EcalLaserAnalyzer.h.

Referenced by endJob().

_noise

const double EcalLaserAnalyzer::_noise

private

Definition at line 107 of file EcalLaserAnalyzer.h.

Referenced by beginJob().

_nsamples

const unsigned int EcalLaserAnalyzer::_nsamples

private

Definition at line 86 of file EcalLaserAnalyzer.h.

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

_nsamplesPN

const unsigned int EcalLaserAnalyzer::_nsamplesPN

private

Definition at line 90 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and EcalLaserAnalyzer().

_presample

unsigned int EcalLaserAnalyzer::_presample

private

Definition at line 87 of file EcalLaserAnalyzer.h.

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

_presamplecut

const double EcalLaserAnalyzer::_presamplecut

private

Definition at line 101 of file EcalLaserAnalyzer.h.

Referenced by endJob().

_presamplePN

const unsigned int EcalLaserAnalyzer::_presamplePN

private

Definition at line 91 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer().

_qualpercent

const double EcalLaserAnalyzer::_qualpercent

private

Definition at line 113 of file EcalLaserAnalyzer.h.

Referenced by endJob().

_ratiomaxcutlow

const double EcalLaserAnalyzer::_ratiomaxcutlow

private

Definition at line 100 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer().

_ratiomincuthigh

const double EcalLaserAnalyzer::_ratiomincuthigh

private

Definition at line 99 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer().

_ratiomincutlow

const double EcalLaserAnalyzer::_ratiomincutlow

private

Definition at line 98 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer().

_saveallevents

const bool EcalLaserAnalyzer::_saveallevents

private

Definition at line 112 of file EcalLaserAnalyzer.h.

Referenced by endJob().

_timingcuthigh

const unsigned int EcalLaserAnalyzer::_timingcuthigh

private

Definition at line 95 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer().

_timingcutlow

const unsigned int EcalLaserAnalyzer::_timingcutlow

private

Definition at line 94 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer().

_timingqualhigh

const unsigned int EcalLaserAnalyzer::_timingqualhigh

private

Definition at line 97 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer().

_timingquallow

const unsigned int EcalLaserAnalyzer::_timingquallow

private

Definition at line 96 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer().

adc

double EcalLaserAnalyzer::adc[10]

private

Definition at line 192 of file EcalLaserAnalyzer.h.

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

ADCfile

std::string EcalLaserAnalyzer::ADCfile

private

Definition at line 131 of file EcalLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

ADCFile

TFile* EcalLaserAnalyzer::ADCFile

private

Definition at line 171 of file EcalLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

adcG

int EcalLaserAnalyzer::adcG[10]

private

Definition at line 193 of file EcalLaserAnalyzer.h.

Referenced by analyze().

adcNoPed

double* EcalLaserAnalyzer::adcNoPed

private

Definition at line 206 of file EcalLaserAnalyzer.h.

Referenced by endJob().

ADCtrees

TTree* EcalLaserAnalyzer::ADCtrees[1700]

private

Definition at line 172 of file EcalLaserAnalyzer.h.

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

alphafile

std::string EcalLaserAnalyzer::alphafile

private

Definition at line 129 of file EcalLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

alphainitfile

std::string EcalLaserAnalyzer::alphainitfile

private

Definition at line 130 of file EcalLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

APD

double EcalLaserAnalyzer::APD[6]

private

Definition at line 224 of file EcalLaserAnalyzer.h.

Referenced by endJob().

apdAmpl

double EcalLaserAnalyzer::apdAmpl

private

Definition at line 197 of file EcalLaserAnalyzer.h.

Referenced by endJob().

apdAmplA

double EcalLaserAnalyzer::apdAmplA

private

Definition at line 198 of file EcalLaserAnalyzer.h.

Referenced by endJob().

apdAmplB

double EcalLaserAnalyzer::apdAmplB

private

Definition at line 199 of file EcalLaserAnalyzer.h.

Referenced by endJob().

APDAnal

TAPD* EcalLaserAnalyzer::APDAnal[1700][nColor]

private

Definition at line 215 of file EcalLaserAnalyzer.h.

Referenced by endJob().

APDfile

std::string EcalLaserAnalyzer::APDfile

private

Definition at line 132 of file EcalLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

APDFile

TFile* EcalLaserAnalyzer::APDFile

private

Definition at line 174 of file EcalLaserAnalyzer.h.

Referenced by endJob().

APDFirstAnal

TAPD* EcalLaserAnalyzer::APDFirstAnal[1700][nColor]

private

Definition at line 214 of file EcalLaserAnalyzer.h.

Referenced by endJob().

APDoAPDA

double EcalLaserAnalyzer::APDoAPDA[6]

private

Definition at line 224 of file EcalLaserAnalyzer.h.

Referenced by endJob().

APDoAPDB

double EcalLaserAnalyzer::APDoAPDB[6]

private

Definition at line 224 of file EcalLaserAnalyzer.h.

Referenced by endJob().

APDoPN

double EcalLaserAnalyzer::APDoPN[6]

private

Definition at line 224 of file EcalLaserAnalyzer.h.

Referenced by endJob().

APDoPNA

double EcalLaserAnalyzer::APDoPNA[6]

private

Definition at line 224 of file EcalLaserAnalyzer.h.

Referenced by endJob().

APDoPNB

double EcalLaserAnalyzer::APDoPNB[6]

private

Definition at line 224 of file EcalLaserAnalyzer.h.

Referenced by endJob().

APDPulse

TAPDPulse* EcalLaserAnalyzer::APDPulse

private

Definition at line 116 of file EcalLaserAnalyzer.h.

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

apdRefMap

std::map<int, unsigned int> EcalLaserAnalyzer::apdRefMap[2]

private

Definition at line 185 of file EcalLaserAnalyzer.h.

Referenced by endJob(), setGeomEB(), and setGeomEE().

apdTime

double EcalLaserAnalyzer::apdTime

private

Definition at line 200 of file EcalLaserAnalyzer.h.

Referenced by endJob().

APDtrees

TTree* EcalLaserAnalyzer::APDtrees[1700]

private

Definition at line 175 of file EcalLaserAnalyzer.h.

Referenced by endJob().

channelID

int EcalLaserAnalyzer::channelID

private

Definition at line 155 of file EcalLaserAnalyzer.h.

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

channelIteratorEE

int EcalLaserAnalyzer::channelIteratorEE

private

Definition at line 222 of file EcalLaserAnalyzer.h.

Referenced by analyze().

channelMapEE

std::map<unsigned int, unsigned int> EcalLaserAnalyzer::channelMapEE

private

Definition at line 183 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and endJob().

color

int EcalLaserAnalyzer::color

private

Definition at line 191 of file EcalLaserAnalyzer.h.

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

colorref

int EcalLaserAnalyzer::colorref

private

Definition at line 204 of file EcalLaserAnalyzer.h.

Referenced by endJob().

colors

std::vector<int> EcalLaserAnalyzer::colors

private

Definition at line 182 of file EcalLaserAnalyzer.h.

dccID

int EcalLaserAnalyzer::dccID

private

Definition at line 157 of file EcalLaserAnalyzer.h.

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

Delta01

TMom* EcalLaserAnalyzer::Delta01

private

Definition at line 119 of file EcalLaserAnalyzer.h.

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

Delta12

TMom* EcalLaserAnalyzer::Delta12

private

Definition at line 120 of file EcalLaserAnalyzer.h.

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

digiCollection_

const std::string EcalLaserAnalyzer::digiCollection_

private

Definition at line 74 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and EcalLaserAnalyzer().

digiPNCollection_

const std::string EcalLaserAnalyzer::digiPNCollection_

private

Definition at line 76 of file EcalLaserAnalyzer.h.

Referenced by analyze().

digiProducer_

const std::string EcalLaserAnalyzer::digiProducer_

private

Definition at line 75 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer().

doesABTreeExist

bool EcalLaserAnalyzer::doesABTreeExist

private

Definition at line 122 of file EcalLaserAnalyzer.h.

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

ebDigiToken_

edm::EDGetTokenT<EBDigiCollection> EcalLaserAnalyzer::ebDigiToken_

private

Definition at line 79 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and EcalLaserAnalyzer().

eeDigiToken_

edm::EDGetTokenT<EEDigiCollection> EcalLaserAnalyzer::eeDigiToken_

private

Definition at line 80 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and EcalLaserAnalyzer().

eta

int EcalLaserAnalyzer::eta

private

Definition at line 189 of file EcalLaserAnalyzer.h.

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

event

int EcalLaserAnalyzer::event

private

Definition at line 190 of file EcalLaserAnalyzer.h.

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

eventHeaderCollection_

const std::string EcalLaserAnalyzer::eventHeaderCollection_

private

Definition at line 72 of file EcalLaserAnalyzer.h.

Referenced by analyze().

eventHeaderProducer_

const std::string EcalLaserAnalyzer::eventHeaderProducer_

private

Definition at line 73 of file EcalLaserAnalyzer.h.

Referenced by analyze().

eventref

int EcalLaserAnalyzer::eventref

private

Definition at line 203 of file EcalLaserAnalyzer.h.

Referenced by endJob().

fedID

int EcalLaserAnalyzer::fedID

private

Definition at line 156 of file EcalLaserAnalyzer.h.

Referenced by analyze().

firstChanMod

unsigned int EcalLaserAnalyzer::firstChanMod[22]

private

Definition at line 233 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), endJob(), setGeomEB(), and setGeomEE().

flag

int EcalLaserAnalyzer::flag

private

Definition at line 221 of file EcalLaserAnalyzer.h.

Referenced by endJob().

flagAB

int EcalLaserAnalyzer::flagAB

private

Definition at line 221 of file EcalLaserAnalyzer.h.

Referenced by endJob().

iChannelID

int EcalLaserAnalyzer::iChannelID[1700]

private

Definition at line 232 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), endJob(), setGeomEB(), and setGeomEE().

idccID

int EcalLaserAnalyzer::idccID[1700]

private

Definition at line 232 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), endJob(), setGeomEB(), and setGeomEE().

iEta

int EcalLaserAnalyzer::iEta[1700]

private

Definition at line 230 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), endJob(), setGeomEB(), and setGeomEE().

iEvent

int EcalLaserAnalyzer::iEvent

private

Definition at line 70 of file EcalLaserAnalyzer.h.

Referenced by analyze().

iModule

unsigned int EcalLaserAnalyzer::iModule[1700]

private

Definition at line 231 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), endJob(), setGeomEB(), and setGeomEE().

iPhi

int EcalLaserAnalyzer::iPhi[1700]

private

Definition at line 230 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), endJob(), setGeomEB(), and setGeomEE().

isFirstChanModFilled

unsigned int EcalLaserAnalyzer::isFirstChanModFilled[22]

private

Definition at line 234 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), setGeomEB(), and setGeomEE().

isGainOK

bool EcalLaserAnalyzer::isGainOK

private

Definition at line 246 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), and endJob().

iside

int EcalLaserAnalyzer::iside[1700]

private

Definition at line 232 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), endJob(), setGeomEB(), and setGeomEE().

IsThereDataADC

int EcalLaserAnalyzer::IsThereDataADC[1700][nColor]

private

Definition at line 217 of file EcalLaserAnalyzer.h.

Referenced by endJob().

isTimingOK

bool EcalLaserAnalyzer::isTimingOK

private

Definition at line 247 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), and endJob().

iTowerID

int EcalLaserAnalyzer::iTowerID[1700]

private

Definition at line 232 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), endJob(), setGeomEB(), and setGeomEE().

iZ

int EcalLaserAnalyzer::iZ

private

Definition at line 160 of file EcalLaserAnalyzer.h.

Referenced by analyze(), EcalLaserAnalyzer(), and setGeomEE().

laserEvents

int EcalLaserAnalyzer::laserEvents

private

Definition at line 163 of file EcalLaserAnalyzer.h.

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

lightside

int EcalLaserAnalyzer::lightside

private

Definition at line 159 of file EcalLaserAnalyzer.h.

Referenced by analyze().

mappingToken_

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

private

Definition at line 82 of file EcalLaserAnalyzer.h.

Referenced by analyze().

Mem

TMem* EcalLaserAnalyzer::Mem

private

Definition at line 118 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and EcalLaserAnalyzer().

moduleID

int EcalLaserAnalyzer::moduleID

private

Definition at line 221 of file EcalLaserAnalyzer.h.

Referenced by endJob().

modules

std::vector<int> EcalLaserAnalyzer::modules

private

Definition at line 184 of file EcalLaserAnalyzer.h.

Referenced by SequenceVisitors.NodeVisitor::enter().

nCrys

unsigned int EcalLaserAnalyzer::nCrys

private

Definition at line 140 of file EcalLaserAnalyzer.h.

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

nevtAB

unsigned int EcalLaserAnalyzer::nevtAB[1700]

private

Definition at line 136 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and beginJob().

nEvtBadGain

int EcalLaserAnalyzer::nEvtBadGain[1700]

private

Definition at line 238 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and endJob().

nEvtBadTiming

int EcalLaserAnalyzer::nEvtBadTiming[1700]

private

Definition at line 239 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and endJob().

nEvtTot

int EcalLaserAnalyzer::nEvtTot[1700]

private

Definition at line 240 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and endJob().

nMod

unsigned int EcalLaserAnalyzer::nMod

private

Definition at line 144 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), and endJob().

nPNPerMod

unsigned int EcalLaserAnalyzer::nPNPerMod

private

Definition at line 141 of file EcalLaserAnalyzer.h.

Referenced by endJob().

nRefChan

unsigned int EcalLaserAnalyzer::nRefChan

private

Definition at line 142 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), endJob(), setGeomEB(), and setGeomEE().

nRefTrees

unsigned int EcalLaserAnalyzer::nRefTrees

private

Definition at line 143 of file EcalLaserAnalyzer.h.

nSides

unsigned int EcalLaserAnalyzer::nSides

private

Definition at line 145 of file EcalLaserAnalyzer.h.

phi

int EcalLaserAnalyzer::phi

private

Definition at line 189 of file EcalLaserAnalyzer.h.

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

pn

double EcalLaserAnalyzer::pn[50]

private

Definition at line 195 of file EcalLaserAnalyzer.h.

Referenced by analyze().

PN

double EcalLaserAnalyzer::PN[6]

private

Definition at line 224 of file EcalLaserAnalyzer.h.

Referenced by endJob().

pn0

double EcalLaserAnalyzer::pn0

private

Definition at line 194 of file EcalLaserAnalyzer.h.

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

pn1

double EcalLaserAnalyzer::pn1

private

Definition at line 194 of file EcalLaserAnalyzer.h.

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

pnAmpl

double EcalLaserAnalyzer::pnAmpl

private

Definition at line 201 of file EcalLaserAnalyzer.h.

Referenced by analyze().

PNAnal

TPN* EcalLaserAnalyzer::PNAnal[22][2][nColor]

private

Definition at line 211 of file EcalLaserAnalyzer.h.

Referenced by endJob().

pncorfile_

const std::string EcalLaserAnalyzer::pncorfile_

private

Definition at line 125 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer().

pnCorrector

TPNCor* EcalLaserAnalyzer::pnCorrector

private

Definition at line 167 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and EcalLaserAnalyzer().

pnDiodeDigiToken_

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

private

Definition at line 81 of file EcalLaserAnalyzer.h.

Referenced by analyze().

PNFirstAnal

TPN* EcalLaserAnalyzer::PNFirstAnal[22][2][nColor]

private

Definition at line 210 of file EcalLaserAnalyzer.h.

Referenced by endJob().

pnG

int EcalLaserAnalyzer::pnG[50]

private

Definition at line 196 of file EcalLaserAnalyzer.h.

Referenced by analyze().

pnID

int EcalLaserAnalyzer::pnID

private

Definition at line 221 of file EcalLaserAnalyzer.h.

Referenced by endJob().

pnNoPed

double* EcalLaserAnalyzer::pnNoPed

private

Definition at line 207 of file EcalLaserAnalyzer.h.

Referenced by analyze().

PNoPN

double EcalLaserAnalyzer::PNoPN[6]

private

Definition at line 224 of file EcalLaserAnalyzer.h.

Referenced by endJob().

PNoPNA

double EcalLaserAnalyzer::PNoPNA[6]

private

Definition at line 224 of file EcalLaserAnalyzer.h.

Referenced by endJob().

PNoPNB

double EcalLaserAnalyzer::PNoPNB[6]

private

Definition at line 224 of file EcalLaserAnalyzer.h.

Referenced by endJob().

PNPulse

TPNPulse* EcalLaserAnalyzer::PNPulse

private

Definition at line 117 of file EcalLaserAnalyzer.h.

Referenced by analyze(), and EcalLaserAnalyzer().

rawDataToken_

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

private

Definition at line 78 of file EcalLaserAnalyzer.h.

Referenced by analyze().

RefAPDtrees

TTree* EcalLaserAnalyzer::RefAPDtrees[2][22]

private

Definition at line 176 of file EcalLaserAnalyzer.h.

Referenced by endJob().

resdir_

const std::string EcalLaserAnalyzer::resdir_

private

Definition at line 124 of file EcalLaserAnalyzer.h.

Referenced by beginJob().

resfile

std::string EcalLaserAnalyzer::resfile

private

Definition at line 133 of file EcalLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

resFile

TFile* EcalLaserAnalyzer::resFile

private

Definition at line 178 of file EcalLaserAnalyzer.h.

Referenced by endJob().

respntrees

TTree* EcalLaserAnalyzer::respntrees[nColor]

private

Definition at line 180 of file EcalLaserAnalyzer.h.

Referenced by endJob().

restrees

TTree* EcalLaserAnalyzer::restrees[nColor]

private

Definition at line 179 of file EcalLaserAnalyzer.h.

Referenced by endJob().

runNum

int EcalLaserAnalyzer::runNum

private

Definition at line 150 of file EcalLaserAnalyzer.h.

Referenced by analyze().

runType

int EcalLaserAnalyzer::runType

private

Definition at line 149 of file EcalLaserAnalyzer.h.

Referenced by analyze().

shapana

TShapeAnalysis* EcalLaserAnalyzer::shapana

private

Definition at line 135 of file EcalLaserAnalyzer.h.

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

side

int EcalLaserAnalyzer::side

private

Definition at line 158 of file EcalLaserAnalyzer.h.

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

Time

double EcalLaserAnalyzer::Time[6]

private

Definition at line 224 of file EcalLaserAnalyzer.h.

Referenced by endJob().

towerID

int EcalLaserAnalyzer::towerID

private

Definition at line 154 of file EcalLaserAnalyzer.h.

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

wasABCalcOK

bool EcalLaserAnalyzer::wasABCalcOK[1700]

private

Definition at line 244 of file EcalLaserAnalyzer.h.

wasGainOK

bool EcalLaserAnalyzer::wasGainOK[1700]

private

Definition at line 242 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), and endJob().

wasTimingOK

bool EcalLaserAnalyzer::wasTimingOK[1700]

private

Definition at line 243 of file EcalLaserAnalyzer.h.

Referenced by EcalLaserAnalyzer(), and endJob().