EcalPerEvtMatacqAnalyzer Class Reference

#include <EcalPerEvtMatacqAnalyzer.h>

Inheritance diagram for EcalPerEvtMatacqAnalyzer:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	beginJob () override
	EcalPerEvtMatacqAnalyzer (const edm::ParameterSet &iConfig)
void	endJob () override
	~EcalPerEvtMatacqAnalyzer () 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 unsigned int	_highlev
const unsigned int	_lowlev
const unsigned int	_nevlasers
const unsigned int	_noiseCut
const unsigned int	_nsamplesaftmax
const unsigned int	_nsamplesbefmax
const unsigned int	_parabnaftmax
const unsigned int	_parabnbefmax
const double	_presample
const unsigned int	_thres
double	ampl
const edm::EDGetTokenT< EcalRawDataCollection >	dccToken_
const std::string	digiCollection_
const std::string	digiProducer_
int	event
const std::string	eventHeaderCollection_
const std::string	eventHeaderProducer_
double	fit
int	iEvent
int	IsFileCreated
int	IsTreeCreated
int	laser_color
double	matacq [2560]
unsigned int	maxsamp
TTree *	mtqShape
unsigned int	nsamples
std::string	outfile
TFile *	outFile
double	peak
const edm::EDGetTokenT< EcalMatacqDigiCollection >	pmatToken_
const std::string	resdir_
int	runNum
int	runType
TFile *	sampFile
double	sigma
int	status
TTree *	tree
double	trise
double	ttrig

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
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 14 of file EcalPerEvtMatacqAnalyzer.h.

Constructor & Destructor Documentation

EcalPerEvtMatacqAnalyzer()

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

explicit

Definition at line 29 of file EcalPerEvtMatacqAnalyzer.cc.

    : iEvent(0),
      // framework parameters with default values

      _presample(iConfig.getUntrackedParameter<double>("nPresamples", 6.7)),
      _nsamplesaftmax(iConfig.getUntrackedParameter<unsigned int>("nSamplesAftMax", 80)),
      _nsamplesbefmax(iConfig.getUntrackedParameter<unsigned int>("nSamplesBefMax", 20)),
      _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)),
      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(digiProducer_)))

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

~EcalPerEvtMatacqAnalyzer()

EcalPerEvtMatacqAnalyzer::~EcalPerEvtMatacqAnalyzer ( )

overridedefault

Member Function Documentation

analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 67 of file EcalPerEvtMatacqAnalyzer.cc.

References EcalMatacqDigi::adcCount(), edm::SortedCollection< T, SORT >::begin(), dccToken_, digiCollection_, MillePedeFileConverter_cfg::e, edm::SortedCollection< T, SORT >::end(), eventHeaderCollection_, mps_fire::i, iEvent, createfilelist::int, IsFileCreated, IsTreeCreated, edm::HandleBase::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, laser_color, EcalDCCHeaderBlock::LASER_DELAY_SCAN, EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_POWER_SCAN, EcalDCCHeaderBlock::LASER_STD, matacq, WZElectronSkims53X_cff::max, maxsamp, N_channels, nsamples, outfile, outFile, pmatToken_, edm::Handle< T >::product(), resdir_, runNum, runType, EcalMatacqDigi::size(), ttrig, EcalMatacqDigi::tTrig(), and EcalDCCHeaderBlock::EcalDCCEventSettings::wavelength.

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

  ++iEvent;

  // retrieving MATACQ :
  const edm::Handle<EcalMatacqDigiCollection>& pmatacqDigi = e.getHandle(pmatToken_);
  const EcalMatacqDigiCollection* matacqDigi = (pmatacqDigi.isValid()) ? pmatacqDigi.product() : nullptr;
  if (!(pmatacqDigi.isValid()))
    edm::LogWarning("EcalPerEvtMatacqAnalyzer") << "Error! can't get the product " << digiCollection_.c_str();

  // retrieving DCC header

  const edm::Handle<EcalRawDataCollection>& pDCCHeader = e.getHandle(dccToken_);
  const EcalRawDataCollection* DCCHeader = (pDCCHeader.isValid()) ? pDCCHeader.product() : nullptr;
  if (!pDCCHeader.isValid())
    edm::LogWarning("EcalPerEvtMatacqAnalyzer") << "Error! can't get the product " << eventHeaderCollection_.c_str();

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

  for (EcalRawDataCollection::const_iterator headerItr = DCCHeader->begin(); headerItr != DCCHeader->end();
       ++headerItr) {
    event = headerItr->getLV1();
    EcalDCCHeaderBlock::EcalDCCEventSettings settings = headerItr->getEventSettings();

    laser_color = (int)settings.wavelength;

    // Get run type and run number

    runType = headerItr->getRunType();
    runNum = headerItr->getRunNumber();

    // Cut on runType

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

    // Define output results files' names

    if (IsFileCreated == 0) {
      std::stringstream namefile;

      namefile << resdir_ << "/MatacqPerEvt-Run" << runNum << ".root";
      outfile = namefile.str();

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

      // declaration of the tree to fill

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

      IsFileCreated = 1;
    }
  }

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

  int iCh = 0;
  double max = 0;

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

    //
    const EcalMatacqDigi& digis = *it;

    if (digis.size() == 0 || iCh >= N_channels)
      continue;

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

    if (IsTreeCreated == 0) {
      //List of branches

      std::stringstream mat;
      mat << "matacq[" << nsamples << "]/D" << std::endl;

      tree->Branch("event", &event, "event/I");
      tree->Branch("laser_color", &laser_color, "laser_color/I");
      //      tree->Branch( "nsamples",    &nsamples ,    "nsamples/I"       );
      tree->Branch("matacq", &matacq, mat.str().c_str());
      tree->Branch("maxsamp", &maxsamp, "maxsamp/I");
      tree->Branch("ttrig", &ttrig, "ttrig/D");

      tree->SetBranchAddress("event", &event);
      tree->SetBranchAddress("laser_color", &laser_color);
      tree->SetBranchAddress("matacq", matacq);
      //      tree->SetBranchAddress( "nsamples",    &nsamples    );
      tree->SetBranchAddress("maxsamp", &maxsamp);
      tree->SetBranchAddress("ttrig", &ttrig);

      IsTreeCreated = 1;
    }

    ttrig = 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;
      }
    }

    iCh++;
  }

  tree->Fill();

}  // analyze

beginJob()

void EcalPerEvtMatacqAnalyzer::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 57 of file EcalPerEvtMatacqAnalyzer.cc.

References IsFileCreated, and IsTreeCreated.

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

  //Initializations

  IsFileCreated = 0;
  IsTreeCreated = 0;
}

endJob()

void EcalPerEvtMatacqAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 186 of file EcalPerEvtMatacqAnalyzer.cc.

References outFile.

Referenced by o2olib.O2ORunMgr::executeJob().

                                      {
  tree->Write();
  outFile->Close();
}

Member Data Documentation

_highlev

const unsigned int EcalPerEvtMatacqAnalyzer::_highlev

private

Definition at line 64 of file EcalPerEvtMatacqAnalyzer.h.

_lowlev

const unsigned int EcalPerEvtMatacqAnalyzer::_lowlev

private

Definition at line 63 of file EcalPerEvtMatacqAnalyzer.h.

_nevlasers

const unsigned int EcalPerEvtMatacqAnalyzer::_nevlasers

private

Definition at line 65 of file EcalPerEvtMatacqAnalyzer.h.

_noiseCut

const unsigned int EcalPerEvtMatacqAnalyzer::_noiseCut

private

Definition at line 59 of file EcalPerEvtMatacqAnalyzer.h.

_nsamplesaftmax

const unsigned int EcalPerEvtMatacqAnalyzer::_nsamplesaftmax

private

Definition at line 57 of file EcalPerEvtMatacqAnalyzer.h.

_nsamplesbefmax

const unsigned int EcalPerEvtMatacqAnalyzer::_nsamplesbefmax

private

Definition at line 58 of file EcalPerEvtMatacqAnalyzer.h.

_parabnaftmax

const unsigned int EcalPerEvtMatacqAnalyzer::_parabnaftmax

private

Definition at line 61 of file EcalPerEvtMatacqAnalyzer.h.

_parabnbefmax

const unsigned int EcalPerEvtMatacqAnalyzer::_parabnbefmax

private

Definition at line 60 of file EcalPerEvtMatacqAnalyzer.h.

_presample

const double EcalPerEvtMatacqAnalyzer::_presample

private

Definition at line 56 of file EcalPerEvtMatacqAnalyzer.h.

_thres

const unsigned int EcalPerEvtMatacqAnalyzer::_thres

private

Definition at line 62 of file EcalPerEvtMatacqAnalyzer.h.

ampl

double EcalPerEvtMatacqAnalyzer::ampl

private

Definition at line 48 of file EcalPerEvtMatacqAnalyzer.h.

dccToken_

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

private

Definition at line 74 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

digiCollection_

const std::string EcalPerEvtMatacqAnalyzer::digiCollection_

private

Definition at line 68 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

digiProducer_

const std::string EcalPerEvtMatacqAnalyzer::digiProducer_

private

Definition at line 69 of file EcalPerEvtMatacqAnalyzer.h.

event

int EcalPerEvtMatacqAnalyzer::event

private

Definition at line 35 of file EcalPerEvtMatacqAnalyzer.h.

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

eventHeaderCollection_

const std::string EcalPerEvtMatacqAnalyzer::eventHeaderCollection_

private

Definition at line 70 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

eventHeaderProducer_

const std::string EcalPerEvtMatacqAnalyzer::eventHeaderProducer_

private

Definition at line 71 of file EcalPerEvtMatacqAnalyzer.h.

fit

double EcalPerEvtMatacqAnalyzer::fit

private

Definition at line 48 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by trackingPlots.Iteration::modules().

iEvent

int EcalPerEvtMatacqAnalyzer::iEvent

private

Definition at line 26 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

IsFileCreated

int EcalPerEvtMatacqAnalyzer::IsFileCreated

private

Definition at line 44 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze(), and beginJob().

IsTreeCreated

int EcalPerEvtMatacqAnalyzer::IsTreeCreated

private

Definition at line 45 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze(), and beginJob().

laser_color

int EcalPerEvtMatacqAnalyzer::laser_color

private

Definition at line 36 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

matacq

double EcalPerEvtMatacqAnalyzer::matacq[2560]

private

Definition at line 37 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

maxsamp

unsigned int EcalPerEvtMatacqAnalyzer::maxsamp

private

Definition at line 38 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

mtqShape

TTree* EcalPerEvtMatacqAnalyzer::mtqShape

private

Definition at line 49 of file EcalPerEvtMatacqAnalyzer.h.

nsamples

unsigned int EcalPerEvtMatacqAnalyzer::nsamples

private

Definition at line 39 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

outfile

std::string EcalPerEvtMatacqAnalyzer::outfile

private

Definition at line 24 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

outFile

TFile* EcalPerEvtMatacqAnalyzer::outFile

private

Definition at line 46 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by Formatter.SimpleHTMLFormatter::__del__(), analyze(), 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 EcalPerEvtMatacqAnalyzer::peak

private

Definition at line 48 of file EcalPerEvtMatacqAnalyzer.h.

pmatToken_

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

private

Definition at line 73 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

resdir_

const std::string EcalPerEvtMatacqAnalyzer::resdir_

private

Definition at line 67 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

runNum

int EcalPerEvtMatacqAnalyzer::runNum

private

Definition at line 31 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

runType

int EcalPerEvtMatacqAnalyzer::runType

private

Definition at line 30 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

sampFile

TFile* EcalPerEvtMatacqAnalyzer::sampFile

private

Definition at line 41 of file EcalPerEvtMatacqAnalyzer.h.

sigma

double EcalPerEvtMatacqAnalyzer::sigma

private

Definition at line 48 of file EcalPerEvtMatacqAnalyzer.h.

status

int EcalPerEvtMatacqAnalyzer::status

private

Definition at line 47 of file EcalPerEvtMatacqAnalyzer.h.

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

tree

TTree* EcalPerEvtMatacqAnalyzer::tree

private

Definition at line 42 of file EcalPerEvtMatacqAnalyzer.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 EcalPerEvtMatacqAnalyzer::trise

private

Definition at line 48 of file EcalPerEvtMatacqAnalyzer.h.

ttrig

double EcalPerEvtMatacqAnalyzer::ttrig

private

Definition at line 48 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().