EcalLaserAnalyzer2 Class Reference

#include <EcalLaserAnalyzer2.h>

Inheritance diagram for EcalLaserAnalyzer2:

Public Types
enum	VarCol { iBlue, iRed, nColor }
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	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
	EcalLaserAnalyzer2 (const edm::ParameterSet &iConfig)
void	endJob () override
void	setGeomEB (int etaG, int phiG, int module, int tower, int strip, int xtal, int apdRefTT, int channel, int lmr)
void	setGeomEE (int etaG, int phiG, int iX, int iY, int iZ, int module, int tower, int ch, int apdRefTT, int channel, int lmr)
	~EcalLaserAnalyzer2 () 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 Member Functions
bool	getShapes ()

Private Attributes
int	_debug
bool	_docorpn
std::string	_ecalPart
int	_fedid
unsigned int	_firstsample
unsigned int	_firstsamplePN
unsigned int	_lastsample
unsigned int	_lastsamplePN
unsigned int	_niter
double	_noise
unsigned int	_nsamples
unsigned int	_nsamplesPN
unsigned int	_presample
double	_presamplecut
unsigned int	_presamplePN
double	_qualpercent
double	_ratiomaxcutlow
double	_ratiomincuthigh
double	_ratiomincutlow
unsigned int	_samplemax
unsigned int	_samplemin
bool	_saveallevents
bool	_saveshapes
unsigned int	_timingcuthigh
unsigned int	_timingcutlow
unsigned int	_timingqualhigh
unsigned int	_timingquallow
double	adc [10]
int	adcC [1700]
std::string	ADCfile
TFile *	ADCFile
int	adcG [10]
double	adcMean [1700][10]
double *	adcNoPed
TTree *	ADCtrees [1700]
double	APD [6]
double	apdAmpl
double	apdAmplA
double	apdAmplB
TAPD *	APDAnal [1700][nColor]
std::string	APDfile
TFile *	APDFile
TAPD *	APDFirstAnal [1700][nColor]
double	APDoAPDA [6]
double	APDoAPDB [6]
double	APDoPN [6]
double	APDoPNA [6]
double	APDoPNB [6]
TAPDPulse *	APDPulse
std::map< int, unsigned int >	apdRefMap [2]
double	apdTime
TTree *	APDtrees [1700]
int	channelID
int	channelIteratorEE
std::map< int, int >	channelMapEE
int	color
int	colorref
std::vector< int >	colors
int	dccID
std::vector< int >	dccMEM
TMom *	Delta01
TMom *	Delta12
std::string	digiCollection_
std::string	digiPNCollection_
std::string	digiProducer_
std::string	elecfile_
int	eta
int	event
std::string	eventHeaderCollection_
std::string	eventHeaderProducer_
int	eventref
int	fedID
unsigned int	firstChanMod [21]
int	flag
int	iChannelID [1700]
int	idccID [1700]
int	iEta [1700]
int	iEvent
unsigned int	iModule [1700]
int	iPhi [1700]
unsigned int	isFirstChanModFilled [21]
bool	isGainOK
int	iside [1700]
bool	IsMatacqOK
bool	isMatacqOK
bool	isSPRFine
int	IsThereDataADC [1700][nColor]
bool	isTimingOK
int	iTowerID [1700]
int	iZ
int	laserEvents
int	lightside
std::string	matfile
TMem *	Mem
int	moduleID
std::vector< int >	modules
unsigned int	nCrys
int	nEvtBadGain [1700]
int	nEvtBadTiming [1700]
int	nEvtTot [1700]
unsigned int	nMod
unsigned int	nPNPerMod
unsigned int	nRefChan
unsigned int	nRefTrees
unsigned int	nSamplesShapes
unsigned int	nSides
int	phi
double	pn [50]
double	PN [6]
double	pn0
double	pn1
double	pnAmpl
TPN *	PNAnal [9][2][nColor]
std::string	pncorfile_
TPNCor *	pnCorrector
TPN *	PNFirstAnal [9][2][nColor]
int	pnG [50]
int	pnID
double *	pnNoPed
double	PNoPN [6]
double	PNoPNA [6]
double	PNoPNB [6]
TPNPulse *	PNPulse
TProfile *	PulseShape
TTree *	RefAPDtrees [2][21]
std::string	resdir_
std::string	resfile
TFile *	resFile
TTree *	respntrees [nColor]
TTree *	restrees [nColor]
int	runNum
int	runType
double	ShapeCor
double	shapeCorrection
std::string	shapefile
TFile *	ShapeFile
double	shapes [250]
std::vector< double >	shapesVec
int	side
double	Time [6]
int	towerID
bool	wasGainOK [1700]
bool	wasTimingOK [1700]

Additional Inherited Members
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 43 of file EcalLaserAnalyzer2.h.

Member Enumeration Documentation

enum EcalLaserAnalyzer2::VarCol

Enumerator
iBlue
iRed
nColor

Definition at line 56 of file EcalLaserAnalyzer2.h.

{ iBlue, iRed, nColor };

Constructor & Destructor Documentation

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

explicit

Definition at line 56 of file EcalLaserAnalyzer2.cc.

References _ecalPart, _fedid, _firstsample, _lastsample, _nsamples, _nsamplesPN, _presample, _presamplePN, _ratiomaxcutlow, _ratiomincuthigh, _ratiomincutlow, _timingcuthigh, _timingcutlow, _timingqualhigh, _timingquallow, APDPulse, Delta01, Delta12, digiCollection_, digiPNCollection_, digiProducer_, elecfile_, eventHeaderCollection_, eventHeaderProducer_, firstChanMod, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), iChannelID, idccID, iEta, cuy::ii, iModule, iPhi, isFirstChanModFilled, isGainOK, iside, isTimingOK, iTowerID, iZ, dqmiolumiharvest::j, ME::lmmodFromDcc(), Mem, nCrys, NCRYSEB, NCRYSEE, nMod, NREFCHAN, nRefChan, pncorfile_, pnCorrector, PNPulse, resdir_, AlCaHLTBitMon_QueryRunRegistry::string, wasGainOK, and wasTimingOK.

    : iEvent(0),

      // 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)),
      _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)),
      _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", 5)),
      _noise(iConfig.getUntrackedParameter<double>("noise", 2.0)),
      _ecalPart(iConfig.getUntrackedParameter<std::string>("ecalPart", "EB")),
      _saveshapes(iConfig.getUntrackedParameter<bool>("saveShapes", true)),
      _docorpn(iConfig.getUntrackedParameter<bool>("doCorPN", false)),
      _fedid(iConfig.getUntrackedParameter<int>("fedID", -999)),
      _saveallevents(iConfig.getUntrackedParameter<bool>("saveAllEvents", false)),
      _qualpercent(iConfig.getUntrackedParameter<double>("qualPercent", 0.2)),
      _debug(iConfig.getUntrackedParameter<int>("debug", 0)),
      nCrys(NCRYSEB),
      nPNPerMod(NPNPERMOD),
      nMod(NMODEB),
      nSides(NSIDES),
      nSamplesShapes(NSAMPSHAPES),
      IsMatacqOK(false),
      runType(-1),
      runNum(0),
      towerID(-1),
      channelID(-1),
      fedID(-1),
      dccID(-1),
      side(2),
      lightside(2),
      iZ(1),
      phi(-1),
      eta(-1),
      event(0),
      color(0),
      pn0(0),
      pn1(0),
      apdAmpl(0),
      apdTime(0),
      pnAmpl(0),
      pnID(-1),
      moduleID(-1),
      flag(0),
      channelIteratorEE(0),
      ShapeCor(0)

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

{
  // Initialization from cfg file

  resdir_ = iConfig.getUntrackedParameter<std::string>("resDir");
  elecfile_ = iConfig.getUntrackedParameter<std::string>("elecFile");
  pncorfile_ = iConfig.getUntrackedParameter<std::string>("pnCorFile");

  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");

  // Geometrical constants initialization
  if (_ecalPart == "EB") {
    nCrys = NCRYSEB;
  } else {
    nCrys = NCRYSEE;
  }
  iZ = 1;
  if (_fedid <= 609)
    iZ = -1;
  modules = ME::lmmodFromDcc(_fedid);
  nMod = modules.size();
  nRefChan = NREFCHAN;

  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;
    wasTimingOK[j] = true;
    wasGainOK[j] = true;
  }

  for (unsigned int j = 0; j < nMod; j++) {
    int ii = modules[j];
    firstChanMod[ii - 1] = 0;
    isFirstChanModFilled[ii - 1] = 0;
  }

  // Quality check flags

  isGainOK = true;
  isTimingOK = true;

  // PN linearity corrector

  pnCorrector = new TPNCor(pncorfile_);

  // Objects dealing with pulses

  APDPulse = new TAPDPulse(_nsamples,
                           _presample,
                           _firstsample,
                           _lastsample,
                           _timingcutlow,
                           _timingcuthigh,
                           _timingquallow,
                           _timingqualhigh,
                           _ratiomincutlow,
                           _ratiomincuthigh,
                           _ratiomaxcutlow);
  PNPulse = new TPNPulse(_nsamplesPN, _presamplePN);

  // Object dealing with MEM numbering

  Mem = new TMem(_fedid);

  // Objects needed for npresample calculation

  Delta01 = new TMom();
  Delta12 = new TMom();
}

