|
|
|
#include <CSCNoiseMatrixConditions.h>
Public Types | |
| typedef const CSCNoiseMatrix * | ReturnType |
Public Member Functions | |
| CSCNoiseMatrixConditions (const edm::ParameterSet &) | |
| ReturnType | produceNoiseMatrix (const CSCNoiseMatrixRcd &) |
| ~CSCNoiseMatrixConditions () | |
Static Public Member Functions | |
| static CSCNoiseMatrix * | prefillNoiseMatrix () |
Private Member Functions | |
| void | setIntervalFor (const edm::eventsetup::EventSetupRecordKey &, const edm::IOVSyncValue &, edm::ValidityInterval &) |
Private Attributes | |
| CSCNoiseMatrix * | cnMatrix |
Definition at line 19 of file CSCNoiseMatrixConditions.h.
| typedef const CSCNoiseMatrix* CSCNoiseMatrixConditions::ReturnType |
Definition at line 27 of file CSCNoiseMatrixConditions.h.
| CSCNoiseMatrixConditions::CSCNoiseMatrixConditions | ( | const edm::ParameterSet & | iConfig | ) |
Definition at line 407 of file CSCNoiseMatrixConditions.cc.
References cnMatrix, prefillNoiseMatrix(), produceNoiseMatrix(), and edm::ESProducer::setWhatProduced().
{
//the following line is needed to tell the framework what
// data is being produced
cnMatrix = prefillNoiseMatrix();
// added by Zhen (changed since 1_2_0)
setWhatProduced(this,&CSCNoiseMatrixConditions::produceNoiseMatrix);
findingRecord<CSCNoiseMatrixRcd>();
//now do what ever other initialization is needed
}
| CSCNoiseMatrixConditions::~CSCNoiseMatrixConditions | ( | ) |
Definition at line 419 of file CSCNoiseMatrixConditions.cc.
References cnMatrix.
{
// do anything here that needs to be done at desctruction time
// (e.g. close files, deallocate resources etc.)
delete cnMatrix;
}
| CSCNoiseMatrix * CSCNoiseMatrixConditions::prefillNoiseMatrix | ( | ) | [static] |
Definition at line 5 of file CSCNoiseMatrixConditions.cc.
References benchmark_cfg::cerr, CSCDetId, cmsRelvalreport::exit, recoMuon::in, CSCNoiseMatrix::matrix, CSCDetId::maxChamberId(), CSCDetId::maxEndcapId(), CSCDetId::maxLayerId(), CSCDetId::maxRingId(), CSCDetId::maxStationId(), CSCDetId::minChamberId(), CSCDetId::minEndcapId(), CSCDetId::minLayerId(), CSCDetId::minRingId(), and CSCDetId::minStationId().
Referenced by CSCNoiseMatrixConditions().
{
int old_chamber_id,old_strip,new_chamber_id,new_strip;
float old_elm33,old_elm34, old_elm44, old_elm35, old_elm45, old_elm55;
float old_elm46, old_elm56, old_elm66, old_elm57, old_elm67, old_elm77;
std::vector<int> old_cham_id;
std::vector<int> old_strips;
std::vector<float> old_elem33;
std::vector<float> old_elem34;
std::vector<float> old_elem44;
std::vector<float> old_elem45;
std::vector<float> old_elem35;
std::vector<float> old_elem55;
std::vector<float> old_elem46;
std::vector<float> old_elem56;
std::vector<float> old_elem66;
std::vector<float> old_elem57;
std::vector<float> old_elem67;
std::vector<float> old_elem77;
float new_elm33,new_elm34, new_elm44, new_elm35, new_elm45, new_elm55;
float new_elm46, new_elm56, new_elm66, new_elm57, new_elm67, new_elm77;
std::vector<int> new_cham_id;
std::vector<int> new_strips;
std::vector<float> new_elem33;
std::vector<float> new_elem34;
std::vector<float> new_elem44;
std::vector<float> new_elem45;
std::vector<float> new_elem35;
std::vector<float> new_elem55;
std::vector<float> new_elem46;
std::vector<float> new_elem56;
std::vector<float> new_elem66;
std::vector<float> new_elem57;
std::vector<float> new_elem67;
std::vector<float> new_elem77;
const CSCDetId& detId = CSCDetId();
CSCNoiseMatrix * cnmatrix = new CSCNoiseMatrix();
int max_istrip,id_layer,max_ring,max_cham;
unsigned int old_nrlines=0;
unsigned int new_nrlines=0;
std::ifstream olddata;
olddata.open("old_matrix.dat",std::ios::in);
if(!olddata) {
std::cerr <<"Error: old_matrix.dat -> no such file!"<< std::endl;
exit(1);
}
while (!olddata.eof() ) {
olddata >> old_chamber_id >> old_strip >> old_elm33 >> old_elm34 >> old_elm44 >> old_elm35 >> old_elm45 >> old_elm55 >> old_elm46 >> old_elm56 >> old_elm66 >> old_elm57 >> old_elm67 >> old_elm77;
old_cham_id.push_back(old_chamber_id);
old_strips.push_back(old_strip);
old_elem33.push_back(old_elm33);
old_elem34.push_back(old_elm34);
old_elem44.push_back(old_elm44);
old_elem35.push_back(old_elm35);
old_elem45.push_back(old_elm45);
old_elem55.push_back(old_elm55);
old_elem46.push_back(old_elm46);
old_elem56.push_back(old_elm56);
old_elem66.push_back(old_elm66);
old_elem57.push_back(old_elm57);
old_elem67.push_back(old_elm67);
old_elem77.push_back(old_elm77);
old_nrlines++;
}
olddata.close();
std::ifstream newdata;
newdata.open("new_matrix.txt",std::ios::in);
if(!newdata) {
std::cerr <<"Error: new_matrix.txt -> no such file!"<< std::endl;
exit(1);
}
while (!newdata.eof() ) {
newdata >> new_chamber_id >> new_strip >> new_elm33 >> new_elm34 >> new_elm44 >> new_elm35 >> new_elm45 >> new_elm55 >> new_elm46 >> new_elm56 >> new_elm66 >> new_elm57 >> new_elm67 >> new_elm77 ;
new_cham_id.push_back(new_chamber_id);
new_strips.push_back(new_strip);
new_elem33.push_back(new_elm33);
new_elem34.push_back(new_elm34);
new_elem44.push_back(new_elm44);
new_elem35.push_back(new_elm35);
new_elem45.push_back(new_elm45);
new_elem55.push_back(new_elm55);
new_elem46.push_back(new_elm46);
new_elem56.push_back(new_elm56);
new_elem66.push_back(new_elm66);
new_elem57.push_back(new_elm57);
new_elem67.push_back(new_elm67);
new_elem77.push_back(new_elm77);
new_nrlines++;
}
newdata.close();
//endcap=1 to 2,station=1 to 4, ring=1 to 4,chamber=1 to 36,layer=1 to 6
for(int iendcap=detId.minEndcapId(); iendcap<=detId.maxEndcapId(); iendcap++){
for(int istation=detId.minStationId() ; istation<=detId.maxStationId(); istation++){
max_ring=detId.maxRingId();
//station 4 ring 4 not there(36 chambers*2 missing)
//3 rings max this way of counting (ME1a & b)
if(istation==1) max_ring=3;
if(istation==2) max_ring=2;
if(istation==3) max_ring=2;
if(istation==4) max_ring=1;
for(int iring=detId.minRingId(); iring<=max_ring; iring++){
max_istrip=80;
max_cham=detId.maxChamberId();
if(istation==1 && iring==1) max_cham=36;
if(istation==1 && iring==2) max_cham=36;
if(istation==1 && iring==3) max_cham=36;
if(istation==2 && iring==1) max_cham=18;
if(istation==2 && iring==2) max_cham=36;
if(istation==3 && iring==1) max_cham=18;
if(istation==3 && iring==2) max_cham=36;
if(istation==4 && iring==1) max_cham=18;
for(int ichamber=detId.minChamberId(); ichamber<=max_cham; ichamber++){
for(int ilayer=detId.minLayerId(); ilayer<=detId.maxLayerId(); ilayer++){
//station 1 ring 3 has 64 strips per layer instead of 80
if(istation==1 && iring==3) max_istrip=64;
std::vector<CSCNoiseMatrix::Item> itemvector;
itemvector.resize(max_istrip);
id_layer = 100000*iendcap + 10000*istation + 1000*iring + 10*ichamber + ilayer;
for(int istrip=0;istrip<max_istrip;istrip++){
if(istation==1 && iring==1){
itemvector[istrip].elem33 = 7.86675;
itemvector[istrip].elem34 = 2.07075;
itemvector[istrip].elem44 = 6.93875;
itemvector[istrip].elem35 = 1.42525;
itemvector[istrip].elem45 = 2.51025;
itemvector[istrip].elem55 = 7.93975;
itemvector[istrip].elem46 = 0.94725;
itemvector[istrip].elem56 = 2.39275;
itemvector[istrip].elem66 = 6.46475;
itemvector[istrip].elem57 = 1.86325;
itemvector[istrip].elem67 = 2.08025;
itemvector[istrip].elem77 = 6.67975;
cnmatrix->matrix[id_layer]=itemvector;
}
if(istation==1 && iring==2){
itemvector[istrip].elem33 = 9.118;
itemvector[istrip].elem34 = 3.884;
itemvector[istrip].elem44 = 7.771;
itemvector[istrip].elem35 = 1.8225;
itemvector[istrip].elem45 = 3.7505;
itemvector[istrip].elem55 = 8.597;
itemvector[istrip].elem46 = 1.651;
itemvector[istrip].elem56 = 2.5225;
itemvector[istrip].elem66 = 6.583;
itemvector[istrip].elem57 = 1.5055;
itemvector[istrip].elem67 = 2.733;
itemvector[istrip].elem77 = 6.988;
cnmatrix->matrix[id_layer]=itemvector;
}
if(istation==1 && iring==3){
itemvector[istrip].elem33 = 9.5245;
itemvector[istrip].elem34 = 3.2415;
itemvector[istrip].elem44 = 7.6265;
itemvector[istrip].elem35 = 1.7225;
itemvector[istrip].elem45 = 3.6075;
itemvector[istrip].elem55 = 8.7275;
itemvector[istrip].elem46 = 1.663;
itemvector[istrip].elem56 = 2.592;
itemvector[istrip].elem66 = 7.5685;
itemvector[istrip].elem57 = 1.7905;
itemvector[istrip].elem67 = 2.409;
itemvector[istrip].elem77 = 7.1495;
cnmatrix->matrix[id_layer]=itemvector;
}
if(istation==2 && iring==1){
itemvector[istrip].elem33 = 9.06825;
itemvector[istrip].elem34 = 3.32025;
itemvector[istrip].elem44 = 7.52925;
itemvector[istrip].elem35 = 3.66125;
itemvector[istrip].elem45 = 3.39125;
itemvector[istrip].elem55 = 9.97625;
itemvector[istrip].elem46 = 1.32725;
itemvector[istrip].elem56 = 3.99025;
itemvector[istrip].elem66 = 8.10125;
itemvector[istrip].elem57 = 2.56456;
itemvector[istrip].elem67 = 2.96625;
itemvector[istrip].elem77 = 7.30925;
cnmatrix->matrix[id_layer]=itemvector;
}
if(istation==2 &&iring==2){
itemvector[istrip].elem33 = 16.7442;
itemvector[istrip].elem34 = 7.96925;
itemvector[istrip].elem44 = 14.1643;
itemvector[istrip].elem35 = 4.67975;
itemvector[istrip].elem45 = 8.44075;
itemvector[istrip].elem55 = 17.2243;
itemvector[istrip].elem46 = 3.68575;
itemvector[istrip].elem56 = 7.48825;
itemvector[istrip].elem66 = 14.4902;
itemvector[istrip].elem57 = 4.4482;
itemvector[istrip].elem67 = 6.47875;
itemvector[istrip].elem77 = 14.6733;
cnmatrix->matrix[id_layer]=itemvector;
}
if(istation==3 && iring==1){
itemvector[istrip].elem33 = 9.3495;
itemvector[istrip].elem34 = 3.529;
itemvector[istrip].elem44 = 7.8715;
itemvector[istrip].elem35 = 3.8155;
itemvector[istrip].elem45 = 3.858;
itemvector[istrip].elem55 = 10.8205;
itemvector[istrip].elem46 = 1.8585;
itemvector[istrip].elem56 = 4.445;
itemvector[istrip].elem66 = 8.0175;
itemvector[istrip].elem57 = 3.29479;
itemvector[istrip].elem67 = 3.625;
itemvector[istrip].elem77 = 8.3895;
cnmatrix->matrix[id_layer]=itemvector;
}
if(istation==3 && iring==2){
itemvector[istrip].elem33 = 13.6193;
itemvector[istrip].elem34 = 5.91025;
itemvector[istrip].elem44 = 11.3842;
itemvector[istrip].elem35 = 3.31775;
itemvector[istrip].elem45 = 5.69775;
itemvector[istrip].elem55 = 11.6652;
itemvector[istrip].elem46 = 2.46175;
itemvector[istrip].elem56 = 4.48325;
itemvector[istrip].elem66 = 9.95725;
itemvector[istrip].elem57 = 2.10561;
itemvector[istrip].elem67 = 4.04625;
itemvector[istrip].elem77 = 9.51625;
cnmatrix->matrix[id_layer]=itemvector;
}
if(istation==4 && iring==1){
itemvector[istrip].elem33 = 10.0;
itemvector[istrip].elem34 = 4.0;
itemvector[istrip].elem44 = 10.0;
itemvector[istrip].elem35 = 3.0;
itemvector[istrip].elem45 = 8.0;
itemvector[istrip].elem55 = 10.0;
itemvector[istrip].elem46 = 2.0;
itemvector[istrip].elem56 = 5.0;
itemvector[istrip].elem66 = 10.0;
itemvector[istrip].elem57 = 3.0;
itemvector[istrip].elem67 = 4.0;
itemvector[istrip].elem77 = 10.0;
cnmatrix->matrix[id_layer]=itemvector;
}
}
}
}
}
}
}
//overwrite fakes with old values from DB
int istrip = 0;
std::vector<CSCNoiseMatrix::Item> itemvector;
itemvector.resize(80);
for(unsigned int mystrip=0; mystrip<old_nrlines-1; mystrip++){
if(old_strips[mystrip]==0) istrip = 0;
itemvector[istrip].elem33 = old_elem33[mystrip];
itemvector[istrip].elem34 = old_elem34[mystrip];
itemvector[istrip].elem44 = old_elem44[mystrip];
itemvector[istrip].elem35 = old_elem35[mystrip];
itemvector[istrip].elem45 = old_elem45[mystrip];
itemvector[istrip].elem55 = old_elem55[mystrip];
itemvector[istrip].elem46 = old_elem46[mystrip];
itemvector[istrip].elem56 = old_elem56[mystrip];
itemvector[istrip].elem66 = old_elem66[mystrip];
itemvector[istrip].elem57 = old_elem57[mystrip];
itemvector[istrip].elem67 = old_elem67[mystrip];
itemvector[istrip].elem77 = old_elem77[mystrip];
cnmatrix->matrix[old_cham_id[mystrip]]=itemvector;
istrip++;
}
itemvector.resize(64);
for(unsigned int mystrip=0; mystrip<old_nrlines-1; mystrip++){
if(old_strips[mystrip]==0) istrip = 0;
if(old_cham_id[mystrip] >= 113000 && old_cham_id[mystrip] <= 113999){
itemvector[istrip].elem33 = old_elem33[mystrip];
itemvector[istrip].elem34 = old_elem34[mystrip];
itemvector[istrip].elem44 = old_elem44[mystrip];
itemvector[istrip].elem35 = old_elem35[mystrip];
itemvector[istrip].elem45 = old_elem45[mystrip];
itemvector[istrip].elem55 = old_elem55[mystrip];
itemvector[istrip].elem46 = old_elem46[mystrip];
itemvector[istrip].elem56 = old_elem56[mystrip];
itemvector[istrip].elem66 = old_elem66[mystrip];
itemvector[istrip].elem57 = old_elem57[mystrip];
itemvector[istrip].elem67 = old_elem67[mystrip];
itemvector[istrip].elem77 = old_elem77[mystrip];
cnmatrix->matrix[old_cham_id[mystrip]]=itemvector;
istrip++;
}
}
itemvector.resize(64);
for(unsigned int mystrip=0; mystrip<old_nrlines-1; mystrip++){
if(old_strips[mystrip]==0) istrip = 0;
if(old_cham_id[mystrip] >= 213000 && old_cham_id[mystrip] <= 213999){
itemvector[istrip].elem33 = old_elem33[mystrip];
itemvector[istrip].elem34 = old_elem34[mystrip];
itemvector[istrip].elem44 = old_elem44[mystrip];
itemvector[istrip].elem35 = old_elem35[mystrip];
itemvector[istrip].elem45 = old_elem45[mystrip];
itemvector[istrip].elem55 = old_elem55[mystrip];
itemvector[istrip].elem46 = old_elem46[mystrip];
itemvector[istrip].elem56 = old_elem56[mystrip];
itemvector[istrip].elem66 = old_elem66[mystrip];
itemvector[istrip].elem57 = old_elem57[mystrip];
itemvector[istrip].elem67 = old_elem67[mystrip];
itemvector[istrip].elem77 = old_elem77[mystrip];
cnmatrix->matrix[old_cham_id[mystrip]]=itemvector;
istrip++;
}
}
//overwrite old values with ones from new runs
itemvector.resize(80);
for(unsigned int mystrip=0; mystrip<new_nrlines-1; mystrip++){
if(new_strips[mystrip]==0) istrip = 0;
itemvector[istrip].elem33 = new_elem33[mystrip];
itemvector[istrip].elem34 = new_elem34[mystrip];
itemvector[istrip].elem44 = new_elem44[mystrip];
itemvector[istrip].elem35 = new_elem35[mystrip];
itemvector[istrip].elem45 = new_elem45[mystrip];
itemvector[istrip].elem55 = new_elem55[mystrip];
itemvector[istrip].elem46 = new_elem46[mystrip];
itemvector[istrip].elem56 = new_elem56[mystrip];
itemvector[istrip].elem66 = new_elem66[mystrip];
itemvector[istrip].elem57 = new_elem57[mystrip];
itemvector[istrip].elem67 = new_elem67[mystrip];
itemvector[istrip].elem77 = new_elem77[mystrip];
cnmatrix->matrix[new_cham_id[mystrip]]=itemvector;
istrip++;
}
itemvector.resize(64);
for(unsigned int mystrip=0; mystrip<new_nrlines-1; mystrip++){
if(new_strips[mystrip]==0) istrip = 0;
if(new_cham_id[mystrip] >= 113000 && new_cham_id[mystrip] <= 113999){
itemvector[istrip].elem33 = new_elem33[mystrip];
itemvector[istrip].elem34 = new_elem34[mystrip];
itemvector[istrip].elem44 = new_elem44[mystrip];
itemvector[istrip].elem35 = new_elem35[mystrip];
itemvector[istrip].elem45 = new_elem45[mystrip];
itemvector[istrip].elem55 = new_elem55[mystrip];
itemvector[istrip].elem46 = new_elem46[mystrip];
itemvector[istrip].elem56 = new_elem56[mystrip];
itemvector[istrip].elem66 = new_elem66[mystrip];
itemvector[istrip].elem57 = new_elem57[mystrip];
itemvector[istrip].elem67 = new_elem67[mystrip];
itemvector[istrip].elem77 = new_elem77[mystrip];
cnmatrix->matrix[new_cham_id[mystrip]]=itemvector;
istrip++;
}
}
itemvector.resize(64);
for(unsigned int mystrip=0; mystrip<new_nrlines-1; mystrip++){
if(new_strips[mystrip]==0) istrip = 0;
if(new_cham_id[mystrip] >= 213000 && new_cham_id[mystrip] <= 213999){
itemvector[istrip].elem33 = new_elem33[mystrip];
itemvector[istrip].elem34 = new_elem34[mystrip];
itemvector[istrip].elem44 = new_elem44[mystrip];
itemvector[istrip].elem35 = new_elem35[mystrip];
itemvector[istrip].elem45 = new_elem45[mystrip];
itemvector[istrip].elem55 = new_elem55[mystrip];
itemvector[istrip].elem46 = new_elem46[mystrip];
itemvector[istrip].elem56 = new_elem56[mystrip];
itemvector[istrip].elem66 = new_elem66[mystrip];
itemvector[istrip].elem57 = new_elem57[mystrip];
itemvector[istrip].elem67 = new_elem67[mystrip];
itemvector[istrip].elem77 = new_elem77[mystrip];
cnmatrix->matrix[new_cham_id[mystrip]]=itemvector;
istrip++;
}
}
return cnmatrix;
}
| CSCNoiseMatrixConditions::ReturnType CSCNoiseMatrixConditions::produceNoiseMatrix | ( | const CSCNoiseMatrixRcd & | iRecord | ) |
Definition at line 434 of file CSCNoiseMatrixConditions.cc.
References cnMatrix.
Referenced by CSCNoiseMatrixConditions().
{
// Added by Zhen, need a new object so to not be deleted at exit
CSCNoiseMatrix* mydata=new CSCNoiseMatrix( *cnMatrix );
return mydata;
}
| void CSCNoiseMatrixConditions::setIntervalFor | ( | const edm::eventsetup::EventSetupRecordKey & | , |
| const edm::IOVSyncValue & | , | ||
| edm::ValidityInterval & | oValidity | ||
| ) | [private, virtual] |
Implements edm::EventSetupRecordIntervalFinder.
Definition at line 443 of file CSCNoiseMatrixConditions.cc.
References edm::IOVSyncValue::beginOfTime(), and edm::IOVSyncValue::endOfTime().
{
oValidity = edm::ValidityInterval(edm::IOVSyncValue::beginOfTime(),edm::IOVSyncValue::endOfTime());
}
CSCNoiseMatrix* CSCNoiseMatrixConditions::cnMatrix [private] |
Definition at line 34 of file CSCNoiseMatrixConditions.h.
Referenced by CSCNoiseMatrixConditions(), produceNoiseMatrix(), and ~CSCNoiseMatrixConditions().
1.7.3