Inheritance diagram for WriteCTPPSPixGainCalibrations:

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

	~WriteCTPPSPixGainCalibrations () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () 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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

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 ( void )

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.

                                            {
   getHistos();
   fillDB();
 }

References fillDB(), and getHistos().

◆ fillDB()

void WriteCTPPSPixGainCalibrations::fillDB ( )

private

Definition at line 190 of file WriteCTPPSPixGainCalibrations.cc.

                                            {
   CTPPSPixelGainCalibrations* gainCalibsTest = new CTPPSPixelGainCalibrations();
   CTPPSPixelGainCalibrations* gainCalibsTest1 = new CTPPSPixelGainCalibrations();
  
   //  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->writeOne(gainCalibsTest, mydbservice->currentTime(), m_record);
 }

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

Referenced by endJob().

◆ fillDescriptions()

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

static

Definition at line 97 of file WriteCTPPSPixGainCalibrations.cc.

                                                                                                {
   //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);
 }

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

◆ 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.

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

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

Referenced by fillDB().

◆ getHistos()

void WriteCTPPSPixGainCalibrations::getHistos ( )

private

Definition at line 105 of file WriteCTPPSPixGainCalibrations.cc.

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

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

Referenced by endJob().

◆ setDummyFullPlane()

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

private

Definition at line 242 of file WriteCTPPSPixGainCalibrations.cc.

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