EcalLaserAnalyzer2::~EcalLaserAnalyzer2 ( )

override

Definition at line 198 of file EcalLaserAnalyzer2.cc.

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

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

Member Function Documentation

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

override

Definition at line 280 of file EcalLaserAnalyzer2.cc.

References _debug, _docorpn, _ecalPart, _fedid, _firstsamplePN, _lastsamplePN, _nsamplesPN, adc, adcG, ADCtrees, TMom::addEntry(), APDPulse, MEEBGeom::apdRefTower(), MEEEGeom::apdRefTower(), edm::SortedCollection< T, SORT >::begin(), edm::DataFrameContainer::begin(), beam_dqm_sourceclient-live_cfg::cerr, EcalElectronicsId::channelId(), channelID, channelIteratorEE, channelMapEE, color, corr, gather_cfg::cout, dccID, MEEEGeom::dee(), Delta01, Delta12, 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, spr::find(), edm::EventSetup::get(), TPNPulse::getAdcWithoutPedestal(), TPNFit::getAmpl(), edm::Event::getByLabel(), TAPDPulse::getDelta(), EcalElectronicsMapping::getElectronicsId(), TPNPulse::getMaxSample(), TPNCor::getPNCorrectionFactor(), ecalpyutils::hashedIndex(), mps_fire::i, EcalPnDiodeDetId::iDCCId(), iEvent, TPNFit::init(), createfilelist::int, EcalPnDiodeDetId::iPnId(), IsMatacqOK, TMem::isMemRelevant(), TAPDPulse::isPulseOK(), TAPDPulse::isTimingQualOK(), iZ, EcalDCCHeaderBlock::LASER_DELAY_SCAN, EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_POWER_SCAN, EcalDCCHeaderBlock::LASER_STD, laserEvents, lightside, MEEEGeom::lmmod(), MEEBGeom::lmmod(), MEEBGeom::lmr(), MEEEGeom::lmr(), MEEBGeom::localCoord(), Max(), TMem::Mem(), Mem, nCrys, nEvtBadGain, nEvtBadTiming, nEvtTot, phi, MEEBGeom::pn(), MEEEGeom::pn(), pn, pn0, pn1, pnAmpl, pnCorrector, pnG, pnNoPed, PNPulse, edm::Handle< T >::product(), edm::ESHandle< T >::product(), runNum, runType, setGeomEB(), setGeomEE(), TPNPulse::setPulse(), TAPDPulse::setPulse(), side, findQualityFiles::size, digitizers_cfi::strip, EcalElectronicsId::stripId(), hgcalTowerProducer_cfi::tower, EcalElectronicsId::towerId(), towerID, EcalDCCHeaderBlock::EcalDCCEventSettings::wavelength, and EcalElectronicsId::xtalId().

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

  ++iEvent;

  // retrieving 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 data from Event
  edm::Handle<EBDigiCollection> pEBDigi;
  const EBDigiCollection* EBDigi = nullptr;
  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 if (_ecalPart == "EE") {
    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;
    }
  } else {
    cout << " Wrong ecalPart in cfg file " << endl;
    return;
  }

  // retrieving 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 " << digiPNCollection_.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) {
    // 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;

    vector<int>::iterator iter = find(colors.begin(), colors.end(), color);
    if (iter == colors.end()) {
      colors.push_back(color);
      cout << " new color found " << color << " " << colors.size() << endl;
    }
  }

  // Check Matacq shape exists

  if (!IsMatacqOK)
    return;

  // Cut on fedID

  if (fedID != _fedid && _fedid != -999)
    return;

  // Count laser events

  laserEvents++;

  // ======================
  // Decode PN Information
  // ======================

  TPNFit* pnfit = new TPNFit();
  pnfit->init(_nsamplesPN, _firstsamplePN, _lastsamplePN);

  double chi2pn = 0;
  unsigned int samplemax = 0;
  int pnGain = 0;

  map<int, vector<double> > allPNAmpl;
  map<int, vector<double> > allPNGain;

  // Loop on PNs digis

  for (EcalPnDiodeDigiCollection::const_iterator pnItr = PNDigi->begin(); pnItr != PNDigi->end(); ++pnItr) {
    EcalPnDiodeDetId pnDetId = EcalPnDiodeDetId((*pnItr).id());

    if (_debug == 1)
      cout << "-- debug test -- Inside PNDigi - pnID=" << pnDetId.iPnId() << ", dccID=" << pnDetId.iDCCId() << endl;

    // Skip MEM DCC without relevant data

    bool isMemRelevant = Mem->isMemRelevant(pnDetId.iDCCId());
    if (!isMemRelevant)
      continue;

    // Loop on PN samples

    for (int samId = 0; samId < (*pnItr).size(); samId++) {
      pn[samId] = (*pnItr).sample(samId).adc();
      pnG[samId] = (*pnItr).sample(samId).gainId();
      if (samId == 0)
        pnGain = pnG[samId];
      if (samId > 0)
        pnGain = int(TMath::Max(pnG[samId], pnGain));
    }

    if (pnGain != 1)
      cout << "PN gain different from 1" << endl;

    // Calculate amplitude from pulse

    PNPulse->setPulse(pn);
    pnNoPed = PNPulse->getAdcWithoutPedestal();
    samplemax = PNPulse->getMaxSample();
    chi2pn = pnfit->doFit(samplemax, &pnNoPed[0]);
    if (chi2pn == 101 || chi2pn == 102 || chi2pn == 103)
      pnAmpl = 0.;
    else
      pnAmpl = pnfit->getAmpl();

    // Apply linearity correction

    double corr = 1.0;
    if (_docorpn)
      corr = pnCorrector->getPNCorrectionFactor(pnAmpl, pnGain);
    pnAmpl *= corr;

    // Fill PN ampl vector

    allPNAmpl[pnDetId.iDCCId()].push_back(pnAmpl);

    if (_debug == 1)
      cout << "-- debug -- Inside PNDigi - PNampl=" << pnAmpl << ", PNgain=" << pnGain << endl;
  }

  // ===========================
  // 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);
      EcalElectronicsId elecid_crystal = TheMapping->getElectronicsId(id_crystal);

      int etaG = id_crystal.ieta();  // global
      int phiG = id_crystal.iphi();  // global

      std::pair<int, int> LocalCoord = MEEBGeom::localCoord(etaG, phiG);

      int etaL = LocalCoord.first;   // local
      int phiL = LocalCoord.second;  // local

      int strip = elecid_crystal.stripId();
      int xtal = elecid_crystal.xtalId();

      int module = MEEBGeom::lmmod(etaG, phiG);
      int tower = elecid_crystal.towerId();

      int apdRefTT = MEEBGeom::apdRefTower(module);

      std::pair<int, int> pnpair = MEEBGeom::pn(module);
      unsigned int MyPn0 = pnpair.first;
      unsigned int MyPn1 = pnpair.second;

      int lmr = MEEBGeom::lmr(etaG, phiG);
      unsigned int channel = MEEBGeom::electronic_channel(etaL, phiL);
      assert(channel < nCrys);

      setGeomEB(etaG, phiG, module, tower, strip, xtal, apdRefTT, channel, lmr);

      if (_debug == 1)
        cout << "-- debug -- Inside EBDigi - towerID:" << towerID << " channelID:" << channelID << " module:" << module
             << " modules:" << modules.size() << endl;

      // 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]++;

      // Associate PN ampl
      //===================

      int mem0 = Mem->Mem(lmr, 0);
      int mem1 = Mem->Mem(lmr, 1);

      if (allPNAmpl[mem0].size() > MyPn0)
        pn0 = allPNAmpl[mem0][MyPn0];
      else
        pn0 = 0;
      if (allPNAmpl[mem1].size() > MyPn1)
        pn1 = allPNAmpl[mem1][MyPn1];
      else
        pn1 = 0;

      // Fill if Pulse is fine
      //=======================

      if (APDPulse->isPulseOK() && lightside == side) {
        ADCtrees[channel]->Fill();

        Delta01->addEntry(APDPulse->getDelta(0, 1));
        Delta12->addEntry(APDPulse->getDelta(1, 2));
      }
    }
  } else if (EEDigi) {
    // Loop on crystals
    //===================

    for (EEDigiCollection::const_iterator digiItr = EEDigi->begin(); digiItr != EEDigi->end(); ++digiItr) {
      // Retrieve geometry
      //===================

      EEDetId id_crystal(digiItr->id());
      EEDataFrame df(*digiItr);
      EcalElectronicsId elecid_crystal = TheMapping->getElectronicsId(id_crystal);

      int etaG = id_crystal.iy();
      int phiG = id_crystal.ix();

      int iX = (phiG - 1) / 5 + 1;
      int iY = (etaG - 1) / 5 + 1;

      int tower = elecid_crystal.towerId();
      int ch = elecid_crystal.channelId() - 1;

      int module = MEEEGeom::lmmod(iX, iY);
      if (module >= 18 && side == 1)
        module += 2;
      int lmr = MEEEGeom::lmr(iX, iY, iZ);
      int dee = MEEEGeom::dee(lmr);
      int apdRefTT = MEEEGeom::apdRefTower(lmr, module);

      std::pair<int, int> pnpair = MEEEGeom::pn(dee, module);
      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];
      assert(channel < nCrys);

      setGeomEE(etaG, phiG, iX, iY, iZ, module, tower, ch, apdRefTT, channel, lmr);

      if (_debug == 1)
        cout << "-- debug -- Inside EEDigi - towerID:" << towerID << " channelID:" << channelID << " module:" << module
             << " modules:" << modules.size() << endl;

      // APD Pulse
      //===========

      if ((*digiItr).size() > 10)
        cout << "SAMPLES SIZE > 10!" << (*digiItr).size() << endl;

      // 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]++;

      // Associate PN ampl
      //===================

      int mem0 = Mem->Mem(lmr, 0);
      int mem1 = Mem->Mem(lmr, 1);

      if (allPNAmpl[mem0].size() > MyPn0)
        pn0 = allPNAmpl[mem0][MyPn0];
      else
        pn0 = 0;
      if (allPNAmpl[mem1].size() > MyPn1)
        pn1 = allPNAmpl[mem1][MyPn1];
      else
        pn1 = 0;

      // Fill if Pulse is fine
      //=======================

      if (APDPulse->isPulseOK() && lightside == side) {
        ADCtrees[channel]->Fill();

        Delta01->addEntry(APDPulse->getDelta(0, 1));
        Delta12->addEntry(APDPulse->getDelta(1, 2));
      }
    }
  }
}

