CSCChamberTimeCorrectionsValues Class Reference

#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 ()
ESProxyIndex const *	getTokenIndices (unsigned int iIndex) const
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &) final
	~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
	ESProducer (const ESProducer &)=delete
ESProducer const &	operator= (const ESProducer &)=delete
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 TArg >
ESConsumesCollectorT< TRecord >	setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel={})
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, 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().