CSCCrosstalkDBConditions Class Reference

#include <CSCCrosstalkDBConditions.h>

Inheritance diagram for CSCCrosstalkDBConditions:

List of all members.

Public Types
typedef const CSCDBCrosstalk *	ReturnType
Public Member Functions
	CSCCrosstalkDBConditions (const edm::ParameterSet &)
ReturnType	produceDBCrosstalk (const CSCDBCrosstalkRcd &)
	~CSCCrosstalkDBConditions ()
Static Public Member Functions
static CSCDBCrosstalk *	prefillDBCrosstalk ()
Private Member Functions
void	setIntervalFor (const edm::eventsetup::EventSetupRecordKey &, const edm::IOVSyncValue &, edm::ValidityInterval &)
Private Attributes
CSCDBCrosstalk *	cndbCrosstalk

---

Detailed Description

Definition at line 20 of file CSCCrosstalkDBConditions.h.

---

Member Typedef Documentation

typedef const CSCDBCrosstalk* CSCCrosstalkDBConditions::ReturnType

Definition at line 28 of file CSCCrosstalkDBConditions.h.

---

Constructor & Destructor Documentation

CSCCrosstalkDBConditions::CSCCrosstalkDBConditions

(

const edm::ParameterSet &

iConfig

)

Definition at line 10 of file CSCCrosstalkDBConditions.cc.

References cndbCrosstalk, prefillDBCrosstalk(), produceDBCrosstalk(), and edm::ESProducer::setWhatProduced().

{
  //the following line is needed to tell the framework what
  // data is being produced
  cndbCrosstalk = prefillDBCrosstalk();
  // added by Zhen (changed since 1_2_0)
  setWhatProduced(this,&CSCCrosstalkDBConditions::produceDBCrosstalk);
  findingRecord<CSCDBCrosstalkRcd>();
  //now do what ever other initialization is needed
}

CSCCrosstalkDBConditions::~CSCCrosstalkDBConditions

(

)

Definition at line 22 of file CSCCrosstalkDBConditions.cc.

References cndbCrosstalk.

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

---

Member Function Documentation

CSCDBCrosstalk * CSCCrosstalkDBConditions::prefillDBCrosstalk

(

)

[inline, static]

Definition at line 46 of file CSCCrosstalkDBConditions.h.

References dtNoiseDBValidation_cfg::cerr, CSCDBCrosstalk::crosstalk, cmsRelvalreport::exit, CSCDBCrosstalk::factor_intercept, CSCDBCrosstalk::factor_slope, i, recoMuon::in, and gen::k.

Referenced by CSCCrosstalkDBConditions().

{
  // const int MAX_SIZE = 273024; //for ME1a unganged
  const int MAX_SIZE = 252288;
  const int SLOPE_FACTOR=10000000;
  const int INTERCEPT_FACTOR=100000;
  const int MAX_SHORT= 32767;
  CSCDBCrosstalk * cndbcrosstalk = new CSCDBCrosstalk();

  int db_index,new_index;
  float db_slope_right,db_slope_left,db_intercept_right;
  float db_intercept_left;
  std::vector<int> db_index_id;
  std::vector<float> db_slope_r;
  std::vector<float> db_intercept_r;
  std::vector<float> db_slope_l;
  std::vector<float> db_intercept_l;
  float new_slope_right,new_slope_left,new_intercept_right;
  float new_intercept_left;
   std::vector<int> new_index_id;
   std::vector<float> new_slope_r;
  std::vector<float> new_intercept_r;
  std::vector<float> new_slope_l;
  std::vector<float> new_intercept_l;
    
  int counter;
  int db_nrlines=0;
  int new_nrlines=0;
    
  std::ifstream dbdata; 
  dbdata.open("old_dbxtalk.dat",std::ios::in); 
  if(!dbdata) {
    std::cerr <<"Error: old_dbxtalk.dat -> no such file!"<< std::endl;
    exit(1);
  }
  
  while (!dbdata.eof() ) { 
    dbdata >> db_index >> db_slope_right >> db_intercept_right >> db_slope_left >> db_intercept_left ;
    db_index_id.push_back(db_index);
    db_slope_r.push_back(db_slope_right);
    db_slope_l.push_back(db_slope_left);
    db_intercept_r.push_back(db_intercept_right);
    db_intercept_l.push_back(db_intercept_left);
    db_nrlines++;
  }
  dbdata.close();

  std::ifstream newdata;
  newdata.open("xtalk.dat",std::ios::in); 
  if(!newdata) {
    std::cerr <<"Error: xtalk.dat -> no such file!"<< std::endl;
    exit(1);
  }
  
  while (!newdata.eof() ) { 
    newdata >> new_index >> new_slope_right >> new_intercept_right >> new_slope_left >> new_intercept_left; 
    new_index_id.push_back(new_index);
    new_slope_r.push_back(new_slope_right);
    new_slope_l.push_back(new_slope_left);
    new_intercept_r.push_back(new_intercept_right);
    new_intercept_l.push_back(new_intercept_left);
    new_nrlines++;
  }
  newdata.close();
  
  CSCDBCrosstalk::CrosstalkContainer & itemvector = cndbcrosstalk->crosstalk;
  itemvector.resize(MAX_SIZE);
  cndbcrosstalk->factor_slope = int (SLOPE_FACTOR);
  cndbcrosstalk->factor_intercept = int (INTERCEPT_FACTOR);

  for(int i=0; i<MAX_SIZE;++i){
    itemvector[i].xtalk_slope_right=(short int) (db_slope_r[i]*SLOPE_FACTOR+0.5);
    itemvector[i].xtalk_intercept_right= (short int) (db_intercept_r[i]*INTERCEPT_FACTOR+0.5); 
    itemvector[i].xtalk_slope_left= (short int) (db_slope_l[i]*SLOPE_FACTOR+0.5);  
    itemvector[i].xtalk_intercept_left= (short int) (db_intercept_l[i]*INTERCEPT_FACTOR+0.5);  
  }

  for(int i=0; i<MAX_SIZE;++i){
    counter=db_index_id[i];  
    for (unsigned int k=0;k<new_index_id.size()-1;k++){
      if(counter==new_index_id[k]){
        if ((short int) (fabs(new_slope_r[k]*SLOPE_FACTOR+0.5))<MAX_SHORT)         itemvector[counter].xtalk_slope_right= int (new_slope_r[k]*SLOPE_FACTOR+0.5);
        if ((short int) (fabs(new_intercept_r[k]*INTERCEPT_FACTOR+0.5))<MAX_SHORT) itemvector[counter].xtalk_intercept_right= int (new_intercept_r[k]*INTERCEPT_FACTOR+0.5); 
        if ((short int) (fabs(new_slope_l[k]*SLOPE_FACTOR+0.5))<MAX_SHORT)         itemvector[counter].xtalk_slope_left= int (new_slope_l[k]*SLOPE_FACTOR+0.5);  
        if ((short int) (fabs(new_intercept_l[k]*INTERCEPT_FACTOR+0.5))<MAX_SHORT) itemvector[counter].xtalk_intercept_left= int (new_intercept_l[k]*INTERCEPT_FACTOR+0.5);  
        itemvector[i] = itemvector[counter];
        //std::cout<<" counter "<<counter <<" dbindex "<<new_index_id[k]<<" dbslope " <<db_slope_r[k]<<" new slope "<<new_slope_r[k]<<std::endl;
      }  
    }

     if (counter>223968){
      itemvector[counter].xtalk_slope_right = int (db_slope_r[i]);
      itemvector[counter].xtalk_slope_left = int (db_slope_l[i]);
      itemvector[counter].xtalk_intercept_right = int (db_intercept_r[i]);
      itemvector[counter].xtalk_intercept_left = int (db_intercept_l[i]);
     }
  }
  
  return cndbcrosstalk;

}

CSCCrosstalkDBConditions::ReturnType CSCCrosstalkDBConditions::produceDBCrosstalk

(

const CSCDBCrosstalkRcd &

iRecord

)

Definition at line 37 of file CSCCrosstalkDBConditions.cc.

References cndbCrosstalk.

Referenced by CSCCrosstalkDBConditions().

{
  //need a new object so to not be deleted at exit
  CSCDBCrosstalk* mydata=new CSCDBCrosstalk( *cndbCrosstalk );
  return mydata;
  
}

void CSCCrosstalkDBConditions::setIntervalFor	(	const edm::eventsetup::EventSetupRecordKey &	,
		const edm::IOVSyncValue &	,
		edm::ValidityInterval &	oValidity
	)		[private, virtual]

Implements edm::EventSetupRecordIntervalFinder.

Definition at line 45 of file CSCCrosstalkDBConditions.cc.

References edm::IOVSyncValue::beginOfTime(), and edm::IOVSyncValue::endOfTime().

 {
 oValidity = edm::ValidityInterval(edm::IOVSyncValue::beginOfTime(),edm::IOVSyncValue::endOfTime());
 
 }

---

Member Data Documentation

CSCDBCrosstalk* CSCCrosstalkDBConditions::cndbCrosstalk [private]

Definition at line 37 of file CSCCrosstalkDBConditions.h.

Referenced by CSCCrosstalkDBConditions(), produceDBCrosstalk(), and ~CSCCrosstalkDBConditions().