#include <EcalMipGraphs.h>

Inheritance diagram for EcalMipGraphs:

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

	~EcalMipGraphs () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	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 (edm::Event const &, edm::EventSetup const &) override

void	beginRun (edm::Run const &, edm::EventSetup const &) override

void	endJob () override

void	endRun (edm::Run const &, edm::EventSetup const &) override

std::string	floatToString (float num)

int	getEEIndex (EcalElectronicsId elecId)

void	initHists (int)

std::string	intToString (int num)

TGraph *	selectDigi (DetId det, int ievt)

void	selectHits (edm::Handle< EcalRecHitCollection > hits, int ievt)

void	writeGraphs ()

Private Attributes
int	abscissa [10]

TH1F *	allFedsTimingHist_

const CaloTopology *	caloTopo_

TTree *	canvasNames_

std::map< int, float >	crysAndAmplitudesMap_

const edm::InputTag	EBDigis_

edm::Handle< EBDigiCollection >	EBdigisHandle

const edm::EDGetTokenT< EBDigiCollection >	ebDigiToken_

const edm::InputTag	EBRecHitCollection_

const edm::EDGetTokenT< EcalRecHitCollection >	ebRecHitToken_

const EcalElectronicsMapping *	ecalElectronicsMap_

const edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	ecalMappingToken_

const edm::InputTag	EEDigis_

edm::Handle< EEDigiCollection >	EEdigisHandle

const edm::EDGetTokenT< EEDigiCollection >	eeDigiToken_

const edm::InputTag	EERecHitCollection_

const edm::EDGetTokenT< EcalRecHitCollection >	eeRecHitToken_

EcalFedMap *	fedMap_

std::map< int, EcalDCCHeaderBlock >	FEDsAndDCCHeaders_

std::map< int, TH1F * >	FEDsAndTimingHists_

TFile *	file_

const edm::InputTag	headerProducer_

std::set< EBDetId >	listEBChannels

std::set< EEDetId >	listEEChannels

std::vector< int >	maskedChannels_

std::vector< std::string >	maskedEBs_

std::vector< int >	maskedFEDs_

const double	minTimingAmp_

int	naiveEvtNum_

std::vector< std::string > *	names

int	ordinate [10]

const edm::EDGetTokenT< EcalRawDataCollection >	rawDataToken_

int	runNum_

std::vector< int >	seedCrys_

std::map< std::string, int >	seedFrequencyMap_

const int	side_

const double	threshold_

const edm::ESGetToken< CaloTopology, CaloTopologyRecord >	topologyToken_

Static Private Attributes
static edm::Service< TFileService >	fileService

static float	gainRatio [3] = {1., 2., 12.}

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 59 of file EcalMipGraphs.h.

Constructor & Destructor Documentation

◆ EcalMipGraphs()

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

explicit

Definition at line 40 of file EcalMipGraphs.cc.

References abscissa, allFedsTimingHist_, fedMap_, fileService, EcalFedMap::getFedFromSlice(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, TFileService::make(), maskedChannels_, maskedEBs_, maskedFEDs_, naiveEvtNum_, and seedCrys_.

     : EBRecHitCollection_(iConfig.getParameter<edm::InputTag>("EcalRecHitCollectionEB")),
       EERecHitCollection_(iConfig.getParameter<edm::InputTag>("EcalRecHitCollectionEE")),
       EBDigis_(iConfig.getParameter<edm::InputTag>("EBDigiCollection")),
       EEDigis_(iConfig.getParameter<edm::InputTag>("EEDigiCollection")),
       headerProducer_(iConfig.getParameter<edm::InputTag>("headerProducer")),
       rawDataToken_(consumes<EcalRawDataCollection>(headerProducer_)),
       ebRecHitToken_(consumes<EcalRecHitCollection>(EBRecHitCollection_)),
       eeRecHitToken_(consumes<EcalRecHitCollection>(EERecHitCollection_)),
       ebDigiToken_(consumes<EBDigiCollection>(EBDigis_)),
       eeDigiToken_(consumes<EEDigiCollection>(EEDigis_)),
       ecalMappingToken_(esConsumes<edm::Transition::BeginRun>()),
       topologyToken_(esConsumes()),
       runNum_(-1),
       side_(iConfig.getUntrackedParameter<int>("side", 3)),
       threshold_(iConfig.getUntrackedParameter<double>("amplitudeThreshold", 12.0)),
       minTimingAmp_(iConfig.getUntrackedParameter<double>("minimumTimingAmplitude", 0.100)) {
   vector<int> listDefaults;
   listDefaults.push_back(-1);
 
   maskedChannels_ = iConfig.getUntrackedParameter<vector<int> >("maskedChannels", listDefaults);
   maskedFEDs_ = iConfig.getUntrackedParameter<vector<int> >("maskedFEDs", listDefaults);
   seedCrys_ = iConfig.getUntrackedParameter<vector<int> >("seedCrys", vector<int>());
 
   vector<string> defaultMaskedEBs;
   defaultMaskedEBs.push_back("none");
   maskedEBs_ = iConfig.getUntrackedParameter<vector<string> >("maskedEBs", defaultMaskedEBs);
 
   fedMap_ = new EcalFedMap();
 
   string title1 = "Jitter for all FEDs";
   string name1 = "JitterAllFEDs";
   allFedsTimingHist_ = fileService->make<TH1F>(name1.c_str(), title1.c_str(), 150, -7, 7);
 
   // load up the maskedFED list with the proper FEDids
   if (maskedFEDs_[0] == -1) {
     //if "actual" EB id given, then convert to FEDid and put in listFEDs_
     if (maskedEBs_[0] != "none") {
       maskedFEDs_.clear();
       for (vector<string>::const_iterator ebItr = maskedEBs_.begin(); ebItr != maskedEBs_.end(); ++ebItr) {
         maskedFEDs_.push_back(fedMap_->getFedFromSlice(*ebItr));
       }
     }
   }
 
   for (int i = 0; i < 10; i++)
     abscissa[i] = i;
 
   naiveEvtNum_ = 0;
 }

◆ ~EcalMipGraphs()

EcalMipGraphs::~EcalMipGraphs ( )

override

Definition at line 91 of file EcalMipGraphs.cc.

91 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 98 of file EcalMipGraphs.cc.

References edm::SortedCollection< T, SORT >::begin(), caloTopo_, canvasNames_, EBdigisHandle, ebDigiToken_, ebRecHitToken_, EEdigisHandle, eeDigiToken_, eeRecHitToken_, edm::SortedCollection< T, SORT >::end(), FEDsAndDCCHeaders_, fileService, edm::EventSetup::getData(), iEvent, listEBChannels, listEEChannels, TFileService::make(), naiveEvtNum_, names, rawDataToken_, runNum_, selectHits(), and topologyToken_.

                                                                                {
   // get the headers
   // (one header for each supermodule)
   edm::Handle<EcalRawDataCollection> DCCHeaders;
   iEvent.getByToken(rawDataToken_, DCCHeaders);
 
   for (EcalRawDataCollection::const_iterator headerItr = DCCHeaders->begin(); headerItr != DCCHeaders->end();
        ++headerItr) {
     FEDsAndDCCHeaders_[headerItr->id() + 600] = *headerItr;
   }
 
   int ievt = iEvent.id().event();
   naiveEvtNum_++;
 
   if (runNum_ == -1) {
     runNum_ = iEvent.id().run();
     canvasNames_ = fileService->make<TTree>("canvasNames", "Names of written canvases");
     names = new std::vector<string>();
     canvasNames_->Branch("names", "vector<string>", &names);
   }
 
   //We only want the 3x3's for this event...
   listEBChannels.clear();
   listEEChannels.clear();
   caloTopo_ = &iSetup.getData(topologyToken_);
 
   Handle<EcalRecHitCollection> EBhits;
   Handle<EcalRecHitCollection> EEhits;
   iEvent.getByToken(ebRecHitToken_, EBhits);
   iEvent.getByToken(eeRecHitToken_, EEhits);
   // Now, retrieve the crystal digi from the event
   iEvent.getByToken(ebDigiToken_, EBdigisHandle);
   iEvent.getByToken(eeDigiToken_, EEdigisHandle);
 
   selectHits(EBhits, ievt);
   selectHits(EEhits, ievt);
 }

◆ beginRun()

void EcalMipGraphs::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	c
	)

overrideprivate

Definition at line 332 of file EcalMipGraphs.cc.

References DummyCfis::c, ecalElectronicsMap_, and ecalMappingToken_.

                                                                   {
   ecalElectronicsMap_ = &c.getData(ecalMappingToken_);
 }

◆ endJob()

void EcalMipGraphs::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 339 of file EcalMipGraphs.cc.

References canvasNames_, ewkTauDQM_cfi::channels, l1tstage2emulator_dqm_sourceclient-live_cfg::feds, intToString(), maskedChannels_, maskedFEDs_, seedFrequencyMap_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by o2olib.O2ORunMgr::executeJob().

                            {
   canvasNames_->Fill();
 
   string frequencies = "";
   for (std::map<std::string, int>::const_iterator itr = seedFrequencyMap_.begin(); itr != seedFrequencyMap_.end();
        ++itr) {
     if (itr->second > 1) {
       frequencies += itr->first;
       frequencies += " Frequency: ";
       frequencies += intToString(itr->second);
       frequencies += "\n";
     }
   }
   LogWarning("EcalMipGraphs") << "Found seeds with frequency > 1: "
                               << "\n\n"
                               << frequencies;
 
   std::string channels;
   for (std::vector<int>::const_iterator itr = maskedChannels_.begin(); itr != maskedChannels_.end(); ++itr) {
     channels += intToString(*itr);
     channels += ",";
   }
 
   std::string feds;
   for (std::vector<int>::const_iterator itr = maskedFEDs_.begin(); itr != maskedFEDs_.end(); ++itr) {
     feds += intToString(*itr);
     feds += ",";
   }
 
   LogWarning("EcalMipGraphs") << "Masked channels are: " << channels;
   LogWarning("EcalMipGraphs") << "Masked FEDs are: " << feds << " and that is all!";
 }

◆ endRun()

void EcalMipGraphs::endRun	(	edm::Run const &	,
		edm::EventSetup const &	c
	)

overrideprivate

Definition at line 336 of file EcalMipGraphs.cc.

336 {}

◆ floatToString()

std::string EcalMipGraphs::floatToString ( float num )

private

Definition at line 379 of file EcalMipGraphs.cc.

References EgammaValidation_cff::num.

Referenced by selectDigi().

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

◆ getEEIndex()

int EcalMipGraphs::getEEIndex ( EcalElectronicsId elecId )

private

Definition at line 314 of file EcalMipGraphs.cc.

References EcalElectronicsId::dccId(), EcalElectronicsId::stripId(), EcalElectronicsId::towerId(), and EcalElectronicsId::xtalId().

Referenced by selectDigi(), and selectHits().

                                                       {
   int FEDid = 600 + elecId.dccId();
   return 10000 * FEDid + 100 * elecId.towerId() + 5 * (elecId.stripId() - 1) + elecId.xtalId();
 }

◆ initHists()

void EcalMipGraphs::initHists ( int FED )

private

Definition at line 320 of file EcalMipGraphs.cc.

References fedMap_, FEDsAndTimingHists_, fileService, EcalFedMap::getSliceFromFed(), intToString(), and TFileService::make().

Referenced by selectHits().

                                      {
   using namespace std;
 
   string title1 = "Jitter for ";
   title1.append(fedMap_->getSliceFromFed(FED));
   string name1 = "JitterFED";
   name1.append(intToString(FED));
   TH1F* timingHist = fileService->make<TH1F>(name1.c_str(), title1.c_str(), 150, -7, 7);
   FEDsAndTimingHists_[FED] = timingHist;
 }

◆ intToString()

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

private

Definition at line 372 of file EcalMipGraphs.cc.

References EgammaValidation_cff::num.

Referenced by endJob(), initHists(), selectDigi(), and selectHits().

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

◆ selectDigi()

TGraph * EcalMipGraphs::selectDigi	(	DetId	det,
		int	ievt
	)

private

Definition at line 136 of file EcalMipGraphs.cc.

References abscissa, edm::DataFrameContainer::begin(), crysAndAmplitudesMap_, EcalElectronicsId::dccId(), DetId::det(), hgcalPerformanceValidation::df, EBdigisHandle, DetId::Ecal, ecalElectronicsMap_, EEdigisHandle, edm::DataFrameContainer::end(), hcalRecHitTable_cff::energy, fedMap_, FEDsAndDCCHeaders_, fileService, nano_mu_digi_cff::float, floatToString(), ecalLiteDTU::gainId(), gainRatio, getEEIndex(), EcalElectronicsMapping::getElectronicsId(), EcalFedMap::getSliceFromFed(), mps_fire::i, intToString(), PixelPluginsPhase0_cfi::isBarrel, TFileService::make(), naiveEvtNum_, Skims_PA_cff::name, ordinate, EcalCondDBWriter_cfi::pedestal, and runGCPTkAlMap::title.

Referenced by selectHits().

                                                          {
   int emptyY[10];
   for (int i = 0; i < 10; i++)
     emptyY[i] = 0;
   TGraph* emptyGraph = fileService->make<TGraph>(10, abscissa, emptyY);
   emptyGraph->SetTitle("NOT ECAL");
 
   //If the DetId is not from Ecal, return
   if (thisDet.det() != DetId::Ecal)
     return emptyGraph;
 
   emptyGraph->SetTitle("NO DIGIS");
   //find digi we need  -- can't get find() to work; need DataFrame(DetId det) to work?
   EcalElectronicsId elecId = ecalElectronicsMap_->getElectronicsId(thisDet);
   int FEDid = 600 + elecId.dccId();
   bool isBarrel = true;
   if (FEDid < 610 || FEDid > 645)
     isBarrel = false;
   int cryIndex = isBarrel ? ((EBDetId)thisDet).hashedIndex() : getEEIndex(elecId);
   int ic = isBarrel ? ((EBDetId)thisDet).ic() : cryIndex;
 
   string sliceName = fedMap_->getSliceFromFed(FEDid);
   EcalDataFrame df;
   if (isBarrel) {
     EBDigiCollection::const_iterator digiItr = EBdigisHandle->begin();
     while (digiItr != EBdigisHandle->end() && ((*digiItr).id() != (EBDetId)thisDet)) {
       ++digiItr;
     }
     if (digiItr == EBdigisHandle->end()) {
       //LogWarning("EcalMipGraphs") << "Cannot find digi for ic:" << ic
       //  << " FED:" << FEDid << " evt:" << naiveEvtNum_;
       return emptyGraph;
     } else
       df = *digiItr;
   } else {
     EEDigiCollection::const_iterator digiItr = EEdigisHandle->begin();
     while (digiItr != EEdigisHandle->end() && ((*digiItr).id() != (EEDetId)thisDet)) {
       ++digiItr;
     }
     if (digiItr == EEdigisHandle->end()) {
       //LogWarning("EcalMipGraphs") << "Cannot find digi for ic:" << ic
       //  << " FED:" << FEDid << " evt:" << naiveEvtNum_;
       return emptyGraph;
     } else
       df = *digiItr;
   }
 
   int gainId = FEDsAndDCCHeaders_[FEDid].getMgpaGain();
   int gainHuman;
   if (gainId == 1)
     gainHuman = 12;
   else if (gainId == 2)
     gainHuman = 6;
   else if (gainId == 3)
     gainHuman = 1;
   else
     gainHuman = -1;
 
   double pedestal = 200;
 
   emptyGraph->SetTitle("FIRST TWO SAMPLES NOT GAIN12");
   if (df.sample(0).gainId() != 1 || df.sample(1).gainId() != 1)
     return emptyGraph;  //goes to the next digi
   else {
     ordinate[0] = df.sample(0).adc();
     ordinate[1] = df.sample(1).adc();
     pedestal = (double)(ordinate[0] + ordinate[1]) / (double)2;
   }
 
   for (int i = 0; i < df.size(); ++i) {
     if (df.sample(i).gainId() != 0)
       ordinate[i] = (int)(pedestal + (df.sample(i).adc() - pedestal) * gainRatio[df.sample(i).gainId() - 1]);
     else
       ordinate[i] = 49152;  //Saturation of gain1
   }
 
   TGraph* oneGraph = fileService->make<TGraph>(10, abscissa, ordinate);
   string name = "Graph_ev" + intToString(naiveEvtNum_) + "_ic" + intToString(ic) + "_FED" + intToString(FEDid);
   string gainString = (gainId == 1) ? "Free" : intToString(gainHuman);
   string title = "Event" + intToString(naiveEvtNum_) + "_lv1a" + intToString(ievt) + "_ic" + intToString(ic) + "_" +
                  sliceName + "_gain" + gainString;
 
   float energy = 0;
   map<int, float>::const_iterator itr;
   itr = crysAndAmplitudesMap_.find(cryIndex);
   if (itr != crysAndAmplitudesMap_.end()) {
     //edm::LogWarning("EcalMipGraphs")<< "itr->second(ampli)="<< itr->second;
     energy = (float)itr->second;
   }
   //else
   //edm::LogWarning("EcalMipGraphs") << "cry " << ic << "not found in ampMap";
 
   title += "_Energy" + floatToString(round(energy * 1000));
 
   oneGraph->SetTitle(title.c_str());
   oneGraph->SetName(name.c_str());
   oneGraph->GetXaxis()->SetTitle("sample");
   oneGraph->GetYaxis()->SetTitle("ADC");
   return oneGraph;
 }

◆ selectHits()

void EcalMipGraphs::selectHits	(	edm::Handle< EcalRecHitCollection >	hits,
		int	ievt
	)

private

Definition at line 237 of file EcalMipGraphs.cc.

References allFedsTimingHist_, caloTopo_, timingPdfMaker::can, crysAndAmplitudesMap_, EcalElectronicsId::dccId(), ecalElectronicsMap_, FEDsAndTimingHists_, spr::find(), getEEIndex(), EcalElectronicsMapping::getElectronicsId(), CaloTopology::getSubdetectorTopology(), hfClusterShapes_cfi::hits, CaloNavigator< T, TOPO >::home(), mps_fire::i, hit::id, initHists(), intToString(), PixelPluginsPhase0_cfi::isBarrel, dqmiolumiharvest::j, maskedChannels_, maskedFEDs_, minTimingAmp_, naiveEvtNum_, Skims_PA_cff::name, names, CaloNavigator< T, TOPO >::offsetBy(), mps_fire::result, fileCollector::seed, seedCrys_, seedFrequencyMap_, selectDigi(), side_, threshold_, and runGCPTkAlMap::title.

Referenced by analyze().

                                                                           {
   for (EcalRecHitCollection::const_iterator hitItr = hits->begin(); hitItr != hits->end(); ++hitItr) {
     EcalRecHit hit = (*hitItr);
     DetId det = hit.id();
     EcalElectronicsId elecId = ecalElectronicsMap_->getElectronicsId(det);
     int FEDid = 600 + elecId.dccId();
     bool isBarrel = true;
     if (FEDid < 610 || FEDid > 645)
       isBarrel = false;
     int cryIndex = isBarrel ? ((EBDetId)det).hashedIndex() : ((EEDetId)det).hashedIndex();
     int ic = isBarrel ? ((EBDetId)det).ic() : getEEIndex(elecId);
 
     float ampli = hit.energy();
 
     vector<int>::iterator result;
     result = find(maskedFEDs_.begin(), maskedFEDs_.end(), FEDid);
     if (result != maskedFEDs_.end()) {
       //LogWarning("EcalMipGraphs") << "skipping uncalRecHit for FED " << FEDid << " ; amplitude " << ampli;
       continue;
     }
     result = find(maskedChannels_.begin(), maskedChannels_.end(), cryIndex);
     if (result != maskedChannels_.end()) {
       //LogWarning("EcalMipGraphs") << "skipping uncalRecHit for channel: " << cryIndex << " in fed: " << FEDid << " with amplitude " << ampli ;
       continue;
     }
     bool cryIsInList = false;
     result = find(seedCrys_.begin(), seedCrys_.end(), cryIndex);
     if (result != seedCrys_.end())
       cryIsInList = true;
 
     // Either we must have a user-requested cry (in which case there is no amplitude selection)
     // Or we pick all crys that pass the amplitude cut (in which case there is no fixed crystal selection)
     if (cryIsInList || (seedCrys_.empty() && ampli > threshold_)) {
       // We have a winner!
       crysAndAmplitudesMap_[cryIndex] = ampli;
       string name = "Event" + intToString(naiveEvtNum_) + "_ic" + intToString(ic) + "_FED" + intToString(FEDid);
       string title = "Digis";
       string seed = "ic" + intToString(ic) + "_FED" + intToString(FEDid);
       int freq = 1;
       pair<map<string, int>::iterator, bool> pair = seedFrequencyMap_.insert(make_pair(seed, freq));
       if (!pair.second) {
         ++(pair.first->second);
       }
 
       TCanvas can(name.c_str(), title.c_str(), 200, 50, 900, 900);
       can.Divide(side_, side_);
       TGraph* myGraph;
       int canvasNum = 1;
 
       CaloNavigator<DetId> cursor = CaloNavigator<DetId>(det, caloTopo_->getSubdetectorTopology(det));
       //Now put each graph in one by one
       for (int j = side_ / 2; j >= -side_ / 2; --j) {
         for (int i = -side_ / 2; i <= side_ / 2; ++i) {
           cursor.home();
           cursor.offsetBy(i, j);
           can.cd(canvasNum);
           myGraph = selectDigi(*cursor, ievt);
           myGraph->Draw("A*");
           canvasNum++;
         }
       }
       can.Write();
       names->push_back(name);
     }
 
     TH1F* timingHist = FEDsAndTimingHists_[FEDid];
     if (timingHist == nullptr) {
       initHists(FEDid);
       timingHist = FEDsAndTimingHists_[FEDid];
     }
     if (ampli > minTimingAmp_) {
       timingHist->Fill(hit.time());
       allFedsTimingHist_->Fill(hit.time());
     }
   }
 }

◆ writeGraphs()

void EcalMipGraphs::writeGraphs ( )

private

Member Data Documentation

◆ abscissa

int EcalMipGraphs::abscissa[10]

private

Definition at line 105 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), and selectDigi().

◆ allFedsTimingHist_

TH1F* EcalMipGraphs::allFedsTimingHist_

private

Definition at line 121 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), and selectHits().

◆ caloTopo_

const CaloTopology* EcalMipGraphs::caloTopo_

private

Definition at line 127 of file EcalMipGraphs.h.

Referenced by analyze(), and selectHits().

◆ canvasNames_

TTree* EcalMipGraphs::canvasNames_

private

Definition at line 124 of file EcalMipGraphs.h.

Referenced by analyze(), and endJob().

◆ crysAndAmplitudesMap_

std::map<int, float> EcalMipGraphs::crysAndAmplitudesMap_

private

Definition at line 117 of file EcalMipGraphs.h.

Referenced by selectDigi(), and selectHits().

◆ EBDigis_

const edm::InputTag EcalMipGraphs::EBDigis_

private

Definition at line 81 of file EcalMipGraphs.h.

◆ EBdigisHandle

edm::Handle<EBDigiCollection> EcalMipGraphs::EBdigisHandle

private

Definition at line 85 of file EcalMipGraphs.h.

Referenced by analyze(), and selectDigi().

◆ ebDigiToken_

const edm::EDGetTokenT<EBDigiCollection> EcalMipGraphs::ebDigiToken_

private

Definition at line 91 of file EcalMipGraphs.h.

Referenced by analyze().

◆ EBRecHitCollection_

const edm::InputTag EcalMipGraphs::EBRecHitCollection_

private

Definition at line 79 of file EcalMipGraphs.h.

◆ ebRecHitToken_

const edm::EDGetTokenT<EcalRecHitCollection> EcalMipGraphs::ebRecHitToken_

private

Definition at line 89 of file EcalMipGraphs.h.

Referenced by analyze().

◆ ecalElectronicsMap_

const EcalElectronicsMapping* EcalMipGraphs::ecalElectronicsMap_

private

Definition at line 126 of file EcalMipGraphs.h.

Referenced by beginRun(), selectDigi(), and selectHits().

◆ ecalMappingToken_

const edm::ESGetToken<EcalElectronicsMapping, EcalMappingRcd> EcalMipGraphs::ecalMappingToken_

private

Definition at line 94 of file EcalMipGraphs.h.

Referenced by beginRun().

◆ EEDigis_

const edm::InputTag EcalMipGraphs::EEDigis_

private

Definition at line 82 of file EcalMipGraphs.h.

◆ EEdigisHandle

edm::Handle<EEDigiCollection> EcalMipGraphs::EEdigisHandle

private

Definition at line 86 of file EcalMipGraphs.h.

Referenced by analyze(), and selectDigi().

◆ eeDigiToken_

const edm::EDGetTokenT<EEDigiCollection> EcalMipGraphs::eeDigiToken_

private

Definition at line 92 of file EcalMipGraphs.h.

Referenced by analyze().

◆ EERecHitCollection_

const edm::InputTag EcalMipGraphs::EERecHitCollection_

private

Definition at line 80 of file EcalMipGraphs.h.

◆ eeRecHitToken_

const edm::EDGetTokenT<EcalRecHitCollection> EcalMipGraphs::eeRecHitToken_

private

Definition at line 90 of file EcalMipGraphs.h.

Referenced by analyze().

◆ fedMap_

EcalFedMap* EcalMipGraphs::fedMap_

private

Definition at line 125 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), initHists(), and selectDigi().

◆ FEDsAndDCCHeaders_

std::map<int, EcalDCCHeaderBlock> EcalMipGraphs::FEDsAndDCCHeaders_

private

Definition at line 118 of file EcalMipGraphs.h.

Referenced by analyze(), and selectDigi().

◆ FEDsAndTimingHists_

std::map<int, TH1F*> EcalMipGraphs::FEDsAndTimingHists_

private

Definition at line 116 of file EcalMipGraphs.h.

Referenced by initHists(), and selectHits().

◆ file_

TFile* EcalMipGraphs::file_

private

Definition at line 123 of file EcalMipGraphs.h.

◆ fileService

edm::Service< TFileService > EcalMipGraphs::fileService

staticprivate

Definition at line 109 of file EcalMipGraphs.h.

Referenced by analyze(), EcalMipGraphs(), initHists(), and selectDigi().

◆ gainRatio

float EcalMipGraphs::gainRatio = {1., 2., 12.}

staticprivate

Definition at line 108 of file EcalMipGraphs.h.

Referenced by selectDigi().

◆ headerProducer_

const edm::InputTag EcalMipGraphs::headerProducer_

private

Definition at line 83 of file EcalMipGraphs.h.

◆ listEBChannels

std::set<EBDetId> EcalMipGraphs::listEBChannels

private

Definition at line 102 of file EcalMipGraphs.h.

Referenced by analyze().

◆ listEEChannels

std::set<EEDetId> EcalMipGraphs::listEEChannels

private

Definition at line 103 of file EcalMipGraphs.h.

Referenced by analyze().

◆ maskedChannels_

std::vector<int> EcalMipGraphs::maskedChannels_

private

Definition at line 112 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), endJob(), and selectHits().

◆ maskedEBs_

std::vector<std::string> EcalMipGraphs::maskedEBs_

private

Definition at line 115 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs().

◆ maskedFEDs_

std::vector<int> EcalMipGraphs::maskedFEDs_

private

Definition at line 113 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), endJob(), and selectHits().

◆ minTimingAmp_

const double EcalMipGraphs::minTimingAmp_

private

Definition at line 100 of file EcalMipGraphs.h.

Referenced by selectHits().

◆ naiveEvtNum_

int EcalMipGraphs::naiveEvtNum_

private

Definition at line 129 of file EcalMipGraphs.h.

Referenced by analyze(), EcalMipGraphs(), selectDigi(), and selectHits().

◆ names

std::vector<std::string>* EcalMipGraphs::names

private

Definition at line 111 of file EcalMipGraphs.h.

Referenced by analyze(), and selectHits().

◆ ordinate

int EcalMipGraphs::ordinate[10]

private

Definition at line 106 of file EcalMipGraphs.h.

Referenced by selectDigi().

◆ rawDataToken_

const edm::EDGetTokenT<EcalRawDataCollection> EcalMipGraphs::rawDataToken_

private

Definition at line 88 of file EcalMipGraphs.h.

Referenced by analyze().

◆ runNum_

int EcalMipGraphs::runNum_

private

Definition at line 97 of file EcalMipGraphs.h.

Referenced by analyze().

◆ seedCrys_

std::vector<int> EcalMipGraphs::seedCrys_

private

Definition at line 114 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), and selectHits().

◆ seedFrequencyMap_

std::map<std::string, int> EcalMipGraphs::seedFrequencyMap_

private

Definition at line 119 of file EcalMipGraphs.h.

Referenced by endJob(), and selectHits().

◆ side_

const int EcalMipGraphs::side_

private

Definition at line 98 of file EcalMipGraphs.h.

Referenced by selectHits().

◆ threshold_

const double EcalMipGraphs::threshold_

private

Definition at line 99 of file EcalMipGraphs.h.

Referenced by selectHits().

◆ topologyToken_

const edm::ESGetToken<CaloTopology, CaloTopologyRecord> EcalMipGraphs::topologyToken_

private

Definition at line 95 of file EcalMipGraphs.h.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalMipGraphs()

◆ ~EcalMipGraphs()

Member Function Documentation

◆ analyze()

◆ beginRun()

◆ endJob()

◆ endRun()

◆ floatToString()

◆ getEEIndex()

◆ initHists()

◆ intToString()

◆ selectDigi()

◆ selectHits()

◆ writeGraphs()

Member Data Documentation

◆ abscissa

◆ allFedsTimingHist_

◆ caloTopo_

◆ canvasNames_

◆ crysAndAmplitudesMap_

◆ EBDigis_

◆ EBdigisHandle

◆ ebDigiToken_

◆ EBRecHitCollection_

◆ ebRecHitToken_

◆ ecalElectronicsMap_

◆ ecalMappingToken_

◆ EEDigis_

◆ EEdigisHandle

◆ eeDigiToken_

◆ EERecHitCollection_

◆ eeRecHitToken_

◆ fedMap_

◆ FEDsAndDCCHeaders_

◆ FEDsAndTimingHists_

◆ file_

◆ fileService

◆ gainRatio

◆ headerProducer_

◆ listEBChannels

◆ listEEChannels

◆ maskedChannels_

◆ maskedEBs_

◆ maskedFEDs_

◆ minTimingAmp_

◆ naiveEvtNum_

◆ names

◆ ordinate

◆ rawDataToken_

◆ runNum_

◆ seedCrys_

◆ seedFrequencyMap_

◆ side_

◆ threshold_

◆ topologyToken_