#include <EcalABAnalyzer.h>

Inheritance diagram for EcalABAnalyzer:

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

	EcalABAnalyzer (const edm::ParameterSet &iConfig)

void	endJob () override

	~EcalABAnalyzer () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
const double	_alpha

const double	_beta

const double	_chi2cut

const int	_debug

const std::string	_ecalPart

const int	_fedid

const unsigned int	_firstsample

bool	_fitab

const unsigned int	_lastsample

const unsigned int	_nevtmax

const unsigned int	_niter

const double	_noise

const unsigned int	_nsamples

unsigned int	_presample

const double	_presamplecut

const double	_qualpercent

const double	_ratiomaxcutlow

const double	_ratiomincuthigh

const double	_ratiomincutlow

const unsigned int	_timingcuthigh

const unsigned int	_timingcutlow

const unsigned int	_timingqualhigh

const unsigned int	_timingquallow

double	adc [10]

int	adcG [10]

std::string	alphafile

std::string	alphainitfile

TAPDPulse *	APDPulse

int	channelIteratorEE

std::map< int, int >	channelMapEE

int	color

std::vector< int >	colors

int	dccID

std::vector< int >	dccMEM

TMom *	Delta01

TMom *	Delta12

const std::string	digiCollection_

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	fedID

int	iChannelID [1700]

int	idccID [1700]

int	iEta [1700]

int	iEvent

int	iPhi [1700]

bool	isGainOK

int	iside [1700]

bool	isTimingOK

int	iTowerID [1700]

int	iZ

int	lightside

const edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	mappingToken_

std::vector< int >	modules

unsigned int	nCrys

unsigned int	nevtAB [1700]

int	nEvtBadGain [1700]

int	nEvtBadTiming [1700]

int	nEvtTot [1700]

int	phi

const edm::EDGetTokenT< EcalRawDataCollection >	rawDataToken_

const std::string	resdir_

int	runNum

int	runType

TShapeAnalysis *	shapana

int	side

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 40 of file EcalABAnalyzer.h.

Member Enumeration Documentation

◆ VarCol

enum EcalABAnalyzer::VarCol

Enumerator
iBlue
iRed
nColor

Definition at line 49 of file EcalABAnalyzer.h.

49 { iBlue, iRed, nColor };

EcalABAnalyzer::iBlue

Definition: EcalABAnalyzer.h:49

EcalABAnalyzer::nColor

Definition: EcalABAnalyzer.h:49

EcalABAnalyzer::iRed

Definition: EcalABAnalyzer.h:49

Constructor & Destructor Documentation

◆ EcalABAnalyzer()

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

explicit

Definition at line 41 of file EcalABAnalyzer.cc.

References _ecalPart, _fedid, _firstsample, _fitab, _lastsample, _nsamples, _presample, _ratiomaxcutlow, _ratiomincuthigh, _ratiomincutlow, _timingcuthigh, _timingcutlow, _timingqualhigh, _timingquallow, APDPulse, Delta01, Delta12, digiCollection_, digiProducer_, ebDigiToken_, eeDigiToken_, iChannelID, idccID, iEta, HLT_2022v11_cff::InputTag, iPhi, isGainOK, iside, isTimingOK, iTowerID, iZ, dqmiolumiharvest::j, nCrys, NCRYSEB, NCRYSEE, nevtAB, 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")),
       rawDataToken_(consumes<EcalRawDataCollection>(edm::InputTag(eventHeaderProducer_, eventHeaderCollection_))),
       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)),
       _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)),
       _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", 100.0)),
       _ecalPart(iConfig.getUntrackedParameter<std::string>("ecalPart", "EB")),
       _fedid(iConfig.getUntrackedParameter<int>("fedId", -999)),
       _qualpercent(iConfig.getUntrackedParameter<double>("qualPercent", 0.2)),
       _debug(iConfig.getUntrackedParameter<int>("debug", 0)),
       resdir_(iConfig.getUntrackedParameter<std::string>("resDir")),
       nCrys(NCRYSEB),
       runType(-1),
       runNum(0),
       fedID(-1),
       dccID(-1),
       side(2),
       lightside(2),
       iZ(1),
       phi(-1),
       eta(-1),
       event(0),
       color(-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;
 
   for (unsigned int j = 0; j < nCrys; j++) {
     iEta[j] = -1;
     iPhi[j] = -1;
     iTowerID[j] = -1;
     iChannelID[j] = -1;
     idccID[j] = -1;
     iside[j] = -1;
     wasTimingOK[j] = true;
     wasGainOK[j] = true;
     nevtAB[j] = 0;
   }
 
   // Quality check flags
 
   isGainOK = true;
   isTimingOK = true;
 
   // Objects dealing with pulses
 
   APDPulse = new TAPDPulse(_nsamples,
                            _presample,
                            _firstsample,
                            _lastsample,
                            _timingcutlow,
                            _timingcuthigh,
                            _timingquallow,
                            _timingqualhigh,
                            _ratiomincutlow,
                            _ratiomincuthigh,
                            _ratiomaxcutlow);
 
   // Objects needed for npresample calculation
 
   Delta01 = new TMom();
   Delta12 = new TMom();
   _fitab = true;
 }

