user includes More...

#include <SiStripGainsCalibTreeWorker.cc>

Inheritance diagram for SiStripGainsCalibTreeWorker:

Public Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &, APVGain::APVGainHistograms &) const override

void	dqmAnalyze (edm::Event const &, edm::EventSetup const &, APVGain::APVGainHistograms const &) const override

	SiStripGainsCalibTreeWorker (const edm::ParameterSet &)

Public Member Functions inherited from DQMGlobalEDAnalyzer< APVGain::APVGainHistograms >
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &, APVGain::APVGainHistogramsconst &) const

void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup) const final

Public Member Functions inherited from DQMGlobalEDAnalyzerBase< APVGain::APVGainHistograms, Args...>
void	accumulate (edm::StreamID id, edm::Event const &event, edm::EventSetup const &setup) const final

	DQMGlobalEDAnalyzerBase ()

std::shared_ptr < APVGain::APVGainHistograms >	globalBeginRun (edm::Run const &run, edm::EventSetup const &setup) const final

void	globalEndRun (edm::Run const &, edm::EventSetup const &) const final

Public Member Functions inherited from edm::global::EDProducer< T >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

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

std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	checkBookAPVColls (const TrackerGeometry *bareTkGeomPtr, APVGain::APVGainHistograms &histograms) const

void	dqmBeginRun (edm::Run const &, edm::EventSetup const &, APVGain::APVGainHistograms &) const override

int	statCollectionFromMode (const char *tag) const

Private Attributes
bool	AllowSaturation

std::string	CalibPrefix_

std::string	CalibSuffix_

std::vector< std::string >	calibTreeFileNames_

std::string	calibTreeName_

bool	doChargeMonitorPerPlane

std::vector< std::string >	dqm_tag_

std::string	EventPrefix_

std::string	EventSuffix_

bool	FirstSetOfConstants

edm::ESGetToken< SiStripGain, SiStripGainRcd >	gainToken_

bool	hasProcessed_

std::string	m_calibrationMode

std::string	m_DQMdir

unsigned int	MaxNrStrips

double	MaxTrackChiOverNdf

double	MaxTrackEta

int	MaxTrackingIteration

double	MaxTrackMomentum

double	MinTrackEta

unsigned int	MinTrackHits

double	MinTrackMomentum

bool	OldGainRemoving

edm::ESGetToken < SiStripQuality, SiStripQualityRcd >	qualityToken_

std::map< unsigned int, APVloc >	theTopologyMap

edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_

edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomTokenBR_

std::string	TrackPrefix_

std::string	TrackSuffix_

edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	tTopoToken_

bool	useCalibration

bool	Validation

std::vector< std::string >	VChargeHisto

Additional Inherited Members
Public Types inherited from DQMGlobalEDAnalyzerBase< APVGain::APVGainHistograms, Args...>
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from DQMGlobalEDAnalyzerBase< APVGain::APVGainHistograms, Args...>
uint64_t	meId (edm::Run const &run) const

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Protected Attributes inherited from DQMGlobalEDAnalyzerBase< APVGain::APVGainHistograms, Args...>
DQMStore *	dqmstore_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

user includes

Description: Fill DQM histograms with SiStrip Charge normalized to path length

Definition at line 50 of file SiStripGainsCalibTreeWorker.cc.

Constructor & Destructor Documentation

SiStripGainsCalibTreeWorker::SiStripGainsCalibTreeWorker ( const edm::ParameterSet & iConfig )

explicit

Definition at line 127 of file SiStripGainsCalibTreeWorker.cc.

