#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 noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
const 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 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, ProducerED_cfi::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(), DummyCfis::c, EcalElectronicsId::channelId(), channelIteratorEE, channelMapEE, color, dccID, Delta01, Delta12, das-up-to-nevents::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, nano_mu_digi_cff::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, 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 tier0.Tier0Handler::_curlQueryTier0(), FrontierCondition_GT_autoExpress_cfi.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(), datamodel.Object::DeltaR(), Jet.Jet::jetID(), datamodel.Object::p4(), 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(), and postprocessor.PostProcessor::run().

◆ 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(), datamodel.Object::DeltaR(), datamodel.Object::p4(), 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