CMS 3D CMS Logo
List of all members | Public Types | Public Member Functions | Private Member Functions | Private Attributes
EcalLaserAnalyzer2 Class Reference

#include <EcalLaserAnalyzer2.h>

Inheritance diagram for EcalLaserAnalyzer2:
edm::EDAnalyzer edm::EDConsumerBase

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 &&)=default
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESProxyIndex > 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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
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 wantsInputProcessBlocks ()
 
static bool wantsProcessBlocks ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
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<Transition Tr = Transition::Event>
constexpr auto esConsumes () noexcept
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag) noexcept
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 

Detailed Description

Definition at line 43 of file EcalLaserAnalyzer2.h.

Member Enumeration Documentation

◆ VarCol

enum EcalLaserAnalyzer2::VarCol
Enumerator
iBlue 
iRed 
nColor 

Definition at line 56 of file EcalLaserAnalyzer2.h.

56 { iBlue, iRed, nColor };

Constructor & Destructor Documentation

◆ EcalLaserAnalyzer2()

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

Definition at line 56 of file EcalLaserAnalyzer2.cc.

58  : iEvent(0),
59 
60  // framework parameters with default values
61  _nsamples(iConfig.getUntrackedParameter<unsigned int>("nSamples", 10)),
62  _presample(iConfig.getUntrackedParameter<unsigned int>("nPresamples", 2)),
63  _firstsample(iConfig.getUntrackedParameter<unsigned int>("firstSample", 1)),
64  _lastsample(iConfig.getUntrackedParameter<unsigned int>("lastSample", 2)),
65  _nsamplesPN(iConfig.getUntrackedParameter<unsigned int>("nSamplesPN", 50)),
66  _presamplePN(iConfig.getUntrackedParameter<unsigned int>("nPresamplesPN", 6)),
67  _firstsamplePN(iConfig.getUntrackedParameter<unsigned int>("firstSamplePN", 7)),
68  _lastsamplePN(iConfig.getUntrackedParameter<unsigned int>("lastSamplePN", 8)),
69  _timingcutlow(iConfig.getUntrackedParameter<unsigned int>("timingCutLow", 2)),
70  _timingcuthigh(iConfig.getUntrackedParameter<unsigned int>("timingCutHigh", 9)),
71  _timingquallow(iConfig.getUntrackedParameter<unsigned int>("timingQualLow", 3)),
72  _timingqualhigh(iConfig.getUntrackedParameter<unsigned int>("timingQualHigh", 8)),
73  _ratiomincutlow(iConfig.getUntrackedParameter<double>("ratioMinCutLow", 0.4)),
74  _ratiomincuthigh(iConfig.getUntrackedParameter<double>("ratioMinCutHigh", 0.95)),
75  _ratiomaxcutlow(iConfig.getUntrackedParameter<double>("ratioMaxCutLow", 0.8)),
76  _presamplecut(iConfig.getUntrackedParameter<double>("presampleCut", 5.0)),
77  _niter(iConfig.getUntrackedParameter<unsigned int>("nIter", 5)),
78  _noise(iConfig.getUntrackedParameter<double>("noise", 2.0)),
79  _ecalPart(iConfig.getUntrackedParameter<std::string>("ecalPart", "EB")),
80  _saveshapes(iConfig.getUntrackedParameter<bool>("saveShapes", true)),
81  _docorpn(iConfig.getUntrackedParameter<bool>("doCorPN", false)),
82  _fedid(iConfig.getUntrackedParameter<int>("fedID", -999)),
83  _saveallevents(iConfig.getUntrackedParameter<bool>("saveAllEvents", false)),
84  _qualpercent(iConfig.getUntrackedParameter<double>("qualPercent", 0.2)),
85  _debug(iConfig.getUntrackedParameter<int>("debug", 0)),
86  nCrys(NCRYSEB),
87  nPNPerMod(NPNPERMOD),
88  nMod(NMODEB),
89  nSides(NSIDES),
90  nSamplesShapes(NSAMPSHAPES),
91  IsMatacqOK(false),
92  runType(-1),
93  runNum(0),
94  towerID(-1),
95  channelID(-1),
96  fedID(-1),
97  dccID(-1),
98  side(2),
99  lightside(2),
100  iZ(1),
101  phi(-1),
102  eta(-1),
103  event(0),
104  color(0),
105  pn0(0),
106  pn1(0),
107  apdAmpl(0),
108  apdTime(0),
109  pnAmpl(0),
110  pnID(-1),
111  moduleID(-1),
112  flag(0),
113  channelIteratorEE(0),
114  ShapeCor(0)
115 
116 //========================================================================
117 
118 {
119  // Initialization from cfg file
120 
121  resdir_ = iConfig.getUntrackedParameter<std::string>("resDir");
122  elecfile_ = iConfig.getUntrackedParameter<std::string>("elecFile");
123  pncorfile_ = iConfig.getUntrackedParameter<std::string>("pnCorFile");
124 
125  digiCollection_ = iConfig.getParameter<std::string>("digiCollection");
126  digiPNCollection_ = iConfig.getParameter<std::string>("digiPNCollection");
127  digiProducer_ = iConfig.getParameter<std::string>("digiProducer");
128 
129  eventHeaderCollection_ = iConfig.getParameter<std::string>("eventHeaderCollection");
130  eventHeaderProducer_ = iConfig.getParameter<std::string>("eventHeaderProducer");
131 
132  // Geometrical constants initialization
133  if (_ecalPart == "EB") {
134  nCrys = NCRYSEB;
135  } else {
136  nCrys = NCRYSEE;
137  }
138  iZ = 1;
139  if (_fedid <= 609)
140  iZ = -1;
141  modules = ME::lmmodFromDcc(_fedid);
142  nMod = modules.size();
143  nRefChan = NREFCHAN;
144 
145  for (unsigned int j = 0; j < nCrys; j++) {
146  iEta[j] = -1;
147  iPhi[j] = -1;
148  iModule[j] = 10;
149  iTowerID[j] = -1;
150  iChannelID[j] = -1;
151  idccID[j] = -1;
152  iside[j] = -1;
153  wasTimingOK[j] = true;
154  wasGainOK[j] = true;
155  }
156 
157  for (unsigned int j = 0; j < nMod; j++) {
158  int ii = modules[j];
159  firstChanMod[ii - 1] = 0;
160  isFirstChanModFilled[ii - 1] = 0;
161  }
162 
163  // Quality check flags
164 
165  isGainOK = true;
166  isTimingOK = true;
167 
168  // PN linearity corrector
169 
170  pnCorrector = new TPNCor(pncorfile_);
171 
172  // Objects dealing with pulses
173 
174  APDPulse = new TAPDPulse(_nsamples,
175  _presample,
176  _firstsample,
177  _lastsample,
178  _timingcutlow,
179  _timingcuthigh,
180  _timingquallow,
181  _timingqualhigh,
182  _ratiomincutlow,
183  _ratiomincuthigh,
184  _ratiomaxcutlow);
185  PNPulse = new TPNPulse(_nsamplesPN, _presamplePN);
186 
187  // Object dealing with MEM numbering
188 
189  Mem = new TMem(_fedid);
190 
191  // Objects needed for npresample calculation
192 
193  Delta01 = new TMom();
194  Delta12 = new TMom();
195 }

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.

◆ ~EcalLaserAnalyzer2()

EcalLaserAnalyzer2::~EcalLaserAnalyzer2 ( )
override

Definition at line 198 of file EcalLaserAnalyzer2.cc.

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

Member Function Documentation

◆ analyze()

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

Implements edm::EDAnalyzer.

Definition at line 280 of file EcalLaserAnalyzer2.cc.

280  {
281  //========================================================================
282 
283  ++iEvent;
284 
285  // retrieving DCC header
286  edm::Handle<EcalRawDataCollection> pDCCHeader;
287  const EcalRawDataCollection* DCCHeader = nullptr;
288  try {
289  e.getByLabel(eventHeaderProducer_, eventHeaderCollection_, pDCCHeader);
290  DCCHeader = pDCCHeader.product();
291  } catch (std::exception& ex) {
292  std::cerr << "Error! can't get the product retrieving DCC header" << eventHeaderCollection_.c_str() << std::endl;
293  }
294 
295  //retrieving crystal data from Event
296  edm::Handle<EBDigiCollection> pEBDigi;
297  const EBDigiCollection* EBDigi = nullptr;
298  edm::Handle<EEDigiCollection> pEEDigi;
299  const EEDigiCollection* EEDigi = nullptr;
300 
301  if (_ecalPart == "EB") {
302  try {
303  e.getByLabel(digiProducer_, digiCollection_, pEBDigi);
304  EBDigi = pEBDigi.product();
305  } catch (std::exception& ex) {
306  std::cerr << "Error! can't get the product retrieving EB crystal data " << digiCollection_.c_str() << std::endl;
307  }
308  } else if (_ecalPart == "EE") {
309  try {
310  e.getByLabel(digiProducer_, digiCollection_, pEEDigi);
311  EEDigi = pEEDigi.product();
312  } catch (std::exception& ex) {
313  std::cerr << "Error! can't get the product retrieving EE crystal data " << digiCollection_.c_str() << std::endl;
314  }
315  } else {
316  cout << " Wrong ecalPart in cfg file " << endl;
317  return;
318  }
319 
320  // retrieving crystal PN diodes from Event
321 
322  edm::Handle<EcalPnDiodeDigiCollection> pPNDigi;
323  const EcalPnDiodeDigiCollection* PNDigi = nullptr;
324  try {
325  e.getByLabel(digiProducer_, digiPNCollection_, pPNDigi);
326  PNDigi = pPNDigi.product();
327  } catch (std::exception& ex) {
328  std::cerr << "Error! can't get the product " << digiPNCollection_.c_str() << std::endl;
329  }
330 
331  // retrieving electronics mapping
332  edm::ESHandle<EcalElectronicsMapping> ecalmapping;
333  const EcalElectronicsMapping* TheMapping = nullptr;
334  try {
335  c.get<EcalMappingRcd>().get(ecalmapping);
336  TheMapping = ecalmapping.product();
337  } catch (std::exception& ex) {
338  std::cerr << "Error! can't get the product EcalMappingRcd" << std::endl;
339  }
340 
341  // ====================================
342  // Decode Basic DCCHeader Information
343  // ====================================
344 
345  for (EcalRawDataCollection::const_iterator headerItr = DCCHeader->begin(); headerItr != DCCHeader->end();
346  ++headerItr) {
347  // Get run type and run number
348 
349  int fed = headerItr->fedId();
350  if (fed != _fedid && _fedid != -999)
351  continue;
352 
353  runType = headerItr->getRunType();
354  runNum = headerItr->getRunNumber();
355  event = headerItr->getLV1();
356 
357  dccID = headerItr->getDccInTCCCommand();
358  fedID = headerItr->fedId();
359  lightside = headerItr->getRtHalf();
360 
361  // Check fed corresponds to the DCC in TCC
362 
363  if (600 + dccID != fedID)
364  continue;
365 
366  // Cut on runType
367 
368  if (runType != EcalDCCHeaderBlock::LASER_STD && runType != EcalDCCHeaderBlock::LASER_GAP &&
369  runType != EcalDCCHeaderBlock::LASER_POWER_SCAN && runType != EcalDCCHeaderBlock::LASER_DELAY_SCAN)
370  return;
371 
372  // Retrieve laser color and event number
373 
374  EcalDCCHeaderBlock::EcalDCCEventSettings settings = headerItr->getEventSettings();
375  color = settings.wavelength;
376  if (color < 0)
377  return;
378 
379  vector<int>::iterator iter = find(colors.begin(), colors.end(), color);
380  if (iter == colors.end()) {
381  colors.push_back(color);
382  cout << " new color found " << color << " " << colors.size() << endl;
383  }
384  }
385 
386  // Check Matacq shape exists
387 
388  if (!IsMatacqOK)
389  return;
390 
391  // Cut on fedID
392 
393  if (fedID != _fedid && _fedid != -999)
394  return;
395 
396  // Count laser events
397 
398  laserEvents++;
399 
400  // ======================
401  // Decode PN Information
402  // ======================
403 
404  TPNFit* pnfit = new TPNFit();
405  pnfit->init(_nsamplesPN, _firstsamplePN, _lastsamplePN);
406 
407  double chi2pn = 0;
408  unsigned int samplemax = 0;
409  int pnGain = 0;
410 
411  map<int, vector<double> > allPNAmpl;
412  map<int, vector<double> > allPNGain;
413 
414  // Loop on PNs digis
415 
416  for (EcalPnDiodeDigiCollection::const_iterator pnItr = PNDigi->begin(); pnItr != PNDigi->end(); ++pnItr) {
417  EcalPnDiodeDetId pnDetId = EcalPnDiodeDetId((*pnItr).id());
418 
419  if (_debug == 1)
420  cout << "-- debug test -- Inside PNDigi - pnID=" << pnDetId.iPnId() << ", dccID=" << pnDetId.iDCCId() << endl;
421 
422  // Skip MEM DCC without relevant data
423 
424  bool isMemRelevant = Mem->isMemRelevant(pnDetId.iDCCId());
425  if (!isMemRelevant)
426  continue;
427 
428  // Loop on PN samples
429 
430  for (int samId = 0; samId < (*pnItr).size(); samId++) {
431  pn[samId] = (*pnItr).sample(samId).adc();
432  pnG[samId] = (*pnItr).sample(samId).gainId();
433  if (samId == 0)
434  pnGain = pnG[samId];
435  if (samId > 0)
436  pnGain = int(TMath::Max(pnG[samId], pnGain));
437  }
438 
439  if (pnGain != 1)
440  cout << "PN gain different from 1" << endl;
441 
442  // Calculate amplitude from pulse
443 
444  PNPulse->setPulse(pn);
445  pnNoPed = PNPulse->getAdcWithoutPedestal();
446  samplemax = PNPulse->getMaxSample();
447  chi2pn = pnfit->doFit(samplemax, &pnNoPed[0]);
448  if (chi2pn == 101 || chi2pn == 102 || chi2pn == 103)
449  pnAmpl = 0.;
450  else
451  pnAmpl = pnfit->getAmpl();
452 
453  // Apply linearity correction
454 
455  double corr = 1.0;
456  if (_docorpn)
457  corr = pnCorrector->getPNCorrectionFactor(pnAmpl, pnGain);
458  pnAmpl *= corr;
459 
460  // Fill PN ampl vector
461 
462  allPNAmpl[pnDetId.iDCCId()].push_back(pnAmpl);
463 
464  if (_debug == 1)
465  cout << "-- debug -- Inside PNDigi - PNampl=" << pnAmpl << ", PNgain=" << pnGain << endl;
466  }
467 
468  // ===========================
469  // Decode EBDigis Information
470  // ===========================
471 
472  int adcGain = 0;
473 
474  if (EBDigi) {
475  // Loop on crystals
476  //===================
477 
478  for (EBDigiCollection::const_iterator digiItr = EBDigi->begin(); digiItr != EBDigi->end();
479  ++digiItr) { // Loop on crystals
480 
481  // Retrieve geometry
482  //===================
483 
484  EBDetId id_crystal(digiItr->id());
485  EBDataFrame df(*digiItr);
486  EcalElectronicsId elecid_crystal = TheMapping->getElectronicsId(id_crystal);
487 
488  int etaG = id_crystal.ieta(); // global
489  int phiG = id_crystal.iphi(); // global
490 
491  std::pair<int, int> LocalCoord = MEEBGeom::localCoord(etaG, phiG);
492 
493  int etaL = LocalCoord.first; // local
494  int phiL = LocalCoord.second; // local
495 
496  int strip = elecid_crystal.stripId();
497  int xtal = elecid_crystal.xtalId();
498 
499  int module = MEEBGeom::lmmod(etaG, phiG);
500  int tower = elecid_crystal.towerId();
501 
502  int apdRefTT = MEEBGeom::apdRefTower(module);
503 
504  std::pair<int, int> pnpair = MEEBGeom::pn(module);
505  unsigned int MyPn0 = pnpair.first;
506  unsigned int MyPn1 = pnpair.second;
507 
508  int lmr = MEEBGeom::lmr(etaG, phiG);
509  unsigned int channel = MEEBGeom::electronic_channel(etaL, phiL);
510  assert(channel < nCrys);
511 
512  setGeomEB(etaG, phiG, module, tower, strip, xtal, apdRefTT, channel, lmr);
513 
514  if (_debug == 1)
515  cout << "-- debug -- Inside EBDigi - towerID:" << towerID << " channelID:" << channelID << " module:" << module
516  << " modules:" << modules.size() << endl;
517 
518  // APD Pulse
519  //===========
520 
521  // Loop on adc samples
522 
523  for (unsigned int i = 0; i < (*digiItr).size(); ++i) {
524  EcalMGPASample samp_crystal(df.sample(i));
525  adc[i] = samp_crystal.adc();
526  adcG[i] = samp_crystal.gainId();
527  adc[i] *= adcG[i];
528  if (i == 0)
529  adcGain = adcG[i];
530  if (i > 0)
531  adcGain = TMath::Max(adcG[i], adcGain);
532  }
533 
534  APDPulse->setPulse(adc);
535 
536  // Quality checks
537  //================
538 
539  if (adcGain != 1)
540  nEvtBadGain[channel]++;
541  if (!APDPulse->isTimingQualOK())
542  nEvtBadTiming[channel]++;
543  nEvtTot[channel]++;
544 
545  // Associate PN ampl
546  //===================
547 
548  int mem0 = Mem->Mem(lmr, 0);
549  int mem1 = Mem->Mem(lmr, 1);
550 
551  if (allPNAmpl[mem0].size() > MyPn0)
552  pn0 = allPNAmpl[mem0][MyPn0];
553  else
554  pn0 = 0;
555  if (allPNAmpl[mem1].size() > MyPn1)
556  pn1 = allPNAmpl[mem1][MyPn1];
557  else
558  pn1 = 0;
559 
560  // Fill if Pulse is fine
561  //=======================
562 
563  if (APDPulse->isPulseOK() && lightside == side) {
564  ADCtrees[channel]->Fill();
565 
566  Delta01->addEntry(APDPulse->getDelta(0, 1));
567  Delta12->addEntry(APDPulse->getDelta(1, 2));
568  }
569  }
570  } else if (EEDigi) {
571  // Loop on crystals
572  //===================
573 
574  for (EEDigiCollection::const_iterator digiItr = EEDigi->begin(); digiItr != EEDigi->end(); ++digiItr) {
575  // Retrieve geometry
576  //===================
577 
578  EEDetId id_crystal(digiItr->id());
579  EEDataFrame df(*digiItr);
580  EcalElectronicsId elecid_crystal = TheMapping->getElectronicsId(id_crystal);
581 
582  int etaG = id_crystal.iy();
583  int phiG = id_crystal.ix();
584 
585  int iX = (phiG - 1) / 5 + 1;
586  int iY = (etaG - 1) / 5 + 1;
587 
588  int tower = elecid_crystal.towerId();
589  int ch = elecid_crystal.channelId() - 1;
590 
591  int module = MEEEGeom::lmmod(iX, iY);
592  if (module >= 18 && side == 1)
593  module += 2;
594  int lmr = MEEEGeom::lmr(iX, iY, iZ);
595  int dee = MEEEGeom::dee(lmr);
596  int apdRefTT = MEEEGeom::apdRefTower(lmr, module);
597 
598  std::pair<int, int> pnpair = MEEEGeom::pn(dee, module);
599  unsigned int MyPn0 = pnpair.first;
600  unsigned int MyPn1 = pnpair.second;
601 
602  int hashedIndex = 100000 * eta + phi;
603  if (channelMapEE.count(hashedIndex) == 0) {
604  channelMapEE[hashedIndex] = channelIteratorEE;
605  channelIteratorEE++;
606  }
607  unsigned int channel = channelMapEE[hashedIndex];
608  assert(channel < nCrys);
609 
610  setGeomEE(etaG, phiG, iX, iY, iZ, module, tower, ch, apdRefTT, channel, lmr);
611 
612  if (_debug == 1)
613  cout << "-- debug -- Inside EEDigi - towerID:" << towerID << " channelID:" << channelID << " module:" << module
614  << " modules:" << modules.size() << endl;
615 
616  // APD Pulse
617  //===========
618 
619  if ((*digiItr).size() > 10)
620  cout << "SAMPLES SIZE > 10!" << (*digiItr).size() << endl;
621 
622  // Loop on adc samples
623 
624  for (unsigned int i = 0; i < (*digiItr).size(); ++i) {
625  EcalMGPASample samp_crystal(df.sample(i));
626  adc[i] = samp_crystal.adc();
627  adcG[i] = samp_crystal.gainId();
628  adc[i] *= adcG[i];
629 
630  if (i == 0)
631  adcGain = adcG[i];
632  if (i > 0)
633  adcGain = TMath::Max(adcG[i], adcGain);
634  }
635 
636  APDPulse->setPulse(adc);
637 
638  // Quality checks
639  //================
640 
641  if (adcGain != 1)
642  nEvtBadGain[channel]++;
643  if (!APDPulse->isTimingQualOK())
644  nEvtBadTiming[channel]++;
645  nEvtTot[channel]++;
646 
647  // Associate PN ampl
648  //===================
649 
650  int mem0 = Mem->Mem(lmr, 0);
651  int mem1 = Mem->Mem(lmr, 1);
652 
653  if (allPNAmpl[mem0].size() > MyPn0)
654  pn0 = allPNAmpl[mem0][MyPn0];
655  else
656  pn0 = 0;
657  if (allPNAmpl[mem1].size() > MyPn1)
658  pn1 = allPNAmpl[mem1][MyPn1];
659  else
660  pn1 = 0;
661 
662  // Fill if Pulse is fine
663  //=======================
664 
665  if (APDPulse->isPulseOK() && lightside == side) {
666  ADCtrees[channel]->Fill();
667 
668  Delta01->addEntry(APDPulse->getDelta(0, 1));
669  Delta12->addEntry(APDPulse->getDelta(1, 2));
670  }
671  }
672  }
673 }

References _debug, _docorpn, _ecalPart, _fedid, _firstsamplePN, _lastsamplePN, _nsamplesPN, adc, adcG, ADCtrees, TMom::addEntry(), APDPulse, MEEBGeom::apdRefTower(), MEEEGeom::apdRefTower(), cms::cuda::assert(), edm::SortedCollection< T, SORT >::begin(), HltBtagPostValidation_cff::c, EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, EcalElectronicsId::channelId(), channelID, channelIteratorEE, channelMapEE, color, alignCSCRings::corr, gather_cfg::cout, dccID, MEEEGeom::dee(), Delta01, Delta12, hgcalPerformanceValidation::df, digiCollection_, digiPNCollection_, digiProducer_, TPNFit::doFit(), MillePedeFileConverter_cfg::e, CollectionTags_cfi::EBDigi, CollectionTags_cfi::EEDigi, MEEBGeom::electronic_channel(), edm::SortedCollection< T, SORT >::end(), eta, eventHeaderCollection_, eventHeaderProducer_, cppFunctionSkipper::exception, fedID, spr::find(), get, TPNPulse::getAdcWithoutPedestal(), TPNFit::getAmpl(), 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().

◆ beginJob()

void EcalLaserAnalyzer2::beginJob ( void  )
overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 206 of file EcalLaserAnalyzer2.cc.

206  {
207  //========================================================================
208 
209  // Create temporary files and trees to save adc samples
210  //======================================================
211 
212  ADCfile = resdir_;
213  ADCfile += "/APDSamplesLaser.root";
214 
215  APDfile = resdir_;
216  APDfile += "/APDPNLaserAllEvents.root";
217 
218  ADCFile = new TFile(ADCfile.c_str(), "RECREATE");
219 
220  for (unsigned int i = 0; i < nCrys; i++) {
221  stringstream name;
222  name << "ADCTree" << i + 1;
223  ADCtrees[i] = new TTree(name.str().c_str(), name.str().c_str());
224 
225  ADCtrees[i]->Branch("iphi", &phi, "phi/I");
226  ADCtrees[i]->Branch("ieta", &eta, "eta/I");
227  ADCtrees[i]->Branch("towerID", &towerID, "towerID/I");
228  ADCtrees[i]->Branch("channelID", &channelID, "channelID/I");
229  ADCtrees[i]->Branch("dccID", &dccID, "dccID/I");
230  ADCtrees[i]->Branch("side", &side, "side/I");
231  ADCtrees[i]->Branch("event", &event, "event/I");
232  ADCtrees[i]->Branch("color", &color, "color/I");
233  ADCtrees[i]->Branch("adc", &adc, "adc[10]/D");
234  ADCtrees[i]->Branch("pn0", &pn0, "pn0/D");
235  ADCtrees[i]->Branch("pn1", &pn1, "pn1/D");
236 
237  ADCtrees[i]->SetBranchAddress("ieta", &eta);
238  ADCtrees[i]->SetBranchAddress("iphi", &phi);
239  ADCtrees[i]->SetBranchAddress("towerID", &towerID);
240  ADCtrees[i]->SetBranchAddress("channelID", &channelID);
241  ADCtrees[i]->SetBranchAddress("dccID", &dccID);
242  ADCtrees[i]->SetBranchAddress("side", &side);
243  ADCtrees[i]->SetBranchAddress("event", &event);
244  ADCtrees[i]->SetBranchAddress("color", &color);
245  ADCtrees[i]->SetBranchAddress("adc", adc);
246  ADCtrees[i]->SetBranchAddress("pn0", &pn0);
247  ADCtrees[i]->SetBranchAddress("pn1", &pn1);
248  }
249 
250  // Define output results filenames
251  //==================================
252  stringstream namefile1;
253  namefile1 << resdir_ << "/SHAPE_LASER.root";
254  shapefile = namefile1.str();
255 
256  stringstream namefile2;
257  namefile2 << resdir_ << "/APDPN_LASER.root";
258  resfile = namefile2.str();
259 
260  stringstream namefile3;
261  namefile3 << resdir_ << "/MATACQ.root";
262 
263  matfile = namefile3.str();
264 
265  // Get Pulse Shapes
266  //==================
267 
268  IsMatacqOK = getShapes();
269  if (!IsMatacqOK) {
270  cout << " ERROR! No matacq shape available: analysis aborted !" << endl;
271  return;
272  }
273 
274  // Laser events counter
275 
276  laserEvents = 0;
277 }

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.

◆ endJob()

void EcalLaserAnalyzer2::endJob ( void  )
overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 677 of file EcalLaserAnalyzer2.cc.

677  {
678  //========================================================================
679 
680  if (!IsMatacqOK) {
681  cout << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+" << endl;
682  cout << "\t+=+ WARNING! NO MATACQ +=+" << endl;
683  cout << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+" << endl;
684  return;
685  }
686 
687  // Adjust channel numbers for EE
688  //===============================
689 
690  if (_ecalPart == "EE") {
691  nCrys = channelMapEE.size();
692  }
693 
694  // Set presamples number
695  //======================
696 
697  double delta01 = Delta01->getMean();
698  double delta12 = Delta12->getMean();
699  if (delta12 > _presamplecut) {
700  _presample = 2;
701  if (delta01 > _presamplecut)
702  _presample = 1;
703  }
704 
705  APDPulse->setPresamples(_presample);
706 
707  // Don't do anything if there is no events
708  //=========================================
709 
710  if (laserEvents == 0) {
711  ADCFile->Close();
712  stringstream del;
713  del << "rm " << ADCfile;
714  system(del.str().c_str());
715  cout << " No Laser Events " << endl;
716  return;
717  }
718 
719  // Set quality flags for gains and timing
720  //=========================================
721 
722  double BadGainEvtPercentage = 0.0;
723  double BadTimingEvtPercentage = 0.0;
724 
725  int nChanBadGain = 0;
726  int nChanBadTiming = 0;
727 
728  for (unsigned int i = 0; i < nCrys; i++) {
729  if (nEvtTot[i] != 0) {
730  BadGainEvtPercentage = double(nEvtBadGain[i]) / double(nEvtTot[i]);
731  BadTimingEvtPercentage = double(nEvtBadTiming[i]) / double(nEvtTot[i]);
732  }
733  if (BadGainEvtPercentage > _qualpercent) {
734  wasGainOK[i] = false;
735  nChanBadGain++;
736  }
737  if (BadTimingEvtPercentage > _qualpercent) {
738  wasTimingOK[i] = false;
739  nChanBadTiming++;
740  }
741  }
742 
743  double BadGainChanPercentage = double(nChanBadGain) / double(nCrys);
744  double BadTimingChanPercentage = double(nChanBadTiming) / double(nCrys);
745 
746  if (BadGainChanPercentage > _qualpercent)
747  isGainOK = false;
748  if (BadTimingChanPercentage > _qualpercent)
749  isTimingOK = false;
750 
751  // Analyze adc samples to get amplitudes
752  //=======================================
753 
754  cout << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+" << endl;
755  cout << "\t+=+ Analyzing laser data: getting APD, PN, APD/PN, PN/PN +=+" << endl;
756 
757  if (!isGainOK)
758  cout << "\t+=+ ............................ WARNING! APD GAIN WAS NOT 1 +=+" << endl;
759  if (!isTimingOK)
760  cout << "\t+=+ ............................ WARNING! TIMING WAS BAD +=+" << endl;
761 
762  APDFile = new TFile(APDfile.c_str(), "RECREATE");
763 
764  int ieta, iphi;
765 
766  int flagfit;
767  for (unsigned int i = 0; i < nCrys; i++) {
768  stringstream name;
769  name << "APDTree" << i + 1;
770 
771  APDtrees[i] = new TTree(name.str().c_str(), name.str().c_str());
772 
773  //List of branches
774 
775  APDtrees[i]->Branch("event", &event, "event/I");
776  APDtrees[i]->Branch("color", &color, "color/I");
777  APDtrees[i]->Branch("iphi", &iphi, "iphi/I");
778  APDtrees[i]->Branch("ieta", &ieta, "ieta/I");
779  APDtrees[i]->Branch("side", &side, "side/I");
780  APDtrees[i]->Branch("dccID", &dccID, "dccID/I");
781  APDtrees[i]->Branch("towerID", &towerID, "towerID/I");
782  APDtrees[i]->Branch("channelID", &channelID, "channelID/I");
783  APDtrees[i]->Branch("apdAmpl", &apdAmpl, "apdAmpl/D");
784  APDtrees[i]->Branch("apdTime", &apdTime, "apdTime/D");
785  if (_saveallevents)
786  APDtrees[i]->Branch("adc", &adc, "adc[10]/D");
787  APDtrees[i]->Branch("flagfit", &flagfit, "flagfit/I");
788  APDtrees[i]->Branch("pn0", &pn0, "pn0/D");
789  APDtrees[i]->Branch("pn1", &pn1, "pn1/D");
790 
791  APDtrees[i]->SetBranchAddress("event", &event);
792  APDtrees[i]->SetBranchAddress("color", &color);
793  APDtrees[i]->SetBranchAddress("iphi", &iphi);
794  APDtrees[i]->SetBranchAddress("ieta", &ieta);
795  APDtrees[i]->SetBranchAddress("side", &side);
796  APDtrees[i]->SetBranchAddress("dccID", &dccID);
797  APDtrees[i]->SetBranchAddress("towerID", &towerID);
798  APDtrees[i]->SetBranchAddress("channelID", &channelID);
799  APDtrees[i]->SetBranchAddress("apdAmpl", &apdAmpl);
800  APDtrees[i]->SetBranchAddress("apdTime", &apdTime);
801  if (_saveallevents)
802  APDtrees[i]->SetBranchAddress("adc", adc);
803  APDtrees[i]->SetBranchAddress("flagfit", &flagfit);
804  APDtrees[i]->SetBranchAddress("pn0", &pn0);
805  APDtrees[i]->SetBranchAddress("pn1", &pn1);
806  }
807 
808  for (unsigned int iref = 0; iref < nRefChan; iref++) {
809  for (unsigned int imod = 0; imod < nMod; imod++) {
810  int jmod = modules[imod];
811 
812  stringstream nameref;
813  nameref << "refAPDTree" << imod << "_" << iref;
814 
815  RefAPDtrees[iref][jmod] = new TTree(nameref.str().c_str(), nameref.str().c_str());
816 
817  RefAPDtrees[iref][jmod]->Branch("eventref", &eventref, "eventref/I");
818  RefAPDtrees[iref][jmod]->Branch("colorref", &colorref, "colorref/I");
819  if (iref == 0)
820  RefAPDtrees[iref][jmod]->Branch("apdAmplA", &apdAmplA, "apdAmplA/D");
821  if (iref == 1)
822  RefAPDtrees[iref][jmod]->Branch("apdAmplB", &apdAmplB, "apdAmplB/D");
823 
824  RefAPDtrees[iref][jmod]->SetBranchAddress("eventref", &eventref);
825  RefAPDtrees[iref][jmod]->SetBranchAddress("colorref", &colorref);
826  if (iref == 0)
827  RefAPDtrees[iref][jmod]->SetBranchAddress("apdAmplA", &apdAmplA);
828  if (iref == 1)
829  RefAPDtrees[iref][jmod]->SetBranchAddress("apdAmplB", &apdAmplB);
830  }
831  }
832 
833  assert(colors.size() <= nColor);
834  unsigned int nCol = colors.size();
835 
836  // Declare PN stuff
837  //===================
838 
839  for (unsigned int iM = 0; iM < nMod; iM++) {
840  unsigned int iMod = modules[iM] - 1;
841  for (unsigned int ich = 0; ich < nPNPerMod; ich++) {
842  for (unsigned int icol = 0; icol < nCol; icol++) {
843  PNFirstAnal[iMod][ich][icol] = new TPN(ich);
844  PNAnal[iMod][ich][icol] = new TPN(ich);
845  }
846  }
847  }
848 
849  // Declare function for APD ampl fit
850  //===================================
851 
852  PulseFitWithShape* psfit = new PulseFitWithShape();
853 
854  for (unsigned int iCry = 0; iCry < nCrys; iCry++) {
855  for (unsigned int icol = 0; icol < nCol; icol++) {
856  // Declare APD stuff
857  //===================
858 
859  APDFirstAnal[iCry][icol] = new TAPD();
860  IsThereDataADC[iCry][icol] = 1;
861  stringstream cut;
862  cut << "color==" << colors.at(icol);
863  if (ADCtrees[iCry]->GetEntries(cut.str().c_str()) < 10)
864  IsThereDataADC[iCry][icol] = 0;
865  }
866 
867  unsigned int iMod = iModule[iCry] - 1;
868  assert(iMod <= nMod);
869 
870  if (isSPRFine)
871  psfit->init(_nsamples, _firstsample, _lastsample, _niter, nSamplesShapes, shapesVec, _noise);
872 
873  // Loop on events
874  //================
875 
876  Long64_t nbytes = 0, nb = 0;
877  for (Long64_t jentry = 0; jentry < ADCtrees[iCry]->GetEntriesFast(); jentry++) { // Loop on events
878  nb = ADCtrees[iCry]->GetEntry(jentry);
879  nbytes += nb;
880 
881  flagfit = 1;
882  apdAmpl = 0.0;
883  apdTime = 0.0;
884  ieta = eta;
885  iphi = phi;
886 
887  // Get back color
888 
889  unsigned int iCol = 0;
890  for (unsigned int i = 0; i < nCol; i++) {
891  if (color == colors[i]) {
892  iCol = i;
893  i = colors.size();
894  }
895  }
896 
897  // Amplitude calculation
898 
899  APDPulse->setPulse(adc);
900  adcNoPed = APDPulse->getAdcWithoutPedestal();
901 
902  if (isSPRFine && APDPulse->isPulseOK()) {
903  psfit->doFit(&adcNoPed[0]);
904  apdAmpl = psfit->getAmpl();
905  apdTime = psfit->getTime();
906 
907  } else {
908  apdAmpl = 0;
909  apdTime = 0;
910  flagfit = 0;
911  }
912 
913  if (_debug >= 1)
914  cout << "-- debug test -- endJob -- apdAmpl:" << apdAmpl << " apdTime:" << apdTime << endl;
915 
916  double pnmean;
917  if (pn0 < 10 && pn1 > 10) {
918  pnmean = pn1;
919  } else if (pn1 < 10 && pn0 > 10) {
920  pnmean = pn0;
921  } else
922  pnmean = 0.5 * (pn0 + pn1);
923 
924  if (_debug >= 1)
925  cout << "-- debug test -- endJob -- pnMean:" << pnmean << endl;
926 
927  // Fill PN stuff
928  //===============
929 
930  if (firstChanMod[iMod] == iCry && IsThereDataADC[iCry][iCol] == 1) {
931  for (unsigned int ichan = 0; ichan < nPNPerMod; ichan++) {
932  PNFirstAnal[iMod][ichan][iCol]->addEntry(pnmean, pn0, pn1);
933  }
934  }
935 
936  // Fill APD stuff
937  //================
938 
939  if (apdAmpl != 0.0)
940  APDFirstAnal[iCry][iCol]->addEntry(apdAmpl, pnmean, pn0, pn1, apdTime);
941  if (_debug >= 1)
942  cout << "-- debug test -- endJob -- filling APDTree" << endl;
943 
944  APDtrees[iCry]->Fill();
945 
946  // Fill reference trees
947  //=====================
948 
949  if (apdRefMap[0][iMod + 1] == iCry || apdRefMap[1][iMod + 1] == iCry) {
950  apdAmplA = 0.0;
951  apdAmplB = 0.0;
952  eventref = event;
953  colorref = color;
954 
955  for (unsigned int ir = 0; ir < nRefChan; ir++) {
956  if (_debug >= 1)
957  cout << "-- debug test -- ir:" << ir << " tt:" << towerID << " refmap:" << apdRefMap[ir][iMod + 1]
958  << " iCry:" << iCry << endl;
959 
960  if (apdRefMap[ir][iMod + 1] == iCry) {
961  if (_debug >= 1)
962  cout << "-- debug test -- cut passed " << endl;
963  if (ir == 0)
964  apdAmplA = apdAmpl;
965  else if (ir == 1)
966  apdAmplB = apdAmpl;
967  if (_debug >= 1)
968  cout << "-- debug test -- apdAmplA=" << apdAmplA << endl;
969  if (_debug >= 1)
970  cout << "-- debug test -- apdAmplB=" << apdAmplB << endl;
971  if (_debug >= 1)
972  cout << "-- debug test -- color=" << color << ", event:" << event << ", ir:" << ir
973  << " tt-1:" << towerID - 1 << endl;
974 
975  RefAPDtrees[ir][iMod + 1]->Fill();
976 
977  if (_debug >= 1)
978  cout << "-- debug test -- tree filled" << event << endl;
979  }
980  }
981  }
982  }
983  }
984 
985  delete psfit;
986 
987  ADCFile->Close();
988 
989  if (_debug == 1)
990  cout << "-- debug test -- endJob -- after apdAmpl Loop" << endl;
991 
992  // Remove temporary file
993  //=======================
994 
995  stringstream del;
996  del << "rm " << ADCfile;
997  system(del.str().c_str());
998 
999  // Create output trees
1000  //=====================
1001 
1002  resFile = new TFile(resfile.c_str(), "RECREATE");
1003 
1004  for (unsigned int iColor = 0; iColor < nCol; iColor++) {
1005  stringstream nametree;
1006  nametree << "APDCol" << colors.at(iColor);
1007  stringstream nametree2;
1008  nametree2 << "PNCol" << colors.at(iColor);
1009 
1010  restrees[iColor] = new TTree(nametree.str().c_str(), nametree.str().c_str());
1011  respntrees[iColor] = new TTree(nametree2.str().c_str(), nametree2.str().c_str());
1012 
1013  restrees[iColor]->Branch("iphi", &iphi, "iphi/I");
1014  restrees[iColor]->Branch("ieta", &ieta, "ieta/I");
1015  restrees[iColor]->Branch("side", &side, "side/I");
1016  restrees[iColor]->Branch("dccID", &dccID, "dccID/I");
1017  restrees[iColor]->Branch("moduleID", &moduleID, "moduleID/I");
1018  restrees[iColor]->Branch("towerID", &towerID, "towerID/I");
1019  restrees[iColor]->Branch("channelID", &channelID, "channelID/I");
1020  restrees[iColor]->Branch("APD", &APD, "APD[6]/D");
1021  restrees[iColor]->Branch("Time", &Time, "Time[6]/D");
1022  restrees[iColor]->Branch("APDoPN", &APDoPN, "APDoPN[6]/D");
1023  restrees[iColor]->Branch("APDoPNA", &APDoPNA, "APDoPNA[6]/D");
1024  restrees[iColor]->Branch("APDoPNB", &APDoPNB, "APDoPNB[6]/D");
1025  restrees[iColor]->Branch("APDoAPDA", &APDoAPDA, "APDoAPDA[6]/D");
1026  restrees[iColor]->Branch("APDoAPDB", &APDoAPDB, "APDoAPDB[6]/D");
1027  restrees[iColor]->Branch("ShapeCor", &ShapeCor, "ShapeCor/D");
1028  restrees[iColor]->Branch("flag", &flag, "flag/I");
1029 
1030  respntrees[iColor]->Branch("moduleID", &moduleID, "moduleID/I");
1031  respntrees[iColor]->Branch("pnID", &pnID, "pnID/I");
1032  respntrees[iColor]->Branch("PN", &PN, "PN[6]/D");
1033  respntrees[iColor]->Branch("PNoPN", &PNoPN, "PNoPN[6]/D");
1034  respntrees[iColor]->Branch("PNoPNA", &PNoPNA, "PNoPNA[6]/D");
1035  respntrees[iColor]->Branch("PNoPNB", &PNoPNB, "PNoPNB[6]/D");
1036 
1037  restrees[iColor]->SetBranchAddress("iphi", &iphi);
1038  restrees[iColor]->SetBranchAddress("ieta", &ieta);
1039  restrees[iColor]->SetBranchAddress("dccID", &dccID);
1040  restrees[iColor]->SetBranchAddress("moduleID", &moduleID);
1041  restrees[iColor]->SetBranchAddress("towerID", &towerID);
1042  restrees[iColor]->SetBranchAddress("channelID", &channelID);
1043  restrees[iColor]->SetBranchAddress("APD", APD);
1044  restrees[iColor]->SetBranchAddress("Time", Time);
1045  restrees[iColor]->SetBranchAddress("APDoPN", APDoPN);
1046  restrees[iColor]->SetBranchAddress("APDoPNA", APDoPNA);
1047  restrees[iColor]->SetBranchAddress("APDoPNB", APDoPNB);
1048  restrees[iColor]->SetBranchAddress("APDoAPDA", APDoAPDA);
1049  restrees[iColor]->SetBranchAddress("APDoAPDB", APDoAPDB);
1050  restrees[iColor]->SetBranchAddress("ShapeCor", &ShapeCor);
1051  restrees[iColor]->SetBranchAddress("flag", &flag);
1052 
1053  respntrees[iColor]->SetBranchAddress("moduleID", &moduleID);
1054  respntrees[iColor]->SetBranchAddress("pnID", &pnID);
1055  respntrees[iColor]->SetBranchAddress("PN", PN);
1056  respntrees[iColor]->SetBranchAddress("PNoPN", PNoPN);
1057  respntrees[iColor]->SetBranchAddress("PNoPNA", PNoPNA);
1058  respntrees[iColor]->SetBranchAddress("PNoPNB", PNoPNB);
1059  }
1060 
1061  // Set Cuts for PN stuff
1062  //=======================
1063 
1064  for (unsigned int iM = 0; iM < nMod; iM++) {
1065  unsigned int iMod = modules[iM] - 1;
1066 
1067  for (unsigned int ich = 0; ich < nPNPerMod; ich++) {
1068  for (unsigned int icol = 0; icol < nCol; icol++) {
1069  PNAnal[iMod][ich][icol]->setPNCut(PNFirstAnal[iMod][ich][icol]->getPN().at(0),
1070  PNFirstAnal[iMod][ich][icol]->getPN().at(1));
1071  }
1072  }
1073  }
1074 
1075  // Build ref trees indexes
1076  //========================
1077 
1078  for (unsigned int imod = 0; imod < nMod; imod++) {
1079  int jmod = modules[imod];
1080  if (RefAPDtrees[0][jmod]->GetEntries() != 0 && RefAPDtrees[1][jmod]->GetEntries() != 0) {
1081  RefAPDtrees[0][jmod]->BuildIndex("eventref");
1082  RefAPDtrees[1][jmod]->BuildIndex("eventref");
1083  }
1084  }
1085 
1086  // Final loop on crystals
1087  //=======================
1088 
1089  for (unsigned int iCry = 0; iCry < nCrys; iCry++) {
1090  unsigned int iMod = iModule[iCry] - 1;
1091 
1092  // Set cuts on APD stuff
1093  //=======================
1094 
1095  for (unsigned int iCol = 0; iCol < nCol; iCol++) {
1096  std::vector<double> lowcut;
1097  std::vector<double> highcut;
1098  double cutMin;
1099  double cutMax;
1100 
1101  cutMin = APDFirstAnal[iCry][iCol]->getAPD().at(0) - 2.0 * APDFirstAnal[iCry][iCol]->getAPD().at(1);
1102  if (cutMin < 0)
1103  cutMin = 0;
1104  cutMax = APDFirstAnal[iCry][iCol]->getAPD().at(0) + 2.0 * APDFirstAnal[iCry][iCol]->getAPD().at(1);
1105 
1106  lowcut.push_back(cutMin);
1107  highcut.push_back(cutMax);
1108 
1109  cutMin = APDFirstAnal[iCry][iCol]->getTime().at(0) - 2.0 * APDFirstAnal[iCry][iCol]->getTime().at(1);
1110  cutMax = APDFirstAnal[iCry][iCol]->getTime().at(0) + 2.0 * APDFirstAnal[iCry][iCol]->getTime().at(1);
1111  lowcut.push_back(cutMin);
1112  highcut.push_back(cutMax);
1113 
1114  APDAnal[iCry][iCol] = new TAPD();
1115  APDAnal[iCry][iCol]->setAPDCut(APDFirstAnal[iCry][iCol]->getAPD().at(0),
1116  APDFirstAnal[iCry][iCol]->getAPD().at(1));
1117  APDAnal[iCry][iCol]->setAPDoPNCut(APDFirstAnal[iCry][iCol]->getAPDoPN().at(0),
1118  APDFirstAnal[iCry][iCol]->getAPDoPN().at(1));
1119  APDAnal[iCry][iCol]->setAPDoPN0Cut(APDFirstAnal[iCry][iCol]->getAPDoPN0().at(0),
1120  APDFirstAnal[iCry][iCol]->getAPDoPN0().at(1));
1121  APDAnal[iCry][iCol]->setAPDoPN1Cut(APDFirstAnal[iCry][iCol]->getAPDoPN1().at(0),
1122  APDFirstAnal[iCry][iCol]->getAPDoPN1().at(1));
1123  APDAnal[iCry][iCol]->setTimeCut(APDFirstAnal[iCry][iCol]->getTime().at(0),
1124  APDFirstAnal[iCry][iCol]->getTime().at(1));
1125 
1126  APDAnal[iCry][iCol]->set2DAPDoAPD0Cut(lowcut, highcut);
1127  APDAnal[iCry][iCol]->set2DAPDoAPD1Cut(lowcut, highcut);
1128  }
1129 
1130  // Final loop on events
1131  //=======================
1132 
1133  Long64_t nbytes = 0, nb = 0;
1134  for (Long64_t jentry = 0; jentry < APDtrees[iCry]->GetEntriesFast(); jentry++) {
1135  nb = APDtrees[iCry]->GetEntry(jentry);
1136  nbytes += nb;
1137 
1138  double pnmean;
1139  if (pn0 < 10 && pn1 > 10) {
1140  pnmean = pn1;
1141  } else if (pn1 < 10 && pn0 > 10) {
1142  pnmean = pn0;
1143  } else
1144  pnmean = 0.5 * (pn0 + pn1);
1145 
1146  // Get back color
1147  //===============
1148 
1149  unsigned int iCol = 0;
1150  for (unsigned int i = 0; i < nCol; i++) {
1151  if (color == colors[i]) {
1152  iCol = i;
1153  i = colors.size();
1154  }
1155  }
1156 
1157  // Fill PN stuff
1158  //===============
1159 
1160  if (firstChanMod[iMod] == iCry && IsThereDataADC[iCry][iCol] == 1) {
1161  for (unsigned int ichan = 0; ichan < nPNPerMod; ichan++) {
1162  PNAnal[iMod][ichan][iCol]->addEntry(pnmean, pn0, pn1);
1163  }
1164  }
1165 
1166  // Get ref amplitudes
1167  //===================
1168 
1169  if (_debug >= 1)
1170  cout << "-- debug test -- LastLoop event:" << event << " apdAmpl:" << apdAmpl << endl;
1171  apdAmplA = 0.0;
1172  apdAmplB = 0.0;
1173 
1174  for (unsigned int iRef = 0; iRef < nRefChan; iRef++) {
1175  RefAPDtrees[iRef][iMod + 1]->GetEntryWithIndex(event);
1176  }
1177 
1178  if (_debug == 1)
1179  cout << "-- debug test -- LastLoop apdAmplA:" << apdAmplA << " apdAmplB:" << apdAmplB << ", event:" << event
1180  << ", eventref:" << eventref << endl;
1181 
1182  // Fill APD stuff
1183  //===============
1184 
1185  APDAnal[iCry][iCol]->addEntry(apdAmpl, pnmean, pn0, pn1, apdTime, apdAmplA, apdAmplB);
1186  }
1187 
1188  moduleID = iMod + 1;
1189 
1190  if (moduleID >= 20)
1191  moduleID -= 2; // Trick to fix endcap specificity
1192 
1193  // Get final results for APD
1194  //===========================
1195 
1196  for (unsigned int iColor = 0; iColor < nCol; iColor++) {
1197  std::vector<double> apdvec = APDAnal[iCry][iColor]->getAPD();
1198  std::vector<double> apdpnvec = APDAnal[iCry][iColor]->getAPDoPN();
1199  std::vector<double> apdpn0vec = APDAnal[iCry][iColor]->getAPDoPN0();
1200  std::vector<double> apdpn1vec = APDAnal[iCry][iColor]->getAPDoPN1();
1201  std::vector<double> timevec = APDAnal[iCry][iColor]->getTime();
1202  std::vector<double> apdapd0vec = APDAnal[iCry][iColor]->getAPDoAPD0();
1203  std::vector<double> apdapd1vec = APDAnal[iCry][iColor]->getAPDoAPD1();
1204 
1205  for (unsigned int i = 0; i < apdvec.size(); i++) {
1206  APD[i] = apdvec.at(i);
1207  APDoPN[i] = apdpnvec.at(i);
1208  APDoPNA[i] = apdpn0vec.at(i);
1209  APDoPNB[i] = apdpn1vec.at(i);
1210  APDoAPDA[i] = apdapd0vec.at(i);
1211  APDoAPDB[i] = apdapd1vec.at(i);
1212  Time[i] = timevec.at(i);
1213  ShapeCor = shapeCorrection;
1214  }
1215 
1216  // Fill APD results trees
1217  //========================
1218 
1219  iphi = iPhi[iCry];
1220  ieta = iEta[iCry];
1221  dccID = idccID[iCry];
1222  towerID = iTowerID[iCry];
1223  side = iside[iCry];
1224  channelID = iChannelID[iCry];
1225 
1226  if (!wasGainOK[iCry] || !wasTimingOK[iCry] || IsThereDataADC[iCry][iColor] == 0)
1227  flag = 1;
1228  else
1229  flag = 0;
1230 
1231  if (_debug >= 1)
1232  cout << "-- debug test -- endJob -- APD[0]" << APD[0] << " APDoPN[0] " << APDoPN[0] << " APDoAPDA[0] "
1233  << APDoAPDA[0] << " flag " << flag << endl;
1234  restrees[iColor]->Fill();
1235  }
1236  }
1237 
1238  // Get final results for PN
1239  //==========================
1240 
1241  for (unsigned int iM = 0; iM < nMod; iM++) {
1242  unsigned int iMod = modules[iM] - 1;
1243 
1244  side = iside[firstChanMod[iMod]];
1245 
1246  for (unsigned int ch = 0; ch < nPNPerMod; ch++) {
1247  pnID = ch;
1248  moduleID = iMod + 1;
1249 
1250  if (moduleID >= 20)
1251  moduleID -= 2; // Trick to fix endcap specificity
1252 
1253  for (unsigned int iColor = 0; iColor < nCol; iColor++) {
1254  std::vector<double> pnvec = PNAnal[iMod][ch][iColor]->getPN();
1255  std::vector<double> pnopnvec = PNAnal[iMod][ch][iColor]->getPNoPN();
1256  std::vector<double> pnopn0vec = PNAnal[iMod][ch][iColor]->getPNoPN0();
1257  std::vector<double> pnopn1vec = PNAnal[iMod][ch][iColor]->getPNoPN1();
1258 
1259  for (unsigned int i = 0; i < pnvec.size(); i++) {
1260  PN[i] = pnvec.at(i);
1261  PNoPN[i] = pnopnvec.at(i);
1262  PNoPNA[i] = pnopn0vec.at(i);
1263  PNoPNB[i] = pnopn1vec.at(i);
1264  }
1265 
1266  if (_debug >= 1)
1267  cout << "-- debug test -- endJob -- filling pn results'tree: PN[0]:" << PN[0] << " iModule:" << iMod
1268  << " iColor:" << iColor << " ch:" << ch << endl;
1269 
1270  // Fill PN results trees
1271  //========================
1272 
1273  respntrees[iColor]->Fill();
1274  }
1275  }
1276  }
1277 
1278  // Remove temporary files
1279  //========================
1280 
1281  if (!_saveallevents) {
1282  APDFile->Close();
1283  stringstream del2;
1284  del2 << "rm " << APDfile;
1285  system(del2.str().c_str());
1286 
1287  } else {
1288  APDFile->cd();
1289  APDtrees[0]->Write();
1290  APDFile->Close();
1291  resFile->cd();
1292  }
1293 
1294  // Save results
1295  //===============
1296 
1297  for (unsigned int i = 0; i < nCol; i++) {
1298  restrees[i]->Write();
1299  respntrees[i]->Write();
1300  }
1301 
1302  resFile->Close();
1303 
1304  cout << "\t+=+ .................................................. done +=+" << endl;
1305  cout << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+" << endl;
1306 }

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, cms::cuda::assert(), 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().

