---

DTVDriftWriter Class Reference

#include <CalibMuon/DTCalibration/plugins/DTVDriftWriter.h>

Inheritance diagram for DTVDriftWriter:

List of all members.

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup)
	DTVDriftWriter (const edm::ParameterSet &pset)
	Constructor.
void	endJob ()
virtual	~DTVDriftWriter ()
	Destructor.
Private Attributes
bool	debug
edm::ParameterSet	theCalibFilePar
TFile *	theFile
DTMeanTimerFitter *	theFitter
std::string	theGranularity
DTMtime *	theMTime
std::string	theRootInputFile
std::string	theVDriftOutputFile

---

Detailed Description

Definition at line 32 of file DTVDriftWriter.h.

---

Constructor & Destructor Documentation

DTVDriftWriter::DTVDriftWriter

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 42 of file DTVDriftWriter.cc.

References GenMuonPlsPt100GeV_cfg::cout, debug, lat::endl(), edm::ParameterSet::getUntrackedParameter(), theCalibFilePar, theFile, theFitter, theGranularity, theMTime, theRootInputFile, and theVDriftOutputFile.

                                                       {
  // get selected debug option
  debug = pset.getUntrackedParameter<bool>("debug", "false");

  // Open the root file which contains the histos
  theRootInputFile = pset.getUntrackedParameter<string>("rootFileName");
  theFile = new TFile(theRootInputFile.c_str(), "READ");
  
  theFitter = new DTMeanTimerFitter(theFile);
  if(debug)
    theFitter->setVerbosity(1);

  // the text file which will contain the histos
  theVDriftOutputFile = pset.getUntrackedParameter<string>("vDriftFileName");

  // get parameter set for DTCalibrationMap constructor
  theCalibFilePar =  pset.getUntrackedParameter<ParameterSet>("calibFileConfig");

  // the granularity to be used for calib constants evaluation
  theGranularity = pset.getUntrackedParameter<string>("calibGranularity","bySL");
  
  theMTime = new DTMtime();

  if(debug)
    cout << "[DTVDriftWriter]Constructor called!" << endl;
}

DTVDriftWriter::~DTVDriftWriter

(

)

[virtual]

Destructor.

Definition at line 70 of file DTVDriftWriter.cc.

References GenMuonPlsPt100GeV_cfg::cout, debug, lat::endl(), theFile, and theFitter.

                               {
  if(debug)
    cout << "[DTVDriftWriter]Destructor called!" << endl;
  theFile->Close();
  delete theFitter;
}

---

Member Function Documentation

void DTVDriftWriter::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)			[virtual]

Implements edm::EDAnalyzer.

Definition at line 77 of file DTVDriftWriter.cc.

