#include <DTVDriftSegmentCalibration.h>

Inheritance diagram for DTVDriftSegmentCalibration:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override

void	beginJob () override

void	beginRun (const edm::Run &run, const edm::EventSetup &setup) override

	DTVDriftSegmentCalibration (const edm::ParameterSet &pset)

void	endJob () override

	~DTVDriftSegmentCalibration () override

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Types
typedef std::map< DTChamberId, std::vector< TH1F * > >	ChamberHistosMapTH1F

typedef std::map< DTChamberId, std::vector< TH2F * > >	ChamberHistosMapTH2F

Private Member Functions
void	bookHistos (DTChamberId)

Private Attributes
TFile *	rootFile_

DTSegmentSelector	select_

std::string	theCalibChamber_

edm::InputTag	theRecHits4DLabel_

ChamberHistosMapTH1F	theVDriftHistoMapTH1F_

ChamberHistosMapTH2F	theVDriftHistoMapTH2F_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Produces histograms from v-drift computation in segment fit to be used for v-drift calibration

Author: A. Vilela Pereira

Definition at line 22 of file DTVDriftSegmentCalibration.h.

Member Typedef Documentation

typedef std::map<DTChamberId, std::vector<TH1F*> > DTVDriftSegmentCalibration::ChamberHistosMapTH1F

private

Definition at line 35 of file DTVDriftSegmentCalibration.h.

typedef std::map<DTChamberId, std::vector<TH2F*> > DTVDriftSegmentCalibration::ChamberHistosMapTH2F

private

Definition at line 36 of file DTVDriftSegmentCalibration.h.

Constructor & Destructor Documentation

DTVDriftSegmentCalibration::DTVDriftSegmentCalibration ( const edm::ParameterSet & pset )

Definition at line 34 of file DTVDriftSegmentCalibration.cc.

References edm::ParameterSet::getUntrackedParameter(), rootFile_, DTAnalyzerDetailed_cfi::rootFileName, and theRecHits4DLabel_.

                                                                               :
   select_(pset),
   theRecHits4DLabel_(pset.getParameter<InputTag>("recHits4DLabel")),
   //writeVDriftDB_(pset.getUntrackedParameter<bool>("writeVDriftDB", false)),
   theCalibChamber_(pset.getUntrackedParameter<string>("calibChamber", "All")) {
 
   LogVerbatim("Calibration") << "[DTVDriftSegmentCalibration] Constructor called!";
   consumes< DTRecSegment4DCollection >(edm::InputTag(theRecHits4DLabel_));
   // the root file which will contain the histos
   string rootFileName = pset.getUntrackedParameter<string>("rootFileName","DTVDriftHistos.root");
   rootFile_ = new TFile(rootFileName.c_str(), "RECREATE");
   rootFile_->cd();
 }

DTVDriftSegmentCalibration::~DTVDriftSegmentCalibration ( )

override

Definition at line 54 of file DTVDriftSegmentCalibration.cc.

References rootFile_.

                                                        {
   rootFile_->Close();
   LogVerbatim("Calibration") << "[DTVDriftSegmentCalibration] Destructor called!";
 }

Member Function Documentation

void DTVDriftSegmentCalibration::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

override

Definition at line 59 of file DTVDriftSegmentCalibration.cc.

References bookHistos(), relativeConstraints::chamber, DTGeometry::chamber(), HcalObjRepresent::Fill(), edm::EventSetup::get(), LogTrace, rootFile_, select_, theCalibChamber_, theRecHits4DLabel_, theVDriftHistoMapTH1F_, theVDriftHistoMapTH2F_, and GeomDet::toGlobal().

                                                                                           {
   rootFile_->cd();
 
   // Get the DT Geometry
   ESHandle<DTGeometry> dtGeom;
   eventSetup.get<MuonGeometryRecord>().get(dtGeom);
 
   // Get the rechit collection from the event
   Handle<DTRecSegment4DCollection> all4DSegments;
   event.getByLabel(theRecHits4DLabel_, all4DSegments); 
 
   DTChamberId chosenChamberId;
   if(theCalibChamber_ != "All") {
     stringstream linestr;
     int selWheel, selStation, selSector;
     linestr << theCalibChamber_;
     linestr >> selWheel >> selStation >> selSector;
     chosenChamberId = DTChamberId(selWheel, selStation, selSector);
     LogVerbatim("Calibration") << " Chosen chamber: " << chosenChamberId << endl;
   }
   // Loop over segments by chamber
   DTRecSegment4DCollection::id_iterator chamberIdIt;
   for(chamberIdIt = all4DSegments->id_begin(); chamberIdIt != all4DSegments->id_end(); ++chamberIdIt){
 
     // Calibrate just the chosen chamber/s    
     if((theCalibChamber_ != "All") && ((*chamberIdIt) != chosenChamberId)) continue;
 
     // Book histos
     if(theVDriftHistoMapTH1F_.find(*chamberIdIt) == theVDriftHistoMapTH1F_.end()){
       LogTrace("Calibration") << "   Booking histos for Chamber: " << *chamberIdIt;
       bookHistos(*chamberIdIt);
     }
    
     // Get the chamber from the setup
     const DTChamber* chamber = dtGeom->chamber(*chamberIdIt);
     // Get the range for the corresponding ChamberId
     DTRecSegment4DCollection::range range = all4DSegments->get((*chamberIdIt));
     // Loop over the rechits of this DetUnit
     for(DTRecSegment4DCollection::const_iterator segment  = range.first;
                                                  segment != range.second; ++segment){
 
       
       LogTrace("Calibration") << "Segment local pos (in chamber RF): " << (*segment).localPosition()
                               << "\nSegment global pos: " << chamber->toGlobal((*segment).localPosition());
       
       if( !select_(*segment, event, eventSetup) ) continue;
 
       // Fill v-drift values
       if( (*segment).hasPhi() ) {
          //if( segment->phiSegment()->ist0Valid() ){
          double segmentVDrift = segment->phiSegment()->vDrift();
          if( segmentVDrift != 0.00 ){   
         (theVDriftHistoMapTH1F_[*chamberIdIt])[0]->Fill(segmentVDrift);
             (theVDriftHistoMapTH2F_[*chamberIdIt])[0]->Fill(segment->localPosition().x(),segmentVDrift);
             (theVDriftHistoMapTH2F_[*chamberIdIt])[1]->Fill(segment->localPosition().y(),segmentVDrift);
     }
       }
       // Probably not meaningful 
       if( (*segment).hasZed() ){
          //if( segment->zSegment()->ist0Valid() ){
          double segmentVDrift = segment->zSegment()->vDrift();
          if( segmentVDrift != 0.00 ){
         (theVDriftHistoMapTH1F_[*chamberIdIt])[1]->Fill(segmentVDrift);
     }
       }
     } // DTRecSegment4DCollection::const_iterator segment
   } // DTRecSegment4DCollection::id_iterator chamberIdIt
 } // DTVDriftSegmentCalibration::analyze

