CastorPedestalsAnalysis Class Reference

Inheritance diagram for CastorPedestalsAnalysis:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override
	CastorPedestalsAnalysis (const edm::ParameterSet &ps)
	~CastorPedestalsAnalysis () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	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 Attributes
std::vector< NewPedBunch >	Bunches
const edm::EDGetTokenT< CastorDigiCollection >	castorDigiCollectionToken_
TH1F *	CASTORMeans
TH1F *	CASTORWidths
TH2F *	dephist
const bool	dumpXML
bool	firsttime
const int	firstTS
const bool	hiSaveFlag
const int	lastTS
std::string	pedsADCfilename
std::string	pedsfCfilename
std::string	ROOTfilename
int	runnum
const edm::ESGetToken< CastorDbService, CastorDbRecord >	tok_cond_
const edm::ESGetToken< CastorElectronicsMap, CastorElectronicsMapRcd >	tok_map_
const bool	verboseflag
std::string	widthsADCfilename
std::string	widthsfCfilename
std::string	XMLfilename
std::string	XMLtag
std::string	ZSfilename

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

Definition at line 74 of file CastorPedestalsAnalysis.cc.

Constructor & Destructor Documentation

CastorPedestalsAnalysis()

CastorPedestalsAnalysis::CastorPedestalsAnalysis

(

const edm::ParameterSet &

ps

)

Definition at line 113 of file CastorPedestalsAnalysis.cc.

References firsttime, and TFileService::kSharedResource.

    : castorDigiCollectionToken_(
          consumes<CastorDigiCollection>(ps.getParameter<edm::InputTag>("castorDigiCollectionTag"))),
      tok_cond_(esConsumes<CastorDbService, CastorDbRecord>()),
      tok_map_(esConsumes<CastorElectronicsMap, CastorElectronicsMapRcd>()),
      hiSaveFlag(ps.getUntrackedParameter<bool>("hiSaveFlag", false)),
      dumpXML(ps.getUntrackedParameter<bool>("dumpXML", false)),
      verboseflag(ps.getUntrackedParameter<bool>("verbose", false)),
      firstTS(ps.getUntrackedParameter<int>("firstTS", 0)),
      lastTS(ps.getUntrackedParameter<int>("lastTS", 9)) {
  usesResource(TFileService::kSharedResource);

  firsttime = true;
}

~CastorPedestalsAnalysis()

CastorPedestalsAnalysis::~CastorPedestalsAnalysis ( )

override

Definition at line 128 of file CastorPedestalsAnalysis.cc.

References CastorCondObjectContainer< Item >::addValues(), Bunches, CASTORMeans, CASTORWidths, gather_cfg::cout, dephist, CastorDbASCIIIO::dumpObject(), dumpXML, mps_fire::i, B2GTnPMonitor_cfi::item, hfnoseParametersInitialization_cfi::name2, pedsADCfilename, pedsfCfilename, runnum, CastorPedestalWidth::setSigma(), AlCaHLTBitMon_QueryRunRegistry::string, verboseflag, widthsADCfilename, widthsfCfilename, and XMLfilename.

                                                  {
  CastorPedestals* rawPedsItem = new CastorPedestals(true);
  CastorPedestalWidths* rawWidthsItem = new CastorPedestalWidths(true);
  CastorPedestals* rawPedsItemfc = new CastorPedestals(false);
  CastorPedestalWidths* rawWidthsItemfc = new CastorPedestalWidths(false);

  //Calculate pedestal constants
  std::cout << "Calculating Pedestal constants...\n";
  std::vector<NewPedBunch>::iterator bunch_it;
  for (bunch_it = Bunches.begin(); bunch_it != Bunches.end(); ++bunch_it) {
    if (bunch_it->usedflag) {
      if (verboseflag)
        std::cout << "Analyzing channel sector= " << bunch_it->detid.sector()
                  << " module = " << bunch_it->detid.module() << std::endl;
      //pedestal constant is the mean
      bunch_it->cap[0] /= bunch_it->num[0][0];
      bunch_it->cap[1] /= bunch_it->num[1][1];
      bunch_it->cap[2] /= bunch_it->num[2][2];
      bunch_it->cap[3] /= bunch_it->num[3][3];
      bunch_it->capfc[0] /= bunch_it->num[0][0];
      bunch_it->capfc[1] /= bunch_it->num[1][1];
      bunch_it->capfc[2] /= bunch_it->num[2][2];
      bunch_it->capfc[3] /= bunch_it->num[3][3];
      //widths are the covariance matrix--assumed symmetric
      bunch_it->sig[0][0] = (bunch_it->prod[0][0] / bunch_it->num[0][0]) - (bunch_it->cap[0] * bunch_it->cap[0]);
      bunch_it->sig[0][1] = (bunch_it->prod[0][1] / bunch_it->num[0][1]) - (bunch_it->cap[0] * bunch_it->cap[1]);
      bunch_it->sig[0][2] = (bunch_it->prod[0][2] / bunch_it->num[0][2]) - (bunch_it->cap[0] * bunch_it->cap[2]);
      bunch_it->sig[0][3] = (bunch_it->prod[0][3] / bunch_it->num[0][3]) - (bunch_it->cap[0] * bunch_it->cap[3]);
      bunch_it->sig[1][0] = (bunch_it->prod[1][0] / bunch_it->num[1][0]) - (bunch_it->cap[1] * bunch_it->cap[0]);
      bunch_it->sig[1][1] = (bunch_it->prod[1][1] / bunch_it->num[1][1]) - (bunch_it->cap[1] * bunch_it->cap[1]);
      bunch_it->sig[1][2] = (bunch_it->prod[1][2] / bunch_it->num[1][2]) - (bunch_it->cap[1] * bunch_it->cap[2]);
      bunch_it->sig[1][3] = (bunch_it->prod[1][3] / bunch_it->num[1][3]) - (bunch_it->cap[1] * bunch_it->cap[3]);
      bunch_it->sig[2][0] = (bunch_it->prod[2][0] / bunch_it->num[2][0]) - (bunch_it->cap[2] * bunch_it->cap[0]);
      bunch_it->sig[2][1] = (bunch_it->prod[2][1] / bunch_it->num[2][1]) - (bunch_it->cap[2] * bunch_it->cap[1]);
      bunch_it->sig[2][2] = (bunch_it->prod[2][2] / bunch_it->num[2][2]) - (bunch_it->cap[2] * bunch_it->cap[2]);
      bunch_it->sig[2][3] = (bunch_it->prod[2][3] / bunch_it->num[2][3]) - (bunch_it->cap[2] * bunch_it->cap[3]);
      bunch_it->sig[3][0] = (bunch_it->prod[3][0] / bunch_it->num[3][0]) - (bunch_it->cap[3] * bunch_it->cap[0]);
      bunch_it->sig[3][1] = (bunch_it->prod[3][1] / bunch_it->num[3][1]) - (bunch_it->cap[3] * bunch_it->cap[1]);
      bunch_it->sig[3][2] = (bunch_it->prod[3][2] / bunch_it->num[3][2]) - (bunch_it->cap[3] * bunch_it->cap[2]);
      bunch_it->sig[3][3] = (bunch_it->prod[3][3] / bunch_it->num[3][3]) - (bunch_it->cap[3] * bunch_it->cap[3]);

      bunch_it->sigfc[0][0] =
          (bunch_it->prodfc[0][0] / bunch_it->num[0][0]) - (bunch_it->capfc[0] * bunch_it->capfc[0]);
      bunch_it->sigfc[0][1] =
          (bunch_it->prodfc[0][1] / bunch_it->num[0][1]) - (bunch_it->capfc[0] * bunch_it->capfc[1]);
      bunch_it->sigfc[0][2] =
          (bunch_it->prodfc[0][2] / bunch_it->num[0][2]) - (bunch_it->capfc[0] * bunch_it->capfc[2]);
      bunch_it->sigfc[0][3] =
          (bunch_it->prodfc[0][3] / bunch_it->num[0][3]) - (bunch_it->capfc[0] * bunch_it->capfc[3]);
      bunch_it->sigfc[1][0] =
          (bunch_it->prodfc[1][0] / bunch_it->num[1][0]) - (bunch_it->capfc[1] * bunch_it->capfc[0]);
      bunch_it->sigfc[1][1] =
          (bunch_it->prodfc[1][1] / bunch_it->num[1][1]) - (bunch_it->capfc[1] * bunch_it->capfc[1]);
      bunch_it->sigfc[1][2] =
          (bunch_it->prodfc[1][2] / bunch_it->num[1][2]) - (bunch_it->capfc[1] * bunch_it->capfc[2]);
      bunch_it->sigfc[1][3] =
          (bunch_it->prodfc[1][3] / bunch_it->num[1][3]) - (bunch_it->capfc[1] * bunch_it->capfc[3]);
      bunch_it->sigfc[2][0] =
          (bunch_it->prodfc[2][0] / bunch_it->num[2][0]) - (bunch_it->capfc[2] * bunch_it->capfc[0]);
      bunch_it->sigfc[2][1] =
          (bunch_it->prodfc[2][1] / bunch_it->num[2][1]) - (bunch_it->capfc[2] * bunch_it->capfc[1]);
      bunch_it->sigfc[2][2] =
          (bunch_it->prodfc[2][2] / bunch_it->num[2][2]) - (bunch_it->capfc[2] * bunch_it->capfc[2]);
      bunch_it->sigfc[2][3] =
          (bunch_it->prodfc[2][3] / bunch_it->num[2][3]) - (bunch_it->capfc[2] * bunch_it->capfc[3]);
      bunch_it->sigfc[3][0] =
          (bunch_it->prodfc[3][0] / bunch_it->num[3][0]) - (bunch_it->capfc[3] * bunch_it->capfc[0]);
      bunch_it->sigfc[3][1] =
          (bunch_it->prodfc[3][1] / bunch_it->num[3][1]) - (bunch_it->capfc[3] * bunch_it->capfc[1]);
      bunch_it->sigfc[3][2] =
          (bunch_it->prodfc[3][2] / bunch_it->num[3][2]) - (bunch_it->capfc[3] * bunch_it->capfc[2]);
      bunch_it->sigfc[3][3] =
          (bunch_it->prodfc[3][3] / bunch_it->num[3][3]) - (bunch_it->capfc[3] * bunch_it->capfc[3]);

      for (int i = 0; i != 3; i++) {
        CASTORMeans->Fill(bunch_it->cap[i]);
        CASTORWidths->Fill(bunch_it->sig[i][i]);
      }

      //if(bunch_it->detid.subdet() == 1){

      int fillphi = bunch_it->detid.sector();
      //if (bunch_it->detid.depth()==4) fillphi++;

      //    dephist[bunch_it->detid.module()-1]->Fill(bunch_it->detid.ieta(),fillphi,
      //             (bunch_it->cap[0]+bunch_it->cap[1]+bunch_it->cap[2]+bunch_it->cap[3])/4);
      dephist->Fill(bunch_it->detid.module(),
                    fillphi,
                    (bunch_it->cap[0] + bunch_it->cap[1] + bunch_it->cap[2] + bunch_it->cap[3]) / 4);

      const CastorPedestal item(bunch_it->detid,
                                bunch_it->cap[0],
                                bunch_it->cap[1],
                                bunch_it->cap[2],
                                bunch_it->cap[3],
                                bunch_it->sig[0][0],
                                bunch_it->sig[1][1],
                                bunch_it->sig[2][2],
                                bunch_it->sig[3][3]);
      rawPedsItem->addValues(item);
      CastorPedestalWidth widthsp(bunch_it->detid);
      widthsp.setSigma(0, 0, bunch_it->sig[0][0]);
      widthsp.setSigma(0, 1, bunch_it->sig[0][1]);
      widthsp.setSigma(0, 2, bunch_it->sig[0][2]);
      widthsp.setSigma(0, 3, bunch_it->sig[0][3]);
      widthsp.setSigma(1, 0, bunch_it->sig[1][0]);
      widthsp.setSigma(1, 1, bunch_it->sig[1][1]);
      widthsp.setSigma(1, 2, bunch_it->sig[1][2]);
      widthsp.setSigma(1, 3, bunch_it->sig[1][3]);
      widthsp.setSigma(2, 0, bunch_it->sig[2][0]);
      widthsp.setSigma(2, 1, bunch_it->sig[2][1]);
      widthsp.setSigma(2, 2, bunch_it->sig[2][2]);
      widthsp.setSigma(2, 3, bunch_it->sig[2][3]);
      widthsp.setSigma(3, 0, bunch_it->sig[3][0]);
      widthsp.setSigma(3, 1, bunch_it->sig[3][1]);
      widthsp.setSigma(3, 2, bunch_it->sig[3][2]);
      widthsp.setSigma(3, 3, bunch_it->sig[3][3]);
      rawWidthsItem->addValues(widthsp);

      const CastorPedestal itemfc(bunch_it->detid,
                                  bunch_it->capfc[0],
                                  bunch_it->capfc[1],
                                  bunch_it->capfc[2],
                                  bunch_it->capfc[3],
                                  bunch_it->sigfc[0][0],
                                  bunch_it->sigfc[1][1],
                                  bunch_it->sigfc[2][2],
                                  bunch_it->sigfc[3][3]);
      rawPedsItemfc->addValues(itemfc);
      CastorPedestalWidth widthspfc(bunch_it->detid);
      widthspfc.setSigma(0, 0, bunch_it->sigfc[0][0]);
      widthspfc.setSigma(0, 1, bunch_it->sigfc[0][1]);
      widthspfc.setSigma(0, 2, bunch_it->sigfc[0][2]);
      widthspfc.setSigma(0, 3, bunch_it->sigfc[0][3]);
      widthspfc.setSigma(1, 0, bunch_it->sigfc[1][0]);
      widthspfc.setSigma(1, 1, bunch_it->sigfc[1][1]);
      widthspfc.setSigma(1, 2, bunch_it->sigfc[1][2]);
      widthspfc.setSigma(1, 3, bunch_it->sigfc[1][3]);
      widthspfc.setSigma(2, 0, bunch_it->sigfc[2][0]);
      widthspfc.setSigma(2, 1, bunch_it->sigfc[2][1]);
      widthspfc.setSigma(2, 2, bunch_it->sigfc[2][2]);
      widthspfc.setSigma(2, 3, bunch_it->sigfc[2][3]);
      widthspfc.setSigma(3, 0, bunch_it->sigfc[3][0]);
      widthspfc.setSigma(3, 1, bunch_it->sigfc[3][1]);
      widthspfc.setSigma(3, 2, bunch_it->sigfc[3][2]);
      widthspfc.setSigma(3, 3, bunch_it->sigfc[3][3]);
      rawWidthsItemfc->addValues(widthspfc);
    }
  }

  // dump the resulting list of pedestals into a file
  std::ofstream outStream1(pedsADCfilename.c_str());
  CastorDbASCIIIO::dumpObject(outStream1, (*rawPedsItem));
  std::ofstream outStream2(widthsADCfilename.c_str());
  CastorDbASCIIIO::dumpObject(outStream2, (*rawWidthsItem));

  std::ofstream outStream3(pedsfCfilename.c_str());
  CastorDbASCIIIO::dumpObject(outStream3, (*rawPedsItemfc));
  std::ofstream outStream4(widthsfCfilename.c_str());
  CastorDbASCIIIO::dumpObject(outStream4, (*rawWidthsItemfc));

  if (dumpXML) {
    std::ofstream outStream5(XMLfilename.c_str());
    //  CastorCondXML::dumpObject (outStream5, runnum, runnum, runnum, XMLtag, 1, (*rawPedsItem), (*rawWidthsItem));
  }

  dephist->SetDrawOption("colz");
  dephist->GetXaxis()->SetTitle("module");
  dephist->GetYaxis()->SetTitle("sector");

  std::stringstream tempstringout;
  tempstringout << runnum;
  std::string name1 = tempstringout.str() + "_pedplots_1d.png";
  std::string name2 = tempstringout.str() + "_pedplots_2d.png";

  TStyle* theStyle = new TStyle("style", "null");
  theStyle->SetPalette(1, nullptr);
  theStyle->SetCanvasDefH(1200);  //Height of canvas
  theStyle->SetCanvasDefW(1600);  //Width of canvas

  gStyle = theStyle;
  /*
    TCanvas * c1 = new TCanvas("c1","graph",1);
    c1->Divide(2,2);
    c1->cd(1);
    CASTORMeans->Draw();
    c1->SaveAs(name1.c_str());   

    theStyle->SetOptStat("n");
    gStyle = theStyle;

    TCanvas * c2 = new TCanvas("c2","graph",1);
 //   c2->Divide(2,2);
    c2->cd(1);
    dephist->Draw();
    dephist->SetDrawOption("colz");
    //c2->cd(2);
    //dephist[1]->Draw();
    //dephist[1]->SetDrawOption("colz");
    //c2->cd(3);
    //dephist[2]->Draw();
    //dephist[2]->SetDrawOption("colz");
    //c2->cd(4);
    //dephist[3]->Draw();
    //dephist[3]->SetDrawOption("colz");
    c2->SaveAs(name2.c_str());
*/
  /*
  std::cout << "Writing ROOT file... ";
  theFile->Close();
  std::cout << "ROOT file closed.\n";
  */
}

Member Function Documentation

analyze()

void CastorPedestalsAnalysis::analyze ( const edm::Event &  event, const edm::EventSetup &  eventSetup  )

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 343 of file CastorPedestalsAnalysis.cc.

References a, HcalQIESample::adc(), CastorElectronicsMap::allPrecisionId(), Bunches, HcalQIESample::capid(), castorDigitizer_cfi::castor, castorDigiCollectionToken_, CASTORMeans, CASTORWidths, CastorQIECoder::charge(), PixelBaryCentreAnalyzer_cfg::conditions, dephist, MillePedeFileConverter_cfg::e, firsttime, firstTS, compareTotals::fs, edm::EventSetup::getData(), mps_fire::i, CastorDataFrame::id(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmiolumiharvest::j, lastTS, pedsADCfilename, pedsfCfilename, DetId::rawId(), ROOTfilename, runnum, CastorDataFrame::sample(), l1trig_cff::shape, CastorDataFrame::size(), AlCaHLTBitMon_QueryRunRegistry::string, tok_cond_, tok_map_, widthsADCfilename, widthsfCfilename, XMLfilename, and XMLtag.

                                                                                    {
  const edm::Handle<CastorDigiCollection>& castor = e.getHandle(castorDigiCollectionToken_);

  auto conditions = &iSetup.getData(tok_cond_);
  const CastorQIEShape* shape = conditions->getCastorShape();

  if (firsttime) {
    runnum = e.id().run();
    std::string runnum_string;
    std::stringstream tempstringout;
    tempstringout << runnum;
    runnum_string = tempstringout.str();
    ROOTfilename = runnum_string + "-peds_ADC.root";
    pedsADCfilename = runnum_string + "-peds_ADC.txt";
    pedsfCfilename = runnum_string + "-peds_fC.txt";
    widthsADCfilename = runnum_string + "-widths_ADC.txt";
    widthsfCfilename = runnum_string + "-widths_fC.txt";
    XMLfilename = runnum_string + "-peds_ADC_complete.xml";
    XMLtag = "Castor_pedestals_" + runnum_string;

    edm::Service<TFileService> fs;
    fs->mkdir("CASTOR");
    CASTORMeans = fs->make<TH1F>("All Ped Means CASTOR", "All Ped Means CASTOR", 100, 0, 9);
    CASTORWidths = fs->make<TH1F>("All Ped Widths CASTOR", "All Ped Widths CASTOR", 100, 0, 3);

    dephist = fs->make<TH2F>("Pedestals (ADC)", "All Castor", 14, 0., 14.5, 16, .5, 16.5);

    const CastorElectronicsMap* myRefEMap = &iSetup.getData(tok_map_);
    std::vector<HcalGenericDetId> listEMap = myRefEMap->allPrecisionId();
    for (std::vector<HcalGenericDetId>::const_iterator it = listEMap.begin(); it != listEMap.end(); ++it) {
      HcalGenericDetId mygenid(it->rawId());
      if (mygenid.isHcalCastorDetId()) {
        NewPedBunch a;
        HcalCastorDetId chanid(mygenid.rawId());
        a.detid = chanid;
        a.usedflag = false;
        std::string type = "CASTOR";
        for (int i = 0; i != 4; i++) {
          a.cap[i] = 0;
          a.capfc[i] = 0;
          for (int j = 0; j != 4; j++) {
            a.sig[i][j] = 0;
            a.sigfc[i][j] = 0;
            a.prod[i][j] = 0;
            a.prodfc[i][j] = 0;
            a.num[i][j] = 0;
          }
        }
        Bunches.push_back(a);
      }
    }
    firsttime = false;
  }

  std::vector<NewPedBunch>::iterator bunch_it;

  for (CastorDigiCollection::const_iterator j = castor->begin(); j != castor->end(); ++j) {
    const CastorDataFrame digi = (const CastorDataFrame)(*j);
    for (bunch_it = Bunches.begin(); bunch_it != Bunches.end(); ++bunch_it)
      if (bunch_it->detid.rawId() == digi.id().rawId())
        break;
    bunch_it->usedflag = true;
    for (int ts = firstTS; ts != lastTS + 1; ts++) {
      const CastorQIECoder* coder = conditions->getCastorCoder(digi.id().rawId());
      bunch_it->num[digi.sample(ts).capid()][digi.sample(ts).capid()] += 1;
      bunch_it->cap[digi.sample(ts).capid()] += digi.sample(ts).adc();
      double charge1 = coder->charge(*shape, digi.sample(ts).adc(), digi.sample(ts).capid());
      bunch_it->capfc[digi.sample(ts).capid()] += charge1;
      bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts).capid()] +=
          (digi.sample(ts).adc() * digi.sample(ts).adc());
      bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts).capid()] += charge1 * charge1;
      if ((ts + 1 < digi.size()) && (ts + 1 < lastTS)) {
        bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts + 1).capid()] +=
            digi.sample(ts).adc() * digi.sample(ts + 1).adc();
        double charge2 = coder->charge(*shape, digi.sample(ts + 1).adc(), digi.sample(ts + 1).capid());
        bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts + 1).capid()] += charge1 * charge2;
        bunch_it->num[digi.sample(ts).capid()][digi.sample(ts + 1).capid()] += 1;
      }
      if ((ts + 2 < digi.size()) && (ts + 2 < lastTS)) {
        bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts + 2).capid()] +=
            digi.sample(ts).adc() * digi.sample(ts + 2).adc();
        double charge2 = coder->charge(*shape, digi.sample(ts + 2).adc(), digi.sample(ts + 2).capid());
        bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts + 2).capid()] += charge1 * charge2;
        bunch_it->num[digi.sample(ts).capid()][digi.sample(ts + 2).capid()] += 1;
      }
      if ((ts + 3 < digi.size()) && (ts + 3 < lastTS)) {
        bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts + 3).capid()] +=
            digi.sample(ts).adc() * digi.sample(ts + 3).adc();
        double charge2 = coder->charge(*shape, digi.sample(ts + 3).adc(), digi.sample(ts + 3).capid());
        bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts + 3).capid()] += charge1 * charge2;
        bunch_it->num[digi.sample(ts).capid()][digi.sample(ts + 3).capid()] += 1;
      }
    }
  }

  //this is the last brace
}

Member Data Documentation

Bunches

std::vector<NewPedBunch> CastorPedestalsAnalysis::Bunches

private

Definition at line 85 of file CastorPedestalsAnalysis.cc.

Referenced by analyze(), and ~CastorPedestalsAnalysis().

castorDigiCollectionToken_

const edm::EDGetTokenT<CastorDigiCollection> CastorPedestalsAnalysis::castorDigiCollectionToken_

private

Definition at line 86 of file CastorPedestalsAnalysis.cc.

Referenced by analyze().

CASTORMeans

TH1F* CastorPedestalsAnalysis::CASTORMeans

private

Definition at line 106 of file CastorPedestalsAnalysis.cc.

Referenced by analyze(), and ~CastorPedestalsAnalysis().

CASTORWidths

TH1F* CastorPedestalsAnalysis::CASTORWidths

private

Definition at line 107 of file CastorPedestalsAnalysis.cc.

Referenced by analyze(), and ~CastorPedestalsAnalysis().

dephist

TH2F* CastorPedestalsAnalysis::dephist

private

Definition at line 110 of file CastorPedestalsAnalysis.cc.

Referenced by analyze(), and ~CastorPedestalsAnalysis().

dumpXML

const bool CastorPedestalsAnalysis::dumpXML

private

Definition at line 91 of file CastorPedestalsAnalysis.cc.

Referenced by ~CastorPedestalsAnalysis().

firsttime

bool CastorPedestalsAnalysis::firsttime

private

Definition at line 95 of file CastorPedestalsAnalysis.cc.

Referenced by analyze(), and CastorPedestalsAnalysis().

firstTS

const int CastorPedestalsAnalysis::firstTS

private

Definition at line 93 of file CastorPedestalsAnalysis.cc.

Referenced by analyze().

hiSaveFlag

const bool CastorPedestalsAnalysis::hiSaveFlag

private

Definition at line 90 of file CastorPedestalsAnalysis.cc.

lastTS

const int CastorPedestalsAnalysis::lastTS

private

Definition at line 94 of file CastorPedestalsAnalysis.cc.

Referenced by analyze().

pedsADCfilename

std::string CastorPedestalsAnalysis::pedsADCfilename

private

Definition at line 98 of file CastorPedestalsAnalysis.cc.

Referenced by analyze(), and ~CastorPedestalsAnalysis().

pedsfCfilename

std::string CastorPedestalsAnalysis::pedsfCfilename

private

Definition at line 99 of file CastorPedestalsAnalysis.cc.

Referenced by analyze(), and ~CastorPedestalsAnalysis().

ROOTfilename

std::string CastorPedestalsAnalysis::ROOTfilename

private

Definition at line 97 of file CastorPedestalsAnalysis.cc.

Referenced by analyze().

runnum

int CastorPedestalsAnalysis::runnum

private

Definition at line 96 of file CastorPedestalsAnalysis.cc.

Referenced by analyze(), and ~CastorPedestalsAnalysis().

tok_cond_

const edm::ESGetToken<CastorDbService, CastorDbRecord> CastorPedestalsAnalysis::tok_cond_

private

Definition at line 87 of file CastorPedestalsAnalysis.cc.

Referenced by analyze().

tok_map_

const edm::ESGetToken<CastorElectronicsMap, CastorElectronicsMapRcd> CastorPedestalsAnalysis::tok_map_

private

Definition at line 88 of file CastorPedestalsAnalysis.cc.

Referenced by analyze().

verboseflag

const bool CastorPedestalsAnalysis::verboseflag

private

Definition at line 92 of file CastorPedestalsAnalysis.cc.

Referenced by ~CastorPedestalsAnalysis().

widthsADCfilename

std::string CastorPedestalsAnalysis::widthsADCfilename

private

Definition at line 100 of file CastorPedestalsAnalysis.cc.

Referenced by analyze(), and ~CastorPedestalsAnalysis().

widthsfCfilename

std::string CastorPedestalsAnalysis::widthsfCfilename

private

Definition at line 101 of file CastorPedestalsAnalysis.cc.

Referenced by analyze(), and ~CastorPedestalsAnalysis().

XMLfilename

std::string CastorPedestalsAnalysis::XMLfilename

private

Definition at line 102 of file CastorPedestalsAnalysis.cc.

Referenced by analyze(), and ~CastorPedestalsAnalysis().

XMLtag

std::string CastorPedestalsAnalysis::XMLtag

private

Definition at line 103 of file CastorPedestalsAnalysis.cc.

Referenced by analyze().

ZSfilename

std::string CastorPedestalsAnalysis::ZSfilename

private

Definition at line 104 of file CastorPedestalsAnalysis.cc.