#include <EcalTestPulseAnalyzer.h>

Inheritance diagram for EcalTestPulseAnalyzer:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	beginJob () override

	EcalTestPulseAnalyzer (const edm::ParameterSet &iConfig)

void	endJob () override

	~EcalTestPulseAnalyzer () override

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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

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

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
const double	_chi2max

const std::string	_ecalPart

const int	_fedid

const unsigned int	_firstsample

const unsigned int	_firstsamplePN

const unsigned int	_lastsample

const unsigned int	_lastsamplePN

const unsigned int	_niter

const unsigned int	_nsamples

const unsigned int	_nsamplesPN

const unsigned int	_presample

const unsigned int	_presamplePN

const unsigned int	_samplemax

const unsigned int	_samplemin

const double	_timeofmax

double	adc [10]

double	APD [6]

double	apdAmpl

int	apdGain

double	apdTime

int	channelID

int	channelIteratorEE

std::map< int, int >	channelMapEE

int	dccID

std::vector< int >	dccMEM

const std::string	digiCollection_

const std::string	digiPNCollection_

const std::string	digiProducer_

edm::EDGetTokenT< EBDigiCollection >	ebDigiToken_

edm::EDGetTokenT< EEDigiCollection >	eeDigiToken_

int	eta

int	event

const std::string	eventHeaderCollection_

const std::string	eventHeaderProducer_

int	fedID

unsigned int	firstChanMod [9]

int	flag

int	gain

int	iChannelID [1700]

int	idccID [1700]

int	ieta

int	iEta [1700]

int	iEvent

unsigned int	iModule [1700]

int	iphi

int	iPhi [1700]

unsigned int	isFirstChanModFilled [9]

int	iside [1700]

int	iTowerID [1700]

int	iZ

const edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	mappingToken_

int	moduleID

std::vector< int >	modules

unsigned int	nCrys

unsigned int	nGainAPD

unsigned int	nGainPN

unsigned int	nMod

TFile *	outFile

int	phi

double	pn [50]

double	PN [6]

double	pnAmpl

double	pnAmpl0

double	pnAmpl1

const edm::EDGetTokenT< EcalPnDiodeDigiCollection >	pnDiodeDigiToken_

int	pnG

int	pnGain

int	pnID

const edm::EDGetTokenT< EcalRawDataCollection >	rawDataToken_

const std::string	resdir_

std::string	resfile

TFile *	resFile

TTree *	respntrees

TTree *	restrees

double	ret_data [20]

std::string	rootfile

int	runNum

int	runType

int	side

int	towerID

int	TPEvents

TTree *	trees [1700]

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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 44 of file EcalTestPulseAnalyzer.h.

Constructor & Destructor Documentation

◆ EcalTestPulseAnalyzer()

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

explicit

Definition at line 36 of file EcalTestPulseAnalyzer.cc.

References _ecalPart, _fedid, dccMEM, digiCollection_, digiProducer_, ebDigiToken_, eeDigiToken_, ProducerED_cfi::InputTag, iZ, ME::lmmodFromDcc(), ME::memFromDcc(), nCrys, NCRYSEB, NCRYSEE, and nMod.

     : iEvent(0),
       eventHeaderCollection_(iConfig.getParameter<std::string>("eventHeaderCollection")),
       eventHeaderProducer_(iConfig.getParameter<std::string>("eventHeaderProducer")),
       digiCollection_(iConfig.getParameter<std::string>("digiCollection")),
       digiProducer_(iConfig.getParameter<std::string>("digiProducer")),
       digiPNCollection_(iConfig.getParameter<std::string>("digiPNCollection")),
       rawDataToken_(consumes<EcalRawDataCollection>(edm::InputTag(eventHeaderProducer_, eventHeaderCollection_))),
       pnDiodeDigiToken_(consumes<EcalPnDiodeDigiCollection>(edm::InputTag(digiProducer_, digiPNCollection_))),
       mappingToken_(esConsumes()),
       // framework parameters with default values
       _nsamples(iConfig.getUntrackedParameter<unsigned int>("nSamples", 10)),
       _presample(iConfig.getUntrackedParameter<unsigned int>("nPresamples", 3)),
       _firstsample(iConfig.getUntrackedParameter<unsigned int>("firstSample", 1)),
       _lastsample(iConfig.getUntrackedParameter<unsigned int>("lastSample", 2)),
       _samplemin(iConfig.getUntrackedParameter<unsigned int>("sampleMin", 3)),
       _samplemax(iConfig.getUntrackedParameter<unsigned int>("sampleMax", 9)),
       _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)),
       _niter(iConfig.getUntrackedParameter<unsigned int>("nIter", 3)),
       _chi2max(iConfig.getUntrackedParameter<double>("chi2Max", 10.0)),
       _timeofmax(iConfig.getUntrackedParameter<double>("timeOfMax", 4.5)),
       _ecalPart(iConfig.getUntrackedParameter<std::string>("ecalPart", "EB")),
       _fedid(iConfig.getUntrackedParameter<int>("fedID", -999)),
       resdir_(iConfig.getUntrackedParameter<std::string>("resDir")),
       nCrys(NCRYSEB),
       nMod(NMODEB),
       nGainPN(NGAINPN),
       nGainAPD(NGAINAPD),
       towerID(-1),
       channelID(-1),
       runType(-1),
       runNum(0),
       fedID(-1),
       dccID(-1),
       side(-1),
       iZ(1),
       phi(-1),
       eta(-1),
       event(0),
       apdAmpl(0),
       apdTime(0),
       pnAmpl(0),
       pnID(-1),
       moduleID(-1),
       channelIteratorEE(0)
 
 //========================================================================
 
 {
   if (_ecalPart == "EB") {
     ebDigiToken_ = consumes<EBDigiCollection>(edm::InputTag(digiProducer_, digiCollection_));
   } else if (_ecalPart == "EE") {
     eeDigiToken_ = consumes<EEDigiCollection>(edm::InputTag(digiProducer_, digiCollection_));
   }
 
   // Define geometrical constants
 
   if (_ecalPart == "EB") {
     nCrys = NCRYSEB;
   } else {
     nCrys = NCRYSEE;
   }
 
   iZ = 1;
   if (_fedid <= 609)
     iZ = -1;
 
   dccMEM = ME::memFromDcc(_fedid);
   modules = ME::lmmodFromDcc(_fedid);
   nMod = modules.size();
 }

