#include <EcalPerEvtLaserAnalyzer.h>

Inheritance diagram for EcalPerEvtLaserAnalyzer:

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

	EcalPerEvtLaserAnalyzer (const edm::ParameterSet &iConfig)

void	endJob () override

	~EcalPerEvtLaserAnalyzer () 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::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	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 Attributes
unsigned int	_channel

std::string	_ecalPart

int	_fedid

unsigned int	_firstsample

unsigned int	_firstsamplePN

unsigned int	_lastsample

unsigned int	_lastsamplePN

unsigned int	_niter

unsigned int	_nsamples

unsigned int	_nsamplesPN

unsigned int	_presample

unsigned int	_presamplePN

unsigned int	_samplemax

unsigned int	_samplemin

unsigned int	_timingcuthigh

unsigned int	_timingcutlow

unsigned int	_tower

double	adc [10]

std::string	ADCfile

TFile *	ADCFile

int	adcG [10]

TTree *	ADCtrees

std::string	alphafile

double	apdAmpl

std::string	APDfile

TFile *	APDFile

double	apdTime

TTree *	APDtrees

int	channelNumber

int	colMat

int	color

std::vector< int >	colors

int	dccID

std::string	digiCollection_

std::string	digiPNCollection_

std::string	digiProducer_

int	doesRefFileExist

int	eta

int	event

std::string	eventHeaderCollection_

std::string	eventHeaderProducer_

int	evtMat

int	fedID

TTree *	header [2]

int	iEvent

int	IsFileCreated

int	IsHeaderFilled [nColor]

int	IsThereDataADC [nColor]

int	lightside

TFile *	matacqFile

TTree *	matacqTree

unsigned int	nCh

unsigned int	nCrys

unsigned int	nSides

unsigned int	nTT

double	peak

double	peakMat

int	phi

double	pn [50]

double	pn0

double	pn1

double	pnAmpl

std::string	refalphabeta_

std::string	resdir_

std::string	resfile

int	runNum

int	runType

double	tt

double	ttMat

double	ttrig

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 33 of file EcalPerEvtLaserAnalyzer.h.

Member Enumeration Documentation

◆ VarCol

enum EcalPerEvtLaserAnalyzer::VarCol

Enumerator
iBlue
iRed
nColor

Definition at line 42 of file EcalPerEvtLaserAnalyzer.h.

42 { iBlue, iRed, nColor };

Constructor & Destructor Documentation

◆ EcalPerEvtLaserAnalyzer()

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

explicit

Definition at line 46 of file EcalPerEvtLaserAnalyzer.cc.

     : iEvent(0),
  
       // framework parameters with default values
       _nsamples(iConfig.getUntrackedParameter<unsigned int>("nSamples", 10)),
       _presample(iConfig.getUntrackedParameter<unsigned int>("nPresamples", 3)),
       _firstsample(iConfig.getUntrackedParameter<unsigned int>("firstSample", 1)),
       _lastsample(iConfig.getUntrackedParameter<unsigned int>("lastSample", 2)),
       _samplemin(iConfig.getUntrackedParameter<unsigned int>("sampleMin", 3)),
       _samplemax(iConfig.getUntrackedParameter<unsigned int>("sampleMax", 9)),
       _nsamplesPN(iConfig.getUntrackedParameter<unsigned int>("nSamplesPN", 50)),
       _presamplePN(iConfig.getUntrackedParameter<unsigned int>("nPresamplesPN", 6)),
       _firstsamplePN(iConfig.getUntrackedParameter<unsigned int>("firstSamplePN", 7)),
       _lastsamplePN(iConfig.getUntrackedParameter<unsigned int>("lastSamplePN", 8)),
       _timingcutlow(iConfig.getUntrackedParameter<unsigned int>("timingCutLow", 3)),
       _timingcuthigh(iConfig.getUntrackedParameter<unsigned int>("timingCutHigh", 7)),
       _niter(iConfig.getUntrackedParameter<unsigned int>("nIter", 3)),
       _fedid(iConfig.getUntrackedParameter<unsigned int>("fedID", 0)),
       _tower(iConfig.getUntrackedParameter<unsigned int>("tower", 1)),
       _channel(iConfig.getUntrackedParameter<unsigned int>("channel", 1)),
       _ecalPart(iConfig.getUntrackedParameter<std::string>("ecalPart", "EB")),
       nCrys(NCRYSEB),
       nTT(NTTEB),
       nSides(NSIDES),
       IsFileCreated(0),
       nCh(0),
       runType(-1),
       runNum(0),
       dccID(-1),
       lightside(2),
       doesRefFileExist(0),
       ttMat(-1),
       peakMat(-1),
       peak(-1),
       evtMat(-1),
       colMat(-1)
 //========================================================================
  
 {
   //now do what ever initialization is needed
  
   resdir_ = iConfig.getUntrackedParameter<std::string>("resDir");
   refalphabeta_ = iConfig.getUntrackedParameter<std::string>("refAlphaBeta");
  
   digiCollection_ = iConfig.getParameter<std::string>("digiCollection");
   digiPNCollection_ = iConfig.getParameter<std::string>("digiPNCollection");
   digiProducer_ = iConfig.getParameter<std::string>("digiProducer");
  
   eventHeaderCollection_ = iConfig.getParameter<std::string>("eventHeaderCollection");
   eventHeaderProducer_ = iConfig.getParameter<std::string>("eventHeaderProducer");
  
   // Define geometrical constants
  
   if (_ecalPart == "EB") {
     nCrys = NCRYSEB;
     nTT = NTTEB;
   } else {
     nCrys = NCRYSEE;
     nTT = NTTEE;
   }
 }