void EcalLaserAnalyzer2::beginJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 206 of file EcalLaserAnalyzer2.cc.

References adc, ADCfile, ADCFile, ADCtrees, APDfile, channelID, color, gather_cfg::cout, dccID, eta, getShapes(), mps_fire::i, IsMatacqOK, laserEvents, matfile, Skims_PA_cff::name, nCrys, phi, pn0, pn1, resdir_, resfile, shapefile, side, and towerID.

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

  // Create temporary files and trees to save adc samples
  //======================================================

  ADCfile = resdir_;
  ADCfile += "/APDSamplesLaser.root";

  APDfile = resdir_;
  APDfile += "/APDPNLaserAllEvents.root";

  ADCFile = new TFile(ADCfile.c_str(), "RECREATE");

  for (unsigned int i = 0; i < nCrys; i++) {
    stringstream name;
    name << "ADCTree" << i + 1;
    ADCtrees[i] = new TTree(name.str().c_str(), name.str().c_str());

    ADCtrees[i]->Branch("iphi", &phi, "phi/I");
    ADCtrees[i]->Branch("ieta", &eta, "eta/I");
    ADCtrees[i]->Branch("towerID", &towerID, "towerID/I");
    ADCtrees[i]->Branch("channelID", &channelID, "channelID/I");
    ADCtrees[i]->Branch("dccID", &dccID, "dccID/I");
    ADCtrees[i]->Branch("side", &side, "side/I");
    ADCtrees[i]->Branch("event", &event, "event/I");
    ADCtrees[i]->Branch("color", &color, "color/I");
    ADCtrees[i]->Branch("adc", &adc, "adc[10]/D");
    ADCtrees[i]->Branch("pn0", &pn0, "pn0/D");
    ADCtrees[i]->Branch("pn1", &pn1, "pn1/D");

    ADCtrees[i]->SetBranchAddress("ieta", &eta);
    ADCtrees[i]->SetBranchAddress("iphi", &phi);
    ADCtrees[i]->SetBranchAddress("towerID", &towerID);
    ADCtrees[i]->SetBranchAddress("channelID", &channelID);
    ADCtrees[i]->SetBranchAddress("dccID", &dccID);
    ADCtrees[i]->SetBranchAddress("side", &side);
    ADCtrees[i]->SetBranchAddress("event", &event);
    ADCtrees[i]->SetBranchAddress("color", &color);
    ADCtrees[i]->SetBranchAddress("adc", adc);
    ADCtrees[i]->SetBranchAddress("pn0", &pn0);
    ADCtrees[i]->SetBranchAddress("pn1", &pn1);
  }

  // Define output results filenames
  //==================================
  stringstream namefile1;
  namefile1 << resdir_ << "/SHAPE_LASER.root";
  shapefile = namefile1.str();

  stringstream namefile2;
  namefile2 << resdir_ << "/APDPN_LASER.root";
  resfile = namefile2.str();

  stringstream namefile3;
  namefile3 << resdir_ << "/MATACQ.root";

  matfile = namefile3.str();

  // Get Pulse Shapes
  //==================

  IsMatacqOK = getShapes();
  if (!IsMatacqOK) {
    cout << " ERROR! No matacq shape available: analysis aborted !" << endl;
    return;
  }

  // Laser events counter

  laserEvents = 0;
}

void EcalLaserAnalyzer2::endJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 677 of file EcalLaserAnalyzer2.cc.

References _debug, _ecalPart, _firstsample, _lastsample, _niter, _noise, _nsamples, _presample, _presamplecut, _qualpercent, _saveallevents, adc, ADCfile, ADCFile, adcNoPed, ADCtrees, TAPD::addEntry(), TPN::addEntry(), APD, apdAmpl, apdAmplA, apdAmplB, APDAnal, APDfile, APDFile, APDFirstAnal, APDoAPDA, APDoAPDB, APDoPN, APDoPNA, APDoPNB, APDPulse, apdRefMap, apdTime, APDtrees, channelID, channelMapEE, color, colorref, gather_cfg::cout, TkAlMuonSelectors_cfi::cut, dccID, Delta01, Delta12, PulseFitWithShape::doFit(), eta, event, eventref, firstChanMod, flag, TAPDPulse::getAdcWithoutPedestal(), PulseFitWithShape::getAmpl(), TAPD::getAPD(), TAPD::getAPDoAPD0(), TAPD::getAPDoAPD1(), TAPD::getAPDoPN(), TAPD::getAPDoPN0(), TAPD::getAPDoPN1(), TMom::getMean(), TPN::getPN(), TPN::getPNoPN(), TPN::getPNoPN0(), TPN::getPNoPN1(), PulseFitWithShape::getTime(), TAPD::getTime(), edmStreamStallGrapher::getTime(), mps_fire::i, iChannelID, idccID, LEDCalibrationChannels::ieta, iEta, iModule, PulseFitWithShape::init(), LEDCalibrationChannels::iphi, iPhi, isGainOK, iside, IsMatacqOK, TAPDPulse::isPulseOK(), isSPRFine, IsThereDataADC, isTimingOK, iTowerID, laserEvents, moduleID, Skims_PA_cff::name, nColor, nCrys, nEvtBadGain, nEvtBadTiming, nEvtTot, nMod, nPNPerMod, nRefChan, nSamplesShapes, phi, PN, pn0, pn1, PNAnal, PNFirstAnal, pnID, PNoPN, PNoPNA, PNoPNB, RefAPDtrees, resfile, resFile, respntrees, restrees, TAPD::set2DAPDoAPD0Cut(), TAPD::set2DAPDoAPD1Cut(), TAPD::setAPDCut(), TAPD::setAPDoPN0Cut(), TAPD::setAPDoPN1Cut(), TAPD::setAPDoPNCut(), TPN::setPNCut(), TAPDPulse::setPresamples(), TAPDPulse::setPulse(), TAPD::setTimeCut(), ShapeCor, shapeCorrection, shapesVec, side, Time, towerID, wasGainOK, and wasTimingOK.

Referenced by o2olib.O2ORunMgr::executeJob().

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

  if (!IsMatacqOK) {
    cout << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+" << endl;
    cout << "\t+=+     WARNING! NO MATACQ        +=+" << endl;
    cout << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+" << endl;
    return;
  }

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

  if (_ecalPart == "EE") {
    nCrys = channelMapEE.size();
  }

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

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

  APDPulse->setPresamples(_presample);

  // Don't do anything if there is no events
  //=========================================

  if (laserEvents == 0) {
    ADCFile->Close();
    stringstream del;
    del << "rm " << ADCfile;
    system(del.str().c_str());
    cout << " No Laser Events " << endl;
    return;
  }

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

  // Analyze adc samples to get amplitudes
  //=======================================

  cout << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+" << endl;
  cout << "\t+=+     Analyzing laser data: getting APD, PN, APD/PN, PN/PN    +=+" << endl;

  if (!isGainOK)
    cout << "\t+=+ ............................ WARNING! APD GAIN WAS NOT 1    +=+" << endl;
  if (!isTimingOK)
    cout << "\t+=+ ............................ WARNING! TIMING WAS BAD        +=+" << endl;

  APDFile = new TFile(APDfile.c_str(), "RECREATE");

  int ieta, iphi;

  int flagfit;
  for (unsigned int i = 0; i < nCrys; i++) {
    stringstream name;
    name << "APDTree" << i + 1;

    APDtrees[i] = new TTree(name.str().c_str(), name.str().c_str());

    //List of branches

    APDtrees[i]->Branch("event", &event, "event/I");
    APDtrees[i]->Branch("color", &color, "color/I");
    APDtrees[i]->Branch("iphi", &iphi, "iphi/I");
    APDtrees[i]->Branch("ieta", &ieta, "ieta/I");
    APDtrees[i]->Branch("side", &side, "side/I");
    APDtrees[i]->Branch("dccID", &dccID, "dccID/I");
    APDtrees[i]->Branch("towerID", &towerID, "towerID/I");
    APDtrees[i]->Branch("channelID", &channelID, "channelID/I");
    APDtrees[i]->Branch("apdAmpl", &apdAmpl, "apdAmpl/D");
    APDtrees[i]->Branch("apdTime", &apdTime, "apdTime/D");
    if (_saveallevents)
      APDtrees[i]->Branch("adc", &adc, "adc[10]/D");
    APDtrees[i]->Branch("flagfit", &flagfit, "flagfit/I");
    APDtrees[i]->Branch("pn0", &pn0, "pn0/D");
    APDtrees[i]->Branch("pn1", &pn1, "pn1/D");

    APDtrees[i]->SetBranchAddress("event", &event);
    APDtrees[i]->SetBranchAddress("color", &color);
    APDtrees[i]->SetBranchAddress("iphi", &iphi);
    APDtrees[i]->SetBranchAddress("ieta", &ieta);
    APDtrees[i]->SetBranchAddress("side", &side);
    APDtrees[i]->SetBranchAddress("dccID", &dccID);
    APDtrees[i]->SetBranchAddress("towerID", &towerID);
    APDtrees[i]->SetBranchAddress("channelID", &channelID);
    APDtrees[i]->SetBranchAddress("apdAmpl", &apdAmpl);
    APDtrees[i]->SetBranchAddress("apdTime", &apdTime);
    if (_saveallevents)
      APDtrees[i]->SetBranchAddress("adc", adc);
    APDtrees[i]->SetBranchAddress("flagfit", &flagfit);
    APDtrees[i]->SetBranchAddress("pn0", &pn0);
    APDtrees[i]->SetBranchAddress("pn1", &pn1);
  }

  for (unsigned int iref = 0; iref < nRefChan; iref++) {
    for (unsigned int imod = 0; imod < nMod; imod++) {
      int jmod = modules[imod];

      stringstream nameref;
      nameref << "refAPDTree" << imod << "_" << iref;

      RefAPDtrees[iref][jmod] = new TTree(nameref.str().c_str(), nameref.str().c_str());

      RefAPDtrees[iref][jmod]->Branch("eventref", &eventref, "eventref/I");
      RefAPDtrees[iref][jmod]->Branch("colorref", &colorref, "colorref/I");
      if (iref == 0)
        RefAPDtrees[iref][jmod]->Branch("apdAmplA", &apdAmplA, "apdAmplA/D");
      if (iref == 1)
        RefAPDtrees[iref][jmod]->Branch("apdAmplB", &apdAmplB, "apdAmplB/D");

      RefAPDtrees[iref][jmod]->SetBranchAddress("eventref", &eventref);
      RefAPDtrees[iref][jmod]->SetBranchAddress("colorref", &colorref);
      if (iref == 0)
        RefAPDtrees[iref][jmod]->SetBranchAddress("apdAmplA", &apdAmplA);
      if (iref == 1)
        RefAPDtrees[iref][jmod]->SetBranchAddress("apdAmplB", &apdAmplB);
    }
  }

  assert(colors.size() <= nColor);
  unsigned int nCol = colors.size();

  // Declare PN stuff
  //===================

  for (unsigned int iM = 0; iM < nMod; iM++) {
    unsigned int iMod = modules[iM] - 1;
    for (unsigned int ich = 0; ich < nPNPerMod; ich++) {
      for (unsigned int icol = 0; icol < nCol; icol++) {
        PNFirstAnal[iMod][ich][icol] = new TPN(ich);
        PNAnal[iMod][ich][icol] = new TPN(ich);
      }
    }
  }

  // Declare function for APD ampl fit
  //===================================

  PulseFitWithShape* psfit = new PulseFitWithShape();

  for (unsigned int iCry = 0; iCry < nCrys; iCry++) {
    for (unsigned int icol = 0; icol < nCol; icol++) {
      // Declare APD stuff
      //===================

      APDFirstAnal[iCry][icol] = new TAPD();
      IsThereDataADC[iCry][icol] = 1;
      stringstream cut;
      cut << "color==" << colors.at(icol);
      if (ADCtrees[iCry]->GetEntries(cut.str().c_str()) < 10)
        IsThereDataADC[iCry][icol] = 0;
    }

    unsigned int iMod = iModule[iCry] - 1;
    assert(iMod <= nMod);

    if (isSPRFine)
      psfit->init(_nsamples, _firstsample, _lastsample, _niter, nSamplesShapes, shapesVec, _noise);

    // Loop on events
    //================

    Long64_t nbytes = 0, nb = 0;
    for (Long64_t jentry = 0; jentry < ADCtrees[iCry]->GetEntriesFast(); jentry++) {  // Loop on events
      nb = ADCtrees[iCry]->GetEntry(jentry);
      nbytes += nb;

      flagfit = 1;
      apdAmpl = 0.0;
      apdTime = 0.0;
      ieta = eta;
      iphi = phi;

      // Get back color

      unsigned int iCol = 0;
      for (unsigned int i = 0; i < nCol; i++) {
        if (color == colors[i]) {
          iCol = i;
          i = colors.size();
        }
      }

      // Amplitude calculation

      APDPulse->setPulse(adc);
      adcNoPed = APDPulse->getAdcWithoutPedestal();

      if (isSPRFine && APDPulse->isPulseOK()) {
        psfit->doFit(&adcNoPed[0]);
        apdAmpl = psfit->getAmpl();
        apdTime = psfit->getTime();

      } else {
        apdAmpl = 0;
        apdTime = 0;
        flagfit = 0;
      }

      if (_debug >= 1)
        cout << "-- debug test -- endJob -- apdAmpl:" << apdAmpl << " apdTime:" << apdTime << endl;

      double pnmean;
      if (pn0 < 10 && pn1 > 10) {
        pnmean = pn1;
      } else if (pn1 < 10 && pn0 > 10) {
        pnmean = pn0;
      } else
        pnmean = 0.5 * (pn0 + pn1);

      if (_debug >= 1)
        cout << "-- debug test -- endJob -- pnMean:" << pnmean << endl;

      // Fill PN stuff
      //===============

      if (firstChanMod[iMod] == iCry && IsThereDataADC[iCry][iCol] == 1) {
        for (unsigned int ichan = 0; ichan < nPNPerMod; ichan++) {
          PNFirstAnal[iMod][ichan][iCol]->addEntry(pnmean, pn0, pn1);
        }
      }

      // Fill APD stuff
      //================

      if (apdAmpl != 0.0)
        APDFirstAnal[iCry][iCol]->addEntry(apdAmpl, pnmean, pn0, pn1, apdTime);
      if (_debug >= 1)
        cout << "-- debug test -- endJob -- filling APDTree" << endl;

      APDtrees[iCry]->Fill();

      // Fill reference trees
      //=====================

      if (apdRefMap[0][iMod + 1] == iCry || apdRefMap[1][iMod + 1] == iCry) {
        apdAmplA = 0.0;
        apdAmplB = 0.0;
        eventref = event;
        colorref = color;

        for (unsigned int ir = 0; ir < nRefChan; ir++) {
          if (_debug >= 1)
            cout << "-- debug test -- ir:" << ir << " tt:" << towerID << " refmap:" << apdRefMap[ir][iMod + 1]
                 << " iCry:" << iCry << endl;

          if (apdRefMap[ir][iMod + 1] == iCry) {
            if (_debug >= 1)
              cout << "-- debug test -- cut passed " << endl;
            if (ir == 0)
              apdAmplA = apdAmpl;
            else if (ir == 1)
              apdAmplB = apdAmpl;
            if (_debug >= 1)
              cout << "-- debug test -- apdAmplA=" << apdAmplA << endl;
            if (_debug >= 1)
              cout << "-- debug test -- apdAmplB=" << apdAmplB << endl;
            if (_debug >= 1)
              cout << "-- debug test -- color=" << color << ", event:" << event << ", ir:" << ir
                   << " tt-1:" << towerID - 1 << endl;

            RefAPDtrees[ir][iMod + 1]->Fill();

            if (_debug >= 1)
              cout << "-- debug test -- tree filled" << event << endl;
          }
        }
      }
    }
  }

  delete psfit;

  ADCFile->Close();

  if (_debug == 1)
    cout << "-- debug test -- endJob -- after apdAmpl Loop" << endl;

  // Remove temporary file
  //=======================

  stringstream del;
  del << "rm " << ADCfile;
  system(del.str().c_str());

  // Create output trees
  //=====================

  resFile = new TFile(resfile.c_str(), "RECREATE");

  for (unsigned int iColor = 0; iColor < nCol; iColor++) {
    stringstream nametree;
    nametree << "APDCol" << colors.at(iColor);
    stringstream nametree2;
    nametree2 << "PNCol" << colors.at(iColor);

    restrees[iColor] = new TTree(nametree.str().c_str(), nametree.str().c_str());
    respntrees[iColor] = new TTree(nametree2.str().c_str(), nametree2.str().c_str());

    restrees[iColor]->Branch("iphi", &iphi, "iphi/I");
    restrees[iColor]->Branch("ieta", &ieta, "ieta/I");
    restrees[iColor]->Branch("side", &side, "side/I");
    restrees[iColor]->Branch("dccID", &dccID, "dccID/I");
    restrees[iColor]->Branch("moduleID", &moduleID, "moduleID/I");
    restrees[iColor]->Branch("towerID", &towerID, "towerID/I");
    restrees[iColor]->Branch("channelID", &channelID, "channelID/I");
    restrees[iColor]->Branch("APD", &APD, "APD[6]/D");
    restrees[iColor]->Branch("Time", &Time, "Time[6]/D");
    restrees[iColor]->Branch("APDoPN", &APDoPN, "APDoPN[6]/D");
    restrees[iColor]->Branch("APDoPNA", &APDoPNA, "APDoPNA[6]/D");
    restrees[iColor]->Branch("APDoPNB", &APDoPNB, "APDoPNB[6]/D");
    restrees[iColor]->Branch("APDoAPDA", &APDoAPDA, "APDoAPDA[6]/D");
    restrees[iColor]->Branch("APDoAPDB", &APDoAPDB, "APDoAPDB[6]/D");
    restrees[iColor]->Branch("ShapeCor", &ShapeCor, "ShapeCor/D");
    restrees[iColor]->Branch("flag", &flag, "flag/I");

    respntrees[iColor]->Branch("moduleID", &moduleID, "moduleID/I");
    respntrees[iColor]->Branch("pnID", &pnID, "pnID/I");
    respntrees[iColor]->Branch("PN", &PN, "PN[6]/D");
    respntrees[iColor]->Branch("PNoPN", &PNoPN, "PNoPN[6]/D");
    respntrees[iColor]->Branch("PNoPNA", &PNoPNA, "PNoPNA[6]/D");
    respntrees[iColor]->Branch("PNoPNB", &PNoPNB, "PNoPNB[6]/D");

    restrees[iColor]->SetBranchAddress("iphi", &iphi);
    restrees[iColor]->SetBranchAddress("ieta", &ieta);
    restrees[iColor]->SetBranchAddress("dccID", &dccID);
    restrees[iColor]->SetBranchAddress("moduleID", &moduleID);
    restrees[iColor]->SetBranchAddress("towerID", &towerID);
    restrees[iColor]->SetBranchAddress("channelID", &channelID);
    restrees[iColor]->SetBranchAddress("APD", APD);
    restrees[iColor]->SetBranchAddress("Time", Time);
    restrees[iColor]->SetBranchAddress("APDoPN", APDoPN);
    restrees[iColor]->SetBranchAddress("APDoPNA", APDoPNA);
    restrees[iColor]->SetBranchAddress("APDoPNB", APDoPNB);
    restrees[iColor]->SetBranchAddress("APDoAPDA", APDoAPDA);
    restrees[iColor]->SetBranchAddress("APDoAPDB", APDoAPDB);
    restrees[iColor]->SetBranchAddress("ShapeCor", &ShapeCor);
    restrees[iColor]->SetBranchAddress("flag", &flag);

    respntrees[iColor]->SetBranchAddress("moduleID", &moduleID);
    respntrees[iColor]->SetBranchAddress("pnID", &pnID);
    respntrees[iColor]->SetBranchAddress("PN", PN);
    respntrees[iColor]->SetBranchAddress("PNoPN", PNoPN);
    respntrees[iColor]->SetBranchAddress("PNoPNA", PNoPNA);
    respntrees[iColor]->SetBranchAddress("PNoPNB", PNoPNB);
  }

  // Set Cuts for PN stuff
  //=======================

  for (unsigned int iM = 0; iM < nMod; iM++) {
    unsigned int iMod = modules[iM] - 1;

    for (unsigned int ich = 0; ich < nPNPerMod; ich++) {
      for (unsigned int icol = 0; icol < nCol; icol++) {
        PNAnal[iMod][ich][icol]->setPNCut(PNFirstAnal[iMod][ich][icol]->getPN().at(0),
                                          PNFirstAnal[iMod][ich][icol]->getPN().at(1));
      }
    }
  }

  // Build ref trees indexes
  //========================

  for (unsigned int imod = 0; imod < nMod; imod++) {
    int jmod = modules[imod];
    if (RefAPDtrees[0][jmod]->GetEntries() != 0 && RefAPDtrees[1][jmod]->GetEntries() != 0) {
      RefAPDtrees[0][jmod]->BuildIndex("eventref");
      RefAPDtrees[1][jmod]->BuildIndex("eventref");
    }
  }

  // Final loop on crystals
  //=======================

  for (unsigned int iCry = 0; iCry < nCrys; iCry++) {
    unsigned int iMod = iModule[iCry] - 1;

    // Set cuts on APD stuff
    //=======================

    for (unsigned int iCol = 0; iCol < nCol; iCol++) {
      std::vector<double> lowcut;
      std::vector<double> highcut;
      double cutMin;
      double cutMax;

      cutMin = APDFirstAnal[iCry][iCol]->getAPD().at(0) - 2.0 * APDFirstAnal[iCry][iCol]->getAPD().at(1);
      if (cutMin < 0)
        cutMin = 0;
      cutMax = APDFirstAnal[iCry][iCol]->getAPD().at(0) + 2.0 * APDFirstAnal[iCry][iCol]->getAPD().at(1);

      lowcut.push_back(cutMin);
      highcut.push_back(cutMax);

      cutMin = APDFirstAnal[iCry][iCol]->getTime().at(0) - 2.0 * APDFirstAnal[iCry][iCol]->getTime().at(1);
      cutMax = APDFirstAnal[iCry][iCol]->getTime().at(0) + 2.0 * APDFirstAnal[iCry][iCol]->getTime().at(1);
      lowcut.push_back(cutMin);
      highcut.push_back(cutMax);

      APDAnal[iCry][iCol] = new TAPD();
      APDAnal[iCry][iCol]->setAPDCut(APDFirstAnal[iCry][iCol]->getAPD().at(0),
                                     APDFirstAnal[iCry][iCol]->getAPD().at(1));
      APDAnal[iCry][iCol]->setAPDoPNCut(APDFirstAnal[iCry][iCol]->getAPDoPN().at(0),
                                        APDFirstAnal[iCry][iCol]->getAPDoPN().at(1));
      APDAnal[iCry][iCol]->setAPDoPN0Cut(APDFirstAnal[iCry][iCol]->getAPDoPN0().at(0),
                                         APDFirstAnal[iCry][iCol]->getAPDoPN0().at(1));
      APDAnal[iCry][iCol]->setAPDoPN1Cut(APDFirstAnal[iCry][iCol]->getAPDoPN1().at(0),
                                         APDFirstAnal[iCry][iCol]->getAPDoPN1().at(1));
      APDAnal[iCry][iCol]->setTimeCut(APDFirstAnal[iCry][iCol]->getTime().at(0),
                                      APDFirstAnal[iCry][iCol]->getTime().at(1));

      APDAnal[iCry][iCol]->set2DAPDoAPD0Cut(lowcut, highcut);
      APDAnal[iCry][iCol]->set2DAPDoAPD1Cut(lowcut, highcut);
    }

    // Final loop on events
    //=======================

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

      double pnmean;
      if (pn0 < 10 && pn1 > 10) {
        pnmean = pn1;
      } else if (pn1 < 10 && pn0 > 10) {
        pnmean = pn0;
      } else
        pnmean = 0.5 * (pn0 + pn1);

      // Get back color
      //===============

      unsigned int iCol = 0;
      for (unsigned int i = 0; i < nCol; i++) {
        if (color == colors[i]) {
          iCol = i;
          i = colors.size();
        }
      }

      // Fill PN stuff
      //===============

      if (firstChanMod[iMod] == iCry && IsThereDataADC[iCry][iCol] == 1) {
        for (unsigned int ichan = 0; ichan < nPNPerMod; ichan++) {
          PNAnal[iMod][ichan][iCol]->addEntry(pnmean, pn0, pn1);
        }
      }

      // Get ref amplitudes
      //===================

      if (_debug >= 1)
        cout << "-- debug test -- LastLoop event:" << event << " apdAmpl:" << apdAmpl << endl;
      apdAmplA = 0.0;
      apdAmplB = 0.0;

      for (unsigned int iRef = 0; iRef < nRefChan; iRef++) {
        RefAPDtrees[iRef][iMod + 1]->GetEntryWithIndex(event);
      }

      if (_debug == 1)
        cout << "-- debug test -- LastLoop apdAmplA:" << apdAmplA << " apdAmplB:" << apdAmplB << ", event:" << event
             << ", eventref:" << eventref << endl;

      // Fill APD stuff
      //===============

      APDAnal[iCry][iCol]->addEntry(apdAmpl, pnmean, pn0, pn1, apdTime, apdAmplA, apdAmplB);
    }

    moduleID = iMod + 1;

    if (moduleID >= 20)
      moduleID -= 2;  // Trick to fix endcap specificity

    // Get final results for APD
    //===========================

    for (unsigned int iColor = 0; iColor < nCol; iColor++) {
      std::vector<double> apdvec = APDAnal[iCry][iColor]->getAPD();
      std::vector<double> apdpnvec = APDAnal[iCry][iColor]->getAPDoPN();
      std::vector<double> apdpn0vec = APDAnal[iCry][iColor]->getAPDoPN0();
      std::vector<double> apdpn1vec = APDAnal[iCry][iColor]->getAPDoPN1();
      std::vector<double> timevec = APDAnal[iCry][iColor]->getTime();
      std::vector<double> apdapd0vec = APDAnal[iCry][iColor]->getAPDoAPD0();
      std::vector<double> apdapd1vec = APDAnal[iCry][iColor]->getAPDoAPD1();

      for (unsigned int i = 0; i < apdvec.size(); i++) {
        APD[i] = apdvec.at(i);
        APDoPN[i] = apdpnvec.at(i);
        APDoPNA[i] = apdpn0vec.at(i);
        APDoPNB[i] = apdpn1vec.at(i);
        APDoAPDA[i] = apdapd0vec.at(i);
        APDoAPDB[i] = apdapd1vec.at(i);
        Time[i] = timevec.at(i);
        ShapeCor = shapeCorrection;
      }

      // Fill APD results trees
      //========================

      iphi = iPhi[iCry];
      ieta = iEta[iCry];
      dccID = idccID[iCry];
      towerID = iTowerID[iCry];
      side = iside[iCry];
      channelID = iChannelID[iCry];

      if (!wasGainOK[iCry] || !wasTimingOK[iCry] || IsThereDataADC[iCry][iColor] == 0)
        flag = 1;
      else
        flag = 0;

      if (_debug >= 1)
        cout << "-- debug test -- endJob -- APD[0]" << APD[0] << " APDoPN[0] " << APDoPN[0] << " APDoAPDA[0] "
             << APDoAPDA[0] << " flag " << flag << endl;
      restrees[iColor]->Fill();
    }
  }

  // Get final results for PN
  //==========================

  for (unsigned int iM = 0; iM < nMod; iM++) {
    unsigned int iMod = modules[iM] - 1;

    side = iside[firstChanMod[iMod]];

    for (unsigned int ch = 0; ch < nPNPerMod; ch++) {
      pnID = ch;
      moduleID = iMod + 1;

      if (moduleID >= 20)
        moduleID -= 2;  // Trick to fix endcap specificity

      for (unsigned int iColor = 0; iColor < nCol; iColor++) {
        std::vector<double> pnvec = PNAnal[iMod][ch][iColor]->getPN();
        std::vector<double> pnopnvec = PNAnal[iMod][ch][iColor]->getPNoPN();
        std::vector<double> pnopn0vec = PNAnal[iMod][ch][iColor]->getPNoPN0();
        std::vector<double> pnopn1vec = PNAnal[iMod][ch][iColor]->getPNoPN1();

        for (unsigned int i = 0; i < pnvec.size(); i++) {
          PN[i] = pnvec.at(i);
          PNoPN[i] = pnopnvec.at(i);
          PNoPNA[i] = pnopn0vec.at(i);
          PNoPNB[i] = pnopn1vec.at(i);
        }

        if (_debug >= 1)
          cout << "-- debug test -- endJob -- filling pn results'tree: PN[0]:" << PN[0] << " iModule:" << iMod
               << " iColor:" << iColor << " ch:" << ch << endl;

        // Fill PN results trees
        //========================

        respntrees[iColor]->Fill();
      }
    }
  }

  // Remove temporary files
  //========================

  if (!_saveallevents) {
    APDFile->Close();
    stringstream del2;
    del2 << "rm " << APDfile;
    system(del2.str().c_str());

  } else {
    APDFile->cd();
    APDtrees[0]->Write();
    APDFile->Close();
    resFile->cd();
  }

  // Save results
  //===============

  for (unsigned int i = 0; i < nCol; i++) {
    restrees[i]->Write();
    respntrees[i]->Write();
  }

  resFile->Close();

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