void DTVDriftSegmentCalibration::beginJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 48 of file DTVDriftSegmentCalibration.cc.

                                          {
   TH1::SetDefaultSumw2(true);
 }

void DTVDriftSegmentCalibration::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

override

Definition at line 52 of file DTVDriftSegmentCalibration.cc.

52 {}

void DTVDriftSegmentCalibration::bookHistos ( DTChamberId chId )

private

Definition at line 149 of file DTVDriftSegmentCalibration.cc.

References DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, theVDriftHistoMapTH1F_, theVDriftHistoMapTH2F_, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by analyze().

                                                             {
 
   // Compose the chamber name
   std::string wheel = std::to_string(chId.wheel());
   std::string station = std::to_string(chId.station());
   std::string sector = std::to_string(chId.sector());
 
   string chHistoName =
     "_W" + wheel +
     "_St" + station +
     "_Sec" +  sector;
 
   vector<TH1F*> histosTH1F;
   histosTH1F.push_back(new TH1F(("hRPhiVDriftCorr" + chHistoName).c_str(), "v-drift corr. from Phi segments", 200, -0.4, 0.4));
   if(chId.station() != 4) histosTH1F.push_back(new TH1F(("hRZVDriftCorr" + chHistoName).c_str(), "v-drift corr. from Z segments", 200, -0.4, 0.4));
   
   vector<TH2F*> histosTH2F;
   histosTH2F.push_back(new TH2F(("hRPhiVDriftCorrVsSegmPosX" + chHistoName).c_str(), "v-drift corr. vs. segment x position", 250, -125., 125., 200, -0.4, 0.4));
   histosTH2F.push_back(new TH2F(("hRPhiVDriftCorrVsSegmPosY" + chHistoName).c_str(), "v-drift corr. vs. segment y position", 250, -125., 125., 200, -0.4, 0.4));
   //if(chId.station() != 4) ...
 
   theVDriftHistoMapTH1F_[chId] = histosTH1F;
   theVDriftHistoMapTH2F_[chId] = histosTH2F;
 }

void DTVDriftSegmentCalibration::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 128 of file DTVDriftSegmentCalibration.cc.

References rootFile_, theVDriftHistoMapTH1F_, and theVDriftHistoMapTH2F_.

Referenced by o2olib.O2ORunMgr::executeJob().

                                         {
   rootFile_->cd();
   
   LogVerbatim("Calibration") << "[DTVDriftSegmentCalibration] Writing histos to file!" << endl;
 
   for(ChamberHistosMapTH1F::const_iterator itChHistos = theVDriftHistoMapTH1F_.begin(); itChHistos != theVDriftHistoMapTH1F_.end(); ++itChHistos){
      vector<TH1F*>::const_iterator itHistTH1F = (*itChHistos).second.begin();
      vector<TH1F*>::const_iterator itHistTH1F_end = (*itChHistos).second.end();
      for(; itHistTH1F != itHistTH1F_end; ++itHistTH1F) (*itHistTH1F)->Write();
 
      vector<TH2F*>::const_iterator itHistTH2F = theVDriftHistoMapTH2F_[(*itChHistos).first].begin();
      vector<TH2F*>::const_iterator itHistTH2F_end = theVDriftHistoMapTH2F_[(*itChHistos).first].end();
      for(; itHistTH2F != itHistTH2F_end; ++itHistTH2F) (*itHistTH2F)->Write();
   }
 
   /*if(writeVDriftDB_){
      // ...
   }*/ 
 }