EcalMatacqAnalyzer Class Reference

#include <EcalMatacqAnalyzer.h>

Inheritance diagram for EcalMatacqAnalyzer:

Public Types
enum	VarCol { iBlue, iRed, nColor }
enum	VarSide { iSide0, iSide1, nSide }
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	beginJob () override
	EcalMatacqAnalyzer (const edm::ParameterSet &iConfig)
void	endJob () override
	~EcalMatacqAnalyzer () override=default
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const noexcept
bool	wantsStreamRuns () const noexcept
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
const double	_cutwindow
const int	_debug
const int	_fedid
const unsigned int	_highlev
const unsigned int	_lowlev
const unsigned int	_nevlasers
const unsigned int	_noiseCut
const unsigned int	_nsamplesaftmax
const unsigned int	_nsamplesshape
const unsigned int	_parabnaftmax
const unsigned int	_parabnbefmax
const double	_presample
const unsigned int	_presampleshape
const unsigned int	_slide
const unsigned int	_thres
const unsigned int	_timeaftmax
const unsigned int	_timebefmax
double	ampl
int	color
std::vector< int >	colors
int	dccID
const edm::EDGetTokenT< EcalRawDataCollection >	dccToken_
const std::string	digiCollection_
const std::string	digiProducer_
int	event
const std::string	eventHeaderCollection_
const std::string	eventHeaderProducer_
int	fedID
double	fit
double	fw20
double	fw80
double	fwhm
int	iEvent
bool	isThereMatacq
int	laserEvents
int	lightside
double	matacq [2560]
int	maxsamp
TTree *	meanTree [nColor]
TMTQ *	MTQ [nColor][nSide]
TTree *	mtqShape
int	nsamples
unsigned int	nSides
std::string	outfile
TFile *	outFile
double	peak
double	ped
double	pedsig
const edm::EDGetTokenT< EcalMatacqDigiCollection >	pmatToken_
const std::string	resdir_
int	runNum
int	runType
std::string	sampfile
TFile *	sampFile
double	sigma
double	sliding
int	status
TTree *	tree
double	trise
double	tt
double	ttrig

Additional Inherited Members
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 16 of file EcalMatacqAnalyzer.h.

Member Enumeration Documentation

VarCol

enum EcalMatacqAnalyzer::VarCol

Enumerator
iBlue
iRed
nColor

Definition at line 25 of file EcalMatacqAnalyzer.h.

{ iBlue, iRed, nColor };

VarSide

enum EcalMatacqAnalyzer::VarSide

Enumerator
iSide0
iSide1
nSide

Definition at line 26 of file EcalMatacqAnalyzer.h.

{ iSide0, iSide1, nSide };

Constructor & Destructor Documentation

EcalMatacqAnalyzer()

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

explicit

Definition at line 34 of file EcalMatacqAnalyzer.cc.

    :

      iEvent(0),

      // framework parameters with default values

      _presample(iConfig.getUntrackedParameter<double>("nPresamples", 6.7)),
      _nsamplesaftmax(iConfig.getUntrackedParameter<unsigned int>("nSamplesAftMax", 80)),
      _noiseCut(iConfig.getUntrackedParameter<unsigned int>("noiseCut", 7)),
      _parabnbefmax(iConfig.getUntrackedParameter<unsigned int>("paraBeforeMax", 8)),
      _parabnaftmax(iConfig.getUntrackedParameter<unsigned int>("paraAfterMax", 7)),
      _thres(iConfig.getUntrackedParameter<unsigned int>("threshold", 10)),
      _lowlev(iConfig.getUntrackedParameter<unsigned int>("lowLevel", 20)),
      _highlev(iConfig.getUntrackedParameter<unsigned int>("highLevel", 80)),
      _nevlasers(iConfig.getUntrackedParameter<unsigned int>("nEventLaser", 600)),
      _timebefmax(iConfig.getUntrackedParameter<unsigned int>("timeBefMax", 100)),
      _timeaftmax(iConfig.getUntrackedParameter<unsigned int>("timeAftMax", 250)),
      _cutwindow(iConfig.getUntrackedParameter<double>("cutWindow", 0.1)),
      _nsamplesshape(iConfig.getUntrackedParameter<unsigned int>("nSamplesShape", 250)),
      _presampleshape(iConfig.getUntrackedParameter<unsigned int>("nPresamplesShape", 50)),
      _slide(iConfig.getUntrackedParameter<unsigned int>("nSlide", 100)),
      _fedid(iConfig.getUntrackedParameter<int>("fedID", -999)),
      _debug(iConfig.getUntrackedParameter<int>("debug", 0)),
      resdir_(iConfig.getUntrackedParameter<std::string>("resDir")),
      digiCollection_(iConfig.getParameter<std::string>("digiCollection")),
      digiProducer_(iConfig.getParameter<std::string>("digiProducer")),
      eventHeaderCollection_(iConfig.getParameter<std::string>("eventHeaderCollection")),
      eventHeaderProducer_(iConfig.getParameter<std::string>("eventHeaderProducer")),
      pmatToken_(consumes<EcalMatacqDigiCollection>(edm::InputTag(digiProducer_, digiCollection_))),
      dccToken_(consumes<EcalRawDataCollection>(edm::InputTag(eventHeaderProducer_, eventHeaderCollection_))),
      nSides(NSIDES),
      lightside(0),
      runType(-1),
      runNum(0),
      event(0),
      color(-1),
      maxsamp(0),
      nsamples(0),
      tt(0)

