#include <EcalLaserAnalyzer2.h>

Inheritance diagram for EcalLaserAnalyzer2:

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

	EcalLaserAnalyzer2 (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)

	~EcalLaserAnalyzer2 () 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 Member Functions
bool	getShapes ()

Private Attributes
const int	_debug

const bool	_docorpn

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 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 bool	_saveshapes

const unsigned int	_timingcuthigh

const unsigned int	_timingcutlow

const unsigned int	_timingqualhigh

const unsigned int	_timingquallow

double	adc [10]

int	adcC [1700]

std::string	ADCfile

TFile *	ADCFile

int	adcG [10]

double	adcMean [1700][10]

double *	adcNoPed

TTree *	ADCtrees [1700]

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< int, int >	channelMapEE

int	color

int	colorref

std::vector< int >	colors

int	dccID

std::vector< int >	dccMEM

TMom *	Delta01

TMom *	Delta12

const std::string	digiCollection_

const std::string	digiPNCollection_

const std::string	digiProducer_

edm::EDGetTokenT< EBDigiCollection >	ebDigiToken_

edm::EDGetTokenT< EEDigiCollection >	eeDigiToken_

const std::string	elecfile_

int	eta

int	event

const std::string	eventHeaderCollection_

const std::string	eventHeaderProducer_

int	eventref

int	fedID

unsigned int	firstChanMod [21]

int	flag

int	iChannelID [1700]

int	idccID [1700]

int	iEta [1700]

int	iEvent

unsigned int	iModule [1700]

int	iPhi [1700]

unsigned int	isFirstChanModFilled [21]

bool	isGainOK

int	iside [1700]

bool	IsMatacqOK

bool	isMatacqOK

bool	isSPRFine

int	IsThereDataADC [1700][nColor]

bool	isTimingOK

int	iTowerID [1700]

int	iZ

int	laserEvents

int	lightside

const edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	mappingToken_

std::string	matfile

TMem *	Mem

int	moduleID

std::vector< int >	modules

unsigned int	nCrys

int	nEvtBadGain [1700]

int	nEvtBadTiming [1700]

int	nEvtTot [1700]

unsigned int	nMod

unsigned int	nPNPerMod

unsigned int	nRefChan

unsigned int	nRefTrees

unsigned int	nSamplesShapes

unsigned int	nSides

int	phi

double	pn [50]

double	PN [6]

double	pn0

double	pn1

double	pnAmpl

TPN *	PNAnal [9][2][nColor]

const std::string	pncorfile_

TPNCor *	pnCorrector

const edm::EDGetTokenT< EcalPnDiodeDigiCollection >	pnDiodeDigiToken_

TPN *	PNFirstAnal [9][2][nColor]

int	pnG [50]

int	pnID

double *	pnNoPed

double	PNoPN [6]

double	PNoPNA [6]

double	PNoPNB [6]

TPNPulse *	PNPulse

TProfile *	PulseShape

const edm::EDGetTokenT< EcalRawDataCollection >	rawDataToken_

TTree *	RefAPDtrees [2][21]

const std::string	resdir_

std::string	resfile

TFile *	resFile

TTree *	respntrees [nColor]

TTree *	restrees [nColor]

int	runNum

int	runType

double	ShapeCor

double	shapeCorrection

std::string	shapefile

TFile *	ShapeFile

double	shapes [250]

std::vector< double >	shapesVec

int	side

double	Time [6]

int	towerID

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 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 ()

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 EcalLaserAnalyzer2.h.

Member Enumeration Documentation

◆ VarCol

enum EcalLaserAnalyzer2::VarCol

Enumerator
iBlue
iRed
nColor

Definition at line 67 of file EcalLaserAnalyzer2.h.

67 { iBlue, iRed, nColor };

EcalLaserAnalyzer2::iRed

Definition: EcalLaserAnalyzer2.h:67

EcalLaserAnalyzer2::nColor

Definition: EcalLaserAnalyzer2.h:67

EcalLaserAnalyzer2::iBlue

Definition: EcalLaserAnalyzer2.h:67

Constructor & Destructor Documentation

◆ EcalLaserAnalyzer2()

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

explicit

Definition at line 48 of file EcalLaserAnalyzer2.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, HLT_2022v12_cff::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_, digiPNCollection_))),
       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", 5)),
       _noise(iConfig.getUntrackedParameter<double>("noise", 2.0)),
       _ecalPart(iConfig.getUntrackedParameter<std::string>("ecalPart", "EB")),
       _saveshapes(iConfig.getUntrackedParameter<bool>("saveShapes", true)),
       _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")),
       elecfile_(iConfig.getUntrackedParameter<std::string>("elecFile")),
       pncorfile_(iConfig.getUntrackedParameter<std::string>("pnCorFile")),
       nCrys(NCRYSEB),
       nPNPerMod(NPNPERMOD),
       nMod(NMODEB),
       nSides(NSIDES),
       nSamplesShapes(NSAMPSHAPES),
       IsMatacqOK(false),
       runType(-1),
       runNum(0),
       towerID(-1),
       channelID(-1),
       fedID(-1),
       dccID(-1),
       side(2),
       lightside(2),
       iZ(1),
       phi(-1),
       eta(-1),
       event(0),
       color(0),
       pn0(0),
       pn1(0),
       apdAmpl(0),
       apdTime(0),
       pnAmpl(0),
       pnID(-1),
       moduleID(-1),
       flag(0),
       channelIteratorEE(0),
       ShapeCor(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();
 }

◆ ~EcalLaserAnalyzer2()

EcalLaserAnalyzer2::~EcalLaserAnalyzer2 ( )

override

Definition at line 193 of file EcalLaserAnalyzer2.cc.

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

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 275 of file EcalLaserAnalyzer2.cc.

