|
|
|
#include <StoreEcalCondition.h>
Definition at line 38 of file StoreEcalCondition.h.
| 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_, and convertXMLtoSQLite_cfg::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 | ( | ) |
Definition at line 142 of file StoreEcalCondition.cc.
{
}
| void StoreEcalCondition::analyze | ( | const edm::Event & | evt, |
| const edm::EventSetup & | evtSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 147 of file StoreEcalCondition.cc.
{
//-------------------------------------------------------------
}
| int StoreEcalCondition::convertFromConstructionSMToSlot | ( | int | sm_constr, |
| int | sm_slot | ||
| ) |
Definition at line 844 of file StoreEcalCondition.cc.
References query::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 | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 41 of file StoreEcalCondition.cc.
References gather_cfg::cout, i, inpFileName_, inpFileNameEE_, edm::Service< T >::isAvailable(), python::rootplot::argparse::message, objectName_, readEcalADCToGeVConstantFromFile(), readEcalChannelStatusFromFile(), readEcalGainRatiosFromFile(), readEcalIntercalibConstantsFromFile(), readEcalIntercalibConstantsMCFromFile(), readEcalTBWeightsFromFile(), readEcalWeightXtalGroupsFromFile(), since_, 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] == "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 194 of file StoreEcalCondition.cc.
References cond::rpcobgas::time, and hcal_dqm_sourceclient-file_cfg::user.
Referenced by writeToLogFile(), and writeToLogFileResults().
| EcalADCToGeVConstant * StoreEcalCondition::readEcalADCToGeVConstantFromFile | ( | const char * | inputFile | ) |
Definition at line 451 of file StoreEcalCondition.cc.
References analyzePatCleaning_cfg::inputFile, geometryCSVtoXML::line, LogDebug, sm_constr_, 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 0;
}
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 1024 of file StoreEcalCondition.cc.
References comment, gather_cfg::cout, Exception, ecalpyutils::hashedIndex(), compare_using_db::ifile, 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, ntuplemaker::status, EBDetId::unhashIndex(), EEDetId::unhashIndex(), EEDetId::validDetId(), and EBDetId::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 865 of file StoreEcalCondition.cc.
References gather_cfg::cout, eta(), EBDetId::ETAPHIMODE, i, 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, 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 0;
}
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 513 of file StoreEcalCondition.cc.
References calib, convertFromConstructionSMToSlot(), gather_cfg::cout, i, 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, 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 0;
}
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 679 of file StoreEcalCondition.cc.
References calib, convertFromConstructionSMToSlot(), gather_cfg::cout, i, 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, 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 0;
}
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;
}
| EcalTBWeights * StoreEcalCondition::readEcalTBWeightsFromFile | ( | const char * | inputFile | ) |
Definition at line 305 of file StoreEcalCondition.cc.
References EcalWeightSet::getChi2WeightsAfterGainSwitch(), EcalWeightSet::getChi2WeightsBeforeGainSwitch(), EcalWeightSet::getWeightsAfterGainSwitch(), EcalWeightSet::getWeightsBeforeGainSwitch(), i, analyzePatCleaning_cfg::inputFile, j, gen::k, LogDebug, EcalTBWeights::setValue(), sm_constr_, 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 0;
}
int smnumber=-99999;
std::ostringstream str;
WeightsFileTB >> smnumber;
if (smnumber == -99999)
return 0;
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 239 of file StoreEcalCondition.cc.
References i, analyzePatCleaning_cfg::inputFile, LogDebug, DetId::rawId(), EcalCondObjectContainer< T >::setValue(), sm_constr_, sm_slot_, EBDetId::SMCRYSTALMODE, 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 0;
}
int smnumber=-99999;
std::ostringstream str;
groupid_in >> smnumber;
if (smnumber == -99999) {
edm::LogError("StoreEcalCondition") << "ERROR: SM number not found in file" << endl;
return 0;
}
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 0;
}
edm::LogInfo("StoreEcalCondition") << "Groups for " << xtals << " xtals written into DB" ;
sm_constr_ = smnumber;
return xtalGroups;
}
| std::string StoreEcalCondition::to_string | ( | char | value[] | ) | [inline, private] |
Definition at line 72 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 154 of file StoreEcalCondition.cc.
References fillHeader(), errorMatrix2Lands::header, logfile_, 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 178 of file StoreEcalCondition.cc.
References fillHeader(), errorMatrix2Lands::header, and logfile_.
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
}
std::vector< std::string > StoreEcalCondition::inpFileName_ [private] |
Definition at line 64 of file StoreEcalCondition.h.
Referenced by endJob(), and StoreEcalCondition().
std::vector< std::string > StoreEcalCondition::inpFileNameEE_ [private] |
Definition at line 65 of file StoreEcalCondition.h.
Referenced by endJob(), and StoreEcalCondition().
std::string StoreEcalCondition::logfile_ [private] |
Definition at line 70 of file StoreEcalCondition.h.
Referenced by StoreEcalCondition(), writeToLogFile(), and writeToLogFileResults().
std::vector< std::string > StoreEcalCondition::objectName_ [private] |
Definition at line 62 of file StoreEcalCondition.h.
Referenced by endJob(), and StoreEcalCondition().
std::string StoreEcalCondition::prog_name_ [private] |
Definition at line 66 of file StoreEcalCondition.h.
Referenced by StoreEcalCondition().
std::vector< unsigned long long > StoreEcalCondition::since_ [private] |
Definition at line 69 of file StoreEcalCondition.h.
Referenced by endJob(), and StoreEcalCondition().
int StoreEcalCondition::sm_constr_ [private] |
Definition at line 67 of file StoreEcalCondition.h.
Referenced by readEcalADCToGeVConstantFromFile(), readEcalGainRatiosFromFile(), readEcalIntercalibConstantsFromFile(), readEcalIntercalibConstantsMCFromFile(), readEcalTBWeightsFromFile(), readEcalWeightXtalGroupsFromFile(), StoreEcalCondition(), and writeToLogFile().
int StoreEcalCondition::sm_slot_ [private] |
Definition at line 68 of file StoreEcalCondition.h.
Referenced by readEcalGainRatiosFromFile(), readEcalWeightXtalGroupsFromFile(), StoreEcalCondition(), and writeToLogFile().
1.7.3