#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 noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept 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 noexcept

bool	wantsStreamRuns () const noexcept

	~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
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 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;
 
   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);
     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);
     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 216 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 236 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_