calculate the best DAC value to obtain a pedestal = 200 More...

#include <testChannel.h>

Inheritance diagram for testChannel:

Public Member Functions
void	analyze (edm::Event const &event, edm::EventSetup const &eventSetup) override
	! Analyze More...

void	beginJob () override
	BeginJob. More...

void	endJob (void) override
	EndJob. More...

void	subscribe (void)
	Subscribe/Unsubscribe to Monitoring Elements. More...

void	subscribeNew (void)

	testChannel (const edm::ParameterSet &ps)
	Constructor. More...

void	unsubscribe (void)

	~testChannel () override
	Destructor. More...

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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
int	getHeaderSMId (const int headerId)

Private Attributes
int	m_bestPed

int	m_DACmax

int	m_DACmin
	name of the xml file to be saved More...

std::string	m_digiCollection

std::string	m_digiProducer
	secondary name given to collection of digis More...

std::string	m_headerProducer
	name of module/plugin/producer making digis More...

TH2F	m_pedVSDAC

double	m_RMSmax

TH2F	m_singlePedVSDAC_1

TH2F	m_singlePedVSDAC_2

TH2F	m_singlePedVSDAC_3

std::string	m_xmlFile
	name of module/plugin/producer making headers More...

int	m_xtal

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

void	consumesMany (const TypeToGet &id)

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

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

calculate the best DAC value to obtain a pedestal = 200

Author: P. Govoni (testC.nosp@m.hann.nosp@m.el.go.nosp@m.voni.nosp@m.@cern.nosp@m.NOSP.nosp@m.AM.ch)

Definition at line 36 of file testChannel.h.

Constructor & Destructor Documentation

testChannel::testChannel ( const edm::ParameterSet & ps )

Constructor.

ctor

Definition at line 17 of file testChannel.cc.

References m_bestPed, m_DACmax, m_DACmin, and m_RMSmax.

     : m_digiCollection(paramSet.getParameter<std::string>("digiCollection")),
       m_digiProducer(paramSet.getParameter<std::string>("digiProducer")),
       m_headerProducer(paramSet.getParameter<std::string>("headerProducer")),
       m_xmlFile(paramSet.getParameter<std::string>("xmlFile")),
       m_DACmin(paramSet.getParameter<int>("DACmin")),
       m_DACmax(paramSet.getParameter<int>("DACmax")),
       m_RMSmax(paramSet.getParameter<double>("RMSmax")),
       m_bestPed(paramSet.getParameter<int>("bestPed")),
       m_xtal(paramSet.getParameter<int>("xtal")),
       m_pedVSDAC("pedVSDAC", "pedVSDAC", 100, 150, 250, m_DACmax - m_DACmin, m_DACmin, m_DACmax),
       m_singlePedVSDAC_1("singlePedVSDAC_1",
                          "pedVSDAC (g1) for xtal " + TString(m_xtal),
                          100,
                          150,
                          250,
                          m_DACmax - m_DACmin,
                          m_DACmin,
                          m_DACmax),
       m_singlePedVSDAC_2("singlePedVSDAC_2",
                          "pedVSDAC (g2) for xtal " + TString(m_xtal),
                          100,
                          150,
                          250,
                          m_DACmax - m_DACmin,
                          m_DACmin,
                          m_DACmax),
       m_singlePedVSDAC_3("singlePedVSDAC_3",
                          "pedVSDAC (g3) for xtal " + TString(m_xtal),
                          100,
                          150,
                          250,
                          m_DACmax - m_DACmin,
                          m_DACmin,
                          m_DACmax) {
   edm::LogInfo("testChannel") << " reading "
                               << " m_DACmin: " << m_DACmin << " m_DACmax: " << m_DACmax << " m_RMSmax: " << m_RMSmax
                               << " m_bestPed: " << m_bestPed;
 }

testChannel::~testChannel ( )

override

Destructor.

dtor

Definition at line 58 of file testChannel.cc.

58 {}

Member Function Documentation

void testChannel::analyze	(	edm::Event const &	event,
		edm::EventSetup const &	eventSetup
	)

overridevirtual

! Analyze

perform te analysis

loop over the headers

Implements edm::one::EDAnalyzerBase.

Definition at line 64 of file testChannel.cc.