bool EcalLaserAnalyzer2::getShapes ( )

private

Definition at line 1309 of file EcalLaserAnalyzer2.cc.

References _saveshapes, gather_cfg::cout, elecfile_, cuy::ii, createfilelist::int, IsMatacqOK, isSPRFine, findQualityFiles::jj, matfile, Skims_PA_cff::name, seedmultiplicitymonitor_newtracking_cfi::nBins, nSamplesShapes, PulseShape, Scenarios_cff::scale, shapeCorrection, shapefile, ShapeFile, shapes, shapesVec, ctpps_dqm_sourceclient-live_cfg::test, x, and y.

Referenced by beginJob().

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

  // Get Pulse From Matacq Analysis:
  //================================

  bool IsMatacqOK = false;

  int doesMatFileExist = 0;
  int doesMatShapeExist = 0;
  FILE* test2;
  TProfile* laserShape = nullptr;
  test2 = fopen(matfile.c_str(), "r");
  if (test2)
    doesMatFileExist = 1;

  TFile* MatShapeFile;
  if (doesMatFileExist == 1) {
    MatShapeFile = new TFile(matfile.c_str());
    laserShape = (TProfile*)MatShapeFile->Get("shapeLaser");
    if (laserShape) {
      doesMatShapeExist = 1;
      double y = laserShape->Integral("w");
      if (y != 0)
        laserShape->Scale(1.0 / y);
    }
  } else {
    cout << " ERROR! Matacq shape file not found !" << endl;
  }
  if (doesMatShapeExist)
    IsMatacqOK = true;

  // Get SPR from the average elec shape in SM6:
  //============================================

  int doesElecFileExist = 0;
  FILE* test;
  test = fopen(elecfile_.c_str(), "r");
  if (test)
    doesElecFileExist = 1;

  TFile* ElecShapesFile;
  TH1D* elecShape = nullptr;

  if (doesElecFileExist == 1) {
    ElecShapesFile = new TFile(elecfile_.c_str());
    stringstream name;
    name << "MeanElecShape";
    elecShape = (TH1D*)ElecShapesFile->Get(name.str().c_str());
    if (elecShape && doesMatShapeExist == 1) {
      double x = elecShape->GetMaximum();
      if (x != 0)
        elecShape->Scale(1.0 / x);
      isSPRFine = true;
    } else {
      isSPRFine = false;
    }

  } else {
    cout << " ERROR! Elec shape file not found !" << endl;
  }

  if (IsMatacqOK) {
    ShapeFile = new TFile(shapefile.c_str(), "RECREATE");

    unsigned int nBins = int(laserShape->GetEntries());
    assert(nSamplesShapes == nBins);
    double elec_jj, laser_iiMinusjj;
    double sum_jj;

    if (isSPRFine == true) {
      unsigned int nBins2 = int(elecShape->GetNbinsX());

      if (nBins2 < nBins) {
        cout << "EcalLaserAnalyzer2::getShapes: wrong configuration of the shapes' number of bins" << std::endl;
        isSPRFine = false;
      }
      assert(nSamplesShapes == nBins2);

      stringstream name;
      name << "PulseShape";

      PulseShape = new TProfile(name.str().c_str(), name.str().c_str(), nBins, -0.5, double(nBins) - 0.5);

      // shift shapes to have max close to the real APD max

      for (int ii = 0; ii < 50; ii++) {
        shapes[ii] = 0.0;
        PulseShape->Fill(ii, 0.0);
      }

      for (unsigned int ii = 0; ii < nBins - 50; ii++) {
        sum_jj = 0.0;
        for (unsigned int jj = 0; jj < ii; jj++) {
          elec_jj = elecShape->GetBinContent(jj + 1);
          laser_iiMinusjj = laserShape->GetBinContent(ii - jj + 1);
          sum_jj += elec_jj * laser_iiMinusjj;
        }
        PulseShape->Fill(ii + 50, sum_jj);
        shapes[ii + 50] = sum_jj;
      }

      double scale = PulseShape->GetMaximum();
      shapeCorrection = scale;

      if (scale != 0) {
        PulseShape->Scale(1.0 / scale);
        for (unsigned int ii = 0; ii < nBins; ii++) {
          shapesVec.push_back(shapes[ii] / scale);
        }
      }

      if (_saveshapes)
        PulseShape->Write();
    }
  }
  ShapeFile->Close();

  if (!_saveshapes) {
    stringstream del;
    del << "rm " << shapefile;
    system(del.str().c_str());
  }

  return IsMatacqOK;
}

void EcalLaserAnalyzer2::setGeomEB	(	int	etaG,
		int	phiG,
		int	module,
		int	tower,
		int	strip,
		int	xtal,
		int	apdRefTT,
		int	channel,
		int	lmr
	)

Definition at line 1436 of file EcalLaserAnalyzer2.cc.

References ME::apdRefChannels(), apdRefMap, channelID, KineDebug3::count(), dccID, eta, firstChanMod, iChannelID, idccID, iEta, iModule, iPhi, isFirstChanModFilled, iside, iTowerID, nRefChan, phi, MEEBGeom::side(), side, hgcalTowerProducer_cfi::tower, and towerID.

Referenced by analyze().

                                                                                                        {
  side = MEEBGeom::side(etaG, phiG);

  assert(module >= *min_element(modules.begin(), modules.end()) &&
         module <= *max_element(modules.begin(), modules.end()));

  eta = etaG;
  phi = phiG;
  channelID = 5 * (strip - 1) + xtal - 1;
  towerID = tower;

  vector<int> apdRefChan = ME::apdRefChannels(module, lmr);
  for (unsigned int iref = 0; iref < nRefChan; iref++) {
    if (channelID == apdRefChan[iref] && towerID == apdRefTT && apdRefMap[iref].count(module) == 0) {
      apdRefMap[iref][module] = channel;
    }
  }

  if (isFirstChanModFilled[module - 1] == 0) {
    firstChanMod[module - 1] = channel;
    isFirstChanModFilled[module - 1] = 1;
  }

  iEta[channel] = eta;
  iPhi[channel] = phi;
  iModule[channel] = module;
  iTowerID[channel] = towerID;
  iChannelID[channel] = channelID;
  idccID[channel] = dccID;
  iside[channel] = side;
}

