EcalSimple2007H4TBAnalyzer Class Reference

#include <EcalSimple2007H4TBAnalyzer.h>

Inheritance diagram for EcalSimple2007H4TBAnalyzer:

Public Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
void	beginRun (edm::Run const &, edm::EventSetup const &) override
	EcalSimple2007H4TBAnalyzer (const edm::ParameterSet &)
void	endJob () override
void	endRun (edm::Run const &, edm::EventSetup const &) override
	~EcalSimple2007H4TBAnalyzer () 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 Attributes
const std::string	digiCollection_
const std::string	digiProducer_
const edm::EDGetTokenT < EEDigiCollection >	eeDigiToken_
const edm::EDGetTokenT < EEUncalibratedRecHitCollection >	eeUncalibratedRecHitToken_
const std::string	eventHeaderCollection_
const std::string	eventHeaderProducer_
const edm::ESGetToken < CaloGeometry, CaloGeometryRecord >	geometryToken_
TH2F *	h_ampltdc
TH1F *	h_bprofx
TH1F *	h_bprofy
TH1F *	h_e1e25
TH2F *	h_e1e25_mapx
TH2F *	h_e1e25_mapy
TH1F *	h_e1e9
TH2F *	h_e1e9_mapx
TH2F *	h_e1e9_mapy
TH1F *	h_e1x1
TH1F *	h_e1x1_center
TH1F *	h_e3x3
TH1F *	h_e3x3_center
TH1F *	h_e5x5
TH1F *	h_e5x5_center
TH1F *	h_e9e25
TH2F *	h_e9e25_mapx
TH2F *	h_e9e25_mapy
TH2F *	h_mapx [25]
TH2F *	h_mapy [25]
TH1F *	h_qualx
TH1F *	h_qualy
TH1F *	h_S6
TH2F *	h_Shape_
TH1F *	h_slopex
TH1F *	h_slopey
TH1F *	h_tableIsMoving
const std::string	hitCollection_
const std::string	hitProducer_
const std::string	hodoRecInfoCollection_
const std::string	hodoRecInfoProducer_
const std::string	rootfile_
const edm::EDGetTokenT < EcalTBEventHeader >	tbEventHeaderToken_
const edm::EDGetTokenT < EcalTBHodoscopeRecInfo >	tbHodoscopeRecInfoToken_
const edm::EDGetTokenT < EcalTBTDCRecInfo >	tbTDCRecInfoToken_
const std::string	tdcRecInfoCollection_
const std::string	tdcRecInfoProducer_
const CaloGeometry *	theTBGeometry_
EEDetId	xtalInBeam_
EBDetId	xtalInBeamTmp
EEDetId	Xtals5x5 [25]

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

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

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

Definition at line 43 of file EcalSimple2007H4TBAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 36 of file EcalSimple2007H4TBAnalyzer.cc.

References hitCollection_, and hitProducer_.

    : rootfile_(iConfig.getUntrackedParameter<std::string>("rootfile", "ecalSimpleTBanalysis.root")),
      digiCollection_(iConfig.getParameter<std::string>("digiCollection")),
      digiProducer_(iConfig.getParameter<std::string>("digiProducer")),
      hitCollection_(iConfig.getParameter<std::string>("hitCollection")),
      hitProducer_(iConfig.getParameter<std::string>("hitProducer")),
      hodoRecInfoCollection_(iConfig.getParameter<std::string>("hodoRecInfoCollection")),
      hodoRecInfoProducer_(iConfig.getParameter<std::string>("hodoRecInfoProducer")),
      tdcRecInfoCollection_(iConfig.getParameter<std::string>("tdcRecInfoCollection")),
      tdcRecInfoProducer_(iConfig.getParameter<std::string>("tdcRecInfoProducer")),
      eventHeaderCollection_(iConfig.getParameter<std::string>("eventHeaderCollection")),
      eventHeaderProducer_(iConfig.getParameter<std::string>("eventHeaderProducer")),
      eeDigiToken_(consumes<>(edm::InputTag(digiProducer_, digiCollection_))),
      eeUncalibratedRecHitToken_(consumes<>(edm::InputTag(hitProducer_, hitCollection_))),
      tbHodoscopeRecInfoToken_(consumes<>(edm::InputTag(hodoRecInfoProducer_, hodoRecInfoCollection_))),
      tbTDCRecInfoToken_(consumes<>(edm::InputTag(tdcRecInfoProducer_, tdcRecInfoCollection_))),
      tbEventHeaderToken_(consumes<>(edm::InputTag(eventHeaderProducer_))),
      geometryToken_(esConsumes<edm::Transition::BeginRun>()),
      xtalInBeam_(0)
//========================================================================
{
  //now do what ever initialization is needed
  edm::LogVerbatim("EcalSimple2007H4TBAnalyzer")
      << "EcalSimple2007H4TBAnalyzer: fetching hitCollection: " << hitCollection_.c_str() << " produced by "
      << hitProducer_.c_str();
}

EcalSimple2007H4TBAnalyzer::~EcalSimple2007H4TBAnalyzer ( )

override

Definition at line 64 of file EcalSimple2007H4TBAnalyzer.cc.

{
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
  // Amplitude vs TDC offset
  //   if (h_ampltdc)
  //   delete h_ampltdc;

  //   // Reconstructed energies
  //   delete h_e1x1;
  //   delete h_e3x3;
  //   delete h_e5x5;

  //   delete h_bprofx;
  //   delete h_bprofy;

  //   delete h_qualx;
  //   delete h_qualy;

  //   delete h_slopex;
  //   delete h_slopey;

  //   delete h_mapx;
  //   delete h_mapy;
}

Member Function Documentation

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 211 of file EcalSimple2007H4TBAnalyzer.cc.

References EcalMGPASample::adc(), digiCollection_, eeDigiToken_, eeUncalibratedRecHitToken_, eventHeaderProducer_, edm::Event::getByToken(), CaloGeometry::getGeometry(), h_ampltdc, h_bprofx, h_bprofy, h_e1e25, h_e1e25_mapx, h_e1e25_mapy, h_e1e9, h_e1e9_mapx, h_e1e9_mapy, h_e1x1, h_e1x1_center, h_e3x3, h_e3x3_center, h_e5x5, h_e5x5_center, h_e9e25, h_e9e25_mapx, h_e9e25_mapy, h_mapx, h_mapy, h_qualx, h_qualy, h_S6, h_Shape_, h_slopex, h_slopey, h_tableIsMoving, hitCollection_, hodoRecInfoCollection_, mps_fire::i, EBDetId::ic(), DetId::null(), phits, makeMEIFBenchmarkPlots::sample, EcalDataFrame::sample(), EcalDataFrame::size(), EBDetId::SMCRYSTALMODE, tbEventHeaderToken_, tbHodoscopeRecInfoToken_, tbTDCRecInfoToken_, tdcRecInfoCollection_, theTBGeometry_, x, xtalInBeam_, xtalInBeamTmp, Xtals5x5, y, and EBDetId::zside().

                                                                                            {
  //========================================================================

  using namespace edm;
  using namespace cms;

  Handle<EEDigiCollection> pdigis;
  const EEDigiCollection* digis = nullptr;
  iEvent.getByToken(eeDigiToken_, pdigis);
  if (pdigis.isValid()) {
    digis = pdigis.product();  // get a ptr to the product
  } else {
    edm::LogError("EcalSimple2007H4TBAnalyzerError") << "Error! can't get the product " << digiCollection_;
  }

  // fetch the digis and compute signal amplitude
  Handle<EEUncalibratedRecHitCollection> phits;
  const EEUncalibratedRecHitCollection* hits = nullptr;
  iEvent.getByToken(eeUncalibratedRecHitToken_, phits);
  if (phits.isValid()) {
    hits = phits.product();  // get a ptr to the product
  } else {
    edm::LogError("EcalSimple2007H4TBAnalyzerError") << "Error! can't get the product " << hitCollection_;
  }

  Handle<EcalTBHodoscopeRecInfo> pHodo;
  const EcalTBHodoscopeRecInfo* recHodo = nullptr;
  iEvent.getByToken(tbHodoscopeRecInfoToken_, pHodo);
  if (pHodo.isValid()) {
    recHodo = pHodo.product();  // get a ptr to the product
  } else {
    edm::LogError("EcalSimple2007H4TBAnalyzerError") << "Error! can't get the product " << hodoRecInfoCollection_;
  }

  Handle<EcalTBTDCRecInfo> pTDC;
  const EcalTBTDCRecInfo* recTDC = nullptr;
  iEvent.getByToken(tbTDCRecInfoToken_, pTDC);
  if (pTDC.isValid()) {
    recTDC = pTDC.product();  // get a ptr to the product
  } else {
    edm::LogError("EcalSimple2007H4TBAnalyzerError") << "Error! can't get the product " << tdcRecInfoCollection_;
  }

  Handle<EcalTBEventHeader> pEventHeader;
  const EcalTBEventHeader* evtHeader = nullptr;
  iEvent.getByToken(tbEventHeaderToken_, pEventHeader);
  if (pEventHeader.isValid()) {
    evtHeader = pEventHeader.product();  // get a ptr to the product
  } else {
    edm::LogError("EcalSimple2007H4TBAnalyzerError") << "Error! can't get the product " << eventHeaderProducer_;
  }

  if (!hits) {
    return;
  }

  if (!recTDC) {
    return;
  }

  if (!recHodo) {
    return;
  }

  if (!evtHeader) {
    return;
  }

  if (hits->empty()) {
    return;
  }

  //Accessing various event information
  if (evtHeader->tableIsMoving())
    h_tableIsMoving->Fill(evtHeader->eventNumber());

  //S6 beam scintillator
  h_S6->Fill(evtHeader->S6ADC());

  if (xtalInBeamTmp.null()) {
    xtalInBeamTmp = EBDetId(1, evtHeader->crystalInBeam(), EBDetId::SMCRYSTALMODE);
    xtalInBeam_ = EEDetId(35 - ((xtalInBeamTmp.ic() - 1) % 20), int(int(xtalInBeamTmp.ic()) / int(20)) + 51, -1);
    edm::LogVerbatim("EcalSimple2007H4TBAnalyzer") << "Xtal In Beam is " << xtalInBeam_.ic() << xtalInBeam_;
    for (unsigned int icry = 0; icry < 25; icry++) {
      unsigned int row = icry / 5;
      unsigned int column = icry % 5;
      int ix = xtalInBeam_.ix() + row - 2;
      int iy = xtalInBeam_.iy() + column - 2;
      EEDetId tempId(ix, iy, xtalInBeam_.zside());
      //Selecting matrix of xtals used in 2007H4TB
      if (tempId.ix() < 16 || tempId.ix() > 35 || tempId.iy() < 51 || tempId.iy() > 75)
        Xtals5x5[icry] = EEDetId(0);
      else {
        Xtals5x5[icry] = tempId;
        auto cell = theTBGeometry_->getGeometry(Xtals5x5[icry]);
        if (!cell)
          continue;
        edm::LogVerbatim("EcalSimple2007H4TBAnalyzer")
            << "** Xtal in the matrix **** row " << row << ", column " << column << ", xtal " << Xtals5x5[icry]
            << " Position " << cell->getPosition(0.);
      }
    }
  } else if (xtalInBeamTmp !=
             EBDetId(1, evtHeader->crystalInBeam(), EBDetId::SMCRYSTALMODE))  //run analysis only on first xtal in beam
    return;

  //Avoid moving table events
  if (evtHeader->tableIsMoving()) {
    edm::LogVerbatim("EcalSimple2007H4TBAnalyzer") << "Table is moving";
    return;
  }

  // Searching for max amplitude xtal alternative to use xtalInBeam_
  EEDetId maxHitId(0);
  float maxHit = -999999.;
  for (EEUncalibratedRecHitCollection::const_iterator ithit = hits->begin(); ithit != hits->end(); ++ithit) {
    if (ithit->amplitude() >= maxHit) {
      maxHit = ithit->amplitude();
      maxHitId = ithit->id();
    }
  }
  if (maxHitId == EEDetId(0)) {
    edm::LogVerbatim("EcalSimple2007H4TBAnalyzer") << "No maxHit found";
    return;
  }

  //Filling the digis shape for the xtalInBeam
  double samples_save[10];
  for (int i = 0; i < 10; ++i)
    samples_save[i] = 0.0;

  double eMax = 0.;
  for (EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr) {
    if (EEDetId((*digiItr).id()) != xtalInBeam_)
      continue;

    EEDataFrame myDigi = (*digiItr);
    for (int sample = 0; sample < myDigi.size(); ++sample) {
      double analogSample = myDigi.sample(sample).adc();
      samples_save[sample] = analogSample;
      if (eMax < analogSample) {
        eMax = analogSample;
      }
    }
  }

  for (int i = 0; i < 10; ++i)
    h_Shape_->Fill(double(i) + recTDC->offset(), samples_save[i]);

  // Taking amplitudes in 5x5
  double amplitude[25];
  double amplitude3x3 = 0;
  double amplitude5x5 = 0;
  for (unsigned int icry = 0; icry < 25; icry++) {
    if (!Xtals5x5[icry].null()) {
      amplitude[icry] = (hits->find(Xtals5x5[icry]))->amplitude();
      amplitude5x5 += amplitude[icry];
      // Is in 3x3?
      if (icry == 6 || icry == 7 || icry == 8 || icry == 11 || icry == 12 || icry == 13 || icry == 16 || icry == 17 ||
          icry == 18) {
        amplitude3x3 += amplitude[icry];
      }
    }
  }

  //Filling amplitudes
  h_e1x1->Fill(amplitude[12]);
  h_e3x3->Fill(amplitude3x3);
  h_e5x5->Fill(amplitude5x5);

  h_e1e9->Fill(amplitude[12] / amplitude3x3);
  h_e1e25->Fill(amplitude[12] / amplitude5x5);
  h_e9e25->Fill(amplitude3x3 / amplitude5x5);

  //Checking stability of amplitude vs TDC
  if (recTDC)
    h_ampltdc->Fill(recTDC->offset(), amplitude[12]);

  //Various amplitudes as a function of hodoscope coordinates
  if (recHodo) {
    float x = recHodo->posX();
    float y = recHodo->posY();
    float xslope = recHodo->slopeX();
    float yslope = recHodo->slopeY();
    float xqual = recHodo->qualX();
    float yqual = recHodo->qualY();

    //Filling beam profiles
    h_bprofx->Fill(x);
    h_bprofy->Fill(y);
    h_qualx->Fill(xqual);
    h_qualy->Fill(yqual);
    h_slopex->Fill(xslope);
    h_slopey->Fill(yslope);

    //Fill central events

    if (fabs(x + 2.5) < 2.5 && fabs(y + 0.5) < 2.5) {
      h_e1x1_center->Fill(amplitude[12]);
      h_e3x3_center->Fill(amplitude3x3);
      h_e5x5_center->Fill(amplitude5x5);
    }

    for (unsigned int icry = 0; icry < 25; icry++) {
      h_mapx[icry]->Fill(x, amplitude[icry]);
      h_mapy[icry]->Fill(y, amplitude[icry]);
    }

    h_e1e9_mapx->Fill(x, amplitude[12] / amplitude3x3);
    h_e1e9_mapy->Fill(y, amplitude[12] / amplitude3x3);

    h_e1e25_mapx->Fill(x, amplitude[12] / amplitude5x5);
    h_e1e25_mapy->Fill(y, amplitude[12] / amplitude5x5);

    h_e9e25_mapx->Fill(x, amplitude3x3 / amplitude5x5);
    h_e9e25_mapy->Fill(y, amplitude3x3 / amplitude5x5);
  }
}

void EcalSimple2007H4TBAnalyzer::beginRun ( edm::Run const &  , edm::EventSetup const &  iSetup  )

override

Definition at line 92 of file EcalSimple2007H4TBAnalyzer.cc.

References alignCSCRings::e, geometryToken_, edm::EventSetup::getData(), h_ampltdc, h_bprofx, h_bprofy, h_e1e25, h_e1e25_mapx, h_e1e25_mapy, h_e1e9, h_e1e9_mapx, h_e1e9_mapy, h_e1x1, h_e1x1_center, h_e3x3, h_e3x3_center, h_e5x5, h_e5x5_center, h_e9e25, h_e9e25_mapx, h_e9e25_mapy, h_mapx, h_mapy, h_qualx, h_qualy, h_S6, h_Shape_, h_slopex, h_slopey, h_tableIsMoving, and theTBGeometry_.

                                                                                    {
  //========================================================================

  theTBGeometry_ = &iSetup.getData(geometryToken_);

  // Amplitude vs TDC offset
  h_ampltdc = new TH2F("h_ampltdc", "Max Amplitude vs TDC offset", 100, 0., 1., 1000, 0., 4000.);

  // Reconstructed energies
  h_tableIsMoving = new TH1F("h_tableIsMoving", "TableIsMoving", 100000, 0., 100000.);

  h_e1x1 = new TH1F("h_e1x1", "E1x1 energy", 1000, 0., 4000.);
  h_e3x3 = new TH1F("h_e3x3", "E3x3 energy", 1000, 0., 4000.);
  h_e5x5 = new TH1F("h_e5x5", "E5x5 energy", 1000, 0., 4000.);

  h_e1x1_center = new TH1F("h_e1x1_center", "E1x1 energy", 1000, 0., 4000.);
  h_e3x3_center = new TH1F("h_e3x3_center", "E3x3 energy", 1000, 0., 4000.);
  h_e5x5_center = new TH1F("h_e5x5_center", "E5x5 energy", 1000, 0., 4000.);

  h_e1e9 = new TH1F("h_e1e9", "E1/E9 ratio", 600, 0., 1.2);
  h_e1e25 = new TH1F("h_e1e25", "E1/E25 ratio", 600, 0., 1.2);
  h_e9e25 = new TH1F("h_e9e25", "E9/E25 ratio", 600, 0., 1.2);

  h_S6 = new TH1F("h_S6", "Amplitude S6", 1000, 0., 4000.);

  h_bprofx = new TH1F("h_bprofx", "Beam Profile X", 100, -20., 20.);
  h_bprofy = new TH1F("h_bprofy", "Beam Profile Y", 100, -20., 20.);

  h_qualx = new TH1F("h_qualx", "Beam Quality X", 5000, 0., 5.);
  h_qualy = new TH1F("h_qualy", "Beam Quality X", 5000, 0., 5.);

  h_slopex = new TH1F("h_slopex", "Beam Slope X", 500, -5e-4, 5e-4);
  h_slopey = new TH1F("h_slopey", "Beam Slope Y", 500, -5e-4, 5e-4);

  char hname[50];
  char htitle[50];
  for (unsigned int icry = 0; icry < 25; icry++) {
    sprintf(hname, "h_mapx_%d", icry);
    sprintf(htitle, "Max Amplitude vs X %d", icry);
    h_mapx[icry] = new TH2F(hname, htitle, 80, -20, 20, 1000, 0., 4000.);
    sprintf(hname, "h_mapy_%d", icry);
    sprintf(htitle, "Max Amplitude vs Y %d", icry);
    h_mapy[icry] = new TH2F(hname, htitle, 80, -20, 20, 1000, 0., 4000.);
  }

  h_e1e9_mapx = new TH2F("h_e1e9_mapx", "E1/E9 vs X", 80, -20, 20, 600, 0., 1.2);
  h_e1e9_mapy = new TH2F("h_e1e9_mapy", "E1/E9 vs Y", 80, -20, 20, 600, 0., 1.2);

  h_e1e25_mapx = new TH2F("h_e1e25_mapx", "E1/E25 vs X", 80, -20, 20, 600, 0., 1.2);
  h_e1e25_mapy = new TH2F("h_e1e25_mapy", "E1/E25 vs Y", 80, -20, 20, 600, 0., 1.2);

  h_e9e25_mapx = new TH2F("h_e9e25_mapx", "E9/E25 vs X", 80, -20, 20, 600, 0., 1.2);
  h_e9e25_mapy = new TH2F("h_e9e25_mapy", "E9/E25 vs Y", 80, -20, 20, 600, 0., 1.2);

  h_Shape_ = new TH2F("h_Shape_", "Xtal in Beam Shape", 250, 0, 10, 350, 0, 3500);
}

void EcalSimple2007H4TBAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 154 of file EcalSimple2007H4TBAnalyzer.cc.

References validate-o2o-wbm::f, h_ampltdc, h_bprofx, h_bprofy, h_e1e25, h_e1e25_mapx, h_e1e25_mapy, h_e1e9, h_e1e9_mapx, h_e1e9_mapy, h_e1x1, h_e1x1_center, h_e3x3, h_e3x3_center, h_e5x5, h_e5x5_center, h_e9e25, h_e9e25_mapx, h_e9e25_mapy, h_mapx, h_mapy, h_qualx, h_qualy, h_S6, h_Shape_, h_slopex, h_slopey, h_tableIsMoving, and rootfile_.

                                        {
  //========================================================================

  TFile f(rootfile_.c_str(), "RECREATE");

  // Amplitude vs TDC offset
  h_ampltdc->Write();

  // Reconstructed energies
  h_e1x1->Write();
  h_e3x3->Write();
  h_e5x5->Write();

  h_e1x1_center->Write();
  h_e3x3_center->Write();
  h_e5x5_center->Write();

  h_e1e9->Write();
  h_e1e25->Write();
  h_e9e25->Write();

  h_S6->Write();
  h_bprofx->Write();
  h_bprofy->Write();

  h_qualx->Write();
  h_qualy->Write();

  h_slopex->Write();
  h_slopey->Write();

  h_Shape_->Write();

  for (unsigned int icry = 0; icry < 25; icry++) {
    h_mapx[icry]->Write();
    h_mapy[icry]->Write();
  }

  h_e1e9_mapx->Write();
  h_e1e9_mapy->Write();

  h_e1e25_mapx->Write();
  h_e1e25_mapy->Write();

  h_e9e25_mapx->Write();
  h_e9e25_mapy->Write();

  h_tableIsMoving->Write();

  f.Close();
}

void EcalSimple2007H4TBAnalyzer::endRun ( edm::Run const &  , edm::EventSetup const &  iSetup  )

override

Definition at line 150 of file EcalSimple2007H4TBAnalyzer.cc.

{}

Member Data Documentation

const std::string EcalSimple2007H4TBAnalyzer::digiCollection_

private

Definition at line 55 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

const std::string EcalSimple2007H4TBAnalyzer::digiProducer_

private

Definition at line 56 of file EcalSimple2007H4TBAnalyzer.h.

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

private

Definition at line 66 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

const edm::EDGetTokenT<EEUncalibratedRecHitCollection> EcalSimple2007H4TBAnalyzer::eeUncalibratedRecHitToken_

private

Definition at line 67 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

const std::string EcalSimple2007H4TBAnalyzer::eventHeaderCollection_

private

Definition at line 63 of file EcalSimple2007H4TBAnalyzer.h.

const std::string EcalSimple2007H4TBAnalyzer::eventHeaderProducer_

private

Definition at line 64 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> EcalSimple2007H4TBAnalyzer::geometryToken_

private

Definition at line 71 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by beginRun().

TH2F* EcalSimple2007H4TBAnalyzer::h_ampltdc

private

Definition at line 74 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_bprofx

private

Definition at line 93 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_bprofy

private

Definition at line 94 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_e1e25

private

Definition at line 89 of file EcalSimple2007H4TBAnalyzer.h.

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

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e25_mapx

private

Definition at line 108 of file EcalSimple2007H4TBAnalyzer.h.

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

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e25_mapy

private

Definition at line 109 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_e1e9

private

Definition at line 88 of file EcalSimple2007H4TBAnalyzer.h.

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

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e9_mapx

private

Definition at line 105 of file EcalSimple2007H4TBAnalyzer.h.

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

TH2F* EcalSimple2007H4TBAnalyzer::h_e1e9_mapy

private

Definition at line 106 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_e1x1

private

Definition at line 80 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_e1x1_center

private

Definition at line 84 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_e3x3

private

Definition at line 81 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_e3x3_center

private

Definition at line 85 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_e5x5

private

Definition at line 82 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_e5x5_center

private

Definition at line 86 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_e9e25

private

Definition at line 90 of file EcalSimple2007H4TBAnalyzer.h.

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

TH2F* EcalSimple2007H4TBAnalyzer::h_e9e25_mapx

private

Definition at line 111 of file EcalSimple2007H4TBAnalyzer.h.

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

TH2F* EcalSimple2007H4TBAnalyzer::h_e9e25_mapy

private

Definition at line 112 of file EcalSimple2007H4TBAnalyzer.h.

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

TH2F* EcalSimple2007H4TBAnalyzer::h_mapx[25]

private

Definition at line 102 of file EcalSimple2007H4TBAnalyzer.h.

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

TH2F* EcalSimple2007H4TBAnalyzer::h_mapy[25]

private

Definition at line 103 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_qualx

private

Definition at line 96 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_qualy

private

Definition at line 97 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_S6

private

Definition at line 92 of file EcalSimple2007H4TBAnalyzer.h.

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

TH2F* EcalSimple2007H4TBAnalyzer::h_Shape_

private

Definition at line 76 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_slopex

private

Definition at line 99 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_slopey

private

Definition at line 100 of file EcalSimple2007H4TBAnalyzer.h.

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

TH1F* EcalSimple2007H4TBAnalyzer::h_tableIsMoving

private

Definition at line 79 of file EcalSimple2007H4TBAnalyzer.h.

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

const std::string EcalSimple2007H4TBAnalyzer::hitCollection_

private

Definition at line 57 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and EcalSimple2007H4TBAnalyzer().

const std::string EcalSimple2007H4TBAnalyzer::hitProducer_

private

Definition at line 58 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by EcalSimple2007H4TBAnalyzer().

const std::string EcalSimple2007H4TBAnalyzer::hodoRecInfoCollection_

private

Definition at line 59 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

const std::string EcalSimple2007H4TBAnalyzer::hodoRecInfoProducer_

private

Definition at line 60 of file EcalSimple2007H4TBAnalyzer.h.

const std::string EcalSimple2007H4TBAnalyzer::rootfile_

private

Definition at line 54 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by endJob().

const edm::EDGetTokenT<EcalTBEventHeader> EcalSimple2007H4TBAnalyzer::tbEventHeaderToken_

private

Definition at line 70 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

const edm::EDGetTokenT<EcalTBHodoscopeRecInfo> EcalSimple2007H4TBAnalyzer::tbHodoscopeRecInfoToken_

private

Definition at line 68 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

const edm::EDGetTokenT<EcalTBTDCRecInfo> EcalSimple2007H4TBAnalyzer::tbTDCRecInfoToken_

private

Definition at line 69 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

const std::string EcalSimple2007H4TBAnalyzer::tdcRecInfoCollection_

private

Definition at line 61 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

const std::string EcalSimple2007H4TBAnalyzer::tdcRecInfoProducer_

private

Definition at line 62 of file EcalSimple2007H4TBAnalyzer.h.

const CaloGeometry* EcalSimple2007H4TBAnalyzer::theTBGeometry_

private

Definition at line 118 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze(), and beginRun().

EEDetId EcalSimple2007H4TBAnalyzer::xtalInBeam_

private

Definition at line 114 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

EBDetId EcalSimple2007H4TBAnalyzer::xtalInBeamTmp

private

Definition at line 115 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().

EEDetId EcalSimple2007H4TBAnalyzer::Xtals5x5[25]

private

Definition at line 116 of file EcalSimple2007H4TBAnalyzer.h.

Referenced by analyze().