#include <CSCChipSpeedCorrectionDBConditions.h>

Inheritance diagram for CSCChipSpeedCorrectionDBConditions:

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

ReturnType	produceDBChipSpeedCorrection (const CSCDBChipSpeedCorrectionRcd &)

	~CSCChipSpeedCorrectionDBConditions () 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 CSCDBChipSpeedCorrection *	prefillDBChipSpeedCorrection (bool isForMC, std::string dataCorrFileName, float dataOffse)

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
CSCDBChipSpeedCorrection *	cndbChipCorr

std::string	dataCorrFileName

float	dataOffset

bool	isForMC

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 21 of file CSCChipSpeedCorrectionDBConditions.h.

Member Typedef Documentation

typedef std::unique_ptr<CSCDBChipSpeedCorrection> CSCChipSpeedCorrectionDBConditions::ReturnType

Definition at line 28 of file CSCChipSpeedCorrectionDBConditions.h.

Constructor & Destructor Documentation

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

Definition at line 9 of file CSCChipSpeedCorrectionDBConditions.cc.

References dataCorrFileName, dataOffset, edm::ParameterSet::getUntrackedParameter(), isForMC, produceDBChipSpeedCorrection(), edm::ESProducer::setWhatProduced(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   //the following line is needed to tell the framework what
   // data is being produced
   isForMC = iConfig.getUntrackedParameter<bool>("isForMC",true);
   dataCorrFileName= iConfig.getUntrackedParameter<std::string>("dataCorrFileName","empty.txt");
   dataOffset=170.;
   // added by Zhen (changed since 1_2_0)
   setWhatProduced(this,&CSCChipSpeedCorrectionDBConditions::produceDBChipSpeedCorrection);
   findingRecord<CSCDBChipSpeedCorrectionRcd>();
   //now do what ever other initialization is needed
 
 }

CSCChipSpeedCorrectionDBConditions::~CSCChipSpeedCorrectionDBConditions ( )

override

Definition at line 24 of file CSCChipSpeedCorrectionDBConditions.cc.

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

Member Function Documentation

CSCDBChipSpeedCorrection * CSCChipSpeedCorrectionDBConditions::prefillDBChipSpeedCorrection	(	bool	isForMC,
		std::string	dataCorrFileName,
		float	dataOffse
	)

inlinestatic

Definition at line 56 of file CSCChipSpeedCorrectionDBConditions.h.

References relativeConstraints::chamber, CSCDetId::chamber(), trackerTree::check(), CSCIndexer::chipIndex(), CSCDBChipSpeedCorrection::chipSpeedCorr, cndbChipCorr, CSCIndexer::detIdFromChamberIndex(), dt, CSCDetId::endcap(), makeMuonMisalignmentScenario::endcap, CSCDBChipSpeedCorrection::factor_speedCorr, groupFilesInBlocks::fin, mps_fire::i, createfilelist::int, relativeConstraints::ring, CSCDetId::ring(), relativeConstraints::station, CSCDetId::station(), and lumiQTWidget::t.

Referenced by produceDBChipSpeedCorrection().

 {
   if (isMC)
     printf("\n Generating fake DB constants for MC\n");
   else {
     printf("\n Reading chip corrections from file %s \n",filename.data());
     printf("my data offset value is %f \n",dataOffset);
   }
   
   CSCIndexer indexer;
 
   const int CHIP_FACTOR=100;
   const int MAX_SIZE = 15768;
   const int MAX_SHORT= 32767;
   CSCDBChipSpeedCorrection * cndbChipCorr = new CSCDBChipSpeedCorrection();
 
   CSCDBChipSpeedCorrection::ChipSpeedContainer & itemvector = cndbChipCorr->chipSpeedCorr;
   itemvector.resize(MAX_SIZE);
   cndbChipCorr->factor_speedCorr= int (CHIP_FACTOR);
 
   //Fill chip corrections for MC is very simple
   if (isMC){
     for (int i=0;i<MAX_SIZE;i++){
       itemvector[i].speedCorr = 0;
     }
     return cndbChipCorr;
   }
 
   //Filling for data takes a little more time
   FILE *fin = fopen(filename.data(),"r");
 
   int serialChamber,endcap,station,ring,chamber,chip;
   float t,dt;
   int nPulses;
 
   std::vector<int> new_index_id;
   std::vector<float> new_chipPulse;
   double runningTotal = 0;
   int numNonZero = 0;
   while (!feof(fin)){
     //note space at end of format string to convert last \n
     // int check = fscanf(fin,"%d %d %f %f \n",&serialChamber,&chip,&t,&dt);
     int check = fscanf(fin,"%d %d %d %d %d %d %f %f %d \n",&serialChamber,&endcap,&station,&ring,&chamber,&chip,&t,&dt,&nPulses);
     if (check != 9){
       printf("The input file format is not as expected\n");
       assert(0);
     }
    
     int serialChamber_safecopy = serialChamber;
     // Now to map from S. Durkin's serialChamber index convention 
     // ME+1/1-ME+41 = 1  -234
     // ME+4/2       = 235-270*
     // ME-1/1-ME-41 = 271-504*
     // ME-4/2       = 505-540
     // To the convention used in /DataFormats/MuonDetId/interface/CSCIndexer.h
     // ME+1/1-ME+41 = 1  -234
     // ME+4/2       = 469-504*
     // ME-1/1-ME-41 = 235-468*
     // ME-4/2       = 505-540
     // Only the chambers marked * need to be remapped
     
     if (serialChamber >=235 && serialChamber <=270)
       serialChamber += 234;
     else { // not in ME+4\2
       if (serialChamber >=271 && serialChamber <=504)
     serialChamber -= 36;
     }
 
     CSCDetId chamberId = indexer.detIdFromChamberIndex(serialChamber);
     // Now to map from S. Durkin's chip index convention 0-29 (with 4,9,14,19,24,29 unused in ME1/3)
     // To the convention used in /DataFormats/MuonDetId/interface/CSCIndexer.h 
     // Layer 1-6 and chip 1-5 for all chambers except ME1/3 which is chip 1-4
     int layer = (chip)/5+1;
     CSCDetId cscId(chamberId.endcap(),chamberId.station(),chamberId.ring(),chamberId.chamber(),layer);
     // This should yield the same CSCDetId as decoding the chamber serial does
     // If not, some debugging is needed
     CSCDetId cscId_doubleCheck(endcap,station,ring,chamber,layer);
     if (cscId != cscId_doubleCheck){
       printf("Why doesn't chamberSerial %d map to e: %d s: %d r: %d c: %d ? \n",serialChamber_safecopy, endcap, station, ring, chamber);
       assert(0);
     }
 
     chip =(chip%5)+1;
     // The file produced by Stan starts from the geometrical strip number stored in the CSCStripDigi 
     // When the strip digis are built, the conversion from electronics channel to geometrical strip (reversing ME+1/1a and, more importantly, ME-1/1b) 
     // is done before the digi is filled (EventFilter/CSCRawToDigi/src/CSCCFEBData.cc).
     // Since we're filling an electronics channel database, I'll flip again ME-1/1b
     if (endcap == 2 && station ==1 && ring==1 && chip <5){
       chip = 5 - chip; //change 1-4 to 4-1
     }
 
 
     new_index_id.push_back(indexer.chipIndex(cscId,chip));
     new_chipPulse.push_back(t);
     if (t!=0){
       runningTotal += t;
       numNonZero++;
     }
   }
   fclose(fin);    
 
 
   // Fill the chip corrections with zeros to start
   for (int i=0;i<MAX_SIZE;i++){
     itemvector[i].speedCorr = 0;
   }
  
   for (unsigned int i=0;i<new_index_id.size();i++){
     if((short int)  (fabs((dataOffset-new_chipPulse[i])*CHIP_FACTOR+0.5))<MAX_SHORT) 
       itemvector[new_index_id[i]-1].speedCorr = (short int) ((dataOffset-new_chipPulse[i])*CHIP_FACTOR+0.5*(dataOffset>=new_chipPulse[i])-0.5*(dataOffset<new_chipPulse[i]));
     //printf("i= %d \t new index id = %d \t corr = %f \n",i,new_index_id[i], new_chipPulse[i]);
   }
 
   //For now, calculate the mean chip correction and use it for all chambers that don't have calibration pulse data (speedCorr ==0)
   //or had values of zero (speedCorr == dataOffset)
   //This should be a temporary fix until all chips that will read out in data have calibration information
   //Since there is only a handful out of 15K chips with values more than 3 ns away from the average, this is probably very safe 
   //to first order
   float ave = runningTotal/numNonZero;
   for (int i=0;i<MAX_SIZE;i++){
     if( itemvector[i].speedCorr == 0 ||itemvector[i].speedCorr == (short int)(dataOffset*CHIP_FACTOR+0.5) )
       itemvector[i].speedCorr = (short int) ((dataOffset-ave)*CHIP_FACTOR+0.5*(dataOffset>=ave)-0.5*(dataOffset<ave));
   }
 
  
   return cndbChipCorr;
 }

CSCChipSpeedCorrectionDBConditions::ReturnType CSCChipSpeedCorrectionDBConditions::produceDBChipSpeedCorrection ( const CSCDBChipSpeedCorrectionRcd & iRecord )

Definition at line 38 of file CSCChipSpeedCorrectionDBConditions.cc.

References dataCorrFileName, dataOffset, isForMC, and prefillDBChipSpeedCorrection().

Referenced by CSCChipSpeedCorrectionDBConditions().

 {
   //need a new object so to not be deleted at exit
   return CSCChipSpeedCorrectionDBConditions::ReturnType( prefillDBChipSpeedCorrection(isForMC,dataCorrFileName,dataOffset) );
   
 }

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

overrideprivatevirtual

Implements edm::EventSetupRecordIntervalFinder.

Definition at line 45 of file CSCChipSpeedCorrectionDBConditions.cc.

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

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

Member Data Documentation

CSCDBChipSpeedCorrection* CSCChipSpeedCorrectionDBConditions::cndbChipCorr

private

Definition at line 35 of file CSCChipSpeedCorrectionDBConditions.h.

Referenced by prefillDBChipSpeedCorrection().

std::string CSCChipSpeedCorrectionDBConditions::dataCorrFileName

private

Definition at line 40 of file CSCChipSpeedCorrectionDBConditions.h.

Referenced by CSCChipSpeedCorrectionDBConditions(), and produceDBChipSpeedCorrection().

float CSCChipSpeedCorrectionDBConditions::dataOffset

private

Definition at line 41 of file CSCChipSpeedCorrectionDBConditions.h.

Referenced by CSCChipSpeedCorrectionDBConditions(), and produceDBChipSpeedCorrection().

bool CSCChipSpeedCorrectionDBConditions::isForMC

private

Definition at line 38 of file CSCChipSpeedCorrectionDBConditions.h.

Referenced by CSCChipSpeedCorrectionDBConditions(), and produceDBChipSpeedCorrection().

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