References _ecalPart, digiCollection_, digiPNCollection_, digiProducer_, eventHeaderCollection_, eventHeaderProducer_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), nCrys, NCRYSEB, NCRYSEE, nTT, NTTEB, NTTEE, refalphabeta_, resdir_, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~EcalPerEvtLaserAnalyzer()

EcalPerEvtLaserAnalyzer::~EcalPerEvtLaserAnalyzer ( )

override

Definition at line 110 of file EcalPerEvtLaserAnalyzer.cc.

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

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 164 of file EcalPerEvtLaserAnalyzer.cc.

                                                                                {
   //========================================================================
  
   ++iEvent;
  
   // retrieving DCC header
   edm::Handle<EcalRawDataCollection> pDCCHeader;
   const EcalRawDataCollection* DCCHeader = nullptr;
   try {
     e.getByLabel(eventHeaderProducer_, eventHeaderCollection_, pDCCHeader);
     DCCHeader = pDCCHeader.product();
   } catch (std::exception& ex) {
     std::cerr << "Error! can't get the product  retrieving DCC header" << eventHeaderCollection_.c_str() << std::endl;
   }
  
   // retrieving crystal data from Event
   edm::Handle<EBDigiCollection> pEBDigi;
   const EBDigiCollection* EBDigi = nullptr;
   edm::Handle<EEDigiCollection> pEEDigi;
   const EEDigiCollection* EEDigi = nullptr;
  
   if (_ecalPart == "EB") {
     try {
       e.getByLabel(digiProducer_, digiCollection_, pEBDigi);
       EBDigi = pEBDigi.product();
     } catch (std::exception& ex) {
       std::cerr << "Error! can't get the product retrieving EB crystal data " << digiCollection_.c_str() << std::endl;
     }
   } else {
     try {
       e.getByLabel(digiProducer_, digiCollection_, pEEDigi);
       EEDigi = pEEDigi.product();
     } catch (std::exception& ex) {
       std::cerr << "Error! can't get the product retrieving EE crystal data " << digiCollection_.c_str() << std::endl;
     }
   }
  
   // retrieving crystal PN diodes from Event
   edm::Handle<EcalPnDiodeDigiCollection> pPNDigi;
   const EcalPnDiodeDigiCollection* PNDigi = nullptr;
   try {
     e.getByLabel(digiProducer_, digiPNCollection_, pPNDigi);
     PNDigi = pPNDigi.product();
   } catch (std::exception& ex) {
     std::cerr << "Error! can't get the product " << digiPNCollection_.c_str() << std::endl;
   }
  
   // retrieving electronics mapping
   edm::ESHandle<EcalElectronicsMapping> ecalmapping;
   const EcalElectronicsMapping* TheMapping = nullptr;
   try {
     c.get<EcalMappingRcd>().get(ecalmapping);
     TheMapping = ecalmapping.product();
   } catch (std::exception& ex) {
     std::cerr << "Error! can't get the product EcalMappingRcd" << std::endl;
   }
  
   // ====================================
   // Decode Basic DCCHeader Information
   // ====================================
  
   for (EcalRawDataCollection::const_iterator headerItr = DCCHeader->begin(); headerItr != DCCHeader->end();
        ++headerItr) {
     // Get run type and run number
  
     int fed = headerItr->fedId();
  
     if (fed != _fedid && _fedid != -999)
       continue;
  
     runType = headerItr->getRunType();
     runNum = headerItr->getRunNumber();
     event = headerItr->getLV1();
     dccID = headerItr->getDccInTCCCommand();
     fedID = headerItr->fedId();
  
     // take event only if the fed corresponds to the DCC in TCC
     if (600 + dccID != fedID)
       continue;
  
     // Cut on runType
     if (runType != EcalDCCHeaderBlock::LASER_STD && runType != EcalDCCHeaderBlock::LASER_GAP)
       return;
  
     // Define output results files' names
  
     if (IsFileCreated == 0) {
       stringstream namefile2;
       namefile2 << resdir_ << "/APDAmpl_Run" << runNum << "_" << _fedid << "_" << _tower << "_" << _channel << ".root";
       resfile = namefile2.str();
  
       // Get Matacq ttrig
  
       stringstream namefile;
       namefile << resdir_ << "/Matacq-Run" << runNum << ".root";
  
       doesRefFileExist = 0;
  
       FILE* test;
       test = fopen(namefile.str().c_str(), "r");
       if (test)
         doesRefFileExist = 1;
  
       if (doesRefFileExist == 1) {
         matacqFile = new TFile((namefile.str().c_str()));
         matacqTree = (TTree*)matacqFile->Get("MatacqShape");
  
         matacqTree->SetBranchAddress("event", &evtMat);
         matacqTree->SetBranchAddress("color", &colMat);
         matacqTree->SetBranchAddress("peak", &peakMat);
         matacqTree->SetBranchAddress("ttrig", &ttMat);
       }
  
       IsFileCreated = 1;
     }
  
     // Retrieve laser color and event number
  
     EcalDCCHeaderBlock::EcalDCCEventSettings settings = headerItr->getEventSettings();
     int color = settings.wavelength;
  
     vector<int>::iterator iter = find(colors.begin(), colors.end(), color);
     if (iter == colors.end()) {
       colors.push_back(color);
       cout << " new color found " << color << " " << colors.size() << endl;
     }
   }
  
   // cut on fedID
  
   if (fedID != _fedid && _fedid != -999)
     return;
  
   // ======================
   // Decode PN Information
   // ======================
  
   TPNFit* pnfit = new TPNFit();
   pnfit->init(_nsamplesPN, _firstsamplePN, _lastsamplePN);
  
   double chi2pn = 0;
   double ypnrange[50];
   double dsum = 0.;
   double dsum1 = 0.;
   double bl = 0.;
   double bl1 = 0.;
   double val_max = 0.;
   unsigned int samplemax = 0;
   unsigned int k;
  
   std::vector<double> allPNAmpl;
  
   for (EcalPnDiodeDigiCollection::const_iterator pnItr = PNDigi->begin(); pnItr != PNDigi->end();
        ++pnItr) {  // Loop on PNs
  
     for (int samId = 0; samId < (*pnItr).size(); samId++) {  // Loop on PN samples
       pn[samId] = (*pnItr).sample(samId).adc();
     }
  
     for (dsum = 0., k = 0; k < _presamplePN; k++) {
       dsum += pn[k];
     }
     bl = dsum / ((double)_presamplePN);
  
     for (val_max = 0., k = 0; k < _nsamplesPN; k++) {
       ypnrange[k] = pn[k] - bl;
  
       if (ypnrange[k] > val_max) {
         val_max = ypnrange[k];
         samplemax = k;
       }
     }
  
     chi2pn = pnfit->doFit(samplemax, &ypnrange[0]);
  
     if (chi2pn == 101 || chi2pn == 102 || chi2pn == 103)
       pnAmpl = 0.;
     else
       pnAmpl = pnfit->getAmpl();
  
     allPNAmpl.push_back(pnAmpl);
   }
  
   // ===========
   // Get Matacq
   // ===========
  
   ttrig = -1.;
   peak = -1.;
  
   if (doesRefFileExist == 1) {
     // FIXME
     if (color == 0)
       matacqTree->GetEntry(event - 1);
     else if (color == 3)
       matacqTree->GetEntry(matacqTree->GetEntries("color==0") + event - 1);
     ttrig = ttMat;
     peak = peakMat;
   }
  
   // ===========================
   // Decode EBDigis Information
   // ===========================
  
   double yrange[10];
   int adcGain = 0;
   int side = 0;
  
   if (EBDigi) {
     for (EBDigiCollection::const_iterator digiItr = EBDigi->begin(); digiItr != EBDigi->end();
          ++digiItr) {  // Loop on crystals
  
       EBDetId id_crystal(digiItr->id());
       EBDataFrame df(*digiItr);
  
       int etaG = id_crystal.ieta();  // global
       int phiG = id_crystal.iphi();  // global
  
       int etaL;  // local
       int phiL;  // local
       std::pair<int, int> LocalCoord = MEEBGeom::localCoord(etaG, phiG);
  
       etaL = LocalCoord.first;
       phiL = LocalCoord.second;
  
       eta = etaG;
       phi = phiG;
  
       side = MEEBGeom::side(etaG, phiG);
  
       EcalElectronicsId elecid_crystal = TheMapping->getElectronicsId(id_crystal);
  
       int towerID = elecid_crystal.towerId();
       // int channelID=elecid_crystal.channelId()-1;  // FIXME so far for endcap only
       int strip = elecid_crystal.stripId();
       int xtal = elecid_crystal.xtalId();
       int channelID = 5 * (strip - 1) + xtal - 1;  // FIXME
  
       int module = MEEBGeom::lmmod(etaG, phiG);
  
       std::pair<int, int> pnpair = MEEBGeom::pn(module);
       unsigned int MyPn0 = pnpair.first;
       unsigned int MyPn1 = pnpair.second;
  
       unsigned int channel = MEEBGeom::electronic_channel(etaL, phiL);
       assert(channel < nCrys);
  
       double adcmax = 0.0;
  
       if (towerID != int(_tower) || channelID != int(_channel) || dccID != int(_fedid - 600))
         continue;
       else
         channelNumber = channel;
  
       for (unsigned int i = 0; i < (*digiItr).size(); ++i) {  // Loop on adc samples
  
         EcalMGPASample samp_crystal(df.sample(i));
         adc[i] = samp_crystal.adc();
         adcG[i] = samp_crystal.gainId();
  
         if (i == 0)
           adcGain = adcG[i];
         if (i > 0)
           adcGain = TMath::Max(adcG[i], adcGain);
  
         if (adc[i] > adcmax) {
           adcmax = adc[i];
         }
       }
  
       for (dsum = 0., dsum1 = 0., k = 0; k < _presample; k++) {
         dsum += adc[k];
         if (k < _presample - 1)
           dsum1 += adc[k];
       }
  
       bl = dsum / ((double)_presample);
       bl1 = dsum1 / ((double)_presample - 1);
  
       for (val_max = 0., k = 0; k < _nsamples; k++) {
         yrange[k] = adc[k] - bl;
         if (yrange[k] > val_max) {
           val_max = yrange[k];
           samplemax = k;
         }
       }
  
       if (samplemax == 4 || samplemax == 5) {
         val_max = val_max + bl - bl1;
         for (k = 0; k < _nsamples; k++) {
           yrange[k] = yrange[k] + bl - bl1;
         }
       }
  
       for (unsigned int k = 0; k < _nsamples; k++) {
         adc[k] = yrange[k];
       }
  
       pn0 = allPNAmpl[MyPn0];
       pn1 = allPNAmpl[MyPn1];
  
       if (samplemax >= _timingcutlow && samplemax <= _timingcuthigh && lightside == side)
         ADCtrees->Fill();
     }
  
   } else {
     for (EEDigiCollection::const_iterator digiItr = EEDigi->begin(); digiItr != EEDigi->end();
          ++digiItr) {  // Loop on crystals
  
       EEDetId id_crystal(digiItr->id());
       EEDataFrame df(*digiItr);
  
       phi = id_crystal.ix() - 1;
       eta = id_crystal.iy() - 1;
       side = 0;  // FIXME
  
       // Recover the TT id and the electronic crystal numbering from EcalElectronicsMapping
  
       EcalElectronicsId elecid_crystal = TheMapping->getElectronicsId(id_crystal);
  
       int towerID = elecid_crystal.towerId();
       int channelID = elecid_crystal.channelId() - 1;
  
       int module = MEEEGeom::lmmod(phi, eta);
  
       std::pair<int, int> pnpair = MEEEGeom::pn(module, _fedid);
       unsigned int MyPn0 = pnpair.first;
       unsigned int MyPn1 = pnpair.second;
  
       unsigned int channel = MEEEGeom::crystal(phi, eta);
       assert(channel < nCrys);
  
       double adcmax = 0.0;
  
       if (towerID != int(_tower) || channelID != int(_channel) || dccID != int(_fedid - 600))
         continue;
       else
         channelNumber = channel;
  
       for (unsigned int i = 0; i < (*digiItr).size(); ++i) {  // Loop on adc samples
  
         EcalMGPASample samp_crystal(df.sample(i));
         adc[i] = samp_crystal.adc();
         adcG[i] = samp_crystal.gainId();
  
         if (i == 0)
           adcGain = adcG[i];
         if (i > 0)
           adcGain = TMath::Max(adcG[i], adcGain);
  
         if (adc[i] > adcmax) {
           adcmax = adc[i];
         }
       }
  
       for (dsum = 0., dsum1 = 0., k = 0; k < _presample; k++) {
         dsum += adc[k];
         if (k < _presample - 1)
           dsum1 += adc[k];
       }
  
       bl = dsum / ((double)_presample);
       bl1 = dsum1 / ((double)_presample - 1);
  
       for (val_max = 0., k = 0; k < _nsamples; k++) {
         yrange[k] = adc[k] - bl;
         if (yrange[k] > val_max) {
           val_max = yrange[k];
           samplemax = k;
         }
       }
  
       if (samplemax == 4 || samplemax == 5) {
         val_max = val_max + bl - bl1;
         for (k = 0; k < _nsamples; k++) {
           yrange[k] = yrange[k] + bl - bl1;
         }
       }
  
       for (unsigned int k = 0; k < _nsamples; k++) {
         adc[k] = yrange[k];
       }
  
       pn0 = allPNAmpl[MyPn0];
       pn1 = allPNAmpl[MyPn1];
  
       if (samplemax >= _timingcutlow && samplemax <= _timingcuthigh && lightside == side)
         ADCtrees->Fill();
     }
   }
  
 }  // analyze

