testChannel Class Reference

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

{
}

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.

{}

void testChannel::subscribeNew

(

void

)

Definition at line 162 of file testChannel.cc.

{}

void testChannel::unsubscribe

(

void

)

Definition at line 165 of file testChannel.cc.

{}

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