References AllowSaturation, CalibPrefix_, CalibSuffix_, calibTreeFileNames_, calibTreeName_, doChargeMonitorPerPlane, dqm_tag_, edm::EDConsumerBase::esConsumes(), EventPrefix_, EventSuffix_, first, FirstSetOfConstants, gainToken_, edm::ParameterSet::getUntrackedParameter(), hasProcessed_, mps_fire::i, m_calibrationMode, m_DQMdir, MaxNrStrips, MaxTrackChiOverNdf, MaxTrackEta, MaxTrackingIteration, MaxTrackMomentum, MinTrackEta, MinTrackHits, MinTrackMomentum, APVGain::monHnames(), OldGainRemoving, qualityToken_, alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, APVGain::subdetectorId(), APVGain::subdetectorPlane(), APVGain::subdetectorSide(), theTopologyMap, APVGain::thickness(), tkGeomTokenBR_, TrackPrefix_, TrackSuffix_, tTopoToken_, useCalibration, Validation, and VChargeHisto.

                                                                                        {
   MinTrackMomentum = iConfig.getUntrackedParameter<double>("minTrackMomentum", 3.0);
   MaxTrackMomentum = iConfig.getUntrackedParameter<double>("maxTrackMomentum", 99999.0);
   MinTrackEta = iConfig.getUntrackedParameter<double>("minTrackEta", -5.0);
   MaxTrackEta = iConfig.getUntrackedParameter<double>("maxTrackEta", 5.0);
   MaxNrStrips = iConfig.getUntrackedParameter<unsigned>("maxNrStrips", 2);
   MinTrackHits = iConfig.getUntrackedParameter<unsigned>("MinTrackHits", 8);
   MaxTrackChiOverNdf = iConfig.getUntrackedParameter<double>("MaxTrackChiOverNdf", 3);
   MaxTrackingIteration = iConfig.getUntrackedParameter<int>("MaxTrackingIteration", 7);
   AllowSaturation = iConfig.getUntrackedParameter<bool>("AllowSaturation", false);
   FirstSetOfConstants = iConfig.getUntrackedParameter<bool>("FirstSetOfConstants", true);
   Validation = iConfig.getUntrackedParameter<bool>("Validation", false);
   OldGainRemoving = iConfig.getUntrackedParameter<bool>("OldGainRemoving", false);
   useCalibration = iConfig.getUntrackedParameter<bool>("UseCalibration", false);
   doChargeMonitorPerPlane = iConfig.getUntrackedParameter<bool>("doChargeMonitorPerPlane", false);
   m_DQMdir = iConfig.getUntrackedParameter<std::string>("DQMdir", "AlCaReco/SiStripGains");
   m_calibrationMode = iConfig.getUntrackedParameter<std::string>("calibrationMode", "StdBunch");
   VChargeHisto = iConfig.getUntrackedParameter<std::vector<std::string>>("ChargeHisto");
 
   // fill in the mapping between the histogram indices and the (id,side,plane) tuple
   std::vector<std::pair<std::string, std::string>> hnames =
       APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "");
   for (unsigned int i = 0; i < hnames.size(); i++) {
     int id = APVGain::subdetectorId((hnames[i]).first);
     int side = APVGain::subdetectorSide((hnames[i]).first);
     int plane = APVGain::subdetectorPlane((hnames[i]).first);
     int thick = APVGain::thickness((hnames[i]).first);
     std::string s = hnames[i].first;
 
     auto loc = APVloc(thick, id, side, plane, s);
     theTopologyMap.insert(std::make_pair(i, loc));
   }
 
   //Set the monitoring element tag and store
   dqm_tag_.reserve(7);
   dqm_tag_.clear();
   dqm_tag_.push_back("StdBunch");    // statistic collection from Standard Collision Bunch @ 3.8 T
   dqm_tag_.push_back("StdBunch0T");  // statistic collection from Standard Collision Bunch @ 0 T
   dqm_tag_.push_back("AagBunch");    // statistic collection from First Collision After Abort Gap @ 3.8 T
   dqm_tag_.push_back("AagBunch0T");  // statistic collection from First Collision After Abort Gap @ 0 T
   dqm_tag_.push_back("IsoMuon");     // statistic collection from Isolated Muon @ 3.8 T
   dqm_tag_.push_back("IsoMuon0T");   // statistic collection from Isolated Muon @ 0 T
   dqm_tag_.push_back("Harvest");     // statistic collection: Harvest
 
   hasProcessed_ = false;
 
   edm::ParameterSet swhallowgain_pset = iConfig.getUntrackedParameter<edm::ParameterSet>("gain");
   CalibPrefix_ = swhallowgain_pset.getUntrackedParameter<std::string>("prefix");
   CalibSuffix_ = swhallowgain_pset.getUntrackedParameter<std::string>("suffix");
   edm::ParameterSet evtinfo_pset = iConfig.getUntrackedParameter<edm::ParameterSet>("evtinfo");
   EventPrefix_ = evtinfo_pset.getUntrackedParameter<std::string>("prefix");
   EventSuffix_ = evtinfo_pset.getUntrackedParameter<std::string>("suffix");
   edm::ParameterSet track_pset = iConfig.getUntrackedParameter<edm::ParameterSet>("tracks");
   TrackPrefix_ = track_pset.getUntrackedParameter<std::string>("prefix");
   TrackSuffix_ = track_pset.getUntrackedParameter<std::string>("suffix");
 
   calibTreeName_ = iConfig.getUntrackedParameter<std::string>("CalibTreeName");
   calibTreeFileNames_ = iConfig.getUntrackedParameter<std::vector<std::string>>("CalibTreeFiles");
 
   tTopoToken_ = esConsumes();
   tkGeomTokenBR_ = esConsumes<edm::Transition::BeginRun>();
   gainToken_ = esConsumes<edm::Transition::BeginRun>();
   qualityToken_ = esConsumes<edm::Transition::BeginRun>();
 }

Member Function Documentation

void SiStripGainsCalibTreeWorker::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	run,
		edm::EventSetup const &	setup,
		APVGain::APVGainHistograms &	histograms
	)		const

overridevirtual

Implements DQMGlobalEDAnalyzerBase< APVGain::APVGainHistograms, Args...>.