◆ getShapes()

bool EcalLaserAnalyzer2::getShapes ( )
private

Definition at line 1309 of file EcalLaserAnalyzer2.cc.

1309  {
1310  //========================================================================
1311 
1312  // Get Pulse From Matacq Analysis:
1313  //================================
1314 
1315  bool IsMatacqOK = false;
1316 
1317  int doesMatFileExist = 0;
1318  int doesMatShapeExist = 0;
1319  FILE* test2;
1320  TProfile* laserShape = nullptr;
1321  test2 = fopen(matfile.c_str(), "r");
1322  if (test2)
1323  doesMatFileExist = 1;
1324 
1325  TFile* MatShapeFile;
1326  if (doesMatFileExist == 1) {
1327  MatShapeFile = new TFile(matfile.c_str());
1328  laserShape = (TProfile*)MatShapeFile->Get("shapeLaser");
1329  if (laserShape) {
1330  doesMatShapeExist = 1;
1331  double y = laserShape->Integral("w");
1332  if (y != 0)
1333  laserShape->Scale(1.0 / y);
1334  }
1335  } else {
1336  cout << " ERROR! Matacq shape file not found !" << endl;
1337  }
1338  if (doesMatShapeExist)
1339  IsMatacqOK = true;
1340 
1341  // Get SPR from the average elec shape in SM6:
1342  //============================================
1343 
1344  int doesElecFileExist = 0;
1345  FILE* test;
1346  test = fopen(elecfile_.c_str(), "r");
1347  if (test)
1348  doesElecFileExist = 1;
1349 
1350  TFile* ElecShapesFile;
1351  TH1D* elecShape = nullptr;
1352 
1353  if (doesElecFileExist == 1) {
1354  ElecShapesFile = new TFile(elecfile_.c_str());
1355  stringstream name;
1356  name << "MeanElecShape";
1357  elecShape = (TH1D*)ElecShapesFile->Get(name.str().c_str());
1358  if (elecShape && doesMatShapeExist == 1) {
1359  double x = elecShape->GetMaximum();
1360  if (x != 0)
1361  elecShape->Scale(1.0 / x);
1362  isSPRFine = true;
1363  } else {
1364  isSPRFine = false;
1365  }
1366 
1367  } else {
1368  cout << " ERROR! Elec shape file not found !" << endl;
1369  }
1370 
1371  if (IsMatacqOK) {
1372  ShapeFile = new TFile(shapefile.c_str(), "RECREATE");
1373 
1374  unsigned int nBins = int(laserShape->GetEntries());
1375  assert(nSamplesShapes == nBins);
1376  double elec_jj, laser_iiMinusjj;
1377  double sum_jj;
1378 
1379  if (isSPRFine == true) {
1380  unsigned int nBins2 = int(elecShape->GetNbinsX());
1381 
1382  if (nBins2 < nBins) {
1383  cout << "EcalLaserAnalyzer2::getShapes: wrong configuration of the shapes' number of bins" << std::endl;
1384  isSPRFine = false;
1385  }
1386  assert(nSamplesShapes == nBins2);
1387 
1388  stringstream name;
1389  name << "PulseShape";
1390 
1391  PulseShape = new TProfile(name.str().c_str(), name.str().c_str(), nBins, -0.5, double(nBins) - 0.5);
1392 
1393  // shift shapes to have max close to the real APD max
1394 
1395  for (int ii = 0; ii < 50; ii++) {
1396  shapes[ii] = 0.0;
1397  PulseShape->Fill(ii, 0.0);
1398  }
1399 
1400  for (unsigned int ii = 0; ii < nBins - 50; ii++) {
1401  sum_jj = 0.0;
1402  for (unsigned int jj = 0; jj < ii; jj++) {
1403  elec_jj = elecShape->GetBinContent(jj + 1);
1404  laser_iiMinusjj = laserShape->GetBinContent(ii - jj + 1);
1405  sum_jj += elec_jj * laser_iiMinusjj;
1406  }
1407  PulseShape->Fill(ii + 50, sum_jj);
1408  shapes[ii + 50] = sum_jj;
1409  }
1410 
1411  double scale = PulseShape->GetMaximum();
1412  shapeCorrection = scale;
1413 
1414  if (scale != 0) {
1415  PulseShape->Scale(1.0 / scale);
1416  for (unsigned int ii = 0; ii < nBins; ii++) {
1417  shapesVec.push_back(shapes[ii] / scale);
1418  }
1419  }
1420 
1421  if (_saveshapes)
1422  PulseShape->Write();
1423  }
1424  }
1425  ShapeFile->Close();
1426 
1427  if (!_saveshapes) {
1428  stringstream del;
1429  del << "rm " << shapefile;
1430  system(del.str().c_str());
1431  }
1432 
1433  return IsMatacqOK;
1434 }

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