References DTCalibrationMap::addCell(), GenMuonPlsPt100GeV_cfg::cout, debug, lat::endl(), DTMeanTimerFitter::evaluateVDriftAndReso(), edm::EventSetup::get(), DTCalibrationMap::getConsts(), DTCalibrationMap::getKey(), N, DTTimeUnits::ns, DTChamberId::sector(), DTMtime::set(), DTChamberId::station(), DTSuperLayerId::superLayer(), theCalibFilePar, theFitter, theGranularity, theMTime, theVDriftOutputFile, DTChamberId::wheel(), and DTCalibrationMap::writeConsts().

                                                                              {
  if(debug)
    cout << "[DTVDriftWriter]Analyzer called!" << endl;

  // Instantiate a DTCalibrationMap object 
  DTCalibrationMap calibValuesFile(theCalibFilePar);  

  // Get the DT Geometry
  ESHandle<DTGeometry> dtGeom;
  eventSetup.get<MuonGeometryRecord>().get(dtGeom);

  if(theGranularity == "bySL") {    
    // Get all the sls from the setup
    const vector<DTSuperLayer*> superLayers = dtGeom->superLayers(); 
    
    // Loop over all SLs
    for(vector<DTSuperLayer*>::const_iterator  slCell = superLayers.begin();
        slCell != superLayers.end(); slCell++) {
      
      DTSuperLayerId slId = (*slCell)->id();
      // evaluate v_drift and sigma from the TMax histograms
      DTWireId wireId(slId, 0, 0);
      vector<float> newConstants;
      TString N=(((((TString) "TMax"+(long) wireId.wheel()) +(long) wireId.station())
                  +(long) wireId.sector())+(long) wireId.superLayer());
      vector<float> vDriftAndReso = theFitter->evaluateVDriftAndReso(N);

      // Don't write the constants for the SL if the vdrift was not computed
      if(vDriftAndReso.front() == -1)
        continue;

      const DTCalibrationMap::CalibConsts* oldConstants = calibValuesFile.getConsts(wireId);
      
      if(oldConstants != 0) {
        newConstants.push_back((*oldConstants)[0]);
        newConstants.push_back((*oldConstants)[1]);
      } else {
        newConstants.push_back(-1);
        newConstants.push_back(-1);
      }
      for(int ivd=0; ivd<=5;ivd++) { 
        // 0=vdrift, 1=reso, 2=(3deltat0-2deltat0), 3=(2deltat0-1deltat0),
        //  4=(1deltat0-0deltat0), 5=deltat0 from hists with max entries,
        newConstants.push_back(vDriftAndReso[ivd]); 
      }
      calibValuesFile.addCell(calibValuesFile.getKey(wireId), newConstants);

      theMTime->set(slId,
                    vDriftAndReso[0],
                    vDriftAndReso[1],
                    DTTimeUnits::ns);
      if(debug) {
        cout << " SL: " << slId
             << " vDrift = " << vDriftAndReso[0]
             << " reso = " << vDriftAndReso[1] << endl;
      }
    }
  }
  // to be implemented: granularity different from bySL

  //   if(theGranularity == "byChamber") {
  //     const vector<DTChamber*> chambers = dMap.chambers();
    
  //     // Loop over all chambers
  //     for(vector<MuBarChamber*>::const_iterator chamber = chambers.begin();
  //    chamber != chambers.end(); chamber ++) {
  //       MuBarChamberId chamber_id = (*chamber)->id();
  //       MuBarDigiParameters::Key wire_id(chamber_id, 0, 0, 0);
  //       vector<float> newConstants;
  //       vector<float> vDriftAndReso = evaluateVDriftAndReso(wire_id, f);
  //       const CalibConsts* oldConstants = digiParams.getConsts(wire_id);
  //       if(oldConstants !=0) {
  //    newConstants = *oldConstants;
  //    newConstants.push_back(vDriftAndReso[0]);
  //    newConstants.push_back(vDriftAndReso[1]);
  //    newConstants.push_back(vDriftAndReso[2]);
  //    newConstants.push_back(vDriftAndReso[3]);
  //       } else {
  //    newConstants.push_back(-1);
  //    newConstants.push_back(-1);
  //    newConstants.push_back(vDriftAndReso[0]);
  //    newConstants.push_back(vDriftAndReso[1]);
  //    newConstants.push_back(vDriftAndReso[2]);
  //    newConstants.push_back(vDriftAndReso[3]);
  //       }
  //       digiParams.addCell(wire_id, newConstants);
  //     }
  //   }
  //write values to a table  
  calibValuesFile.writeConsts(theVDriftOutputFile);
}

void DTVDriftWriter::endJob

(

void

)

[virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 171 of file DTVDriftWriter.cc.

References GenMuonPlsPt100GeV_cfg::cout, debug, lat::endl(), ecalRecalibSequence_cff::record, and theMTime.

                            {

  if(debug) 
    cout << "[DTVDriftWriter]Writing vdrift object to DB!" << endl;

  // Write the ttrig object to DB
  string record = "DTMtimeRcd";
  DTCalibDBUtils::writeToDB<DTMtime>(record, theMTime);
}

---

Member Data Documentation

bool DTVDriftWriter::debug [private]

Definition at line 49 of file DTVDriftWriter.h.

Referenced by analyze(), DTVDriftWriter(), endJob(), and ~DTVDriftWriter().

edm::ParameterSet DTVDriftWriter::theCalibFilePar [private]

Definition at line 61 of file DTVDriftWriter.h.

Referenced by analyze(), and DTVDriftWriter().

TFile* DTVDriftWriter::theFile [private]

Definition at line 52 of file DTVDriftWriter.h.

Referenced by DTVDriftWriter(), and ~DTVDriftWriter().

DTMeanTimerFitter* DTVDriftWriter::theFitter [private]

Definition at line 67 of file DTVDriftWriter.h.

Referenced by analyze(), DTVDriftWriter(), and ~DTVDriftWriter().

std::string DTVDriftWriter::theGranularity [private]

Definition at line 64 of file DTVDriftWriter.h.

Referenced by analyze(), and DTVDriftWriter().

DTMtime* DTVDriftWriter::theMTime [private]

Definition at line 70 of file DTVDriftWriter.h.

Referenced by analyze(), DTVDriftWriter(), and endJob().

std::string DTVDriftWriter::theRootInputFile [private]

Definition at line 55 of file DTVDriftWriter.h.

Referenced by DTVDriftWriter().

std::string DTVDriftWriter::theVDriftOutputFile [private]

Definition at line 58 of file DTVDriftWriter.h.

Referenced by analyze(), and DTVDriftWriter().

---

The documentation for this class was generated from the following files:

• CalibMuon/DTCalibration/plugins/DTVDriftWriter.h
• CalibMuon/DTCalibration/plugins/DTVDriftWriter.cc

---