Definition at line 622 of file SiStripGainsCalibTreeWorker.cc.

                                                                                              {
   ibooker.cd();
   std::string dqm_dir = m_DQMdir;
   const char* tag = dqm_tag_[statCollectionFromMode(m_calibrationMode.c_str())].c_str();
 
   edm::LogInfo("SiStripGainsCalibTreeWorker")
       << "Setting " << dqm_dir << " in DQM and booking histograms for tag " << tag << std::endl;
 
   ibooker.setCurrentFolder(dqm_dir);
 
   std::string stag(tag);
   if (!stag.empty() && stag[0] != '_')
     stag.insert(0, 1, '_');
 
   std::string cvi = std::string("Charge_Vs_Index") + stag;
   std::string cvpTIB = std::string("Charge_Vs_PathlengthTIB") + stag;
   std::string cvpTOB = std::string("Charge_Vs_PathlengthTOB") + stag;
   std::string cvpTIDP = std::string("Charge_Vs_PathlengthTIDP") + stag;
   std::string cvpTIDM = std::string("Charge_Vs_PathlengthTIDM") + stag;
   std::string cvpTECP1 = std::string("Charge_Vs_PathlengthTECP1") + stag;
   std::string cvpTECP2 = std::string("Charge_Vs_PathlengthTECP2") + stag;
   std::string cvpTECM1 = std::string("Charge_Vs_PathlengthTECM1") + stag;
   std::string cvpTECM2 = std::string("Charge_Vs_PathlengthTECM2") + stag;
 
   int elepos = statCollectionFromMode(tag);
 
   histograms.Charge_Vs_Index.reserve(dqm_tag_.size());
   histograms.Charge_Vs_PathlengthTIB.reserve(dqm_tag_.size());
   histograms.Charge_Vs_PathlengthTOB.reserve(dqm_tag_.size());
   histograms.Charge_Vs_PathlengthTIDP.reserve(dqm_tag_.size());
   histograms.Charge_Vs_PathlengthTIDM.reserve(dqm_tag_.size());
   histograms.Charge_Vs_PathlengthTECP1.reserve(dqm_tag_.size());
   histograms.Charge_Vs_PathlengthTECP2.reserve(dqm_tag_.size());
   histograms.Charge_Vs_PathlengthTECM1.reserve(dqm_tag_.size());
   histograms.Charge_Vs_PathlengthTECM2.reserve(dqm_tag_.size());
 
   // The cluster charge is stored by exploiting a non uniform binning in order
   // reduce the histogram memory size. The bin width is relaxed with a falling
   // exponential function and the bin boundaries are stored in the binYarray.
   // The binXarray is used to provide as many bins as the APVs.
   //
   // More details about this implementations are here:
   // https://indico.cern.ch/event/649344/contributions/2672267/attachments/1498323/2332518/OptimizeChHisto.pdf
 
   std::vector<float> binXarray;
   binXarray.reserve(histograms.NStripAPVs + 1);
   for (unsigned int a = 0; a <= histograms.NStripAPVs; a++) {
     binXarray.push_back((float)a);
   }
 
   std::array<float, 688> binYarray;
   double p0 = 5.445;
   double p1 = 0.002113;
   double p2 = 69.01576;
   double y = 0.;
   for (int b = 0; b < 687; b++) {
     binYarray[b] = y;
     if (y <= 902.)
       y = y + 2.;
     else
       y = (p0 - log(exp(p0 - p1 * y) - p2 * p1)) / p1;
   }
   binYarray[687] = 4000.;
 
   histograms.Charge_1[elepos].clear();
   histograms.Charge_2[elepos].clear();
   histograms.Charge_3[elepos].clear();
   histograms.Charge_4[elepos].clear();
 
   auto it = histograms.Charge_Vs_Index.begin();
   histograms.Charge_Vs_Index.insert(
       it + elepos,
       ibooker.book2S(cvi.c_str(), cvi.c_str(), histograms.NStripAPVs, &binXarray[0], 687, binYarray.data()));
 
   it = histograms.Charge_Vs_PathlengthTIB.begin();
   histograms.Charge_Vs_PathlengthTIB.insert(it + elepos,
                                             ibooker.book2S(cvpTIB.c_str(), cvpTIB.c_str(), 20, 0.3, 1.3, 250, 0, 2000));
 
   it = histograms.Charge_Vs_PathlengthTOB.begin();
   histograms.Charge_Vs_PathlengthTOB.insert(it + elepos,
                                             ibooker.book2S(cvpTOB.c_str(), cvpTOB.c_str(), 20, 0.3, 1.3, 250, 0, 2000));
 
   it = histograms.Charge_Vs_PathlengthTIDP.begin();
   histograms.Charge_Vs_PathlengthTIDP.insert(
       it + elepos, ibooker.book2S(cvpTIDP.c_str(), cvpTIDP.c_str(), 20, 0.3, 1.3, 250, 0, 2000));
 
   it = histograms.Charge_Vs_PathlengthTIDM.begin();
   histograms.Charge_Vs_PathlengthTIDM.insert(
       it + elepos, ibooker.book2S(cvpTIDM.c_str(), cvpTIDM.c_str(), 20, 0.3, 1.3, 250, 0, 2000));
 
   it = histograms.Charge_Vs_PathlengthTECP1.begin();
   histograms.Charge_Vs_PathlengthTECP1.insert(
       it + elepos, ibooker.book2S(cvpTECP1.c_str(), cvpTECP1.c_str(), 20, 0.3, 1.3, 250, 0, 2000));
 
   it = histograms.Charge_Vs_PathlengthTECP2.begin();
   histograms.Charge_Vs_PathlengthTECP2.insert(
       it + elepos, ibooker.book2S(cvpTECP2.c_str(), cvpTECP2.c_str(), 20, 0.3, 1.3, 250, 0, 2000));
 
   it = histograms.Charge_Vs_PathlengthTECM1.begin();
   histograms.Charge_Vs_PathlengthTECM1.insert(
       it + elepos, ibooker.book2S(cvpTECM1.c_str(), cvpTECM1.c_str(), 20, 0.3, 1.3, 250, 0, 2000));
 
   it = histograms.Charge_Vs_PathlengthTECM2.begin();
   histograms.Charge_Vs_PathlengthTECM2.insert(
       it + elepos, ibooker.book2S(cvpTECM2.c_str(), cvpTECM2.c_str(), 20, 0.3, 1.3, 250, 0, 2000));
 
   std::vector<std::pair<std::string, std::string>> hnames =
       APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "");
   for (unsigned int i = 0; i < hnames.size(); i++) {
     std::string htag = (hnames[i]).first + stag;
     histograms.Charge_1[elepos].push_back(ibooker.book1DD(htag.c_str(), (hnames[i]).second.c_str(), 100, 0., 1000.));
   }
 
   hnames = APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "woG2");
   for (unsigned int i = 0; i < hnames.size(); i++) {
     std::string htag = (hnames[i]).first + stag;
     histograms.Charge_2[elepos].push_back(ibooker.book1DD(htag.c_str(), (hnames[i]).second.c_str(), 100, 0., 1000.));
   }
 
   hnames = APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "woG1");
   for (unsigned int i = 0; i < hnames.size(); i++) {
     std::string htag = (hnames[i]).first + stag;
     histograms.Charge_3[elepos].push_back(ibooker.book1DD(htag.c_str(), (hnames[i]).second.c_str(), 100, 0., 1000.));
   }
 
   hnames = APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "woG1G2");
   for (unsigned int i = 0; i < hnames.size(); i++) {
     std::string htag = (hnames[i]).first + stag;
     histograms.Charge_4[elepos].push_back(ibooker.book1DD(htag.c_str(), (hnames[i]).second.c_str(), 100, 0., 1000.));
   }
 }

void SiStripGainsCalibTreeWorker::checkBookAPVColls	(	const TrackerGeometry *	bareTkGeomPtr,
		APVGain::APVGainHistograms &	histograms
	)		const

private

Definition at line 507 of file SiStripGainsCalibTreeWorker.cc.

References sistrip::APV, APVGain::APVGainHistograms::APVsColl, APVGain::APVGainHistograms::APVsCollOrdered, TrackerGeometry::dets(), mps_fire::i, dqmiolumiharvest::j, PixelTopology::ncolumns(), APVGain::APVGainHistograms::NPixelDets, PixelTopology::nrows(), APVGain::APVGainHistograms::NStripAPVs, StripTopology::nstrips(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::rawId(), ntupleEnum::SubDet, DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, and StripSubdetector::TOB.

Referenced by dqmBeginRun().

                                                                                                 {
   if (bareTkGeomPtr) {  // pointer not yet set: called the first time => fill the APVColls
     auto const& Det = bareTkGeomPtr->dets();
 
     edm::LogInfo("SiStripGainsCalibTreeWorker") << " Resetting APV struct" << std::endl;
 
     unsigned int Index = 0;
 
     for (unsigned int i = 0; i < Det.size(); i++) {
       DetId Detid = Det[i]->geographicalId();
       int SubDet = Detid.subdetId();
 
       if (SubDet == StripSubdetector::TIB || SubDet == StripSubdetector::TID || SubDet == StripSubdetector::TOB ||
           SubDet == StripSubdetector::TEC) {
         auto DetUnit = dynamic_cast<const StripGeomDetUnit*>(Det[i]);
         if (!DetUnit)
           continue;
 
         const StripTopology& Topo = DetUnit->specificTopology();
         unsigned int NAPV = Topo.nstrips() / 128;
 
         for (unsigned int j = 0; j < NAPV; j++) {
           auto APV = std::make_shared<stAPVGain>();
           APV->Index = Index;
           APV->Bin = -1;
           APV->DetId = Detid.rawId();
           APV->APVId = j;
           APV->SubDet = SubDet;
           APV->FitMPV = -1;
           APV->FitMPVErr = -1;
           APV->FitWidth = -1;
           APV->FitWidthErr = -1;
           APV->FitChi2 = -1;
           APV->FitNorm = -1;
           APV->Gain = -1;
           APV->PreviousGain = 1;
           APV->PreviousGainTick = 1;
           APV->x = DetUnit->position().basicVector().x();
           APV->y = DetUnit->position().basicVector().y();
           APV->z = DetUnit->position().basicVector().z();
           APV->Eta = DetUnit->position().basicVector().eta();
           APV->Phi = DetUnit->position().basicVector().phi();
           APV->R = DetUnit->position().basicVector().transverse();
           APV->Thickness = DetUnit->surface().bounds().thickness();
           APV->NEntries = 0;
           APV->isMasked = false;
 
           histograms.APVsCollOrdered.push_back(APV);
           histograms.APVsColl[(APV->DetId << 4) | APV->APVId] = APV;
           Index++;
           histograms.NStripAPVs++;
         }  // loop on APVs
       }    // if is Strips
     }      // loop on dets
 
     for (unsigned int i = 0; i < Det.size();
          i++) {  //Make two loop such that the Pixel information is added at the end --> make transition simpler
       DetId Detid = Det[i]->geographicalId();
       int SubDet = Detid.subdetId();
       if (SubDet == PixelSubdetector::PixelBarrel || SubDet == PixelSubdetector::PixelEndcap) {
         auto DetUnit = dynamic_cast<const PixelGeomDetUnit*>(Det[i]);
         if (!DetUnit)
           continue;
 
         const PixelTopology& Topo = DetUnit->specificTopology();
         unsigned int NROCRow = Topo.nrows() / (80.);
         unsigned int NROCCol = Topo.ncolumns() / (52.);
 
         for (unsigned int j = 0; j < NROCRow; j++) {
           for (unsigned int i = 0; i < NROCCol; i++) {
             auto APV = std::make_shared<stAPVGain>();
             APV->Index = Index;
             APV->Bin = -1;
             APV->DetId = Detid.rawId();
             APV->APVId = (j << 3 | i);
             APV->SubDet = SubDet;
             APV->FitMPV = -1;
             APV->FitMPVErr = -1;
             APV->FitWidth = -1;
             APV->FitWidthErr = -1;
             APV->FitChi2 = -1;
             APV->Gain = -1;
             APV->PreviousGain = 1;
             APV->PreviousGainTick = 1;
             APV->x = DetUnit->position().basicVector().x();
             APV->y = DetUnit->position().basicVector().y();
             APV->z = DetUnit->position().basicVector().z();
             APV->Eta = DetUnit->position().basicVector().eta();
             APV->Phi = DetUnit->position().basicVector().phi();
             APV->R = DetUnit->position().basicVector().transverse();
             APV->Thickness = DetUnit->surface().bounds().thickness();
             APV->isMasked = false;  //SiPixelQuality_->IsModuleBad(Detid.rawId());
             APV->NEntries = 0;
 
             histograms.APVsCollOrdered.push_back(APV);
             histograms.APVsColl[(APV->DetId << 4) | APV->APVId] = APV;
             Index++;
             histograms.NPixelDets++;
 
           }  // loop on ROC cols
         }    // loop on ROC rows
       }      // if Pixel
     }        // loop on Dets
   }          //if (!bareTkGeomPtr_) ...
 }

void SiStripGainsCalibTreeWorker::dqmAnalyze	(	edm::Event const &	iEvent,
		edm::EventSetup const &	iSetup,
		APVGain::APVGainHistograms const &	histograms
	)		const

overridevirtual

Implements DQMGlobalEDAnalyzerBase< APVGain::APVGainHistograms, Args...>.

Definition at line 243 of file SiStripGainsCalibTreeWorker.cc.

                                                                                                {
   if (!hasProcessed_) {
     const TrackerTopology* topo = &iSetup.getData(tTopoToken_);
 
     for (const auto& elem : theTopologyMap) {
       LogDebug("SiStripGainsCalibTreeWorker")
           << elem.first << " - " << elem.second.m_string << " " << elem.second.m_subdetectorId << " "
           << elem.second.m_subdetectorSide << " " << elem.second.m_subdetectorPlane << std::endl;
     }
 
     LogDebug("SiStripGainsCalibTreeWorker") << "for mode" << m_calibrationMode << std::endl;
 
     // data members
     unsigned int eventnumber = 0;
     unsigned int runnumber = 0;
 #ifdef ExtendedCALIBTree
     const std::vector<bool>* TrigTech = nullptr;
     const std::vector<double>* chargeoverpath = nullptr;
 #endif
     // Track data
     const std::vector<double>* trackchi2ndof = nullptr;
     const std::vector<float>* trackp = nullptr;
     const std::vector<float>* trackpt = nullptr;
     const std::vector<double>* tracketa = nullptr;
     const std::vector<double>* trackphi = nullptr;
     const std::vector<unsigned int>* trackhitsvalid = nullptr;
     const std::vector<int>* trackalgo = nullptr;
     // CalibTree data
     const std::vector<int>* trackindex = nullptr;
     const std::vector<unsigned int>* rawid = nullptr;
     const std::vector<double>* localdirx = nullptr;
     const std::vector<double>* localdiry = nullptr;
     const std::vector<double>* localdirz = nullptr;
     const std::vector<unsigned short>* firststrip = nullptr;
     const std::vector<unsigned short>* nstrips = nullptr;
     const std::vector<bool>* saturation = nullptr;
     const std::vector<bool>* overlapping = nullptr;
     const std::vector<bool>* farfromedge = nullptr;
     const std::vector<unsigned int>* charge = nullptr;
     const std::vector<double>* path = nullptr;
     const std::vector<unsigned char>* amplitude = nullptr;
     const std::vector<double>* gainused = nullptr;
     const std::vector<double>* gainusedTick = nullptr;
 
     TChain tree{calibTreeName_.c_str()};
     for (const auto& ctFn : calibTreeFileNames_) {
       tree.Add(ctFn.c_str());
     }
 
     tree.SetBranchAddress((EventPrefix_ + "event" + EventSuffix_).c_str(), &eventnumber, nullptr);
     tree.SetBranchAddress((EventPrefix_ + "run" + EventSuffix_).c_str(), &runnumber, nullptr);
 #ifdef ExtendedCALIBTree
     tree.SetBranchAddress((EventPrefix_ + "TrigTech" + EventSuffix_).c_str(), &TrigTech, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "chargeoverpath" + CalibSuffix_).c_str(), &chargeoverpath, nullptr);
 #endif
     tree.SetBranchAddress((TrackPrefix_ + "chi2ndof" + TrackSuffix_).c_str(), &trackchi2ndof, nullptr);
     tree.SetBranchAddress((TrackPrefix_ + "momentum" + TrackSuffix_).c_str(), &trackp, nullptr);
     tree.SetBranchAddress((TrackPrefix_ + "pt" + TrackSuffix_).c_str(), &trackpt, nullptr);
     tree.SetBranchAddress((TrackPrefix_ + "eta" + TrackSuffix_).c_str(), &tracketa, nullptr);
     tree.SetBranchAddress((TrackPrefix_ + "phi" + TrackSuffix_).c_str(), &trackphi, nullptr);
     tree.SetBranchAddress((TrackPrefix_ + "hitsvalid" + TrackSuffix_).c_str(), &trackhitsvalid, nullptr);
     tree.SetBranchAddress((TrackPrefix_ + "algo" + TrackSuffix_).c_str(), &trackalgo, nullptr);
 
     tree.SetBranchAddress((CalibPrefix_ + "trackindex" + CalibSuffix_).c_str(), &trackindex, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "rawid" + CalibSuffix_).c_str(), &rawid, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "localdirx" + CalibSuffix_).c_str(), &localdirx, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "localdiry" + CalibSuffix_).c_str(), &localdiry, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "localdirz" + CalibSuffix_).c_str(), &localdirz, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "firststrip" + CalibSuffix_).c_str(), &firststrip, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "nstrips" + CalibSuffix_).c_str(), &nstrips, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "saturation" + CalibSuffix_).c_str(), &saturation, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "overlapping" + CalibSuffix_).c_str(), &overlapping, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "farfromedge" + CalibSuffix_).c_str(), &farfromedge, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "charge" + CalibSuffix_).c_str(), &charge, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "path" + CalibSuffix_).c_str(), &path, nullptr);
 
     tree.SetBranchAddress((CalibPrefix_ + "amplitude" + CalibSuffix_).c_str(), &amplitude, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "gainused" + CalibSuffix_).c_str(), &gainused, nullptr);
     tree.SetBranchAddress((CalibPrefix_ + "gainusedTick" + CalibSuffix_).c_str(), &gainusedTick, nullptr);
 
     int elepos = statCollectionFromMode(m_calibrationMode.c_str());
     const auto nEntries = tree.GetEntries();
     printf("Progressing Bar              :0%%       20%%       40%%       60%%       80%%       100%%\n");
     printf("Looping on the Tree          :");
     const auto treeProgStep = std::max(nEntries / 50, Long64_t(1));
     for (Long64_t i{0}; i != nEntries; ++i) {
       if ((i % treeProgStep) == 0) {
         printf(".");
         fflush(stdout);
       }
       tree.GetEntry(i);
 
       unsigned int FirstAmplitude = 0;
       for (unsigned int i = 0; i < charge->size(); i++) {
         FirstAmplitude += (*nstrips)[i];
         int TI = (*trackindex)[i];
 
         if ((*tracketa)[TI] < MinTrackEta)
           continue;
         if ((*tracketa)[TI] > MaxTrackEta)
           continue;
         if ((*trackp)[TI] < MinTrackMomentum)
           continue;
         if ((*trackp)[TI] > MaxTrackMomentum)
           continue;
         if ((*trackhitsvalid)[TI] < MinTrackHits)
           continue;
         if ((*trackchi2ndof)[TI] > MaxTrackChiOverNdf)
           continue;
         if ((*trackalgo)[TI] > MaxTrackingIteration)
           continue;
 
         std::shared_ptr<stAPVGain> APV = histograms.APVsColl.at(
             ((*rawid)[i] << 4) |
             ((*firststrip)[i] /
              128));  //works for both strip and pixel thanks to firstStrip encoding for pixel in the calibTree
 
         if (APV->SubDet > 2 && (*farfromedge)[i] == false)
           continue;
         if (APV->SubDet > 2 && (*overlapping)[i] == true)
           continue;
         if (APV->SubDet > 2 && (*saturation)[i] && !AllowSaturation)
           continue;
         if (APV->SubDet > 2 && (*nstrips)[i] > MaxNrStrips)
           continue;
 
         int Charge = 0;
         if (APV->SubDet > 2 && (useCalibration || !FirstSetOfConstants)) {
           bool Saturation = false;
           for (unsigned int s = 0; s < (*nstrips)[i]; s++) {
             int StripCharge = (*amplitude)[FirstAmplitude - (*nstrips)[i] + s];
             if (useCalibration && !FirstSetOfConstants) {
               StripCharge = (int)(StripCharge * (APV->PreviousGain / APV->CalibGain));
             } else if (useCalibration) {
               StripCharge = (int)(StripCharge / APV->CalibGain);
             } else if (!FirstSetOfConstants) {
               StripCharge = (int)(StripCharge * APV->PreviousGain);
             }
             if (StripCharge > 1024) {
               StripCharge = 255;
               Saturation = true;
             } else if (StripCharge > 254) {
               StripCharge = 254;
               Saturation = true;
             }
             Charge += StripCharge;
           }
           if (Saturation && !AllowSaturation)
             continue;
         } else if (APV->SubDet > 2) {
           Charge = (*charge)[i];
         } else {
           Charge = (*charge)[i] / 265.0;  //expected scale factor between pixel and strip charge
         }
 
         double ClusterChargeOverPath = ((double)Charge) / (*path)[i];
         if (APV->SubDet > 2) {
           if (Validation) {
             ClusterChargeOverPath /= (*gainused)[i];
           }
           if (OldGainRemoving) {
             ClusterChargeOverPath *= (*gainused)[i];
           }
         }
 
         // keep processing of pixel cluster charge until here
         if (APV->SubDet <= 2)
           continue;
 
         // real histogram for calibration
         histograms.Charge_Vs_Index[elepos]->Fill(APV->Index, ClusterChargeOverPath);
         LogDebug("SiStripGainsCalibTreeWorker")
             << " for mode " << m_calibrationMode << "\n"
             << " i " << i << " useCalibration " << useCalibration << " FirstSetOfConstants " << FirstSetOfConstants
             << " APV->PreviousGain " << APV->PreviousGain << " APV->CalibGain " << APV->CalibGain << " APV->DetId "
             << APV->DetId << " APV->Index " << APV->Index << " Charge " << Charge << " Path " << (*path)[i]
             << " ClusterChargeOverPath " << ClusterChargeOverPath << std::endl;
 
         // Fill monitoring histograms
         int mCharge1 = 0;
         int mCharge2 = 0;
         int mCharge3 = 0;
         int mCharge4 = 0;
         if (APV->SubDet > 2) {
           for (unsigned int s = 0; s < (*nstrips)[i]; s++) {
             int StripCharge = (*amplitude)[FirstAmplitude - (*nstrips)[i] + s];
             if (StripCharge > 1024)
               StripCharge = 255;
             else if (StripCharge > 254)
               StripCharge = 254;
             mCharge1 += StripCharge;
             mCharge2 += StripCharge;
             mCharge3 += StripCharge;
             mCharge4 += StripCharge;
           }
           // Revome gains for monitoring
           mCharge2 *= (*gainused)[i];                         // remove G2
           mCharge3 *= (*gainusedTick)[i];                     // remove G1
           mCharge4 *= ((*gainused)[i] * (*gainusedTick)[i]);  // remove G1 and G2
         }
 
         LogDebug("SiStripGainsCalibTreeWorker")
             << " full charge " << mCharge1 << " remove G2 " << mCharge2 << " remove G1 " << mCharge3 << " remove G1*G2 "
             << mCharge4 << std::endl;
 
         auto indices = APVGain::FetchIndices(theTopologyMap, (*rawid)[i], topo);
 
         for (auto m : indices)
           histograms.Charge_1[elepos][m]->Fill(((double)mCharge1) / (*path)[i]);
         for (auto m : indices)
           histograms.Charge_2[elepos][m]->Fill(((double)mCharge2) / (*path)[i]);
         for (auto m : indices)
           histograms.Charge_3[elepos][m]->Fill(((double)mCharge3) / (*path)[i]);
         for (auto m : indices)
           histograms.Charge_4[elepos][m]->Fill(((double)mCharge4) / (*path)[i]);
 
         if (APV->SubDet == StripSubdetector::TIB) {
           histograms.Charge_Vs_PathlengthTIB[elepos]->Fill((*path)[i], Charge);  // TIB
 
         } else if (APV->SubDet == StripSubdetector::TOB) {
           histograms.Charge_Vs_PathlengthTOB[elepos]->Fill((*path)[i], Charge);  // TOB
 
         } else if (APV->SubDet == StripSubdetector::TID) {
           if (APV->Eta < 0) {
             histograms.Charge_Vs_PathlengthTIDM[elepos]->Fill((*path)[i], Charge);
           }  // TID minus
           else if (APV->Eta > 0) {
             histograms.Charge_Vs_PathlengthTIDP[elepos]->Fill((*path)[i], Charge);
           }  // TID plus
 
         } else if (APV->SubDet == StripSubdetector::TEC) {
           if (APV->Eta < 0) {
             if (APV->Thickness < 0.04) {
               histograms.Charge_Vs_PathlengthTECM1[elepos]->Fill((*path)[i], Charge);
             }  // TEC minus, type 1
             else if (APV->Thickness > 0.04) {
               histograms.Charge_Vs_PathlengthTECM2[elepos]->Fill((*path)[i], Charge);
             }  // TEC minus, type 2
           } else if (APV->Eta > 0) {
             if (APV->Thickness < 0.04) {
               histograms.Charge_Vs_PathlengthTECP1[elepos]->Fill((*path)[i], Charge);
             }  // TEC plus, type 1
             else if (APV->Thickness > 0.04) {
               histograms.Charge_Vs_PathlengthTECP2[elepos]->Fill((*path)[i], Charge);
             }  // TEC plus, type 2
           }
         }
 
       }  // END OF ON-CLUSTER LOOP
 
       //LogDebug("SiStripGainsCalibTreeWorker")<<" for mode"<< m_calibrationMode
       //                   <<" entries in histogram:"<< histograms.Charge_Vs_Index[elepos].getEntries()
 
       //                   <<std::endl;
     }
 
     hasProcessed_ = true;
   }
 }

void SiStripGainsCalibTreeWorker::dqmBeginRun	(	edm::Run const &	run,
		edm::EventSetup const &	iSetup,
		APVGain::APVGainHistograms &	histograms
	)		const

overrideprivatevirtual

Reimplemented from DQMGlobalEDAnalyzerBase< APVGain::APVGainHistograms, Args...>.

Definition at line 193 of file SiStripGainsCalibTreeWorker.cc.

References a, sistrip::APV, APVGain::APVGainHistograms::APVsCollOrdered, checkBookAPVColls(), beamvalidation::exit(), gainToken_, edm::EventSetup::getData(), edm::EventSetup::getHandle(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, qualityToken_, submitPVResolutionJobs::stdout, and tkGeomTokenBR_.

                                                                                           {
   using namespace edm;
 
   // fills the APV collections at each begin run
   const TrackerGeometry* bareTkGeomPtr = &iSetup.getData(tkGeomTokenBR_);
   checkBookAPVColls(bareTkGeomPtr, histograms);
 
   const auto gainHandle = iSetup.getHandle(gainToken_);
   if (!gainHandle.isValid()) {
     edm::LogError("SiStripGainPCLWorker") << "gainHandle is not valid\n";
     exit(0);
   }
 
   const auto& SiStripQuality_ = &iSetup.getData(qualityToken_);
 
   for (unsigned int a = 0; a < histograms.APVsCollOrdered.size(); a++) {
     std::shared_ptr<stAPVGain> APV = histograms.APVsCollOrdered[a];
 
     if (APV->SubDet == PixelSubdetector::PixelBarrel || APV->SubDet == PixelSubdetector::PixelEndcap)
       continue;
 
     APV->isMasked = SiStripQuality_->IsApvBad(APV->DetId, APV->APVId);
 
     if (gainHandle->getNumberOfTags() != 2) {
       edm::LogError("SiStripGainPCLWorker") << "NUMBER OF GAIN TAG IS EXPECTED TO BE 2\n";
       fflush(stdout);
       exit(0);
     };
     float newPreviousGain = gainHandle->getApvGain(APV->APVId, gainHandle->getRange(APV->DetId, 1), 1);
     if (APV->PreviousGain != 1 and newPreviousGain != APV->PreviousGain)
       edm::LogWarning("SiStripGainPCLWorker") << "WARNING: ParticleGain in the global tag changed\n";
     APV->PreviousGain = newPreviousGain;
 
     float newPreviousGainTick = gainHandle->getApvGain(APV->APVId, gainHandle->getRange(APV->DetId, 0), 0);
     if (APV->PreviousGainTick != 1 and newPreviousGainTick != APV->PreviousGainTick) {
       edm::LogWarning("SiStripGainPCLWorker")
           << "WARNING: TickMarkGain in the global tag changed\n"
           << std::endl
           << " APV->SubDet: " << APV->SubDet << " APV->APVId:" << APV->APVId << std::endl
           << " APV->PreviousGainTick: " << APV->PreviousGainTick << " newPreviousGainTick: " << newPreviousGainTick
           << std::endl;
     }
     APV->PreviousGainTick = newPreviousGainTick;
   }
 }

void SiStripGainsCalibTreeWorker::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 615 of file SiStripGainsCalibTreeWorker.cc.

References edm::ConfigurationDescriptions::addDefault(), submitPVResolutionJobs::desc, and edm::ParameterSetDescription::setUnknown().

                                                                                              {
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

int SiStripGainsCalibTreeWorker::statCollectionFromMode ( const char * tag ) const

inlineprivate

Definition at line 109 of file SiStripGainsCalibTreeWorker.cc.

References dqm_tag_, relativeConstraints::empty, None, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by bookHistograms(), and dqmAnalyze().

                                                                                     {
   std::vector<std::string>::const_iterator it = dqm_tag_.begin();
   while (it != dqm_tag_.end()) {
     if (*it == std::string(tag))
       return it - dqm_tag_.begin();
     it++;
   }
 
   if (std::string(tag).empty())
     return 0;  // return StdBunch calibration mode for backward compatibility
 
   return None;
 }