//========================================================================
{
  //now do what ever initialization is needed
}

~EcalMatacqAnalyzer()

EcalMatacqAnalyzer::~EcalMatacqAnalyzer ( )

overridedefault

Member Function Documentation

analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 126 of file EcalMatacqAnalyzer.cc.

References _debug, _fedid, EcalMatacqDigi::adcCount(), edm::SortedCollection< T, SORT >::begin(), color, dccID, dccToken_, digiCollection_, digiProducer_, MillePedeFileConverter_cfg::e, edm::SortedCollection< T, SORT >::end(), eventHeaderCollection_, eventHeaderProducer_, fedID, spr::find(), mps_fire::i, iEvent, ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::if(), createfilelist::int, isThereMatacq, edm::HandleBase::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, EcalDCCHeaderBlock::LASER_DELAY_SCAN, EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_POWER_SCAN, EcalDCCHeaderBlock::LASER_STD, laserEvents, lightside, matacq, SiStripPI::max, maxsamp, nsamples, pmatToken_, edm::Handle< T >::product(), runNum, runType, EcalMatacqDigi::size(), EcalMatacqDigi::tTrig(), and EcalDCCHeaderBlock::EcalDCCEventSettings::wavelength.

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

  ++iEvent;

  if (_debug == 1)
    edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- Entering Analyze -- event= " << iEvent;

  // retrieving MATACQ :
  const edm::Handle<EcalMatacqDigiCollection>& pmatacqDigi = e.getHandle(pmatToken_);
  const EcalMatacqDigiCollection* matacqDigi = (pmatacqDigi.isValid()) ? pmatacqDigi.product() : nullptr;
  if (pmatacqDigi.isValid()) {
    if (_debug == 1)
      edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- Matacq Digis Found -- ";
  } else {
    edm::LogError("EcalMatacqAnalyzzer") << "Error! can't get the product EcalMatacqDigi producer:"
                                         << digiProducer_.c_str() << " collection:" << digiCollection_.c_str();
    return;
  }

  // retrieving DCC header

  const edm::Handle<EcalRawDataCollection>& pDCCHeader = e.getHandle(dccToken_);
  const EcalRawDataCollection* DCCHeader = (pDCCHeader.isValid()) ? pDCCHeader.product() : nullptr;
  if (!pDCCHeader.isValid()) {
    edm::LogError("EcalMatacqAnalyzzer") << "Error! can't get the product EcalRawData producer:"
                                         << eventHeaderProducer_.c_str()
                                         << " collection:" << eventHeaderCollection_.c_str();
    return;
  }

  // ====================================
  // Decode Basic DCCHeader Information
  // ====================================

  if (_debug == 1)
    edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- Before header -- ";

  for (EcalRawDataCollection::const_iterator headerItr = DCCHeader->begin(); headerItr != DCCHeader->end();
       ++headerItr) {
    EcalDCCHeaderBlock::EcalDCCEventSettings settings = headerItr->getEventSettings();
    color = (int)settings.wavelength;
    if (color < 0)
      return;

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

    if (_debug == 1)
      edm::LogVerbatim("EcalMatacqAnalyzzer")
          << "-- debug test -- runtype:" << runType << " event:" << event << " runNum:" << runNum;

    dccID = headerItr->getDccInTCCCommand();
    fedID = headerItr->fedId();
    lightside = headerItr->getRtHalf();

    //assert (lightside<2 && lightside>=0);

    if (lightside != 1 && lightside != 0) {
      edm::LogVerbatim("EcalMatacqAnalyzzer") << "Unexpected lightside: " << lightside << " for event " << iEvent;
      return;
    }
    if (_debug == 1) {
      edm::LogVerbatim("EcalMatacqAnalyzzer")
          << "-- debug test -- Inside header before fed cut -- color=" << color << ", dcc=" << dccID
          << ", fed=" << fedID << ",  lightside=" << lightside << ", runType=" << runType;
    }

    // take event only if the fed corresponds to the DCC in TCC
    if (600 + dccID != fedID)
      continue;

    if (_debug == 1) {
      edm::LogVerbatim("EcalMatacqAnalyzzer")
          << "-- debug test -- Inside header after fed cut -- color=" << color << ", dcc=" << dccID << ", fed=" << fedID
          << ",  lightside=" << lightside << ", runType=" << runType;
    }

    // Cut on runType

    if (runType != EcalDCCHeaderBlock::LASER_STD && runType != EcalDCCHeaderBlock::LASER_GAP &&
        runType != EcalDCCHeaderBlock::LASER_POWER_SCAN && runType != EcalDCCHeaderBlock::LASER_DELAY_SCAN)
      return;

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

  if (_debug == 1)
    edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- Before digis -- Event:" << iEvent;

  // Count laser events
  laserEvents++;

  // ===========================
  // Decode Matacq Information
  // ===========================

  double max = 0;

  for (EcalMatacqDigiCollection::const_iterator it = matacqDigi->begin(); it != matacqDigi->end();
       ++it) {  // Loop on matacq channel

    //
    const EcalMatacqDigi& digis = *it;

    if (_debug == 1) {
      edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- Inside digis -- digi size=" << digis.size();
    }

    if (digis.size() == 0)
      continue;
    else
      isThereMatacq = true;

    max = 0;
    maxsamp = 0;
    nsamples = digis.size();
    tt = digis.tTrig();

    for (int i = 0; i < digis.size(); ++i) {  // Loop on matacq samples
      matacq[i] = -digis.adcCount(i);
      if (matacq[i] > max) {
        max = matacq[i];
        maxsamp = i;
      }
    }
    if (_debug == 1) {
      edm::LogVerbatim("EcalMatacqAnalyzzer")
          << "-- debug test -- Inside digis -- nsamples=" << nsamples << ", max=" << max;
    }
  }

  if (_debug == 1)
    edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- After digis -- Event: " << iEvent;
  tree->Fill();

}  // analyze