References _channel, _ecalPart, _fedid, _firstsamplePN, _lastsamplePN, _nsamples, _nsamplesPN, _presample, _presamplePN, _timingcuthigh, _timingcutlow, _tower, adc, adcG, ADCtrees, cms::cuda::assert(), edm::SortedCollection< T, SORT >::begin(), c, EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, EcalElectronicsId::channelId(), channelNumber, colMat, color, gather_cfg::cout, MEEEGeom::crystal(), dccID, hgcalPerformanceValidation::df, digiCollection_, digiPNCollection_, digiProducer_, doesRefFileExist, TPNFit::doFit(), MillePedeFileConverter_cfg::e, CollectionTags_cfi::EBDigi, CollectionTags_cfi::EEDigi, MEEBGeom::electronic_channel(), edm::SortedCollection< T, SORT >::end(), eta, eventHeaderCollection_, eventHeaderProducer_, evtMat, cppFunctionSkipper::exception, fedID, spr::find(), get, TPNFit::getAmpl(), EcalElectronicsMapping::getElectronicsId(), mps_fire::i, iEvent, TPNFit::init(), IsFileCreated, dqmdumpme::k, EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_STD, lightside, MEEEGeom::lmmod(), MEEBGeom::lmmod(), MEEBGeom::localCoord(), matacqFile, matacqTree, Max(), nCrys, peak, peakMat, phi, MEEBGeom::pn(), MEEEGeom::pn(), pn, pn0, pn1, pnAmpl, edm::Handle< T >::product(), edm::ESHandle< T >::product(), resdir_, resfile, runNum, runType, MEEBGeom::side(), digitizers_cfi::strip, EcalElectronicsId::stripId(), ctpps_dqm_sourceclient-live_cfg::test, EcalElectronicsId::towerId(), ecaldqm::towerID(), ttMat, ttrig, EcalDCCHeaderBlock::EcalDCCEventSettings::wavelength, and EcalElectronicsId::xtalId().

◆ beginJob()

void EcalPerEvtLaserAnalyzer::beginJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 118 of file EcalPerEvtLaserAnalyzer.cc.

                                        {
   //========================================================================
  
   // Define temporary files' names
  
   stringstream namefile1;
   namefile1 << resdir_ << "/ADCSamples.root";
  
   ADCfile = namefile1.str();
  
   // Create temporary file and trees to save adc samples
  
   ADCFile = new TFile(ADCfile.c_str(), "RECREATE");
  
   stringstream name;
   name << "ADCTree";
  
   ADCtrees = new TTree(name.str().c_str(), name.str().c_str());
  
   ADCtrees->Branch("iphi", &phi, "phi/I");
   ADCtrees->Branch("ieta", &eta, "eta/I");
   ADCtrees->Branch("dccID", &dccID, "dccID/I");
   ADCtrees->Branch("event", &event, "event/I");
   ADCtrees->Branch("color", &color, "color/I");
   ADCtrees->Branch("adc", &adc, "adc[10]/D");
   ADCtrees->Branch("ttrigMatacq", &ttrig, "ttrig/D");
   ADCtrees->Branch("peakMatacq", &peak, "peak/D");
   ADCtrees->Branch("pn0", &pn0, "pn0/D");
   ADCtrees->Branch("pn1", &pn1, "pn1/D");
  
   ADCtrees->SetBranchAddress("ieta", &eta);
   ADCtrees->SetBranchAddress("iphi", &phi);
   ADCtrees->SetBranchAddress("dccID", &dccID);
   ADCtrees->SetBranchAddress("event", &event);
   ADCtrees->SetBranchAddress("color", &color);
   ADCtrees->SetBranchAddress("adc", adc);
   ADCtrees->SetBranchAddress("ttrigMatacq", &ttrig);
   ADCtrees->SetBranchAddress("peakMatacq", &peak);
   ADCtrees->SetBranchAddress("pn0", &pn0);
   ADCtrees->SetBranchAddress("pn1", &pn1);
  
   IsFileCreated = 0;
   nCh = nCrys;
 }

