#include <CSCChamberTimeCorrectionsValues.h>

Inheritance diagram for CSCChamberTimeCorrectionsValues:

Public Types
typedef std::unique_ptr< CSCChamberTimeCorrections >	ReturnType

Public Member Functions
	CSCChamberTimeCorrectionsValues (const edm::ParameterSet &)

ReturnType	produceChamberTimeCorrections (const CSCChamberTimeCorrectionsRcd &)

	~CSCChamberTimeCorrectionsValues () override

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

	ESProducer (const ESProducer &)=delete

	ESProducer (ESProducer &&)=delete

ESResolverIndex const *	getTokenIndices (unsigned int iIndex) const

ESRecordIndex const *	getTokenRecordIndices (unsigned int iIndex) const

bool	hasMayConsumes () const noexcept

size_t	numberOfTokenIndices (unsigned int iIndex) const

ESProducer &	operator= (const ESProducer &)=delete

ESProducer &	operator= (ESProducer &&)=delete

SerialTaskQueueChain &	queue ()

template<typename Record >
std::optional< std::vector< ESResolverIndex > >	updateFromMayConsumes (unsigned int iIndex, const Record &iRecord) const

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &) final

	~ESProducer () noexcept(false) override

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

	ESProductResolverFactoryProducer (const ESProductResolverFactoryProducer &)=delete

const ESProductResolverFactoryProducer &	operator= (const ESProductResolverFactoryProducer &)=delete

	~ESProductResolverFactoryProducer () noexcept(false) override

Public Member Functions inherited from edm::eventsetup::ESProductResolverProvider
void	createKeyedResolvers (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)

const ComponentDescription &	description () const

	ESProductResolverProvider ()

	ESProductResolverProvider (const ESProductResolverProvider &)=delete

void	fillRecordsNotAllowingConcurrentIOVs (std::set< EventSetupRecordKey > &recordsNotAllowingConcurrentIOVs) const

virtual void	initConcurrentIOVs (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)

bool	isUsingRecord (const EventSetupRecordKey &key) const

KeyedResolvers &	keyedResolvers (const EventSetupRecordKey &iRecordKey, unsigned int iovIndex=0)

const ESProductResolverProvider &	operator= (const ESProductResolverProvider &)=delete

void	setAppendToDataLabel (const edm::ParameterSet &)

void	setDescription (const ComponentDescription &iDescription)

std::set< EventSetupRecordKey >	usingRecords () const

virtual	~ESProductResolverProvider () noexcept(false)

Public Member Functions inherited from edm::EventSetupRecordIntervalFinder
bool	concurrentFinder () const

const eventsetup::ComponentDescription &	descriptionForFinder () const

	EventSetupRecordIntervalFinder ()

	EventSetupRecordIntervalFinder (const EventSetupRecordIntervalFinder &)=delete

std::set< eventsetup::EventSetupRecordKey >	findingForRecords () const

const ValidityInterval &	findIntervalFor (const eventsetup::EventSetupRecordKey &, const IOVSyncValue &)

bool	nonconcurrentAndIOVNeedsUpdate (const eventsetup::EventSetupRecordKey &key, const IOVSyncValue &syncValue) const

const EventSetupRecordIntervalFinder &	operator= (const EventSetupRecordIntervalFinder &)=delete

void	resetInterval (const eventsetup::EventSetupRecordKey &)

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::ESProductResolverProvider
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 Types inherited from edm::ESProductResolverFactoryProducer
using	EventSetupRecordKey = eventsetup::EventSetupRecordKey

Protected Types inherited from edm::eventsetup::ESProductResolverProvider
using	KeyedResolversVector = std::vector< std::pair< DataKey, std::shared_ptr< ESProductResolver > >>

Protected Member Functions inherited from edm::ESProducer
ESConsumesInfo *	consumesInfoPushBackNew ()

unsigned int	consumesInfoSize () const

template<typename CallbackT , typename TList , typename TRecord >
void	registerProducts (std::shared_ptr< std::pair< unsigned int, std::shared_ptr< CallbackT >>> iCallback, const TList , const TRecord iRecord, const es::Label &iLabel)

template<typename CallbackT , typename TRecord >
void	registerProducts (std::shared_ptr< std::pair< unsigned int, std::shared_ptr< CallbackT >>>, const eventsetup::produce::Null , const TRecord , const es::Label &)

template<typename T >
auto	setWhatProduced (T *iThis, const es::Label &iLabel={})

template<typename T >
auto	setWhatProduced (T iThis, const char iLabel)

template<typename T >
auto	setWhatProduced (T *iThis, const std::string &iLabel)

template<typename T , typename TDecorator >
auto	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel={})

template<typename T , typename TReturn , typename TRecord >
auto	setWhatProduced (T iThis, TReturn(T::iMethod)(const TRecord &), const es::Label &iLabel={})

template<typename T , typename TReturn , typename TRecord , typename TDecorator >
auto	setWhatProduced (T iThis, TReturn(T ::iMethod)(const TRecord &), const TDecorator &iDec, const es::Label &iLabel={})

template<typename TFunc >
auto	setWhatProduced (TFunc &&func, const es::Label &iLabel={})

template<typename TReturn , typename TRecord , typename TFunc , typename TDecorator >
ESConsumesCollectorT< TRecord >	setWhatProduced (TFunc &&func, TDecorator &&iDec, const es::Label &iLabel={})

void	usesResources (std::vector< std::string > const &)

Protected Member Functions inherited from edm::ESProductResolverFactoryProducer
template<class TFactory >
void	registerFactory (std::unique_ptr< TFactory > iFactory, const std::string &iLabel=std::string())

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

KeyedResolversVector	registerResolvers (const EventSetupRecordKey &, unsigned int iovIndex) override

Protected Member Functions inherited from edm::eventsetup::ESProductResolverProvider
template<class T >
void	usingRecord ()

void	usingRecordWithKey (const EventSetupRecordKey &key)

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

◆ ReturnType

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

Definition at line 28 of file CSCChamberTimeCorrectionsValues.h.

Constructor & Destructor Documentation

◆ CSCChamberTimeCorrectionsValues()

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::~CSCChamberTimeCorrectionsValues ( )

override

Definition at line 23 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

◆ prefill()

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

inlinestatic

Definition at line 53 of file CSCChamberTimeCorrectionsValues.h.

References cms::cuda::assert(), csccableread::cable_read(), relativeConstraints::chamber, CSCChamberTimeCorrections::chamberCorrections, CSCIndexer::chamberIndex(), RPCNoise_example::check, alignCSCRings::corr, submitPVResolutionJobs::count, phase2TrackerDigitizer_cfi::delay, CSCChamberTimeCorrections::factor_precision, groupFilesInBlocks::fin, nano_mu_digi_cff::float, mps_fire::i, createfilelist::int, PV_cfg::isMC, label, pixelgpudetails::MAX_SIZE, hltrates_dqm_sourceclient-live_cfg::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");
   //must add space for null terminator
   char label[1024 + 1];
   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;
 }

◆ produceChamberTimeCorrections()

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

Definition at line 33 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));
 }

◆ setIntervalFor()

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

overrideprivatevirtual

Implements edm::EventSetupRecordIntervalFinder.

Definition at line 40 of file CSCChamberTimeCorrectionsValues.cc.

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

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