#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

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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
unsigned int	_highlev

unsigned int	_lowlev

unsigned int	_nevlasers

unsigned int	_noiseCut

unsigned int	_nsamplesaftmax

unsigned int	_nsamplesbefmax

unsigned int	_parabnaftmax

unsigned int	_parabnbefmax

double	_presample

unsigned int	_thres

double	ampl

std::string	digiCollection_

std::string	digiProducer_

int	event

std::string	eventHeaderCollection_

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

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::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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)

Detailed Description

Definition at line 15 of file EcalPerEvtMatacqAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 34 of file EcalPerEvtMatacqAnalyzer.cc.

References digiCollection_, digiProducer_, eventHeaderCollection_, eventHeaderProducer_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), resdir_, and AlCaHLTBitMon_QueryRunRegistry::string.

   :
   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 ) )
 
 //========================================================================
 {
 
 
   //now do what ever initialization is needed
 
   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");
 
 }

EcalPerEvtMatacqAnalyzer::~EcalPerEvtMatacqAnalyzer ( )

override

Definition at line 67 of file EcalPerEvtMatacqAnalyzer.cc.

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

Member Function Documentation

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

override

Definition at line 90 of file EcalPerEvtMatacqAnalyzer.cc.

References EcalMatacqDigi::adcCount(), edm::SortedCollection< T, SORT >::begin(), MessageLogger_cfi::cerr, digiCollection_, digiProducer_, edm::SortedCollection< T, SORT >::end(), eventHeaderCollection_, cppFunctionSkipper::exception, edm::Event::getByLabel(), mps_fire::i, iEvent, createfilelist::int, IsFileCreated, IsTreeCreated, 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, edm::Handle< T >::product(), resdir_, runNum, runType, EcalMatacqDigi::size(), ttrig, EcalMatacqDigi::tTrig(), and EcalDCCHeaderBlock::EcalDCCEventSettings::wavelength.

                                                                                    {
 //========================================================================
 
   ++iEvent;
 
   // retrieving MATACQ :
   edm::Handle<EcalMatacqDigiCollection> pmatacqDigi;
   const EcalMatacqDigiCollection* matacqDigi=nullptr;
   try {
     e.getByLabel(digiProducer_,digiCollection_, pmatacqDigi); 
     matacqDigi=pmatacqDigi.product();
   }catch ( std::exception& ex ) {
     std::cerr << "Error! can't get the product " << digiCollection_.c_str() << std::endl;
 
   }
   
   // retrieving DCC header
 
   edm::Handle<EcalRawDataCollection> pDCCHeader;
   const  EcalRawDataCollection* DCCHeader=nullptr;
   try {
      e.getByLabel(digiProducer_, pDCCHeader); 
      //e.getByLabel(eventHeaderProducer_,eventHeaderCollection_, pDCCHeader);
     DCCHeader=pDCCHeader.product();
   }catch ( std::exception& ex ) {
     std::cerr << "Error! can't get the product " << eventHeaderCollection_.c_str() << std::endl;
   }
 
 
   // ====================================
   // 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){
 
       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
       
       stringstream mat;
       mat << "matacq["<<nsamples<<"]/D"<< 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

void EcalPerEvtMatacqAnalyzer::beginJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 78 of file EcalPerEvtMatacqAnalyzer.cc.

References IsFileCreated, and IsTreeCreated.

                                         {
 //========================================================================
 
     //Initializations
 
     IsFileCreated=0;
     IsTreeCreated=0;
 
 }

void EcalPerEvtMatacqAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 224 of file EcalPerEvtMatacqAnalyzer.cc.

References DEFINE_FWK_MODULE, and outFile.

Referenced by o2olib.O2ORunMgr::executeJob().

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

Member Data Documentation

unsigned int EcalPerEvtMatacqAnalyzer::_highlev

private

Definition at line 77 of file EcalPerEvtMatacqAnalyzer.h.

unsigned int EcalPerEvtMatacqAnalyzer::_lowlev

private

Definition at line 76 of file EcalPerEvtMatacqAnalyzer.h.

unsigned int EcalPerEvtMatacqAnalyzer::_nevlasers

private

Definition at line 78 of file EcalPerEvtMatacqAnalyzer.h.

unsigned int EcalPerEvtMatacqAnalyzer::_noiseCut

private

Definition at line 72 of file EcalPerEvtMatacqAnalyzer.h.

unsigned int EcalPerEvtMatacqAnalyzer::_nsamplesaftmax

private

Definition at line 70 of file EcalPerEvtMatacqAnalyzer.h.

unsigned int EcalPerEvtMatacqAnalyzer::_nsamplesbefmax

private

Definition at line 71 of file EcalPerEvtMatacqAnalyzer.h.

unsigned int EcalPerEvtMatacqAnalyzer::_parabnaftmax

private

Definition at line 74 of file EcalPerEvtMatacqAnalyzer.h.

unsigned int EcalPerEvtMatacqAnalyzer::_parabnbefmax

private

Definition at line 73 of file EcalPerEvtMatacqAnalyzer.h.

double EcalPerEvtMatacqAnalyzer::_presample

private

Definition at line 69 of file EcalPerEvtMatacqAnalyzer.h.

unsigned int EcalPerEvtMatacqAnalyzer::_thres

private

Definition at line 75 of file EcalPerEvtMatacqAnalyzer.h.

double EcalPerEvtMatacqAnalyzer::ampl

private

Definition at line 60 of file EcalPerEvtMatacqAnalyzer.h.

std::string EcalPerEvtMatacqAnalyzer::digiCollection_

private

Definition at line 31 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze(), and EcalPerEvtMatacqAnalyzer().

std::string EcalPerEvtMatacqAnalyzer::digiProducer_

private

Definition at line 32 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze(), and EcalPerEvtMatacqAnalyzer().

int EcalPerEvtMatacqAnalyzer::event

private

Definition at line 47 of file EcalPerEvtMatacqAnalyzer.h.

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

std::string EcalPerEvtMatacqAnalyzer::eventHeaderCollection_

private

Definition at line 33 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze(), and EcalPerEvtMatacqAnalyzer().

std::string EcalPerEvtMatacqAnalyzer::eventHeaderProducer_

private

Definition at line 34 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by EcalPerEvtMatacqAnalyzer().

double EcalPerEvtMatacqAnalyzer::fit

private

Definition at line 60 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by trackingPlots.Iteration::modules().

int EcalPerEvtMatacqAnalyzer::iEvent

private

Definition at line 38 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

int EcalPerEvtMatacqAnalyzer::IsFileCreated

private

Definition at line 56 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze(), and beginJob().

int EcalPerEvtMatacqAnalyzer::IsTreeCreated

private

Definition at line 57 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze(), and beginJob().

int EcalPerEvtMatacqAnalyzer::laser_color

private

Definition at line 48 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

double EcalPerEvtMatacqAnalyzer::matacq[2560]

private

Definition at line 49 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

unsigned int EcalPerEvtMatacqAnalyzer::maxsamp

private

Definition at line 50 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

TTree* EcalPerEvtMatacqAnalyzer::mtqShape

private

Definition at line 61 of file EcalPerEvtMatacqAnalyzer.h.

unsigned int EcalPerEvtMatacqAnalyzer::nsamples

private

Definition at line 51 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

std::string EcalPerEvtMatacqAnalyzer::outfile

private

Definition at line 36 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

TFile* EcalPerEvtMatacqAnalyzer::outFile

private

Definition at line 58 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze(), and endJob().

double EcalPerEvtMatacqAnalyzer::peak

private

Definition at line 60 of file EcalPerEvtMatacqAnalyzer.h.

std::string EcalPerEvtMatacqAnalyzer::resdir_

private

Definition at line 30 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze(), and EcalPerEvtMatacqAnalyzer().

int EcalPerEvtMatacqAnalyzer::runNum

private

Definition at line 43 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

int EcalPerEvtMatacqAnalyzer::runType

private

Definition at line 42 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

TFile* EcalPerEvtMatacqAnalyzer::sampFile

private

Definition at line 53 of file EcalPerEvtMatacqAnalyzer.h.

double EcalPerEvtMatacqAnalyzer::sigma

private

Definition at line 60 of file EcalPerEvtMatacqAnalyzer.h.

int EcalPerEvtMatacqAnalyzer::status

private

Definition at line 59 of file EcalPerEvtMatacqAnalyzer.h.

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

TTree* EcalPerEvtMatacqAnalyzer::tree

private

Definition at line 54 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by eventstfile.Events::__iter__(), svgfig.SVG::__str__(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::beginLoop(), MetTreeProducer.MetTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::fillTree(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::process(), MetTreeProducer.MetTreeProducer::process(), python.cmstools.EventTree::SetAlias(), and eventstfile.Events::to().

double EcalPerEvtMatacqAnalyzer::trise

private

Definition at line 60 of file EcalPerEvtMatacqAnalyzer.h.

double EcalPerEvtMatacqAnalyzer::ttrig

private

Definition at line 60 of file EcalPerEvtMatacqAnalyzer.h.

Referenced by analyze().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation