EcalTestPulseAnalyzer Class Reference

#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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () 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
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::vector< ModuleDescription const * > &modules, 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
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
double	_chi2max
std::string	_ecalPart
int	_fedid
unsigned int	_firstsample
unsigned int	_firstsamplePN
unsigned int	_lastsample
unsigned int	_lastsamplePN
unsigned int	_niter
unsigned int	_nsamples
unsigned int	_nsamplesPN
unsigned int	_presample
unsigned int	_presamplePN
unsigned int	_samplemax
unsigned int	_samplemin
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
std::string	digiCollection_
std::string	digiPNCollection_
std::string	digiProducer_
int	eta
int	event
std::string	eventHeaderCollection_
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
int	moduleID
std::vector< int >	modules
unsigned int	nCrys
unsigned int	nGainAPD
unsigned int	nGainPN
unsigned int	nMod
unsigned int	nTT
TFile *	outFile
int	phi
double	pn [50]
double	PN [6]
double	pnAmpl
double	pnAmpl0
double	pnAmpl1
int	pnG
int	pnGain
int	pnID
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::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
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<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)

Detailed Description

Definition at line 35 of file EcalTestPulseAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 45 of file EcalTestPulseAnalyzer.cc.

References _ecalPart, _fedid, dccMEM, digiCollection_, digiPNCollection_, digiProducer_, eventHeaderCollection_, eventHeaderProducer_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), iZ, ME::lmmodFromDcc(), ME::memFromDcc(), nCrys, NCRYSEB, NCRYSEE, nMod, nTT, NTTEB, NTTEE, resdir_, and AlCaHLTBitMon_QueryRunRegistry::string.

    : iEvent(0),

      // 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)),
      nCrys(NCRYSEB),
      nTT(NTTEB),
      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)

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

{
  //now do what ever initialization is needed

  resdir_ = iConfig.getUntrackedParameter<std::string>("resDir");

  digiCollection_ = iConfig.getParameter<std::string>("digiCollection");
  digiPNCollection_ = iConfig.getParameter<std::string>("digiPNCollection");
  digiProducer_ = iConfig.getParameter<std::string>("digiProducer");

  eventHeaderCollection_ = iConfig.getParameter<std::string>("eventHeaderCollection");
  eventHeaderProducer_ = iConfig.getParameter<std::string>("eventHeaderProducer");

  // Define geometrical constants

  if (_ecalPart == "EB") {
    nCrys = NCRYSEB;
    nTT = NTTEB;
  } else {
    nCrys = NCRYSEE;
    nTT = NTTEE;
  }

  iZ = 1;
  if (_fedid <= 609)
    iZ = -1;

  dccMEM = ME::memFromDcc(_fedid);
  modules = ME::lmmodFromDcc(_fedid);
  nMod = modules.size();
}

EcalTestPulseAnalyzer::~EcalTestPulseAnalyzer ( )

override

Definition at line 123 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

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

override

Definition at line 204 of file EcalTestPulseAnalyzer.cc.

References _chi2max, _ecalPart, _fedid, _firstsample, _firstsamplePN, _lastsample, _lastsamplePN, _niter, _nsamples, _nsamplesPN, _presample, _presamplePN, _samplemax, _samplemin, _timeofmax, adc, apdAmpl, apdGain, apdTime, edm::SortedCollection< T, SORT >::begin(), edm::DataFrameContainer::begin(), beam_dqm_sourceclient-live_cfg::cerr, EcalElectronicsId::channelId(), channelID, channelIteratorEE, channelMapEE, hltPixelTracks_cff::chi2, gather_cfg::cout, dccID, dccMEM, digiCollection_, digiPNCollection_, digiProducer_, TPNFit::doFit(), CollectionTags_cfi::EBDigi, CollectionTags_cfi::EEDigi, MEEBGeom::electronic_channel(), edm::SortedCollection< T, SORT >::end(), edm::DataFrameContainer::end(), eta, eventHeaderCollection_, eventHeaderProducer_, cppFunctionSkipper::exception, fedID, firstChanMod, TSFit::fit_third_degree_polynomial(), edm::EventSetup::get(), TPNFit::getAmpl(), edm::Event::getByLabel(), EcalElectronicsMapping::getElectronicsId(), ecalpyutils::hashedIndex(), mps_fire::i, iChannelID, EcalPnDiodeDetId::iDCCId(), idccID, iEta, iEvent, iModule, TPNFit::init(), TSFit::init_errmat(), iPhi, isFirstChanModFilled, iside, iTowerID, iZ, dqmdumpme::k, MEEEGeom::lmmod(), MEEBGeom::lmmod(), MEEBGeom::localCoord(), Max(), nCrys, phi, MEEBGeom::pn(), MEEEGeom::pn(), pn, pnAmpl, pnAmpl0, pnAmpl1, pnG, PedestalClient_cfi::pngain, pnGain, edm::Handle< T >::product(), edm::ESHandle< T >::product(), ret_data, runNum, runType, TSFit::set_params(), MEEBGeom::side(), MEEEGeom::side(), side, findQualityFiles::size, digitizers_cfi::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;
  try {
    e.getByLabel(eventHeaderProducer_, eventHeaderCollection_, pDCCHeader);
    DCCHeader = pDCCHeader.product();
  } catch (std::exception& ex) {
    std::cerr << "Error! can't get the product  retrieving DCC header" << eventHeaderCollection_.c_str() << std::endl;
  }

  // retrieving crystal EB data from Event
  edm::Handle<EBDigiCollection> pEBDigi;
  const EBDigiCollection* EBDigi = nullptr;

  // retrieving crystal EE data from Event
  edm::Handle<EEDigiCollection> pEEDigi;
  const EEDigiCollection* EEDigi = nullptr;

  if (_ecalPart == "EB") {
    try {
      e.getByLabel(digiProducer_, digiCollection_, pEBDigi);
      EBDigi = pEBDigi.product();
    } catch (std::exception& ex) {
      std::cerr << "Error! can't get the product retrieving EB crystal data " << digiCollection_.c_str() << std::endl;
    }
  } else {
    try {
      e.getByLabel(digiProducer_, digiCollection_, pEEDigi);
      EEDigi = pEEDigi.product();
    } catch (std::exception& ex) {
      std::cerr << "Error! can't get the product retrieving EE crystal data " << digiCollection_.c_str() << std::endl;
    }
  }

  // Retrieve crystal PN diodes from Event
  edm::Handle<EcalPnDiodeDigiCollection> pPNDigi;
  const EcalPnDiodeDigiCollection* PNDigi = nullptr;
  try {
    e.getByLabel(digiProducer_, digiPNCollection_, pPNDigi);
    PNDigi = pPNDigi.product();
  } catch (std::exception& ex) {
    std::cerr << "Error! can't get the product " << digiCollection_.c_str() << std::endl;
  }

  // retrieving electronics mapping
  edm::ESHandle<EcalElectronicsMapping> ecalmapping;
  const EcalElectronicsMapping* TheMapping = nullptr;
  try {
    c.get<EcalMappingRcd>().get(ecalmapping);
    TheMapping = ecalmapping.product();
  } catch (std::exception& ex) {
    std::cerr << "Error! can't get the product EcalMappingRcd" << std::endl;
  }

  // ====================================
  // 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)
        std::cout << "PN gain different from 1 for sample " << samId << std::endl;
      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];
          samplemax = 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];
          samplemax = 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

void EcalTestPulseAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 131 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;
}

void EcalTestPulseAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 642 of file EcalTestPulseAnalyzer.cc.

