#include <StoreEcalCondition.h>

Inheritance diagram for StoreEcalCondition:

Public Member Functions
void	analyze (const edm::Event &evt, const edm::EventSetup &evtSetup) override

int	convertFromConstructionSMToSlot (int, int)

void	endJob () override

EcalADCToGeVConstant *	readEcalADCToGeVConstantFromFile (const char *)

EcalChannelStatus *	readEcalChannelStatusFromFile (const char *)

EcalGainRatios *	readEcalGainRatiosFromFile (const char *)

EcalIntercalibConstants *	readEcalIntercalibConstantsFromFile (const char , const char )

EcalIntercalibConstantsMC *	readEcalIntercalibConstantsMCFromFile (const char , const char )

EcalPFRecHitThresholds *	readEcalPFRecHitThresholdsFromFile (const char , const char )

EcalTBWeights *	readEcalTBWeightsFromFile (const char *)

EcalWeightXtalGroups *	readEcalWeightXtalGroupsFromFile (const char *)

	StoreEcalCondition (const edm::ParameterSet &iConfig)

void	writeToLogFile (std::string, std::string, unsigned long long)

void	writeToLogFileResults (char *)

	~StoreEcalCondition () override

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	fillHeader (char *)

std::string	to_string (char value[])

Private Attributes
std::vector< std::string >	inpFileName_

std::vector< std::string >	inpFileNameEE_

std::string	logfile_

std::vector< std::string >	objectName_

std::string	prog_name_

std::vector< unsigned long long >	since_

int	sm_constr_

int	sm_slot_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

Detailed Description

Definition at line 39 of file StoreEcalCondition.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 21 of file StoreEcalCondition.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), inpFileName_, inpFileNameEE_, logfile_, objectName_, prog_name_, since_, sm_constr_, sm_slot_, AlCaHLTBitMon_QueryRunRegistry::string, and DBConfiguration_cff::toPut.

                                                                      {
 
   prog_name_ = "StoreEcalCondition";
 
   logfile_ = iConfig.getParameter< std::string >("logfile");
   sm_slot_ = iConfig.getUntrackedParameter< unsigned int >("smSlot", 1);
 
   typedef std::vector< edm::ParameterSet > Parameters;
   Parameters toPut=iConfig.getParameter<Parameters>("toPut");
   for(Parameters::iterator itToPut = toPut.begin(); itToPut != toPut.end(); ++itToPut) 
     {
       inpFileName_.push_back(itToPut->getUntrackedParameter<std::string>("inputFile"));
       inpFileNameEE_.push_back(itToPut->getUntrackedParameter<std::string>("inputFileEE"));
       objectName_.push_back(itToPut->getUntrackedParameter<std::string>("conditionType"));
       since_.push_back(itToPut->getUntrackedParameter<unsigned int>("since"));
     }
 
   sm_constr_ = -1;
 }

StoreEcalCondition::~StoreEcalCondition ( )

override

Definition at line 149 of file StoreEcalCondition.cc.

                                         {
 
 }

Member Function Documentation

void StoreEcalCondition::analyze	(	const edm::Event &	evt,
		const edm::EventSetup &	evtSetup
	)

override

Definition at line 154 of file StoreEcalCondition.cc.

                                                                                     {
 //-------------------------------------------------------------
 
 }

int StoreEcalCondition::convertFromConstructionSMToSlot	(	int	sm_constr,
		int	sm_slot
	)

Definition at line 929 of file StoreEcalCondition.cc.

References mps_fire::result.

Referenced by readEcalIntercalibConstantsFromFile(), and readEcalIntercalibConstantsMCFromFile().

                                                                                 {
   // input either cosntruction number or slot number and returns the other
   // the slots are numbered first EB+ slot 1 ...18 then EB- 1... 18 
   // the slots start at 1 and the SM start at 0
 //-------------------------------------------------------------
   int slot_to_constr[37]={-1,12,17,10,1,8,4,27,20,23,25,6,34,35,15,18,30,21,9 
               ,24,22,13,31,26,16,2,11,5,0,29,28,14,33,32,3,7,19};
   int constr_to_slot[36]={28,4,25,34,6,27,11,35,5,18,3,26,1,21,31,14,24,2,15,
               36,8,17,20,9,19,10,23,7,30,29,16,22,33,32,12,13  };
   
   int result=0;
   if(sm_constr!=-1) {
     result=constr_to_slot[sm_constr];
   } else if(sm_slot!=-1) {
     result=slot_to_constr[sm_slot];
   }
   return result;
 }

void StoreEcalCondition::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 41 of file StoreEcalCondition.cc.

References cond::service::PoolDBOutputService::appendSinceTime(), cond::service::PoolDBOutputService::beginOfTime(), gather_cfg::cout, cond::service::PoolDBOutputService::createNewIOV(), cond::service::PoolDBOutputService::endOfTime(), mps_fire::i, inpFileName_, inpFileNameEE_, edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), python.rootplot.argparse::message, objectName_, readEcalADCToGeVConstantFromFile(), readEcalChannelStatusFromFile(), readEcalGainRatiosFromFile(), readEcalIntercalibConstantsFromFile(), readEcalIntercalibConstantsMCFromFile(), readEcalPFRecHitThresholdsFromFile(), readEcalTBWeightsFromFile(), readEcalWeightXtalGroupsFromFile(), since_, AlCaHLTBitMon_QueryRunRegistry::string, and writeToLogFileResults().

                                 {
   
   edm::Service<cond::service::PoolDBOutputService> mydbservice;
   if( !mydbservice.isAvailable() ){
     edm::LogError("StoreEcalCondition")<<"PoolDBOutputService is unavailable"<<std::endl;
     return;
   }
 
   bool toAppend=false;
   // copy a string to the char *
   std::string message = "finished OK\n";
   size_t messageSize = message.size() + 1;
   char * messChar = new char [messageSize];
   strncpy(messChar, message.c_str(), messageSize);
       
   for (unsigned int i=0;i<objectName_.size();i++) {
       cond::Time_t newTime;
       
       if ( mydbservice->isNewTagRequest( objectName_[i]+std::string("Rcd") ) ) {
     // This is the first object for this tag.
     // Append mode should be off.
     // newTime is the end of this new objects IOV.
     newTime = mydbservice->beginOfTime();
       } else {
     // There should already be an object in the DB for this tag.
     // Append IOV mode should be on.
     // newTime is the beginning of this new objects IOV.
     toAppend=true;
     newTime = (cond::Time_t)since_[i];
       }
       edm::LogInfo("StoreEcalCondition") << "Reading " << objectName_[i] 
                      << " from file and writing to DB with newTime " << newTime << endl;
       std::cout << "Reading " << objectName_[i] 
                      << " from file and writing to DB with newTime " << newTime << endl;
       
       if (objectName_[i] == "EcalWeightXtalGroups") {
     EcalWeightXtalGroups* mycali = readEcalWeightXtalGroupsFromFile(inpFileName_[i].c_str());
     if(!toAppend){
       mydbservice->createNewIOV<EcalWeightXtalGroups>(mycali,newTime,mydbservice->endOfTime(),"EcalWeightXtalGroupsRcd");
     }else{
       mydbservice->appendSinceTime<EcalWeightXtalGroups>(mycali,newTime,"EcalWeightXtalGroupsRcd");
     }
       }else if (objectName_[i]  =="EcalTBWeights") {
     EcalTBWeights* mycali=readEcalTBWeightsFromFile(inpFileName_[i].c_str());
     if(!toAppend){
       mydbservice->createNewIOV<EcalTBWeights>(mycali,newTime,mydbservice->endOfTime(),"EcalTBWeightsRcd");
     }else{
       mydbservice->appendSinceTime<EcalTBWeights>(mycali,newTime,"EcalTBWeightsRcd");
     }
       } else if (objectName_[i]  == "EcalADCToGeVConstant") {
     EcalADCToGeVConstant* mycali=readEcalADCToGeVConstantFromFile(inpFileName_[i].c_str());
     if(!toAppend){
       mydbservice->createNewIOV<EcalADCToGeVConstant>(mycali,newTime,mydbservice->endOfTime(),"EcalADCToGeVConstantRcd");
     }else{
       mydbservice->appendSinceTime<EcalADCToGeVConstant>(mycali,newTime,"EcalADCToGeVConstantRcd");
     }
       } else if (objectName_[i]  ==  "EcalIntercalibConstants") {
     EcalIntercalibConstants* mycali=readEcalIntercalibConstantsFromFile(inpFileName_[i].c_str(),inpFileNameEE_[i].c_str());
     if(!toAppend){
       mydbservice->createNewIOV<EcalIntercalibConstants>(mycali,newTime,mydbservice->endOfTime(),"EcalIntercalibConstantsRcd");
     }else{
       mydbservice->appendSinceTime<EcalIntercalibConstants>(mycali,newTime,"EcalIntercalibConstantsRcd");
     }
       } else if (objectName_[i]  ==  "EcalPFRecHitThresholds") {
     EcalPFRecHitThresholds* mycali=readEcalPFRecHitThresholdsFromFile(inpFileName_[i].c_str(),inpFileNameEE_[i].c_str());
     if(!toAppend){
       mydbservice->createNewIOV<EcalPFRecHitThresholds>(mycali,newTime,mydbservice->endOfTime(),"EcalPFRecHitThresholdsRcd");
     }else{
       mydbservice->appendSinceTime<EcalPFRecHitThresholds>(mycali,newTime,"EcalPFRecHitThresholdsRcd");
     }
       } else if (objectName_[i]  ==  "EcalIntercalibConstantsMC") {
     EcalIntercalibConstantsMC* mycali=readEcalIntercalibConstantsMCFromFile(inpFileName_[i].c_str(),inpFileNameEE_[i].c_str());
     if(!toAppend){
       mydbservice->createNewIOV<EcalIntercalibConstantsMC>(mycali,newTime,mydbservice->endOfTime(),"EcalIntercalibConstantsMCRcd");
     }else{
       mydbservice->appendSinceTime<EcalIntercalibConstantsMC>(mycali,newTime,"EcalIntercalibConstantsMCRcd");
     }
       } else if (objectName_[i]  ==  "EcalGainRatios") {
     EcalGainRatios* mycali=readEcalGainRatiosFromFile(inpFileName_[i].c_str());
     if(!toAppend){
       mydbservice->createNewIOV<EcalGainRatios>(mycali,newTime,mydbservice->endOfTime(),"EcalGainRatiosRcd");
     }else{
       mydbservice->appendSinceTime<EcalGainRatios>(mycali,newTime,"EcalGainRatiosRcd");
     }
       } else if (objectName_[i]  ==  "EcalChannelStatus") 
     {
       EcalChannelStatus* mycali=readEcalChannelStatusFromFile(inpFileName_[i].c_str());
       if(!toAppend){
         mydbservice->createNewIOV<EcalChannelStatus>(mycali,newTime,mydbservice->endOfTime(),"EcalChannelStatusRcd");
       }else{
         mydbservice->appendSinceTime<EcalChannelStatus>(mycali,newTime,"EcalChannelStatusRcd");
       } 
     } else {
     edm::LogError("StoreEcalCondition")<< "Object " << objectName_[i]  << " is not supported by this program." << endl;
       }
       
 
       //      writeToLogFile(objectName_[i], inpFileName_[i], since_[i]);
       //writeToLogFileResults("finished OK\n");
       writeToLogFileResults(messChar);
    
       edm::LogInfo("StoreEcalCondition") << "Finished endJob" << endl;
   }
   
   delete [] messChar;
 }

void StoreEcalCondition::fillHeader ( char * header )

private

Definition at line 201 of file StoreEcalCondition.cc.

References ntuplemaker::time, and dataset::user.

Referenced by writeToLogFile(), and writeToLogFileResults().

 {
   time_t rawtime;
   struct tm * timeinfo;
   time ( &rawtime );
   timeinfo = localtime ( &rawtime );
   char user[50];
   sprintf(user,"%s",getlogin());
   sprintf(header,"%s %s:",asctime(timeinfo),user);
 }

EcalADCToGeVConstant * StoreEcalCondition::readEcalADCToGeVConstantFromFile ( const char * inputFile )

Definition at line 458 of file StoreEcalCondition.cc.

References analyzePatCleaning_cfg::inputFile, geometryCSVtoXML::line, LogDebug, sm_constr_, harvestTrackValidationPlots::str, and to_string().

Referenced by endJob().

                                                                           {
 //-------------------------------------------------------------
   
     
     FILE *inpFile; // input file
     inpFile = fopen(inputFile,"r");
     if(!inpFile) {
       edm::LogError("StoreEcalCondition")<<"*** Can not open file: "<<inputFile;
       return nullptr;
     }
 
     char line[256];
     
     std::ostringstream str;
 
     fgets(line,255,inpFile);
     int sm_number=atoi(line);
     str << "sm= " << sm_number << endl ;  
 
     fgets(line,255,inpFile);
     //int nevents=atoi(line); // not necessary here just for online conddb
     
     fgets(line,255,inpFile);
     string gen_tag=to_string(line);
     str << "gen tag " << gen_tag << endl ;  // should I use this? 
 
     fgets(line,255,inpFile);
     string cali_method=to_string(line);
     str << "cali method " << cali_method << endl ; // not important 
 
     fgets(line,255,inpFile);
     string cali_version=to_string(line);
     str << "cali version " << cali_version << endl ; // not important 
 
     fgets(line,255,inpFile);
     string cali_type=to_string(line);
     str << "cali type " << cali_type << endl ; // not important
 
     edm::LogInfo("StoreEcalCondition") << "ADCToGeV file " << str.str() ;
 
     fgets(line,255,inpFile);
     float adc_to_gev=0;
     sscanf(line, "%f", &adc_to_gev );
     LogDebug("StoreEcalCondition") <<" calib="<< adc_to_gev ;
     fgets(line,255,inpFile);
     float adc_to_gev_ee=0;
     sscanf(line, "%f", &adc_to_gev_ee );
     LogDebug("StoreEcalCondition") <<" calib="<< adc_to_gev_ee ;
     
     fclose(inpFile);           // close inp. file
 
     sm_constr_ = sm_number;
 
     // barrel and endcaps the same 
     EcalADCToGeVConstant* agc = new EcalADCToGeVConstant(adc_to_gev,adc_to_gev_ee );
     edm::LogInfo("StoreEcalCondition") << "ADCtoGeV scale written into the DB";
     return agc;
 }

EcalChannelStatus * StoreEcalCondition::readEcalChannelStatusFromFile ( const char * inputFile )

Definition at line 1109 of file StoreEcalCondition.cc.

References comment, gather_cfg::cout, Exception, ecalpyutils::hashedIndex(), compare_using_db::ifile, cuy::ii, ecalpyutils::ism(), EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, geometryCSVtoXML::line, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EcalCondObjectContainer< T >::setValue(), EBDetId::SMCRYSTALMODE, mps_update::status, AlCaHLTBitMon_QueryRunRegistry::string, EBDetId::unhashIndex(), EEDetId::unhashIndex(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by endJob().

 {
         EcalChannelStatus* status = new EcalChannelStatus();
         // barrel
         for(int ieta=-EBDetId::MAX_IETA; ieta<=EBDetId::MAX_IETA; ++ieta) 
       {
         if(ieta==0) continue;
         for(int iphi=EBDetId::MIN_IPHI; iphi<=EBDetId::MAX_IPHI; ++iphi) {
           if (EBDetId::validDetId(ieta,iphi)) {
         EBDetId ebid(ieta,iphi);
         status->setValue( ebid, 0 );
           }
         }
       }
         // endcap
         for(int iX=EEDetId::IX_MIN; iX<=EEDetId::IX_MAX ;++iX) 
       {
         for(int iY=EEDetId::IY_MIN; iY<=EEDetId::IY_MAX; ++iY) 
           {
         // make an EEDetId since we need EEDetId::rawId() to be used as the key for the pedestals
         if (EEDetId::validDetId(iX,iY,1)) 
           {
           EEDetId eedetidpos(iX,iY,1);
           status->setValue( eedetidpos, 0 );
           }
         if (EEDetId::validDetId(iX,iY,-1)) 
           {
             EEDetId eedetidneg(iX,iY,-1);
             status->setValue( eedetidneg, 0 );
           }
           }
       }
 
     
     //        edm::LogInfo("EcalTrivialConditionRetriever") << "Reading channel status from file " << edm::FileInPath(channelStatusFile_).fullPath().c_str() ;
         std::cout << "Reading channel status from file " << inputFile << std::endl;
         FILE *ifile = fopen( inputFile ,"r" );
         if ( !ifile ) 
       throw cms::Exception ("Cannot open ECAL channel status file") ;
 
     char line[256];
 
     fgets(line,255,ifile);
     std::string gen_tag= line;
     std::cout << "Gen tag " << gen_tag << std::endl ;  
 
     fgets(line,255,ifile);
     std::string comment = line ;
     std::cout << "Gen comment " << comment << std::endl ; 
 
     int iovRunStart(0);
     fgets(line,255,ifile);
     sscanf (line, "%d", &iovRunStart);
     std::cout << "IOV START " << iovRunStart << std::endl;
     //if -1 start of time
 
     int iovRunEnd(0);
     fgets(line,255,ifile);
     sscanf (line, "%d", &iovRunEnd);
     std::cout << "IOV END " << iovRunEnd << std::endl;
     //if -1 end of time
 
     int ii = 0;
     while(fgets(line,255,ifile)) 
       {
         std::string EBorEE;
         int hashedIndex(0);
         int chStatus(0);
         std::stringstream aStrStream;
         aStrStream << line;
         aStrStream >> EBorEE >> hashedIndex >> chStatus;
         //      if(ii==0) 
         std::cout << EBorEE << " hashedIndex " << hashedIndex << " status " <<  chStatus << std::endl;
         
         if (EBorEE == "EB") 
           {
         EBDetId aEBDetId=EBDetId::unhashIndex(hashedIndex);
         status->setValue( aEBDetId, chStatus );
           }
         else if (EBorEE == "EE")
           {
         //      chStatus=1;
         EEDetId aEEDetId=EEDetId::unhashIndex(hashedIndex);
         status->setValue( aEEDetId, chStatus );
           }
         else if (EBorEE == "EBTT")
           {
         int ism=hashedIndex;
         int itt=chStatus;
 
         int ixtt=(itt-1)%4;
         int iytt=(itt-1)/4;
         int ixmin=ixtt*5;
         int iymin=iytt*5;
         int ixmax=(ixtt+1)*5-1;
         int iymax=(iytt+1)*5-1;
         for(int ieta=iymin; ieta<=iymax; ieta++) {
           for(int iphi=ixmin; iphi<=ixmax; iphi++) {
             int ixt=ieta*20+iphi+1;
             std::cout<< "killing crystal "<< ism << "/" << ixt << endl;
             EBDetId ebid(ism,ixt,EBDetId::SMCRYSTALMODE);
             status->setValue( ebid, 1 );
           }
         }
           }
         
         ii++ ;
         
       }
     
         fclose(ifile);
 
     
     /*
     std::cout <<"KILLING CHANNELS FOR CRAFT EB+16 AND EB+7"<<endl; 
 
     int ism=7;
     for(int ixt=1; ixt<=500; ixt++) {
       EBDetId ebid(ism,ixt,EBDetId::SMCRYSTALMODE);
       status->setValue( ebid, 1 );
     }
     for(int ixt=501; ixt<=900; ixt++) {
       EBDetId ebid(ism,ixt,EBDetId::SMCRYSTALMODE);
       if( ((ixt)%20==0) || ((ixt)%20>10) ){  
         status->setValue( ebid, 1 );
       }
     }
     ism=16;
     for(int ixt=501; ixt<=900; ixt++) {
       EBDetId ebid(ism,ixt,EBDetId::SMCRYSTALMODE);
       if( ((ixt)%20==0) || ((ixt)%20>10) ){  
         status->setValue( ebid, 1 );
       }
     }
 
     */
 
 
         return status;
 
 
 }

EcalGainRatios * StoreEcalCondition::readEcalGainRatiosFromFile ( const char * inputFile )

Definition at line 950 of file StoreEcalCondition.cc.

References gather_cfg::cout, PVValHelper::eta, EBDetId::ETAPHIMODE, mps_fire::i, cuy::ii, analyzePatCleaning_cfg::inputFile, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, geometryCSVtoXML::line, phi, DetId::rawId(), EcalMGPAGainRatio::setGain12Over6(), EcalMGPAGainRatio::setGain6Over1(), EcalCondObjectContainer< T >::setValue(), sm_constr_, sm_slot_, EBDetId::SMCRYSTALMODE, harvestTrackValidationPlots::str, to_string(), and EEDetId::validDetId().

Referenced by endJob().

                                                                     {
 //-------------------------------------------------------------
 
   // create gain ratios
   EcalGainRatios* gratio = new EcalGainRatios;
     
 
     FILE *inpFile; // input file
     inpFile = fopen(inputFile,"r");
     if(!inpFile) {
       edm::LogError("StoreEcalCondition")<<"*** Can not open file: "<<inputFile;
       return nullptr;
     }
 
     char line[256];
     std::ostringstream str;
               
     fgets(line,255,inpFile);
     string sm_or_all=to_string(line);
     int sm_number=0;
     sm_number=atoi(line);
     str << "sm= " << sm_number << endl ;  
 
     fgets(line,255,inpFile);
     //int nevents=atoi(line);
     
     fgets(line,255,inpFile);
     string gen_tag=to_string(line);
     str << "gen tag " << gen_tag << endl ;    
 
     fgets(line,255,inpFile);
     string cali_method=to_string(line);
     str << "cali method " << cali_method << endl ;
 
     fgets(line,255,inpFile);
     string cali_version=to_string(line);
     str << "cali version " << cali_version << endl ;
 
 
     fgets(line,255,inpFile);
     string cali_type=to_string(line);
 
     str << "cali type " << cali_type << endl ;
 
     edm::LogInfo("StoreEcalCondition") << "GainRatio file " << str.str() ;
 
 
 
     int cry_num[61200]={0};
     float g1_g12[61200]={0};
     float g6_g12[61200]={0};
     int calib_status[61200]={0};
     int dummy1=0;
     int dummy2=0;
     int hash1=0;
 
     int ii = 0;
 
     if(sm_number!=-1){
       while(fgets(line,255,inpFile)) {
     sscanf(line, "%d %d %d %f %f %d", &dummy1, &dummy2, &cry_num[ii], &g1_g12[ii], &g6_g12[ii], &calib_status[ii] );
     ii++ ;
       } 
 
     
       fclose(inpFile);           // close inp. file
 
 
 
       edm::LogInfo("StoreEcalCondition") << "Read gainRatios for " << ii << " xtals " ;
       if(ii!=1700) edm::LogWarning("StoreEcalCondition") << " Missing crystals:: missing channels will be set to 0" << endl;
       
       // Get channel ID 
       sm_constr_ = sm_number;
       
       for(int i=0; i<1700; i++){
     // EBDetId(int index1, int index2, int mode = ETAPHIMODE)
     // sm and crys index SMCRYSTALMODE index1 is SM index2 is crystal number a la H4
     EBDetId ebid(sm_slot_,cry_num[i],EBDetId::SMCRYSTALMODE);
     // cout << "ebid.rawId()"<< ebid.rawId()<< endl;
     EcalMGPAGainRatio gr;
     gr.setGain12Over6( g6_g12[i] );
     gr.setGain6Over1(g1_g12[i]/g6_g12[i]);
     gratio->setValue( ebid.rawId(), gr );
       } // loop over channels 
       
       
 
       
 
 
     } else {
       // this is for the whole Barrel 
       cout<<"mode ALL BARREL" <<endl; 
       while(fgets(line,255,inpFile)) {
     int eta=0; 
     int phi=0;
     sscanf(line, "%d %d %d %f %f",&hash1, &eta, &phi, &g1_g12[ii], &g6_g12[ii]);
     if(ii<20) cout<<"crystal eta/phi="<<eta<<"/"<<phi<<" g1_12/g6_12= "<< g1_g12[ii]<<"/"<<g6_g12[ii]<<endl;
 
     if(g1_g12[ii]<9 || g1_g12[ii]>15 ) g1_g12[ii]=12.0;
     if(g6_g12[ii]<1 || g6_g12[ii]>3 ) g6_g12[ii]=2.0;
 
     if(eta<-85|| eta>85 || eta==0) std::cout<<"error!!!"<<endl;
     if(phi<1 || phi>360) std::cout<<"error!!!"<<endl;
 
     EBDetId ebid(eta, phi,EBDetId::ETAPHIMODE);
     EcalMGPAGainRatio gr;
     gr.setGain12Over6( g6_g12[ii] );
     gr.setGain6Over1(g1_g12[ii]/g6_g12[ii]);
     gratio->setValue( ebid.rawId(), gr );
     
     ii++ ;
       }
     
       fclose(inpFile);           // close inp. file
       if(ii!=61200) edm::LogWarning("StoreEcalCondition") << " Missing crystals !!!!!!!" << endl;
       
       std::cout<< "number of crystals read:"<<ii<<endl;
 
     }
 
 
 
     for(int iX=EEDetId::IX_MIN; iX<=EEDetId::IX_MAX ;++iX) {
       for(int iY=EEDetId::IY_MIN; iY<=EEDetId::IY_MAX; ++iY) {
     // make an EEDetId since we need EEDetId::rawId() to be used as the key for the pedestals
     EcalMGPAGainRatio gr;
     gr.setGain12Over6( 2. );
     gr.setGain6Over1( 6. );
 
 
     if (EEDetId::validDetId(iX,iY,1))
       {
         EEDetId eedetidpos(iX,iY,1);
         gratio->setValue( eedetidpos.rawId(), gr );
       }
     if (EEDetId::validDetId(iX,iY,-1))
       {
         EEDetId eedetidneg(iX,iY,-1);
         gratio->setValue( eedetidneg.rawId(), gr );
       }
       }
     }
   
 
     std::cout <<" gratio pointer="<<gratio<<endl;
 
 
 
     std::cout<< "now leaving"<<endl;
 
 
 
     return gratio;
 
 }

EcalIntercalibConstants * StoreEcalCondition::readEcalIntercalibConstantsFromFile	(	const char *	inputFile,
		const char *	inputFileEE
	)

Definition at line 598 of file StoreEcalCondition.cc.

References convertFromConstructionSMToSlot(), gather_cfg::cout, mps_fire::i, cuy::ii, analyzePatCleaning_cfg::inputFile, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, geometryCSVtoXML::line, nchan, DetId::rawId(), EcalCondObjectContainer< T >::setValue(), sm_constr_, EBDetId::SMCRYSTALMODE, harvestTrackValidationPlots::str, to_string(), and EEDetId::validDetId().

Referenced by endJob().

                                                                                                      {
 //-------------------------------------------------------------
 
   EcalIntercalibConstants* ical = new EcalIntercalibConstants();
 
     
     FILE *inpFile; // input file
     inpFile = fopen(inputFile,"r");
     if(!inpFile) {
       edm::LogError("StoreEcalCondition")<<"*** Can not open file: "<<inputFile;
       return nullptr;
     }
 
     char line[256];
  
     std::ostringstream str;
 
     fgets(line,255,inpFile);
     string sm_or_all=to_string(line);
     int sm_number=0;
     int nchan=1700;
     sm_number=atoi(line);
     str << "sm= " << sm_number << endl ;  
     if(sm_number!=-1){
       nchan=1700;
     } else {
       nchan=61200;
     }
     
 
     fgets(line,255,inpFile);
     //int nevents=atoi(line); // not necessary here just for online conddb
     
     fgets(line,255,inpFile);
     string gen_tag=to_string(line);
     str << "gen tag " << gen_tag << endl ;  // should I use this? 
 
     fgets(line,255,inpFile);
     string cali_method=to_string(line);
     str << "cali method " << cali_method << endl ; // not important 
 
     fgets(line,255,inpFile);
     string cali_version=to_string(line);
     str << "cali version " << cali_version << endl ; // not important 
 
     fgets(line,255,inpFile);
     string cali_type=to_string(line);
     str << "cali type " << cali_type << endl ; // not important
 
     edm::LogInfo("StoreEcalCondition") << "Intercalibration file " << str.str() ;
 
     int sm_num[61200]={0};
     int cry_num[61200]={0};
     float calib[61200]={0};
     float calib_rms[61200]={0};
     int calib_nevents[61200]={0};
     int calib_status[61200]={0};
 
     int ii = 0;
     if(sm_number!=-1){
       while(fgets(line,255,inpFile)) {
     sscanf(line, "%d %f %f %d %d", &cry_num[ii], &calib[ii], &calib_rms[ii], &calib_nevents[ii], &calib_status[ii] );
 //       if(ii<10) { // print out only the first ten channels 
 //  cout << "cry="<<cry_num[ii]<<" calib="<<calib[ii]<<endl;
 //       }
     sm_num[ii]=sm_number;
     ii++ ;
       }
     } else {
       // this is for the whole Barrel 
       cout<<"mode ALL BARREL" <<endl; 
       while(fgets(line,255,inpFile)) {
     sscanf(line, "%d %d %f %f %d", &sm_num[ii], &cry_num[ii], &calib[ii], &calib_rms[ii], &calib_status[ii] );
     if(ii==0) cout<<"crystal "<<cry_num[ii]<<" of sm "<<sm_num[ii]<<" cali= "<< calib[ii]<<endl;
     ii++ ;
       }
     }
 
     //    inf.close();           // close inp. file
     fclose(inpFile);           // close inp. file
 
     edm::LogInfo("StoreEcalCondition") << "Read intercalibrations for " << ii << " xtals " ; 
 
     cout << " I read the calibrations for "<< ii<< " crystals " << endl;
     if(ii!=nchan) edm::LogWarning("StoreEcalCondition") << "Some crystals missing. Missing channels will be set to 0" << endl;
 
     // Get channel ID 
     
     sm_constr_ = sm_number;
 
 
     // Set the data
     for(int i=0; i<nchan; i++){
     
     // EBDetId(int index1, int index2, int mode = ETAPHIMODE)
     // sm and crys index SMCRYSTALMODE index1 is SM index2 is crystal number a la H4
     
       int slot_num=convertFromConstructionSMToSlot(sm_num[i],-1);
       EBDetId ebid(slot_num,cry_num[i],EBDetId::SMCRYSTALMODE);
 
       ical->setValue( ebid.rawId(), calib[i]  );
 
       if(i==0) cout<<"crystal "<<cry_num[i]<<" of sm "<<sm_num[i]<< " in slot " <<slot_num<<" calib= "<< calib[i]<<endl;
 
     } // loop over channels 
 
     cout<<"loop on channels done" << endl;
 
     FILE *inpFileEE; // input file
     inpFileEE = fopen(inputFileEE,"r");
     if(!inpFileEE) {
       edm::LogError("StoreEcalCondition")<<"*** Can not open file: "<<inputFile;
 
     // dummy endcap data
 
     for(int iX=EEDetId::IX_MIN; iX<=EEDetId::IX_MAX ;++iX) 
       {
     for(int iY=EEDetId::IY_MIN; iY<=EEDetId::IY_MAX; ++iY) 
       {
         // make an EEDetId since we need EEDetId::rawId() to be used as the key for the pedestals
         if (EEDetId::validDetId(iX,iY,1)) 
           {
         EEDetId eedetidpos(iX,iY,1);
         ical->setValue( eedetidpos, 1.0 );
           }
         if (EEDetId::validDetId(iX,iY,-1)) 
           {
         EEDetId eedetidneg(iX,iY,-1);
         ical->setValue( eedetidneg, 1.0 );
           }
       }
       }
 
     } else {
       cout<<"... now reading EE file ..." <<endl; 
 
       int ii=0;
       while(fgets(line,255,inpFileEE)) {
     int iz,ix,iy;
     float calibee;
     sscanf(line, "%d %d %d %f", &iz, &ix, &iy, &calibee );
     if(ii<=0) cout<<"crystal "<<iz<<"/"<<ix<<"/"<<iy<<" cali="<< calibee<<endl;
 
     if (EEDetId::validDetId(ix,iy,iz)) 
       {
         EEDetId eedetid(ix,iy,iz);
         ical->setValue( eedetid, calibee );
       }
     
     ii++ ;
       }
 
 
       fclose(inpFileEE);           // close inp. file
 
     }
 
     cout<<"loop on EE channels done" << endl;
 
 
 
   return ical;
 }

EcalIntercalibConstantsMC * StoreEcalCondition::readEcalIntercalibConstantsMCFromFile	(	const char *	inputFile,
		const char *	inputFileEE
	)

Definition at line 764 of file StoreEcalCondition.cc.

References convertFromConstructionSMToSlot(), gather_cfg::cout, mps_fire::i, cuy::ii, analyzePatCleaning_cfg::inputFile, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, geometryCSVtoXML::line, nchan, DetId::rawId(), EcalCondObjectContainer< T >::setValue(), sm_constr_, EBDetId::SMCRYSTALMODE, harvestTrackValidationPlots::str, to_string(), and EEDetId::validDetId().

Referenced by endJob().

                                                                                                        {
 //-------------------------------------------------------------
 
   EcalIntercalibConstantsMC* ical = new EcalIntercalibConstantsMC();
 
     
     FILE *inpFile; // input file
     inpFile = fopen(inputFile,"r");
     if(!inpFile) {
       edm::LogError("StoreEcalCondition")<<"*** Can not open file: "<<inputFile;
       return nullptr;
     }
 
     char line[256];
  
     std::ostringstream str;
 
     fgets(line,255,inpFile);
     string sm_or_all=to_string(line);
     int sm_number=0;
     int nchan=1700;
     sm_number=atoi(line);
     str << "sm= " << sm_number << endl ;  
     if(sm_number!=-1){
       nchan=1700;
     } else {
       nchan=61200;
     }
     
 
     fgets(line,255,inpFile);
     //int nevents=atoi(line); // not necessary here just for online conddb
     
     fgets(line,255,inpFile);
     string gen_tag=to_string(line);
     str << "gen tag " << gen_tag << endl ;  // should I use this? 
 
     fgets(line,255,inpFile);
     string cali_method=to_string(line);
     str << "cali method " << cali_method << endl ; // not important 
 
     fgets(line,255,inpFile);
     string cali_version=to_string(line);
     str << "cali version " << cali_version << endl ; // not important 
 
     fgets(line,255,inpFile);
     string cali_type=to_string(line);
     str << "cali type " << cali_type << endl ; // not important
 
     edm::LogInfo("StoreEcalCondition") << "Intercalibration file " << str.str() ;
 
     int sm_num[61200]={0};
     int cry_num[61200]={0};
     float calib[61200]={0};
     float calib_rms[61200]={0};
     int calib_nevents[61200]={0};
     int calib_status[61200]={0};
 
     int ii = 0;
     if(sm_number!=-1){
       while(fgets(line,255,inpFile)) {
     sscanf(line, "%d %f %f %d %d", &cry_num[ii], &calib[ii], &calib_rms[ii], &calib_nevents[ii], &calib_status[ii] );
 //       if(ii<10) { // print out only the first ten channels 
 //  cout << "cry="<<cry_num[ii]<<" calib="<<calib[ii]<<endl;
 //       }
     sm_num[ii]=sm_number;
     ii++ ;
       }
     } else {
       // this is for the whole Barrel 
       cout<<"mode ALL BARREL" <<endl; 
       while(fgets(line,255,inpFile)) {
     sscanf(line, "%d %d %f %f %d", &sm_num[ii], &cry_num[ii], &calib[ii], &calib_rms[ii], &calib_status[ii] );
     if(ii==0) cout<<"crystal "<<cry_num[ii]<<" of sm "<<sm_num[ii]<<" cali= "<< calib[ii]<<endl;
     ii++ ;
       }
     }
 
     //    inf.close();           // close inp. file
     fclose(inpFile);           // close inp. file
 
     edm::LogInfo("StoreEcalCondition") << "Read intercalibrations for " << ii << " xtals " ; 
 
     cout << " I read the calibrations for "<< ii<< " crystals " << endl;
     if(ii!=nchan) edm::LogWarning("StoreEcalCondition") << "Some crystals missing. Missing channels will be set to 0" << endl;
 
     // Get channel ID 
     
     sm_constr_ = sm_number;
 
 
     // Set the data
     for(int i=0; i<nchan; i++){
     
     // EBDetId(int index1, int index2, int mode = ETAPHIMODE)
     // sm and crys index SMCRYSTALMODE index1 is SM index2 is crystal number a la H4
     
       int slot_num=convertFromConstructionSMToSlot(sm_num[i],-1);
       EBDetId ebid(slot_num,cry_num[i],EBDetId::SMCRYSTALMODE);
 
       ical->setValue( ebid.rawId(), calib[i]  );
 
       if(i==0) cout<<"crystal "<<cry_num[i]<<" of sm "<<sm_num[i]<< " in slot " <<slot_num<<" calib= "<< calib[i]<<endl;
 
     } // loop over channels 
 
     cout<<"loop on channels done" << endl;
 
     FILE *inpFileEE; // input file
     inpFileEE = fopen(inputFileEE,"r");
     if(!inpFileEE) {
       edm::LogError("StoreEcalCondition")<<"*** Can not open file: "<<inputFile;
 
     // dummy endcap data
 
     for(int iX=EEDetId::IX_MIN; iX<=EEDetId::IX_MAX ;++iX) 
       {
     for(int iY=EEDetId::IY_MIN; iY<=EEDetId::IY_MAX; ++iY) 
       {
         // make an EEDetId since we need EEDetId::rawId() to be used as the key for the pedestals
         if (EEDetId::validDetId(iX,iY,1)) 
           {
         EEDetId eedetidpos(iX,iY,1);
         ical->setValue( eedetidpos, 1.0 );
           }
         if (EEDetId::validDetId(iX,iY,-1)) 
           {
         EEDetId eedetidneg(iX,iY,-1);
         ical->setValue( eedetidneg, 1.0 );
           }
       }
       }
 
     } else {
       cout<<"... now reading EE file ..." <<endl; 
 
       int ii=0;
       while(fgets(line,255,inpFileEE)) {
     int iz,ix,iy;
     float calibee;
     sscanf(line, "%d %d %d %f", &iz, &ix, &iy, &calibee );
     if(ii<=0) cout<<"crystal "<<iz<<"/"<<ix<<"/"<<iy<<" cali="<< calibee<<endl;
 
     if (EEDetId::validDetId(ix,iy,iz)) 
       {
         EEDetId eedetid(ix,iy,iz);
         ical->setValue( eedetid, calibee );
       }
     
     ii++ ;
       }
 
 
       fclose(inpFileEE);           // close inp. file
 
     }
 
     cout<<"loop on EE channels done" << endl;
 
 
 
   return ical;
 }

EcalPFRecHitThresholds * StoreEcalCondition::readEcalPFRecHitThresholdsFromFile	(	const char *	inputFile,
		const char *	inputFileEE
	)

Definition at line 520 of file StoreEcalCondition.cc.

References gather_cfg::cout, cuy::ii, analyzePatCleaning_cfg::inputFile, geometryCSVtoXML::line, DetId::rawId(), EcalCondObjectContainer< T >::setValue(), scrapingFilter_cfi::thresh, EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by endJob().

                                                                                                     {
 //-------------------------------------------------------------
 
   EcalPFRecHitThresholds* ical = new EcalPFRecHitThresholds();
 
     
     FILE *inpFile; // input file
     inpFile = fopen(inputFile,"r");
     if(!inpFile) {
       edm::LogError("StoreEcalCondition")<<"*** Can not open file: "<<inputFile;
       return nullptr;
     }
 
     char line[256];
  
 
     int ieta=0;
     int iphi=0;
     int ix=0;
     int iy=0;
     int iz=0;
     
     float thresh=0;
 
 
     int ii = 0;
     while(fgets(line,255,inpFile)) {
       sscanf(line, "%d %d %f ", &ieta, &iphi, &thresh);
       if(ii==0) cout<<"crystal "<<ieta<<"/"<<iphi<<" Thresh= "<< thresh<<endl;
       
       if (EBDetId::validDetId(ieta,iphi)) {
     EBDetId ebid(ieta,iphi);
     ical->setValue( ebid.rawId(), thresh  );
     ii++ ;
       }
     }
     
 
     //    inf.close();           // close inp. file
     fclose(inpFile);           // close inp. file
     
     edm::LogInfo("StoreEcalCondition") << "Read PF RecHits for " << ii << " xtals " ; 
     
     cout << " I read the thresholds for "<< ii<< " crystals " << endl;
     
     
     
     FILE *inpFileEE; // input file
     inpFileEE = fopen(inputFileEE,"r");
     if(!inpFileEE) {
       edm::LogError("StoreEcalCondition")<<"*** Can not open file: "<<inputFileEE;
       return nullptr;
     } 
     ii=0;
     while(fgets(line,255,inpFileEE)) {
       sscanf(line, "%d %d %d %f ", &ix,&iy,&iz,  &thresh);
       if(ii==0) cout<<"crystal "<<ix<<"/"<<iy<<"/"<<iz<<" Thresh= "<< thresh<<endl;
       if (EEDetId::validDetId(ix,iy,iz)) {
     EEDetId eeid(ix,iy,iz);
     ical->setValue( eeid.rawId(), thresh  );
     ii++ ;
       }
     }
     
     
     //    inf.close();           // close inp. file
     fclose(inpFileEE);           // close inp. file
     
     
     cout<<"loop on EE channels done - number of crystals =" <<ii<< std::endl;
     
     
     
     return ical;
     
 }

EcalTBWeights * StoreEcalCondition::readEcalTBWeightsFromFile ( const char * inputFile )

Definition at line 312 of file StoreEcalCondition.cc.

References EcalWeightSet::getChi2WeightsAfterGainSwitch(), EcalWeightSet::getChi2WeightsBeforeGainSwitch(), EcalWeightSet::getWeightsAfterGainSwitch(), EcalWeightSet::getWeightsBeforeGainSwitch(), mps_fire::i, analyzePatCleaning_cfg::inputFile, gen::k, LogDebug, hgc_digi::nSamples, EcalTBWeights::setValue(), sm_constr_, harvestTrackValidationPlots::str, and groupFilesInBlocks::temp.

Referenced by endJob().

                                                                    {
   //-------------------------------------------------------------
   
   // Zabi code to be written here 
   
   EcalTBWeights* tbwgt = new EcalTBWeights();
   
   std::ifstream WeightsFileTB(inputFile);
   if(!WeightsFileTB.is_open()) {
     edm::LogError("StoreEcalCondition")<< "*** Can not open file: "<< inputFile ;
     return nullptr;
   }
 
   int smnumber=-99999;
 
   std::ostringstream str;
   WeightsFileTB >> smnumber;
   if (smnumber == -99999)
     return nullptr;
 
   str << "sm= " << smnumber << endl;
 
   char temp[256];
   //Reading the other 5 header lines containing various informations
   for (int i=0;i<=5;i++)
     {
       WeightsFileTB.getline(temp,255);
       str << temp << endl;
     }
 
   edm::LogInfo("StoreEcalCondition") << "Weights file " << str.str() ;
 
   int ngroups=0;
   while (WeightsFileTB.good())
     {
       int igroup_ID = -99999;
       int nSamples  = -99999;
       int nTdcBins  = -99999;
       
       WeightsFileTB >> igroup_ID >> nSamples >> nTdcBins;
       if (igroup_ID == -99999 || nSamples == -99999 || nTdcBins == -99999)
     break;
       
       std::ostringstream str;
       str << "Igroup=" << igroup_ID << " Nsamples=" << nSamples << " NTdcBins=" << nTdcBins <<"\n" ;
       
       for (int iTdcBin = 0; iTdcBin < nTdcBins; iTdcBin++) {
     
     EcalWeightSet wgt; // one set of weights
     EcalWeightSet::EcalWeightMatrix& wgt1   = wgt.getWeightsBeforeGainSwitch();
     EcalWeightSet::EcalWeightMatrix& wgt2   = wgt.getWeightsAfterGainSwitch();  
     EcalWeightSet::EcalChi2WeightMatrix& chisq1 = wgt.getChi2WeightsBeforeGainSwitch();
     EcalWeightSet::EcalChi2WeightMatrix& chisq2 = wgt.getChi2WeightsAfterGainSwitch();
     
 //  std::vector<EcalWeight> wamp, wped, wtime; //weights before gain switch
 //  std::vector<EcalWeight> wamp2, wped2, wtime2; //weights after gain switch
     
     //WEIGHTS BEFORE GAIN SWITCH
     //Amplitude weights 
     for(int j = 0; j < nSamples; ++j) {
       double ww = 0.0; WeightsFileTB >> ww;
       wgt1(0,j)=ww;
       str << ww << " ";
     }// loop Samples
     str << std::endl;
     
     //Pedestal weights
     for(int j = 0; j < nSamples; ++j) {
       double ww = 0.0; WeightsFileTB >> ww;
       wgt1(1,j)=ww;
       str << ww << " ";
     }//loop Samples
     str << std::endl;
     
     //Timing weights
     for(int j = 0; j < nSamples; ++j) {
       double ww = 0.0; WeightsFileTB >> ww;
       wgt1(2,j)=ww;
       str << ww << " ";
     }//loop Samples
     str << std::endl;
     
     for(int j = 0; j < nSamples; ++j) {
       // fill chi2 matrix
       //std::vector<EcalWeight> vChi2; // row of chi2 matrix
       for(int k = 0; k < nSamples; ++k) {
         double ww = 0.0; WeightsFileTB >> ww;
         chisq1(j,k)=ww;
         str << ww << " ";
       } //loop samples
       str << std::endl;
     }//loop lines
     
     //WEIGHTS AFTER GAIN SWITCH
     for(int j = 0; j < nSamples; ++j) {
       double ww = 0.0; WeightsFileTB >> ww;
       wgt2(0,j)=ww;
       str << ww << " ";
     }// loop Samples
     str << std::endl;
     
     //Pedestal weights
     for(int j = 0; j < nSamples; ++j) {
       double ww = 0.0; WeightsFileTB >> ww;
       wgt2(1,j)=ww;
       str << ww << " ";
     }//loop Samples
     str << std::endl;
     
     //Timing weights
     for(int j = 0; j < nSamples; ++j) {
       double ww = 0.0; WeightsFileTB >> ww;
       wgt2(2,j)=ww;
       str << ww << " ";
     }//loop Samples
     str << std::endl;
     
     for(int j = 0; j < nSamples; ++j) {
       // fill chi2 matrix
       //std::vector<EcalWeight> vChi2; // row of chi2 matrix
       for(int k = 0; k < nSamples; ++k) {
         double ww = 0.0; WeightsFileTB >> ww;
         chisq2(j,k)=ww;
         str << ww << " ";
       } //loop samples
       str << std::endl;
     }//loop lines
 
     LogDebug("StoreEcalCondition") << str.str();
       
     //modif-27-07-06 tdcid should start from 1
     tbwgt->setValue(std::make_pair( igroup_ID , iTdcBin+1 ), wgt);
       }//loop Tdc bins
       ngroups++;
     }//loop groupID
 
   sm_constr_ = smnumber;
 
 
   edm::LogInfo("StoreEcalCondition") << "Weights for " << ngroups << " groups written into DB" ;
   return tbwgt;
 }

EcalWeightXtalGroups * StoreEcalCondition::readEcalWeightXtalGroupsFromFile ( const char * inputFile )

Definition at line 246 of file StoreEcalCondition.cc.

References mps_fire::i, analyzePatCleaning_cfg::inputFile, LogDebug, DetId::rawId(), EcalCondObjectContainer< T >::setValue(), sm_constr_, sm_slot_, EBDetId::SMCRYSTALMODE, harvestTrackValidationPlots::str, and groupFilesInBlocks::temp.

Referenced by endJob().

                                                                           {
 //-------------------------------------------------------------
 
 // Code taken from EcalWeightTools/test/MakeOfflineDbFromAscii.cpp
 
   EcalWeightXtalGroups* xtalGroups = new EcalWeightXtalGroups();
   std::ifstream groupid_in(inputFile);
 
   if(!groupid_in.is_open()) {
     edm::LogError("StoreEcalCondition")<< "*** Can not open file: "<< inputFile ;
     return nullptr;
   }  
 
   int smnumber=-99999;
 
   std::ostringstream str;
   groupid_in >> smnumber;
   if (smnumber == -99999) {
     edm::LogError("StoreEcalCondition") << "ERROR: SM number not found in file" << endl;
     return nullptr;
   }
   str << "sm= " << smnumber << endl;
   sm_constr_ = smnumber;
 
   char temp[256];
   //Reading the other 5 header lines containing various informations
   for (int i=0;i<=5;i++) {
     groupid_in.getline(temp,255);
     str << temp << endl;
   }
 
   // Skip the nGroup/Mean line
   groupid_in.getline(temp, 255);
   str << temp << endl;
 
   edm::LogInfo("StoreEcalCondition") << "GROUPID file " << str.str() ;
 
   int xtals = 0;
   int xtal, ietaf, iphif, groupID;
   while (groupid_in.good()) {
     groupid_in >> xtal >> ietaf >> iphif >> groupID;
     if (groupid_in.eof()) { break; }
     
     LogDebug("StoreEcalCondition") << "XTAL=" << xtal << " ETA=" << ietaf << " PHI=" << iphif 
                        << " GROUP=" << groupID ;
 
     //EBDetId ebid(ieta,iphi);  
     EBDetId ebid(sm_slot_,xtal,EBDetId::SMCRYSTALMODE); 
     // xtalGroups->setValue(ebid.rawId(), EcalXtalGroupId( ebid.hashedIndex()) );
     xtalGroups->setValue(ebid.rawId(), EcalXtalGroupId( groupID ) );
     xtals++;
   }//loop iphi
 
   if (xtals != 1700) {
     edm::LogError("StoreEcalCondition") << "ERROR:  GROUPID file did not contain data for 1700 crystals" << endl;
     return nullptr;
   }
   
   edm::LogInfo("StoreEcalCondition") << "Groups for " << xtals << " xtals written into DB" ;
   sm_constr_ = smnumber;
   
   return xtalGroups;
 }

std::string StoreEcalCondition::to_string ( char value[] )

inlineprivate

Definition at line 74 of file StoreEcalCondition.h.

References relativeConstraints::value.

Referenced by readEcalADCToGeVConstantFromFile(), readEcalGainRatiosFromFile(), readEcalIntercalibConstantsFromFile(), and readEcalIntercalibConstantsMCFromFile().

                                      {
     std::ostringstream streamOut;
     streamOut << value;
     return streamOut.str();
   }

void StoreEcalCondition::writeToLogFile	(	std::string	,
		std::string	,
		unsigned long	long
	)

Definition at line 161 of file StoreEcalCondition.cc.

References fillHeader(), RecoTauValidation_cfi::header, createfilelist::int, logfile_, GetRecoTauVFromDQM_MC_cff::outFile, sm_constr_, and sm_slot_.

                                                                                     {
 //-------------------------------------------------------------
   
     FILE *outFile; // output log file for appending 
     outFile = fopen(logfile_.c_str(),"a");  
     if(!outFile) {
       edm::LogError("StoreEcalCondition") <<"*** Can not open file: " << logfile_;
       return;
     }
     char header[256];
     fillHeader(header);
     char appendMode[10];
     if (since != 0)
       sprintf(appendMode,"append");
     else
       sprintf(appendMode,"create");
 
     fprintf(outFile, "%s %s condition from file %s written into DB for SM %d (mapped to SM %d) in %s mode (since run %u)\n", 
         header, a.c_str(),b .c_str(), sm_constr_, sm_slot_, appendMode, (unsigned int)since);
 
     fclose(outFile);           // close out file
 
 }

void StoreEcalCondition::writeToLogFileResults ( char * arg )

Definition at line 185 of file StoreEcalCondition.cc.

References fillHeader(), RecoTauValidation_cfi::header, logfile_, and GetRecoTauVFromDQM_MC_cff::outFile.

Referenced by endJob().

                                                         {
 //-------------------------------------------------------------
   
     FILE *outFile; // output log file for appending 
     outFile = fopen(logfile_.c_str(),"a");  
     if(!outFile) {
       edm::LogError("StoreEcalCondition")<<"*** Can not open file: " << logfile_;
       return;
     }
     char header[256];
     fillHeader(header);
     fprintf(outFile, "%s %s\n", header,arg);
     fclose(outFile);           // close out file
 }