|
|
|
#include <CSCChamberTimeCorrectionsValues.h>
Public Types | |
| typedef const CSCChamberTimeCorrections * | ReturnType |
Public Member Functions | |
| CSCChamberTimeCorrectionsValues (const edm::ParameterSet &) | |
| ReturnType | produceChamberTimeCorrections (const CSCChamberTimeCorrectionsRcd &) |
| ~CSCChamberTimeCorrectionsValues () | |
Static Public Member Functions | |
| static CSCChamberTimeCorrections * | prefill (bool isMC, float ME11offset, float nonME11offset) |
Private Member Functions | |
| void | setIntervalFor (const edm::eventsetup::EventSetupRecordKey &, const edm::IOVSyncValue &, edm::ValidityInterval &) |
Private Attributes | |
| CSCChamberTimeCorrections * | chamberObj |
| bool | isForMC |
| float | ME11offsetData |
| float | ME11offsetMC |
| float | nonME11offsetData |
| float | nonME11offsetMC |
Definition at line 23 of file CSCChamberTimeCorrectionsValues.h.
Definition at line 30 of file CSCChamberTimeCorrectionsValues.h.
| CSCChamberTimeCorrectionsValues::CSCChamberTimeCorrectionsValues | ( | const edm::ParameterSet & | iConfig | ) |
Definition at line 10 of file CSCChamberTimeCorrectionsValues.cc.
References chamberObj, edm::ParameterSet::getUntrackedParameter(), isForMC, ME11offsetData, ME11offsetMC, nonME11offsetData, nonME11offsetMC, prefill(), 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;
chamberObj = prefill(isForMC, isForMC ? ME11offsetMC : ME11offsetData, isForMC ? nonME11offsetMC : nonME11offsetData);
setWhatProduced(this,&CSCChamberTimeCorrectionsValues::produceChamberTimeCorrections);
findingRecord<CSCChamberTimeCorrectionsRcd>();
//now do what ever other initialization is needed
}
| CSCChamberTimeCorrectionsValues::~CSCChamberTimeCorrectionsValues | ( | ) |
Definition at line 26 of file CSCChamberTimeCorrectionsValues.cc.
References chamberObj.
{
// do anything here that needs to be done at desctruction time
// (e.g. close files, deallocate resources etc.)
delete chamberObj;
}
| CSCChamberTimeCorrections * CSCChamberTimeCorrectionsValues::prefill | ( | bool | isMC, |
| float | ME11offset, | ||
| float | nonME11offset | ||
| ) | [inline, static] |
Definition at line 53 of file CSCChamberTimeCorrectionsValues.h.
References csccableread::cable_read(), CSCChamberTimeCorrections::chamberCorrections, CSCIndexer::chamberIndex(), chamberObj, CastorDataFrameFilter_impl::check(), corr, prof2calltree::count, CSCChamberTimeCorrections::factor_precision, i, label, evf::evtn::offset(), and cond::rpcobtemp::temp.
Referenced by CSCChamberTimeCorrectionsValues().
{
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 = new csccableread ();
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 41 of file CSCChamberTimeCorrectionsValues.cc.
References chamberObj.
Referenced by CSCChamberTimeCorrectionsValues().
{
//need a new object so to not be deleted at exit
CSCChamberTimeCorrections* mydata=new CSCChamberTimeCorrections( *chamberObj );
return mydata;
}
| void CSCChamberTimeCorrectionsValues::setIntervalFor | ( | const edm::eventsetup::EventSetupRecordKey & | , |
| const edm::IOVSyncValue & | , | ||
| edm::ValidityInterval & | oValidity | ||
| ) | [private, virtual] |
Implements edm::EventSetupRecordIntervalFinder.
Definition at line 49 of file CSCChamberTimeCorrectionsValues.cc.
References edm::IOVSyncValue::beginOfTime(), and edm::IOVSyncValue::endOfTime().
{
oValidity = edm::ValidityInterval(edm::IOVSyncValue::beginOfTime(),edm::IOVSyncValue::endOfTime());
}
Definition at line 37 of file CSCChamberTimeCorrectionsValues.h.
Referenced by CSCChamberTimeCorrectionsValues(), prefill(), produceChamberTimeCorrections(), and ~CSCChamberTimeCorrectionsValues().
bool CSCChamberTimeCorrectionsValues::isForMC [private] |
Definition at line 40 of file CSCChamberTimeCorrectionsValues.h.
Referenced by CSCChamberTimeCorrectionsValues().
float CSCChamberTimeCorrectionsValues::ME11offsetData [private] |
Definition at line 42 of file CSCChamberTimeCorrectionsValues.h.
Referenced by CSCChamberTimeCorrectionsValues().
float CSCChamberTimeCorrectionsValues::ME11offsetMC [private] |
Definition at line 41 of file CSCChamberTimeCorrectionsValues.h.
Referenced by CSCChamberTimeCorrectionsValues().
float CSCChamberTimeCorrectionsValues::nonME11offsetData [private] |
Definition at line 44 of file CSCChamberTimeCorrectionsValues.h.
Referenced by CSCChamberTimeCorrectionsValues().
float CSCChamberTimeCorrectionsValues::nonME11offsetMC [private] |
Definition at line 43 of file CSCChamberTimeCorrectionsValues.h.
Referenced by CSCChamberTimeCorrectionsValues().
1.7.3