References adc, ADCfile, ADCFile, ADCtrees, color, dccID, eta, IsFileCreated, Skims_PA_cff::name, nCh, nCrys, peak, phi, pn0, pn1, resdir_, and ttrig.

◆ endJob()

void EcalPerEvtLaserAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 558 of file EcalPerEvtLaserAnalyzer.cc.

                                      {
   //========================================================================
  
   assert(colors.size() <= nColor);
   unsigned int nCol = colors.size();
  
   ADCtrees->Write();
  
   cout << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+" << endl;
   cout << "\t+=+       Analyzing laser data: getting per event               +=+" << endl;
   cout << "\t+=+            APD Amplitudes and ADC samples                   +=+" << endl;
   cout << "\t+=+    for fed:" << _fedid << ", tower:" << _tower << ", and channel:" << _channel << endl;
  
   // Define temporary tree to save APD amplitudes
  
   APDFile = new TFile(resfile.c_str(), "RECREATE");
  
   int ieta, iphi, channelID, towerID, flag;
   double alpha, beta;
  
   colors.push_back(color);
  
   for (unsigned int i = 0; i < nCol; i++) {
     stringstream name1;
     name1 << "headerCol" << colors[i];
  
     header[i] = new TTree(name1.str().c_str(), name1.str().c_str());
  
     header[i]->Branch("alpha", &alpha, "alpha/D");
     header[i]->Branch("beta", &beta, "beta/D");
     header[i]->Branch("iphi", &iphi, "iphi/I");
     header[i]->Branch("ieta", &ieta, "ieta/I");
     header[i]->Branch("dccID", &dccID, "dccID/I");
     header[i]->Branch("towerID", &towerID, "towerID/I");
     header[i]->Branch("channelID", &channelID, "channelID/I");
  
     header[i]->SetBranchAddress("alpha", &alpha);
     header[i]->SetBranchAddress("beta", &beta);
     header[i]->SetBranchAddress("iphi", &iphi);
     header[i]->SetBranchAddress("ieta", &ieta);
     header[i]->SetBranchAddress("dccID", &dccID);
     header[i]->SetBranchAddress("towerID", &towerID);
     header[i]->SetBranchAddress("channelID", &channelID);
   }
  
   stringstream name2;
   name2 << "APDTree";
   APDtrees = new TTree(name2.str().c_str(), name2.str().c_str());
  
   //List of branches
  
   APDtrees->Branch("event", &event, "event/I");
   APDtrees->Branch("color", &color, "color/I");
   APDtrees->Branch("adc", &adc, "adc[10]/D");
   APDtrees->Branch("peakMatacq", &peak, "peak/D");
   APDtrees->Branch("ttrigMatacq", &ttrig, "ttrig/D");
   APDtrees->Branch("apdAmpl", &apdAmpl, "apdAmpl/D");
   APDtrees->Branch("apdTime", &apdTime, "apdTime/D");
   APDtrees->Branch("flag", &flag, "flag/I");
   APDtrees->Branch("pn0", &pn0, "pn0/D");
   APDtrees->Branch("pn1", &pn1, "pn1/D");
  
   APDtrees->SetBranchAddress("event", &event);
   APDtrees->SetBranchAddress("color", &color);
   APDtrees->SetBranchAddress("adc", adc);
   APDtrees->SetBranchAddress("peakMatacq", &peak);
   APDtrees->SetBranchAddress("ttrigMatacq", &ttrig);
   APDtrees->SetBranchAddress("apdAmpl", &apdAmpl);
   APDtrees->SetBranchAddress("apdTime", &apdTime);
   APDtrees->SetBranchAddress("flag", &flag);
   APDtrees->SetBranchAddress("pn0", &pn0);
   APDtrees->SetBranchAddress("pn1", &pn1);
  
   // Retrieve alpha and beta for APD amplitudes calculation
  
   TFile* alphaFile = new TFile(refalphabeta_.c_str());
   TTree* alphaTree[2];
  
   Double_t alphaRun, betaRun;
   int ietaRun, iphiRun, channelIDRun, towerIDRun, dccIDRun, flagRun;
  
   for (unsigned int i = 0; i < nCol; i++) {
     stringstream name3;
     name3 << "ABCol" << i;
     alphaTree[i] = (TTree*)alphaFile->Get(name3.str().c_str());
     alphaTree[i]->SetBranchStatus("*", false);
     alphaTree[i]->SetBranchStatus("alpha", true);
     alphaTree[i]->SetBranchStatus("beta", true);
     alphaTree[i]->SetBranchStatus("iphi", true);
     alphaTree[i]->SetBranchStatus("ieta", true);
     alphaTree[i]->SetBranchStatus("dccID", true);
     alphaTree[i]->SetBranchStatus("towerID", true);
     alphaTree[i]->SetBranchStatus("channelID", true);
     alphaTree[i]->SetBranchStatus("flag", true);
  
     alphaTree[i]->SetBranchAddress("alpha", &alphaRun);
     alphaTree[i]->SetBranchAddress("beta", &betaRun);
     alphaTree[i]->SetBranchAddress("iphi", &iphiRun);
     alphaTree[i]->SetBranchAddress("ieta", &ietaRun);
     alphaTree[i]->SetBranchAddress("dccID", &dccIDRun);
     alphaTree[i]->SetBranchAddress("towerID", &towerIDRun);
     alphaTree[i]->SetBranchAddress("channelID", &channelIDRun);
     alphaTree[i]->SetBranchAddress("flag", &flagRun);
   }
  
   PulseFitWithFunction* pslsfit = new PulseFitWithFunction();
  
   double chi2;
  
   for (unsigned int icol = 0; icol < nCol; icol++) {
     IsThereDataADC[icol] = 1;
     stringstream cut;
     cut << "color==" << colors.at(icol);
     if (ADCtrees->GetEntries(cut.str().c_str()) < 10)
       IsThereDataADC[icol] = 0;
     IsHeaderFilled[icol] = 0;
   }
  
   // Define submodule and channel number inside the submodule (as Patrice)
  
   Long64_t nbytes = 0, nb = 0;
   for (Long64_t jentry = 0; jentry < ADCtrees->GetEntriesFast(); jentry++) {  // Loop on events
     nb = ADCtrees->GetEntry(jentry);
     nbytes += nb;
  
     int iCry = channelNumber;
  
     // get back color
  
     unsigned int iCol = 0;
     for (unsigned int i = 0; i < nCol; i++) {
       if (color == colors[i]) {
         iCol = i;
         i = colors.size();
       }
     }
  
     alphaTree[iCol]->GetEntry(iCry);
  
     flag = flagRun;
     iphi = iphiRun;
     ieta = ietaRun;
     towerID = towerIDRun;
     channelID = channelIDRun;
     alpha = alphaRun;
     beta = betaRun;
  
     if (IsHeaderFilled[iCol] == 0) {
       header[iCol]->Fill();
       IsHeaderFilled[iCol] = 1;
     }
     // Amplitude calculation
  
     apdAmpl = 0;
     apdTime = 0;
  
     pslsfit->init(_nsamples, _firstsample, _lastsample, _niter, alphaRun, betaRun);
     chi2 = pslsfit->doFit(&adc[0]);
  
     if (chi2 < 0. || chi2 == 102 || chi2 == 101) {
       apdAmpl = 0;
       apdTime = 0;
  
     } else {
       apdAmpl = pslsfit->getAmpl();
       apdTime = pslsfit->getTime();
     }
  
     APDtrees->Fill();
   }
  
   alphaFile->Close();
  
   ADCFile->Close();
  
   APDFile->Write();
   APDFile->Close();
  
   // Remove unwanted files
  
   stringstream del;
   del << "rm " << ADCfile;
   system(del.str().c_str());
  
   cout << "\t+=+    .................................................. done  +=+" << endl;
   cout << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+" << endl;
 }

References _channel, _fedid, _firstsample, _lastsample, _niter, _nsamples, _tower, adc, ADCfile, ADCFile, ADCtrees, alpha, apdAmpl, APDFile, apdTime, APDtrees, cms::cuda::assert(), HLT_FULL_cff::beta, channelNumber, hltPixelTracks_cff::chi2, color, gather_cfg::cout, TkAlMuonSelectors_cfi::cut, dccID, PulseFitWithFunction::doFit(), RemoveAddSevLevel::flag, PulseFitWithFunction::getAmpl(), PulseFitWithFunction::getTime(), header, mps_fire::i, LEDCalibrationChannels::ieta, PulseFitWithFunction::init(), LEDCalibrationChannels::iphi, IsHeaderFilled, IsThereDataADC, hfnoseParametersInitialization_cfi::name2, nColor, peak, pn0, pn1, refalphabeta_, resfile, ecaldqm::towerID(), and ttrig.

Referenced by o2olib.O2ORunMgr::executeJob().

Member Data Documentation

◆ _channel

unsigned int EcalPerEvtLaserAnalyzer::_channel

private

Definition at line 64 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

◆ _ecalPart

std::string EcalPerEvtLaserAnalyzer::_ecalPart

private

Definition at line 65 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

◆ _fedid

int EcalPerEvtLaserAnalyzer::_fedid

private

Definition at line 62 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

◆ _firstsample

unsigned int EcalPerEvtLaserAnalyzer::_firstsample

private

Definition at line 51 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

◆ _firstsamplePN

unsigned int EcalPerEvtLaserAnalyzer::_firstsamplePN

private

Definition at line 57 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ _lastsample

unsigned int EcalPerEvtLaserAnalyzer::_lastsample

private

Definition at line 52 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

◆ _lastsamplePN

unsigned int EcalPerEvtLaserAnalyzer::_lastsamplePN

private

Definition at line 58 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ _niter

unsigned int EcalPerEvtLaserAnalyzer::_niter

private

Definition at line 61 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

◆ _nsamples

unsigned int EcalPerEvtLaserAnalyzer::_nsamples

private

Definition at line 49 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

◆ _nsamplesPN

unsigned int EcalPerEvtLaserAnalyzer::_nsamplesPN

private

Definition at line 55 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ _presample

unsigned int EcalPerEvtLaserAnalyzer::_presample

private

Definition at line 50 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ _presamplePN

unsigned int EcalPerEvtLaserAnalyzer::_presamplePN

private

Definition at line 56 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ _samplemax

unsigned int EcalPerEvtLaserAnalyzer::_samplemax

private

Definition at line 54 of file EcalPerEvtLaserAnalyzer.h.

◆ _samplemin

unsigned int EcalPerEvtLaserAnalyzer::_samplemin

private

Definition at line 53 of file EcalPerEvtLaserAnalyzer.h.

◆ _timingcuthigh

unsigned int EcalPerEvtLaserAnalyzer::_timingcuthigh

private

Definition at line 60 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ _timingcutlow

unsigned int EcalPerEvtLaserAnalyzer::_timingcutlow

private

Definition at line 59 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ _tower

unsigned int EcalPerEvtLaserAnalyzer::_tower

private

Definition at line 63 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

◆ adc

double EcalPerEvtLaserAnalyzer::adc[10]

private

Definition at line 130 of file EcalPerEvtLaserAnalyzer.h.

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

◆ ADCfile

std::string EcalPerEvtLaserAnalyzer::ADCfile

private

Definition at line 78 of file EcalPerEvtLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

◆ ADCFile

TFile* EcalPerEvtLaserAnalyzer::ADCFile

private

Definition at line 111 of file EcalPerEvtLaserAnalyzer.h.

Referenced by beginJob(), and endJob().

◆ adcG

int EcalPerEvtLaserAnalyzer::adcG[10]

private

Definition at line 131 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ ADCtrees

TTree* EcalPerEvtLaserAnalyzer::ADCtrees

private

Definition at line 112 of file EcalPerEvtLaserAnalyzer.h.

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

◆ alphafile

std::string EcalPerEvtLaserAnalyzer::alphafile

private

Definition at line 77 of file EcalPerEvtLaserAnalyzer.h.

◆ apdAmpl

double EcalPerEvtLaserAnalyzer::apdAmpl

private

Definition at line 136 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

◆ APDfile

std::string EcalPerEvtLaserAnalyzer::APDfile

private

Definition at line 79 of file EcalPerEvtLaserAnalyzer.h.

◆ APDFile

TFile* EcalPerEvtLaserAnalyzer::APDFile

private

Definition at line 114 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

◆ apdTime

double EcalPerEvtLaserAnalyzer::apdTime

private

Definition at line 137 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

◆ APDtrees

TTree* EcalPerEvtLaserAnalyzer::APDtrees

private

Definition at line 116 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

◆ channelNumber

int EcalPerEvtLaserAnalyzer::channelNumber

private

Definition at line 103 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

◆ colMat

int EcalPerEvtLaserAnalyzer::colMat

private

Definition at line 123 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ color

int EcalPerEvtLaserAnalyzer::color

private

Definition at line 129 of file EcalPerEvtLaserAnalyzer.h.

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

◆ colors

std::vector<int> EcalPerEvtLaserAnalyzer::colors

private

Definition at line 105 of file EcalPerEvtLaserAnalyzer.h.

◆ dccID

int EcalPerEvtLaserAnalyzer::dccID

private

Definition at line 99 of file EcalPerEvtLaserAnalyzer.h.

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

◆ digiCollection_

std::string EcalPerEvtLaserAnalyzer::digiCollection_

private

Definition at line 69 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

◆ digiPNCollection_

std::string EcalPerEvtLaserAnalyzer::digiPNCollection_

private

Definition at line 70 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

◆ digiProducer_

std::string EcalPerEvtLaserAnalyzer::digiProducer_

private

Definition at line 71 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

◆ doesRefFileExist

int EcalPerEvtLaserAnalyzer::doesRefFileExist

private

Definition at line 118 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ eta

int EcalPerEvtLaserAnalyzer::eta

private

Definition at line 127 of file EcalPerEvtLaserAnalyzer.h.

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

◆ event

int EcalPerEvtLaserAnalyzer::event

private

Definition at line 128 of file EcalPerEvtLaserAnalyzer.h.

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

◆ eventHeaderCollection_

std::string EcalPerEvtLaserAnalyzer::eventHeaderCollection_

private

Definition at line 72 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

◆ eventHeaderProducer_

std::string EcalPerEvtLaserAnalyzer::eventHeaderProducer_

private

Definition at line 73 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and EcalPerEvtLaserAnalyzer().

◆ evtMat

int EcalPerEvtLaserAnalyzer::evtMat

private

Definition at line 123 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ fedID

int EcalPerEvtLaserAnalyzer::fedID

private

Definition at line 100 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ header

TTree* EcalPerEvtLaserAnalyzer::header[2]

private

Definition at line 115 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

◆ iEvent

int EcalPerEvtLaserAnalyzer::iEvent

private

Definition at line 45 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ IsFileCreated

int EcalPerEvtLaserAnalyzer::IsFileCreated

private

Definition at line 89 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ IsHeaderFilled

int EcalPerEvtLaserAnalyzer::IsHeaderFilled[nColor]

private

Definition at line 120 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

◆ IsThereDataADC

int EcalPerEvtLaserAnalyzer::IsThereDataADC[nColor]

private

Definition at line 119 of file EcalPerEvtLaserAnalyzer.h.

Referenced by endJob().

◆ lightside

int EcalPerEvtLaserAnalyzer::lightside

private

Definition at line 101 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ matacqFile

TFile* EcalPerEvtLaserAnalyzer::matacqFile

private

Definition at line 108 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ matacqTree

TTree* EcalPerEvtLaserAnalyzer::matacqTree

private

Definition at line 109 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ nCh

unsigned int EcalPerEvtLaserAnalyzer::nCh

private

Definition at line 93 of file EcalPerEvtLaserAnalyzer.h.

Referenced by beginJob().

◆ nCrys

unsigned int EcalPerEvtLaserAnalyzer::nCrys

private

Definition at line 85 of file EcalPerEvtLaserAnalyzer.h.

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

◆ nSides

unsigned int EcalPerEvtLaserAnalyzer::nSides

private

Definition at line 87 of file EcalPerEvtLaserAnalyzer.h.

◆ nTT

unsigned int EcalPerEvtLaserAnalyzer::nTT

private

Definition at line 86 of file EcalPerEvtLaserAnalyzer.h.

Referenced by EcalPerEvtLaserAnalyzer().

◆ peak

double EcalPerEvtLaserAnalyzer::peak

private

Definition at line 122 of file EcalPerEvtLaserAnalyzer.h.

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

◆ peakMat

double EcalPerEvtLaserAnalyzer::peakMat

private

Definition at line 122 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ phi

int EcalPerEvtLaserAnalyzer::phi

private

Definition at line 127 of file EcalPerEvtLaserAnalyzer.h.

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

◆ pn

double EcalPerEvtLaserAnalyzer::pn[50]

private

Definition at line 135 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ pn0

double EcalPerEvtLaserAnalyzer::pn0

private

Definition at line 134 of file EcalPerEvtLaserAnalyzer.h.

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

◆ pn1

double EcalPerEvtLaserAnalyzer::pn1

private

Definition at line 134 of file EcalPerEvtLaserAnalyzer.h.

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

◆ pnAmpl

double EcalPerEvtLaserAnalyzer::pnAmpl

private

Definition at line 138 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ refalphabeta_

std::string EcalPerEvtLaserAnalyzer::refalphabeta_

private

Definition at line 68 of file EcalPerEvtLaserAnalyzer.h.

Referenced by EcalPerEvtLaserAnalyzer(), and endJob().

◆ resdir_

std::string EcalPerEvtLaserAnalyzer::resdir_

private

Definition at line 67 of file EcalPerEvtLaserAnalyzer.h.

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

◆ resfile

std::string EcalPerEvtLaserAnalyzer::resfile

private

Definition at line 80 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze(), and endJob().

◆ runNum

int EcalPerEvtLaserAnalyzer::runNum

private

Definition at line 98 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ runType

int EcalPerEvtLaserAnalyzer::runType

private

Definition at line 97 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ tt

double EcalPerEvtLaserAnalyzer::tt

private

Definition at line 132 of file EcalPerEvtLaserAnalyzer.h.

◆ ttMat

double EcalPerEvtLaserAnalyzer::ttMat

private

Definition at line 122 of file EcalPerEvtLaserAnalyzer.h.

Referenced by analyze().

◆ ttrig

double EcalPerEvtLaserAnalyzer::ttrig

private

Definition at line 133 of file EcalPerEvtLaserAnalyzer.h.

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

Public Types

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ VarCol

Constructor & Destructor Documentation

◆ EcalPerEvtLaserAnalyzer()

◆ ~EcalPerEvtLaserAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

Member Data Documentation

◆ _channel

◆ _ecalPart

◆ _fedid

◆ _firstsample

◆ _firstsamplePN

◆ _lastsample

◆ _lastsamplePN

◆ _niter

◆ _nsamples

◆ _nsamplesPN

◆ _presample

◆ _presamplePN

◆ _samplemax

◆ _samplemin

◆ _timingcuthigh

◆ _timingcutlow

◆ _tower

◆ adc

◆ ADCfile

◆ ADCFile

◆ adcG

◆ ADCtrees

◆ alphafile

◆ apdAmpl

◆ APDfile

◆ APDFile

◆ apdTime

◆ APDtrees

◆ channelNumber

◆ colMat

◆ color

◆ colors

◆ dccID

◆ digiCollection_

◆ digiPNCollection_

◆ digiProducer_

◆ doesRefFileExist

◆ eta

◆ event

◆ eventHeaderCollection_

◆ eventHeaderProducer_

◆ evtMat

◆ fedID

◆ header

◆ iEvent

◆ IsFileCreated

◆ IsHeaderFilled

◆ IsThereDataADC

◆ lightside

◆ matacqFile

◆ matacqTree

◆ nCh

◆ nCrys

◆ nSides

◆ nTT

◆ peak

◆ peakMat

◆ phi

◆ pn

◆ pn0

◆ pn1

◆ pnAmpl

◆ refalphabeta_

◆ resdir_

◆ resfile