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)
	! Analyze More...

void	beginJob ()
	BeginJob. More...

void	endJob (void)
	EndJob. More...

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

void	subscribeNew (void)

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

void	unsubscribe (void)

virtual	~testChannel ()
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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::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 &)

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

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 "+m_xtal,100,150,250,m_DACmax-m_DACmin,m_DACmin,m_DACmax) ,
   m_singlePedVSDAC_2 ("singlePedVSDAC_2","pedVSDAC (g2) for xtal "+m_xtal,100,150,250,m_DACmax-m_DACmin,m_DACmin,m_DACmax) ,
   m_singlePedVSDAC_3 ("singlePedVSDAC_3","pedVSDAC (g3) for xtal "+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 ( )

virtual

Destructor.

dtor

Definition at line 41 of file testChannel.cc.

42 {

43 }

Member Function Documentation

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

virtual

! Analyze

perform te analysis

loop over the headers

Implements edm::EDAnalyzer.

Definition at line 54 of file testChannel.cc.

References EcalMGPASample::adc(), ecalMGPA::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 )

virtual

BeginJob.

begin the job

Reimplemented from edm::EDAnalyzer.

Definition at line 47 of file testChannel.cc.

References LogDebug.

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

void testChannel::endJob ( void )

virtual

EndJob.

perform the minimiation and write results

Reimplemented from edm::EDAnalyzer.

Definition at line 124 of file testChannel.cc.

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

 {
   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 149 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 159 of file testChannel.cc.

160 {}

void testChannel::subscribeNew ( void )

Definition at line 162 of file testChannel.cc.

163 {}

void testChannel::unsubscribe ( void )

Definition at line 165 of file testChannel.cc.

166 {}

Member Data Documentation

int testChannel::m_bestPed

private

Definition at line 74 of file testChannel.h.

Referenced by testChannel().

int testChannel::m_DACmax

private

Definition at line 72 of file testChannel.h.

Referenced by testChannel().

int testChannel::m_DACmin

private

name of the xml file to be saved

Definition at line 71 of file testChannel.h.

Referenced by testChannel().

std::string testChannel::m_digiCollection

private

Definition at line 65 of file testChannel.h.

Referenced by analyze().

std::string testChannel::m_digiProducer

private

secondary name given to collection of digis

Definition at line 66 of file testChannel.h.

Referenced by analyze().

std::string testChannel::m_headerProducer

private

name of module/plugin/producer making digis

Definition at line 67 of file testChannel.h.

Referenced by analyze().

TH2F testChannel::m_pedVSDAC

private

Definition at line 78 of file testChannel.h.

Referenced by analyze(), and endJob().

double testChannel::m_RMSmax

private

Definition at line 73 of file testChannel.h.

Referenced by testChannel().

TH2F testChannel::m_singlePedVSDAC_1

private

Definition at line 79 of file testChannel.h.

Referenced by analyze(), and endJob().

TH2F testChannel::m_singlePedVSDAC_2

private

Definition at line 80 of file testChannel.h.

Referenced by analyze(), and endJob().

TH2F testChannel::m_singlePedVSDAC_3

private

Definition at line 81 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 69 of file testChannel.h.

int testChannel::m_xtal

private

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