References _debug, _docorpn, _ecalPart, _fedid, _firstsamplePN, _lastsamplePN, _nsamplesPN, adc, adcG, ADCtrees, TMom::addEntry(), APDPulse, MEEBGeom::apdRefTower(), MEEEGeom::apdRefTower(), cms::cuda::assert(), edm::SortedCollection< T, SORT >::begin(), c, EcalElectronicsId::channelId(), channelID, channelIteratorEE, channelMapEE, color, alignCSCRings::corr, dccID, MEEEGeom::dee(), Delta01, Delta12, hgcalPerformanceValidation::df, digiCollection_, digiPNCollection_, TPNFit::doFit(), MillePedeFileConverter_cfg::e, CollectionTags_cfi::EBDigi, ebDigiToken_, CollectionTags_cfi::EEDigi, eeDigiToken_, MEEBGeom::electronic_channel(), edm::SortedCollection< T, SORT >::end(), eta, eventHeaderCollection_, 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(), IsMatacqOK, 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, callgraph::module, nCrys, nEvtBadGain, nEvtBadTiming, nEvtTot, phi, MEEBGeom::pn(), MEEEGeom::pn(), pn, pn0, pn1, pnAmpl, pnCorrector, pnDiodeDigiToken_, pnG, pnNoPed, PNPulse, edm::Handle< T >::product(), rawDataToken_, runNum, runType, setGeomEB(), setGeomEE(), TPNPulse::setPulse(), TAPDPulse::setPulse(), side, findQualityFiles::size, digitizers_cfi::strip, EcalElectronicsId::stripId(), hgcalTowerProducer_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();
   } 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 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();
     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;
 
     vector<int>::iterator iter = find(colors.begin(), colors.end(), color);
     if (iter == colors.end()) {
       colors.push_back(color);
       edm::LogVerbatim("EcalLaserAnalyzer2") << " new color found " << color << " " << colors.size();
     }
   }
 
   // Check Matacq shape exists
 
   if (!IsMatacqOK)
     return;
 
   // 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;
 
   map<int, vector<double> > allPNAmpl;
   map<int, 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("EcalLaserAnalyzer2")
           << "-- debug test -- 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("EcalLaserAnalyzer2") << "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("EcalLaserAnalyzer2")
           << "-- 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) {  // Loop on crystals
 
       // 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("EcalLaserAnalyzer2")
             << "-- 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));
       }
     }
 
   } 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("EcalLaserAnalyzer2")
             << "-- debug -- Inside EEDigi - towerID:" << towerID << " channelID:" << channelID << " module:" << module
             << " modules:" << modules.size();
 
       // APD Pulse
       //===========
 
       if ((*digiItr).size() > 10)
         edm::LogVerbatim("EcalLaserAnalyzer2") << "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));
       }
     }
   }
 }

◆ beginJob()

void EcalLaserAnalyzer2::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 201 of file EcalLaserAnalyzer2.cc.

References adc, ADCfile, ADCFile, ADCtrees, APDfile, channelID, color, dccID, eta, getShapes(), mps_fire::i, IsMatacqOK, laserEvents, matfile, Skims_PA_cff::name, nCrys, phi, pn0, pn1, resdir_, resfile, shapefile, side, 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++) {
     stringstream name;
     name << "ADCTree" << i + 1;
     ADCtrees[i] = new TTree(name.str().c_str(), name.str().c_str());
 
     ADCtrees[i]->Branch("iphi", &phi, "phi/I");
     ADCtrees[i]->Branch("ieta", &eta, "eta/I");
     ADCtrees[i]->Branch("towerID", &towerID, "towerID/I");
     ADCtrees[i]->Branch("channelID", &channelID, "channelID/I");
     ADCtrees[i]->Branch("dccID", &dccID, "dccID/I");
     ADCtrees[i]->Branch("side", &side, "side/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("towerID", &towerID);
     ADCtrees[i]->SetBranchAddress("channelID", &channelID);
     ADCtrees[i]->SetBranchAddress("dccID", &dccID);
     ADCtrees[i]->SetBranchAddress("side", &side);
     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);
   }
 
   // Define output results filenames
   //==================================
   stringstream namefile1;
   namefile1 << resdir_ << "/SHAPE_LASER.root";
   shapefile = namefile1.str();
 
   stringstream namefile2;
   namefile2 << resdir_ << "/APDPN_LASER.root";
   resfile = namefile2.str();
 
   stringstream namefile3;
   namefile3 << resdir_ << "/MATACQ.root";
 
   matfile = namefile3.str();
 
   // Get Pulse Shapes
   //==================
 
   IsMatacqOK = getShapes();
   if (!IsMatacqOK) {
     edm::LogError("noshape") << " ERROR! No matacq shape available: analysis aborted !";
     return;
   }
 
   // Laser events counter
 
   laserEvents = 0;
 }

◆ endJob()

void EcalLaserAnalyzer2::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 668 of file EcalLaserAnalyzer2.cc.

References _debug, _ecalPart, _firstsample, _lastsample, _niter, _noise, _nsamples, _presample, _presamplecut, _qualpercent, _saveallevents, adc, ADCfile, ADCFile, adcNoPed, ADCtrees, TAPD::addEntry(), TPN::addEntry(), APD, apdAmpl, apdAmplA, apdAmplB, APDAnal, APDfile, APDFile, APDFirstAnal, APDoAPDA, APDoAPDB, APDoPN, APDoPNA, APDoPNB, APDPulse, apdRefMap, apdTime, APDtrees, cms::cuda::assert(), channelID, channelMapEE, color, colorref, PA_MinBiasSkim_cff::cut, dccID, Delta01, Delta12, PulseFitWithShape::doFit(), eta, event, eventref, firstChanMod, flag, TAPDPulse::getAdcWithoutPedestal(), PulseFitWithShape::getAmpl(), TAPD::getAPD(), TAPD::getAPDoAPD0(), TAPD::getAPDoAPD1(), TAPD::getAPDoPN(), TAPD::getAPDoPN0(), TAPD::getAPDoPN1(), TMom::getMean(), TPN::getPN(), TPN::getPNoPN(), TPN::getPNoPN0(), TPN::getPNoPN1(), PulseFitWithShape::getTime(), TAPD::getTime(), edmStreamStallGrapher::getTime(), mps_fire::i, iChannelID, idccID, LEDCalibrationChannels::ieta, iEta, iModule, PulseFitWithShape::init(), LEDCalibrationChannels::iphi, iPhi, isGainOK, iside, IsMatacqOK, TAPDPulse::isPulseOK(), isSPRFine, IsThereDataADC, isTimingOK, iTowerID, laserEvents, moduleID, Skims_PA_cff::name, nColor, nCrys, nEvtBadGain, nEvtBadTiming, nEvtTot, nMod, nPNPerMod, nRefChan, nSamplesShapes, phi, PN, pn0, pn1, PNAnal, PNFirstAnal, pnID, PNoPN, PNoPNA, PNoPNB, RefAPDtrees, resfile, resFile, respntrees, restrees, TAPD::set2DAPDoAPD0Cut(), TAPD::set2DAPDoAPD1Cut(), TAPD::setAPDCut(), TAPD::setAPDoPN0Cut(), TAPD::setAPDoPN1Cut(), TAPD::setAPDoPNCut(), TPN::setPNCut(), TAPDPulse::setPresamples(), TAPDPulse::setPulse(), TAPD::setTimeCut(), ShapeCor, shapeCorrection, shapesVec, side, Time, towerID, wasGainOK, and wasTimingOK.

Referenced by o2olib.O2ORunMgr::executeJob().

                                 {
   //========================================================================
 
   if (!IsMatacqOK) {
     edm::LogVerbatim("EcalLaserAnalyzer2") << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
     edm::LogVerbatim("EcalLaserAnalyzer2") << "\t+=+     WARNING! NO MATACQ        +=+";
     edm::LogVerbatim("EcalLaserAnalyzer2") << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
     return;
   }
 
   // Adjust channel numbers for EE
   //===============================
 
   if (_ecalPart == "EE") {
     nCrys = channelMapEE.size();
   }
 
   // Set presamples number
   //======================
 
   double delta01 = Delta01->getMean();
   double delta12 = Delta12->getMean();
   if (delta12 > _presamplecut) {
     _presample = 2;
     if (delta01 > _presamplecut)
       _presample = 1;
   }
 
   APDPulse->setPresamples(_presample);
 
   // Don't do anything if there is no events
   //=========================================
 
   if (laserEvents == 0) {
     ADCFile->Close();
     stringstream del;
     del << "rm " << ADCfile;
     system(del.str().c_str());
     edm::LogVerbatim("EcalLaserAnalyzer2") << " 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("EcalLaserAnalyzer2") << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
   edm::LogVerbatim("EcalLaserAnalyzer2") << "\t+=+     Analyzing laser data: getting APD, PN, APD/PN, PN/PN    +=+";
 
   if (!isGainOK)
     edm::LogVerbatim("EcalLaserAnalyzer2") << "\t+=+ ............................ WARNING! APD GAIN WAS NOT 1    +=+";
   if (!isTimingOK)
     edm::LogVerbatim("EcalLaserAnalyzer2") << "\t+=+ ............................ WARNING! TIMING WAS BAD        +=+";
 
   APDFile = new TFile(APDfile.c_str(), "RECREATE");
 
   int ieta, iphi;
 
   int flagfit;
   for (unsigned int i = 0; i < nCrys; i++) {
     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("flagfit", &flagfit, "flagfit/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("flagfit", &flagfit);
     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];
 
       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
   //===================================
 
   PulseFitWithShape* psfit = new PulseFitWithShape();
 
   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;
       stringstream cut;
       cut << "color==" << colors.at(icol);
       if (ADCtrees[iCry]->GetEntries(cut.str().c_str()) < 10)
         IsThereDataADC[iCry][icol] = 0;
     }
 
     unsigned int iMod = iModule[iCry] - 1;
     assert(iMod <= nMod);
 
     if (isSPRFine)
       psfit->init(_nsamples, _firstsample, _lastsample, _niter, nSamplesShapes, shapesVec, _noise);
 
     // Loop on events
     //================
 
     Long64_t nbytes = 0, nb = 0;
     for (Long64_t jentry = 0; jentry < ADCtrees[iCry]->GetEntriesFast(); jentry++) {  // Loop on events
       nb = ADCtrees[iCry]->GetEntry(jentry);
       nbytes += nb;
 
       flagfit = 1;
       apdAmpl = 0.0;
       apdTime = 0.0;
       ieta = eta;
       iphi = phi;
 
       // Get back color
 
       unsigned int iCol = 0;
       for (unsigned int i = 0; i < nCol; i++) {
         if (color == colors[i]) {
           iCol = i;
           i = colors.size();
         }
       }
 
       // Amplitude calculation
 
       APDPulse->setPulse(adc);
       adcNoPed = APDPulse->getAdcWithoutPedestal();
 
       if (isSPRFine && APDPulse->isPulseOK()) {
         psfit->doFit(&adcNoPed[0]);
         apdAmpl = psfit->getAmpl();
         apdTime = psfit->getTime();
 
       } else {
         apdAmpl = 0;
         apdTime = 0;
         flagfit = 0;
       }
 
       if (_debug >= 1)
         edm::LogVerbatim("EcalLaserAnalyzer2")
             << "-- debug test -- endJob -- 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("EcalLaserAnalyzer2") << "-- debug test -- endJob -- pnMean:" << pnmean;
 
       // 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 (apdAmpl != 0.0)
         APDFirstAnal[iCry][iCol]->addEntry(apdAmpl, pnmean, pn0, pn1, apdTime);
       if (_debug >= 1)
         edm::LogVerbatim("EcalLaserAnalyzer2") << "-- debug test -- endJob -- filling APDTree";
 
       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 (_debug >= 1)
             edm::LogVerbatim("EcalLaserAnalyzer2") << "-- debug test -- ir:" << ir << " tt:" << towerID
                                                    << " refmap:" << apdRefMap[ir][iMod + 1] << " iCry:" << iCry;
 
           if (apdRefMap[ir][iMod + 1] == iCry) {
             if (_debug >= 1)
               edm::LogVerbatim("EcalLaserAnalyzer2") << "-- debug test -- cut passed ";
             if (ir == 0)
               apdAmplA = apdAmpl;
             else if (ir == 1)
               apdAmplB = apdAmpl;
             if (_debug >= 1)
               edm::LogVerbatim("EcalLaserAnalyzer2") << "-- debug test -- apdAmplA=" << apdAmplA;
             if (_debug >= 1)
               edm::LogVerbatim("EcalLaserAnalyzer2") << "-- debug test -- apdAmplB=" << apdAmplB;
             if (_debug >= 1)
               edm::LogVerbatim("EcalLaserAnalyzer2") << "-- debug test -- color=" << color << ", event:" << event
                                                      << ", ir:" << ir << " tt-1:" << towerID - 1;
 
             RefAPDtrees[ir][iMod + 1]->Fill();
 
             if (_debug >= 1)
               edm::LogVerbatim("EcalLaserAnalyzer2") << "-- debug test -- tree filled" << event;
           }
         }
       }
     }
   }
 
   delete psfit;
 
   ADCFile->Close();
 
   if (_debug == 1)
     edm::LogVerbatim("EcalLaserAnalyzer2") << "-- debug test -- endJob -- after apdAmpl Loop";
 
   // Remove temporary file
   //=======================
   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++) {
     stringstream nametree;
     nametree << "APDCol" << colors.at(iColor);
     stringstream nametree2;
     nametree2 << "PNCol" << colors.at(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("ShapeCor", &ShapeCor, "ShapeCor/D");
     restrees[iColor]->Branch("flag", &flag, "flag/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("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("ShapeCor", &ShapeCor);
     restrees[iColor]->SetBranchAddress("flag", &flag);
 
     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
     //=======================
 
     Long64_t nbytes = 0, nb = 0;
     for (Long64_t jentry = 0; jentry < APDtrees[iCry]->GetEntriesFast(); jentry++) {
       nb = APDtrees[iCry]->GetEntry(jentry);
       nbytes += nb;
 
       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("EcalLaserAnalyzer2") << "-- debug test -- LastLoop 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("EcalLaserAnalyzer2")
             << "-- debug test -- LastLoop 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);
         ShapeCor = shapeCorrection;
       }
 
       // Fill APD results trees
       //========================
 
       iphi = iPhi[iCry];
       ieta = iEta[iCry];
       dccID = idccID[iCry];
       towerID = iTowerID[iCry];
       side = iside[iCry];
       channelID = iChannelID[iCry];
 
       if (!wasGainOK[iCry] || !wasTimingOK[iCry] || IsThereDataADC[iCry][iColor] == 0)
         flag = 1;
       else
         flag = 0;
 
       if (_debug >= 1)
         edm::LogVerbatim("EcalLaserAnalyzer2") << "-- debug test -- endJob -- APD[0]" << APD[0] << " APDoPN[0] "
                                                << APDoPN[0] << " APDoAPDA[0] " << APDoAPDA[0] << " flag " << flag;
       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);
         }
 
         if (_debug >= 1)
           edm::LogVerbatim("EcalLaserAnalyzer2")
               << "-- debug test -- endJob -- filling pn results'tree: PN[0]:" << PN[0] << " iModule:" << iMod
               << " iColor:" << iColor << " ch:" << ch;
 
         // Fill PN results trees
         //========================
 
         respntrees[iColor]->Fill();
       }
     }
   }
 
   // Remove temporary files
   //========================
   if (!_saveallevents) {
     APDFile->Close();
     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("EcalLaserAnalyzer2") << "\t+=+    .................................................. done  +=+";
   edm::LogVerbatim("EcalLaserAnalyzer2") << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
 }

◆ getShapes()

bool EcalLaserAnalyzer2::getShapes ( )

private

Definition at line 1300 of file EcalLaserAnalyzer2.cc.

References _saveshapes, cms::cuda::assert(), elecfile_, Exception, cuy::ii, createfilelist::int, IsMatacqOK, isSPRFine, findQualityFiles::jj, matfile, Skims_PA_cff::name, seedmultiplicitymonitor_newtracking_cfi::nBins, nSamplesShapes, PulseShape, L1EGammaClusterEmuProducer_cfi::scale, shapeCorrection, shapefile, ShapeFile, shapes, shapesVec, ctpps_dqm_sourceclient-live_cfg::test, x, and y.

Referenced by beginJob().

                                    {
   //========================================================================
 
   // Get Pulse From Matacq Analysis:
   //================================
 
   bool IsMatacqOK = false;
 
   int doesMatFileExist = 0;
   int doesMatShapeExist = 0;
   FILE* test2;
   TProfile* laserShape = nullptr;
   test2 = fopen(matfile.c_str(), "r");
   if (test2)
     doesMatFileExist = 1;
 
   TFile* MatShapeFile;
   if (doesMatFileExist == 1) {
     MatShapeFile = new TFile(matfile.c_str());
     laserShape = (TProfile*)MatShapeFile->Get("shapeLaser");
     if (laserShape) {
       doesMatShapeExist = 1;
       double y = laserShape->Integral("w");
       if (y != 0)
         laserShape->Scale(1.0 / y);
     }
   } else {
     edm::LogError("file_not_found") << " ERROR! Matacq shape file not found !";
   }
   if (doesMatShapeExist)
     IsMatacqOK = true;
 
   // Get SPR from the average elec shape in SM6:
   //============================================
 
   int doesElecFileExist = 0;
   FILE* test;
   test = fopen(elecfile_.c_str(), "r");
   if (test)
     doesElecFileExist = 1;
 
   TFile* ElecShapesFile;
   TH1D* elecShape = nullptr;
 
   if (doesElecFileExist == 1) {
     ElecShapesFile = new TFile(elecfile_.c_str());
     stringstream name;
     name << "MeanElecShape";
     elecShape = (TH1D*)ElecShapesFile->Get(name.str().c_str());
     if (elecShape && doesMatShapeExist == 1) {
       double x = elecShape->GetMaximum();
       if (x != 0)
         elecShape->Scale(1.0 / x);
       isSPRFine = true;
     } else {
       isSPRFine = false;
     }
 
   } else {
     throw cms::Exception("file_not_found") << " ERROR! Elec shape file not found !";
   }
 
   if (IsMatacqOK) {
     ShapeFile = new TFile(shapefile.c_str(), "RECREATE");
 
     unsigned int nBins = int(laserShape->GetEntries());
     assert(nSamplesShapes == nBins);
     double elec_jj, laser_iiMinusjj;
     double sum_jj;
 
     if (isSPRFine == true) {
       unsigned int nBins2 = int(elecShape->GetNbinsX());
 
       if (nBins2 < nBins) {
         edm::LogError("cfg_error")
             << "EcalLaserAnalyzer2::getShapes: wrong configuration of the shapes' number of bins";
         isSPRFine = false;
       }
       assert(nSamplesShapes == nBins2);
 
       stringstream name;
       name << "PulseShape";
 
       PulseShape = new TProfile(name.str().c_str(), name.str().c_str(), nBins, -0.5, double(nBins) - 0.5);
 
       // shift shapes to have max close to the real APD max
 
       for (int ii = 0; ii < 50; ii++) {
         shapes[ii] = 0.0;
         PulseShape->Fill(ii, 0.0);
       }
 
       for (unsigned int ii = 0; ii < nBins - 50; ii++) {
         sum_jj = 0.0;
         for (unsigned int jj = 0; jj < ii; jj++) {
           elec_jj = elecShape->GetBinContent(jj + 1);
           laser_iiMinusjj = laserShape->GetBinContent(ii - jj + 1);
           sum_jj += elec_jj * laser_iiMinusjj;
         }
         PulseShape->Fill(ii + 50, sum_jj);
         shapes[ii + 50] = sum_jj;
       }
 
       double scale = PulseShape->GetMaximum();
       shapeCorrection = scale;
 
       if (scale != 0) {
         PulseShape->Scale(1.0 / scale);
         for (unsigned int ii = 0; ii < nBins; ii++) {
           shapesVec.push_back(shapes[ii] / scale);
         }
       }
 
       if (_saveshapes)
         PulseShape->Write();
     }
   }
   ShapeFile->Close();
 
   if (!_saveshapes) {
     stringstream del;
     del << "rm " << shapefile;
     system(del.str().c_str());
   }
 
   return IsMatacqOK;
 }

◆ setGeomEB()

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

Definition at line 1428 of file EcalLaserAnalyzer2.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, digitizers_cfi::strip, hgcalTowerProducer_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;
 
   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 EcalLaserAnalyzer2::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 1461 of file EcalLaserAnalyzer2.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, hgcalTowerProducer_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;
 
   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

◆ _debug

const int EcalLaserAnalyzer2::_debug

private

Definition at line 110 of file EcalLaserAnalyzer2.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 EcalLaserAnalyzer2::_docorpn

private

Definition at line 106 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ _ecalPart

const std::string EcalLaserAnalyzer2::_ecalPart

private

Definition at line 104 of file EcalLaserAnalyzer2.h.

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

◆ _fedid

const int EcalLaserAnalyzer2::_fedid

private

Definition at line 107 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ _firstsample

const unsigned int EcalLaserAnalyzer2::_firstsample

private

Definition at line 88 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ _firstsamplePN

const unsigned int EcalLaserAnalyzer2::_firstsamplePN

private

Definition at line 92 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ _lastsample

const unsigned int EcalLaserAnalyzer2::_lastsample

private

Definition at line 89 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ _lastsamplePN

const unsigned int EcalLaserAnalyzer2::_lastsamplePN

private

Definition at line 93 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ _niter

const unsigned int EcalLaserAnalyzer2::_niter

private

Definition at line 102 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ _noise

const double EcalLaserAnalyzer2::_noise

private

Definition at line 103 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ _nsamples

const unsigned int EcalLaserAnalyzer2::_nsamples

private

Definition at line 86 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ _nsamplesPN

const unsigned int EcalLaserAnalyzer2::_nsamplesPN

private

Definition at line 90 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ _presample

unsigned int EcalLaserAnalyzer2::_presample

private

Definition at line 87 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ _presamplecut

const double EcalLaserAnalyzer2::_presamplecut

private

Definition at line 101 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ _presamplePN

const unsigned int EcalLaserAnalyzer2::_presamplePN

private

Definition at line 91 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _qualpercent

const double EcalLaserAnalyzer2::_qualpercent

