EcalABAnalyzer Class Reference

#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.

{ iBlue, iRed, nColor };

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, hgcalPerformanceValidation::df, digiCollection_, doesABTreeExist, MillePedeFileConverter_cfg::e, CollectionTags_cfi::EBDigi, ebDigiToken_, CollectionTags_cfi::EEDigi, eeDigiToken_, MEEBGeom::electronic_channel(), edm::SortedCollection< T, SORT >::end(), eta, eventHeaderCollection_, eventHeaderProducer_, fedID, spr::find(), TAPDPulse::getDelta(), ecalpyutils::hashedIndex(), mps_fire::i, iChannelID, idccID, iEta, iEvent, iPhi, iside, TAPDPulse::isPulseOK(), TAPDPulse::isTimingQualOK(), edm::HandleBase::isValid(), iTowerID, iZ, EcalDCCHeaderBlock::LASER_DELAY_SCAN, EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_POWER_SCAN, EcalDCCHeaderBlock::LASER_STD, lightside, MEEBGeom::localCoord(), mappingToken_, METSkim_cff::Max, nCrys, nevtAB, nEvtBadGain, nEvtBadTiming, nEvtTot, phi, edm::Handle< T >::product(), TShapeAnalysis::putAllVals(), rawDataToken_, runNum, runType, TAPDPulse::setPulse(), shapana, MEEBGeom::side(), MEEEGeom::side(), side, 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, l1tGTMenu_lepSeeds_cff::delta12, Delta12, doesABTreeExist, TMom::getMean(), mps_fire::i, isGainOK, isTimingOK, nCrys, nEvtBadGain, nEvtBadTiming, nEvtTot, TShapeAnalysis::set_nch(), TShapeAnalysis::set_presample(), TAPDPulse::setPresamples(), shapana, wasGainOK, and wasTimingOK.

Referenced by o2olib.O2ORunMgr::executeJob().

                            {
  //========================================================================

  edm::LogVerbatim("EcalABAnalyzer") << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
  edm::LogVerbatim("EcalABAnalyzer") << "\t+=+     Analyzing data: getting (alpha, beta)     +=+";

  // Adjust channel numbers for EE
  //===============================

  if (_ecalPart == "EE") {
    nCrys = channelMapEE.size();
    shapana->set_nch(nCrys);
  }

  // Set presamples number
  //======================

  double delta01 = Delta01->getMean();
  double delta12 = Delta12->getMean();
  if (delta12 > _presamplecut) {
    _presample = 2;
    if (delta01 > _presamplecut)
      _presample = 1;
  }

  APDPulse->setPresamples(_presample);
  shapana->set_presample(_presample);

  //  Get alpha and beta
  //======================

  if (_fitab) {
    edm::LogVerbatim("EcalABAnalyzer") << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
    edm::LogVerbatim("EcalABAnalyzer") << "\t+=+     Analyzing data: getting (alpha, beta)     +=+";
    TFile* fAB = nullptr;
    TTree* ABInit = nullptr;
    if (doesABTreeExist) {
      fAB = new TFile(alphainitfile.c_str());
    }
    if (doesABTreeExist && fAB) {
      ABInit = (TTree*)fAB->Get("ABCol0");
    }
    shapana->computeShape(alphafile, ABInit);

    // Set quality flags for gains and timing

    double BadGainEvtPercentage = 0.0;
    double BadTimingEvtPercentage = 0.0;

    int nChanBadGain = 0;
    int nChanBadTiming = 0;

    for (unsigned int i = 0; i < nCrys; i++) {
      if (nEvtTot[i] != 0) {
        BadGainEvtPercentage = double(nEvtBadGain[i]) / double(nEvtTot[i]);
        BadTimingEvtPercentage = double(nEvtBadTiming[i]) / double(nEvtTot[i]);
      }
      if (BadGainEvtPercentage > _qualpercent) {
        wasGainOK[i] = false;
        nChanBadGain++;
      }
      if (BadTimingEvtPercentage > _qualpercent) {
        wasTimingOK[i] = false;
        nChanBadTiming++;
      }
    }

    double BadGainChanPercentage = double(nChanBadGain) / double(nCrys);
    double BadTimingChanPercentage = double(nChanBadTiming) / double(nCrys);

    if (BadGainChanPercentage > _qualpercent)
      isGainOK = false;
    if (BadTimingChanPercentage > _qualpercent)
      isTimingOK = false;

    if (!isGainOK)
      edm::LogVerbatim("EcalABAnalyzer") << "\t+=+ ............................ WARNING! APD GAIN WAS NOT 1    +=+";
    if (!isTimingOK)
      edm::LogVerbatim("EcalABAnalyzer") << "\t+=+ ............................ WARNING! TIMING WAS BAD        +=+";

    edm::LogVerbatim("EcalABAnalyzer") << "\t+=+    .................................... done  +=+";
    edm::LogVerbatim("EcalABAnalyzer") << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
  }
}

Member Data Documentation

_alpha

const double EcalABAnalyzer::_alpha

private

Definition at line 79 of file EcalABAnalyzer.h.

Referenced by beginJob().

_beta

const double EcalABAnalyzer::_beta

private

Definition at line 80 of file EcalABAnalyzer.h.

Referenced by beginJob().

_chi2cut

const double EcalABAnalyzer::_chi2cut

private

Definition at line 83 of file EcalABAnalyzer.h.

Referenced by beginJob().

_debug

const int EcalABAnalyzer::_debug

private

Definition at line 87 of file EcalABAnalyzer.h.

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

_ecalPart

const std::string EcalABAnalyzer::_ecalPart

private

Definition at line 84 of file EcalABAnalyzer.h.

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

_fedid

const int EcalABAnalyzer::_fedid

private

Definition at line 85 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

_firstsample

const unsigned int EcalABAnalyzer::_firstsample

private

Definition at line 68 of file EcalABAnalyzer.h.

Referenced by beginJob(), and EcalABAnalyzer().

_fitab

bool EcalABAnalyzer::_fitab

private

Definition at line 109 of file EcalABAnalyzer.h.

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

_lastsample

const unsigned int EcalABAnalyzer::_lastsample

private

Definition at line 69 of file EcalABAnalyzer.h.

Referenced by beginJob(), and EcalABAnalyzer().

_nevtmax

const unsigned int EcalABAnalyzer::_nevtmax

private

Definition at line 81 of file EcalABAnalyzer.h.

Referenced by analyze(), and beginJob().

_niter

const unsigned int EcalABAnalyzer::_niter

private

Definition at line 78 of file EcalABAnalyzer.h.

_noise

const double EcalABAnalyzer::_noise

private

Definition at line 82 of file EcalABAnalyzer.h.

Referenced by beginJob().

_nsamples

const unsigned int EcalABAnalyzer::_nsamples

private

Definition at line 66 of file EcalABAnalyzer.h.

Referenced by beginJob(), and EcalABAnalyzer().

_presample

unsigned int EcalABAnalyzer::_presample

private

Definition at line 67 of file EcalABAnalyzer.h.

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

_presamplecut

const double EcalABAnalyzer::_presamplecut

private

Definition at line 77 of file EcalABAnalyzer.h.

Referenced by endJob().

_qualpercent

const double EcalABAnalyzer::_qualpercent

private

Definition at line 86 of file EcalABAnalyzer.h.

Referenced by endJob().

_ratiomaxcutlow

const double EcalABAnalyzer::_ratiomaxcutlow

private

Definition at line 76 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

_ratiomincuthigh

const double EcalABAnalyzer::_ratiomincuthigh

private

Definition at line 75 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

_ratiomincutlow

const double EcalABAnalyzer::_ratiomincutlow

private

Definition at line 74 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

_timingcuthigh

const unsigned int EcalABAnalyzer::_timingcuthigh

private

Definition at line 71 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

_timingcutlow

const unsigned int EcalABAnalyzer::_timingcutlow

private

Definition at line 70 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

_timingqualhigh

const unsigned int EcalABAnalyzer::_timingqualhigh

private

Definition at line 73 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

_timingquallow

const unsigned int EcalABAnalyzer::_timingquallow

private

Definition at line 72 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

adc

double EcalABAnalyzer::adc[10]

private

Definition at line 132 of file EcalABAnalyzer.h.

Referenced by analyze().

adcG

int EcalABAnalyzer::adcG[10]

private

Definition at line 133 of file EcalABAnalyzer.h.

Referenced by analyze().

alphafile

std::string EcalABAnalyzer::alphafile

private

Definition at line 97 of file EcalABAnalyzer.h.

Referenced by beginJob(), and endJob().

alphainitfile

std::string EcalABAnalyzer::alphainitfile

private

Definition at line 98 of file EcalABAnalyzer.h.

Referenced by beginJob(), and endJob().

APDPulse

TAPDPulse* EcalABAnalyzer::APDPulse

private

Definition at line 89 of file EcalABAnalyzer.h.

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

channelIteratorEE

int EcalABAnalyzer::channelIteratorEE

private

Definition at line 134 of file EcalABAnalyzer.h.

Referenced by analyze().

channelMapEE

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

private

Definition at line 123 of file EcalABAnalyzer.h.

Referenced by analyze(), and endJob().

color

int EcalABAnalyzer::color

private

Definition at line 131 of file EcalABAnalyzer.h.

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

colors

std::vector<int> EcalABAnalyzer::colors

private

Definition at line 122 of file EcalABAnalyzer.h.

dccID

int EcalABAnalyzer::dccID

private

Definition at line 115 of file EcalABAnalyzer.h.

Referenced by analyze().

dccMEM

std::vector<int> EcalABAnalyzer::dccMEM

private

Definition at line 124 of file EcalABAnalyzer.h.

Delta01

TMom* EcalABAnalyzer::Delta01

private

Definition at line 90 of file EcalABAnalyzer.h.

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

Delta12

TMom* EcalABAnalyzer::Delta12

private

Definition at line 91 of file EcalABAnalyzer.h.

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

digiCollection_

const std::string EcalABAnalyzer::digiCollection_

private

Definition at line 56 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

digiProducer_

const std::string EcalABAnalyzer::digiProducer_

private

Definition at line 57 of file EcalABAnalyzer.h.

Referenced by EcalABAnalyzer().

doesABTreeExist

bool EcalABAnalyzer::doesABTreeExist

private

Definition at line 107 of file EcalABAnalyzer.h.

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

ebDigiToken_

edm::EDGetTokenT<EBDigiCollection> EcalABAnalyzer::ebDigiToken_

private

Definition at line 60 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

eeDigiToken_

edm::EDGetTokenT<EEDigiCollection> EcalABAnalyzer::eeDigiToken_

private

Definition at line 61 of file EcalABAnalyzer.h.

Referenced by analyze(), and EcalABAnalyzer().

eta

int EcalABAnalyzer::eta

private

Definition at line 129 of file EcalABAnalyzer.h.

Referenced by Particle.Particle::__str__(), analyze(), 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().