References TMom::addEntry(), APD, apdAmpl, apdGain, channelID, gather_cfg::cout, dccID, DEFINE_FWK_MODULE, 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());

    std::cout << " No TP Events " << std::endl;
    return;
  }

  std::cout << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+" << std::endl;
  std::cout << "\t+=+     Analyzing test pulse data: getting APD, PN  +=+" << std::endl;

  // Create output ntuples:

  //std::cout<< "TP Test Name File "<< resfile.c_str() << std::endl;

  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

    Long64_t nbytes = 0, nb = 0;
    for (Long64_t jentry = 0; jentry < trees[iCry]->GetEntriesFast(); jentry++) {
      nb = trees[iCry]->GetEntry(jentry);
      nbytes += nb;

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

  std::cout << "\t+=+    ...................................... done  +=+" << std::endl;
  std::cout << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+" << std::endl;
}

Member Data Documentation

double EcalTestPulseAnalyzer::_chi2max

private

Definition at line 60 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

std::string EcalTestPulseAnalyzer::_ecalPart

private

Definition at line 62 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

int EcalTestPulseAnalyzer::_fedid

private

Definition at line 63 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

unsigned int EcalTestPulseAnalyzer::_firstsample

private

Definition at line 51 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::_firstsamplePN

private

Definition at line 57 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::_lastsample

private

Definition at line 52 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::_lastsamplePN

private

Definition at line 58 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::_niter

private

Definition at line 59 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::_nsamples

private

Definition at line 49 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::_nsamplesPN

private

Definition at line 55 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::_presample

private

Definition at line 50 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::_presamplePN

private

Definition at line 56 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::_samplemax

private

Definition at line 54 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::_samplemin

private

Definition at line 53 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

double EcalTestPulseAnalyzer::_timeofmax

private

Definition at line 61 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

double EcalTestPulseAnalyzer::adc[10]

private

Definition at line 116 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

double EcalTestPulseAnalyzer::APD[6]

private

Definition at line 142 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

double EcalTestPulseAnalyzer::apdAmpl

private

Definition at line 122 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::apdGain

private

Definition at line 119 of file EcalTestPulseAnalyzer.h.

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

double EcalTestPulseAnalyzer::apdTime

private

Definition at line 123 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

int EcalTestPulseAnalyzer::channelID

private

Definition at line 92 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::channelIteratorEE

private

Definition at line 141 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 133 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

int EcalTestPulseAnalyzer::dccID

private

Definition at line 99 of file EcalTestPulseAnalyzer.h.

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

std::vector<int> EcalTestPulseAnalyzer::dccMEM

private

Definition at line 134 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

std::string EcalTestPulseAnalyzer::digiCollection_

private

Definition at line 66 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

std::string EcalTestPulseAnalyzer::digiPNCollection_

private

Definition at line 67 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

std::string EcalTestPulseAnalyzer::digiProducer_

private

Definition at line 68 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

int EcalTestPulseAnalyzer::eta

private

Definition at line 114 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::event

private

Definition at line 115 of file EcalTestPulseAnalyzer.h.

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

std::string EcalTestPulseAnalyzer::eventHeaderCollection_

private

Definition at line 69 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

std::string EcalTestPulseAnalyzer::eventHeaderProducer_

private

Definition at line 70 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

int EcalTestPulseAnalyzer::fedID

private

Definition at line 98 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::firstChanMod[9]

private

Definition at line 148 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::flag

private

Definition at line 139 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

int EcalTestPulseAnalyzer::gain

private

Definition at line 139 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

int EcalTestPulseAnalyzer::iChannelID[1700]

private

Definition at line 146 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::idccID[1700]

private

Definition at line 146 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::ieta

private

Definition at line 139 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

int EcalTestPulseAnalyzer::iEta[1700]

private

Definition at line 144 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::iEvent

private

Definition at line 45 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

unsigned int EcalTestPulseAnalyzer::iModule[1700]

private

Definition at line 145 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::iphi

private

Definition at line 139 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

int EcalTestPulseAnalyzer::iPhi[1700]

private

Definition at line 144 of file EcalTestPulseAnalyzer.h.

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

unsigned int EcalTestPulseAnalyzer::isFirstChanModFilled[9]

private

Definition at line 149 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and beginJob().

int EcalTestPulseAnalyzer::iside[1700]

private

Definition at line 146 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::iTowerID[1700]

private

Definition at line 146 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::iZ

private

Definition at line 101 of file EcalTestPulseAnalyzer.h.

Referenced by analyze(), and EcalTestPulseAnalyzer().

int EcalTestPulseAnalyzer::moduleID

private

Definition at line 140 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

std::vector<int> EcalTestPulseAnalyzer::modules

private

Definition at line 135 of file EcalTestPulseAnalyzer.h.

unsigned int EcalTestPulseAnalyzer::nCrys

private

Definition at line 80 of file EcalTestPulseAnalyzer.h.

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

unsigned int EcalTestPulseAnalyzer::nGainAPD

private

Definition at line 84 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

unsigned int EcalTestPulseAnalyzer::nGainPN

private

Definition at line 83 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

unsigned int EcalTestPulseAnalyzer::nMod

private

Definition at line 82 of file EcalTestPulseAnalyzer.h.

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

unsigned int EcalTestPulseAnalyzer::nTT

private

Definition at line 81 of file EcalTestPulseAnalyzer.h.

Referenced by EcalTestPulseAnalyzer().

TFile* EcalTestPulseAnalyzer::outFile

private

Definition at line 105 of file EcalTestPulseAnalyzer.h.

Referenced by beginJob(), and endJob().

int EcalTestPulseAnalyzer::phi

private

Definition at line 114 of file EcalTestPulseAnalyzer.h.

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

double EcalTestPulseAnalyzer::pn[50]

private

Definition at line 117 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

double EcalTestPulseAnalyzer::PN[6]

private

Definition at line 142 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

double EcalTestPulseAnalyzer::pnAmpl

private

Definition at line 126 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

double EcalTestPulseAnalyzer::pnAmpl0

private

Definition at line 124 of file EcalTestPulseAnalyzer.h.

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

double EcalTestPulseAnalyzer::pnAmpl1

private

Definition at line 125 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::pnG

private

Definition at line 121 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

int EcalTestPulseAnalyzer::pnGain

private

Definition at line 120 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::pnID

private

Definition at line 140 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

std::string EcalTestPulseAnalyzer::resdir_

private

Definition at line 65 of file EcalTestPulseAnalyzer.h.

Referenced by beginJob(), and EcalTestPulseAnalyzer().

std::string EcalTestPulseAnalyzer::resfile

private

Definition at line 75 of file EcalTestPulseAnalyzer.h.

Referenced by beginJob(), and endJob().

TFile* EcalTestPulseAnalyzer::resFile

private

Definition at line 106 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

TTree* EcalTestPulseAnalyzer::respntrees

private

Definition at line 131 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

TTree* EcalTestPulseAnalyzer::restrees

private

Definition at line 130 of file EcalTestPulseAnalyzer.h.

Referenced by endJob().

double EcalTestPulseAnalyzer::ret_data[20]

private

Definition at line 89 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

std::string EcalTestPulseAnalyzer::rootfile

private

Definition at line 74 of file EcalTestPulseAnalyzer.h.

Referenced by beginJob(), and endJob().

int EcalTestPulseAnalyzer::runNum

private

Definition at line 97 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

int EcalTestPulseAnalyzer::runType

private

Definition at line 96 of file EcalTestPulseAnalyzer.h.

Referenced by analyze().

int EcalTestPulseAnalyzer::side

private

Definition at line 100 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::towerID

private

Definition at line 91 of file EcalTestPulseAnalyzer.h.

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

int EcalTestPulseAnalyzer::TPEvents

private

Definition at line 87 of file EcalTestPulseAnalyzer.h.

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

TTree* EcalTestPulseAnalyzer::trees[1700]

private

Definition at line 110 of file EcalTestPulseAnalyzer.h.

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