private

Definition at line 109 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ _ratiomaxcutlow

const double EcalLaserAnalyzer2::_ratiomaxcutlow

private

Definition at line 100 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _ratiomincuthigh

const double EcalLaserAnalyzer2::_ratiomincuthigh

private

Definition at line 99 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _ratiomincutlow

const double EcalLaserAnalyzer2::_ratiomincutlow

private

Definition at line 98 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _saveallevents

const bool EcalLaserAnalyzer2::_saveallevents

private

Definition at line 108 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ _saveshapes

const bool EcalLaserAnalyzer2::_saveshapes

private

Definition at line 105 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

◆ _timingcuthigh

const unsigned int EcalLaserAnalyzer2::_timingcuthigh

private

Definition at line 95 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _timingcutlow

const unsigned int EcalLaserAnalyzer2::_timingcutlow

private

Definition at line 94 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _timingqualhigh

const unsigned int EcalLaserAnalyzer2::_timingqualhigh

private

Definition at line 97 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _timingquallow

const unsigned int EcalLaserAnalyzer2::_timingquallow

private

Definition at line 96 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ adc

double EcalLaserAnalyzer2::adc[10]

private

Definition at line 197 of file EcalLaserAnalyzer2.h.

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

◆ adcC

int EcalLaserAnalyzer2::adcC[1700]

private

Definition at line 242 of file EcalLaserAnalyzer2.h.

◆ ADCfile

std::string EcalLaserAnalyzer2::ADCfile

private

Definition at line 126 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

◆ ADCFile

TFile* EcalLaserAnalyzer2::ADCFile

private

Definition at line 178 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

◆ adcG

int EcalLaserAnalyzer2::adcG[10]

private

Definition at line 198 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ adcMean

double EcalLaserAnalyzer2::adcMean[1700][10]

private

Definition at line 241 of file EcalLaserAnalyzer2.h.

◆ adcNoPed

double* EcalLaserAnalyzer2::adcNoPed

private

Definition at line 211 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ ADCtrees

TTree* EcalLaserAnalyzer2::ADCtrees[1700]

private

Definition at line 179 of file EcalLaserAnalyzer2.h.

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

◆ APD

double EcalLaserAnalyzer2::APD[6]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ apdAmpl

double EcalLaserAnalyzer2::apdAmpl

private

Definition at line 202 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ apdAmplA

double EcalLaserAnalyzer2::apdAmplA

private

Definition at line 203 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ apdAmplB

double EcalLaserAnalyzer2::apdAmplB

private

Definition at line 204 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDAnal

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

private

Definition at line 220 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDfile

std::string EcalLaserAnalyzer2::APDfile

private

Definition at line 127 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

◆ APDFile

TFile* EcalLaserAnalyzer2::APDFile

private

Definition at line 181 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDFirstAnal

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

private

Definition at line 219 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDoAPDA

double EcalLaserAnalyzer2::APDoAPDA[6]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDoAPDB

double EcalLaserAnalyzer2::APDoAPDB[6]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDoPN

double EcalLaserAnalyzer2::APDoPN[6]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDoPNA

double EcalLaserAnalyzer2::APDoPNA[6]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDoPNB

double EcalLaserAnalyzer2::APDoPNB[6]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDPulse

TAPDPulse* EcalLaserAnalyzer2::APDPulse

private

Definition at line 112 of file EcalLaserAnalyzer2.h.

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

◆ apdRefMap

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

private

Definition at line 166 of file EcalLaserAnalyzer2.h.

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

◆ apdTime

double EcalLaserAnalyzer2::apdTime

private

Definition at line 205 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDtrees

TTree* EcalLaserAnalyzer2::APDtrees[1700]

private

Definition at line 182 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ channelID

int EcalLaserAnalyzer2::channelID

private

Definition at line 152 of file EcalLaserAnalyzer2.h.

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

◆ channelIteratorEE

int EcalLaserAnalyzer2::channelIteratorEE

private

Definition at line 235 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ channelMapEE

std::map<int, int> EcalLaserAnalyzer2::channelMapEE

private

Definition at line 163 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

◆ color

int EcalLaserAnalyzer2::color

private

Definition at line 196 of file EcalLaserAnalyzer2.h.

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

◆ colorref

int EcalLaserAnalyzer2::colorref

private

Definition at line 209 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ colors

std::vector<int> EcalLaserAnalyzer2::colors

private

Definition at line 162 of file EcalLaserAnalyzer2.h.

◆ dccID

int EcalLaserAnalyzer2::dccID

private

Definition at line 154 of file EcalLaserAnalyzer2.h.

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

◆ dccMEM

std::vector<int> EcalLaserAnalyzer2::dccMEM

private

Definition at line 164 of file EcalLaserAnalyzer2.h.

◆ Delta01

TMom* EcalLaserAnalyzer2::Delta01

private

Definition at line 115 of file EcalLaserAnalyzer2.h.

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

◆ Delta12

TMom* EcalLaserAnalyzer2::Delta12

private

Definition at line 116 of file EcalLaserAnalyzer2.h.

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

◆ digiCollection_

const std::string EcalLaserAnalyzer2::digiCollection_

private

Definition at line 74 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ digiPNCollection_

const std::string EcalLaserAnalyzer2::digiPNCollection_

private

Definition at line 76 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ digiProducer_

const std::string EcalLaserAnalyzer2::digiProducer_

private

Definition at line 75 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ ebDigiToken_

edm::EDGetTokenT<EBDigiCollection> EcalLaserAnalyzer2::ebDigiToken_

private

Definition at line 79 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ eeDigiToken_

edm::EDGetTokenT<EEDigiCollection> EcalLaserAnalyzer2::eeDigiToken_

private

Definition at line 80 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ elecfile_

const std::string EcalLaserAnalyzer2::elecfile_

private

Definition at line 119 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

◆ eta

int EcalLaserAnalyzer2::eta

private

Definition at line 194 of file EcalLaserAnalyzer2.h.

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

◆ event

int EcalLaserAnalyzer2::event

private

Definition at line 195 of file EcalLaserAnalyzer2.h.

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

◆ eventHeaderCollection_

const std::string EcalLaserAnalyzer2::eventHeaderCollection_

private

Definition at line 72 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ eventHeaderProducer_

const std::string EcalLaserAnalyzer2::eventHeaderProducer_

private

Definition at line 73 of file EcalLaserAnalyzer2.h.

◆ eventref

int EcalLaserAnalyzer2::eventref

private

Definition at line 208 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ fedID

int EcalLaserAnalyzer2::fedID

private

Definition at line 153 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ firstChanMod

unsigned int EcalLaserAnalyzer2::firstChanMod[21]

private

Definition at line 247 of file EcalLaserAnalyzer2.h.

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

◆ flag

int EcalLaserAnalyzer2::flag

private

Definition at line 234 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ iChannelID

int EcalLaserAnalyzer2::iChannelID[1700]

private

Definition at line 246 of file EcalLaserAnalyzer2.h.

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

◆ idccID

int EcalLaserAnalyzer2::idccID[1700]

private

Definition at line 246 of file EcalLaserAnalyzer2.h.

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

◆ iEta

int EcalLaserAnalyzer2::iEta[1700]

private

Definition at line 244 of file EcalLaserAnalyzer2.h.

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

◆ iEvent

int EcalLaserAnalyzer2::iEvent

private

Definition at line 70 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ iModule

unsigned int EcalLaserAnalyzer2::iModule[1700]

private

Definition at line 245 of file EcalLaserAnalyzer2.h.

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

◆ iPhi

int EcalLaserAnalyzer2::iPhi[1700]

private

Definition at line 244 of file EcalLaserAnalyzer2.h.

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

◆ isFirstChanModFilled

unsigned int EcalLaserAnalyzer2::isFirstChanModFilled[21]

private

Definition at line 248 of file EcalLaserAnalyzer2.h.

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

◆ isGainOK

bool EcalLaserAnalyzer2::isGainOK

private

Definition at line 259 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ iside

int EcalLaserAnalyzer2::iside[1700]

private

Definition at line 246 of file EcalLaserAnalyzer2.h.

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

◆ IsMatacqOK

bool EcalLaserAnalyzer2::IsMatacqOK

private

Definition at line 142 of file EcalLaserAnalyzer2.h.

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

◆ isMatacqOK

bool EcalLaserAnalyzer2::isMatacqOK

private

Definition at line 229 of file EcalLaserAnalyzer2.h.

◆ isSPRFine

bool EcalLaserAnalyzer2::isSPRFine

private

Definition at line 230 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

◆ IsThereDataADC

int EcalLaserAnalyzer2::IsThereDataADC[1700][nColor]

private

Definition at line 222 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ isTimingOK

bool EcalLaserAnalyzer2::isTimingOK

private

Definition at line 260 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ iTowerID

int EcalLaserAnalyzer2::iTowerID[1700]

private

Definition at line 246 of file EcalLaserAnalyzer2.h.

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

◆ iZ

int EcalLaserAnalyzer2::iZ

private

Definition at line 157 of file EcalLaserAnalyzer2.h.

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

◆ laserEvents

int EcalLaserAnalyzer2::laserEvents

private

Definition at line 160 of file EcalLaserAnalyzer2.h.

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

◆ lightside

int EcalLaserAnalyzer2::lightside

private

Definition at line 156 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ mappingToken_

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

private

Definition at line 82 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ matfile

std::string EcalLaserAnalyzer2::matfile

private

Definition at line 125 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and getShapes().

◆ Mem

TMem* EcalLaserAnalyzer2::Mem

private

Definition at line 114 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ moduleID

int EcalLaserAnalyzer2::moduleID

private

Definition at line 234 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ modules

std::vector<int> EcalLaserAnalyzer2::modules

private

Definition at line 165 of file EcalLaserAnalyzer2.h.

Referenced by SequenceVisitors.NodeVisitor::enter().

◆ nCrys

unsigned int EcalLaserAnalyzer2::nCrys

private

Definition at line 133 of file EcalLaserAnalyzer2.h.

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

◆ nEvtBadGain

int EcalLaserAnalyzer2::nEvtBadGain[1700]

private

Definition at line 252 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

◆ nEvtBadTiming

int EcalLaserAnalyzer2::nEvtBadTiming[1700]

private

Definition at line 253 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

◆ nEvtTot

int EcalLaserAnalyzer2::nEvtTot[1700]

private

Definition at line 254 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

◆ nMod

unsigned int EcalLaserAnalyzer2::nMod

private

Definition at line 137 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ nPNPerMod

unsigned int EcalLaserAnalyzer2::nPNPerMod

private

Definition at line 134 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ nRefChan

unsigned int EcalLaserAnalyzer2::nRefChan

private

Definition at line 135 of file EcalLaserAnalyzer2.h.

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

◆ nRefTrees

unsigned int EcalLaserAnalyzer2::nRefTrees

private

Definition at line 136 of file EcalLaserAnalyzer2.h.

◆ nSamplesShapes

unsigned int EcalLaserAnalyzer2::nSamplesShapes

private

Definition at line 140 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

◆ nSides

unsigned int EcalLaserAnalyzer2::nSides

private

Definition at line 138 of file EcalLaserAnalyzer2.h.

◆ phi

int EcalLaserAnalyzer2::phi

private

Definition at line 194 of file EcalLaserAnalyzer2.h.

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

◆ pn

double EcalLaserAnalyzer2::pn[50]

private

Definition at line 200 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ PN

double EcalLaserAnalyzer2::PN[6]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ pn0

double EcalLaserAnalyzer2::pn0

private

Definition at line 199 of file EcalLaserAnalyzer2.h.

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

◆ pn1

double EcalLaserAnalyzer2::pn1

private

Definition at line 199 of file EcalLaserAnalyzer2.h.

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

◆ pnAmpl

double EcalLaserAnalyzer2::pnAmpl

private

Definition at line 206 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ PNAnal

TPN* EcalLaserAnalyzer2::PNAnal[9][2][nColor]

private

Definition at line 216 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ pncorfile_

const std::string EcalLaserAnalyzer2::pncorfile_

private

Definition at line 120 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ pnCorrector

TPNCor* EcalLaserAnalyzer2::pnCorrector

private

Definition at line 170 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ pnDiodeDigiToken_

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

private

Definition at line 81 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ PNFirstAnal

TPN* EcalLaserAnalyzer2::PNFirstAnal[9][2][nColor]

private

Definition at line 215 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ pnG

int EcalLaserAnalyzer2::pnG[50]

private

Definition at line 201 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ pnID

int EcalLaserAnalyzer2::pnID

private

Definition at line 234 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ pnNoPed

double* EcalLaserAnalyzer2::pnNoPed

private

Definition at line 212 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ PNoPN

double EcalLaserAnalyzer2::PNoPN[6]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ PNoPNA

double EcalLaserAnalyzer2::PNoPNA[6]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ PNoPNB

double EcalLaserAnalyzer2::PNoPNB[6]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ PNPulse

TPNPulse* EcalLaserAnalyzer2::PNPulse

private

Definition at line 113 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ PulseShape

TProfile* EcalLaserAnalyzer2::PulseShape

private

Definition at line 190 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

◆ rawDataToken_

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

private

Definition at line 78 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ RefAPDtrees

TTree* EcalLaserAnalyzer2::RefAPDtrees[2][21]

private

Definition at line 183 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ resdir_

const std::string EcalLaserAnalyzer2::resdir_

private

Definition at line 118 of file EcalLaserAnalyzer2.h.

Referenced by beginJob().

◆ resfile

std::string EcalLaserAnalyzer2::resfile

private

Definition at line 128 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

◆ resFile

TFile* EcalLaserAnalyzer2::resFile

private

Definition at line 185 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ respntrees

TTree* EcalLaserAnalyzer2::respntrees[nColor]

private

Definition at line 187 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ restrees

TTree* EcalLaserAnalyzer2::restrees[nColor]

private

Definition at line 186 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ runNum

int EcalLaserAnalyzer2::runNum

private

Definition at line 147 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ runType

int EcalLaserAnalyzer2::runType

private

Definition at line 146 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ ShapeCor

double EcalLaserAnalyzer2::ShapeCor

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ shapeCorrection

double EcalLaserAnalyzer2::shapeCorrection

private

Definition at line 228 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

◆ shapefile

std::string EcalLaserAnalyzer2::shapefile

private

Definition at line 124 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and getShapes().

◆ ShapeFile

TFile* EcalLaserAnalyzer2::ShapeFile

private

Definition at line 189 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

◆ shapes

double EcalLaserAnalyzer2::shapes[250]

private

Definition at line 227 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

◆ shapesVec

std::vector<double> EcalLaserAnalyzer2::shapesVec

private

Definition at line 226 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

◆ side

int EcalLaserAnalyzer2::side

private

Definition at line 155 of file EcalLaserAnalyzer2.h.

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

◆ Time

double EcalLaserAnalyzer2::Time[6]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ towerID

int EcalLaserAnalyzer2::towerID

private

Definition at line 151 of file EcalLaserAnalyzer2.h.

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

◆ wasGainOK

bool EcalLaserAnalyzer2::wasGainOK[1700]

private

Definition at line 256 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ wasTimingOK

bool EcalLaserAnalyzer2::wasTimingOK[1700]

private

Definition at line 257 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ VarCol

Constructor & Destructor Documentation

◆ EcalLaserAnalyzer2()

◆ ~EcalLaserAnalyzer2()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ getShapes()

◆ setGeomEB()

◆ setGeomEE()

Member Data Documentation

◆ _debug

◆ _docorpn

◆ _ecalPart

◆ _fedid

◆ _firstsample

◆ _firstsamplePN

◆ _lastsample

◆ _lastsamplePN

◆ _niter

◆ _noise

◆ _nsamples

◆ _nsamplesPN

◆ _presample

◆ _presamplecut

◆ _presamplePN

◆ _qualpercent

◆ _ratiomaxcutlow

◆ _ratiomincuthigh

◆ _ratiomincutlow

◆ _saveallevents

◆ _saveshapes

◆ _timingcuthigh

◆ _timingcutlow

◆ _timingqualhigh

◆ _timingquallow

◆ adc

◆ adcC

◆ ADCfile

◆ ADCFile

◆ adcG

◆ adcMean

◆ adcNoPed

◆ ADCtrees

◆ APD

◆ apdAmpl

◆ apdAmplA

◆ apdAmplB

◆ APDAnal

◆ APDfile

◆ APDFile

◆ APDFirstAnal

◆ APDoAPDA

◆ APDoAPDB

◆ APDoPN

◆ APDoPNA

◆ APDoPNB

◆ APDPulse

◆ apdRefMap

◆ apdTime

◆ APDtrees

◆ channelID

◆ channelIteratorEE

◆ channelMapEE

◆ color

◆ colorref

◆ colors

◆ dccID

◆ dccMEM

◆ Delta01

◆ Delta12

◆ digiCollection_