#include <EcalMipGraphs.h>

Inheritance diagram for EcalMipGraphs:

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

	~EcalMipGraphs () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

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

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::vector< ModuleDescription const * > &modules, 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

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

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

void	endJob () 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, edm::ESHandle< CaloTopology > caloTopo)

void	writeGraphs ()

Private Attributes
int	abscissa [10]

TH1F *	allFedsTimingHist_

TTree *	canvasNames_

std::map< int, float >	crysAndAmplitudesMap_

edm::InputTag	EBDigis_

edm::Handle< EBDigiCollection >	EBdigisHandle

edm::InputTag	EBRecHitCollection_

const EcalElectronicsMapping *	ecalElectronicsMap_

edm::InputTag	EEDigis_

edm::Handle< EEDigiCollection >	EEdigisHandle

edm::InputTag	EERecHitCollection_

EcalFedMap *	fedMap_

std::map< int, EcalDCCHeaderBlock >	FEDsAndDCCHeaders_

std::map< int, TH1F * >	FEDsAndTimingHists_

TFile *	file_

edm::InputTag	headerProducer_

std::set< EBDetId >	listEBChannels

std::set< EEDetId >	listEEChannels

std::vector< int >	maskedChannels_

std::vector< std::string >	maskedEBs_

std::vector< int >	maskedFEDs_

double	minTimingAmp_

int	naiveEvtNum_

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

int	ordinate [10]

int	runNum_

std::vector< int >	seedCrys_

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

int	side_

double	threshold_

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

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

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

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

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

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

Definition at line 59 of file EcalMipGraphs.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 42 of file EcalMipGraphs.cc.

References abscissa, allFedsTimingHist_, fedMap_, fileService, 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")),
       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 ( )

override

Definition at line 86 of file EcalMipGraphs.cc.

86 {}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 93 of file EcalMipGraphs.cc.

References edm::SortedCollection< T, SORT >::begin(), canvasNames_, EBDigis_, EBdigisHandle, EBRecHitCollection_, EEDigis_, EEdigisHandle, EERecHitCollection_, edm::SortedCollection< T, SORT >::end(), edm::EventID::event(), FEDsAndDCCHeaders_, fileService, edm::EventSetup::get(), edm::Event::getByLabel(), headerProducer_, edm::EventBase::id(), listEBChannels, listEEChannels, TFileService::make(), naiveEvtNum_, names, edm::EventID::run(), runNum_, and selectHits().

                                                                                {
   // get the headers
   // (one header for each supermodule)
   edm::Handle<EcalRawDataCollection> DCCHeaders;
   iEvent.getByLabel(headerProducer_, 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();
   Handle<EcalRecHitCollection> EBhits;
   Handle<EcalRecHitCollection> EEhits;
   ESHandle<CaloTopology> caloTopo;
   iSetup.get<CaloTopologyRecord>().get(caloTopo);
   iEvent.getByLabel(EBRecHitCollection_, EBhits);
   iEvent.getByLabel(EERecHitCollection_, EEhits);
   // Now, retrieve the crystal digi from the event
   iEvent.getByLabel(EBDigis_, EBdigisHandle);
   iEvent.getByLabel(EEDigis_, EEdigisHandle);
   //debug
   //LogWarning("EcalMipGraphs") << "event " << ievt << " EBhits collection size " << EBhits->size();
   //LogWarning("EcalMipGraphs") << "event " << ievt << " EEhits collection size " << EEhits->size();
 
   selectHits(EBhits, ievt, caloTopo);
   selectHits(EEhits, ievt, caloTopo);
 }

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

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 330 of file EcalMipGraphs.cc.

References ecalElectronicsMap_, edm::EventSetup::get(), cmsBatch::handle, and edm::ESHandle< T >::product().

                                                                   {
   edm::ESHandle<EcalElectronicsMapping> handle;
   c.get<EcalMappingRcd>().get(handle);
   ecalElectronicsMap_ = handle.product();
 }

void EcalMipGraphs::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 337 of file EcalMipGraphs.cc.

References canvasNames_, 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!";
 }

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

private

Definition at line 377 of file EcalMipGraphs.cc.

Referenced by selectDigi().

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

int EcalMipGraphs::getEEIndex ( EcalElectronicsId elecId )

private

Definition at line 312 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();
 }

void EcalMipGraphs::initHists ( int FED )

private

Definition at line 318 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;
 }

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

private

Definition at line 370 of file EcalMipGraphs.cc.

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
 }

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

private

Definition at line 134 of file EcalMipGraphs.cc.

References abscissa, EcalMGPASample::adc(), edm::DataFrameContainer::begin(), crysAndAmplitudesMap_, EcalElectronicsId::dccId(), DetId::det(), EBdigisHandle, DetId::Ecal, ecalElectronicsMap_, EEdigisHandle, edm::DataFrameContainer::end(), randomXiThetaGunProducer_cfi::energy, fedMap_, FEDsAndDCCHeaders_, fileService, objects.autophobj::float, floatToString(), ecalMGPA::gainId(), EcalMGPASample::gainId(), gainRatio, getEEIndex(), EcalElectronicsMapping::getElectronicsId(), EcalFedMap::getSliceFromFed(), mps_fire::i, intToString(), GeomDetEnumerators::isBarrel(), TFileService::make(), naiveEvtNum_, dataset::name, ordinate, muonCSCDigis_cfi::pedestal, EcalDataFrame::sample(), EcalDataFrame::size(), and fftjetcommon_cfi::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;
 }

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

private

Definition at line 235 of file EcalMipGraphs.cc.

References allFedsTimingHist_, edm::SortedCollection< T, SORT >::begin(), crysAndAmplitudesMap_, EcalElectronicsId::dccId(), ecalElectronicsMap_, edm::SortedCollection< T, SORT >::end(), EcalRecHit::energy(), FEDsAndTimingHists_, spr::find(), getEEIndex(), EcalElectronicsMapping::getElectronicsId(), CaloTopology::getSubdetectorTopology(), CaloNavigator< T, TOPO >::home(), mps_fire::i, EcalRecHit::id(), initHists(), intToString(), GeomDetEnumerators::isBarrel(), maskedChannels_, maskedFEDs_, minTimingAmp_, naiveEvtNum_, dataset::name, names, CaloNavigator< T, TOPO >::offsetBy(), mps_fire::result, SurveyInfoScenario_cff::seed, seedCrys_, seedFrequencyMap_, selectDigi(), side_, threshold_, EcalRecHit::time(), and fftjetcommon_cfi::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());
     }
   }
 }

void EcalMipGraphs::writeGraphs ( )

private

Member Data Documentation

int EcalMipGraphs::abscissa[10]

private

Definition at line 95 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), and selectDigi().

TH1F* EcalMipGraphs::allFedsTimingHist_

private

Definition at line 111 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), and selectHits().

TTree* EcalMipGraphs::canvasNames_

private

Definition at line 114 of file EcalMipGraphs.h.

Referenced by analyze(), and endJob().

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

private

Definition at line 107 of file EcalMipGraphs.h.

Referenced by selectDigi(), and selectHits().

edm::InputTag EcalMipGraphs::EBDigis_

private

Definition at line 80 of file EcalMipGraphs.h.

Referenced by analyze().

edm::Handle<EBDigiCollection> EcalMipGraphs::EBdigisHandle

private

Definition at line 84 of file EcalMipGraphs.h.

Referenced by analyze(), and selectDigi().

edm::InputTag EcalMipGraphs::EBRecHitCollection_

private

Definition at line 78 of file EcalMipGraphs.h.

Referenced by analyze().

const EcalElectronicsMapping* EcalMipGraphs::ecalElectronicsMap_

private

Definition at line 116 of file EcalMipGraphs.h.

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

edm::InputTag EcalMipGraphs::EEDigis_

private

Definition at line 81 of file EcalMipGraphs.h.

Referenced by analyze().

edm::Handle<EEDigiCollection> EcalMipGraphs::EEdigisHandle

private

Definition at line 85 of file EcalMipGraphs.h.

Referenced by analyze(), and selectDigi().

edm::InputTag EcalMipGraphs::EERecHitCollection_

private

Definition at line 79 of file EcalMipGraphs.h.

Referenced by analyze().

EcalFedMap* EcalMipGraphs::fedMap_

private

Definition at line 115 of file EcalMipGraphs.h.

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

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

private

Definition at line 108 of file EcalMipGraphs.h.

Referenced by analyze(), and selectDigi().

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

private

Definition at line 106 of file EcalMipGraphs.h.

Referenced by initHists(), and selectHits().

TFile* EcalMipGraphs::file_

private

Definition at line 113 of file EcalMipGraphs.h.

edm::Service< TFileService > EcalMipGraphs::fileService

staticprivate

Definition at line 99 of file EcalMipGraphs.h.

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

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

staticprivate

Definition at line 98 of file EcalMipGraphs.h.

Referenced by selectDigi().

edm::InputTag EcalMipGraphs::headerProducer_

private

Definition at line 82 of file EcalMipGraphs.h.

Referenced by analyze().

std::set<EBDetId> EcalMipGraphs::listEBChannels

private

Definition at line 92 of file EcalMipGraphs.h.

Referenced by analyze().

std::set<EEDetId> EcalMipGraphs::listEEChannels

private

Definition at line 93 of file EcalMipGraphs.h.

Referenced by analyze().

std::vector<int> EcalMipGraphs::maskedChannels_

private

Definition at line 102 of file EcalMipGraphs.h.

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

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

private

Definition at line 105 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs().

std::vector<int> EcalMipGraphs::maskedFEDs_

private

Definition at line 103 of file EcalMipGraphs.h.

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

double EcalMipGraphs::minTimingAmp_

private

Definition at line 90 of file EcalMipGraphs.h.

Referenced by selectHits().

int EcalMipGraphs::naiveEvtNum_

private

Definition at line 118 of file EcalMipGraphs.h.

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

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

private

Definition at line 101 of file EcalMipGraphs.h.

Referenced by analyze(), and selectHits().

int EcalMipGraphs::ordinate[10]

private

Definition at line 96 of file EcalMipGraphs.h.

Referenced by selectDigi().

int EcalMipGraphs::runNum_

private

Definition at line 87 of file EcalMipGraphs.h.

Referenced by analyze().

std::vector<int> EcalMipGraphs::seedCrys_

private

Definition at line 104 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), and selectHits().

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

private

Definition at line 109 of file EcalMipGraphs.h.

Referenced by endJob(), and selectHits().

int EcalMipGraphs::side_

private

Definition at line 88 of file EcalMipGraphs.h.

Referenced by selectHits().

double EcalMipGraphs::threshold_

private

Definition at line 89 of file EcalMipGraphs.h.

Referenced by selectHits().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation