#include <EcalPulseShapeGrapher.h>

Inheritance diagram for EcalPulseShapeGrapher:

Public Member Functions
	EcalPulseShapeGrapher (const edm::ParameterSet &)

	~EcalPulseShapeGrapher () override=default

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
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	endJob () override

std::string	intToString (int)

Private Attributes
int	abscissa [10]

int	ampCut_

std::map< int, TH1F * >	ampHistMap_

std::map< int, TH1F * >	cutAmpHistMap_

const edm::EDGetTokenT< EBDigiCollection >	EBDigis_

const edm::EDGetTokenT< EcalUncalibratedRecHitCollection >	EBUncalibratedRecHitCollection_

const edm::EDGetTokenT< EEDigiCollection >	EEDigis_

const edm::EDGetTokenT< EcalUncalibratedRecHitCollection >	EEUncalibratedRecHitCollection_

EcalFedMap *	fedMap_

TFile *	file_

std::map< int, TH1F * >	firstSampleHistMap_

std::vector< int >	listChannels_

int	ordinate [10]

std::map< int, TH2F * >	pulseShapeHistMap_

std::map< int, TH2F * >	rawPulseShapeHistMap_

std::string	rootFilename_

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

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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 46 of file EcalPulseShapeGrapher.h.

Constructor & Destructor Documentation

◆ EcalPulseShapeGrapher()

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

explicit

Definition at line 32 of file EcalPulseShapeGrapher.cc.

References abscissa, ampHistMap_, cutAmpHistMap_, fedMap_, firstSampleHistMap_, edm::ParameterSet::getUntrackedParameter(), mps_fire::i, intToString(), listChannels_, Skims_PA_cff::name, pulseShapeHistMap_, rawPulseShapeHistMap_, AlCaHLTBitMon_QueryRunRegistry::string, and runGCPTkAlMap::title.

     : EBUncalibratedRecHitCollection_(consumes<EcalUncalibratedRecHitCollection>(
           iConfig.getParameter<edm::InputTag>("EBUncalibratedRecHitCollection"))),
       EBDigis_(consumes<EBDigiCollection>(iConfig.getParameter<edm::InputTag>("EBDigiCollection"))),
       EEUncalibratedRecHitCollection_(consumes<EcalUncalibratedRecHitCollection>(
           iConfig.getParameter<edm::InputTag>("EEUncalibratedRecHitCollection"))),
       EEDigis_(consumes<EEDigiCollection>(iConfig.getParameter<edm::InputTag>("EEDigiCollection"))),
       ampCut_(iConfig.getUntrackedParameter<int>("AmplitudeCutADC", 13)),
       rootFilename_(iConfig.getUntrackedParameter<std::string>("rootFilename", "pulseShapeGrapher")) {
   //now do what ever initialization is needed
 
   std::vector<int> listDefaults;
   listDefaults.push_back(-1);
   listChannels_ = iConfig.getUntrackedParameter<std::vector<int> >("listChannels", listDefaults);
 
   for (std::vector<int>::const_iterator itr = listChannels_.begin(); itr != listChannels_.end(); ++itr) {
     std::string title = "Amplitude of cry " + intToString(*itr);
     std::string name = "ampOfCry" + intToString(*itr);
     ampHistMap_[*itr] = new TH1F(name.c_str(), title.c_str(), 100, 0, 100);
     ampHistMap_[*itr]->GetXaxis()->SetTitle("ADC");
 
     title = "Amplitude (over 13 ADC) of cry " + intToString(*itr);
     name = "cutAmpOfCry" + intToString(*itr);
     cutAmpHistMap_[*itr] = new TH1F(name.c_str(), title.c_str(), 100, 0, 100);
     cutAmpHistMap_[*itr]->GetXaxis()->SetTitle("ADC");
 
     title = "Pulse shape of cry " + intToString(*itr);
     name = "PulseShapeCry" + intToString(*itr);
     pulseShapeHistMap_[*itr] = new TH2F(name.c_str(), title.c_str(), 10, 0, 10, 220, -20, 2);
     pulseShapeHistMap_[*itr]->GetXaxis()->SetTitle("sample");
 
     title = "Raw Pulse shape of cry " + intToString(*itr);
     name = "RawPulseShapeCry" + intToString(*itr);
     rawPulseShapeHistMap_[*itr] = new TH2F(name.c_str(), title.c_str(), 10, 0, 10, 500, 0, 500);
     rawPulseShapeHistMap_[*itr]->GetXaxis()->SetTitle("sample");
 
     title = "Amplitude of first sample, cry " + intToString(*itr);
     name = "AmpOfFirstSampleCry" + intToString(*itr);
     firstSampleHistMap_[*itr] = new TH1F(name.c_str(), title.c_str(), 300, 100, 400);
     firstSampleHistMap_[*itr]->GetXaxis()->SetTitle("ADC");
   }
 
   fedMap_ = new EcalFedMap();
 
   for (int i = 0; i < 10; i++)
     abscissa[i] = i;
 }

◆ ~EcalPulseShapeGrapher()

EcalPulseShapeGrapher::~EcalPulseShapeGrapher ( )

overridedefault

Member Function Documentation

◆ analyze()

void EcalPulseShapeGrapher::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 85 of file EcalPulseShapeGrapher.cc.

References ampCut_, ampHistMap_, CustomPhysics_cfi::amplitude, edm::SortedCollection< T, SORT >::begin(), edm::DataFrameContainer::begin(), cutAmpHistMap_, EcalElectronicsId::dccId(), hgcalPerformanceValidation::df, EBDigis_, EBUncalibratedRecHitCollection_, EEDigis_, EEUncalibratedRecHitCollection_, edm::SortedCollection< T, SORT >::end(), edm::DataFrameContainer::end(), fedMap_, firstSampleHistMap_, PedestalClient_cfi::gain, EcalFedMap::getSliceFromFed(), EBDetId::hashedIndex(), EEDetId::hashedIndex(), mps_fire::i, hit::id, iEvent, createfilelist::int, listChannels_, EcalCondDBWriter_cfi::pedestal, pulseShapeHistMap_, rawPulseShapeHistMap_, and findQualityFiles::size.

                                                                                        {
   using namespace edm;
   using namespace std;
 
   int numHitsWithActivity = 0;
 
   const Handle<EcalUncalibratedRecHitCollection>& EBHits = iEvent.getHandle(EBUncalibratedRecHitCollection_);
   const Handle<EcalUncalibratedRecHitCollection>& EEHits = iEvent.getHandle(EEUncalibratedRecHitCollection_);
   const Handle<EBDigiCollection>& EBdigis = iEvent.getHandle(EBDigis_);
   const Handle<EEDigiCollection>& EEdigis = iEvent.getHandle(EEDigis_);
 
   unique_ptr<EcalElectronicsMapping> ecalElectronicsMap(new EcalElectronicsMapping);
 
   // Loop over the hits
   for (EcalUncalibratedRecHitCollection::const_iterator hitItr = EBHits->begin(); hitItr != EBHits->end(); ++hitItr) {
     EcalUncalibratedRecHit hit = (*hitItr);
     float amplitude = hit.amplitude();
     EBDetId hitDetId = hit.id();
 
     // Get the Fedid
     EcalElectronicsId elecId = ecalElectronicsMap->getElectronicsId(hitDetId);
     int FEDid = 600 + elecId.dccId();
     string SMname = fedMap_->getSliceFromFed(FEDid);
 
     vector<int>::const_iterator itr = listChannels_.begin();
     while (itr != listChannels_.end() && (*itr) != hitDetId.hashedIndex()) {
       itr++;
     }
     if (itr == listChannels_.end())
       continue;
 
     ampHistMap_[hitDetId.hashedIndex()]->Fill(amplitude);
     if (amplitude < ampCut_)
       continue;
 
     cutAmpHistMap_[hitDetId.hashedIndex()]->Fill(amplitude);
     numHitsWithActivity++;
     EBDigiCollection::const_iterator digiItr = EBdigis->begin();
     while (digiItr != EBdigis->end() && digiItr->id() != hitItr->id()) {
       digiItr++;
     }
     if (digiItr == EBdigis->end())
       continue;
 
     double sampleADC[10];
     EBDataFrame df(*digiItr);
     double pedestal = 200;
 
     if (df.sample(0).gainId() != 1 || df.sample(1).gainId() != 1)
       continue;  //goes to the next digi
     else {
       sampleADC[0] = df.sample(0).adc();
       sampleADC[1] = df.sample(1).adc();
       pedestal = (double)(sampleADC[0] + sampleADC[1]) / (double)2;
     }
 
     for (int i = 0; (unsigned int)i < digiItr->size(); ++i) {
       EBDataFrame df(*digiItr);
       double gain = 12.;
       if (df.sample(i).gainId() == 1)
         gain = 1.;
       else if (df.sample(i).gainId() == 2)
         gain = 2.;
       sampleADC[i] = pedestal + (df.sample(i).adc() - pedestal) * gain;
     }
 
     for (int i = 0; i < 10; ++i) {
       pulseShapeHistMap_[hitDetId.hashedIndex()]->Fill(i, (float)(sampleADC[i] - pedestal) / amplitude);
       rawPulseShapeHistMap_[hitDetId.hashedIndex()]->Fill(i, (float)(sampleADC[i]));
     }
     firstSampleHistMap_[hitDetId.hashedIndex()]->Fill(sampleADC[0]);
   }
 
   // Now do the same for the EE hits
   for (EcalUncalibratedRecHitCollection::const_iterator hitItr = EEHits->begin(); hitItr != EEHits->end(); ++hitItr) {
     EcalUncalibratedRecHit hit = (*hitItr);
     float amplitude = hit.amplitude();
     EEDetId hitDetId = hit.id();
 
     // Get the Fedid
     EcalElectronicsId elecId = ecalElectronicsMap->getElectronicsId(hitDetId);
     int FEDid = 600 + elecId.dccId();
     string SMname = fedMap_->getSliceFromFed(FEDid);
 
     vector<int>::const_iterator itr = listChannels_.begin();
     while (itr != listChannels_.end() && (*itr) != hitDetId.hashedIndex()) {
       itr++;
     }
     if (itr == listChannels_.end())
       continue;
 
     ampHistMap_[hitDetId.hashedIndex()]->Fill(amplitude);
     if (amplitude < ampCut_)
       continue;
 
     cutAmpHistMap_[hitDetId.hashedIndex()]->Fill(amplitude);
     numHitsWithActivity++;
     EEDigiCollection::const_iterator digiItr = EEdigis->begin();
     while (digiItr != EEdigis->end() && digiItr->id() != hitItr->id()) {
       digiItr++;
     }
     if (digiItr == EEdigis->end())
       continue;
 
     double sampleADC[10];
     EEDataFrame df(*digiItr);
     double pedestal = 200;
 
     if (df.sample(0).gainId() != 1 || df.sample(1).gainId() != 1)
       continue;  //goes to the next digi
     else {
       sampleADC[0] = df.sample(0).adc();
       sampleADC[1] = df.sample(1).adc();
       pedestal = (double)(sampleADC[0] + sampleADC[1]) / (double)2;
     }
 
     for (int i = 0; (unsigned int)i < digiItr->size(); ++i) {
       EEDataFrame df(*digiItr);
       double gain = 12.;
       if (df.sample(i).gainId() == 1)
         gain = 1.;
       else if (df.sample(i).gainId() == 2)
         gain = 2.;
       sampleADC[i] = pedestal + (df.sample(i).adc() - pedestal) * gain;
     }
 
     for (int i = 0; i < 10; ++i) {
       pulseShapeHistMap_[hitDetId.hashedIndex()]->Fill(i, (float)(sampleADC[i] - pedestal) / amplitude);
       rawPulseShapeHistMap_[hitDetId.hashedIndex()]->Fill(i, (float)(sampleADC[i]));
     }
     firstSampleHistMap_[hitDetId.hashedIndex()]->Fill(sampleADC[0]);
   }
 }

◆ endJob()

void EcalPulseShapeGrapher::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 220 of file EcalPulseShapeGrapher.cc.

References ampHistMap_, cutAmpHistMap_, file_, firstSampleHistMap_, listChannels_, rawPulseShapeHistMap_, rootFilename_, and RandomServiceHelper::t2.

Referenced by o2olib.O2ORunMgr::executeJob().

                                    {
   rootFilename_ += ".root";
   file_ = new TFile(rootFilename_.c_str(), "RECREATE");
   TH1::AddDirectory(false);
 
   for (std::vector<int>::const_iterator itr = listChannels_.begin(); itr != listChannels_.end(); ++itr) {
     ampHistMap_[*itr]->Write();
     cutAmpHistMap_[*itr]->Write();
     firstSampleHistMap_[*itr]->Write();
 
     rawPulseShapeHistMap_[*itr]->Write();
     TProfile* t2 = (TProfile*)(rawPulseShapeHistMap_[*itr]->ProfileX());
     t2->Write();
     //TODO: fix the normalization so these are correct
   }
 
   file_->Write();
   file_->Close();
 }

◆ intToString()

std::string EcalPulseShapeGrapher::intToString ( int num )

private

Definition at line 240 of file EcalPulseShapeGrapher.cc.

References EgammaValidation_cff::num.

Referenced by EcalPulseShapeGrapher().

                                                     {
   using namespace std;
   ostringstream myStream;
   myStream << num << flush;
   return (myStream.str());  //returns the string form of the stringstream object
 }

Member Data Documentation

◆ abscissa

int EcalPulseShapeGrapher::abscissa[10]

private

Definition at line 62 of file EcalPulseShapeGrapher.h.

Referenced by EcalPulseShapeGrapher().

◆ ampCut_

int EcalPulseShapeGrapher::ampCut_

private

Definition at line 71 of file EcalPulseShapeGrapher.h.

Referenced by analyze().

◆ ampHistMap_

std::map<int, TH1F*> EcalPulseShapeGrapher::ampHistMap_

private

Definition at line 65 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), EcalPulseShapeGrapher(), and endJob().

◆ cutAmpHistMap_

std::map<int, TH1F*> EcalPulseShapeGrapher::cutAmpHistMap_

private

Definition at line 69 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), EcalPulseShapeGrapher(), and endJob().

◆ EBDigis_

const edm::EDGetTokenT<EBDigiCollection> EcalPulseShapeGrapher::EBDigis_

private

Definition at line 58 of file EcalPulseShapeGrapher.h.

Referenced by analyze().

◆ EBUncalibratedRecHitCollection_

const edm::EDGetTokenT<EcalUncalibratedRecHitCollection> EcalPulseShapeGrapher::EBUncalibratedRecHitCollection_

private

Definition at line 57 of file EcalPulseShapeGrapher.h.

Referenced by analyze().

◆ EEDigis_

const edm::EDGetTokenT<EEDigiCollection> EcalPulseShapeGrapher::EEDigis_

private

Definition at line 60 of file EcalPulseShapeGrapher.h.

Referenced by analyze().

◆ EEUncalibratedRecHitCollection_

const edm::EDGetTokenT<EcalUncalibratedRecHitCollection> EcalPulseShapeGrapher::EEUncalibratedRecHitCollection_

private

Definition at line 59 of file EcalPulseShapeGrapher.h.

Referenced by analyze().

◆ fedMap_

EcalFedMap* EcalPulseShapeGrapher::fedMap_

private

Definition at line 76 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), and EcalPulseShapeGrapher().

◆ file_

TFile* EcalPulseShapeGrapher::file_

private

Definition at line 74 of file EcalPulseShapeGrapher.h.

Referenced by endJob().

◆ firstSampleHistMap_

std::map<int, TH1F*> EcalPulseShapeGrapher::firstSampleHistMap_

private

Definition at line 67 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), EcalPulseShapeGrapher(), and endJob().

◆ listChannels_

std::vector<int> EcalPulseShapeGrapher::listChannels_

private

Definition at line 64 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), EcalPulseShapeGrapher(), and endJob().

◆ ordinate

int EcalPulseShapeGrapher::ordinate[10]

private

Definition at line 63 of file EcalPulseShapeGrapher.h.

◆ pulseShapeHistMap_

std::map<int, TH2F*> EcalPulseShapeGrapher::pulseShapeHistMap_

private

Definition at line 66 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), and EcalPulseShapeGrapher().

◆ rawPulseShapeHistMap_

std::map<int, TH2F*> EcalPulseShapeGrapher::rawPulseShapeHistMap_

private

Definition at line 68 of file EcalPulseShapeGrapher.h.

Referenced by analyze(), EcalPulseShapeGrapher(), and endJob().

◆ rootFilename_

std::string EcalPulseShapeGrapher::rootFilename_

private

Definition at line 72 of file EcalPulseShapeGrapher.h.

Referenced by endJob().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalPulseShapeGrapher()

◆ ~EcalPulseShapeGrapher()

Member Function Documentation

◆ analyze()

◆ endJob()

◆ intToString()

Member Data Documentation

◆ abscissa

◆ ampCut_

◆ ampHistMap_

◆ cutAmpHistMap_

◆ EBDigis_

◆ EBUncalibratedRecHitCollection_

◆ EEDigis_

◆ EEUncalibratedRecHitCollection_

◆ fedMap_

◆ file_

◆ firstSampleHistMap_

◆ listChannels_

◆ ordinate

◆ pulseShapeHistMap_

◆ rawPulseShapeHistMap_

◆ rootFilename_