#include <CSCChamberTimeCorrectionsValues.h>

Inheritance diagram for CSCChamberTimeCorrectionsValues:

Public Types
typedef std::unique_ptr< CSCChamberTimeCorrections >	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
	CSCChamberTimeCorrectionsValues (const edm::ParameterSet &)

ReturnType	produceChamberTimeCorrections (const CSCChamberTimeCorrectionsRcd &)

	~CSCChamberTimeCorrectionsValues () override

Public Member Functions inherited from edm::ESProducer
	ESProducer ()

	~ESProducer ()(false) override

Public Member Functions inherited from edm::ESProxyFactoryProducer
	ESProxyFactoryProducer ()

void	newInterval (const eventsetup::EventSetupRecordKey &iRecordType, const ValidityInterval &iInterval) override
	overrides DataProxyProvider method More...

	~ESProxyFactoryProducer () noexcept(false) override

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 () noexcept(false)

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 () noexcept(false)

Static Public Member Functions
static CSCChamberTimeCorrections *	prefill (bool isMC, float ME11offset, float nonME11offset)

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

Private Attributes
bool	isForMC

float	ME11offsetData

float	ME11offsetMC

float	nonME11offsetData

float	nonME11offsetMC

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::unique_ptr< TFactory > iFactory, const std::string &iLabel=std::string())

virtual void	registerFactoryWithKey (const eventsetup::EventSetupRecordKey &iRecord, std::unique_ptr< eventsetup::ProxyFactoryBase > iFactory, const std::string &iLabel=std::string())

void	registerProxies (const eventsetup::EventSetupRecordKey &iRecord, KeyedProxies &aProxyList) override
	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 23 of file CSCChamberTimeCorrectionsValues.h.

Member Typedef Documentation

typedef std::unique_ptr<CSCChamberTimeCorrections> CSCChamberTimeCorrectionsValues::ReturnType

Definition at line 28 of file CSCChamberTimeCorrectionsValues.h.

Constructor & Destructor Documentation

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

Definition at line 10 of file CSCChamberTimeCorrectionsValues.cc.

References edm::ParameterSet::getUntrackedParameter(), isForMC, ME11offsetData, ME11offsetMC, nonME11offsetData, nonME11offsetMC, produceChamberTimeCorrections(), and edm::ESProducer::setWhatProduced().

 {
   //the following line is needed to tell the framework what
   // data is being produced
   isForMC = iConfig.getUntrackedParameter<bool>("isForMC",true);
   ME11offsetMC = 184;
   ME11offsetData = 205;
   nonME11offsetMC = 174;
   nonME11offsetData = 216;
   setWhatProduced(this,&CSCChamberTimeCorrectionsValues::produceChamberTimeCorrections);
   findingRecord<CSCChamberTimeCorrectionsRcd>();
   //now do what ever other initialization is needed
 }

CSCChamberTimeCorrectionsValues::~CSCChamberTimeCorrectionsValues ( )

override

Definition at line 25 of file CSCChamberTimeCorrectionsValues.cc.

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

Member Function Documentation

CSCChamberTimeCorrections * CSCChamberTimeCorrectionsValues::prefill	(	bool	isMC,
		float	ME11offset,
		float	nonME11offset
	)

inlinestatic

Definition at line 52 of file CSCChamberTimeCorrectionsValues.h.

References csccableread::cable_read(), relativeConstraints::chamber, CSCChamberTimeCorrections::chamberCorrections, CSCIndexer::chamberIndex(), trackerTree::check(), corr, KineDebug3::count(), CSCChamberTimeCorrections::factor_precision, groupFilesInBlocks::fin, objects.autophobj::float, mps_fire::i, createfilelist::int, diffTwoXMLs::label, PFRecoTauDiscriminationByIsolation_cfi::offset, AlCaHLTBitMon_QueryRunRegistry::string, and groupFilesInBlocks::temp.

Referenced by produceChamberTimeCorrections().

 {
   if (isMC)
     printf("\n Generating fake DB constants for MC\n");
   else 
     printf("\n Getting chamber corrections from the cable data base and possibly other files \n");
 
   const int FACTOR=100;
   const int MAX_SIZE = 540;
   //const int MAX_SHORT= 32767;
 
   CSCChamberTimeCorrections * chamberObj = new CSCChamberTimeCorrections();
 
   int i; //i - chamber index.
   int count=0;
   std::string chamber_label, cfeb_rev, alct_rev;
   float cfeb_length=0, alct_length=0, cfeb_tmb_skew_delay=0, cfeb_timing_corr=0;
 
   // Only the first 481 chambers have interesting cable lengths at present
   // The rest of the chambers will be filled with zeros
   chamberObj->factor_precision = FACTOR;
 
   chamberObj->chamberCorrections.resize(MAX_SIZE);
   // fill the database with dummy values
   for(i=1;i<=MAX_SIZE;++i){
     chamberObj->chamberCorrections[i-1].cfeb_length=0;
     chamberObj->chamberCorrections[i-1].cfeb_rev='X';
     chamberObj->chamberCorrections[i-1].alct_length=0;
     chamberObj->chamberCorrections[i-1].alct_rev='X';
     chamberObj->chamberCorrections[i-1].cfeb_tmb_skew_delay=0;
     chamberObj->chamberCorrections[i-1].cfeb_timing_corr=0;
     chamberObj->chamberCorrections[i-1].cfeb_cable_delay=0;
     chamberObj->chamberCorrections[i-1].anode_bx_offset=0;
   }
 
 
   // for MC there will is a different correction for each chamber type
   if (isMC){
     float OffsetByType;
     float anodeOffset;
     for(i=1;i<=MAX_SIZE;++i){   
       if (i<= 36 || (i>= 235 && i<=270))       { OffsetByType=172.;  anodeOffset=6.18; }// 1/1
       else if (i<= 72 || (i>= 271 && i<=306))  { OffsetByType=168.;  anodeOffset=6.22; }// 1/2
       else if (i<= 108 || (i>= 307 && i<=342)) { OffsetByType=177.;  anodeOffset=6.19; }// 1/3
       else if (i<= 126 || (i>= 343 && i<=360)) { OffsetByType=171.;  anodeOffset=6.25; }// 2/1
       else if (i<= 162 || (i>= 361 && i<=396)) { OffsetByType=175.;  anodeOffset=6.21; }// 2/2
       else if (i<= 180 || (i>= 397 && i<=414)) { OffsetByType=171.;  anodeOffset=6.25; }// 3/1
       else if (i<= 216 || (i>= 415 && i<=450)) { OffsetByType=175.;  anodeOffset=6.20; }// 3/2
       else if (i<= 234 || (i>= 451 && i<=468)) { OffsetByType=172.;  anodeOffset=6.19; }// 4/1
       else {OffsetByType=175; anodeOffset=6.21; }// 4/2
 
     chamberObj->chamberCorrections[i-1].cfeb_timing_corr=(short int)(-1*OffsetByType*FACTOR+0.5*(-1*OffsetByType>=0)-0.5*(-1*OffsetByType<0));
     chamberObj->chamberCorrections[i-1].anode_bx_offset =(short int)(anodeOffset *FACTOR+0.5*(anodeOffset >=0)-0.5*(anodeOffset <0));
     }
 
     return chamberObj;
   }
 
   // ***************************************************************************
   // Everything below this point is for setting the chamber corrections for data
   // ***************************************************************************
 
   csccableread cable;
   for(i=1;i<=MAX_SIZE;++i){
     // the anode bx offset is 8.15 bx for chambers in 2/1, 3/1, and 4/1 
     // and 8.18 bx for all other chambers for early runs (8.20 for runs> 149357)
     float anodeOffset;
     if (i<= 36 || (i>= 235 && i<=270))       {   anodeOffset=8.20; }// 1/1
     else if (i<= 72 || (i>= 271 && i<=306))  {   anodeOffset=8.20; }// 1/2
     else if (i<= 108 || (i>= 307 && i<=342)) {   anodeOffset=8.20; }// 1/3
     else if (i<= 126 || (i>= 343 && i<=360)) {   anodeOffset=8.15; }// 2/1
     else if (i<= 162 || (i>= 361 && i<=396)) {   anodeOffset=8.20; }// 2/2
     else if (i<= 180 || (i>= 397 && i<=414)) {   anodeOffset=8.15; }// 3/1
     else if (i<= 216 || (i>= 415 && i<=450)) {   anodeOffset=8.20; }// 3/2
     else if (i<= 234 || (i>= 451 && i<=468)) {   anodeOffset=8.15; }// 4/1
     else {anodeOffset=8.20; }// 4/2
 
     // for data we will read in from Igor's database
     cable.cable_read(i, &chamber_label, &cfeb_length, &cfeb_rev, &alct_length,
               &alct_rev, &cfeb_tmb_skew_delay, &cfeb_timing_corr);
     // If the read of the cable database is useful (if there is information for the chamber there)
     // re-enter the information the cable object
     if(!chamber_label.empty() && !(cfeb_length==0)){
       chamberObj->chamberCorrections[i-1].cfeb_length=(short int)(cfeb_length*FACTOR+0.5);
       chamberObj->chamberCorrections[i-1].cfeb_rev=cfeb_rev[0];
       chamberObj->chamberCorrections[i-1].alct_length=(short int)(alct_length*FACTOR+0.5);
       chamberObj->chamberCorrections[i-1].alct_rev=alct_rev[0];
       chamberObj->chamberCorrections[i-1].cfeb_tmb_skew_delay=(short int)(cfeb_tmb_skew_delay*FACTOR+0.5);
       chamberObj->chamberCorrections[i-1].cfeb_timing_corr=(short int)(cfeb_timing_corr*FACTOR+0.5);
       chamberObj->chamberCorrections[i-1].cfeb_cable_delay=0;
       chamberObj->chamberCorrections[i-1].anode_bx_offset=(short int)(anodeOffset*FACTOR+0.5);
     }
     count=count+1;
   }
 
   //Read in the changes you want to make in the extra chamber variable cfeb_timing_corr
   FILE *fin = fopen("/afs/cern.ch/user/d/deisher/public/TimingCorrections2009/ttcrx_delay_effects_23April_2010.txt","r");
   int chamber;
   float corr;
   while (!feof(fin)){
     //note space at end of format string to convert last \n
     int check = fscanf(fin,"%d %f \n",&chamber,&corr);
     if (check != 2){
       printf("cfeb timing corr file has an unexpected format \n");
       assert(0);  
     }
     //printf("chamber %d corr %f \n",chamber,corr);
     chamberObj->chamberCorrections[chamber-1].cfeb_timing_corr= (short int)(corr*FACTOR+0.5*(corr>=0)-0.5*(corr<0));
   }
   fclose(fin);  
   
   // Remove the offsets inherent to ME11 and non ME11 chambers
   for(i=1;i<=MAX_SIZE;++i){
     float temp= float(chamberObj->chamberCorrections[i-1].cfeb_timing_corr)/FACTOR;
     if (i<= 36 || (i>= 235 && i<=270))
       chamberObj->chamberCorrections[i-1].cfeb_timing_corr=(short int)((temp-1*ME11offset)*FACTOR+0.5*(temp>=ME11offset)-0.5*(temp<ME11offset));
     else
       chamberObj->chamberCorrections[i-1].cfeb_timing_corr=(short int)((temp-1*nonME11offset)*FACTOR+0.5*(temp>=nonME11offset)-0.5*(temp<nonME11offset));
   }
 
   //Read in the cfeb_cable_delay values (0 or 1) and don't use a precision correction factor 
   FILE *fdelay = fopen("/afs/cern.ch/user/d/deisher/public/TimingCorrections2009/cfeb_cable_delay_20100423_both.txt","r");
   char label[1024];
   int delay;
   CSCIndexer indexer;
   while (!feof(fdelay)){
     //note space at end of format string to convert last \n
     int check = fscanf(fdelay,"%1024s %d \n",label,&delay);
     if (check != 2){
       printf("cfeb cable delay file has an unexpected format \n");
       assert(0);  
     }
     int chamberSerial = 0;
     int c_endcap = (label[2]=='+' ? 1 : 2);
     int c_station = atoi(&label[3]);
     int c_ring = atoi(&label[5]);
     if (c_station ==1 && c_ring ==4)
       c_ring =1;
     int c_chamber = (label[7]=='0' ? atoi(&label[8]) : atoi(&label[7]));
     chamberSerial = indexer.chamberIndex(c_endcap, c_station, c_ring, c_chamber);
     //printf("chamberLabel %s (%d %d %d %d) chamberSerial %d delay %d \n",label,c_endcap,c_station, c_ring, c_chamber, chamberSerial,delay);
     chamberObj->chamberCorrections[chamberSerial-1].cfeb_cable_delay= (short int)delay;
   }
   fclose(fdelay);  
 
   //Read in a 2nd order correction for chamber offsets derived from data 
   FILE *foffset = fopen("/afs/cern.ch/user/d/deisher/public/TimingCorrections2009/offset_26July2010_codeOverhaul_slope012.txt","r");
   float offset;
   int iE,iS,iR,iC;
   while (!feof(foffset)){
     //note space at end of format string to convert last \n
     int check = fscanf(foffset,"%d %d %d %d %f \n",&iE,&iS,&iR,&iC,&offset);
     if (check != 5){
       printf("offset file has an unexpected format \n");
       assert(0);  
     }
     int chamberSerial = 0;
     if (iS ==1 && iR ==4)
       iR =1;
     chamberSerial = indexer.chamberIndex(iE,iS,iR,iC);
     //printf("chamberLabel %s (%d %d %d %d) chamberSerial %d delay %d \n",label,c_endcap,c_station, c_ring, c_chamber, chamberSerial,delay);
     float temp= float(chamberObj->chamberCorrections[chamberSerial-1].cfeb_timing_corr)/FACTOR;
     chamberObj->chamberCorrections[chamberSerial-1].cfeb_timing_corr= (short int)((temp-offset)*FACTOR+0.5*(temp>=offset)-0.5*(temp<offset));
     printf("Serial %d old corr  %f change %f newcorr %f \n",chamberSerial,temp,offset,(float)chamberObj->chamberCorrections[chamberSerial-1].cfeb_timing_corr/FACTOR);
   }
   fclose(foffset);  
   
   //Read in a 3rd order correction for chamber offsets derived from data 
   FILE *foffsetAgain = fopen("/afs/cern.ch/user/d/deisher/public/TimingCorrections2009/CathodeTimingCorrection_DB_12082010.txt","r");
   while (!feof(foffsetAgain)){
     //note space at end of format string to convert last \n
     int check = fscanf(foffsetAgain,"%d %d %d %d %f \n",&iE,&iS,&iR,&iC,&offset);
     if (check != 5){
       printf("offsetAgain file has an unexpected format \n");
       assert(0);  
     }
     int chamberSerial = 0;
     if (iS ==1 && iR ==4)
       iR =1;
     chamberSerial = indexer.chamberIndex(iE,iS,iR,iC);
     //printf("chamberLabel %s (%d %d %d %d) chamberSerial %d delay %d \n",label,c_endcap,c_station, c_ring, c_chamber, chamberSerial,delay);
     float temp= float(chamberObj->chamberCorrections[chamberSerial-1].cfeb_timing_corr)/FACTOR;
     chamberObj->chamberCorrections[chamberSerial-1].cfeb_timing_corr= (short int)((temp-offset)*FACTOR+0.5*(temp>=offset)-0.5*(temp<offset));
     printf("Serial %d old corr  %f change %f newcorr %f \n",chamberSerial,temp,offset,(float)chamberObj->chamberCorrections[chamberSerial-1].cfeb_timing_corr/FACTOR);
   }
   fclose(foffsetAgain);  
 
 
 
   return chamberObj;
 }

CSCChamberTimeCorrectionsValues::ReturnType CSCChamberTimeCorrectionsValues::produceChamberTimeCorrections ( const CSCChamberTimeCorrectionsRcd & iRecord )

Definition at line 39 of file CSCChamberTimeCorrectionsValues.cc.

References isForMC, ME11offsetData, ME11offsetMC, nonME11offsetData, nonME11offsetMC, and prefill().

Referenced by CSCChamberTimeCorrectionsValues().

 {
   //need a new object so to not be deleted at exit
   return ReturnType(prefill(isForMC, isForMC ? ME11offsetMC : ME11offsetData, isForMC ? nonME11offsetMC : nonME11offsetData));
 
   
 }

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

overrideprivatevirtual

Implements edm::EventSetupRecordIntervalFinder.

Definition at line 47 of file CSCChamberTimeCorrectionsValues.cc.

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

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

Member Data Documentation

bool CSCChamberTimeCorrectionsValues::isForMC

private

Definition at line 39 of file CSCChamberTimeCorrectionsValues.h.

Referenced by CSCChamberTimeCorrectionsValues(), and produceChamberTimeCorrections().

float CSCChamberTimeCorrectionsValues::ME11offsetData

private

Definition at line 41 of file CSCChamberTimeCorrectionsValues.h.

Referenced by CSCChamberTimeCorrectionsValues(), and produceChamberTimeCorrections().

float CSCChamberTimeCorrectionsValues::ME11offsetMC

private

Definition at line 40 of file CSCChamberTimeCorrectionsValues.h.

Referenced by CSCChamberTimeCorrectionsValues(), and produceChamberTimeCorrections().

float CSCChamberTimeCorrectionsValues::nonME11offsetData

private

Definition at line 43 of file CSCChamberTimeCorrectionsValues.h.

Referenced by CSCChamberTimeCorrectionsValues(), and produceChamberTimeCorrections().

float CSCChamberTimeCorrectionsValues::nonME11offsetMC

private

Definition at line 42 of file CSCChamberTimeCorrectionsValues.h.

Referenced by CSCChamberTimeCorrectionsValues(), and produceChamberTimeCorrections().

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation