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 final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const 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
bool	wantsStreamRuns () const
	~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(), laser_color, EcalDCCHeaderBlock::LASER_DELAY_SCAN, EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_POWER_SCAN, EcalDCCHeaderBlock::LASER_STD, matacq, SiStripPI::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().