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

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 Member Functions
int	getHeaderSMId (const int headerId)

Private Attributes
const int	m_bestPed

const int	m_DACmax

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

const edm::EDGetTokenT< EBDigiCollection >	m_digiProducerToken

const edm::EDGetTokenT< EcalRawDataCollection >	m_headerProducerToken
	Token to access digis. More...

TH2F	m_pedVSDAC

const double	m_RMSmax

TH2F	m_singlePedVSDAC_1

TH2F	m_singlePedVSDAC_2

TH2F	m_singlePedVSDAC_3

const std::string	m_xmlFile
	Token to access headers. More...

const 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 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

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::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_digiProducerToken(
           consumes<EBDigiCollection>(edm::InputTag(paramSet.getParameter<std::string>("digiProducer")))),
       m_headerProducerToken(
           consumes<EcalRawDataCollection>(edm::InputTag(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::~testChannel ( )

override

Destructor.

dtor

Definition at line 59 of file testChannel.cc.

59 {}

Member Function Documentation

◆ analyze()

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 65 of file testChannel.cc.

References edm::SortedCollection< T, SORT >::begin(), edm::DataFrameContainer::begin(), hgcalPerformanceValidation::df, edm::SortedCollection< T, SORT >::end(), edm::DataFrameContainer::end(), ecalLiteDTU::gainId(), getHeaderSMId(), ecalpyutils::ism(), edm::HandleBase::isValid(), LogDebug, m_digiProducerToken, m_headerProducerToken, m_pedVSDAC, m_singlePedVSDAC_1, m_singlePedVSDAC_2, m_singlePedVSDAC_3, m_xtal, EcalDataFrame::MAXSAMPLES, and EcalDCCHeaderBlock::EcalDCCEventSettings::ped_offset.

                                                                                 {
   LogDebug("testChannel") << "entering analyze ...";
 
   // get the headers
   // (one header for each supermodule)
   const edm::Handle<EcalRawDataCollection> &DCCHeaders = event.getHandle(m_headerProducerToken);
   if (!DCCHeaders.isValid()) {
     edm::LogError("testChannel") << "Error! can't get the product for EcalRawDataCollection";
   }
 
   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)
   const edm::Handle<EBDigiCollection> &pDigis = event.getHandle(m_digiProducerToken);
   if (!pDigis.isValid()) {
     edm::LogError("testChannel") << "Error! can't get the product for EBDigiCollection";
   }
 
   // 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
 }

◆ beginJob()

void testChannel::beginJob ( )

overridevirtual

BeginJob.

begin the job

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 62 of file testChannel.cc.

References LogDebug.

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

LogDebug

#define LogDebug(id)

Definition: MessageLogger.h:241

◆ endJob()

void testChannel::endJob ( void )

overridevirtual

EndJob.

perform the minimiation and write results

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 120 of file testChannel.cc.

References m_pedVSDAC, m_singlePedVSDAC_1, m_singlePedVSDAC_2, m_singlePedVSDAC_3, m_xtal, and MillePedeFileConverter_cfg::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();
 }

◆ getHeaderSMId()

int testChannel::getHeaderSMId ( const int headerId )

private

Definition at line 144 of file testChannel.cc.

Referenced by analyze().

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

◆ subscribe()

void testChannel::subscribe ( void )

Subscribe/Unsubscribe to Monitoring Elements.

Definition at line 150 of file testChannel.cc.

150 {}

◆ subscribeNew()

void testChannel::subscribeNew ( void )

Definition at line 152 of file testChannel.cc.

152 {}

◆ unsubscribe()

void testChannel::unsubscribe ( void )

Definition at line 154 of file testChannel.cc.

154 {}

Member Data Documentation

◆ m_bestPed

const int testChannel::m_bestPed

private

Definition at line 69 of file testChannel.h.

Referenced by testChannel().

◆ m_DACmax

const int testChannel::m_DACmax

private

Definition at line 67 of file testChannel.h.

Referenced by testChannel().

◆ m_DACmin

const int testChannel::m_DACmin

private

name of the xml file to be saved

Definition at line 66 of file testChannel.h.

Referenced by testChannel().

◆ m_digiProducerToken

const edm::EDGetTokenT<EBDigiCollection> testChannel::m_digiProducerToken

private

Definition at line 61 of file testChannel.h.

Referenced by analyze().

◆ m_headerProducerToken

const edm::EDGetTokenT<EcalRawDataCollection> testChannel::m_headerProducerToken

private

Token to access digis.

Definition at line 62 of file testChannel.h.

Referenced by analyze().

◆ m_pedVSDAC

TH2F testChannel::m_pedVSDAC

private

Definition at line 73 of file testChannel.h.

Referenced by analyze(), and endJob().

◆ m_RMSmax

const double testChannel::m_RMSmax

private

Definition at line 68 of file testChannel.h.

Referenced by testChannel().

◆ m_singlePedVSDAC_1

TH2F testChannel::m_singlePedVSDAC_1

private

Definition at line 74 of file testChannel.h.

Referenced by analyze(), and endJob().

◆ m_singlePedVSDAC_2

TH2F testChannel::m_singlePedVSDAC_2

private

Definition at line 75 of file testChannel.h.

Referenced by analyze(), and endJob().

◆ m_singlePedVSDAC_3

TH2F testChannel::m_singlePedVSDAC_3

private

Definition at line 76 of file testChannel.h.

Referenced by analyze(), and endJob().

◆ m_xmlFile

const std::string testChannel::m_xmlFile

private

Token to access headers.

Definition at line 64 of file testChannel.h.

◆ m_xtal

const int testChannel::m_xtal

private

Definition at line 71 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

◆ testChannel()

◆ ~testChannel()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ getHeaderSMId()

◆ subscribe()

◆ subscribeNew()

◆ unsubscribe()

Member Data Documentation

◆ m_bestPed

◆ m_DACmax

◆ m_DACmin

◆ m_digiProducerToken

◆ m_headerProducerToken

◆ m_pedVSDAC

◆ m_RMSmax

◆ m_singlePedVSDAC_1

◆ m_singlePedVSDAC_2

◆ m_singlePedVSDAC_3

◆ m_xmlFile

◆ m_xtal