◆ ~EcalTestPulseAnalyzer()

EcalTestPulseAnalyzer::~EcalTestPulseAnalyzer ( )

override

Definition at line 113 of file EcalTestPulseAnalyzer.cc.

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

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 194 of file EcalTestPulseAnalyzer.cc.

References _chi2max, _ecalPart, _fedid, _firstsample, _firstsamplePN, _lastsample, _lastsamplePN, _niter, _nsamples, _nsamplesPN, _presample, _presamplePN, _samplemax, _samplemin, _timeofmax, adc, apdAmpl, apdGain, apdTime, cms::cuda::assert(), edm::SortedCollection< T, SORT >::begin(), HltBtagPostValidation_cff::c, EcalElectronicsId::channelId(), channelID, channelIteratorEE, channelMapEE, nano_mu_local_reco_cff::chi2, dccID, dccMEM, hgcalPerformanceValidation::df, digiCollection_, TPNFit::doFit(), MillePedeFileConverter_cfg::e, CollectionTags_cfi::EBDigi, ebDigiToken_, CollectionTags_cfi::EEDigi, eeDigiToken_, MEEBGeom::electronic_channel(), edm::SortedCollection< T, SORT >::end(), eta, eventHeaderCollection_, fedID, firstChanMod, TSFit::fit_third_degree_polynomial(), TPNFit::getAmpl(), ecalpyutils::hashedIndex(), mps_fire::i, iChannelID, EcalPnDiodeDetId::iDCCId(), idccID, iEta, iEvent, iModule, TPNFit::init(), TSFit::init_errmat(), iPhi, isFirstChanModFilled, iside, edm::HandleBase::isValid(), iTowerID, iZ, dqmdumpme::k, MEEEGeom::lmmod(), MEEBGeom::lmmod(), MEEBGeom::localCoord(), mappingToken_, METSkim_cff::Max, callgraph::module, nCrys, phi, MEEBGeom::pn(), MEEEGeom::pn(), pn, pnAmpl, pnAmpl0, pnAmpl1, pnDiodeDigiToken_, pnG, PedestalClient_cfi::pngain, pnGain, edm::Handle< T >::product(), rawDataToken_, ret_data, runNum, runType, TSFit::set_params(), MEEBGeom::side(), MEEEGeom::side(), side, findQualityFiles::size, nano_mu_digi_cff::strip, EcalElectronicsId::stripId(), EcalDCCHeaderBlock::TESTPULSE_GAP, EcalDCCHeaderBlock::TESTPULSE_MGPA, EcalDCCHeaderBlock::TESTPULSE_SCAN_MEM, EcalElectronicsId::towerId(), towerID, TPEvents, trees, and EcalElectronicsId::xtalId().

                                                                              {
   //========================================================================
 
   ++iEvent;
 
   // Retrieve DCC header
   edm::Handle<EcalRawDataCollection> pDCCHeader;
   const EcalRawDataCollection* DCCHeader = nullptr;
   e.getByToken(rawDataToken_, pDCCHeader);
   if (!pDCCHeader.isValid()) {
     edm::LogError("nodata") << "Error! can't get the product  retrieving DCC header" << eventHeaderCollection_.c_str();
   } else {
     DCCHeader = pDCCHeader.product();
   }
 
   // retrieving crystal data from Event
   edm::Handle<EBDigiCollection> pEBDigi;
   const EBDigiCollection* EBDigi = nullptr;
   edm::Handle<EEDigiCollection> pEEDigi;
   const EEDigiCollection* EEDigi = nullptr;
   if (_ecalPart == "EB") {
     e.getByToken(ebDigiToken_, pEBDigi);
     if (!pEBDigi.isValid()) {
       edm::LogError("nodata") << "Error! can't get the product retrieving EB crystal data " << digiCollection_.c_str();
     } else {
       EBDigi = pEBDigi.product();
     }
   } else {
     e.getByToken(eeDigiToken_, pEEDigi);
     if (!pEEDigi.isValid()) {
       edm::LogError("nodata") << "Error! can't get the product retrieving EE crystal data " << digiCollection_.c_str();
     } else {
       EEDigi = pEEDigi.product();
     }
   }
 
   // Retrieve 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 " << digiCollection_.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) {
     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();
 
     if (600 + dccID != fedID)
       continue;
 
     // Cut on runType
 
     if (runType != EcalDCCHeaderBlock::TESTPULSE_MGPA && runType != EcalDCCHeaderBlock::TESTPULSE_GAP &&
         runType != EcalDCCHeaderBlock::TESTPULSE_SCAN_MEM)
       return;
   }
 
   // Cut on fedID
 
   if (fedID != _fedid && _fedid != -999)
     return;
 
   // Count TP events
   TPEvents++;
 
   // ======================
   // Decode PN Information
   // ======================
 
   TPNFit* pnfit = new TPNFit();
   pnfit->init(_nsamplesPN, _firstsamplePN, _lastsamplePN);
 
   double chi2pn = 0;
   double ypnrange[50];
   double dsum = 0.;
   double dsum1 = 0.;
   double bl = 0.;
   double val_max = 0.;
   int samplemax = 0;
   unsigned int k;
   int pnG[50];
   int pngain = 0;
 
   std::map<int, std::vector<double> > allPNAmpl;
   std::map<int, std::vector<int> > allPNGain;
 
   for (EcalPnDiodeDigiCollection::const_iterator pnItr = PNDigi->begin(); pnItr != PNDigi->end();
        ++pnItr) {  // Loop on PNs
 
     EcalPnDiodeDetId pnDetId = EcalPnDiodeDetId((*pnItr).id());
 
     bool isMemRelevant = false;
     for (unsigned int imem = 0; imem < dccMEM.size(); imem++) {
       if (pnDetId.iDCCId() == dccMEM[imem]) {
         isMemRelevant = true;
       }
     }
 
     // skip mem dcc without relevant data
     if (!isMemRelevant)
       continue;
 
     for (int samId = 0; samId < (*pnItr).size(); samId++) {  // Loop on PN samples
       pn[samId] = (*pnItr).sample(samId).adc();
       pnG[samId] = (*pnItr).sample(samId).gainId();
 
       if (pnG[samId] != 1)
         edm::LogVerbatim("EcalTestPulseAnalyzer") << "PN gain different from 1 for sample " << samId;
       if (samId == 0)
         pngain = pnG[samId];
       if (samId > 0)
         pngain = TMath::Max(pnG[samId], pngain);
     }
 
     for (dsum = 0., k = 0; k < _presamplePN; k++) {
       dsum += pn[k];
     }
     bl = dsum / ((double)_presamplePN);
 
     for (val_max = 0., k = 0; k < _nsamplesPN; k++) {
       ypnrange[k] = pn[k] - bl;
 
       if (ypnrange[k] > val_max) {
         val_max = ypnrange[k];
         samplemax = k;
       }
     }
 
     chi2pn = pnfit->doFit(samplemax, &ypnrange[0]);
 
     if (chi2pn == 101 || chi2pn == 102 || chi2pn == 103)
       pnAmpl = 0.;
     else
       pnAmpl = pnfit->getAmpl();
 
     allPNAmpl[pnDetId.iDCCId()].push_back(pnAmpl);
     allPNGain[pnDetId.iDCCId()].push_back(pngain);
   }
 
   // ===========================
   // Decode EBDigis Information
   // ===========================
 
   TSFit* pstpfit = new TSFit(_nsamples, 650);
   pstpfit->set_params(
       _nsamples, _niter, _presample, _samplemin, _samplemax, _timeofmax, _chi2max, _firstsample, _lastsample);
   pstpfit->init_errmat(10.);
 
   double chi2 = 0;
   double yrange[10];
   int adcgain = 0;
   int adcG[10];
 
   if (EBDigi) {
     for (EBDigiCollection::const_iterator digiItr = EBDigi->begin(); digiItr != EBDigi->end();
          ++digiItr) {  // Loop on EB crystals
 
       EBDetId id_crystal(digiItr->id());
       EBDataFrame df(*digiItr);
 
       int etaG = id_crystal.ieta();  // global
       int phiG = id_crystal.iphi();  // global
 
       int etaL;  // local
       int phiL;  // local
       std::pair<int, int> LocalCoord = MEEBGeom::localCoord(etaG, phiG);
 
       etaL = LocalCoord.first;
       phiL = LocalCoord.second;
 
       eta = etaG;
       phi = phiG;
 
       side = MEEBGeom::side(etaG, phiG);
 
       EcalElectronicsId elecid_crystal = TheMapping.getElectronicsId(id_crystal);
 
       towerID = elecid_crystal.towerId();
       int strip = elecid_crystal.stripId();
       int xtal = elecid_crystal.xtalId();
       channelID = 5 * (strip - 1) + xtal - 1;  // FIXME
 
       int module = MEEBGeom::lmmod(etaG, phiG);
       int iMod = module - 1;
 
       assert(module >= *min_element(modules.begin(), modules.end()) &&
              module <= *max_element(modules.begin(), modules.end()));
 
       std::pair<int, int> pnpair = MEEBGeom::pn(module);
       unsigned int MyPn0 = pnpair.first;
       unsigned int MyPn1 = pnpair.second;
 
       unsigned int channel = MEEBGeom::electronic_channel(etaL, phiL);
 
       if (isFirstChanModFilled[iMod] == 0) {
         firstChanMod[iMod] = channel;
         isFirstChanModFilled[iMod] = 1;
       }
 
       iEta[channel] = eta;
       iPhi[channel] = phi;
       iModule[channel] = module;
       iTowerID[channel] = towerID;
       iChannelID[channel] = channelID;
       idccID[channel] = dccID;
       iside[channel] = side;
 
       // get 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();
 
         if (i == 0)
           adcgain = adcG[i];
         if (i > 0)
           adcgain = TMath::Max(adcG[i], adcgain);
       }
       // Remove pedestal
       //====================
       for (dsum = 0., dsum1 = 0., k = 0; k < _presample; k++) {
         dsum += adc[k];
         if (k < _presample - 1)
           dsum1 += adc[k];
       }
 
       bl = dsum / ((double)_presample);
 
       for (val_max = 0., k = 0; k < _nsamples; k++) {
         yrange[k] = adc[k] - bl;
         if (yrange[k] > val_max) {
           val_max = yrange[k];
         }
       }
 
       apdGain = adcgain;
 
       if (allPNAmpl[dccMEM[0]].size() > MyPn0)
         pnAmpl0 = allPNAmpl[dccMEM[0]][MyPn0];
       else
         pnAmpl0 = 0;
       if (allPNAmpl[dccMEM[0]].size() > MyPn1)
         pnAmpl1 = allPNAmpl[dccMEM[0]][MyPn1];
       else
         pnAmpl1 = 0;
 
       if (allPNGain[dccMEM[0]].size() > MyPn0)
         pnGain = allPNGain[dccMEM[0]][MyPn0];
       else
         pnGain = 0;
 
       // Perform the fit on apd samples
       //================================
 
       chi2 = pstpfit->fit_third_degree_polynomial(&yrange[0], ret_data);
 
       //Retrieve APD amplitude from fit
       //================================
 
       if (val_max > 100000. || chi2 < 0. || chi2 == 102) {
         apdAmpl = 0;
         apdTime = 0;
 
       } else {
         apdAmpl = ret_data[0];
         apdTime = ret_data[1];
       }
 
       trees[channel]->Fill();
     }
 
   } else {
     for (EEDigiCollection::const_iterator digiItr = EEDigi->begin(); digiItr != EEDigi->end();
          ++digiItr) {  // Loop on EE crystals
 
       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);
 
       // Recover the TT id and the electronic crystal numbering from EcalElectronicsMapping
 
       EcalElectronicsId elecid_crystal = TheMapping.getElectronicsId(id_crystal);
 
       towerID = elecid_crystal.towerId();
       channelID = elecid_crystal.channelId() - 1;
 
       int module = MEEEGeom::lmmod(iX, iY);
       if (module >= 18 && side == 1)
         module += 2;  // Trick to fix endcap specificity
       int iMod = module - 1;
 
       assert(module >= *min_element(modules.begin(), modules.end()) &&
              module <= *max_element(modules.begin(), modules.end()));
 
       std::pair<int, int> pnpair = MEEEGeom::pn(module, _fedid);
 
       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];
 
       if (isFirstChanModFilled[iMod] == 0) {
         firstChanMod[iMod] = channel;
         isFirstChanModFilled[iMod] = 1;
       }
 
       iEta[channel] = eta;
       iPhi[channel] = phi;
       iModule[channel] = module;
       iTowerID[channel] = towerID;
       iChannelID[channel] = channelID;
       idccID[channel] = dccID;
       iside[channel] = side;
 
       assert(channel < nCrys);
 
       // Get 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();
 
         if (i == 0)
           adcgain = adcG[i];
         if (i > 0)
           adcgain = TMath::Max(adcG[i], adcgain);
       }
 
       // Remove pedestal
       //====================
       for (dsum = 0., dsum1 = 0., k = 0; k < _presample; k++) {
         dsum += adc[k];
         if (k < _presample - 1)
           dsum1 += adc[k];
       }
 
       bl = dsum / ((double)_presample);
 
       for (val_max = 0., k = 0; k < _nsamples; k++) {
         yrange[k] = adc[k] - bl;
         if (yrange[k] > val_max) {
           val_max = yrange[k];
         }
       }
       apdGain = adcgain;
 
       int dccMEMIndex = 0;
       if (side == 1)
         dccMEMIndex += 2;  // Trick to fix endcap specificity
 
       if (allPNAmpl[dccMEM[dccMEMIndex]].size() > MyPn0)
         pnAmpl0 = allPNAmpl[dccMEM[dccMEMIndex]][MyPn0];
       else
         pnAmpl0 = 0;
       if (allPNAmpl[dccMEM[dccMEMIndex + 1]].size() > MyPn1)
         pnAmpl1 = allPNAmpl[dccMEM[dccMEMIndex + 1]][MyPn1];
       else
         pnAmpl1 = 0;
 
       if (allPNGain[dccMEM[dccMEMIndex]].size() > MyPn0)
         pnGain = allPNGain[dccMEM[dccMEMIndex]][MyPn0];
       else
         pnGain = 0;
 
       // Perform the fit on apd samples
       //=================================
 
       chi2 = pstpfit->fit_third_degree_polynomial(&yrange[0], ret_data);
 
       //Retrieve APD amplitude from fit
       //=================================
 
       if (val_max > 100000. || chi2 < 0. || chi2 == 102) {
         apdAmpl = 0;
         apdTime = 0;
 
       } else {
         apdAmpl = ret_data[0];
         apdTime = ret_data[1];
       }
 
       trees[channel]->Fill();
     }
   }
 
 }  // end of analyze

◆ beginJob()

void EcalTestPulseAnalyzer::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 121 of file EcalTestPulseAnalyzer.cc.

References apdAmpl, apdGain, channelID, dccID, eta, firstChanMod, mps_fire::i, iChannelID, idccID, iEta, iModule, iPhi, isFirstChanModFilled, iside, iTowerID, dqmiolumiharvest::j, Skims_PA_cff::name, nCrys, nMod, outFile, phi, pnAmpl0, pnAmpl1, pnGain, resdir_, resfile, rootfile, side, towerID, TPEvents, and trees.

                                      {
   //========================================================================
 
   // Define temporary file
 
   rootfile = resdir_;
   rootfile += "/TmpTreeTestPulseAnalyzer.root";
 
   outFile = new TFile(rootfile.c_str(), "RECREATE");
 
   for (unsigned int i = 0; i < nCrys; i++) {
     std::stringstream name;
     name << "Tree" << i;
 
     trees[i] = new TTree(name.str().c_str(), name.str().c_str());
 
     //List of branches
 
     trees[i]->Branch("iphi", &phi, "phi/I");
     trees[i]->Branch("ieta", &eta, "eta/I");
     trees[i]->Branch("side", &side, "side/I");
     trees[i]->Branch("dccID", &dccID, "dccID/I");
     trees[i]->Branch("towerID", &towerID, "towerID/I");
     trees[i]->Branch("channelID", &channelID, "channelID/I");
     trees[i]->Branch("event", &event, "event/I");
     trees[i]->Branch("apdGain", &apdGain, "apdGain/I");
     trees[i]->Branch("pnGain", &pnGain, "pnGain/I");
     trees[i]->Branch("apdAmpl", &apdAmpl, "apdAmpl/D");
     trees[i]->Branch("pnAmpl0", &pnAmpl0, "pnAmpl0/D");
     trees[i]->Branch("pnAmpl1", &pnAmpl1, "pnAmpl1/D");
 
     trees[i]->SetBranchAddress("ieta", &eta);
     trees[i]->SetBranchAddress("iphi", &phi);
     trees[i]->SetBranchAddress("side", &side);
     trees[i]->SetBranchAddress("dccID", &dccID);
     trees[i]->SetBranchAddress("towerID", &towerID);
     trees[i]->SetBranchAddress("channelID", &channelID);
     trees[i]->SetBranchAddress("event", &event);
     trees[i]->SetBranchAddress("apdGain", &apdGain);
     trees[i]->SetBranchAddress("pnGain", &pnGain);
     trees[i]->SetBranchAddress("apdAmpl", &apdAmpl);
     trees[i]->SetBranchAddress("pnAmpl0", &pnAmpl0);
     trees[i]->SetBranchAddress("pnAmpl1", &pnAmpl1);
   }
 
   // Initializations
 
   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;
   }
 
   for (unsigned int j = 0; j < nMod; j++) {
     firstChanMod[j] = 0;
     isFirstChanModFilled[j] = 0;
   }
 
   // Define output results file name
 
   std::stringstream namefile;
   namefile << resdir_ << "/APDPN_TESTPULSE.root";
   resfile = namefile.str();
 
   // TP events counter
   TPEvents = 0;
 }