beginJob()

void EcalMatacqAnalyzer::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 82 of file EcalMatacqAnalyzer.cc.

References color, isThereMatacq, laserEvents, lightside, matacq, maxsamp, nsamples, outfile, resdir_, sampfile, and sampFile.

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

  // Define temporary file name

  sampfile = resdir_;
  sampfile += "/TmpTreeMatacqAnalyzer.root";

  sampFile = new TFile(sampfile.c_str(), "RECREATE");

  // declaration of the tree to fill

  tree = new TTree("MatacqTree", "MatacqTree");

  //List of branches

  tree->Branch("event", &event, "event/I");
  tree->Branch("color", &color, "color/I");
  tree->Branch("matacq", &matacq, "matacq[2560]/D");
  tree->Branch("nsamples", &nsamples, "nsamples/I");
  tree->Branch("maxsamp", &maxsamp, "maxsamp/I");
  tree->Branch("tt", &tt, "tt/D");
  tree->Branch("lightside", &lightside, "lightside/I");

  tree->SetBranchAddress("event", &event);
  tree->SetBranchAddress("color", &color);
  tree->SetBranchAddress("matacq", matacq);
  tree->SetBranchAddress("nsamples", &nsamples);
  tree->SetBranchAddress("maxsamp", &maxsamp);
  tree->SetBranchAddress("tt", &tt);
  tree->SetBranchAddress("lightside", &lightside);

  // Define output results files' names

  std::stringstream namefile;
  namefile << resdir_ << "/MATACQ.root";
  outfile = namefile.str();

  // Laser events counter
  laserEvents = 0;
  isThereMatacq = false;
}

endJob()

void EcalMatacqAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 277 of file EcalMatacqAnalyzer.cc.

References _cutwindow, _debug, _highlev, _lowlev, _nevlasers, _noiseCut, _nsamplesaftmax, _nsamplesshape, _parabnaftmax, _parabnbefmax, _presample, _presampleshape, _slide, _thres, _timeaftmax, _timebefmax, a0, testProducerWithPsetDescEmpty_cfi::a2, ampl, cms::cuda::assert(), color, TMatacq::compute_trise(), TMatacq::doFit(), TMatacq::findPeak(), fw20, fw80, fwhm, TMTQ::get(), TMatacq::getAmpl(), TMatacq::getBaseLine(), TMatacq::getFwhm(), TMatacq::getsigBaseLine(), TMatacq::getsigTimpeak(), TMatacq::getSlide(), TMatacq::getTimax(), TMatacq::getTimpeak(), TMatacq::getTrise(), TMatacq::getWidth20(), TMatacq::getWidth80(), mps_fire::i, TMTQ::iAmpl, TMTQ::iFit, TMTQ::iFw20, TMTQ::iFw80, TMTQ::iFwhm, createfilelist::int, TMTQ::iPeak, TMTQ::iPed, TMTQ::iPedsig, TMTQ::iSigma, TMTQ::iSlide, isThereMatacq, TMTQ::iTrise, lightside, matacq, maxsamp, meanTree, MTQ, mtqShape, nColor, TMTQ::nOutVar, nsamples, nSide, nSides, outfile, outFile, peak, ped, pedsig, TMatacq::rawPulseAnalysis(), sampfile, sampFile, sigma, sliding, status, ctpps_dqm_sourceclient-live_cfg::test, trise, tt, ttrig, and heppy_batch::val.

Referenced by o2olib.O2ORunMgr::executeJob().

                                {
  // Don't do anything if there is no events
  if (!isThereMatacq) {
    edm::LogVerbatim("EcalMatacqAnalyzzer") << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
    edm::LogVerbatim("EcalMatacqAnalyzzer") << "\t+=+     WARNING! NO MATACQ        +=+";
    edm::LogVerbatim("EcalMatacqAnalyzzer") << "\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";

    // Remove temporary file
    FILE* test;
    test = fopen(sampfile.c_str(), "r");
    if (test) {
      std::stringstream del2;
      del2 << "rm " << sampfile;
      system(del2.str().c_str());
    }
    return;
  }

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

  for (unsigned int iCol = 0; iCol < nCol; iCol++) {
    for (unsigned int iSide = 0; iSide < nSide; iSide++) {
      MTQ[iCol][iSide] = new TMTQ();
    }
  }

  outFile = new TFile(outfile.c_str(), "RECREATE");

  TProfile* shapeMat = new TProfile("shapeLaser", "shapeLaser", _nsamplesshape, -0.5, double(_nsamplesshape) - 0.5);
  TProfile* shapeMatTmp = new TProfile(
      "shapeLaserTmp", "shapeLaserTmp", _timeaftmax + _timebefmax, -0.5, double(_timeaftmax + _timebefmax) - 0.5);

  edm::LogVerbatim("EcalMatacqAnalyzzer") << "\n\t+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+";
  edm::LogVerbatim("EcalMatacqAnalyzzer") << "\t+=+     Analyzing MATACQ data     +=+";

  //
  // create output ntuple
  //

  mtqShape = new TTree("MatacqShape", "MatacqShape");

  // list of branches
  // keep Patrice's notations

  mtqShape->Branch("event", &event, "event/I");
  mtqShape->Branch("color", &color, "color/I");
  mtqShape->Branch("status", &status, "status/I");
  mtqShape->Branch("peak", &peak, "peak/D");
  mtqShape->Branch("sigma", &sigma, "sigma/D");
  mtqShape->Branch("fit", &fit, "fit/D");
  mtqShape->Branch("ampl", &ampl, "ampl/D");
  mtqShape->Branch("trise", &trise, "trise/D");
  mtqShape->Branch("fwhm", &fwhm, "fwhm/D");
  mtqShape->Branch("fw20", &fw20, "fw20/D");
  mtqShape->Branch("fw80", &fw80, "fw80/D");
  mtqShape->Branch("ped", &ped, "ped/D");
  mtqShape->Branch("pedsig", &pedsig, "pedsig/D");
  mtqShape->Branch("ttrig", &ttrig, "ttrig/D");
  mtqShape->Branch("sliding", &sliding, "sliding/D");

  mtqShape->SetBranchAddress("event", &event);
  mtqShape->SetBranchAddress("color", &color);
  mtqShape->SetBranchAddress("status", &status);
  mtqShape->SetBranchAddress("peak", &peak);
  mtqShape->SetBranchAddress("sigma", &sigma);
  mtqShape->SetBranchAddress("fit", &fit);
  mtqShape->SetBranchAddress("ampl", &ampl);
  mtqShape->SetBranchAddress("fwhm", &fwhm);
  mtqShape->SetBranchAddress("fw20", &fw20);
  mtqShape->SetBranchAddress("fw80", &fw80);
  mtqShape->SetBranchAddress("trise", &trise);
  mtqShape->SetBranchAddress("ped", &ped);
  mtqShape->SetBranchAddress("pedsig", &pedsig);
  mtqShape->SetBranchAddress("ttrig", &ttrig);
  mtqShape->SetBranchAddress("sliding", &sliding);

  unsigned int endsample;
  unsigned int presample;

  // loop over the entries of the tree
  TChain* fChain = (TChain*)tree;
  Long64_t nentries = fChain->GetEntriesFast();

  for (Long64_t jentry = 0; jentry < nentries; jentry++) {
    // load the event
    Long64_t ientry = fChain->LoadTree(jentry);
    if (ientry < 0)
      break;
    fChain->GetEntry(jentry);

    status = 0;
    peak = -1;
    sigma = 0;
    fit = -1;
    ampl = -1;
    trise = -1;
    ttrig = tt;
    fwhm = 0;
    fw20 = 0;
    fw80 = 0;
    ped = 0;
    pedsig = 0;
    sliding = 0;

    if (_debug == 1)
      edm::LogVerbatim("EcalMatacqAnalyzzer")
          << "-- debug test -- inside loop 1  -- jentry:" << jentry << " over nentries=" << nentries;

    // create the object for Matacq data analysis

    endsample = maxsamp + _nsamplesaftmax;
    presample = int(_presample * nsamples / 100.);
    TMatacq* mtq = new TMatacq(nsamples,
                               presample,
                               endsample,
                               _noiseCut,
                               _parabnbefmax,
                               _parabnaftmax,
                               _thres,
                               _lowlev,
                               _highlev,
                               _nevlasers,
                               _slide);

    if (_debug == 1)
      edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- inside loop 2  -- ";

    // analyse the Matacq data
    if (mtq->rawPulseAnalysis(nsamples, &matacq[0]) == 0) {
      status = 1;
      ped = mtq->getBaseLine();
      pedsig = mtq->getsigBaseLine();

      if (_debug == 1)
        edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- inside loop 3  -- ped:" << ped;
      if (mtq->findPeak() == 0) {
        peak = mtq->getTimpeak();
        sigma = mtq->getsigTimpeak();
      }
      if (_debug == 1)
        edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- inside loop 4  -- peak:" << peak;
      if (mtq->doFit() == 0) {
        fit = mtq->getTimax();
        ampl = mtq->getAmpl();
        fwhm = mtq->getFwhm();
        fw20 = mtq->getWidth20();
        fw80 = mtq->getWidth80();
        sliding = mtq->getSlide();
      }
      if (_debug == 1)
        edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- inside loop 4  -- ampl:" << ampl;
      if (mtq->compute_trise() == 0) {
        trise = mtq->getTrise();
      }
      if (_debug == 1)
        edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- inside loop 5  -- trise:" << trise;
    }

    if (_debug == 1)
      edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- inside loop 6  -- status:" << status;

    if (status == 1 && mtq->findPeak() == 0) {
      int firstS = int(peak - double(_timebefmax));
      int lastS = int(peak + double(_timeaftmax));

      // Fill histo if there are enough samples
      if (_debug == 1)
        edm::LogVerbatim("EcalMatacqAnalyzzer")
            << "-- debug test -- inside loop 7  -- firstS:" << firstS << ", nsamples:" << nsamples;

      if (firstS >= 0 && lastS <= nsamples) {
        for (int i = firstS; i < lastS; i++) {
          shapeMatTmp->Fill(double(i) - firstS, matacq[i]);
        }

      } else {  // else extrapolate

        int firstSBis;

        if (firstS < 0) {  // fill first bins with 0

          double thisped;
          thisped = (matacq[0] + matacq[1] + matacq[2] + matacq[4] + matacq[5]) / 5.0;

          for (int i = firstS; i < 0; i++) {
            shapeMatTmp->Fill(double(i) - firstS, thisped);
          }
          firstSBis = 0;

        } else {
          firstSBis = firstS;
        }

        if (lastS > nsamples) {
          for (int i = firstSBis; i < int(nsamples); i++) {
            shapeMatTmp->Fill(double(i) - firstS, matacq[i]);
          }

          //extrapolate with expo tail

          double expb = 0.998;
          double matacqval = expb * matacq[nsamples - 1];

          for (int i = nsamples; i < lastS; i++) {
            shapeMatTmp->Fill(double(i) - firstS, matacqval);
            matacqval *= expb;
          }

        } else {
          for (int i = firstSBis; i < lastS; i++) {
            shapeMatTmp->Fill(double(i) - firstS, matacq[i]);
          }
        }
      }
    }
    if (_debug == 1)
      edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- inside loop 8";

    // get back color

    int iCol = nCol;
    for (unsigned int i = 0; i < nCol; i++) {
      if (color == colors[i]) {
        iCol = i;
        i = nCol;
      }
    }
    if (_debug == 1)
      edm::LogVerbatim("EcalMatacqAnalyzzer")
          << "-- debug test -- inside loop 8bis color:" << color << " iCol:" << iCol << " nCol:" << nCol;

    // fill TMTQ

    if (status == 1 && mtq->findPeak() == 0 && mtq->doFit() == 0 && mtq->compute_trise() == 0)
      MTQ[iCol][lightside]->addEntry(peak, sigma, fit, ampl, trise, fwhm, fw20, fw80, ped, pedsig, sliding);

    // fill the output tree

    if (_debug == 1)
      edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- inside loop 9";
    mtqShape->Fill();

    // clean up
    delete mtq;
  }

  if (_debug == 1)
    edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- after loop ";
  sampFile->Close();

  double Peak[6], Sigma[6], Fit[6], Ampl[6], Trise[6], Fwhm[6], Fw20[6], Fw80[6], Ped[6], Pedsig[6], Sliding[6];
  int Side;

  for (unsigned int iColor = 0; iColor < nCol; iColor++) {
    std::stringstream nametree;
    nametree << "MatacqCol" << colors[iColor];
    meanTree[iColor] = new TTree(nametree.str().c_str(), nametree.str().c_str());
    meanTree[iColor]->Branch("side", &Side, "Side/I");
    meanTree[iColor]->Branch("peak", &Peak, "Peak[6]/D");
    meanTree[iColor]->Branch("sigma", &Sigma, "Sigma[6]/D");
    meanTree[iColor]->Branch("fit", &Fit, "Fit[6]/D");
    meanTree[iColor]->Branch("ampl", &Ampl, "Ampl[6]/D");
    meanTree[iColor]->Branch("trise", &Trise, "Trise[6]/D");
    meanTree[iColor]->Branch("fwhm", &Fwhm, "Fwhm[6]/D");
    meanTree[iColor]->Branch("fw20", &Fw20, "Fw20[6]/D");
    meanTree[iColor]->Branch("fw80", &Fw80, "Fw80[6]/D");
    meanTree[iColor]->Branch("ped", &Ped, "Ped[6]/D");
    meanTree[iColor]->Branch("pedsig", &Pedsig, "Pedsig[6]/D");
    meanTree[iColor]->Branch("sliding", &Sliding, "Sliding[6]/D");

    meanTree[iColor]->SetBranchAddress("side", &Side);
    meanTree[iColor]->SetBranchAddress("peak", Peak);
    meanTree[iColor]->SetBranchAddress("sigma", Sigma);
    meanTree[iColor]->SetBranchAddress("fit", Fit);
    meanTree[iColor]->SetBranchAddress("ampl", Ampl);
    meanTree[iColor]->SetBranchAddress("fwhm", Fwhm);
    meanTree[iColor]->SetBranchAddress("fw20", Fw20);
    meanTree[iColor]->SetBranchAddress("fw80", Fw80);
    meanTree[iColor]->SetBranchAddress("trise", Trise);
    meanTree[iColor]->SetBranchAddress("ped", Ped);
    meanTree[iColor]->SetBranchAddress("pedsig", Pedsig);
    meanTree[iColor]->SetBranchAddress("sliding", Sliding);
  }

  for (unsigned int iCol = 0; iCol < nCol; iCol++) {
    for (unsigned int iSide = 0; iSide < nSides; iSide++) {
      Side = iSide;
      std::vector<double> val[TMTQ::nOutVar];

      for (int iVar = 0; iVar < TMTQ::nOutVar; iVar++) {
        val[iVar] = MTQ[iCol][iSide]->get(iVar);

        for (unsigned int i = 0; i < val[iVar].size(); i++) {
          switch (iVar) {
            case TMTQ::iPeak:
              Peak[i] = val[iVar].at(i);
              break;
            case TMTQ::iSigma:
              Sigma[i] = val[iVar].at(i);
              break;
            case TMTQ::iFit:
              Fit[i] = val[iVar].at(i);
              break;
            case TMTQ::iAmpl:
              Ampl[i] = val[iVar].at(i);
              break;
            case TMTQ::iFwhm:
              Fwhm[i] = val[iVar].at(i);
              break;
            case TMTQ::iFw20:
              Fw20[i] = val[iVar].at(i);
              break;
            case TMTQ::iFw80:
              Fw80[i] = val[iVar].at(i);
              break;
            case TMTQ::iTrise:
              Trise[i] = val[iVar].at(i);
              break;
            case TMTQ::iPed:
              Ped[i] = val[iVar].at(i);
              break;
            case TMTQ::iPedsig:
              Pedsig[i] = val[iVar].at(i);
              break;
            case TMTQ::iSlide:
              Sliding[i] = val[iVar].at(i);
              break;
          }
        }
      }
      meanTree[iCol]->Fill();
      if (_debug == 1)
        edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- inside final loop  ";
    }
  }

  // Calculate maximum with pol 2

  int im = shapeMatTmp->GetMaximumBin();
  double q1 = shapeMatTmp->GetBinContent(im - 1);
  double q2 = shapeMatTmp->GetBinContent(im);
  double q3 = shapeMatTmp->GetBinContent(im + 1);

  double a2 = (q3 + q1) / 2.0 - q2;
  double a1 = q2 - q1 + a2 * (1 - 2 * im);
  double a0 = q2 - a1 * im - a2 * im * im;

  double am = a0 - a1 * a1 / (4 * a2);

  // Compute pedestal

  double bl = 0;
  for (unsigned int i = 1; i < _presampleshape + 1; i++) {
    bl += shapeMatTmp->GetBinContent(i);
  }
  bl /= _presampleshape;

  // Compute and save laser shape

  if (_debug == 1)
    edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- computing shape  ";

  int firstBin = 0;
  double height = 0.0;

  for (unsigned int i = _timebefmax; i > _presampleshape; i--) {
    height = shapeMatTmp->GetBinContent(i) - bl;

    if (height < (am - bl) * _cutwindow) {
      firstBin = i;
      i = _presampleshape;
    }
  }

  unsigned int lastBin = firstBin + _nsamplesshape;

  for (unsigned int i = firstBin; i < lastBin; i++) {
    shapeMat->Fill(i - firstBin, shapeMatTmp->GetBinContent(i) - bl);
  }

  mtqShape->Write();
  for (unsigned int iColor = 0; iColor < nCol; iColor++) {
    meanTree[iColor]->Write();
  }
  if (_debug == 1)
    edm::LogVerbatim("EcalMatacqAnalyzzer") << "-- debug test -- writing  ";
  shapeMat->Write();

  // close the output file
  outFile->Close();

  // Remove temporary file
  FILE* test;
  test = fopen(sampfile.c_str(), "r");
  if (test) {
    std::stringstream del2;
    del2 << "rm " << sampfile;
    system(del2.str().c_str());
  }

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

Member Data Documentation

_cutwindow

const double EcalMatacqAnalyzer::_cutwindow

private

Definition at line 46 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_debug

const int EcalMatacqAnalyzer::_debug

private

Definition at line 51 of file EcalMatacqAnalyzer.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().

_fedid

const int EcalMatacqAnalyzer::_fedid

private

Definition at line 50 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

_highlev

const unsigned int EcalMatacqAnalyzer::_highlev

private

Definition at line 42 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_lowlev

const unsigned int EcalMatacqAnalyzer::_lowlev

private

Definition at line 41 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_nevlasers

const unsigned int EcalMatacqAnalyzer::_nevlasers

private

Definition at line 43 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_noiseCut

const unsigned int EcalMatacqAnalyzer::_noiseCut

private

Definition at line 37 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_nsamplesaftmax

const unsigned int EcalMatacqAnalyzer::_nsamplesaftmax

private

Definition at line 36 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_nsamplesshape

const unsigned int EcalMatacqAnalyzer::_nsamplesshape

private

Definition at line 47 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_parabnaftmax

const unsigned int EcalMatacqAnalyzer::_parabnaftmax

private

Definition at line 39 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_parabnbefmax

const unsigned int EcalMatacqAnalyzer::_parabnbefmax

private

Definition at line 38 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_presample

const double EcalMatacqAnalyzer::_presample

private

Definition at line 35 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_presampleshape

const unsigned int EcalMatacqAnalyzer::_presampleshape

private

Definition at line 48 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_slide

const unsigned int EcalMatacqAnalyzer::_slide

private

Definition at line 49 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_thres

const unsigned int EcalMatacqAnalyzer::_thres

private

Definition at line 40 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_timeaftmax

const unsigned int EcalMatacqAnalyzer::_timeaftmax

private

Definition at line 45 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

_timebefmax

const unsigned int EcalMatacqAnalyzer::_timebefmax

private

Definition at line 44 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

ampl

double EcalMatacqAnalyzer::ampl

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

color

int EcalMatacqAnalyzer::color

private

Definition at line 80 of file EcalMatacqAnalyzer.h.

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

colors

std::vector<int> EcalMatacqAnalyzer::colors

private

Definition at line 92 of file EcalMatacqAnalyzer.h.

dccID

int EcalMatacqAnalyzer::dccID

private

Definition at line 70 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

dccToken_

const edm::EDGetTokenT<EcalRawDataCollection> EcalMatacqAnalyzer::dccToken_

private

Definition at line 59 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

digiCollection_

const std::string EcalMatacqAnalyzer::digiCollection_

private

Definition at line 54 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

digiProducer_

const std::string EcalMatacqAnalyzer::digiProducer_

private

Definition at line 55 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

event

int EcalMatacqAnalyzer::event

private

Definition at line 79 of file EcalMatacqAnalyzer.h.

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

eventHeaderCollection_

const std::string EcalMatacqAnalyzer::eventHeaderCollection_

private

Definition at line 56 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

eventHeaderProducer_

const std::string EcalMatacqAnalyzer::eventHeaderProducer_

private

Definition at line 57 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

fedID

int EcalMatacqAnalyzer::fedID

private

Definition at line 71 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

fit

double EcalMatacqAnalyzer::fit

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by trackingPlots.Iteration::modules().

fw20

double EcalMatacqAnalyzer::fw20

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

fw80

double EcalMatacqAnalyzer::fw80

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

fwhm

double EcalMatacqAnalyzer::fwhm

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

iEvent

int EcalMatacqAnalyzer::iEvent

private

Definition at line 29 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

isThereMatacq

bool EcalMatacqAnalyzer::isThereMatacq

private

Definition at line 75 of file EcalMatacqAnalyzer.h.

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

laserEvents

int EcalMatacqAnalyzer::laserEvents

private

Definition at line 74 of file EcalMatacqAnalyzer.h.

Referenced by analyze(), and beginJob().

lightside

int EcalMatacqAnalyzer::lightside

private

Definition at line 67 of file EcalMatacqAnalyzer.h.

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

matacq

double EcalMatacqAnalyzer::matacq[2560]

private

Definition at line 81 of file EcalMatacqAnalyzer.h.

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

maxsamp

int EcalMatacqAnalyzer::maxsamp

private

Definition at line 82 of file EcalMatacqAnalyzer.h.

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

meanTree

TTree* EcalMatacqAnalyzer::meanTree[nColor]

private

Definition at line 90 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

MTQ

TMTQ* EcalMatacqAnalyzer::MTQ[nColor][nSide]

private

Definition at line 89 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

mtqShape

TTree* EcalMatacqAnalyzer::mtqShape

private

Definition at line 97 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

nsamples

int EcalMatacqAnalyzer::nsamples

private

Definition at line 83 of file EcalMatacqAnalyzer.h.

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

nSides

unsigned int EcalMatacqAnalyzer::nSides

private

Definition at line 66 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

outfile

std::string EcalMatacqAnalyzer::outfile

private

Definition at line 61 of file EcalMatacqAnalyzer.h.

Referenced by beginJob(), and endJob().

outFile

TFile* EcalMatacqAnalyzer::outFile

private

Definition at line 94 of file EcalMatacqAnalyzer.h.

Referenced by Formatter.SimpleHTMLFormatter::__del__(), endJob(), Formatter.SimpleHTMLFormatter::endTable(), Formatter.SimpleHTMLFormatter::headers(), Formatter.SimpleHTMLFormatter::showLine(), Formatter.SimpleHTMLFormatter::startTable(), Formatter.SimpleHTMLFormatter::trailers(), Formatter.SimpleHTMLFormatter::write(), Formatter.SimpleHTMLFormatter::writeAnchor(), Formatter.SimpleHTMLFormatter::writeB(), Formatter.SimpleHTMLFormatter::writeBr(), Formatter.SimpleHTMLFormatter::writeH1(), Formatter.SimpleHTMLFormatter::writeH2(), Formatter.SimpleHTMLFormatter::writeH3(), Formatter.SimpleHTMLFormatter::writeRow(), and Formatter.SimpleHTMLFormatter::writeStyledRow().

peak

double EcalMatacqAnalyzer::peak

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

ped

double EcalMatacqAnalyzer::ped

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

pedsig

double EcalMatacqAnalyzer::pedsig

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

pmatToken_

const edm::EDGetTokenT<EcalMatacqDigiCollection> EcalMatacqAnalyzer::pmatToken_

private

Definition at line 58 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

resdir_

const std::string EcalMatacqAnalyzer::resdir_

private

Definition at line 53 of file EcalMatacqAnalyzer.h.

Referenced by beginJob().

runNum

int EcalMatacqAnalyzer::runNum

private

Definition at line 69 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

runType

int EcalMatacqAnalyzer::runType

private

Definition at line 68 of file EcalMatacqAnalyzer.h.

Referenced by analyze().

sampfile

std::string EcalMatacqAnalyzer::sampfile

private

Definition at line 62 of file EcalMatacqAnalyzer.h.

Referenced by beginJob(), and endJob().

sampFile

TFile* EcalMatacqAnalyzer::sampFile

private

Definition at line 86 of file EcalMatacqAnalyzer.h.

Referenced by beginJob(), and endJob().

sigma

double EcalMatacqAnalyzer::sigma

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

sliding

double EcalMatacqAnalyzer::sliding

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

status

int EcalMatacqAnalyzer::status

private

Definition at line 95 of file EcalMatacqAnalyzer.h.

Referenced by dirstructure.Comparison::__make_image(), dirstructure.Comparison::__repr__(), and endJob().

tree

TTree* EcalMatacqAnalyzer::tree

private

Definition at line 87 of file EcalMatacqAnalyzer.h.

Referenced by eventstfile.Events::__iter__(), svgfig.SVG::__str__(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::beginLoop(), tree.Tree::branch_(), MetTreeProducer.MetTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::fillTree(), SimpleTreeProducer.SimpleTreeProducer::process(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::process(), MetTreeProducer.MetTreeProducer::process(), python.cmstools.EventTree::SetAlias(), eventstfile.Events::size(), eventstfile.Events::to(), tree.Tree::var(), and tree.Tree::vector().

trise

double EcalMatacqAnalyzer::trise

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

tt

double EcalMatacqAnalyzer::tt

private

Definition at line 84 of file EcalMatacqAnalyzer.h.

Referenced by endJob().

ttrig

double EcalMatacqAnalyzer::ttrig

private

Definition at line 96 of file EcalMatacqAnalyzer.h.

Referenced by endJob().