References EcalMGPASample::adc(), isotrackApplyRegressor::df, ecalLiteDTU::gainId(), EcalMGPASample::gainId(), getHeaderSMId(), ecalpyutils::ism(), edm::HandleBase::isValid(), LogDebug, m_digiCollection, m_digiProducer, m_headerProducer, m_pedVSDAC, m_singlePedVSDAC_1, m_singlePedVSDAC_2, m_singlePedVSDAC_3, m_xtal, EcalDataFrame::MAXSAMPLES, EcalDCCHeaderBlock::EcalDCCEventSettings::ped_offset, and EcalDataFrame::sample().

                                                                                 {
   LogDebug("testChannel") << "entering analyze ...";
 
   // get the headers
   // (one header for each supermodule)
   edm::Handle<EcalRawDataCollection> DCCHeaders;
   event.getByLabel(m_headerProducer, DCCHeaders);
   if (!DCCHeaders.isValid()) {
     edm::LogError("testChannel") << "Error! can't get the product " << m_headerProducer.c_str();
   }
 
   std::map<int, int> DACvalues;
 
   // loop over the headers
   for (EcalRawDataCollection::const_iterator headerItr = DCCHeaders->begin(); headerItr != DCCHeaders->end();
        ++headerItr) {
     EcalDCCHeaderBlock::EcalDCCEventSettings settings = headerItr->getEventSettings();
     DACvalues[getHeaderSMId(headerItr->id())] = settings.ped_offset;
     //       std::cout << "DCCid: " << headerItr->id () << "" ;
     //       std::cout << "Ped offset DAC: " << settings.ped_offset << "" ;
   }  
 
   // get the digis
   // (one digi for each crystal)
   edm::Handle<EBDigiCollection> pDigis;
   event.getByLabel(m_digiProducer, pDigis);
   if (!pDigis.isValid()) {
     edm::LogError("testChannel") << "Error! can't get the product " << m_digiCollection.c_str();
   }
 
   // loop over the digis
   for (EBDigiCollection::const_iterator itdigi = pDigis->begin(); itdigi != pDigis->end(); ++itdigi) {
     EBDataFrame df(*itdigi);
     int gainId = df.sample(0).gainId();
     int crystalId = EBDetId(itdigi->id()).ic();
     int smId = EBDetId(itdigi->id()).ism();
 
     edm::LogInfo("testChannel") << "channel " << event.id() << "\tcry: " << crystalId << "\tG: " << gainId
                                 << "\tDAC: " << DACvalues[smId];
 
     // loop over the samples
     for (int iSample = 0; iSample < EBDataFrame::MAXSAMPLES; ++iSample) {
       edm::LogInfo("testChannel") << "\t`-->" << df.sample(iSample).adc();
       m_pedVSDAC.Fill(df.sample(iSample).adc(), DACvalues[smId]);
       if (crystalId == m_xtal) {
         if (gainId == 1)
           m_singlePedVSDAC_1.Fill(df.sample(iSample).adc(), DACvalues[smId]);
         if (gainId == 2)
           m_singlePedVSDAC_2.Fill(df.sample(iSample).adc(), DACvalues[smId]);
         if (gainId == 3)
           m_singlePedVSDAC_3.Fill(df.sample(iSample).adc(), DACvalues[smId]);
       }
     }  // loop over the samples
   }    // loop over the digis
 }

void testChannel::beginJob ( void )

overridevirtual

BeginJob.

begin the job

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 61 of file testChannel.cc.

References LogDebug.

61 { LogDebug("testChannel") << "entering beginJob ..."; }

LogDebug

#define LogDebug(id)

Definition: MessageLogger.h:233

void testChannel::endJob ( void )

overridevirtual

EndJob.

perform the minimiation and write results

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 121 of file testChannel.cc.

References m_pedVSDAC, m_singlePedVSDAC_1, m_singlePedVSDAC_2, m_singlePedVSDAC_3, m_xtal, and submitPVResolutionJobs::out.

Referenced by o2olib.O2ORunMgr::executeJob().

                          {
   char ccout[80];
   sprintf(ccout, "out_%d.root", m_xtal);
   TFile out(ccout, "RECREATE");
   out.cd();
   m_pedVSDAC.Write();
   m_singlePedVSDAC_1.Write();
   m_singlePedVSDAC_2.Write();
   m_singlePedVSDAC_3.Write();
   TProfile *profilo1 = m_singlePedVSDAC_1.ProfileX();
   TProfile *profilo2 = m_singlePedVSDAC_2.ProfileX();
   TProfile *profilo3 = m_singlePedVSDAC_3.ProfileX();
   profilo1->Write("singleProfile_1");
   profilo2->Write("singleProfile_2");
   profilo3->Write("singleProfile_3");
   out.Close();
 }

int testChannel::getHeaderSMId ( const int headerId )

private

Definition at line 145 of file testChannel.cc.

Referenced by analyze().

                                                  {
   // PG FIXME temporary solution
   // PG FIXME check it is consistent with the TB!
   return 1;
 }

void testChannel::subscribe ( void )

Subscribe/Unsubscribe to Monitoring Elements.

Definition at line 151 of file testChannel.cc.

151 {}

void testChannel::subscribeNew ( void )

Definition at line 153 of file testChannel.cc.

153 {}

void testChannel::unsubscribe ( void )

Definition at line 155 of file testChannel.cc.

155 {}

Member Data Documentation

int testChannel::m_bestPed

private

Definition at line 70 of file testChannel.h.

Referenced by testChannel().

int testChannel::m_DACmax

private

Definition at line 68 of file testChannel.h.

Referenced by testChannel().

int testChannel::m_DACmin

private

name of the xml file to be saved

Definition at line 67 of file testChannel.h.

Referenced by testChannel().

std::string testChannel::m_digiCollection

private

Definition at line 61 of file testChannel.h.

Referenced by analyze().

std::string testChannel::m_digiProducer

private

secondary name given to collection of digis

Definition at line 62 of file testChannel.h.

Referenced by analyze().

std::string testChannel::m_headerProducer

private

name of module/plugin/producer making digis

Definition at line 63 of file testChannel.h.

Referenced by analyze().

TH2F testChannel::m_pedVSDAC

private

Definition at line 74 of file testChannel.h.

Referenced by analyze(), and endJob().

double testChannel::m_RMSmax

private

Definition at line 69 of file testChannel.h.

Referenced by testChannel().

TH2F testChannel::m_singlePedVSDAC_1

private

Definition at line 75 of file testChannel.h.

Referenced by analyze(), and endJob().

TH2F testChannel::m_singlePedVSDAC_2

private

Definition at line 76 of file testChannel.h.

Referenced by analyze(), and endJob().

TH2F testChannel::m_singlePedVSDAC_3

private

Definition at line 77 of file testChannel.h.

Referenced by analyze(), and endJob().

std::string testChannel::m_xmlFile

private

name of module/plugin/producer making headers

Definition at line 65 of file testChannel.h.

int testChannel::m_xtal

private

Definition at line 72 of file testChannel.h.

Referenced by analyze(), and endJob().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation