EcalMipGraphs Class Reference

#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 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 (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 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 () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
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 ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 40 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")),
      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 ( )

override

Definition at line 91 of file EcalMipGraphs.cc.

{}

Member Function Documentation

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 caloTopo_, canvasNames_, EBdigisHandle, ebDigiToken_, ebRecHitToken_, EEdigisHandle, eeDigiToken_, eeRecHitToken_, edm::EventID::event(), FEDsAndDCCHeaders_, fileService, edm::Event::getByToken(), edm::EventSetup::getData(), edm::EventBase::id(), listEBChannels, listEEChannels, TFileService::make(), naiveEvtNum_, names, rawDataToken_, edm::EventID::run(), 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);
}

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

overrideprivate

Definition at line 332 of file EcalMipGraphs.cc.

References ecalElectronicsMap_, ecalMappingToken_, and edm::EventSetup::getData().

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

void EcalMipGraphs::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 339 of file EcalMipGraphs.cc.

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

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

overrideprivate

Definition at line 336 of file EcalMipGraphs.cc.

{}

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

private

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

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;
}

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

private

Definition at line 372 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 136 of file EcalMipGraphs.cc.

References abscissa, EcalMGPASample::adc(), crysAndAmplitudesMap_, EcalElectronicsId::dccId(), DetId::det(), isotrackApplyRegressor::df, EBdigisHandle, DetId::Ecal, ecalElectronicsMap_, EEdigisHandle, relval_parameters_module::energy, fedMap_, FEDsAndDCCHeaders_, fileService, floatToString(), ecalLiteDTU::gainId(), EcalMGPASample::gainId(), gainRatio, getEEIndex(), EcalElectronicsMapping::getElectronicsId(), EcalFedMap::getSliceFromFed(), mps_fire::i, intToString(), PixelPluginsPhase0_cfi::isBarrel, TFileService::make(), naiveEvtNum_, mergeVDriftHistosByStation::name, ordinate, EcalCondDBWriter_cfi::pedestal, EcalDataFrame::sample(), EcalDataFrame::size(), 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;
}

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_, EcalRecHit::energy(), FEDsAndTimingHists_, spr::find(), getEEIndex(), EcalElectronicsMapping::getElectronicsId(), CaloTopology::getSubdetectorTopology(), CaloNavigator< T, TOPO >::home(), mps_fire::i, EcalRecHit::id(), initHists(), intToString(), PixelPluginsPhase0_cfi::isBarrel, dqmiolumiharvest::j, maskedChannels_, maskedFEDs_, minTimingAmp_, naiveEvtNum_, mergeVDriftHistosByStation::name, names, CaloNavigator< T, TOPO >::offsetBy(), mps_fire::result, fileCollector::seed, seedCrys_, seedFrequencyMap_, selectDigi(), side_, threshold_, EcalRecHit::time(), 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());
    }
  }
}

void EcalMipGraphs::writeGraphs ( )

private

Member Data Documentation

int EcalMipGraphs::abscissa[10]

private

Definition at line 105 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), and selectDigi().

TH1F* EcalMipGraphs::allFedsTimingHist_

private

Definition at line 121 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), and selectHits().

const CaloTopology* EcalMipGraphs::caloTopo_

private

Definition at line 127 of file EcalMipGraphs.h.

Referenced by analyze(), and selectHits().

TTree* EcalMipGraphs::canvasNames_

private

Definition at line 124 of file EcalMipGraphs.h.

Referenced by analyze(), and endJob().

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

private

Definition at line 117 of file EcalMipGraphs.h.

Referenced by selectDigi(), and selectHits().

const edm::InputTag EcalMipGraphs::EBDigis_

private

Definition at line 81 of file EcalMipGraphs.h.

edm::Handle<EBDigiCollection> EcalMipGraphs::EBdigisHandle

private

Definition at line 85 of file EcalMipGraphs.h.

Referenced by analyze(), and selectDigi().

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

private

Definition at line 91 of file EcalMipGraphs.h.

Referenced by analyze().

const edm::InputTag EcalMipGraphs::EBRecHitCollection_

private

Definition at line 79 of file EcalMipGraphs.h.

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

private

Definition at line 89 of file EcalMipGraphs.h.

Referenced by analyze().

const EcalElectronicsMapping* EcalMipGraphs::ecalElectronicsMap_

private

Definition at line 126 of file EcalMipGraphs.h.

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

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

private

Definition at line 94 of file EcalMipGraphs.h.

Referenced by beginRun().

const edm::InputTag EcalMipGraphs::EEDigis_

private

Definition at line 82 of file EcalMipGraphs.h.

edm::Handle<EEDigiCollection> EcalMipGraphs::EEdigisHandle

private

Definition at line 86 of file EcalMipGraphs.h.

Referenced by analyze(), and selectDigi().

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

private

Definition at line 92 of file EcalMipGraphs.h.

Referenced by analyze().

const edm::InputTag EcalMipGraphs::EERecHitCollection_

private

Definition at line 80 of file EcalMipGraphs.h.

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

private

Definition at line 90 of file EcalMipGraphs.h.

Referenced by analyze().

EcalFedMap* EcalMipGraphs::fedMap_

private

Definition at line 125 of file EcalMipGraphs.h.

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

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

private

Definition at line 118 of file EcalMipGraphs.h.

Referenced by analyze(), and selectDigi().

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

private

Definition at line 116 of file EcalMipGraphs.h.

Referenced by initHists(), and selectHits().

TFile* EcalMipGraphs::file_

private

Definition at line 123 of file EcalMipGraphs.h.

edm::Service< TFileService > EcalMipGraphs::fileService

staticprivate

Definition at line 109 of file EcalMipGraphs.h.

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

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

staticprivate

Definition at line 108 of file EcalMipGraphs.h.

Referenced by selectDigi().

const edm::InputTag EcalMipGraphs::headerProducer_

private

Definition at line 83 of file EcalMipGraphs.h.

std::set<EBDetId> EcalMipGraphs::listEBChannels

private

Definition at line 102 of file EcalMipGraphs.h.

Referenced by analyze().

std::set<EEDetId> EcalMipGraphs::listEEChannels

private

Definition at line 103 of file EcalMipGraphs.h.

Referenced by analyze().

std::vector<int> EcalMipGraphs::maskedChannels_

private

Definition at line 112 of file EcalMipGraphs.h.

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

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

private

Definition at line 115 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs().

std::vector<int> EcalMipGraphs::maskedFEDs_

private

Definition at line 113 of file EcalMipGraphs.h.

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

const double EcalMipGraphs::minTimingAmp_

private

Definition at line 100 of file EcalMipGraphs.h.

Referenced by selectHits().

int EcalMipGraphs::naiveEvtNum_

private

Definition at line 129 of file EcalMipGraphs.h.

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

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

private

Definition at line 111 of file EcalMipGraphs.h.

Referenced by analyze(), and selectHits().

int EcalMipGraphs::ordinate[10]

private

Definition at line 106 of file EcalMipGraphs.h.

Referenced by selectDigi().

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

private

Definition at line 88 of file EcalMipGraphs.h.

Referenced by analyze().

int EcalMipGraphs::runNum_

private

Definition at line 97 of file EcalMipGraphs.h.

Referenced by analyze().

std::vector<int> EcalMipGraphs::seedCrys_

private

Definition at line 114 of file EcalMipGraphs.h.

Referenced by EcalMipGraphs(), and selectHits().

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

private

Definition at line 119 of file EcalMipGraphs.h.

Referenced by endJob(), and selectHits().

const int EcalMipGraphs::side_

private

Definition at line 98 of file EcalMipGraphs.h.

Referenced by selectHits().

const double EcalMipGraphs::threshold_

private

Definition at line 99 of file EcalMipGraphs.h.

Referenced by selectHits().

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

private

Definition at line 95 of file EcalMipGraphs.h.

Referenced by analyze().