Inheritance diagram for WriteCTPPSPixGainCalibrations:

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

	~WriteCTPPSPixGainCalibrations () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	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

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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

void	endJob () override

void	fillDB ()

void	getGainsPedsFromHistos (uint32_t detid, int rocId, int index, std::vector< double > &peds, std::vector< double > &gains, std::map< int, int > &myindxmap, int nrocs)

void	getHistos ()

void	setDummyFullPlane (std::vector< float > &peds, std::vector< float > &gains, int npixplane)

Private Attributes
std::map< uint32_t, std::vector< std::string > >	detidHistoNameMap

std::map< uint32_t, std::vector< int > >	detidROCsPresent

double	gainhigh

double	gainlow

std::string	m_inputHistosFileName

TFile *	m_inputRootFile

std::string	m_record

bool	m_usedummy

int	npfitmin

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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 26 of file WriteCTPPSPixGainCalibrations.cc.

Constructor & Destructor Documentation

◆ WriteCTPPSPixGainCalibrations()

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

explicit

Definition at line 65 of file WriteCTPPSPixGainCalibrations.cc.

     : m_record(iConfig.getUntrackedParameter<std::string>("record", "CTPPSPixelGainCalibrationsRcd")),
       m_inputHistosFileName(iConfig.getUntrackedParameter<std::string>("inputrootfile", "inputfile.root")),
       m_usedummy(iConfig.getUntrackedParameter<bool>("useDummyValues", true)),
       npfitmin(iConfig.getUntrackedParameter<int>("minimumNpfit", 3)),
       gainlow(iConfig.getUntrackedParameter<double>("gainLowLimit", 0.0)),
       gainhigh(iConfig.getUntrackedParameter<double>("gainHighLimit", 100.0)) {}

◆ ~WriteCTPPSPixGainCalibrations()

WriteCTPPSPixGainCalibrations::~WriteCTPPSPixGainCalibrations ( )

override

Definition at line 73 of file WriteCTPPSPixGainCalibrations.cc.

                                                               {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 83 of file WriteCTPPSPixGainCalibrations.cc.

                                                                                                {
   //   using namespace edm;
 }

◆ beginJob()

void WriteCTPPSPixGainCalibrations::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 88 of file WriteCTPPSPixGainCalibrations.cc.

88 {}

◆ endJob()

void WriteCTPPSPixGainCalibrations::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 91 of file WriteCTPPSPixGainCalibrations.cc.

References fillDB(), and getHistos().

                                            {
   getHistos();
   fillDB();
 }

◆ fillDB()

void WriteCTPPSPixGainCalibrations::fillDB ( )

private

Definition at line 190 of file WriteCTPPSPixGainCalibrations.cc.

References cond::service::PoolDBOutputService::currentTime(), detidROCsPresent, getGainsPedsFromHistos(), mps_fire::i, edm::Service< T >::isAvailable(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmdumpme::k, m_record, genParticles_cff::map, DeadROCCounter::nrocs, CTPPSPixelGainCalibration::putData(), PixelMapPlotter::rocs, setDummyFullPlane(), CTPPSPixelGainCalibrations::setGainCalibration(), and cond::service::PoolDBOutputService::writeOneIOV().

Referenced by endJob().

                                            {
   CTPPSPixelGainCalibrations gainCalibsTest;
   CTPPSPixelGainCalibrations gainCalibsTest1;
 
   //  std::cout<<"Here! "<<std::endl;
 
   for (std::map<uint32_t, std::vector<int> >::iterator it = detidROCsPresent.begin(); it != detidROCsPresent.end();
        ++it) {
     uint32_t tempdetid = it->first;
     std::vector<int> rocs = it->second;
     unsigned int nrocs = rocs.size();
     std::map<int, int> mapIPixIndx;
 
     std::vector<double> gainsFromHistos;
     std::vector<double> pedsFromHistos;
 
     CTPPSPixelGainCalibration tempPGCalib;
 
     for (unsigned int i = 0; i < nrocs; i++) {
       getGainsPedsFromHistos(tempdetid, i, rocs[i], pedsFromHistos, gainsFromHistos, mapIPixIndx, nrocs);
     }
 
     std::vector<float> orderedGains;
     std::vector<float> orderedPeds;
     for (unsigned int k = 0; k < nrocs * 52 * 80; k++) {
       int indx = mapIPixIndx[k];
       float tmpped = pedsFromHistos[indx];
       float tmpgain = gainsFromHistos[indx];
       orderedGains.push_back(tmpgain);
       orderedPeds.push_back(tmpped);
       tempPGCalib.putData(k, tmpped, tmpgain);
     }
 
     if (nrocs == 0) {
       edm::LogWarning("CTPPSPixGainsCalibrationWriter") << " plane with detID =" << tempdetid << " is empty";
       setDummyFullPlane(orderedPeds, orderedGains, 6 * 52 * 80);
     }
 
     gainCalibsTest.setGainCalibration(tempdetid, orderedPeds, orderedGains);
     //   std::cout << "Here detid = "<<tempdetid <<std::endl;
     gainCalibsTest1.setGainCalibration(tempdetid, tempPGCalib);
     //   std::cout << "Here again"<<std::endl;
   }
   //  std::cout<<" Here 3!"<<std::endl;
   edm::Service<cond::service::PoolDBOutputService> mydbservice;
   if (!mydbservice.isAvailable()) {
     edm::LogError("CTPPSPixGainsCalibrationWriter") << "Db Service Unavailable";
     return;
   }
   mydbservice->writeOneIOV(gainCalibsTest, mydbservice->currentTime(), m_record);
 }

◆ fillDescriptions()

void WriteCTPPSPixGainCalibrations::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 97 of file WriteCTPPSPixGainCalibrations.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

                                                                                                {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

◆ getGainsPedsFromHistos()

void WriteCTPPSPixGainCalibrations::getGainsPedsFromHistos	(	uint32_t	detid,
		int	rocId,
		int	index,
		std::vector< double > &	peds,
		std::vector< double > &	gains,
		std::map< int, int > &	myindxmap,
		int	nrocs
	)

private

Definition at line 251 of file WriteCTPPSPixGainCalibrations.cc.

References detidHistoNameMap, gainhigh, ecalph2::gains, m_inputRootFile, m_usedummy, hgcalPlots::ncols, npfitmin, DeadROCCounter::nrocs, AlCaHLTBitMon_QueryRunRegistry::string, and CTPPSPixelIndices::transformToModule().

Referenced by fillDB().

                                                                       {
   CTPPSPixelIndices modulepixels(52 * nrocs / 2, 160);
 
   std::string tmpslopename = detidHistoNameMap[detid][ROCindex] + "_Slope2D";
   std::string tmpitcpname = detidHistoNameMap[detid][ROCindex] + "_Intercept2D";
   std::string tmpchi2name = detidHistoNameMap[detid][ROCindex] + "_Chisquare2D";
   std::string tmpnpfitsname = detidHistoNameMap[detid][ROCindex] + "_Npfits2D";
   TH2D* tempslope = (TH2D*)m_inputRootFile->Get(tmpslopename.c_str());
   TH2D* tempintrcpt = (TH2D*)m_inputRootFile->Get(tmpitcpname.c_str());
   // TH2D * tempchi2    = (TH2D*) m_inputRootFile->Get(tmpchi2name.c_str());
   TH2D* tempnpfit = (TH2D*)m_inputRootFile->Get(tmpnpfitsname.c_str());
   int ncols = tempslope->GetNbinsX();
   int nrows = tempslope->GetNbinsY();
   if (nrows != 80 || ncols != 52)
     edm::LogWarning("CTPPSPixGainsCalibrationWriter")
         << "Something wrong ncols = " << ncols << " and nrows = " << nrows;
 
   for (int jrow = 0; jrow < nrows;
        ++jrow)  // when scanning through the 2d histo make sure to avoid underflow bin i or j ==0
     for (int icol = 0; icol < ncols; ++icol) {
       double tmpslp = tempslope->GetBinContent(icol + 1, jrow + 1);
       double tmpgain = (tmpslp == 0.0) ? 0.0 : 1.0 / tmpslp;
       double tmpped = tempintrcpt->GetBinContent(icol + 1, jrow + 1);
       // check for noisy/badly calibrated pixels
       int tmpnpfit = tempnpfit->GetBinContent(icol + 1, jrow + 1);
       //double tmpchi2  = tempchi2 -> GetBinContent(icol+1,jrow+1);
       if (tmpnpfit < npfitmin || tmpgain < gainlow || tmpgain > gainhigh) {
         //  std::cout << " *** Badly calibrated pixel, NPoints = "<<tmpnpfit << " setting dummy values gain = 0.5 and  ped =20.0 ***" <<std::endl;
         //  std::cout << " **** bad Pixel column icol = "<<icol <<" and jrow = "<<jrow <<" Name= "<< tmpslopename <<std::endl;
         if (m_usedummy) {
           tmpgain = 1.0 / 2.0;
           tmpped = 20.0;
         }
         // else  leave as is and set noisy in mask?
       }
 
       gains.push_back(tmpgain);
       peds.push_back(tmpped);
       int modCol = -1;
       int modRow = -1;
       modulepixels.transformToModule(icol, jrow, rocId, modCol, modRow);
       int indx = gains.size() - 1;
       int pixIndx = modCol + modRow * (52 * nrocs / 2);
       mymap[pixIndx] = indx;
     }
 }

◆ getHistos()

void WriteCTPPSPixGainCalibrations::getHistos ( )

private

Definition at line 105 of file WriteCTPPSPixGainCalibrations.cc.

References protons_cff::arm, detidHistoNameMap, detidROCsPresent, mps_fire::i, m_inputHistosFileName, m_inputRootFile, DetId::rawId(), PixelMapPlotter::roc, nano_mu_digi_cff::sector, relativeConstraints::station, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by endJob().

                                               {
   //  std::cout <<"Parsing file " <<m_inputHistosFileName << std::endl;
   m_inputRootFile = new TFile(m_inputHistosFileName.c_str());
   m_inputRootFile->cd();
 
   int sector[2] = {45, 56};  // arm 0 is sector 45 and arm 1 is sector 56
   int nsec = 2;
   int station[2] = {0, 2};  // for each arm
   int nst = 2;
   //int pot[6]={3}; // index of the pot within the 6 pot configuration (vertical top or bottom and horizontal)
   int npt = 6;
 
   for (int i = 0; i < nsec; i++)
     for (int st = 0; st < nst; st++)
       for (int pot = 0; pot < npt; pot++) {
         int arm = i;
 
         //Check which pots present
         char temppathrp[100];
         sprintf(temppathrp, "CTPPS/CTPPS_SEC%d/CTPPS_SEC%d_RP%d%d%d", sector[i], sector[i], arm, station[st], pot);
         if (!m_inputRootFile->Get(temppathrp))
           continue;
 
         for (int plane = 0; plane < 6; plane++) {
           //Check which planes present
           char temppathplane[100];
           sprintf(temppathplane,
                   "CTPPS/CTPPS_SEC%d/CTPPS_SEC%d_RP%d%d%d/CTPPS_SEC%d_RP%d%d%d_PLN%d",
                   sector[i],
                   sector[i],
                   arm,
                   station[st],
                   pot,
                   sector[i],
                   arm,
                   station[st],
                   pot,
                   plane);
 
           // do not skip the missing plane, instead put dummy values
           //      if(!m_inputRootFile->Get(temppathplane)) continue;
 
           CTPPSPixelDetId mytempid(arm, station[st], pot, plane);
           std::vector<std::string> histnamevec;
           std::vector<int> listrocs;
           for (int roc = 0; roc < 6; roc++) {
             char temppathhistos[200];
 
             sprintf(
                 temppathhistos,
                 "CTPPS/CTPPS_SEC%d/CTPPS_SEC%d_RP%d%d%d/CTPPS_SEC%d_RP%d%d%d_PLN%d/CTPPS_SEC%d_RP%d%d%d_PLN%d_ROC%d",
                 sector[i],
                 sector[i],
                 arm,
                 station[st],
                 pot,
                 sector[i],
                 arm,
                 station[st],
                 pot,
                 plane,
                 sector[i],
                 arm,
                 station[st],
                 pot,
                 plane,
                 roc);
 
             std::string pathhistos(temppathhistos);
             std::string pathslope = pathhistos + "_Slope2D";
             std::string pathintercept = pathhistos + "_Intercept2D";
             if (m_inputRootFile->Get(pathslope.c_str()) && m_inputRootFile->Get(pathintercept.c_str())) {
               histnamevec.push_back(pathhistos);
               listrocs.push_back(roc);
             }
           }
           detidHistoNameMap[mytempid.rawId()] = histnamevec;
           detidROCsPresent[mytempid.rawId()] = listrocs;
           edm::LogInfo("CTPPSPixGainsCalibrationWriter")
               << "Raw DetId = " << mytempid.rawId() << " Arm = " << arm << " Sector = " << sector[arm]
               << " Station = " << station[st] << " Pot = " << pot << " Plane = " << plane;
         }
       }
 }

◆ setDummyFullPlane()

void WriteCTPPSPixGainCalibrations::setDummyFullPlane	(	std::vector< float > &	peds,
		std::vector< float > &	gains,
		int	npixplane
	)

private

Definition at line 242 of file WriteCTPPSPixGainCalibrations.cc.

References ecalph2::gains, and mps_fire::i.

Referenced by fillDB().

                                                                      {
   for (int i = 0; i < npixplane; ++i) {
     peds.push_back(20.);
     gains.push_back(0.5);
   }
 }