CSCCrosstalkConditions Class Reference

#include <CSCCrosstalkConditions.h>

Inheritance diagram for CSCCrosstalkConditions:

Public Types
typedef const CSCcrosstalk *	ReturnType
Public Types inherited from edm::eventsetup::DataProxyProvider
typedef std::vector< std::pair < DataKey, edm::propagate_const < std::shared_ptr< DataProxy > > > >	KeyedProxies
typedef std::vector < EventSetupRecordKey >	Keys
typedef std::map < EventSetupRecordKey, KeyedProxies >	RecordProxies

Public Member Functions
	CSCCrosstalkConditions (const edm::ParameterSet &)
ReturnType	produceCrosstalk (const CSCcrosstalkRcd &)
	~CSCCrosstalkConditions ()
Public Member Functions inherited from edm::ESProducer
	ESProducer ()
virtual	~ESProducer ()
Public Member Functions inherited from edm::ESProxyFactoryProducer
	ESProxyFactoryProducer ()
virtual void	newInterval (const eventsetup::EventSetupRecordKey &iRecordType, const ValidityInterval &iInterval)
	overrides DataProxyProvider method More...
virtual	~ESProxyFactoryProducer ()
Public Member Functions inherited from edm::eventsetup::DataProxyProvider
	DataProxyProvider ()
const ComponentDescription &	description () const
bool	isUsingRecord (const EventSetupRecordKey &) const
const KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey) const
void	resetProxies (const EventSetupRecordKey &iRecordType)
void	resetProxiesIfTransient (const EventSetupRecordKey &iRecordType)
void	setAppendToDataLabel (const edm::ParameterSet &)
void	setDescription (const ComponentDescription &iDescription)
std::set< EventSetupRecordKey >	usingRecords () const
virtual	~DataProxyProvider ()
Public Member Functions inherited from edm::EventSetupRecordIntervalFinder
const eventsetup::ComponentDescription &	descriptionForFinder () const
	EventSetupRecordIntervalFinder ()
std::set < eventsetup::EventSetupRecordKey >	findingForRecords () const
const ValidityInterval &	findIntervalFor (const eventsetup::EventSetupRecordKey &, const IOVSyncValue &)
void	setDescriptionForFinder (const eventsetup::ComponentDescription &iDescription)
virtual	~EventSetupRecordIntervalFinder ()

Static Public Member Functions
static CSCcrosstalk *	prefillCrosstalk ()
Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)

Private Member Functions
void	setIntervalFor (const edm::eventsetup::EventSetupRecordKey &, const edm::IOVSyncValue &, edm::ValidityInterval &)

Private Attributes
CSCcrosstalk *	cnCrosstalk

Additional Inherited Members
Protected Member Functions inherited from edm::ESProducer
template<typename T >
void	setWhatProduced (T *iThis, const es::Label &iLabel=es::Label())
template<typename T >
void	setWhatProduced (T *iThis, const char *iLabel)
template<typename T >
void	setWhatProduced (T *iThis, const std::string &iLabel)
template<typename T , typename TDecorator >
void	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel=es::Label())
template<typename T , typename TReturn , typename TRecord >
void	setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const es::Label &iLabel=es::Label())
template<typename T , typename TReturn , typename TRecord , typename TArg >
void	setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel=es::Label())
Protected Member Functions inherited from edm::ESProxyFactoryProducer
template<class TFactory >
void	registerFactory (std::auto_ptr< TFactory > iFactory, const std::string &iLabel=std::string())
virtual void	registerFactoryWithKey (const eventsetup::EventSetupRecordKey &iRecord, std::auto_ptr< eventsetup::ProxyFactoryBase > &iFactory, const std::string &iLabel=std::string())
virtual void	registerProxies (const eventsetup::EventSetupRecordKey &iRecord, KeyedProxies &aProxyList)
	override DataProxyProvider method More...
Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
void	eraseAll (const EventSetupRecordKey &iRecordKey)
	deletes all the Proxies in aStream More...
void	invalidateProxies (const EventSetupRecordKey &iRecordKey)
template<class T >
void	usingRecord ()
void	usingRecordWithKey (const EventSetupRecordKey &)
Protected Member Functions inherited from edm::EventSetupRecordIntervalFinder
template<class T >
void	findingRecord ()
void	findingRecordWithKey (const eventsetup::EventSetupRecordKey &)

Detailed Description

Definition at line 19 of file CSCCrosstalkConditions.h.

Member Typedef Documentation

typedef const CSCcrosstalk* CSCCrosstalkConditions::ReturnType

Definition at line 27 of file CSCCrosstalkConditions.h.

Constructor & Destructor Documentation

CSCCrosstalkConditions::CSCCrosstalkConditions

(

const edm::ParameterSet &

iConfig

)

Definition at line 253 of file CSCCrosstalkConditions.cc.

References cnCrosstalk, prefillCrosstalk(), produceCrosstalk(), and edm::ESProducer::setWhatProduced().

{
  //the following line is needed to tell the framework what
  // data is being produced
  cnCrosstalk = prefillCrosstalk();
  // added by Zhen (changed since 1_2_0)
  setWhatProduced(this,&CSCCrosstalkConditions::produceCrosstalk);
  findingRecord<CSCcrosstalkRcd>();
  //now do what ever other initialization is needed
}

CSCCrosstalkConditions::~CSCCrosstalkConditions

(

)

Definition at line 265 of file CSCCrosstalkConditions.cc.

References cnCrosstalk.

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

Member Function Documentation

CSCcrosstalk * CSCCrosstalkConditions::prefillCrosstalk ( )

static

Definition at line 6 of file CSCCrosstalkConditions.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, CSCcrosstalk::crosstalk, CSCDetId, cmsRelvalreport::exit, recoMuon::in, CSCDetId::maxChamberId(), CSCDetId::maxEndcapId(), CSCDetId::maxLayerId(), CSCDetId::maxRingId(), CSCDetId::maxStationId(), timingPdfMaker::mean, min(), CSCDetId::minChamberId(), CSCDetId::minEndcapId(), CSCDetId::minLayerId(), CSCDetId::minRingId(), CSCDetId::minStationId(), rand(), and fileCollector::seed.

Referenced by CSCCrosstalkConditions().

                                                        {

  float mean,min,minchi;
  int seed;
  int old_chamber_id,old_strip,new_chamber_id,new_strip;
  float old_slope_right,old_slope_left,old_intercept_right;
  float old_intercept_left, old_chi2_right,old_chi2_left;
  std::vector<int> old_cham_id;
  std::vector<int> old_strips;
  std::vector<float> old_slope_r;
  std::vector<float> old_intercept_r;
  std::vector<float> old_chi2_r;
  std::vector<float> old_slope_l;
  std::vector<float> old_intercept_l;
  std::vector<float> old_chi2_l;
  float new_slope_right,new_slope_left,new_intercept_right;
  float new_intercept_left, new_chi2_right,new_chi2_left;
  std::vector<int> new_cham_id;
  std::vector<int> new_strips;
  std::vector<float> new_slope_r;
  std::vector<float> new_intercept_r;
  std::vector<float> new_chi2_r;
  std::vector<float> new_slope_l;
  std::vector<float> new_intercept_l;
  std::vector<float> new_chi2_l;
  


  const CSCDetId& detId = CSCDetId();
  CSCcrosstalk * cncrosstalk = new CSCcrosstalk();
  
  int max_istrip,id_layer,max_ring,max_cham;
  unsigned int old_nrlines=0;
  unsigned int new_nrlines=0;
  seed = 10000; 
  srand(seed);
  mean=-0.0009, min=0.035, minchi=1.5;
 
  //endcap=1 to 2,station=1 to 4, ring=1 to 4,chamber=1 to 36,layer=1 to 6 
  std::ifstream olddata; 
  olddata.open("old_xtalk.dat",std::ios::in); 
  if(!olddata) {
    std::cerr <<"Error: old_xtalk.dat -> no such file!"<< std::endl;
    exit(1);
  }
  
  while (!olddata.eof() ) { 
    olddata >> old_chamber_id >> old_strip >> old_slope_right >> old_intercept_right >> old_chi2_right >> old_slope_left >> old_intercept_left >> old_chi2_left ; 
    old_cham_id.push_back(old_chamber_id);
    old_strips.push_back(old_strip);
    old_slope_r.push_back(old_slope_right);
    old_slope_l.push_back(old_slope_left);
    old_intercept_r.push_back(old_intercept_right);
    old_intercept_l.push_back(old_intercept_left);
    old_chi2_r.push_back(old_chi2_right);
    old_chi2_l.push_back(old_chi2_left);
    old_nrlines++;
  }
  olddata.close();

  std::ifstream newdata;
  newdata.open("new_xtalk.txt",std::ios::in); 
  if(!newdata) {
    std::cerr <<"Error: new_xtalk.dat -> no such file!"<< std::endl;
    exit(1);
  }
  
  while (!newdata.eof() ) { 
    newdata >> new_chamber_id >> new_strip >> new_slope_right >> new_intercept_right >> new_chi2_right >> new_slope_left >> new_intercept_left >> new_chi2_left ; 
    new_cham_id.push_back(new_chamber_id);
    new_strips.push_back(new_strip);
    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_chi2_r.push_back(new_chi2_right);
    new_chi2_l.push_back(new_chi2_left);
    new_nrlines++;
  }
  newdata.close();

  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;
    //station 1 ring 3 has 64 strips per layer instead of 80(minus & plus side!!!)

    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<CSCcrosstalk::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++){
          //create fake values
          itemvector[istrip].xtalk_slope_right=-((double)rand()/((double)(RAND_MAX)+(double)(1)))/10000+mean;
          itemvector[istrip].xtalk_intercept_right=((double)rand()/((double)(RAND_MAX)+(double)(1)))/100+min;
          itemvector[istrip].xtalk_chi2_right=((double)rand()/((double)(RAND_MAX)+(double)(1)))+minchi;
          itemvector[istrip].xtalk_slope_left=-((double)rand()/((double)(RAND_MAX)+(double)(1)))/10000+mean;
          itemvector[istrip].xtalk_intercept_left=((double)rand()/((double)(RAND_MAX)+(double)(1)))/100+min;
          itemvector[istrip].xtalk_chi2_left=((double)rand()/((double)(RAND_MAX)+(double)(1)))+minchi;
          cncrosstalk->crosstalk[id_layer]=itemvector;

          if(istrip==0){
        itemvector[istrip].xtalk_slope_right=-((double)rand()/((double)(RAND_MAX)+(double)(1)))/10000+mean;
        itemvector[istrip].xtalk_intercept_right=((double)rand()/((double)(RAND_MAX)+(double)(1)))/100+min;
        itemvector[istrip].xtalk_chi2_right=((double)rand()/((double)(RAND_MAX)+(double)(1)))+minchi;
        itemvector[istrip].xtalk_slope_left=0.0;
        itemvector[istrip].xtalk_intercept_left=0.0;
        itemvector[istrip].xtalk_chi2_left=0.0;
        cncrosstalk->crosstalk[id_layer]=itemvector;
          }
          
          if(istrip==79){
        itemvector[istrip].xtalk_slope_right=0.0;
        itemvector[istrip].xtalk_intercept_right=0.0;
        itemvector[istrip].xtalk_chi2_right=0.0;
        itemvector[istrip].xtalk_slope_left=-((double)rand()/((double)(RAND_MAX)+(double)(1)))/10000+mean;
        itemvector[istrip].xtalk_intercept_left=((double)rand()/((double)(RAND_MAX)+(double)(1)))/100+min;
        itemvector[istrip].xtalk_chi2_left=((double)rand()/((double)(RAND_MAX)+(double)(1)))+minchi;
        cncrosstalk->crosstalk[id_layer]=itemvector;
          }
        }
      }
    }
      }
    }
    
    //overwrite fakes with old values from DB
    int istrip = 0;
    std::vector<CSCcrosstalk::Item> itemvector;
    itemvector.resize(80);
    
    for(unsigned int mystrip=0; mystrip<old_nrlines-1; mystrip++){
      if(old_strips[mystrip]==0) istrip = 0;
      itemvector[istrip].xtalk_slope_right=old_slope_r[mystrip];
      itemvector[istrip].xtalk_intercept_right=old_intercept_r[mystrip]; 
      itemvector[istrip].xtalk_chi2_right=old_chi2_r[mystrip];
      itemvector[istrip].xtalk_slope_left=old_slope_l[mystrip];  
      itemvector[istrip].xtalk_intercept_left=old_intercept_l[mystrip];  
      itemvector[istrip].xtalk_chi2_left=old_chi2_l[mystrip];
      cncrosstalk->crosstalk[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].xtalk_slope_right=old_slope_r[mystrip];
    itemvector[istrip].xtalk_intercept_right=old_intercept_r[mystrip]; 
    itemvector[istrip].xtalk_chi2_right=old_chi2_r[mystrip];
    itemvector[istrip].xtalk_slope_left=old_slope_l[mystrip];  
    itemvector[istrip].xtalk_intercept_left=old_intercept_l[mystrip];  
    itemvector[istrip].xtalk_chi2_left=old_chi2_l[mystrip];
    cncrosstalk->crosstalk[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].xtalk_slope_right=old_slope_r[mystrip];
    itemvector[istrip].xtalk_intercept_right=old_intercept_r[mystrip]; 
    itemvector[istrip].xtalk_chi2_right=old_chi2_r[mystrip];
    itemvector[istrip].xtalk_slope_left=old_slope_l[mystrip];  
    itemvector[istrip].xtalk_intercept_left=old_intercept_l[mystrip];  
    itemvector[istrip].xtalk_chi2_left=old_chi2_l[mystrip];
    cncrosstalk->crosstalk[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].xtalk_slope_right=new_slope_r[mystrip];
      itemvector[istrip].xtalk_intercept_right=new_intercept_r[mystrip]; 
      itemvector[istrip].xtalk_chi2_right=new_chi2_r[mystrip];
      itemvector[istrip].xtalk_slope_left=new_slope_l[mystrip];  
      itemvector[istrip].xtalk_intercept_left=new_intercept_l[mystrip];  
      itemvector[istrip].xtalk_chi2_left=new_chi2_l[mystrip];
      cncrosstalk->crosstalk[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].xtalk_slope_right=new_slope_r[mystrip];
    itemvector[istrip].xtalk_intercept_right=new_intercept_r[mystrip]; 
    itemvector[istrip].xtalk_chi2_right=new_chi2_r[mystrip];
    itemvector[istrip].xtalk_slope_left=new_slope_l[mystrip];  
    itemvector[istrip].xtalk_intercept_left=new_intercept_l[mystrip];  
    itemvector[istrip].xtalk_chi2_left=new_chi2_l[mystrip];
    cncrosstalk->crosstalk[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].xtalk_slope_right=new_slope_r[mystrip];
    itemvector[istrip].xtalk_intercept_right=new_intercept_r[mystrip]; 
    itemvector[istrip].xtalk_chi2_right=new_chi2_r[mystrip];
    itemvector[istrip].xtalk_slope_left=new_slope_l[mystrip];  
    itemvector[istrip].xtalk_intercept_left=new_intercept_l[mystrip];  
    itemvector[istrip].xtalk_chi2_left=new_chi2_l[mystrip];
    cncrosstalk->crosstalk[new_cham_id[mystrip]]=itemvector;
    istrip++;
      }
    } 
  }

  return  cncrosstalk;

}

CSCCrosstalkConditions::ReturnType CSCCrosstalkConditions::produceCrosstalk

(

const CSCcrosstalkRcd &

iRecord

)

Definition at line 280 of file CSCCrosstalkConditions.cc.

References cnCrosstalk.

Referenced by CSCCrosstalkConditions().

{
  // Added by Zhen, need a new object so to not be deleted at exit
  CSCcrosstalk* mydata=new CSCcrosstalk( *cnCrosstalk );
  
  return mydata;
  
}

void CSCCrosstalkConditions::setIntervalFor ( const edm::eventsetup::EventSetupRecordKey &  , const edm::IOVSyncValue &  , edm::ValidityInterval &  oValidity  )

privatevirtual

Implements edm::EventSetupRecordIntervalFinder.

Definition at line 289 of file CSCCrosstalkConditions.cc.

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

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

Member Data Documentation

CSCcrosstalk* CSCCrosstalkConditions::cnCrosstalk

private

Definition at line 34 of file CSCCrosstalkConditions.h.

Referenced by CSCCrosstalkConditions(), produceCrosstalk(), and ~CSCCrosstalkConditions().