Referenced by beginJob().

◆ setGeomEB()

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.

1437  {
1438  side = MEEBGeom::side(etaG, phiG);
1439 
1440  assert(module >= *min_element(modules.begin(), modules.end()) &&
1441  module <= *max_element(modules.begin(), modules.end()));
1442 
1443  eta = etaG;
1444  phi = phiG;
1445  channelID = 5 * (strip - 1) + xtal - 1;
1446  towerID = tower;
1447 
1448  vector<int> apdRefChan = ME::apdRefChannels(module, lmr);
1449  for (unsigned int iref = 0; iref < nRefChan; iref++) {
1450  if (channelID == apdRefChan[iref] && towerID == apdRefTT && apdRefMap[iref].count(module) == 0) {
1451  apdRefMap[iref][module] = channel;
1452  }
1453  }
1454 
1455  if (isFirstChanModFilled[module - 1] == 0) {
1456  firstChanMod[module - 1] = channel;
1457  isFirstChanModFilled[module - 1] = 1;
1458  }
1459 
1460  iEta[channel] = eta;
1461  iPhi[channel] = phi;
1462  iModule[channel] = module;
1463  iTowerID[channel] = towerID;
1464  iChannelID[channel] = channelID;
1465  idccID[channel] = dccID;
1466  iside[channel] = side;
1467 }

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

Referenced by analyze().

◆ setGeomEE()

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.

1470  {
1471  side = MEEEGeom::side(iX, iY, iZ);
1472 
1473  assert(module >= *min_element(modules.begin(), modules.end()) &&
1474  module <= *max_element(modules.begin(), modules.end()));
1475 
1476  eta = etaG;
1477  phi = phiG;
1478  channelID = ch;
1479  towerID = tower;
1480 
1481  vector<int> apdRefChan = ME::apdRefChannels(module, lmr);
1482  for (unsigned int iref = 0; iref < nRefChan; iref++) {
1483  if (channelID == apdRefChan[iref] && towerID == apdRefTT && apdRefMap[iref].count(module) == 0) {
1484  apdRefMap[iref][module] = channel;
1485  }
1486  }
1487 
1488  if (isFirstChanModFilled[module - 1] == 0) {
1489  firstChanMod[module - 1] = channel;
1490  isFirstChanModFilled[module - 1] = 1;
1491  }
1492 
1493  iEta[channel] = eta;
1494  iPhi[channel] = phi;
1495  iModule[channel] = module;
1496  iTowerID[channel] = towerID;
1497  iChannelID[channel] = channelID;
1498  idccID[channel] = dccID;
1499  iside[channel] = side;
1500 }

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

Referenced by analyze().

Member Data Documentation

◆ _debug

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().

◆ _docorpn

bool EcalLaserAnalyzer2::_docorpn
private

Definition at line 85 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ _ecalPart

std::string EcalLaserAnalyzer2::_ecalPart
private

Definition at line 83 of file EcalLaserAnalyzer2.h.

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

◆ _fedid

int EcalLaserAnalyzer2::_fedid
private

Definition at line 86 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ _firstsample

unsigned int EcalLaserAnalyzer2::_firstsample
private

Definition at line 65 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ _firstsamplePN

unsigned int EcalLaserAnalyzer2::_firstsamplePN
private

Definition at line 71 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ _lastsample

unsigned int EcalLaserAnalyzer2::_lastsample
private

Definition at line 66 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ _lastsamplePN

unsigned int EcalLaserAnalyzer2::_lastsamplePN
private

Definition at line 72 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ _niter

unsigned int EcalLaserAnalyzer2::_niter
private

Definition at line 81 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ _noise

double EcalLaserAnalyzer2::_noise
private

Definition at line 82 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ _nsamples

unsigned int EcalLaserAnalyzer2::_nsamples
private

Definition at line 63 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ _nsamplesPN

unsigned int EcalLaserAnalyzer2::_nsamplesPN
private

Definition at line 69 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ _presample

unsigned int EcalLaserAnalyzer2::_presample
private

Definition at line 64 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ _presamplecut

double EcalLaserAnalyzer2::_presamplecut
private

Definition at line 80 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ _presamplePN

unsigned int EcalLaserAnalyzer2::_presamplePN
private

Definition at line 70 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _qualpercent

double EcalLaserAnalyzer2::_qualpercent
private

Definition at line 88 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ _ratiomaxcutlow

double EcalLaserAnalyzer2::_ratiomaxcutlow
private

Definition at line 79 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _ratiomincuthigh

double EcalLaserAnalyzer2::_ratiomincuthigh
private

Definition at line 78 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _ratiomincutlow

double EcalLaserAnalyzer2::_ratiomincutlow
private

Definition at line 77 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _samplemax

unsigned int EcalLaserAnalyzer2::_samplemax
private

Definition at line 68 of file EcalLaserAnalyzer2.h.

◆ _samplemin

unsigned int EcalLaserAnalyzer2::_samplemin
private

Definition at line 67 of file EcalLaserAnalyzer2.h.

◆ _saveallevents

bool EcalLaserAnalyzer2::_saveallevents
private

Definition at line 87 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ _saveshapes

bool EcalLaserAnalyzer2::_saveshapes
private

Definition at line 84 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

◆ _timingcuthigh

unsigned int EcalLaserAnalyzer2::_timingcuthigh
private

Definition at line 74 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _timingcutlow

unsigned int EcalLaserAnalyzer2::_timingcutlow
private

Definition at line 73 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _timingqualhigh

unsigned int EcalLaserAnalyzer2::_timingqualhigh
private

Definition at line 76 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ _timingquallow

unsigned int EcalLaserAnalyzer2::_timingquallow
private

Definition at line 75 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ adc

double EcalLaserAnalyzer2::adc[10]
private

Definition at line 181 of file EcalLaserAnalyzer2.h.

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

◆ adcC

int EcalLaserAnalyzer2::adcC[1700]
private

Definition at line 226 of file EcalLaserAnalyzer2.h.

◆ ADCfile

std::string EcalLaserAnalyzer2::ADCfile
private

Definition at line 110 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

◆ ADCFile

TFile* EcalLaserAnalyzer2::ADCFile
private

Definition at line 162 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

◆ adcG

int EcalLaserAnalyzer2::adcG[10]
private

Definition at line 182 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ adcMean

double EcalLaserAnalyzer2::adcMean[1700][10]
private

Definition at line 225 of file EcalLaserAnalyzer2.h.

◆ adcNoPed

double* EcalLaserAnalyzer2::adcNoPed
private

Definition at line 195 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ ADCtrees

TTree* EcalLaserAnalyzer2::ADCtrees[1700]
private

Definition at line 163 of file EcalLaserAnalyzer2.h.

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

◆ APD

double EcalLaserAnalyzer2::APD[6]
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ apdAmpl

double EcalLaserAnalyzer2::apdAmpl
private

Definition at line 186 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ apdAmplA

double EcalLaserAnalyzer2::apdAmplA
private

Definition at line 187 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ apdAmplB

double EcalLaserAnalyzer2::apdAmplB
private

Definition at line 188 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDAnal

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

Definition at line 204 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDfile

std::string EcalLaserAnalyzer2::APDfile
private

Definition at line 111 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

◆ APDFile

TFile* EcalLaserAnalyzer2::APDFile
private

Definition at line 165 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDFirstAnal

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

Definition at line 203 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDoAPDA

double EcalLaserAnalyzer2::APDoAPDA[6]
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDoAPDB

double EcalLaserAnalyzer2::APDoAPDB[6]
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDoPN

double EcalLaserAnalyzer2::APDoPN[6]
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDoPNA

double EcalLaserAnalyzer2::APDoPNA[6]
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDoPNB

double EcalLaserAnalyzer2::APDoPNB[6]
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDPulse

TAPDPulse* EcalLaserAnalyzer2::APDPulse
private

Definition at line 91 of file EcalLaserAnalyzer2.h.

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

◆ apdRefMap

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

Definition at line 150 of file EcalLaserAnalyzer2.h.

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

◆ apdTime

double EcalLaserAnalyzer2::apdTime
private

Definition at line 189 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ APDtrees

TTree* EcalLaserAnalyzer2::APDtrees[1700]
private

Definition at line 166 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ channelID

int EcalLaserAnalyzer2::channelID
private

Definition at line 136 of file EcalLaserAnalyzer2.h.

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

◆ channelIteratorEE

int EcalLaserAnalyzer2::channelIteratorEE
private

Definition at line 219 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ channelMapEE

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

Definition at line 147 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

◆ color

int EcalLaserAnalyzer2::color
private

Definition at line 180 of file EcalLaserAnalyzer2.h.

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

◆ colorref

int EcalLaserAnalyzer2::colorref
private

Definition at line 193 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ colors

std::vector<int> EcalLaserAnalyzer2::colors
private

Definition at line 146 of file EcalLaserAnalyzer2.h.

◆ dccID

int EcalLaserAnalyzer2::dccID
private

Definition at line 138 of file EcalLaserAnalyzer2.h.

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

◆ dccMEM

std::vector<int> EcalLaserAnalyzer2::dccMEM
private

Definition at line 148 of file EcalLaserAnalyzer2.h.