◆ endJob()

void EcalTestPulseAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 621 of file EcalTestPulseAnalyzer.cc.

References TMom::addEntry(), APD, apdAmpl, apdGain, channelID, dccID, firstChanMod, flag, gain, TMom::getM3(), TMom::getMax(), TMom::getMean(), TMom::getMin(), TMom::getNevt(), TMom::getRMS(), iChannelID, idccID, ieta, iEta, iModule, iphi, iPhi, iside, iTowerID, dqmiolumiharvest::j, moduleID, nCrys, nGainAPD, nGainPN, nMod, outFile, PN, pnAmpl0, pnAmpl1, pnGain, pnID, resfile, resFile, respntrees, restrees, rootfile, side, towerID, TPEvents, and trees.

Referenced by o2olib.O2ORunMgr::executeJob().

                                    {
   //========================================================================
 
   // Don't do anything if there is no events
   if (TPEvents == 0) {
     outFile->Close();
 
     // Remove temporary file
 
     std::stringstream del;
     del << "rm " << rootfile;
     system(del.str().c_str());
 
     edm::LogVerbatim("EcalTestPulseAnalyzer") << " No TP Events ";
     return;
   }
 
   edm::LogVerbatim("EcalTestPulseAnalyzer") << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
   edm::LogVerbatim("EcalTestPulseAnalyzer") << "\t+=+     Analyzing test pulse data: getting APD, PN  +=+";
 
   // Create output ntuples:
 
   //edm::LogVerbatim("EcalTestPulseAnalyzer")<< "TP Test Name File "<< resfile.c_str();
 
   resFile = new TFile(resfile.c_str(), "RECREATE");
 
   restrees = new TTree("TPAPD", "TPAPD");
   respntrees = new TTree("TPPN", "TPPN");
 
   restrees->Branch("iphi", &iphi, "iphi/I");
   restrees->Branch("ieta", &ieta, "ieta/I");
   restrees->Branch("dccID", &dccID, "dccID/I");
   restrees->Branch("side", &side, "side/I");
   restrees->Branch("towerID", &towerID, "towerID/I");
   restrees->Branch("channelID", &channelID, "channelID/I");
   restrees->Branch("moduleID", &moduleID, "moduleID/I");
   restrees->Branch("flag", &flag, "flag/I");
   restrees->Branch("gain", &gain, "gain/I");
   restrees->Branch("APD", &APD, "APD[6]/D");
 
   respntrees->Branch("pnID", &pnID, "pnID/I");
   respntrees->Branch("moduleID", &moduleID, "moduleID/I");
   respntrees->Branch("gain", &gain, "gain/I");
   respntrees->Branch("PN", &PN, "PN[6]/D");
 
   restrees->SetBranchAddress("iphi", &iphi);
   restrees->SetBranchAddress("ieta", &ieta);
   restrees->SetBranchAddress("dccID", &dccID);
   restrees->SetBranchAddress("side", &side);
   restrees->SetBranchAddress("towerID", &towerID);
   restrees->SetBranchAddress("channelID", &channelID);
   restrees->SetBranchAddress("moduleID", &moduleID);
   restrees->SetBranchAddress("flag", &flag);
   restrees->SetBranchAddress("gain", &gain);
   restrees->SetBranchAddress("APD", APD);
 
   respntrees->SetBranchAddress("pnID", &pnID);
   respntrees->SetBranchAddress("moduleID", &moduleID);
   respntrees->SetBranchAddress("gain", &gain);
   respntrees->SetBranchAddress("PN", PN);
 
   TMom* APDAnal[1700][10];
   TMom* PNAnal[9][2][10];
 
   for (unsigned int iMod = 0; iMod < nMod; iMod++) {
     for (unsigned int ich = 0; ich < 2; ich++) {
       for (unsigned int ig = 0; ig < nGainPN; ig++) {
         PNAnal[iMod][ich][ig] = new TMom();
       }
     }
   }
 
   for (unsigned int iCry = 0; iCry < nCrys; iCry++) {  // Loop on data trees (ie on cristals)
 
     for (unsigned int iG = 0; iG < nGainAPD; iG++) {
       APDAnal[iCry][iG] = new TMom();
     }
 
     // Define submodule and channel number inside the submodule (as Patrice)
 
     unsigned int iMod = iModule[iCry] - 1;
 
     moduleID = iMod + 1;
     if (moduleID >= 20)
       moduleID -= 2;  // Trick to fix endcap specificity
 
     for (Long64_t jentry = 0; jentry < trees[iCry]->GetEntriesFast(); jentry++) {
       // PN Means and RMS
 
       if (firstChanMod[iMod] == iCry) {
         PNAnal[iMod][0][pnGain]->addEntry(pnAmpl0);
         PNAnal[iMod][1][pnGain]->addEntry(pnAmpl1);
       }
 
       // APD means and RMS
 
       APDAnal[iCry][apdGain]->addEntry(apdAmpl);
     }
 
     if (trees[iCry]->GetEntries() < 10) {
       flag = -1;
       for (int j = 0; j < 6; j++) {
         APD[j] = 0.0;
       }
     } else
       flag = 1;
 
     iphi = iPhi[iCry];
     ieta = iEta[iCry];
     dccID = idccID[iCry];
     side = iside[iCry];
     towerID = iTowerID[iCry];
     channelID = iChannelID[iCry];
 
     for (unsigned int ig = 0; ig < nGainAPD; ig++) {
       APD[0] = APDAnal[iCry][ig]->getMean();
       APD[1] = APDAnal[iCry][ig]->getRMS();
       APD[2] = APDAnal[iCry][ig]->getM3();
       APD[3] = APDAnal[iCry][ig]->getNevt();
       APD[4] = APDAnal[iCry][ig]->getMin();
       APD[5] = APDAnal[iCry][ig]->getMax();
       gain = ig;
 
       // Fill APD tree
 
       restrees->Fill();
     }
   }
 
   // Get final results for PN and PN/PN
 
   for (unsigned int ig = 0; ig < nGainPN; ig++) {
     for (unsigned int iMod = 0; iMod < nMod; iMod++) {
       for (int ch = 0; ch < 2; ch++) {
         pnID = ch;
         moduleID = iMod;
         if (moduleID >= 20)
           moduleID -= 2;  // Trick to fix endcap specificity
 
         PN[0] = PNAnal[iMod][ch][ig]->getMean();
         PN[1] = PNAnal[iMod][ch][ig]->getRMS();
         PN[2] = PNAnal[iMod][ch][ig]->getM3();
         PN[3] = PNAnal[iMod][ch][ig]->getNevt();
         PN[4] = PNAnal[iMod][ch][ig]->getMin();
         PN[5] = PNAnal[iMod][ch][ig]->getMax();
         gain = ig;
 
         // Fill PN tree
         respntrees->Fill();
       }
     }
   }
 
   outFile->Close();
 
   // Remove temporary file
 
   std::stringstream del;
   del << "rm " << rootfile;
   system(del.str().c_str());
 
   // Save final results
 
   restrees->Write();
   respntrees->Write();
   resFile->Close();
 
   edm::LogVerbatim("EcalTestPulseAnalyzer") << "\t+=+    ...................................... done  +=+";
   edm::LogVerbatim("EcalTestPulseAnalyzer") << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
 }

Member Data Documentation

◆ _chi2max

const double EcalTestPulseAnalyzer::_chi2max

private

Definition at line 81 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _ecalPart

const std::string EcalTestPulseAnalyzer::_ecalPart

private

Definition at line 83 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

◆ _fedid

const int EcalTestPulseAnalyzer::_fedid

private

Definition at line 84 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

◆ _firstsample

const unsigned int EcalTestPulseAnalyzer::_firstsample

private

Definition at line 72 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _firstsamplePN

const unsigned int EcalTestPulseAnalyzer::_firstsamplePN

private

Definition at line 78 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _lastsample

const unsigned int EcalTestPulseAnalyzer::_lastsample

private

Definition at line 73 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _lastsamplePN

const unsigned int EcalTestPulseAnalyzer::_lastsamplePN

private

Definition at line 79 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _niter

const unsigned int EcalTestPulseAnalyzer::_niter

private

Definition at line 80 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _nsamples

const unsigned int EcalTestPulseAnalyzer::_nsamples

private

Definition at line 70 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _nsamplesPN

const unsigned int EcalTestPulseAnalyzer::_nsamplesPN

private

Definition at line 76 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _presample

const unsigned int EcalTestPulseAnalyzer::_presample

private

Definition at line 71 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _presamplePN

const unsigned int EcalTestPulseAnalyzer::_presamplePN

private

Definition at line 77 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _samplemax

const unsigned int EcalTestPulseAnalyzer::_samplemax

private

Definition at line 75 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _samplemin

const unsigned int EcalTestPulseAnalyzer::_samplemin

private

Definition at line 74 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ _timeofmax

const double EcalTestPulseAnalyzer::_timeofmax

private

Definition at line 82 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ adc

double EcalTestPulseAnalyzer::adc[10]

private

Definition at line 131 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ APD

double EcalTestPulseAnalyzer::APD[6]

private

Definition at line 157 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ apdAmpl

double EcalTestPulseAnalyzer::apdAmpl

private

Definition at line 137 of file EcalTestPulseAnalyzer.h.

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

◆ apdGain

int EcalTestPulseAnalyzer::apdGain

private

Definition at line 134 of file EcalTestPulseAnalyzer.h.

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

◆ apdTime

double EcalTestPulseAnalyzer::apdTime

private

Definition at line 138 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ channelID

int EcalTestPulseAnalyzer::channelID

private

Definition at line 107 of file EcalTestPulseAnalyzer.h.

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

◆ channelIteratorEE

int EcalTestPulseAnalyzer::channelIteratorEE

private

Definition at line 156 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ channelMapEE

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

private

Definition at line 148 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ dccID

int EcalTestPulseAnalyzer::dccID

private

Definition at line 114 of file EcalTestPulseAnalyzer.h.

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

◆ dccMEM

std::vector<int> EcalTestPulseAnalyzer::dccMEM

private

Definition at line 149 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

◆ digiCollection_

const std::string EcalTestPulseAnalyzer::digiCollection_

private

Definition at line 58 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

◆ digiPNCollection_

const std::string EcalTestPulseAnalyzer::digiPNCollection_

private

Definition at line 60 of file EcalTestPulseAnalyzer.h.

◆ digiProducer_

const std::string EcalTestPulseAnalyzer::digiProducer_

private

Definition at line 59 of file EcalTestPulseAnalyzer.h.

Referenced by EcalTestPulseAnalyzer().

◆ ebDigiToken_

edm::EDGetTokenT<EBDigiCollection> EcalTestPulseAnalyzer::ebDigiToken_

private

Definition at line 63 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

◆ eeDigiToken_

edm::EDGetTokenT<EEDigiCollection> EcalTestPulseAnalyzer::eeDigiToken_

private

Definition at line 64 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

◆ eta

int EcalTestPulseAnalyzer::eta

private

Definition at line 129 of file EcalTestPulseAnalyzer.h.

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

◆ event

int EcalTestPulseAnalyzer::event

private

Definition at line 130 of file EcalTestPulseAnalyzer.h.

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

◆ eventHeaderCollection_

const std::string EcalTestPulseAnalyzer::eventHeaderCollection_

private

Definition at line 56 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ eventHeaderProducer_

const std::string EcalTestPulseAnalyzer::eventHeaderProducer_

private

Definition at line 57 of file EcalTestPulseAnalyzer.h.

◆ fedID

int EcalTestPulseAnalyzer::fedID

private

Definition at line 113 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ firstChanMod

unsigned int EcalTestPulseAnalyzer::firstChanMod[9]

private

Definition at line 163 of file EcalTestPulseAnalyzer.h.

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

◆ flag

int EcalTestPulseAnalyzer::flag

private

Definition at line 154 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ gain

int EcalTestPulseAnalyzer::gain

private

Definition at line 154 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ iChannelID

int EcalTestPulseAnalyzer::iChannelID[1700]

private

Definition at line 161 of file EcalTestPulseAnalyzer.h.

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

◆ idccID

int EcalTestPulseAnalyzer::idccID[1700]

private

Definition at line 161 of file EcalTestPulseAnalyzer.h.

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

◆ ieta

int EcalTestPulseAnalyzer::ieta

private

Definition at line 154 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ iEta

int EcalTestPulseAnalyzer::iEta[1700]

private

Definition at line 159 of file EcalTestPulseAnalyzer.h.

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

◆ iEvent

int EcalTestPulseAnalyzer::iEvent

private

Definition at line 54 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ iModule

unsigned int EcalTestPulseAnalyzer::iModule[1700]

private

Definition at line 160 of file EcalTestPulseAnalyzer.h.

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

◆ iphi

int EcalTestPulseAnalyzer::iphi

private

Definition at line 154 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ iPhi

int EcalTestPulseAnalyzer::iPhi[1700]

private

Definition at line 159 of file EcalTestPulseAnalyzer.h.

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

◆ isFirstChanModFilled

unsigned int EcalTestPulseAnalyzer::isFirstChanModFilled[9]

private

Definition at line 164 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ iside

int EcalTestPulseAnalyzer::iside[1700]

private

Definition at line 161 of file EcalTestPulseAnalyzer.h.

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

◆ iTowerID

int EcalTestPulseAnalyzer::iTowerID[1700]

private

Definition at line 161 of file EcalTestPulseAnalyzer.h.

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

◆ iZ

int EcalTestPulseAnalyzer::iZ

private

Definition at line 116 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

◆ mappingToken_

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

private

Definition at line 66 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ moduleID

int EcalTestPulseAnalyzer::moduleID

private

Definition at line 155 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ modules

std::vector<int> EcalTestPulseAnalyzer::modules