void EcalLaserAnalyzer2::setGeomEE	(	int	etaG,
		int	phiG,
		int	iX,
		int	iY,
		int	iZ,
		int	module,
		int	tower,
		int	ch,
		int	apdRefTT,
		int	channel,
		int	lmr
	)

Definition at line 1469 of file EcalLaserAnalyzer2.cc.

References ME::apdRefChannels(), apdRefMap, channelID, KineDebug3::count(), dccID, DEFINE_FWK_MODULE, eta, firstChanMod, iChannelID, idccID, iEta, iModule, iPhi, isFirstChanModFilled, iside, iTowerID, nRefChan, phi, MEEEGeom::side(), side, hgcalTowerProducer_cfi::tower, and towerID.

Referenced by analyze().

                                                                                                                   {
  side = MEEEGeom::side(iX, iY, iZ);

  assert(module >= *min_element(modules.begin(), modules.end()) &&
         module <= *max_element(modules.begin(), modules.end()));

  eta = etaG;
  phi = phiG;
  channelID = ch;
  towerID = tower;

  vector<int> apdRefChan = ME::apdRefChannels(module, lmr);
  for (unsigned int iref = 0; iref < nRefChan; iref++) {
    if (channelID == apdRefChan[iref] && towerID == apdRefTT && apdRefMap[iref].count(module) == 0) {
      apdRefMap[iref][module] = channel;
    }
  }

  if (isFirstChanModFilled[module - 1] == 0) {
    firstChanMod[module - 1] = channel;
    isFirstChanModFilled[module - 1] = 1;
  }

  iEta[channel] = eta;
  iPhi[channel] = phi;
  iModule[channel] = module;
  iTowerID[channel] = towerID;
  iChannelID[channel] = channelID;
  idccID[channel] = dccID;
  iside[channel] = side;
}

Member Data Documentation

int EcalLaserAnalyzer2::_debug

private

Definition at line 89 of file EcalLaserAnalyzer2.h.

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

bool EcalLaserAnalyzer2::_docorpn

private

Definition at line 85 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

std::string EcalLaserAnalyzer2::_ecalPart

private

Definition at line 83 of file EcalLaserAnalyzer2.h.

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

int EcalLaserAnalyzer2::_fedid

private

Definition at line 86 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

unsigned int EcalLaserAnalyzer2::_firstsample

private

Definition at line 65 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

unsigned int EcalLaserAnalyzer2::_firstsamplePN

private

Definition at line 71 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

unsigned int EcalLaserAnalyzer2::_lastsample

private

Definition at line 66 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

unsigned int EcalLaserAnalyzer2::_lastsamplePN

private

Definition at line 72 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

unsigned int EcalLaserAnalyzer2::_niter

private

Definition at line 81 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::_noise

private

Definition at line 82 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

unsigned int EcalLaserAnalyzer2::_nsamples

private

Definition at line 63 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

unsigned int EcalLaserAnalyzer2::_nsamplesPN

private

Definition at line 69 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

unsigned int EcalLaserAnalyzer2::_presample

private

Definition at line 64 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

double EcalLaserAnalyzer2::_presamplecut

private

Definition at line 80 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

unsigned int EcalLaserAnalyzer2::_presamplePN

private

Definition at line 70 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

double EcalLaserAnalyzer2::_qualpercent

private

Definition at line 88 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::_ratiomaxcutlow

private

Definition at line 79 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

double EcalLaserAnalyzer2::_ratiomincuthigh

private

Definition at line 78 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

double EcalLaserAnalyzer2::_ratiomincutlow

private

Definition at line 77 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

unsigned int EcalLaserAnalyzer2::_samplemax

private

Definition at line 68 of file EcalLaserAnalyzer2.h.

unsigned int EcalLaserAnalyzer2::_samplemin

private

Definition at line 67 of file EcalLaserAnalyzer2.h.

bool EcalLaserAnalyzer2::_saveallevents

private

Definition at line 87 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

bool EcalLaserAnalyzer2::_saveshapes

private

Definition at line 84 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

unsigned int EcalLaserAnalyzer2::_timingcuthigh

private

Definition at line 74 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

unsigned int EcalLaserAnalyzer2::_timingcutlow

private

Definition at line 73 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

unsigned int EcalLaserAnalyzer2::_timingqualhigh

private

Definition at line 76 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

unsigned int EcalLaserAnalyzer2::_timingquallow

private

Definition at line 75 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

double EcalLaserAnalyzer2::adc[10]

private

Definition at line 181 of file EcalLaserAnalyzer2.h.

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

int EcalLaserAnalyzer2::adcC[1700]

private

Definition at line 226 of file EcalLaserAnalyzer2.h.

std::string EcalLaserAnalyzer2::ADCfile

private

Definition at line 110 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

TFile* EcalLaserAnalyzer2::ADCFile

private

Definition at line 162 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

int EcalLaserAnalyzer2::adcG[10]

private

Definition at line 182 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

double EcalLaserAnalyzer2::adcMean[1700][10]

private

Definition at line 225 of file EcalLaserAnalyzer2.h.

double* EcalLaserAnalyzer2::adcNoPed

private

Definition at line 195 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

TTree* EcalLaserAnalyzer2::ADCtrees[1700]

private

Definition at line 163 of file EcalLaserAnalyzer2.h.

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

double EcalLaserAnalyzer2::APD[6]

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::apdAmpl

private

Definition at line 186 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::apdAmplA

private

Definition at line 187 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::apdAmplB

private

Definition at line 188 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

TAPD* EcalLaserAnalyzer2::APDAnal[1700][nColor]

private

Definition at line 204 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

std::string EcalLaserAnalyzer2::APDfile

private

Definition at line 111 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

TFile* EcalLaserAnalyzer2::APDFile

private

Definition at line 165 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

TAPD* EcalLaserAnalyzer2::APDFirstAnal[1700][nColor]

private

Definition at line 203 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::APDoAPDA[6]

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::APDoAPDB[6]

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::APDoPN[6]

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::APDoPNA[6]

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::APDoPNB[6]

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

TAPDPulse* EcalLaserAnalyzer2::APDPulse

private

Definition at line 91 of file EcalLaserAnalyzer2.h.

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

std::map<int, unsigned int> EcalLaserAnalyzer2::apdRefMap[2]

private

Definition at line 150 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), setGeomEB(), and setGeomEE().

double EcalLaserAnalyzer2::apdTime

private

Definition at line 189 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

TTree* EcalLaserAnalyzer2::APDtrees[1700]

private

Definition at line 166 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

int EcalLaserAnalyzer2::channelID

private

Definition at line 136 of file EcalLaserAnalyzer2.h.

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

int EcalLaserAnalyzer2::channelIteratorEE

private

Definition at line 219 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

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

private

Definition at line 147 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

int EcalLaserAnalyzer2::color

private

Definition at line 180 of file EcalLaserAnalyzer2.h.

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

int EcalLaserAnalyzer2::colorref

private

Definition at line 193 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

std::vector<int> EcalLaserAnalyzer2::colors

private

Definition at line 146 of file EcalLaserAnalyzer2.h.

int EcalLaserAnalyzer2::dccID

private

Definition at line 138 of file EcalLaserAnalyzer2.h.

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

std::vector<int> EcalLaserAnalyzer2::dccMEM

private

Definition at line 148 of file EcalLaserAnalyzer2.h.

TMom* EcalLaserAnalyzer2::Delta01

private

Definition at line 94 of file EcalLaserAnalyzer2.h.

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

TMom* EcalLaserAnalyzer2::Delta12

private

Definition at line 95 of file EcalLaserAnalyzer2.h.

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

std::string EcalLaserAnalyzer2::digiCollection_

private

Definition at line 98 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

std::string EcalLaserAnalyzer2::digiPNCollection_

private

Definition at line 101 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

std::string EcalLaserAnalyzer2::digiProducer_

private

Definition at line 102 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

std::string EcalLaserAnalyzer2::elecfile_

private

Definition at line 99 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and getShapes().

int EcalLaserAnalyzer2::eta

private

Definition at line 178 of file EcalLaserAnalyzer2.h.

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

int EcalLaserAnalyzer2::event

private

Definition at line 179 of file EcalLaserAnalyzer2.h.

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

std::string EcalLaserAnalyzer2::eventHeaderCollection_

private

Definition at line 103 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

std::string EcalLaserAnalyzer2::eventHeaderProducer_

private

Definition at line 104 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

int EcalLaserAnalyzer2::eventref

private

Definition at line 192 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

int EcalLaserAnalyzer2::fedID

private

Definition at line 137 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

unsigned int EcalLaserAnalyzer2::firstChanMod[21]

private

Definition at line 231 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), endJob(), setGeomEB(), and setGeomEE().

int EcalLaserAnalyzer2::flag

private

Definition at line 218 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

int EcalLaserAnalyzer2::iChannelID[1700]

private

Definition at line 230 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), endJob(), setGeomEB(), and setGeomEE().

int EcalLaserAnalyzer2::idccID[1700]

private

Definition at line 230 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), endJob(), setGeomEB(), and setGeomEE().

int EcalLaserAnalyzer2::iEta[1700]

private

Definition at line 228 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), endJob(), setGeomEB(), and setGeomEE().

int EcalLaserAnalyzer2::iEvent

private

Definition at line 59 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

unsigned int EcalLaserAnalyzer2::iModule[1700]

private

Definition at line 229 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), endJob(), setGeomEB(), and setGeomEE().

int EcalLaserAnalyzer2::iPhi[1700]

private

Definition at line 228 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), endJob(), setGeomEB(), and setGeomEE().

unsigned int EcalLaserAnalyzer2::isFirstChanModFilled[21]

private

Definition at line 232 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), setGeomEB(), and setGeomEE().

bool EcalLaserAnalyzer2::isGainOK

private

Definition at line 243 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

int EcalLaserAnalyzer2::iside[1700]

private

Definition at line 230 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), endJob(), setGeomEB(), and setGeomEE().

bool EcalLaserAnalyzer2::IsMatacqOK

private

Definition at line 126 of file EcalLaserAnalyzer2.h.

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

bool EcalLaserAnalyzer2::isMatacqOK

private

Definition at line 213 of file EcalLaserAnalyzer2.h.

bool EcalLaserAnalyzer2::isSPRFine

private

Definition at line 214 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

int EcalLaserAnalyzer2::IsThereDataADC[1700][nColor]

private

Definition at line 206 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

bool EcalLaserAnalyzer2::isTimingOK

private

Definition at line 244 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

int EcalLaserAnalyzer2::iTowerID[1700]

private

Definition at line 230 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), endJob(), setGeomEB(), and setGeomEE().

int EcalLaserAnalyzer2::iZ

private

Definition at line 141 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

int EcalLaserAnalyzer2::laserEvents

private

Definition at line 144 of file EcalLaserAnalyzer2.h.

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

int EcalLaserAnalyzer2::lightside

private

Definition at line 140 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

std::string EcalLaserAnalyzer2::matfile

private

Definition at line 109 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and getShapes().

TMem* EcalLaserAnalyzer2::Mem

private

Definition at line 93 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

int EcalLaserAnalyzer2::moduleID

private

Definition at line 218 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

std::vector<int> EcalLaserAnalyzer2::modules

private

Definition at line 149 of file EcalLaserAnalyzer2.h.

unsigned int EcalLaserAnalyzer2::nCrys

private

Definition at line 117 of file EcalLaserAnalyzer2.h.

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

int EcalLaserAnalyzer2::nEvtBadGain[1700]

private

Definition at line 236 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

int EcalLaserAnalyzer2::nEvtBadTiming[1700]

private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

int EcalLaserAnalyzer2::nEvtTot[1700]

private

Definition at line 238 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

unsigned int EcalLaserAnalyzer2::nMod

private

Definition at line 121 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

unsigned int EcalLaserAnalyzer2::nPNPerMod

private

Definition at line 118 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

unsigned int EcalLaserAnalyzer2::nRefChan

private

Definition at line 119 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), endJob(), setGeomEB(), and setGeomEE().

unsigned int EcalLaserAnalyzer2::nRefTrees

private

Definition at line 120 of file EcalLaserAnalyzer2.h.

unsigned int EcalLaserAnalyzer2::nSamplesShapes

private

Definition at line 124 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

unsigned int EcalLaserAnalyzer2::nSides

private

Definition at line 122 of file EcalLaserAnalyzer2.h.

int EcalLaserAnalyzer2::phi

private

Definition at line 178 of file EcalLaserAnalyzer2.h.

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

double EcalLaserAnalyzer2::pn[50]

private

Definition at line 184 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

double EcalLaserAnalyzer2::PN[6]

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::pn0

private

Definition at line 183 of file EcalLaserAnalyzer2.h.

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

double EcalLaserAnalyzer2::pn1

private

Definition at line 183 of file EcalLaserAnalyzer2.h.

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

double EcalLaserAnalyzer2::pnAmpl

private

Definition at line 190 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

TPN* EcalLaserAnalyzer2::PNAnal[9][2][nColor]

private

Definition at line 200 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

std::string EcalLaserAnalyzer2::pncorfile_

private

Definition at line 100 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

TPNCor* EcalLaserAnalyzer2::pnCorrector

private

Definition at line 154 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

TPN* EcalLaserAnalyzer2::PNFirstAnal[9][2][nColor]

private

Definition at line 199 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

int EcalLaserAnalyzer2::pnG[50]

private

Definition at line 185 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

int EcalLaserAnalyzer2::pnID

private

Definition at line 218 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double* EcalLaserAnalyzer2::pnNoPed

private

Definition at line 196 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

double EcalLaserAnalyzer2::PNoPN[6]

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::PNoPNA[6]

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::PNoPNB[6]

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

TPNPulse* EcalLaserAnalyzer2::PNPulse

private

Definition at line 92 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

TProfile* EcalLaserAnalyzer2::PulseShape

private

Definition at line 174 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

TTree* EcalLaserAnalyzer2::RefAPDtrees[2][21]

private

Definition at line 167 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

std::string EcalLaserAnalyzer2::resdir_

private

Definition at line 97 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and EcalLaserAnalyzer2().

std::string EcalLaserAnalyzer2::resfile

private

Definition at line 112 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

TFile* EcalLaserAnalyzer2::resFile

private

Definition at line 169 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

TTree* EcalLaserAnalyzer2::respntrees[nColor]

private

Definition at line 171 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

TTree* EcalLaserAnalyzer2::restrees[nColor]

private

Definition at line 170 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

int EcalLaserAnalyzer2::runNum

private

Definition at line 131 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

int EcalLaserAnalyzer2::runType

private

Definition at line 130 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

double EcalLaserAnalyzer2::ShapeCor

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

double EcalLaserAnalyzer2::shapeCorrection

private

Definition at line 212 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

std::string EcalLaserAnalyzer2::shapefile

private

Definition at line 108 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and getShapes().

TFile* EcalLaserAnalyzer2::ShapeFile

private

Definition at line 173 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

double EcalLaserAnalyzer2::shapes[250]

private

Definition at line 211 of file EcalLaserAnalyzer2.h.

Referenced by DOTExport.DotExport::export(), getShapes(), DOTExport.DotProducer::produceLegend(), DOTExport.DotProducer::produceServices(), and DOTExport.DotProducer::produceSource().

std::vector<double> EcalLaserAnalyzer2::shapesVec

private

Definition at line 210 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

int EcalLaserAnalyzer2::side

private

Definition at line 139 of file EcalLaserAnalyzer2.h.

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

double EcalLaserAnalyzer2::Time[6]

private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

int EcalLaserAnalyzer2::towerID

private

Definition at line 135 of file EcalLaserAnalyzer2.h.

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

bool EcalLaserAnalyzer2::wasGainOK[1700]

private

Definition at line 240 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

bool EcalLaserAnalyzer2::wasTimingOK[1700]

private

Definition at line 241 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().