◆ Delta01

TMom* EcalLaserAnalyzer2::Delta01
private

Definition at line 94 of file EcalLaserAnalyzer2.h.

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

◆ Delta12

TMom* EcalLaserAnalyzer2::Delta12
private

Definition at line 95 of file EcalLaserAnalyzer2.h.

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

◆ digiCollection_

std::string EcalLaserAnalyzer2::digiCollection_
private

Definition at line 98 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ digiPNCollection_

std::string EcalLaserAnalyzer2::digiPNCollection_
private

Definition at line 101 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ digiProducer_

std::string EcalLaserAnalyzer2::digiProducer_
private

Definition at line 102 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ elecfile_

std::string EcalLaserAnalyzer2::elecfile_
private

Definition at line 99 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and getShapes().

◆ eta

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().

◆ event

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().

◆ eventHeaderCollection_

std::string EcalLaserAnalyzer2::eventHeaderCollection_
private

Definition at line 103 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ eventHeaderProducer_

std::string EcalLaserAnalyzer2::eventHeaderProducer_
private

Definition at line 104 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ eventref

int EcalLaserAnalyzer2::eventref
private

Definition at line 192 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ fedID

int EcalLaserAnalyzer2::fedID
private

Definition at line 137 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ firstChanMod

unsigned int EcalLaserAnalyzer2::firstChanMod[21]
private

Definition at line 231 of file EcalLaserAnalyzer2.h.

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

◆ flag

int EcalLaserAnalyzer2::flag
private

Definition at line 218 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ iChannelID

int EcalLaserAnalyzer2::iChannelID[1700]
private

Definition at line 230 of file EcalLaserAnalyzer2.h.

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

◆ idccID

int EcalLaserAnalyzer2::idccID[1700]
private

Definition at line 230 of file EcalLaserAnalyzer2.h.

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

◆ iEta

int EcalLaserAnalyzer2::iEta[1700]
private

Definition at line 228 of file EcalLaserAnalyzer2.h.

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

◆ iEvent

int EcalLaserAnalyzer2::iEvent
private

Definition at line 59 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ iModule

unsigned int EcalLaserAnalyzer2::iModule[1700]
private

Definition at line 229 of file EcalLaserAnalyzer2.h.

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

◆ iPhi

int EcalLaserAnalyzer2::iPhi[1700]
private

Definition at line 228 of file EcalLaserAnalyzer2.h.

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

◆ isFirstChanModFilled

unsigned int EcalLaserAnalyzer2::isFirstChanModFilled[21]
private

Definition at line 232 of file EcalLaserAnalyzer2.h.

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

◆ isGainOK

bool EcalLaserAnalyzer2::isGainOK
private

Definition at line 243 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ iside

int EcalLaserAnalyzer2::iside[1700]
private

Definition at line 230 of file EcalLaserAnalyzer2.h.

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

◆ IsMatacqOK

bool EcalLaserAnalyzer2::IsMatacqOK
private

Definition at line 126 of file EcalLaserAnalyzer2.h.

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

◆ isMatacqOK

bool EcalLaserAnalyzer2::isMatacqOK
private

Definition at line 213 of file EcalLaserAnalyzer2.h.

◆ isSPRFine

bool EcalLaserAnalyzer2::isSPRFine
private

Definition at line 214 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

◆ IsThereDataADC

int EcalLaserAnalyzer2::IsThereDataADC[1700][nColor]
private

Definition at line 206 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ isTimingOK

bool EcalLaserAnalyzer2::isTimingOK
private

Definition at line 244 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ iTowerID

int EcalLaserAnalyzer2::iTowerID[1700]
private

Definition at line 230 of file EcalLaserAnalyzer2.h.

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

◆ iZ

int EcalLaserAnalyzer2::iZ
private

Definition at line 141 of file EcalLaserAnalyzer2.h.

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

◆ laserEvents

int EcalLaserAnalyzer2::laserEvents
private

Definition at line 144 of file EcalLaserAnalyzer2.h.

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

◆ lightside

int EcalLaserAnalyzer2::lightside
private

Definition at line 140 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ matfile

std::string EcalLaserAnalyzer2::matfile
private

Definition at line 109 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and getShapes().

◆ Mem

TMem* EcalLaserAnalyzer2::Mem
private

Definition at line 93 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ moduleID

int EcalLaserAnalyzer2::moduleID
private

Definition at line 218 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ modules

std::vector<int> EcalLaserAnalyzer2::modules
private

Definition at line 149 of file EcalLaserAnalyzer2.h.

Referenced by SequenceVisitors.NodeVisitor::enter().

◆ nCrys

unsigned int EcalLaserAnalyzer2::nCrys
private

Definition at line 117 of file EcalLaserAnalyzer2.h.

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

◆ nEvtBadGain

int EcalLaserAnalyzer2::nEvtBadGain[1700]
private

Definition at line 236 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

◆ nEvtBadTiming

int EcalLaserAnalyzer2::nEvtBadTiming[1700]
private

Definition at line 237 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

◆ nEvtTot

int EcalLaserAnalyzer2::nEvtTot[1700]
private

Definition at line 238 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and endJob().

◆ nMod

unsigned int EcalLaserAnalyzer2::nMod
private

Definition at line 121 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ nPNPerMod

unsigned int EcalLaserAnalyzer2::nPNPerMod
private

Definition at line 118 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ nRefChan

unsigned int EcalLaserAnalyzer2::nRefChan
private

Definition at line 119 of file EcalLaserAnalyzer2.h.

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

◆ nRefTrees

unsigned int EcalLaserAnalyzer2::nRefTrees
private

Definition at line 120 of file EcalLaserAnalyzer2.h.

◆ nSamplesShapes

unsigned int EcalLaserAnalyzer2::nSamplesShapes
private

Definition at line 124 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

◆ nSides

unsigned int EcalLaserAnalyzer2::nSides
private

Definition at line 122 of file EcalLaserAnalyzer2.h.

◆ phi

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().

◆ pn

double EcalLaserAnalyzer2::pn[50]
private

Definition at line 184 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ PN

double EcalLaserAnalyzer2::PN[6]
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ pn0

double EcalLaserAnalyzer2::pn0
private

Definition at line 183 of file EcalLaserAnalyzer2.h.

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

◆ pn1

double EcalLaserAnalyzer2::pn1
private

Definition at line 183 of file EcalLaserAnalyzer2.h.

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

◆ pnAmpl

double EcalLaserAnalyzer2::pnAmpl
private

Definition at line 190 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ PNAnal

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

Definition at line 200 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ pncorfile_

std::string EcalLaserAnalyzer2::pncorfile_
private

Definition at line 100 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2().

◆ pnCorrector

TPNCor* EcalLaserAnalyzer2::pnCorrector
private

Definition at line 154 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ PNFirstAnal

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

Definition at line 199 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ pnG

int EcalLaserAnalyzer2::pnG[50]
private

Definition at line 185 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ pnID

int EcalLaserAnalyzer2::pnID
private

Definition at line 218 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ pnNoPed

double* EcalLaserAnalyzer2::pnNoPed
private

Definition at line 196 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ PNoPN

double EcalLaserAnalyzer2::PNoPN[6]
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ PNoPNA

double EcalLaserAnalyzer2::PNoPNA[6]
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ PNoPNB

double EcalLaserAnalyzer2::PNoPNB[6]
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ PNPulse

TPNPulse* EcalLaserAnalyzer2::PNPulse
private

Definition at line 92 of file EcalLaserAnalyzer2.h.

Referenced by analyze(), and EcalLaserAnalyzer2().

◆ PulseShape

TProfile* EcalLaserAnalyzer2::PulseShape
private

Definition at line 174 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

◆ RefAPDtrees

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

Definition at line 167 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ resdir_

std::string EcalLaserAnalyzer2::resdir_
private

Definition at line 97 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and EcalLaserAnalyzer2().

◆ resfile

std::string EcalLaserAnalyzer2::resfile
private

Definition at line 112 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and endJob().

◆ resFile

TFile* EcalLaserAnalyzer2::resFile
private

Definition at line 169 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ respntrees

TTree* EcalLaserAnalyzer2::respntrees[nColor]
private

Definition at line 171 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ restrees

TTree* EcalLaserAnalyzer2::restrees[nColor]
private

Definition at line 170 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ runNum

int EcalLaserAnalyzer2::runNum
private

Definition at line 131 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ runType

int EcalLaserAnalyzer2::runType
private

Definition at line 130 of file EcalLaserAnalyzer2.h.

Referenced by analyze().

◆ ShapeCor

double EcalLaserAnalyzer2::ShapeCor
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ shapeCorrection

double EcalLaserAnalyzer2::shapeCorrection
private

Definition at line 212 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

◆ shapefile

std::string EcalLaserAnalyzer2::shapefile
private

Definition at line 108 of file EcalLaserAnalyzer2.h.

Referenced by beginJob(), and getShapes().

◆ ShapeFile

TFile* EcalLaserAnalyzer2::ShapeFile
private

Definition at line 173 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

◆ shapes

double EcalLaserAnalyzer2::shapes[250]
private

Definition at line 211 of file EcalLaserAnalyzer2.h.

Referenced by getShapes().

◆ shapesVec

std::vector<double> EcalLaserAnalyzer2::shapesVec
private

Definition at line 210 of file EcalLaserAnalyzer2.h.

Referenced by endJob(), and getShapes().

◆ side

int EcalLaserAnalyzer2::side
private

Definition at line 139 of file EcalLaserAnalyzer2.h.

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

◆ Time

double EcalLaserAnalyzer2::Time[6]
private

Definition at line 221 of file EcalLaserAnalyzer2.h.

Referenced by endJob().

◆ towerID

int EcalLaserAnalyzer2::towerID
private

Definition at line 135 of file EcalLaserAnalyzer2.h.

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

◆ wasGainOK

bool EcalLaserAnalyzer2::wasGainOK[1700]
private

Definition at line 240 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

◆ wasTimingOK

bool EcalLaserAnalyzer2::wasTimingOK[1700]
private

Definition at line 241 of file EcalLaserAnalyzer2.h.

Referenced by EcalLaserAnalyzer2(), and endJob().

EcalLaserAnalyzer2::pncorfile_
std::string pncorfile_
Definition: EcalLaserAnalyzer2.h:100
EcalLaserAnalyzer2::eta
int eta
Definition: EcalLaserAnalyzer2.h:178
EcalLaserAnalyzer2::iZ
int iZ
Definition: EcalLaserAnalyzer2.h:141
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
NMODEB
#define NMODEB
Definition: EcalLaserAnalyzer.h:37
EcalElectronicsMapping
Definition: EcalElectronicsMapping.h:28
EcalLaserAnalyzer2::dccID
int dccID
Definition: EcalLaserAnalyzer2.h:138
NCRYSEE
#define NCRYSEE
Definition: EcalABAnalyzer.h:33
TMom::addEntry
void addEntry(double val)
Definition: TMom.cc:88
EcalLaserAnalyzer2::APDoAPDA
double APDoAPDA[6]
Definition: EcalLaserAnalyzer2.h:221
TAPDPulse::getAdcWithoutPedestal
double * getAdcWithoutPedestal()
Definition: TAPDPulse.cc:237
DDAxes::y
EcalLaserAnalyzer2::PNoPN
double PNoPN[6]
Definition: EcalLaserAnalyzer2.h:221
EcalLaserAnalyzer2::channelID
int channelID
Definition: EcalLaserAnalyzer2.h:136
TAPD::getAPDoAPD1
std::vector< double > getAPDoAPD1()
Definition: TAPD.cc:241
EcalLaserAnalyzer2::laserEvents
int laserEvents
Definition: EcalLaserAnalyzer2.h:144
mps_fire.i
i
Definition: mps_fire.py:428
EcalLaserAnalyzer2::_firstsample
unsigned int _firstsample
Definition: EcalLaserAnalyzer2.h:65
EcalLaserAnalyzer2::iPhi
int iPhi[1700]
Definition: EcalLaserAnalyzer2.h:228
EcalLaserAnalyzer2::APDtrees
TTree * APDtrees[1700]
Definition: EcalLaserAnalyzer2.h:166
EcalLaserAnalyzer2::_timingcutlow
unsigned int _timingcutlow
Definition: EcalLaserAnalyzer2.h:73
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:80
EcalLaserAnalyzer2::iModule
unsigned int iModule[1700]
Definition: EcalLaserAnalyzer2.h:229
EcalLaserAnalyzer2::iBlue
Definition: EcalLaserAnalyzer2.h:56
EcalLaserAnalyzer2::pnAmpl
double pnAmpl
Definition: EcalLaserAnalyzer2.h:190
EcalLaserAnalyzer2::lightside
int lightside
Definition: EcalLaserAnalyzer2.h:140
TAPD::getAPD
std::vector< double > getAPD()
Definition: TAPD.cc:216
colors
Definition: colors.py:1
edm::Handle::product
T const * product() const
Definition: Handle.h:70
L1EGammaCrystalsEmulatorProducer_cfi.scale
scale
Definition: L1EGammaCrystalsEmulatorProducer_cfi.py:10
TPNFit
Definition: TPNFit.h:8
TPNPulse::setPulse
bool setPulse(double *)
Definition: TPNPulse.cc:45
EcalLaserAnalyzer2::_timingqualhigh
unsigned int _timingqualhigh
Definition: EcalLaserAnalyzer2.h:76
EcalLaserAnalyzer2::pnCorrector
TPNCor * pnCorrector
Definition: EcalLaserAnalyzer2.h:154
EcalLaserAnalyzer2::matfile
std::string matfile
Definition: EcalLaserAnalyzer2.h:109
TkAlMuonSelectors_cfi.cut
cut
Definition: TkAlMuonSelectors_cfi.py:5
EcalLaserAnalyzer2::_presamplecut
double _presamplecut
Definition: EcalLaserAnalyzer2.h:80
edm::DataFrameContainer::const_iterator
boost::transform_iterator< IterHelp, boost::counting_iterator< int > > const_iterator
Definition: DataFrameContainer.h:61
TPN::getPNoPN0
std::vector< double > getPNoPN0()
Definition: TPN.cc:98
TAPDPulse::setPresamples
void setPresamples(int)
Definition: TAPDPulse.cc:251
MEEBGeom::electronic_channel
static int electronic_channel(EBLocalCoord ix, EBLocalCoord iy)
Definition: MEEBGeom.cc:326
EcalLaserAnalyzer2::_lastsample
unsigned int _lastsample
Definition: EcalLaserAnalyzer2.h:66
EcalLaserAnalyzer2::side
int side
Definition: EcalLaserAnalyzer2.h:139
TPN
Definition: TPN.h:8
modules
Definition: ZHLTMatchFilter.cc:17
EBDetId
Definition: EBDetId.h:17
EcalLaserAnalyzer2::_noise
double _noise
Definition: EcalLaserAnalyzer2.h:82
digitizers_cfi.strip
strip
Definition: digitizers_cfi.py:19
EcalLaserAnalyzer2::iEvent
int iEvent
Definition: EcalLaserAnalyzer2.h:59
EBDataFrame
Definition: EBDataFrame.h:11
ME::apdRefChannels
static std::vector< int > apdRefChannels(ME::LMMid ilmmod, ME::LMRid ilmr)
Definition: ME.cc:545
gather_cfg.cout
cout
Definition: gather_cfg.py:144
EcalLaserAnalyzer2::_ratiomincuthigh
double _ratiomincuthigh
Definition: EcalLaserAnalyzer2.h:78
hgcalTowerProducer_cfi.tower
tower
Definition: hgcalTowerProducer_cfi.py:4
EcalLaserAnalyzer2::APDFile
TFile * APDFile
Definition: EcalLaserAnalyzer2.h:165
EcalLaserAnalyzer2::Delta01
TMom * Delta01
Definition: EcalLaserAnalyzer2.h:94
EcalPnDiodeDetId::iPnId
int iPnId() const
get the PnId
Definition: EcalPnDiodeDetId.h:46
NSAMPSHAPES
#define NSAMPSHAPES
Definition: EcalLaserAnalyzer2.h:41
EcalLaserAnalyzer2::Time
double Time[6]
Definition: EcalLaserAnalyzer2.h:221
MEEBGeom::lmr
static int lmr(EBGlobalCoord ieta, EBGlobalCoord iphi)
Definition: MEEBGeom.cc:110
CollectionTags_cfi.EBDigi
EBDigi
Definition: CollectionTags_cfi.py:27
cms::cuda::assert
assert(be >=bs)
edm::SortedCollection
Definition: SortedCollection.h:49
EcalLaserAnalyzer2::APDoPNA
double APDoPNA[6]
Definition: EcalLaserAnalyzer2.h:221
TPN::getPN
std::vector< double > getPN()
Definition: TPN.cc:90
EcalLaserAnalyzer2::adcNoPed
double * adcNoPed
Definition: EcalLaserAnalyzer2.h:195
EcalLaserAnalyzer2::_qualpercent
double _qualpercent
Definition: EcalLaserAnalyzer2.h:88
EcalLaserAnalyzer2::setGeomEB
void setGeomEB(int etaG, int phiG, int module, int tower, int strip, int xtal, int apdRefTT, int channel, int lmr)
Definition: EcalLaserAnalyzer2.cc:1436
NCRYSEB
#define NCRYSEB
Definition: EcalABAnalyzer.h:30
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
EcalLaserAnalyzer2::eventHeaderCollection_
std::string eventHeaderCollection_
Definition: EcalLaserAnalyzer2.h:103
DDAxes::x
EcalLaserAnalyzer2::pn1
double pn1
Definition: EcalLaserAnalyzer2.h:183
EcalLaserAnalyzer2::_debug
int _debug
Definition: EcalLaserAnalyzer2.h:89
EcalLaserAnalyzer2::_nsamplesPN
unsigned int _nsamplesPN
Definition: EcalLaserAnalyzer2.h:69
EcalLaserAnalyzer2::wasTimingOK
bool wasTimingOK[1700]
Definition: EcalLaserAnalyzer2.h:241
EcalPnDiodeDetId
Definition: EcalPnDiodeDetId.h:22
edmStreamStallGrapher.getTime
def getTime(line)
Definition: edmStreamStallGrapher.py:179
EcalLaserAnalyzer2::_fedid
int _fedid
Definition: EcalLaserAnalyzer2.h:86
MEEBGeom::pn
static std::pair< int, int > pn(int ilmmod)
Definition: MEEBGeom.cc:447
EcalLaserAnalyzer2::elecfile_
std::string elecfile_
Definition: EcalLaserAnalyzer2.h:99
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
edm::Handle
Definition: AssociativeIterator.h:50
ME::lmmodFromDcc
static std::vector< ME::LMMid > lmmodFromDcc(ME::DCCid idcc)
Definition: ME.cc:574
EcalLaserAnalyzer2::phi
int phi
Definition: EcalLaserAnalyzer2.h:178
EcalLaserAnalyzer2::runType
int runType
Definition: EcalLaserAnalyzer2.h:130
EcalLaserAnalyzer2::pnG
int pnG[50]
Definition: EcalLaserAnalyzer2.h:185
MEEEGeom::apdRefTower
static int apdRefTower(int ilmr, int ilmmod)
Definition: MEEEGeom.cc:1208
EcalLaserAnalyzer2::runNum
int runNum
Definition: EcalLaserAnalyzer2.h:131
TAPDPulse
Definition: TAPDPulse.h:7
TPN::getPNoPN
std::vector< double > getPNoPN()
Definition: TPN.cc:94
EcalLaserAnalyzer2::APDAnal
TAPD * APDAnal[1700][nColor]
Definition: EcalLaserAnalyzer2.h:204
EcalLaserAnalyzer2::idccID
int idccID[1700]
Definition: EcalLaserAnalyzer2.h:230
EcalElectronicsId::towerId
int towerId() const
get the tower id
Definition: EcalElectronicsId.h:33
LEDCalibrationChannels.iphi
iphi
Definition: LEDCalibrationChannels.py:64
EcalLaserAnalyzer2::resdir_
std::string resdir_
Definition: EcalLaserAnalyzer2.h:97
EcalDCCHeaderBlock::LASER_GAP
Definition: EcalDCCHeaderBlock.h:42
EcalLaserAnalyzer2::_firstsamplePN
unsigned int _firstsamplePN
Definition: EcalLaserAnalyzer2.h:71
EcalLaserAnalyzer2::iEta
int iEta[1700]
Definition: EcalLaserAnalyzer2.h:228
TAPD::getAPDoPN
std::vector< double > getAPDoPN()
Definition: TAPD.cc:220
MEEBGeom::side
static int side(EBGlobalCoord ieta, EBGlobalCoord iphi)
Definition: MEEBGeom.cc:105
EcalLaserAnalyzer2::pnNoPed
double * pnNoPed
Definition: EcalLaserAnalyzer2.h:196
EcalLaserAnalyzer2::APDFirstAnal
TAPD * APDFirstAnal[1700][nColor]
Definition: EcalLaserAnalyzer2.h:203
EcalLaserAnalyzer2::resFile
TFile * resFile
Definition: EcalLaserAnalyzer2.h:169
EcalLaserAnalyzer2::PNFirstAnal
TPN * PNFirstAnal[9][2][nColor]
Definition: EcalLaserAnalyzer2.h:199
EcalLaserAnalyzer2::nEvtTot
int nEvtTot[1700]
Definition: EcalLaserAnalyzer2.h:238
test
Definition: SmallWORMDict.h:13
EcalLaserAnalyzer2::iRed
Definition: EcalLaserAnalyzer2.h:56
EcalLaserAnalyzer2::APDoPN
double APDoPN[6]
Definition: EcalLaserAnalyzer2.h:221
EcalLaserAnalyzer2::pn0
double pn0
Definition: EcalLaserAnalyzer2.h:183
EcalDCCHeaderBlock::LASER_STD
Definition: EcalDCCHeaderBlock.h:26
EcalLaserAnalyzer2::apdAmpl
double apdAmpl
Definition: EcalLaserAnalyzer2.h:186
ctpps_dqm_sourceclient-live_cfg.test
test
Definition: ctpps_dqm_sourceclient-live_cfg.py:7
EcalLaserAnalyzer2::_lastsamplePN
unsigned int _lastsamplePN
Definition: EcalLaserAnalyzer2.h:72
CollectionTags_cfi.EEDigi
EEDigi
Definition: CollectionTags_cfi.py:28
EcalLaserAnalyzer2::_ratiomaxcutlow
double _ratiomaxcutlow
Definition: EcalLaserAnalyzer2.h:79
EcalLaserAnalyzer2::restrees
TTree * restrees[nColor]
Definition: EcalLaserAnalyzer2.h:170
TMem::Mem
int Mem(int, int)
Definition: TMem.cc:41
TAPDPulse::setPulse
bool setPulse(double *)
Definition: TAPDPulse.cc:86
EcalLaserAnalyzer2::APDoPNB
double APDoPNB[6]
Definition: EcalLaserAnalyzer2.h:221
EcalLaserAnalyzer2::isGainOK
bool isGainOK
Definition: EcalLaserAnalyzer2.h:243
EcalElectronicsId::channelId
int channelId() const
so far for EndCap only :
Definition: EcalElectronicsId.cc:71
TAPD::getTime
std::vector< double > getTime()
Definition: TAPD.cc:232
TPNPulse::getAdcWithoutPedestal
double * getAdcWithoutPedestal()
Definition: TPNPulse.cc:89
EcalLaserAnalyzer2::_saveallevents
bool _saveallevents
Definition: EcalLaserAnalyzer2.h:87
EcalLaserAnalyzer2::nSides
unsigned int nSides
Definition: EcalLaserAnalyzer2.h:122
alignCSCRings.corr
dictionary corr
Definition: alignCSCRings.py:124
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:262
EcalLaserAnalyzer2::ADCfile
std::string ADCfile
Definition: EcalLaserAnalyzer2.h:110
EcalLaserAnalyzer2::adcG
int adcG[10]
Definition: EcalLaserAnalyzer2.h:182
EcalLaserAnalyzer2::_timingcuthigh
unsigned int _timingcuthigh
Definition: EcalLaserAnalyzer2.h:74
TAPD::setAPDoPNCut
void setAPDoPNCut(double, double)
Definition: TAPD.cc:142
EcalLaserAnalyzer2::_niter
unsigned int _niter
Definition: EcalLaserAnalyzer2.h:81
edm::ESHandle< EcalElectronicsMapping >
EcalElectronicsMapping::getElectronicsId
EcalElectronicsId getElectronicsId(const DetId &id) const
Get the electronics id for this det id.
Definition: EcalElectronicsMapping.cc:324
EcalLaserAnalyzer2::_timingquallow
unsigned int _timingquallow
Definition: EcalLaserAnalyzer2.h:75
EcalLaserAnalyzer2::_docorpn
bool _docorpn
Definition: EcalLaserAnalyzer2.h:85
submitPVResolutionJobs.count
count
Definition: submitPVResolutionJobs.py:352
EcalLaserAnalyzer2::iTowerID
int iTowerID[1700]
Definition: EcalLaserAnalyzer2.h:230
NSIDES
#define NSIDES
Definition: EcalLaserAnalyzer.h:44
EcalLaserAnalyzer2::PNAnal
TPN * PNAnal[9][2][nColor]
Definition: EcalLaserAnalyzer2.h:200
seedmultiplicitymonitor_newtracking_cfi.nBins
nBins
Definition: seedmultiplicitymonitor_newtracking_cfi.py:8
PulseFitWithShape::getTime
double getTime()
Definition: PulseFitWithShape.h:29
TAPDPulse::isTimingQualOK
bool isTimingQualOK()
Definition: TAPDPulse.cc:145
EcalLaserAnalyzer2::PN
double PN[6]
Definition: EcalLaserAnalyzer2.h:221
EcalLaserAnalyzer2::isSPRFine
bool isSPRFine
Definition: EcalLaserAnalyzer2.h:214
NREFCHAN
#define NREFCHAN
Definition: EcalLaserAnalyzer.h:45
PulseFitWithShape::init
virtual void init(int, int, int, int, int, const std::vector< double > &, double)
Definition: PulseFitWithShape.cc:31
EcalElectronicsId
Ecal readout channel identification [32:20] Unused (so far) [19:13] DCC id [12:6] tower [5:3] strip [...
Definition: EcalElectronicsId.h:18
MEEEGeom::pn
static std::pair< int, int > pn(int dee, int ilmod)
Definition: MEEEGeom.cc:574
EcalLaserAnalyzer2::shapefile
std::string shapefile
Definition: EcalLaserAnalyzer2.h:108
EEDetId
Definition: EEDetId.h:14
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
EcalDCCHeaderBlock::LASER_POWER_SCAN
Definition: EcalDCCHeaderBlock.h:27
EcalDCCHeaderBlock::EcalDCCEventSettings
Definition: EcalDCCHeaderBlock.h:64
EcalLaserAnalyzer2::APDfile
std::string APDfile
Definition: EcalLaserAnalyzer2.h:111
PulseFitWithShape::doFit
virtual double doFit(double *, double *cova=nullptr)
Definition: PulseFitWithShape.cc:60
MEEEGeom::dee
static int dee(SuperCrysCoord iX, SuperCrysCoord iY, int iz)
Definition: MEEEGeom.cc:292
cppFunctionSkipper.exception
exception
Definition: cppFunctionSkipper.py:10
TAPD::setAPDCut
void setAPDCut(double, double)
Definition: TAPD.cc:141
EcalLaserAnalyzer2::RefAPDtrees
TTree * RefAPDtrees[2][21]
Definition: EcalLaserAnalyzer2.h:167
EcalLaserAnalyzer2::nCrys
unsigned int nCrys
Definition: EcalLaserAnalyzer2.h:117
LEDCalibrationChannels.ieta
ieta
Definition: LEDCalibrationChannels.py:63
NPNPERMOD
#define NPNPERMOD
Definition: EcalLaserAnalyzer.h:38
EcalLaserAnalyzer2::nSamplesShapes
unsigned int nSamplesShapes
Definition: EcalLaserAnalyzer2.h:124
TAPD::setAPDoPN1Cut
void setAPDoPN1Cut(double, double)
Definition: TAPD.cc:144
EcalPnDiodeDetId::iDCCId
int iDCCId() const
get the DCCId
Definition: EcalPnDiodeDetId.h:44
TMom
Definition: TMom.h:7
EcalLaserAnalyzer2::getShapes
bool getShapes()
Definition: EcalLaserAnalyzer2.cc:1309
TPNFit::getAmpl
double getAmpl()
Definition: TPNFit.h:31
TAPDPulse::isPulseOK
bool isPulseOK()
Definition: TAPDPulse.cc:162
EcalLaserAnalyzer2::PNoPNB
double PNoPNB[6]
Definition: EcalLaserAnalyzer2.h:221
TPN::addEntry
void addEntry(double, double, double)
Definition: TPN.cc:33
EcalLaserAnalyzer2::respntrees
TTree * respntrees[nColor]
Definition: EcalLaserAnalyzer2.h:171
EBDigiCollection
Definition: EcalDigiCollections.h:56
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:267
EEDigiCollection
Definition: EcalDigiCollections.h:69
EcalLaserAnalyzer2::fedID
int fedID
Definition: EcalLaserAnalyzer2.h:137
TPNCor::getPNCorrectionFactor
double getPNCorrectionFactor(double val0, int gain)
Definition: TPNCor.cc:64
EcalLaserAnalyzer2::event
int event
Definition: EcalLaserAnalyzer2.h:179
EcalLaserAnalyzer2::ShapeCor
double ShapeCor
Definition: EcalLaserAnalyzer2.h:221
EcalLaserAnalyzer2::resfile
std::string resfile
Definition: EcalLaserAnalyzer2.h:112
EcalLaserAnalyzer2::_saveshapes
bool _saveshapes
Definition: EcalLaserAnalyzer2.h:84
EcalLaserAnalyzer2::channelMapEE
std::map< int, int > channelMapEE
Definition: EcalLaserAnalyzer2.h:147
EcalLaserAnalyzer2::Delta12
TMom * Delta12
Definition: EcalLaserAnalyzer2.h:95
EcalLaserAnalyzer2::digiProducer_
std::string digiProducer_
Definition: EcalLaserAnalyzer2.h:102
EcalLaserAnalyzer2::towerID
int towerID
Definition: EcalLaserAnalyzer2.h:135
createfilelist.int
int
Definition: createfilelist.py:10
EcalLaserAnalyzer2::iside
int iside[1700]
Definition: EcalLaserAnalyzer2.h:230
EcalLaserAnalyzer2::_ecalPart
std::string _ecalPart
Definition: EcalLaserAnalyzer2.h:83
EcalLaserAnalyzer2::wasGainOK
bool wasGainOK[1700]
Definition: EcalLaserAnalyzer2.h:240
EcalLaserAnalyzer2::adc
double adc[10]
Definition: EcalLaserAnalyzer2.h:181
TMem
Definition: TMem.h:7
PulseFitWithShape::getAmpl
double getAmpl()
Definition: PulseFitWithShape.h:28
Max
T Max(T a, T b)
Definition: MathUtil.h:44
TAPD::addEntry
void addEntry(double, double, double, double, double, double, double)
Definition: TAPD.cc:42
EcalLaserAnalyzer2::iChannelID
int iChannelID[1700]
Definition: EcalLaserAnalyzer2.h:230
EcalLaserAnalyzer2::ShapeFile
TFile * ShapeFile
Definition: EcalLaserAnalyzer2.h:173
TPNFit::init
void init(int, int, int)
Definition: TPNFit.cc:24
EcalLaserAnalyzer2::eventHeaderProducer_
std::string eventHeaderProducer_
Definition: EcalLaserAnalyzer2.h:104
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:31
get
#define get
TAPD
Definition: TAPD.h:8
EcalLaserAnalyzer2::flag
int flag
Definition: EcalLaserAnalyzer2.h:218
TAPDPulse::getDelta
double getDelta(int, int)
Definition: TAPDPulse.cc:116
EcalLaserAnalyzer2::isFirstChanModFilled
unsigned int isFirstChanModFilled[21]
Definition: EcalLaserAnalyzer2.h:232
EcalLaserAnalyzer2::PulseShape
TProfile * PulseShape
Definition: EcalLaserAnalyzer2.h:174
MEEBGeom::apdRefTower
static int apdRefTower(int ilmmod)
Definition: MEEBGeom.cc:490
TMom::getMean
double getMean()
Definition: TMom.cc:121
EcalMGPASample
Definition: EcalMGPASample.h:22
EcalLaserAnalyzer2::setGeomEE
void setGeomEE(int etaG, int phiG, int iX, int iY, int iZ, int module, int tower, int ch, int apdRefTT, int channel, int lmr)
Definition: EcalLaserAnalyzer2.cc:1469
TMem::isMemRelevant
bool isMemRelevant(int)
Definition: TMem.cc:30
EcalLaserAnalyzer2::moduleID
int moduleID
Definition: EcalLaserAnalyzer2.h:218
EcalLaserAnalyzer2::_ratiomincutlow
double _ratiomincutlow
Definition: EcalLaserAnalyzer2.h:77
EcalLaserAnalyzer2::shapeCorrection
double shapeCorrection
Definition: EcalLaserAnalyzer2.h:212
MEEEGeom::lmmod
static int lmmod(SuperCrysCoord iX, SuperCrysCoord iY)
Definition: MEEEGeom.cc:112
EcalLaserAnalyzer2::colorref
int colorref
Definition: EcalLaserAnalyzer2.h:193
MEEEGeom::lmr
static int lmr(SuperCrysCoord iX, SuperCrysCoord iY, int iz)
Definition: MEEEGeom.cc:254
EcalLaserAnalyzer2::nMod
unsigned int nMod
Definition: EcalLaserAnalyzer2.h:121
TPNFit::doFit
double doFit(int, double *)
Definition: TPNFit.cc:39
EcalLaserAnalyzer2::pnID
int pnID
Definition: EcalLaserAnalyzer2.h:218
EcalLaserAnalyzer2::_nsamples
unsigned int _nsamples
Definition: EcalLaserAnalyzer2.h:63
TPNCor
Definition: TPNCor.h:7
EcalLaserAnalyzer2::_presample
unsigned int _presample
Definition: EcalLaserAnalyzer2.h:64
TAPD::getAPDoPN0
std::vector< double > getAPDoPN0()
Definition: TAPD.cc:224
EcalLaserAnalyzer2::apdTime
double apdTime
Definition: EcalLaserAnalyzer2.h:189
TAPD::setTimeCut
void setTimeCut(double, double)
Definition: TAPD.cc:145
EcalLaserAnalyzer2::APDoAPDB
double APDoAPDB[6]
Definition: EcalLaserAnalyzer2.h:221
EcalMappingRcd
Definition: EcalMappingRcd.h:15
EcalLaserAnalyzer2::shapesVec
std::vector< double > shapesVec
Definition: EcalLaserAnalyzer2.h:210
hgcalPerformanceValidation.df
df
Definition: hgcalPerformanceValidation.py:640
EcalLaserAnalyzer2::apdRefMap
std::map< int, unsigned int > apdRefMap[2]
Definition: EcalLaserAnalyzer2.h:150
EcalDCCHeaderBlock::EcalDCCEventSettings::wavelength
short wavelength
Definition: EcalDCCHeaderBlock.h:67
TAPD::set2DAPDoAPD1Cut
void set2DAPDoAPD1Cut(const std::vector< double > &, const std::vector< double > &)
Definition: TAPD.cc:181
EcalLaserAnalyzer2::pn
double pn[50]
Definition: EcalLaserAnalyzer2.h:184
EcalLaserAnalyzer2::firstChanMod
unsigned int firstChanMod[21]
Definition: EcalLaserAnalyzer2.h:231
MEEBGeom::lmmod
static int lmmod(EBGlobalCoord ieta, EBGlobalCoord iphi)
Definition: MEEBGeom.cc:90
EcalLaserAnalyzer2::isTimingOK
bool isTimingOK
Definition: EcalLaserAnalyzer2.h:244
EcalLaserAnalyzer2::digiPNCollection_
std::string digiPNCollection_
Definition: EcalLaserAnalyzer2.h:101
EcalLaserAnalyzer2::PNPulse
TPNPulse * PNPulse
Definition: EcalLaserAnalyzer2.h:92
findQualityFiles.jj
string jj
Definition: findQualityFiles.py:188
EcalLaserAnalyzer2::IsMatacqOK
bool IsMatacqOK
Definition: EcalLaserAnalyzer2.h:126
EcalDCCHeaderBlock::LASER_DELAY_SCAN
Definition: EcalDCCHeaderBlock.h:28
EcalLaserAnalyzer2::nRefChan
unsigned int nRefChan
Definition: EcalLaserAnalyzer2.h:119
EcalLaserAnalyzer2::APDPulse
TAPDPulse * APDPulse
Definition: EcalLaserAnalyzer2.h:91
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
EcalLaserAnalyzer2::PNoPNA
double PNoPNA[6]
Definition: EcalLaserAnalyzer2.h:221
EcalLaserAnalyzer2::apdAmplB
double apdAmplB
Definition: EcalLaserAnalyzer2.h:188
TPN::setPNCut
void setPNCut(double, double)
Definition: TPN.cc:70
EcalLaserAnalyzer2::channelIteratorEE
int channelIteratorEE
Definition: EcalLaserAnalyzer2.h:219
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
MEEBGeom::localCoord
static XYCoord localCoord(int icr)
Definition: MEEBGeom.cc:142
EcalLaserAnalyzer2::nColor
Definition: EcalLaserAnalyzer2.h:56
TAPD::getAPDoAPD0
std::vector< double > getAPDoAPD0()
Definition: TAPD.cc:236
EcalLaserAnalyzer2::eventref
int eventref
Definition: EcalLaserAnalyzer2.h:192
EcalLaserAnalyzer2::ADCFile
TFile * ADCFile
Definition: EcalLaserAnalyzer2.h:162
EcalLaserAnalyzer2::nPNPerMod
unsigned int nPNPerMod
Definition: EcalLaserAnalyzer2.h:118
EcalLaserAnalyzer2::color
int color
Definition: EcalLaserAnalyzer2.h:180
MEEEGeom::side
static int side(SuperCrysCoord iX, SuperCrysCoord iY, int iz)
Definition: MEEEGeom.cc:1155
EcalElectronicsId::stripId
int stripId() const
get the tower id
Definition: EcalElectronicsId.h:35
EcalLaserAnalyzer2::ADCtrees
TTree * ADCtrees[1700]
Definition: EcalLaserAnalyzer2.h:163
EcalLaserAnalyzer2::shapes
double shapes[250]
Definition: EcalLaserAnalyzer2.h:211
PulseFitWithShape
Definition: PulseFitWithShape.h:10
EcalLaserAnalyzer2::APD
double APD[6]
Definition: EcalLaserAnalyzer2.h:221
EcalLaserAnalyzer2::_presamplePN
unsigned int _presamplePN
Definition: EcalLaserAnalyzer2.h:70
EcalLaserAnalyzer2::IsThereDataADC
int IsThereDataADC[1700][nColor]
Definition: EcalLaserAnalyzer2.h:206
TAPD::setAPDoPN0Cut
void setAPDoPN0Cut(double, double)
Definition: TAPD.cc:143
EcalLaserAnalyzer2::Mem
TMem * Mem
Definition: EcalLaserAnalyzer2.h:93
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
TPN::getPNoPN1
std::vector< double > getPNoPN1()
Definition: TPN.cc:102
EcalLaserAnalyzer2::digiCollection_
std::string digiCollection_
Definition: EcalLaserAnalyzer2.h:98
TAPD::getAPDoPN1
std::vector< double > getAPDoPN1()
Definition: TAPD.cc:228
EcalLaserAnalyzer2::apdAmplA
double apdAmplA
Definition: EcalLaserAnalyzer2.h:187
event
Definition: event.py:1
EEDataFrame
Definition: EEDataFrame.h:12
TPNPulse::getMaxSample
int getMaxSample()
Definition: TPNPulse.cc:70
EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0.cerr
cerr
Definition: EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0.py:8
EcalElectronicsId::xtalId
int xtalId() const
get the channel id
Definition: EcalElectronicsId.h:37
ecalpyutils::hashedIndex
int hashedIndex(int ieta, int iphi)
Definition: EcalPyUtils.cc:38
cuy.ii
ii
Definition: cuy.py:590
TAPD::set2DAPDoAPD0Cut
void set2DAPDoAPD0Cut(const std::vector< double > &, const std::vector< double > &)
Definition: TAPD.cc:173
EcalLaserAnalyzer2::nEvtBadTiming
int nEvtBadTiming[1700]
Definition: EcalLaserAnalyzer2.h:237
EcalLaserAnalyzer2::nEvtBadGain
int nEvtBadGain[1700]
Definition: EcalLaserAnalyzer2.h:236
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
TPNPulse
Definition: TPNPulse.h:7
Generated for CMSSW Reference Manual by doxygen 1.8.16