◆ ~EcalABAnalyzer()

EcalABAnalyzer::~EcalABAnalyzer ( )

override

Definition at line 147 of file EcalABAnalyzer.cc.

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

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 224 of file EcalABAnalyzer.cc.

References _ecalPart, _fedid, _fitab, _nevtmax, adc, adcG, TMom::addEntry(), APDPulse, cms::cuda::assert(), edm::SortedCollection< T, SORT >::begin(), c, EcalElectronicsId::channelId(), channelIteratorEE, channelMapEE, color, dccID, Delta01, Delta12, hgcalPerformanceValidation::df, digiCollection_, doesABTreeExist, 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(), TAPDPulse::getDelta(), ecalpyutils::hashedIndex(), mps_fire::i, iChannelID, idccID, iEta, iEvent, iPhi, iside, TAPDPulse::isPulseOK(), TAPDPulse::isTimingQualOK(), edm::HandleBase::isValid(), iTowerID, iZ, EcalDCCHeaderBlock::LASER_DELAY_SCAN, EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_POWER_SCAN, EcalDCCHeaderBlock::LASER_STD, lightside, MEEBGeom::localCoord(), mappingToken_, METSkim_cff::Max, nCrys, nevtAB, nEvtBadGain, nEvtBadTiming, nEvtTot, phi, edm::Handle< T >::product(), TShapeAnalysis::putAllVals(), rawDataToken_, runNum, runType, TAPDPulse::setPulse(), shapana, MEEBGeom::side(), MEEEGeom::side(), side, digitizers_cfi::strip, EcalElectronicsId::stripId(), ecaldqm::towerID(), EcalElectronicsId::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 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);
     }
   }
 
   // Cut on fedID
 
   if (fedID != _fedid && _fedid != -999)
     return;
 
   // ===========================
   // 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);
 
       int etaG = id_crystal.ieta();  // global
       int phiG = id_crystal.iphi();  // global
 
       int etaL;  // local
       int phiL;  // local
       std::pair<int, int> LocalCoord = MEEBGeom::localCoord(etaG, phiG);
 
       etaL = LocalCoord.first;
       phiL = LocalCoord.second;
 
       eta = etaG;
       phi = phiG;
 
       side = MEEBGeom::side(etaG, phiG);
 
       // Recover the TT id and the electronic crystal numbering from EcalElectronicsMapping
 
       EcalElectronicsId elecid_crystal = TheMapping.getElectronicsId(id_crystal);
 
       int towerID = elecid_crystal.towerId();
       int strip = elecid_crystal.stripId();
       int xtal = elecid_crystal.xtalId();
       int channelID = 5 * (strip - 1) + xtal - 1;
 
       unsigned int channel = MEEBGeom::electronic_channel(etaL, phiL);
 
       assert(channel < nCrys);
 
       iEta[channel] = eta;
       iPhi[channel] = phi;
       iTowerID[channel] = towerID;
       iChannelID[channel] = channelID;
       idccID[channel] = dccID;
       iside[channel] = side;
 
       // 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]++;
 
       // Fill if Pulse is fine
       //=======================
 
       if (APDPulse->isPulseOK() && lightside == side) {
         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) {  // Loop on crystals
 
       // Retrieve geometry
       //===================
 
       EEDetId id_crystal(digiItr->id());
       EEDataFrame df(*digiItr);
 
       phi = id_crystal.ix();
       eta = id_crystal.iy();
 
       int iX = (phi - 1) / 5 + 1;
       int iY = (eta - 1) / 5 + 1;
 
       side = MEEEGeom::side(iX, iY, iZ);
       EcalElectronicsId elecid_crystal = TheMapping.getElectronicsId(id_crystal);
 
       int towerID = elecid_crystal.towerId();
       int channelID = elecid_crystal.channelId() - 1;
 
       int hashedIndex = 100000 * eta + phi;
 
       if (channelMapEE.count(hashedIndex) == 0) {
         channelMapEE[hashedIndex] = channelIteratorEE;
         channelIteratorEE++;
       }
 
       unsigned int channel = channelMapEE[hashedIndex];
 
       assert(channel < nCrys);
 
       iEta[channel] = eta;
       iPhi[channel] = phi;
       iTowerID[channel] = towerID;
       iChannelID[channel] = channelID;
       idccID[channel] = dccID;
       iside[channel] = side;
 
       // APD Pulse
       //===========
 
       if ((*digiItr).size() > 10)
         edm::LogVerbatim("EcalABAnalyzer") << "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]++;
 
       // Fill if Pulse is fine
       //=======================
 
       if (APDPulse->isPulseOK() && lightside == side) {
         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]++;
         }
       }
     }
   }
 }  // analyze