private

Definition at line 150 of file EcalTestPulseAnalyzer.h.

Referenced by SequenceVisitors.NodeVisitor::enter().

◆ nCrys

unsigned int EcalTestPulseAnalyzer::nCrys

private

Definition at line 96 of file EcalTestPulseAnalyzer.h.

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

◆ nGainAPD

unsigned int EcalTestPulseAnalyzer::nGainAPD

private

Definition at line 99 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ nGainPN

unsigned int EcalTestPulseAnalyzer::nGainPN

private

Definition at line 98 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ nMod

unsigned int EcalTestPulseAnalyzer::nMod

private

Definition at line 97 of file EcalTestPulseAnalyzer.h.

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

◆ outFile

TFile* EcalTestPulseAnalyzer::outFile

private

Definition at line 120 of file EcalTestPulseAnalyzer.h.

◆ phi

int EcalTestPulseAnalyzer::phi

private

Definition at line 129 of file EcalTestPulseAnalyzer.h.

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

◆ pn

double EcalTestPulseAnalyzer::pn[50]

private

Definition at line 132 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ PN

double EcalTestPulseAnalyzer::PN[6]

private

Definition at line 157 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ pnAmpl

double EcalTestPulseAnalyzer::pnAmpl

private

Definition at line 141 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ pnAmpl0

double EcalTestPulseAnalyzer::pnAmpl0

private

Definition at line 139 of file EcalTestPulseAnalyzer.h.

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

◆ pnAmpl1

double EcalTestPulseAnalyzer::pnAmpl1

private

Definition at line 140 of file EcalTestPulseAnalyzer.h.

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

◆ pnDiodeDigiToken_

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

private

Definition at line 65 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ pnG

int EcalTestPulseAnalyzer::pnG

private

Definition at line 136 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ pnGain

int EcalTestPulseAnalyzer::pnGain

private

Definition at line 135 of file EcalTestPulseAnalyzer.h.

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

◆ pnID

int EcalTestPulseAnalyzer::pnID

private

Definition at line 155 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ rawDataToken_

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

private

Definition at line 62 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ resdir_

const std::string EcalTestPulseAnalyzer::resdir_

private

Definition at line 86 of file EcalTestPulseAnalyzer.h.

Referenced by beginJob().

◆ resfile

std::string EcalTestPulseAnalyzer::resfile

private

Definition at line 91 of file EcalTestPulseAnalyzer.h.

Referenced by beginJob(), and endJob().

◆ resFile

TFile* EcalTestPulseAnalyzer::resFile

private

Definition at line 121 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ respntrees

TTree* EcalTestPulseAnalyzer::respntrees

private

Definition at line 146 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ restrees

TTree* EcalTestPulseAnalyzer::restrees

private

Definition at line 145 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

◆ ret_data

double EcalTestPulseAnalyzer::ret_data[20]

private

Definition at line 104 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ rootfile

std::string EcalTestPulseAnalyzer::rootfile

private

Definition at line 90 of file EcalTestPulseAnalyzer.h.

Referenced by beginJob(), endJob(), and SimpleTreeProducer.SimpleTreeProducer::write().

◆ runNum

int EcalTestPulseAnalyzer::runNum

private

Definition at line 112 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ runType

int EcalTestPulseAnalyzer::runType

private

Definition at line 111 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

◆ side

int EcalTestPulseAnalyzer::side

private

Definition at line 115 of file EcalTestPulseAnalyzer.h.

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

◆ towerID

int EcalTestPulseAnalyzer::towerID

private

Definition at line 106 of file EcalTestPulseAnalyzer.h.

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

◆ TPEvents

int EcalTestPulseAnalyzer::TPEvents

private

Definition at line 102 of file EcalTestPulseAnalyzer.h.

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

◆ trees

TTree* EcalTestPulseAnalyzer::trees[1700]

private

Definition at line 125 of file EcalTestPulseAnalyzer.h.

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

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalTestPulseAnalyzer()

◆ ~EcalTestPulseAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

Member Data Documentation

◆ _chi2max

◆ _ecalPart

◆ _fedid

◆ _firstsample

◆ _firstsamplePN

◆ _lastsample

◆ _lastsamplePN

◆ _niter

◆ _nsamples

◆ _nsamplesPN

◆ _presample

◆ _presamplePN

◆ _samplemax

◆ _samplemin

◆ _timeofmax

◆ adc

◆ APD

◆ apdAmpl

◆ apdGain

◆ apdTime

◆ channelID

◆ channelIteratorEE

◆ channelMapEE

◆ dccID

◆ dccMEM

◆ digiCollection_

◆ digiPNCollection_

◆ digiProducer_

◆ ebDigiToken_

◆ eeDigiToken_

◆ eta

◆ event

◆ eventHeaderCollection_

◆ eventHeaderProducer_

◆ fedID

◆ firstChanMod

◆ flag

◆ gain

◆ iChannelID

◆ idccID

◆ ieta

◆ iEta

◆ iEvent

◆ iModule

◆ iphi

◆ iPhi

◆ isFirstChanModFilled

◆ iside

◆ iTowerID

◆ iZ

◆ mappingToken_

◆ moduleID

◆ modules

◆ nCrys

◆ nGainAPD

◆ nGainPN

◆ nMod

◆ outFile

◆ phi

◆ pn

◆ PN

◆ pnAmpl

◆ pnAmpl0

◆ pnAmpl1

◆ pnDiodeDigiToken_

◆ pnG

◆ pnGain

◆ pnID