|
|
|
#include <DTT0Calibration.h>
Public Member Functions | |
| void | analyze (const edm::Event &event, const edm::EventSetup &eventSetup) |
| Fill the maps with t0 (by channel) | |
| DTT0Calibration (const edm::ParameterSet &pset) | |
| Constructor. | |
| void | endJob () |
| Compute the mean and the RMS of the t0 from the maps and write them to the DB with channel granularity. | |
| virtual | ~DTT0Calibration () |
| Destructor. | |
Private Member Functions | |
| std::string | getHistoName (const DTWireId &wId) const |
| std::string | getHistoName (const DTLayerId &lId) const |
Private Attributes | |
| std::vector< std::string > | cellsWithHistos |
| bool | correctByChamberMean_ |
| bool | debug |
| std::string | digiLabel |
| edm::ESHandle< DTGeometry > | dtGeom |
| unsigned int | eventsForLayerT0 |
| unsigned int | eventsForWireT0 |
| TH1D * | hT0SectorHisto |
| std::map< DTWireId, double > | mK |
| std::map< DTWireId, double > | mK_ref |
| std::map< DTWireId, int > | nDigiPerWire |
| std::map< DTWireId, int > | nDigiPerWire_ref |
| unsigned int | nevents |
| std::map< DTWireId, double > | qK |
| unsigned int | rejectDigiFromPeak |
| int | selSector |
| int | selWheel |
| std::map< DTWireId, double > | theAbsoluteT0PerWire |
| std::string | theCalibSector |
| std::string | theCalibWheel |
| TFile * | theFile |
| std::map< DTLayerId, TH1I * > | theHistoLayerMap |
| std::map< DTWireId, TH1I * > | theHistoWireMap |
| std::map< DTWireId, TH1I * > | theHistoWireMap_ref |
| TFile * | theOutputFile |
| std::map< DTWireId, double > | theRelativeT0PerWire |
| std::map< std::string, double > | theSigmaT0LayerMap |
| std::map< DTWireId, double > | theSigmaT0PerWire |
| std::map< std::string, double > | theT0LayerMap |
| double | tpPeakWidth |
| std::vector< DTWireId > | wireIdWithHistos |
Analyzer class computes the mean and RMS of t0 from pulses. Those values are written in the DB with cell granularity. The mean value for each channel is normalized to a reference time common to all the sector. The t0 of wires in odd layers are corrected for the relative difference between odd and even layers
Definition at line 30 of file DTT0Calibration.h.
| DTT0Calibration::DTT0Calibration | ( | const edm::ParameterSet & | pset | ) |
Constructor.
Definition at line 29 of file DTT0Calibration.cc.
References gather_cfg::cout, debug, dtTPAnalyzer_cfg::digiLabel, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), nevents, dtT0Analyzer_cfg::rootFileName, relativeConstraints::station, interactiveExample::theFile, and createXMLFile::wire.
{
// Get the debug parameter for verbose output
debug = pset.getUntrackedParameter<bool>("debug");
if(debug)
cout << "[DTT0Calibration]Constructor called!" << endl;
// Get the label to retrieve digis from the event
digiLabel = pset.getUntrackedParameter<string>("digiLabel");
// The root file which contain the histos per layer
string rootFileName = pset.getUntrackedParameter<string>("rootFileName","DTT0PerLayer.root");
theFile = new TFile(rootFileName.c_str(), "RECREATE");
theCalibWheel = pset.getUntrackedParameter<string>("calibWheel", "All"); //FIXME amke a vector of integer instead of a string
if(theCalibWheel != "All") {
stringstream linestr;
int selWheel;
linestr << theCalibWheel;
linestr >> selWheel;
cout << "[DTT0CalibrationPerLayer] chosen wheel " << selWheel << endl;
}
// Sector/s to calibrate
theCalibSector = pset.getUntrackedParameter<string>("calibSector", "All"); //FIXME amke a vector of integer instead of a string
if(theCalibSector != "All") {
stringstream linestr;
int selSector;
linestr << theCalibSector;
linestr >> selSector;
cout << "[DTT0CalibrationPerLayer] chosen sector " << selSector << endl;
}
vector<string> defaultCell;
defaultCell.push_back("None");
cellsWithHistos = pset.getUntrackedParameter<vector<string> >("cellsWithHisto", defaultCell);
for(vector<string>::const_iterator cell = cellsWithHistos.begin(); cell != cellsWithHistos.end(); cell++){
if((*cell) != "None"){
stringstream linestr;
int wheel,sector,station,sl,layer,wire;
linestr << (*cell);
linestr >> wheel >> sector >> station >> sl >> layer >> wire;
wireIdWithHistos.push_back(DTWireId(wheel,station,sector,sl,layer,wire));
}
}
hT0SectorHisto=0;
nevents=0;
eventsForLayerT0 = pset.getParameter<unsigned int>("eventsForLayerT0");
eventsForWireT0 = pset.getParameter<unsigned int>("eventsForWireT0");
rejectDigiFromPeak = pset.getParameter<unsigned int>("rejectDigiFromPeak");
tpPeakWidth = pset.getParameter<double>("tpPeakWidth");
//useReferenceWireInLayer_ = true;
correctByChamberMean_ = pset.getParameter<bool>("correctByChamberMean");
}
| DTT0Calibration::~DTT0Calibration | ( | ) | [virtual] |
Destructor.
Definition at line 86 of file DTT0Calibration.cc.
References gather_cfg::cout, debug, and interactiveExample::theFile.
| void DTT0Calibration::analyze | ( | const edm::Event & | event, |
| const edm::EventSetup & | eventSetup | ||
| ) | [virtual] |
Fill the maps with t0 (by channel)
Perform the real analysis.
Implements edm::EDAnalyzer.
Definition at line 94 of file DTT0Calibration.cc.
References abs, DTSuperLayerId::chamberId(), gather_cfg::cout, debug, dtTPAnalyzer_cfg::digiLabel, edm::EventID::event(), spr::find(), edm::EventSetup::get(), mergeVDriftHistosByStation::getHistoName(), edm::EventBase::id(), nevents, edm::EventID::run(), DTChamberId::sector(), DTLayerId::superlayerId(), interactiveExample::theFile, and DTChamberId::wheel().
{
if(debug || event.id().event() % 500==0)
cout << "--- [DTT0Calibration] Analysing Event: #Run: " << event.id().run()
<< " #Event: " << event.id().event() << endl;
nevents++;
// Get the digis from the event
Handle<DTDigiCollection> digis;
event.getByLabel(digiLabel, digis);
// Get the DT Geometry
eventSetup.get<MuonGeometryRecord>().get(dtGeom);
// Iterate through all digi collections ordered by LayerId
DTDigiCollection::DigiRangeIterator dtLayerIt;
for (dtLayerIt = digis->begin();
dtLayerIt != digis->end();
++dtLayerIt){
// Get the iterators over the digis associated with this LayerId
const DTDigiCollection::Range& digiRange = (*dtLayerIt).second;
// Get the layerId
const DTLayerId layerId = (*dtLayerIt).first; //FIXME: check to be in the right sector
if((theCalibWheel != "All") && (layerId.superlayerId().chamberId().wheel() != selWheel))
continue;
if((theCalibSector != "All") && (layerId.superlayerId().chamberId().sector() != selSector))
continue;
//if(debug) {
// cout << "Layer " << layerId<<" with "<<distance(digiRange.first, digiRange.second)<<" digi"<<endl;
//}
// Loop over all digis in the given layer
for (DTDigiCollection::const_iterator digi = digiRange.first;
digi != digiRange.second;
digi++) {
double t0 = (*digi).countsTDC();
//Use first bunch of events to fill t0 per layer
if(nevents < eventsForLayerT0){
//Get the per-layer histo from the map
TH1I *hT0LayerHisto = theHistoLayerMap[layerId];
//If it doesn't exist, book it
if(hT0LayerHisto == 0){
theFile->cd();
hT0LayerHisto = new TH1I(getHistoName(layerId).c_str(),
"T0 from pulses by layer (TDC counts, 1 TDC count = 0.781 ns)",
200, t0-100, t0+100);
if(debug)
cout << " New T0 per Layer Histo: " << hT0LayerHisto->GetName() << endl;
theHistoLayerMap[layerId] = hT0LayerHisto;
}
//Fill the histos
theFile->cd();
if(hT0LayerHisto != 0) {
// if(debug)
// cout<<"Filling histo "<<hT0LayerHisto->GetName()<<" with digi "<<t0<<" TDC counts"<<endl;
hT0LayerHisto->Fill(t0);
}
}
//Use all the remaining events to compute t0 per wire
if(nevents>eventsForLayerT0){
// Get the wireId
const DTWireId wireId(layerId, (*digi).wire());
if(debug) {
cout << " Wire: " << wireId << endl
<< " time (TDC counts): " << (*digi).countsTDC()<< endl;
}
//Fill the histos per wire for the chosen cells
vector<DTWireId>::iterator it_wire = find(wireIdWithHistos.begin(),wireIdWithHistos.end(),wireId);
if(it_wire != wireIdWithHistos.end()){
if(theHistoWireMap.find(wireId) == theHistoWireMap.end()){
theHistoWireMap[wireId] = new TH1I(getHistoName(wireId).c_str(),"T0 from pulses by wire (TDC counts, 1 TDC count = 0.781 ns)",7000,0,7000);
if(debug) cout << " New T0 per wire Histo: " << (theHistoWireMap[wireId])->GetName() << endl;
}
if(theHistoWireMap_ref.find(wireId) == theHistoWireMap_ref.end()){
theHistoWireMap_ref[wireId] = new TH1I((getHistoName(wireId) + "_ref").c_str(),"T0 from pulses by wire (TDC counts, 1 TDC count = 0.781 ns)",7000,0,7000);
if(debug) cout << " New T0 per wire Histo: " << (theHistoWireMap_ref[wireId])->GetName() << endl;
}
TH1I* hT0WireHisto = theHistoWireMap[wireId];
//Fill the histos
theFile->cd();
if(hT0WireHisto) hT0WireHisto->Fill(t0);
}
//Check the tzero has reasonable value
if(abs(hT0SectorHisto->GetBinCenter(hT0SectorHisto->GetMaximumBin()) - t0) > rejectDigiFromPeak){
if(debug)
cout<<"digi skipped because t0 per sector "<<hT0SectorHisto->GetBinCenter(hT0SectorHisto->GetMaximumBin())<<endl;
continue;
}
//Use second bunch of events to compute a t0 reference per wire
if(nevents< (eventsForLayerT0 + eventsForWireT0)){
//Fill reference wire histos
if(it_wire != wireIdWithHistos.end()){
TH1I* hT0WireHisto_ref = theHistoWireMap_ref[wireId];
theFile->cd();
if(hT0WireHisto_ref) hT0WireHisto_ref->Fill(t0);
}
if(!nDigiPerWire_ref[wireId]){
mK_ref[wireId] = 0;
}
nDigiPerWire_ref[wireId] = nDigiPerWire_ref[wireId] + 1;
mK_ref[wireId] = mK_ref[wireId] + (t0-mK_ref[wireId])/nDigiPerWire_ref[wireId];
}
//Use last all the remaining events to compute the mean and sigma t0 per wire
else if(nevents>(eventsForLayerT0 + eventsForWireT0)){
if(abs(t0-mK_ref[wireId]) > tpPeakWidth)
continue;
if(!nDigiPerWire[wireId]){
theAbsoluteT0PerWire[wireId] = 0;
qK[wireId] = 0;
mK[wireId] = 0;
}
nDigiPerWire[wireId] = nDigiPerWire[wireId] + 1;
theAbsoluteT0PerWire[wireId] = theAbsoluteT0PerWire[wireId] + t0;
//theSigmaT0PerWire[wireId] = theSigmaT0PerWire[wireId] + (t0*t0);
qK[wireId] = qK[wireId] + ((nDigiPerWire[wireId]-1)*(t0-mK[wireId])*(t0-mK[wireId])/nDigiPerWire[wireId]);
mK[wireId] = mK[wireId] + (t0-mK[wireId])/nDigiPerWire[wireId];
}
}//end if(nevents>1000)
}//end loop on digi
}//end loop on layer
//Use the t0 per layer histos to have an indication about the t0 position
if(nevents == eventsForLayerT0){
for(map<DTLayerId, TH1I*>::const_iterator lHisto = theHistoLayerMap.begin();
lHisto != theHistoLayerMap.end();
lHisto++){
if(debug)
cout<<"Reading histogram "<<(*lHisto).second->GetName()<<" with mean "<<(*lHisto).second->GetMean()<<" and RMS "<<(*lHisto).second->GetRMS();
//Take the mean as a first t0 estimation
if((*lHisto).second->GetRMS()<5.0){
if(hT0SectorHisto == 0){
hT0SectorHisto = new TH1D("hT0AllLayerOfSector","T0 from pulses per layer in sector",
//20, (*lHisto).second->GetMean()-100, (*lHisto).second->GetMean()+100);
700, 0, 7000);
}
if(debug)
cout<<" accepted"<<endl;
hT0SectorHisto->Fill((*lHisto).second->GetMean());
}
//Take the mean of noise + 400ns as a first t0 estimation
// if((*lHisto).second->GetRMS()>10.0 && ((*lHisto).second->GetRMS()<15.0)){
// double t0_estim = (*lHisto).second->GetMean() + 400;
// if(hT0SectorHisto == 0){
// hT0SectorHisto = new TH1D("hT0AllLayerOfSector","T0 from pulses per layer in sector",
// //20, t0_estim-100, t0_estim+100);
// 700, 0, 7000);
// }
// if(debug)
// cout<<" accepted + 400ns"<<endl;
// hT0SectorHisto->Fill((*lHisto).second->GetMean() + 400);
// }
if(debug)
cout<<endl;
theT0LayerMap[(*lHisto).second->GetName()] = (*lHisto).second->GetMean();
theSigmaT0LayerMap[(*lHisto).second->GetName()] = (*lHisto).second->GetRMS();
}
if(!hT0SectorHisto){
cout<<"[DTT0Calibration]: All the t0 per layer are still uncorrect: trying with greater number of events"<<endl;
eventsForLayerT0 = eventsForLayerT0*2;
return;
}
if(debug)
cout<<"[DTT0Calibration] t0 reference for this sector "<<
hT0SectorHisto->GetBinCenter(hT0SectorHisto->GetMaximumBin())<<endl;
}
}
| void DTT0Calibration::endJob | ( | void | ) | [virtual] |
Compute the mean and the RMS of the t0 from the maps and write them to the DB with channel granularity.
Loop on superlayer to correct between even-odd layers (2 different test pulse lines!)
Change t0 absolute reference -> from sector peak to chamber average
Write the t0 map into DB
Reimplemented from edm::EDAnalyzer.
Definition at line 273 of file DTT0Calibration.cc.
References abs, DTT0::begin(), DTSuperLayerId::chamberId(), DTTimeUnits::counts, gather_cfg::cout, debug, DTT0::end(), DTT0::get(), DTT0::set(), mathSSE::sqrt(), interactiveExample::theFile, tzero, and DTCalibDBUtils::writeToDB().
{
DTT0* t0sAbsolute = new DTT0();
DTT0* t0sRelative = new DTT0();
DTT0* t0sWRTChamber = new DTT0();
//if(debug)
cout << "[DTT0CalibrationPerLayer]Writing histos to file!" << endl;
theFile->cd();
hT0SectorHisto->Write();
for(map<DTWireId, TH1I*>::const_iterator wHisto = theHistoWireMap.begin();
wHisto != theHistoWireMap.end();
wHisto++) {
(*wHisto).second->Write();
}
for(map<DTWireId, TH1I*>::const_iterator wHisto = theHistoWireMap_ref.begin();
wHisto != theHistoWireMap_ref.end();
wHisto++) {
(*wHisto).second->Write();
}
for(map<DTLayerId, TH1I*>::const_iterator lHisto = theHistoLayerMap.begin();
lHisto != theHistoLayerMap.end();
lHisto++) {
(*lHisto).second->Write();
}
//if(debug)
cout << "[DTT0Calibration] Compute and store t0 and sigma per wire" << endl;
for(map<DTWireId, double>::const_iterator wiret0 = theAbsoluteT0PerWire.begin();
wiret0 != theAbsoluteT0PerWire.end();
wiret0++){
if(nDigiPerWire[(*wiret0).first]){
double t0 = (*wiret0).second/nDigiPerWire[(*wiret0).first];
theRelativeT0PerWire[(*wiret0).first] = t0 - hT0SectorHisto->GetBinCenter(hT0SectorHisto->GetMaximumBin());
//theSigmaT0PerWire[(*wiret0).first] = sqrt((theSigmaT0PerWire[(*wiret0).first] / nDigiPerWire[(*wiret0).first]) - t0*t0);
theSigmaT0PerWire[(*wiret0).first] = sqrt(qK[(*wiret0).first]/nDigiPerWire[(*wiret0).first]);
cout << "Wire " << (*wiret0).first << " has t0 " << t0 << "(absolute) "
<< theRelativeT0PerWire[(*wiret0).first] << "(relative)"
<< " sigma " << theSigmaT0PerWire[(*wiret0).first] << endl;
t0sAbsolute->set((*wiret0).first, t0, theSigmaT0PerWire[(*wiret0).first],DTTimeUnits::counts);
}
else{
cout<<"[DTT0Calibration] ERROR: no digis in wire "<<(*wiret0).first<<endl;
abort();
}
}
if(correctByChamberMean_){
// Get all the sls from the setup
const vector<DTSuperLayer*> superLayers = dtGeom->superLayers();
// Loop over all SLs
for(vector<DTSuperLayer*>::const_iterator sl = superLayers.begin();
sl != superLayers.end(); sl++) {
//Compute mean for odd and even superlayers
double oddLayersMean=0;
double evenLayersMean=0;
double oddLayersDen=0;
double evenLayersDen=0;
for(map<DTWireId, double>::const_iterator wiret0 = theRelativeT0PerWire.begin();
wiret0 != theRelativeT0PerWire.end();
wiret0++){
if((*wiret0).first.layerId().superlayerId() == (*sl)->id()){
if(debug)
cout<<"[DTT0Calibration] Superlayer "<<(*sl)->id()
<<"layer " <<(*wiret0).first.layerId().layer()<<" with "<<(*wiret0).second<<endl;
if(((*wiret0).first.layerId().layer()) % 2){
oddLayersMean = oddLayersMean + (*wiret0).second;
oddLayersDen++;
}
else{
evenLayersMean = evenLayersMean + (*wiret0).second;
evenLayersDen++;
}
}
}
oddLayersMean = oddLayersMean/oddLayersDen;
evenLayersMean = evenLayersMean/evenLayersDen;
//if(debug && oddLayersMean)
cout<<"[DTT0Calibration] Relative T0 mean for odd layers "<<oddLayersMean<<" even layers"<<evenLayersMean<<endl;
//Compute sigma for odd and even superlayers
double oddLayersSigma=0;
double evenLayersSigma=0;
for(map<DTWireId, double>::const_iterator wiret0 = theRelativeT0PerWire.begin();
wiret0 != theRelativeT0PerWire.end();
wiret0++){
if((*wiret0).first.layerId().superlayerId() == (*sl)->id()){
if(((*wiret0).first.layerId().layer()) % 2){
oddLayersSigma = oddLayersSigma + ((*wiret0).second - oddLayersMean) * ((*wiret0).second - oddLayersMean);
}
else{
evenLayersSigma = evenLayersSigma + ((*wiret0).second - evenLayersMean) * ((*wiret0).second - evenLayersMean);
}
}
}
oddLayersSigma = oddLayersSigma/oddLayersDen;
evenLayersSigma = evenLayersSigma/evenLayersDen;
oddLayersSigma = sqrt(oddLayersSigma);
evenLayersSigma = sqrt(evenLayersSigma);
//if(debug && oddLayersMean)
cout<<"[DTT0Calibration] Relative T0 sigma for odd layers "<<oddLayersSigma<<" even layers"<<evenLayersSigma<<endl;
//Recompute the mean for odd and even superlayers discarding fluctations
double oddLayersFinalMean=0;
double evenLayersFinalMean=0;
for(map<DTWireId, double>::const_iterator wiret0 = theRelativeT0PerWire.begin();
wiret0 != theRelativeT0PerWire.end();
wiret0++){
if((*wiret0).first.layerId().superlayerId() == (*sl)->id()){
if(((*wiret0).first.layerId().layer()) % 2){
if(abs((*wiret0).second - oddLayersMean) < (2*oddLayersSigma))
oddLayersFinalMean = oddLayersFinalMean + (*wiret0).second;
}
else{
if(abs((*wiret0).second - evenLayersMean) < (2*evenLayersSigma))
evenLayersFinalMean = evenLayersFinalMean + (*wiret0).second;
}
}
}
oddLayersFinalMean = oddLayersFinalMean/oddLayersDen;
evenLayersFinalMean = evenLayersFinalMean/evenLayersDen;
//if(debug && oddLayersMean)
cout<<"[DTT0Calibration] Final relative T0 mean for odd layers "<<oddLayersFinalMean<<" even layers"<<evenLayersFinalMean<<endl;
for(map<DTWireId, double>::const_iterator wiret0 = theRelativeT0PerWire.begin();
wiret0 != theRelativeT0PerWire.end();
wiret0++){
if((*wiret0).first.layerId().superlayerId() == (*sl)->id()){
double t0=-999;
if(((*wiret0).first.layerId().layer()) % 2)
t0 = (*wiret0).second + (evenLayersFinalMean - oddLayersFinalMean);
else
t0 = (*wiret0).second;
cout << "[DTT0Calibration] Wire " << (*wiret0).first << " has t0 " << (*wiret0).second
<< " (relative, after even-odd layer corrections) "
<< " sigma " << theSigmaT0PerWire[(*wiret0).first] << endl;
//Store the results into DB
t0sRelative->set((*wiret0).first, t0, theSigmaT0PerWire[(*wiret0).first],DTTimeUnits::counts);
}
}
}
//if(debug)
cout << "[DTT0Calibration]Computing relative t0 wrt to chamber average" << endl;
//Compute the reference for each chamber
map<DTChamberId,double> sumT0ByChamber;
map<DTChamberId,int> countT0ByChamber;
for(DTT0::const_iterator tzero = t0sRelative->begin();
tzero != t0sRelative->end(); tzero++) {
int channelId = tzero->channelId;
if ( channelId == 0 ) continue;
DTWireId wireId(channelId);
DTChamberId chamberId(wireId.chamberId());
//sumT0ByChamber[chamberId] = sumT0ByChamber[chamberId] + tzero->t0mean;
// @@@ better DTT0 usage
float t0mean_f;
float t0rms_f;
t0sRelative->get(wireId,t0mean_f,t0rms_f,DTTimeUnits::counts);
sumT0ByChamber[chamberId] = sumT0ByChamber[chamberId] + t0mean_f;
// @@@ NEW DTT0 END
countT0ByChamber[chamberId]++;
}
//Change reference for each wire and store the new t0s in the new map
for(DTT0::const_iterator tzero = t0sRelative->begin();
tzero != t0sRelative->end(); tzero++) {
int channelId = tzero->channelId;
if ( channelId == 0 ) continue;
DTWireId wireId(channelId);
DTChamberId chamberId(wireId.chamberId());
//double t0mean = (tzero->t0mean) - (sumT0ByChamber[chamberId]/countT0ByChamber[chamberId]);
//double t0rms = tzero->t0rms;
// @@@ better DTT0 usage
float t0mean_f;
float t0rms_f;
t0sRelative->get(wireId,t0mean_f,t0rms_f,DTTimeUnits::counts);
double t0mean = t0mean_f - (sumT0ByChamber[chamberId]/countT0ByChamber[chamberId]);
double t0rms = t0rms_f;
// @@@ NEW DTT0 END
t0sWRTChamber->set(wireId,
t0mean,
t0rms,
DTTimeUnits::counts);
//if(debug)
//cout<<"Chamber "<<chamberId<<" has reference "<<(sumT0ByChamber[chamberId]/countT0ByChamber[chamberId]);
cout << "Changing t0 of wire " << wireId << " from " << t0mean_f
<< " to " << t0mean << endl;
}
}
if(debug)
cout << "[DTT0Calibration]Writing values in DB!" << endl;
// FIXME: to be read from cfg?
string t0Record = "DTT0Rcd";
// Write the t0 map to DB
if( correctByChamberMean_ ) DTCalibDBUtils::writeToDB(t0Record, t0sWRTChamber);
else DTCalibDBUtils::writeToDB(t0Record, t0sAbsolute);
}
| string DTT0Calibration::getHistoName | ( | const DTWireId & | wId | ) | const [private] |
Definition at line 486 of file DTT0Calibration.cc.
References DTLayerId::layer(), DTChamberId::sector(), DTChamberId::station(), DTSuperLayerId::superlayer(), DTChamberId::wheel(), and DTWireId::wire().
| string DTT0Calibration::getHistoName | ( | const DTLayerId & | lId | ) | const [private] |
Definition at line 495 of file DTT0Calibration.cc.
References DTLayerId::layer(), DTChamberId::sector(), DTChamberId::station(), DTSuperLayerId::superlayer(), and DTChamberId::wheel().
{
string histoName;
stringstream theStream;
theStream << "Ch_" << lId.wheel() << "_" << lId.station() << "_" << lId.sector()
<< "_SL" << lId.superlayer() << "_L" << lId.layer() <<"_hT0Histo";
theStream >> histoName;
return histoName;
}
std::vector<std::string> DTT0Calibration::cellsWithHistos [private] |
Definition at line 94 of file DTT0Calibration.h.
bool DTT0Calibration::correctByChamberMean_ [private] |
Definition at line 79 of file DTT0Calibration.h.
bool DTT0Calibration::debug [private] |
Definition at line 55 of file DTT0Calibration.h.
std::string DTT0Calibration::digiLabel [private] |
Definition at line 58 of file DTT0Calibration.h.
edm::ESHandle<DTGeometry> DTT0Calibration::dtGeom [private] |
Definition at line 113 of file DTT0Calibration.h.
unsigned int DTT0Calibration::eventsForLayerT0 [private] |
Definition at line 68 of file DTT0Calibration.h.
unsigned int DTT0Calibration::eventsForWireT0 [private] |
Definition at line 70 of file DTT0Calibration.h.
TH1D* DTT0Calibration::hT0SectorHisto [private] |
Definition at line 90 of file DTT0Calibration.h.
std::map<DTWireId,double> DTT0Calibration::mK [private] |
Definition at line 102 of file DTT0Calibration.h.
std::map<DTWireId,double> DTT0Calibration::mK_ref [private] |
Definition at line 103 of file DTT0Calibration.h.
std::map<DTWireId,int> DTT0Calibration::nDigiPerWire [private] |
Definition at line 100 of file DTT0Calibration.h.
std::map<DTWireId,int> DTT0Calibration::nDigiPerWire_ref [private] |
Definition at line 101 of file DTT0Calibration.h.
unsigned int DTT0Calibration::nevents [private] |
Definition at line 66 of file DTT0Calibration.h.
std::map<DTWireId,double> DTT0Calibration::qK [private] |
Definition at line 104 of file DTT0Calibration.h.
unsigned int DTT0Calibration::rejectDigiFromPeak [private] |
Definition at line 73 of file DTT0Calibration.h.
int DTT0Calibration::selSector [private] |
Definition at line 85 of file DTT0Calibration.h.
int DTT0Calibration::selWheel [private] |
Definition at line 83 of file DTT0Calibration.h.
std::map<DTWireId,double> DTT0Calibration::theAbsoluteT0PerWire [private] |
Definition at line 97 of file DTT0Calibration.h.
std::string DTT0Calibration::theCalibSector [private] |
Definition at line 84 of file DTT0Calibration.h.
std::string DTT0Calibration::theCalibWheel [private] |
Definition at line 82 of file DTT0Calibration.h.
TFile* DTT0Calibration::theFile [private] |
Definition at line 61 of file DTT0Calibration.h.
std::map<DTLayerId, TH1I*> DTT0Calibration::theHistoLayerMap [private] |
Definition at line 88 of file DTT0Calibration.h.
std::map<DTWireId,TH1I*> DTT0Calibration::theHistoWireMap [private] |
Definition at line 106 of file DTT0Calibration.h.
std::map<DTWireId,TH1I*> DTT0Calibration::theHistoWireMap_ref [private] |
Definition at line 107 of file DTT0Calibration.h.
TFile* DTT0Calibration::theOutputFile [private] |
Definition at line 63 of file DTT0Calibration.h.
std::map<DTWireId,double> DTT0Calibration::theRelativeT0PerWire [private] |
Definition at line 98 of file DTT0Calibration.h.
std::map<std::string,double> DTT0Calibration::theSigmaT0LayerMap [private] |
Definition at line 110 of file DTT0Calibration.h.
std::map<DTWireId,double> DTT0Calibration::theSigmaT0PerWire [private] |
Definition at line 99 of file DTT0Calibration.h.
std::map<std::string,double> DTT0Calibration::theT0LayerMap [private] |
Definition at line 109 of file DTT0Calibration.h.
double DTT0Calibration::tpPeakWidth [private] |
Definition at line 76 of file DTT0Calibration.h.
std::vector<DTWireId> DTT0Calibration::wireIdWithHistos [private] |
Definition at line 93 of file DTT0Calibration.h.
1.7.3