◆ beginJob()

void EcalABAnalyzer::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 155 of file EcalABAnalyzer.cc.

References _alpha, _beta, _chi2cut, _firstsample, _fitab, _lastsample, _nevtmax, _noise, _nsamples, _presample, alphafile, alphainitfile, doesABTreeExist, MainPageGenerator::link, resdir_, TShapeAnalysis::set_const(), shapana, and ctpps_dqm_sourceclient-live_cfg::test.

                               {
   //========================================================================
 
   //Calculate alpha and beta
 
   // 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;
   std::stringstream link;
   nameabfile << resdir_ << "/AB.root";
 
   FILE* test;
   test = fopen(nameabinitfile.str().c_str(), "r");
   if (test == nullptr) {
     doesABTreeExist = false;
     _fitab = true;
   } else {
     fclose(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();
 
   if (_fitab) {
     alphafile = nameabfile.str();
   } else {
     alphafile = alphainitfile;
     link << "ln -s " << resdir_ << "/ABInit.root " << resdir_ << "/AB.root";
     system(link.str().c_str());
   }
 
   // Define output results files' names
 
   std::stringstream namefile;
   namefile << resdir_ << "/AB.root";
   alphafile = namefile.str();
 }

◆ endJob()

void EcalABAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 508 of file EcalABAnalyzer.cc.

References _ecalPart, _fitab, _presample, _presamplecut, _qualpercent, alphafile, alphainitfile, APDPulse, channelMapEE, TShapeAnalysis::computeShape(), Delta01, Delta12, doesABTreeExist, TMom::getMean(), mps_fire::i, isGainOK, isTimingOK, nCrys, nEvtBadGain, nEvtBadTiming, nEvtTot, TShapeAnalysis::set_nch(), TShapeAnalysis::set_presample(), TAPDPulse::setPresamples(), shapana, wasGainOK, and wasTimingOK.

Referenced by o2olib.O2ORunMgr::executeJob().

                             {
   //========================================================================
 
   edm::LogVerbatim("EcalABAnalyzer") << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
   edm::LogVerbatim("EcalABAnalyzer") << "\t+=+     Analyzing data: getting (alpha, beta)     +=+";
 
   // 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("EcalABAnalyzer") << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
     edm::LogVerbatim("EcalABAnalyzer") << "\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);
 
     // 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;
 
     if (!isGainOK)
       edm::LogVerbatim("EcalABAnalyzer") << "\t+=+ ............................ WARNING! APD GAIN WAS NOT 1    +=+";
     if (!isTimingOK)
       edm::LogVerbatim("EcalABAnalyzer") << "\t+=+ ............................ WARNING! TIMING WAS BAD        +=+";
 
     edm::LogVerbatim("EcalABAnalyzer") << "\t+=+    .................................... done  +=+";
     edm::LogVerbatim("EcalABAnalyzer") << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
   }
 }

Member Data Documentation

◆ _alpha

const double EcalABAnalyzer::_alpha

private

Definition at line 79 of file EcalABAnalyzer.h.

Referenced by beginJob().

◆ _beta

const double EcalABAnalyzer::_beta

private

Definition at line 80 of file EcalABAnalyzer.h.

Referenced by beginJob().

◆ _chi2cut

const double EcalABAnalyzer::_chi2cut

private

Definition at line 83 of file EcalABAnalyzer.h.

Referenced by beginJob().

◆ _debug

const int EcalABAnalyzer::_debug

private

Definition at line 87 of file EcalABAnalyzer.h.

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

◆ _ecalPart

const std::string EcalABAnalyzer::_ecalPart

private

Definition at line 84 of file EcalABAnalyzer.h.

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

◆ _fedid

const int EcalABAnalyzer::_fedid

private

Definition at line 85 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ _firstsample

const unsigned int EcalABAnalyzer::_firstsample

private

Definition at line 68 of file EcalABAnalyzer.h.

Referenced by beginJob(), and EcalABAnalyzer().

◆ _fitab

bool EcalABAnalyzer::_fitab

private

Definition at line 109 of file EcalABAnalyzer.h.

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

◆ _lastsample

const unsigned int EcalABAnalyzer::_lastsample

private

Definition at line 69 of file EcalABAnalyzer.h.

Referenced by beginJob(), and EcalABAnalyzer().

◆ _nevtmax

const unsigned int EcalABAnalyzer::_nevtmax

private

Definition at line 81 of file EcalABAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ _niter

const unsigned int EcalABAnalyzer::_niter

private

Definition at line 78 of file EcalABAnalyzer.h.

◆ _noise

const double EcalABAnalyzer::_noise

private

Definition at line 82 of file EcalABAnalyzer.h.

Referenced by beginJob().

◆ _nsamples

const unsigned int EcalABAnalyzer::_nsamples

private

Definition at line 66 of file EcalABAnalyzer.h.

Referenced by beginJob(), and EcalABAnalyzer().

◆ _presample

unsigned int EcalABAnalyzer::_presample

private

Definition at line 67 of file EcalABAnalyzer.h.

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

◆ _presamplecut

const double EcalABAnalyzer::_presamplecut

private

Definition at line 77 of file EcalABAnalyzer.h.

Referenced by endJob().

◆ _qualpercent

const double EcalABAnalyzer::_qualpercent

private

Definition at line 86 of file EcalABAnalyzer.h.

Referenced by endJob().

◆ _ratiomaxcutlow

const double EcalABAnalyzer::_ratiomaxcutlow

private

Definition at line 76 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

◆ _ratiomincuthigh

const double EcalABAnalyzer::_ratiomincuthigh

private

Definition at line 75 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

◆ _ratiomincutlow

const double EcalABAnalyzer::_ratiomincutlow

private

Definition at line 74 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

◆ _timingcuthigh

const unsigned int EcalABAnalyzer::_timingcuthigh

private

Definition at line 71 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

◆ _timingcutlow

const unsigned int EcalABAnalyzer::_timingcutlow

private

Definition at line 70 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

◆ _timingqualhigh

const unsigned int EcalABAnalyzer::_timingqualhigh

private

Definition at line 73 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

◆ _timingquallow

const unsigned int EcalABAnalyzer::_timingquallow

private

Definition at line 72 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

◆ adc

double EcalABAnalyzer::adc[10]

private

Definition at line 132 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ adcG

int EcalABAnalyzer::adcG[10]

private

Definition at line 133 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ alphafile

std::string EcalABAnalyzer::alphafile

private

Definition at line 97 of file EcalABAnalyzer.h.

Referenced by beginJob(), and endJob().

◆ alphainitfile

std::string EcalABAnalyzer::alphainitfile

private

Definition at line 98 of file EcalABAnalyzer.h.

Referenced by beginJob(), and endJob().

◆ APDPulse

TAPDPulse* EcalABAnalyzer::APDPulse

private

Definition at line 89 of file EcalABAnalyzer.h.

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

◆ channelIteratorEE

int EcalABAnalyzer::channelIteratorEE

private

Definition at line 134 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ channelMapEE

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

private

Definition at line 123 of file EcalABAnalyzer.h.

Referenced by analyze(), and endJob().

◆ color

int EcalABAnalyzer::color

private

Definition at line 131 of file EcalABAnalyzer.h.

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

◆ colors

std::vector<int> EcalABAnalyzer::colors

private

Definition at line 122 of file EcalABAnalyzer.h.

◆ dccID

int EcalABAnalyzer::dccID

private

Definition at line 115 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ dccMEM

std::vector<int> EcalABAnalyzer::dccMEM

private

Definition at line 124 of file EcalABAnalyzer.h.

◆ Delta01

TMom* EcalABAnalyzer::Delta01

private

Definition at line 90 of file EcalABAnalyzer.h.

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

◆ Delta12

TMom* EcalABAnalyzer::Delta12

private

Definition at line 91 of file EcalABAnalyzer.h.

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

◆ digiCollection_

const std::string EcalABAnalyzer::digiCollection_

private

Definition at line 56 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ digiProducer_

const std::string EcalABAnalyzer::digiProducer_

private

Definition at line 57 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

◆ doesABTreeExist

bool EcalABAnalyzer::doesABTreeExist

private

Definition at line 107 of file EcalABAnalyzer.h.

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

◆ ebDigiToken_

edm::EDGetTokenT<EBDigiCollection> EcalABAnalyzer::ebDigiToken_

private

Definition at line 60 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ eeDigiToken_

edm::EDGetTokenT<EEDigiCollection> EcalABAnalyzer::eeDigiToken_

private

Definition at line 61 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ eta

int EcalABAnalyzer::eta

private

Definition at line 129 of file EcalABAnalyzer.h.

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

◆ event

int EcalABAnalyzer::event

private

Definition at line 130 of file EcalABAnalyzer.h.

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

◆ eventHeaderCollection_

const std::string EcalABAnalyzer::eventHeaderCollection_

private

Definition at line 54 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ eventHeaderProducer_

const std::string EcalABAnalyzer::eventHeaderProducer_

private

Definition at line 55 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ fedID

int EcalABAnalyzer::fedID

private

Definition at line 114 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ iChannelID

int EcalABAnalyzer::iChannelID[1700]

private

Definition at line 137 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ idccID

int EcalABAnalyzer::idccID[1700]

private

Definition at line 137 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ iEta

int EcalABAnalyzer::iEta[1700]

private

Definition at line 136 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ iEvent

int EcalABAnalyzer::iEvent

private

Definition at line 52 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ iPhi

int EcalABAnalyzer::iPhi[1700]

private

Definition at line 136 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ isGainOK

bool EcalABAnalyzer::isGainOK

private

Definition at line 148 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer(), and endJob().

◆ iside

int EcalABAnalyzer::iside[1700]

private

Definition at line 137 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ isTimingOK

bool EcalABAnalyzer::isTimingOK

private

Definition at line 149 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer(), and endJob().

◆ iTowerID

int EcalABAnalyzer::iTowerID[1700]

private

Definition at line 137 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ iZ

int EcalABAnalyzer::iZ

private

Definition at line 118 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ lightside

int EcalABAnalyzer::lightside

private

Definition at line 117 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ mappingToken_

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

private

Definition at line 62 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ modules

std::vector<int> EcalABAnalyzer::modules

private

Definition at line 125 of file EcalABAnalyzer.h.

Referenced by SequenceVisitors.NodeVisitor::enter().

◆ nCrys

unsigned int EcalABAnalyzer::nCrys

private

Definition at line 106 of file EcalABAnalyzer.h.

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

◆ nevtAB

unsigned int EcalABAnalyzer::nevtAB[1700]

private

Definition at line 101 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

◆ nEvtBadGain

int EcalABAnalyzer::nEvtBadGain[1700]

private

Definition at line 141 of file EcalABAnalyzer.h.

Referenced by analyze(), and endJob().

◆ nEvtBadTiming

int EcalABAnalyzer::nEvtBadTiming[1700]

private

Definition at line 142 of file EcalABAnalyzer.h.

Referenced by analyze(), and endJob().

◆ nEvtTot

int EcalABAnalyzer::nEvtTot[1700]

private

Definition at line 143 of file EcalABAnalyzer.h.

Referenced by analyze(), and endJob().

◆ phi

int EcalABAnalyzer::phi

private

Definition at line 129 of file EcalABAnalyzer.h.

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

◆ rawDataToken_

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

private

Definition at line 59 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ resdir_

const std::string EcalABAnalyzer::resdir_

private

Definition at line 93 of file EcalABAnalyzer.h.

Referenced by beginJob().

◆ runNum

int EcalABAnalyzer::runNum

private

Definition at line 113 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ runType

int EcalABAnalyzer::runType

private

Definition at line 112 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ shapana

TShapeAnalysis* EcalABAnalyzer::shapana

private

Definition at line 100 of file EcalABAnalyzer.h.

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

◆ side

int EcalABAnalyzer::side

private

Definition at line 116 of file EcalABAnalyzer.h.

Referenced by analyze().

◆ wasGainOK

bool EcalABAnalyzer::wasGainOK[1700]

private

Definition at line 145 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer(), and endJob().

◆ wasTimingOK

bool EcalABAnalyzer::wasTimingOK[1700]

private

Definition at line 146 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer(), and endJob().

Public Types

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ VarCol

Constructor & Destructor Documentation

◆ EcalABAnalyzer()

◆ ~EcalABAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

Member Data Documentation

◆ _alpha

◆ _beta

◆ _chi2cut

◆ _debug

◆ _ecalPart

◆ _fedid

◆ _firstsample

◆ _fitab

◆ _lastsample

◆ _nevtmax

◆ _niter

◆ _noise

◆ _nsamples

◆ _presample

◆ _presamplecut

◆ _qualpercent

◆ _ratiomaxcutlow

◆ _ratiomincuthigh

◆ _ratiomincutlow

◆ _timingcuthigh

◆ _timingcutlow

◆ _timingqualhigh

◆ _timingquallow

◆ adc

◆ adcG

◆ alphafile

◆ alphainitfile

◆ APDPulse

◆ channelIteratorEE

◆ channelMapEE

◆ color

◆ colors

◆ dccID

◆ dccMEM

◆ Delta01

◆ Delta12

◆ digiCollection_

◆ digiProducer_

◆ doesABTreeExist

◆ ebDigiToken_

◆ eeDigiToken_

◆ eta

◆ event

◆ eventHeaderCollection_

◆ eventHeaderProducer_

◆ fedID

◆ iChannelID

◆ idccID

◆ iEta

◆ iEvent

◆ iPhi

◆ isGainOK

◆ iside

◆ isTimingOK

◆ iTowerID

◆ iZ

◆ lightside

◆ mappingToken_

◆ modules

◆ nCrys

◆ nevtAB

◆ nEvtBadGain

◆ nEvtBadTiming

◆ nEvtTot

◆ phi