Inheritance diagram for SiStripGainFromCalibTree:

Public Types
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	SiStripGainFromCalibTree (const edm::ParameterSet &)

	~SiStripGainFromCalibTree () override

Public Member Functions inherited from ConditionDBWriter< SiStripApvGain >
	ConditionDBWriter (const edm::ParameterSet &iConfig)

	~ConditionDBWriter () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () 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)

Private Member Functions
void	algoAnalyze (const edm::Event &, const edm::EventSetup &) override

void	algoAnalyzeTheTree ()

void	algoBeginJob (const edm::EventSetup &iSetup) override

void	algoBeginRun (const edm::Run &run, const edm::EventSetup &iSetup) override

void	algoComputeMPVandGain ()

void	algoEndJob () override

void	algoEndRun (const edm::Run &run, const edm::EventSetup &iSetup) override

void	bookDQMHistos (const char dqm_dir, const char tag)

template<typename T >
edm::Handle< T >	connect (const T *&ptr, edm::EDGetTokenT< T > token, const edm::Event &evt)

std::unique_ptr< SiStripApvGain >	getNewObject () override

void	getPeakOfLandau (TH1 InputHisto, double FitResults, double LowRange=50, double HighRange=5400)

bool	isBFieldConsistentWithMode (const edm::EventSetup &iSetup) const

bool	IsGoodLandauFit (double *FitResults)

void	MakeCalibrationMap ()

void	merge (TH2 A, TH2 B)

void	processEvent ()

bool	produceTagFilter ()

void	qualityMonitor ()

int	statCollectionFromMode (const char *tag) const

void	storeOnTree (TFileService *tfs)

void	swapBFieldMode (void)

Private Attributes
std::string	AlgoMode

bool	AllowSaturation

const std::vector< unsigned char > *	amplitude = nullptr

edm::EDGetTokenT< std::vector< unsigned char > >	amplitude_token_

std::unordered_map< unsigned int, stAPVGain * >	APVsColl

std::vector< stAPVGain * >	APVsCollOrdered

unsigned int	BAD

std::string	booked_dir_

string	CalibPrefix_

int	CalibrationLevel

string	CalibSuffix_

const std::vector< unsigned int > *	charge = nullptr

std::array< std::vector< APVGain::APVmon >, 7 >	Charge_1

std::array< std::vector< APVGain::APVmon >, 7 >	Charge_2

std::array< std::vector< APVGain::APVmon >, 7 >	Charge_3

std::array< std::vector< APVGain::APVmon >, 7 >	Charge_4

edm::EDGetTokenT< std::vector< unsigned int > >	charge_token_

std::vector< MonitorElement * >	Charge_Vs_Index

std::vector< MonitorElement * >	Charge_Vs_PathlengthTECM1

std::vector< MonitorElement * >	Charge_Vs_PathlengthTECM2

std::vector< MonitorElement * >	Charge_Vs_PathlengthTECP1

std::vector< MonitorElement * >	Charge_Vs_PathlengthTECP2

std::vector< MonitorElement * >	Charge_Vs_PathlengthTIB

std::vector< MonitorElement * >	Charge_Vs_PathlengthTIDM

std::vector< MonitorElement * >	Charge_Vs_PathlengthTIDP

std::vector< MonitorElement * >	Charge_Vs_PathlengthTOB

const std::vector< double > *	chargeoverpath = nullptr

edm::EDGetTokenT< std::vector< double > >	chargeoverpath_token_

DQMStore *	dbe

MonitorElement *	DiffWRTPrevGainTEC

MonitorElement *	DiffWRTPrevGainTIB

MonitorElement *	DiffWRTPrevGainTID

MonitorElement *	DiffWRTPrevGainTOB

bool	doChargeMonitorPerPlane

std::vector< string >	dqm_tag_

unsigned int	ERun

unsigned int	eventnumber = 0

string	EventPrefix_

string	EventSuffix_

const std::vector< bool > *	farfromedge = nullptr

edm::EDGetTokenT< std::vector< bool > >	farfromedge_token_

bool	FirstSetOfConstants

const std::vector< unsigned short > *	firststrip = nullptr

edm::EDGetTokenT< std::vector< unsigned short > >	firststrip_token_

MonitorElement *	Gains

edm::ESGetToken< SiStripGain, SiStripGainRcd >	gainToken_

const std::vector< double > *	gainused = nullptr

edm::EDGetTokenT< std::vector< double > >	gainused_token_

const std::vector< double > *	gainusedTick = nullptr

edm::EDGetTokenT< std::vector< double > >	gainusedTick_token_

MonitorElement *	GainVsPrevGainTEC

MonitorElement *	GainVsPrevGainTIB

MonitorElement *	GainVsPrevGainTID

MonitorElement *	GainVsPrevGainTOB

unsigned int	GOOD

const std::vector< double > *	localdirx = nullptr

edm::EDGetTokenT< std::vector< double > >	localdirx_token_

const std::vector< double > *	localdiry = nullptr

edm::EDGetTokenT< std::vector< double > >	localdiry_token_

const std::vector< double > *	localdirz = nullptr

edm::EDGetTokenT< std::vector< double > >	localdirz_token_

std::string	m_calibrationMode

std::string	m_calibrationPath

std::string	m_DQMdir

bool	m_harvestingMode

bool	m_splitDQMstat

double	MagFieldCurrentTh

unsigned int	MASKED

double	MaxChi2OverNDF

double	MaxMPVError

unsigned int	MaxNrStrips

double	MaxTrackChiOverNdf

double	MaxTrackEta

int	MaxTrackingIteration

double	MaxTrackMomentum

double	MinNrEntries

double	MinTrackEta

unsigned int	MinTrackHits

double	MinTrackMomentum

MonitorElement *	MPV_Vs_EtaTEC

MonitorElement *	MPV_Vs_EtaTECthick

MonitorElement *	MPV_Vs_EtaTECthin

MonitorElement *	MPV_Vs_EtaTIB

MonitorElement *	MPV_Vs_EtaTID

MonitorElement *	MPV_Vs_EtaTOB

MonitorElement *	MPV_Vs_PhiTEC

MonitorElement *	MPV_Vs_PhiTECthick

MonitorElement *	MPV_Vs_PhiTECthin

MonitorElement *	MPV_Vs_PhiTIB

MonitorElement *	MPV_Vs_PhiTID

MonitorElement *	MPV_Vs_PhiTOB

MonitorElement *	MPVError

MonitorElement *	MPVErrorVsEta

MonitorElement *	MPVErrorVsMPV

MonitorElement *	MPVErrorVsN

MonitorElement *	MPVErrorVsPhi

MonitorElement *	MPVs

MonitorElement *	MPVs320

MonitorElement *	MPVs500

MonitorElement *	MPVsTEC

MonitorElement *	MPVsTECM

MonitorElement *	MPVsTECM1

MonitorElement *	MPVsTECM2

MonitorElement *	MPVsTECP

MonitorElement *	MPVsTECP1

MonitorElement *	MPVsTECP2

MonitorElement *	MPVsTECthick

MonitorElement *	MPVsTECthin

MonitorElement *	MPVsTIB

MonitorElement *	MPVsTID

MonitorElement *	MPVsTIDM

MonitorElement *	MPVsTIDP

MonitorElement *	MPVsTOB

unsigned int	NClusterPixel

unsigned int	NClusterStrip

unsigned int	NEvent

std::vector< APVGain::APVmon >	newCharge

MonitorElement *	NoMPVfit

MonitorElement *	NoMPVmasked

int	NPixelDets

int	NStripAPVs

const std::vector< unsigned short > *	nstrips = nullptr

edm::EDGetTokenT< std::vector< unsigned short > >	nstrips_token_

unsigned int	NTrack

bool	OldGainRemoving

std::string	OutputGains

const std::vector< bool > *	overlapping = nullptr

edm::EDGetTokenT< std::vector< bool > >	overlapping_token_

const std::vector< double > *	path = nullptr

edm::EDGetTokenT< std::vector< double > >	path_token_

edm::ESGetToken< SiStripQuality, SiStripQualityRcd >	qualityToken_

const std::vector< unsigned int > *	rawid = nullptr

edm::EDGetTokenT< std::vector< unsigned int > >	rawid_token_

edm::ESGetToken< RunInfo, RunInfoRcd >	runInfoToken_

unsigned int	runnumber = 0

const std::vector< bool > *	saturation = nullptr

edm::EDGetTokenT< std::vector< bool > >	saturation_token_

bool	saveSummary

unsigned int	SRun

double	tagCondition_GoodFrac

double	tagCondition_NClusters

TFileService *	tfs

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_

const std::vector< int > *	trackalgo = nullptr

edm::EDGetTokenT< std::vector< int > >	trackalgo_token_

const std::vector< double > *	trackchi2ndof = nullptr

edm::EDGetTokenT< std::vector< double > >	trackchi2ndof_token_

const std::vector< double > *	tracketa = nullptr

edm::EDGetTokenT< std::vector< double > >	tracketa_token_

const std::vector< unsigned int > *	trackhitsvalid = nullptr

edm::EDGetTokenT< std::vector< unsigned int > >	trackhitsvalid_token_

const std::vector< int > *	trackindex = nullptr

edm::EDGetTokenT< std::vector< int > >	trackindex_token_

const std::vector< float > *	trackp = nullptr

edm::EDGetTokenT< std::vector< float > >	trackp_token_

const std::vector< double > *	trackphi = nullptr

edm::EDGetTokenT< std::vector< double > >	trackphi_token_

string	TrackPrefix_

const std::vector< float > *	trackpt = nullptr

edm::EDGetTokenT< std::vector< float > >	trackpt_token_

string	TrackSuffix_

const std::vector< bool > *	TrigTech = nullptr

edm::EDGetTokenT< std::vector< bool > >	TrigTech_token_

const TrackerTopology *	tTopo_ = nullptr

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

bool	useCalibration

bool	Validation

std::vector< std::string >	VChargeHisto

vector< string >	VInputFiles

Additional Inherited Members
Static Public Member Functions inherited from ConditionDBWriter< SiStripApvGain >
static void	fillPSetDescription (edm::ParameterSetDescription &desc)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from ConditionDBWriter< SiStripApvGain >
void	setDoStore (const bool doStore)
	When set to false the payload will not be written to the db. More...

void	storeOnDbNow ()

cond::Time_t	timeOfLastIOV ()

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)

Detailed Description

Definition at line 88 of file SiStripGainFromCalibTree.cc.

Member Typedef Documentation

◆ DQMStore

typedef dqm::legacy::DQMStore SiStripGainFromCalibTree::DQMStore

Definition at line 91 of file SiStripGainFromCalibTree.cc.

◆ MonitorElement

typedef dqm::legacy::MonitorElement SiStripGainFromCalibTree::MonitorElement

Definition at line 90 of file SiStripGainFromCalibTree.cc.

Constructor & Destructor Documentation

◆ SiStripGainFromCalibTree()

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

explicit

Definition at line 359 of file SiStripGainFromCalibTree.cc.

     : ConditionDBWriter<SiStripApvGain>(iConfig) {
   OutputGains = iConfig.getParameter<std::string>("OutputGains");
 
   AlgoMode = iConfig.getUntrackedParameter<std::string>("AlgoMode", "CalibTree");
   MagFieldCurrentTh = iConfig.getUntrackedParameter<double>("MagFieldCurrentTh", 2000.);
   MinNrEntries = iConfig.getUntrackedParameter<double>("minNrEntries", 20);
   MaxMPVError = iConfig.getUntrackedParameter<double>("maxMPVError", 500.0);
   MaxChi2OverNDF = iConfig.getUntrackedParameter<double>("maxChi2OverNDF", 5.0);
   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);
 
   CalibrationLevel = iConfig.getUntrackedParameter<int>("CalibrationLevel", 0);
   VInputFiles = iConfig.getUntrackedParameter<vector<string>>("InputFiles");
   VChargeHisto = iConfig.getUntrackedParameter<vector<string>>("ChargeHisto");
 
   useCalibration = iConfig.getUntrackedParameter<bool>("UseCalibration", false);
   m_harvestingMode = iConfig.getUntrackedParameter<bool>("harvestingMode", false);
   m_splitDQMstat = iConfig.getUntrackedParameter<bool>("splitDQMstat", false);
   m_calibrationMode = iConfig.getUntrackedParameter<string>("calibrationMode", "StdBunch");
   m_calibrationPath = iConfig.getUntrackedParameter<string>("calibrationPath");
   m_DQMdir = iConfig.getUntrackedParameter<string>("DQMdir", "AlCaReco/SiStripGains");
 
   tagCondition_NClusters = iConfig.getUntrackedParameter<double>("NClustersForTagProd", 2E8);
   tagCondition_GoodFrac = iConfig.getUntrackedParameter<double>("GoodFracForTagProd", 0.95);
 
   saveSummary = iConfig.getUntrackedParameter<bool>("saveSummary", false);
 
   doChargeMonitorPerPlane = iConfig.getUntrackedParameter<bool>("doChargeMonitorPerPlane", false);
 
   // Gather DQM Service
   dbe = edm::Service<DQMStore>().operator->();
 
   //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
 
   Charge_Vs_Index.insert(Charge_Vs_Index.begin(), dqm_tag_.size(), nullptr);
   //Charge_Vs_Index_Absolute.insert( Charge_Vs_Index_Absolute.begin(), dqm_tag_.size(), 0);
   Charge_Vs_PathlengthTIB.insert(Charge_Vs_PathlengthTIB.begin(), dqm_tag_.size(), nullptr);
   Charge_Vs_PathlengthTOB.insert(Charge_Vs_PathlengthTOB.begin(), dqm_tag_.size(), nullptr);
   Charge_Vs_PathlengthTIDP.insert(Charge_Vs_PathlengthTIDP.begin(), dqm_tag_.size(), nullptr);
   Charge_Vs_PathlengthTIDM.insert(Charge_Vs_PathlengthTIDM.begin(), dqm_tag_.size(), nullptr);
   Charge_Vs_PathlengthTECP1.insert(Charge_Vs_PathlengthTECP1.begin(), dqm_tag_.size(), nullptr);
   Charge_Vs_PathlengthTECP2.insert(Charge_Vs_PathlengthTECP2.begin(), dqm_tag_.size(), nullptr);
   Charge_Vs_PathlengthTECM1.insert(Charge_Vs_PathlengthTECM1.begin(), dqm_tag_.size(), nullptr);
   Charge_Vs_PathlengthTECM2.insert(Charge_Vs_PathlengthTECM2.begin(), dqm_tag_.size(), nullptr);
 
   // configure token for gathering the ntuple variables
   edm::ParameterSet swhallowgain_pset = iConfig.getUntrackedParameter<edm::ParameterSet>("gain");
 
   string label = swhallowgain_pset.getUntrackedParameter<string>("label");
   CalibPrefix_ = swhallowgain_pset.getUntrackedParameter<string>("prefix");
   CalibSuffix_ = swhallowgain_pset.getUntrackedParameter<string>("suffix");
 
   trackindex_token_ = consumes<std::vector<int>>(edm::InputTag(label, CalibPrefix_ + "trackindex" + CalibSuffix_));
   rawid_token_ = consumes<std::vector<unsigned int>>(edm::InputTag(label, CalibPrefix_ + "rawid" + CalibSuffix_));
   localdirx_token_ = consumes<std::vector<double>>(edm::InputTag(label, CalibPrefix_ + "localdirx" + CalibSuffix_));
   localdiry_token_ = consumes<std::vector<double>>(edm::InputTag(label, CalibPrefix_ + "localdiry" + CalibSuffix_));
   localdirz_token_ = consumes<std::vector<double>>(edm::InputTag(label, CalibPrefix_ + "localdirz" + CalibSuffix_));
   firststrip_token_ =
       consumes<std::vector<unsigned short>>(edm::InputTag(label, CalibPrefix_ + "firststrip" + CalibSuffix_));
   nstrips_token_ = consumes<std::vector<unsigned short>>(edm::InputTag(label, CalibPrefix_ + "nstrips" + CalibSuffix_));
   saturation_token_ = consumes<std::vector<bool>>(edm::InputTag(label, CalibPrefix_ + "saturation" + CalibSuffix_));
   overlapping_token_ = consumes<std::vector<bool>>(edm::InputTag(label, CalibPrefix_ + "overlapping" + CalibSuffix_));
   farfromedge_token_ = consumes<std::vector<bool>>(edm::InputTag(label, CalibPrefix_ + "farfromedge" + CalibSuffix_));
   charge_token_ = consumes<std::vector<unsigned int>>(edm::InputTag(label, CalibPrefix_ + "charge" + CalibSuffix_));
   path_token_ = consumes<std::vector<double>>(edm::InputTag(label, CalibPrefix_ + "path" + CalibSuffix_));
   chargeoverpath_token_ =
       consumes<std::vector<double>>(edm::InputTag(label, CalibPrefix_ + "chargeoverpath" + CalibSuffix_));
   amplitude_token_ =
       consumes<std::vector<unsigned char>>(edm::InputTag(label, CalibPrefix_ + "amplitude" + CalibSuffix_));
   gainused_token_ = consumes<std::vector<double>>(edm::InputTag(label, CalibPrefix_ + "gainused" + CalibSuffix_));
   gainusedTick_token_ =
       consumes<std::vector<double>>(edm::InputTag(label, CalibPrefix_ + "gainusedTick" + CalibSuffix_));
 
   edm::ParameterSet evtinfo_pset = iConfig.getUntrackedParameter<edm::ParameterSet>("evtinfo");
   label = evtinfo_pset.getUntrackedParameter<string>("label");
   EventPrefix_ = evtinfo_pset.getUntrackedParameter<string>("prefix");
   EventSuffix_ = evtinfo_pset.getUntrackedParameter<string>("suffix");
   TrigTech_token_ = consumes<std::vector<bool>>(edm::InputTag(label, EventPrefix_ + "TrigTech" + EventSuffix_));
 
   edm::ParameterSet track_pset = iConfig.getUntrackedParameter<edm::ParameterSet>("tracks");
   label = track_pset.getUntrackedParameter<string>("label");
   TrackPrefix_ = track_pset.getUntrackedParameter<string>("prefix");
   TrackSuffix_ = track_pset.getUntrackedParameter<string>("suffix");
 
   trackchi2ndof_token_ = consumes<std::vector<double>>(edm::InputTag(label, TrackPrefix_ + "chi2ndof" + TrackSuffix_));
   trackp_token_ = consumes<std::vector<float>>(edm::InputTag(label, TrackPrefix_ + "momentum" + TrackSuffix_));
   trackpt_token_ = consumes<std::vector<float>>(edm::InputTag(label, TrackPrefix_ + "pt" + TrackSuffix_));
   tracketa_token_ = consumes<std::vector<double>>(edm::InputTag(label, TrackPrefix_ + "eta" + TrackSuffix_));
   trackphi_token_ = consumes<std::vector<double>>(edm::InputTag(label, TrackPrefix_ + "phi" + TrackSuffix_));
   trackhitsvalid_token_ =
       consumes<std::vector<unsigned int>>(edm::InputTag(label, TrackPrefix_ + "hitsvalid" + TrackSuffix_));
   trackalgo_token_ = consumes<std::vector<int>>(edm::InputTag(label, TrackPrefix_ + "algo" + TrackSuffix_));
 
   tTopoToken_ = esConsumes<edm::Transition::BeginRun>();
   tkGeomToken_ = esConsumes<edm::Transition::BeginRun>();
   runInfoToken_ = esConsumes<edm::Transition::BeginRun>();
   gainToken_ = esConsumes<edm::Transition::BeginRun>();
   qualityToken_ = esConsumes<edm::Transition::BeginRun>();
 }

◆ ~SiStripGainFromCalibTree()

SiStripGainFromCalibTree::~SiStripGainFromCalibTree ( )

override

Definition at line 1807 of file SiStripGainFromCalibTree.cc.

References a, sistrip::APV, APVsColl, and APVsCollOrdered.

                                                     {
   APVsColl.clear();
   for (unsigned int a = 0; a < APVsCollOrdered.size(); a++) {
     stAPVGain* APV = APVsCollOrdered[a];
     if (APV != nullptr)
       delete APV;
   }
   APVsCollOrdered.clear();
 }

Member Function Documentation

◆ algoAnalyze()

void SiStripGainFromCalibTree::algoAnalyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 1841 of file SiStripGainFromCalibTree.cc.

References AlgoMode, amplitude, amplitude_token_, charge, charge_token_, chargeoverpath, chargeoverpath_token_, connect(), eventnumber, farfromedge, farfromedge_token_, firststrip, firststrip_token_, gainused, gainused_token_, gainusedTick, gainusedTick_token_, iEvent, localdirx, localdirx_token_, localdiry, localdiry_token_, localdirz, localdirz_token_, m_harvestingMode, nstrips, nstrips_token_, overlapping, overlapping_token_, path, path_token_, processEvent(), rawid, rawid_token_, runnumber, saturation, saturation_token_, trackalgo, trackalgo_token_, trackchi2ndof, trackchi2ndof_token_, tracketa, tracketa_token_, trackhitsvalid, trackhitsvalid_token_, trackindex, trackindex_token_, trackp, trackp_token_, trackphi, trackphi_token_, trackpt, trackpt_token_, TrigTech, and TrigTech_token_.

                                                                                               {
   // in AlCaHarvesting mode we just need to run the logic in the endJob step
   if (m_harvestingMode)
     return;
 
   if (AlgoMode == "CalibTree")
     return;
 
   eventnumber = iEvent.id().event();
   runnumber = iEvent.id().run();
   auto handle01 = connect(TrigTech, TrigTech_token_, iEvent);
   auto handle02 = connect(trackchi2ndof, trackchi2ndof_token_, iEvent);
   auto handle03 = connect(trackp, trackp_token_, iEvent);
   auto handle04 = connect(trackpt, trackpt_token_, iEvent);
   auto handle05 = connect(tracketa, tracketa_token_, iEvent);
   auto handle06 = connect(trackphi, trackphi_token_, iEvent);
   auto handle07 = connect(trackhitsvalid, trackhitsvalid_token_, iEvent);
   auto handle08 = connect(trackindex, trackindex_token_, iEvent);
   auto handle09 = connect(rawid, rawid_token_, iEvent);
   auto handle11 = connect(localdirx, localdirx_token_, iEvent);
   auto handle12 = connect(localdiry, localdiry_token_, iEvent);
   auto handle13 = connect(localdirz, localdirz_token_, iEvent);
   auto handle14 = connect(firststrip, firststrip_token_, iEvent);
   auto handle15 = connect(nstrips, nstrips_token_, iEvent);
   auto handle16 = connect(saturation, saturation_token_, iEvent);
   auto handle17 = connect(overlapping, overlapping_token_, iEvent);
   auto handle18 = connect(farfromedge, farfromedge_token_, iEvent);
   auto handle19 = connect(charge, charge_token_, iEvent);
   auto handle21 = connect(path, path_token_, iEvent);
   auto handle22 = connect(chargeoverpath, chargeoverpath_token_, iEvent);
   auto handle23 = connect(amplitude, amplitude_token_, iEvent);
   auto handle24 = connect(gainused, gainused_token_, iEvent);
   auto handle25 = connect(gainusedTick, gainusedTick_token_, iEvent);
   auto handle26 = connect(trackalgo, trackalgo_token_, iEvent);
 
   processEvent();
 }

◆ algoAnalyzeTheTree()

void SiStripGainFromCalibTree::algoAnalyzeTheTree ( )

private

Definition at line 1311 of file SiStripGainFromCalibTree.cc.

References amplitude, CalibPrefix_, CalibSuffix_, charge, chargeoverpath, eventnumber, EventPrefix_, EventSuffix_, farfromedge, firststrip, gainused, gainusedTick, mps_fire::i, localdirx, localdiry, localdirz, m_calibrationMode, NEvent, nstrips, overlapping, path, processEvent(), rawid, runnumber, saturation, mps_setup::stdout, trackalgo, trackchi2ndof, tracketa, trackhitsvalid, trackindex, trackp, trackphi, TrackPrefix_, trackpt, TrackSuffix_, TrigTech, and VInputFiles.

Referenced by algoEndJob().

                                                   {
   for (unsigned int i = 0; i < VInputFiles.size(); i++) {
     printf("Openning file %3i/%3i --> %s\n", i + 1, (int)VInputFiles.size(), (char*)(VInputFiles[i].c_str()));
     fflush(stdout);
     TFile* tfile = TFile::Open(VInputFiles[i].c_str());
     TString tree_path = TString::Format("gainCalibrationTree%s/tree", m_calibrationMode.c_str());
     TTree* tree = dynamic_cast<TTree*>(tfile->Get(tree_path.Data()));
 
     tree->SetBranchAddress((EventPrefix_ + "event" + EventSuffix_).c_str(), &eventnumber, nullptr);
     tree->SetBranchAddress((EventPrefix_ + "run" + EventSuffix_).c_str(), &runnumber, nullptr);
     tree->SetBranchAddress((EventPrefix_ + "TrigTech" + EventSuffix_).c_str(), &TrigTech, nullptr);
 
     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_ + "chargeoverpath" + CalibSuffix_).c_str(), &chargeoverpath, 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);
 
     unsigned int nentries = tree->GetEntries();
     printf("Number of Events = %i + %i = %i\n", NEvent, nentries, (NEvent + nentries));
     printf("Progressing Bar              :0%%       20%%       40%%       60%%       80%%       100%%\n");
     printf("Looping on the Tree          :");
     int TreeStep = nentries / 50;
     if (TreeStep <= 1)
       TreeStep = 1;
     for (unsigned int ientry = 0; ientry < tree->GetEntries(); ientry++) {
       if (ientry % TreeStep == 0) {
         printf(".");
         fflush(stdout);
       }
       tree->GetEntry(ientry);
       processEvent();
     }
     printf("\n");  // END OF EVENT LOOP
   }
 }

◆ algoBeginJob()

void SiStripGainFromCalibTree::algoBeginJob ( const edm::EventSetup & iSetup )

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 661 of file SiStripGainFromCalibTree.cc.

References AlgoMode, sistrip::APV, APVsColl, APVsCollOrdered, BAD, bookDQMHistos(), dqm_tag_, ERun, edm::EventSetup::getData(), GOOD, mps_fire::i, isBFieldConsistentWithMode(), dqmiolumiharvest::j, m_calibrationMode, m_DQMdir, m_harvestingMode, m_splitDQMstat, MakeCalibrationMap(), MASKED, NClusterPixel, NClusterStrip, PixelTopology::ncolumns(), NEvent, NPixelDets, PixelTopology::nrows(), NStripAPVs, StripTopology::nstrips(), NTrack, or, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::rawId(), SRun, statCollectionFromMode(), AlCaHLTBitMon_QueryRunRegistry::string, ntupleEnum::SubDet, DetId::subdetId(), swapBFieldMode(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, tkGeomToken_, StripSubdetector::TOB, tTopo_, and tTopoToken_.

                                                                        {
   edm::LogInfo("SiStripGainFromCalibTree") << "AlgoMode        : " << AlgoMode << "\n"
                                            << "CalibrationMode : " << m_calibrationMode << "\n"
                                            << "HarvestingMode  : " << m_harvestingMode << std::endl;
   //Setup DQM histograms
   if (AlgoMode != "PCL" or m_harvestingMode) {
     const char* dqm_dir = "AlCaReco/SiStripGainsHarvesting/";
     this->bookDQMHistos(dqm_dir, dqm_tag_[statCollectionFromMode(m_calibrationMode.c_str())].c_str());
   } else {
     //Check consistency of calibration Mode and BField only for the ALCAPROMPT in the PCL workflow
     if (!isBFieldConsistentWithMode(iSetup)) {
       string prevMode = m_calibrationMode;
       swapBFieldMode();
       edm::LogInfo("SiStripGainFromCalibTree") << "Switching calibration mode for endorsing BField status: " << prevMode
                                                << " ==> " << m_calibrationMode << std::endl;
     }
     std::string dqm_dir = m_DQMdir + ((m_splitDQMstat) ? m_calibrationMode : "") + "/";
     int elem = statCollectionFromMode(m_calibrationMode.c_str());
     this->bookDQMHistos(dqm_dir.c_str(), dqm_tag_[elem].c_str());
   }
 
   tTopo_ = &iSetup.getData(tTopoToken_);
 
   auto const& Det = iSetup.getData(tkGeomToken_).dets();
 
   NPixelDets = 0;
   NStripAPVs = 0;
   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++) {
         stAPVGain* APV = new 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;
 
         APVsCollOrdered.push_back(APV);
         APVsColl[(APV->DetId << 4) | APV->APVId] = APV;
         Index++;
         NStripAPVs++;
       }
     }
   }
 
   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++) {
           stAPVGain* APV = new 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;
 
           APVsCollOrdered.push_back(APV);
           APVsColl[(APV->DetId << 4) | APV->APVId] = APV;
           Index++;
           NPixelDets++;
         }
       }
     }
   }
 
   MakeCalibrationMap();
 
   NEvent = 0;
   NTrack = 0;
   NClusterStrip = 0;
   NClusterPixel = 0;
   SRun = 1 << 31;
   ERun = 0;
   GOOD = 0;
   BAD = 0;
   MASKED = 0;
 }

◆ algoBeginRun()

void SiStripGainFromCalibTree::algoBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 815 of file SiStripGainFromCalibTree.cc.

References a, AlgoMode, sistrip::APV, APVsCollOrdered, beamvalidation::exit(), gainToken_, edm::EventSetup::getData(), edm::EventSetup::getHandle(), isBFieldConsistentWithMode(), m_calibrationMode, m_harvestingMode, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, qualityToken_, mps_setup::stdout, and swapBFieldMode().

                                                                                           {
   if (!m_harvestingMode && AlgoMode == "PCL") {
     //Check consistency of calibration Mode and BField only for the ALCAPROMPT in the PCL workflow
     if (!isBFieldConsistentWithMode(iSetup)) {
       string prevMode = m_calibrationMode;
       swapBFieldMode();
       edm::LogInfo("SiStripGainFromCalibTree") << "Switching calibration mode for endorsing BField status: " << prevMode
                                                << " ==> " << m_calibrationMode << std::endl;
     }
   }
 
   const auto gainHandle = iSetup.getHandle(gainToken_);
   if (!gainHandle.isValid()) {
     edm::LogError("SiStripGainFromCalibTree") << "gainHandle is not valid\n";
     exit(0);
   }
 
   const auto& siStripQuality = iSetup.getData(qualityToken_);
   for (unsigned int a = 0; a < APVsCollOrdered.size(); a++) {
     stAPVGain* APV = APVsCollOrdered[a];
 
     // MM. 03/03/2017 all of this makes sense only for SiStrip (i.e. get me out of here if I am on a pixel ROC)
     if (APV->SubDet == PixelSubdetector::PixelBarrel || APV->SubDet == PixelSubdetector::PixelEndcap)
       continue;
 
     APV->isMasked = siStripQuality.IsApvBad(APV->DetId, APV->APVId);
     //      if(!FirstSetOfConstants){
 
     if (gainHandle->getNumberOfTags() != 2) {
       edm::LogError("SiStripGainFromCalibTree") << "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("SiStripGainFromCalibTree") << "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("SiStripGainFromCalibTree")
           << "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;
 
     //printf("DETID = %7i APVID=%1i Previous Gain=%8.4f (G1) x %8.4f (G2)\n",APV->DetId,APV->APVId,APV->PreviousGainTick, APV->PreviousGain);
     //      }
   }
 }

◆ algoComputeMPVandGain()

void SiStripGainFromCalibTree::algoComputeMPVandGain ( )

private

Definition at line 1367 of file SiStripGainFromCalibTree.cc.

References AlgoMode, sistrip::APV, stAPVGain::APVId, APVsColl, BAD, stAPVGain::Bin, CalibrationLevel, Charge_Vs_Index, stAPVGain::DetId, getPeakOfLandau(), GOOD, Harvest, Exhume::I, mps_fire::i, stAPVGain::Index, IsGoodLandauFit(), m_calibrationMode, MASKED, statCollectionFromMode(), and mps_setup::stdout.

Referenced by algoEndJob().

                                                      {
   unsigned int I = 0;
   TH1F* Proj = nullptr;
   double FitResults[6];
   double MPVmean = 300;
 
   int elepos = (AlgoMode == "PCL") ? Harvest : statCollectionFromMode(m_calibrationMode.c_str());
 
   if (Charge_Vs_Index[elepos] == nullptr) {
     edm::LogError("SiStripGainFromCalibTree")
         << "Harvesting: could not execute algoComputeMPVandGain method because " << m_calibrationMode.c_str()
         << " statistics cannot be retrieved.\n"
         << "Please check if input contains " << m_calibrationMode.c_str() << " data." << std::endl;
     return;
   }
 
   TH2S* chvsidx = (Charge_Vs_Index[elepos])->getTH2S();
 
   printf("Progressing Bar              :0%%       20%%       40%%       60%%       80%%       100%%\n");
   printf("Fitting Charge Distribution  :");
   int TreeStep = APVsColl.size() / 50;
   for (auto it = APVsColl.begin(); it != APVsColl.end(); it++, I++) {
     if (I % TreeStep == 0) {
       printf(".");
       fflush(stdout);
     }
     stAPVGain* APV = it->second;
     if (APV->Bin < 0)
       APV->Bin = chvsidx->GetXaxis()->FindBin(APV->Index);
 
     if (APV->isMasked) {
       APV->Gain = APV->PreviousGain;
       MASKED++;
       continue;
     }
 
     Proj = (TH1F*)(chvsidx->ProjectionY(
         "", chvsidx->GetXaxis()->FindBin(APV->Index), chvsidx->GetXaxis()->FindBin(APV->Index), "e"));
     if (!Proj)
       continue;
 
     if (CalibrationLevel == 0) {
     } else if (CalibrationLevel == 1) {
       int SecondAPVId = APV->APVId;
       if (SecondAPVId % 2 == 0) {
         SecondAPVId = SecondAPVId + 1;
       } else {
         SecondAPVId = SecondAPVId - 1;
       }
       stAPVGain* APV2 = APVsColl[(APV->DetId << 4) | SecondAPVId];
       if (APV2->Bin < 0)
         APV2->Bin = chvsidx->GetXaxis()->FindBin(APV2->Index);
       TH1F* Proj2 = (TH1F*)(chvsidx->ProjectionY("", APV2->Bin, APV2->Bin, "e"));
       if (Proj2) {
         Proj->Add(Proj2, 1);
         delete Proj2;
       }
     } else if (CalibrationLevel == 2) {
       for (unsigned int i = 0; i < 16; i++) {  //loop up to 6APV for Strip and up to 16 for Pixels
         auto tmpit = APVsColl.find((APV->DetId << 4) | i);
         if (tmpit == APVsColl.end())
           continue;
         stAPVGain* APV2 = tmpit->second;
         if (APV2->DetId != APV->DetId || APV2->APVId == APV->APVId)
           continue;
         if (APV2->Bin < 0)
           APV2->Bin = chvsidx->GetXaxis()->FindBin(APV2->Index);
         TH1F* Proj2 = (TH1F*)(chvsidx->ProjectionY("", APV2->Bin, APV2->Bin, "e"));
         if (Proj2) {
           Proj->Add(Proj2, 1);
           delete Proj2;
         }
       }
     } else {
       CalibrationLevel = 0;
       printf("Unknown Calibration Level, will assume %i\n", CalibrationLevel);
     }
 
     getPeakOfLandau(Proj, FitResults);
     APV->FitMPV = FitResults[0];
     APV->FitMPVErr = FitResults[1];
     APV->FitWidth = FitResults[2];
     APV->FitWidthErr = FitResults[3];
     APV->FitChi2 = FitResults[4];
     APV->FitNorm = FitResults[5];
     APV->NEntries = Proj->GetEntries();
 
     if (IsGoodLandauFit(FitResults)) {
       APV->Gain = APV->FitMPV / MPVmean;
       if (APV->SubDet > 2)
         GOOD++;
     } else {
       APV->Gain = APV->PreviousGain;
       if (APV->SubDet > 2)
         BAD++;
     }
     if (APV->Gain <= 0)
       APV->Gain = 1;
 
     //printf("%5i/%5i:  %6i - %1i  %5E Entries --> MPV = %f +- %f\n",I,APVsColl.size(),APV->DetId, APV->APVId, Proj->GetEntries(), FitResults[0], FitResults[1]);fflush(stdout);
     delete Proj;
   }
   printf("\n");
 }

◆ algoEndJob()

void SiStripGainFromCalibTree::algoEndJob ( )

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 1029 of file SiStripGainFromCalibTree.cc.

                                           {
   if (AlgoMode == "PCL" && !m_harvestingMode)
     return;  //nothing to do in that case
 
   if (AlgoMode == "CalibTree") {
     edm::LogInfo("SiStripGainFromCalibTree") << "Analyzing calibration tree" << std::endl;
     // Loop on calibTrees to fill the 2D histograms
     algoAnalyzeTheTree();
   } else if (m_harvestingMode) {
     NClusterStrip = (Charge_Vs_Index[Harvest])->getTH2S()->Integral(0, NStripAPVs + 1, 0, 99999);
     //NClusterPixel = (Charge_Vs_Index[Harvest])->getTH2S()->Integral(NStripAPVs+2, NStripAPVs+NPixelDets+2, 0, 99999 );
   }
 
   // Now that we have the full statistics we can extract the information of the 2D histograms
   algoComputeMPVandGain();
 
   // Result monitoring
   qualityMonitor();
 
   // Force the DB object writing,
   // thus setting the IOV as the first processed run (if timeFromEndRun is set to false)
   storeOnDbNow();
 
   if (AlgoMode != "PCL" or saveSummary) {
     edm::LogInfo("SiStripGainFromCalibTree") << "Saving summary into root file" << std::endl;
 
     //also save the 2D monitor elements to this file as TH2D tfs
     tfs = edm::Service<TFileService>().operator->();
 
     //save only the statistics for the calibrationTag
     int elepos = statCollectionFromMode(m_calibrationMode.c_str());
 
     if (Charge_Vs_Index[elepos] != nullptr)
       tfs->make<TH2S>(*(Charge_Vs_Index[elepos])->getTH2S());
     //if( Charge_Vs_Index_Absolute[elepos]!=0 )  tfs->make<TH2S> ( *(Charge_Vs_Index_Absolute[elepos])->getTH2S() );
     if (Charge_Vs_PathlengthTIB[elepos] != nullptr)
       tfs->make<TH2S>(*(Charge_Vs_PathlengthTIB[elepos])->getTH2S());
     if (Charge_Vs_PathlengthTOB[elepos] != nullptr)
       tfs->make<TH2S>(*(Charge_Vs_PathlengthTOB[elepos])->getTH2S());
     if (Charge_Vs_PathlengthTIDP[elepos] != nullptr)
       tfs->make<TH2S>(*(Charge_Vs_PathlengthTIDP[elepos])->getTH2S());
     if (Charge_Vs_PathlengthTIDM[elepos] != nullptr)
       tfs->make<TH2S>(*(Charge_Vs_PathlengthTIDM[elepos])->getTH2S());
     if (Charge_Vs_PathlengthTECP1[elepos] != nullptr)
       tfs->make<TH2S>(*(Charge_Vs_PathlengthTECP1[elepos])->getTH2S());
     if (Charge_Vs_PathlengthTECP2[elepos] != nullptr)
       tfs->make<TH2S>(*(Charge_Vs_PathlengthTECP2[elepos])->getTH2S());
     if (Charge_Vs_PathlengthTECM1[elepos] != nullptr)
       tfs->make<TH2S>(*(Charge_Vs_PathlengthTECM1[elepos])->getTH2S());
     if (Charge_Vs_PathlengthTECM2[elepos] != nullptr)
       tfs->make<TH2S>(*(Charge_Vs_PathlengthTECM2[elepos])->getTH2S());
 
     storeOnTree(tfs);
   }
 }

◆ algoEndRun()

void SiStripGainFromCalibTree::algoEndRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 869 of file SiStripGainFromCalibTree.cc.

                                                                                         {
   if (AlgoMode == "PCL" && !m_harvestingMode)
     return;  //nothing to do in that case
 
   if (AlgoMode == "PCL" and m_harvestingMode) {
     // Load the 2D histograms from the DQM objects
     // When running in AlCaHarvesting mode the histos are already booked and should be just retrieved from
     // DQMStore so that they can be used in the fit
 
     edm::LogInfo("SiStripGainFromCalibTree") << "Starting harvesting statistics" << std::endl;
 
     // check the required tag before adding histograms
     int elepos = statCollectionFromMode(m_calibrationMode.c_str());
     if (elepos != Harvest) {
       //collect statistics from DQM into the related monitored elements
       std::string stag = m_calibrationMode;
       if (!stag.empty() && stag[0] != '_')
         stag.insert(0, 1, '_');
 
       if (elepos == -1) {
         //implememt backward compatibility
         elepos = 0;
         stag = "";
       }
 
       std::string DQM_dir = m_DQMdir + ((m_splitDQMstat) ? m_calibrationMode : "");
 
       std::string cvi = DQM_dir + std::string("/Charge_Vs_Index") + stag;
       //std::string cviA     = DQM_dir + std::string("/Charge_Vs_Index_Absolute")  + stag;
       std::string cvpTIB = DQM_dir + std::string("/Charge_Vs_PathlengthTIB") + stag;
       std::string cvpTOB = DQM_dir + std::string("/Charge_Vs_PathlengthTOB") + stag;
       std::string cvpTIDP = DQM_dir + std::string("/Charge_Vs_PathlengthTIDP") + stag;
       std::string cvpTIDM = DQM_dir + std::string("/Charge_Vs_PathlengthTIDM") + stag;
       std::string cvpTECP1 = DQM_dir + std::string("/Charge_Vs_PathlengthTECP1") + stag;
       std::string cvpTECP2 = DQM_dir + std::string("/Charge_Vs_PathlengthTECP2") + stag;
       std::string cvpTECM1 = DQM_dir + std::string("/Charge_Vs_PathlengthTECM1") + stag;
       std::string cvpTECM2 = DQM_dir + std::string("/Charge_Vs_PathlengthTECM2") + stag;
 
       Charge_Vs_Index[elepos] = dbe->get(cvi);
       //Charge_Vs_Index_Absolute[elepos]  = dbe->get(cviA.c_str());
       Charge_Vs_PathlengthTIB[elepos] = dbe->get(cvpTIB);
       Charge_Vs_PathlengthTOB[elepos] = dbe->get(cvpTOB);
       Charge_Vs_PathlengthTIDP[elepos] = dbe->get(cvpTIDP);
       Charge_Vs_PathlengthTIDM[elepos] = dbe->get(cvpTIDM);
       Charge_Vs_PathlengthTECP1[elepos] = dbe->get(cvpTECP1);
       Charge_Vs_PathlengthTECP2[elepos] = dbe->get(cvpTECP2);
       Charge_Vs_PathlengthTECM1[elepos] = dbe->get(cvpTECM1);
       Charge_Vs_PathlengthTECM2[elepos] = dbe->get(cvpTECM2);
 
       if (Charge_Vs_Index[elepos] == nullptr) {
         edm::LogError("SiStripGainFromCalibTree")
             << "Harvesting: could not retrieve " << cvi.c_str() << ", statistics will not be summed!" << std::endl;
       } else {
         merge((Charge_Vs_Index[Harvest])->getTH2S(), (Charge_Vs_Index[elepos])->getTH2S());
         edm::LogInfo("SiStripGainFromCalibTree")
             << "Harvesting " << (Charge_Vs_Index[elepos])->getTH2S()->GetEntries() << " more clusters" << std::endl;
       }
 
       //if (Charge_Vs_Index_Absolute[elepos]==0) {
       //    edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << cviA.c_str()
       //                                              << ", statistics will not be summed!" << std::endl;
       //} else merge( (Charge_Vs_Index_Absolute[Harvest])->getTH2S(), (Charge_Vs_Index_Absolute[elepos])->getTH2S() );
 
       if (Charge_Vs_PathlengthTIB[elepos] == nullptr) {
         edm::LogError("SiStripGainFromCalibTree")
             << "Harvesting: could not retrieve " << cvpTIB.c_str() << ", statistics will not be summed!" << std::endl;
       } else
         (Charge_Vs_PathlengthTIB[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTIB[elepos])->getTH2S());
 
       if (Charge_Vs_PathlengthTOB[elepos] == nullptr) {
         edm::LogError("SiStripGainFromCalibTree")
             << "Harvesting: could not retrieve " << cvpTOB.c_str() << ", statistics will not be summed!" << std::endl;
       } else
         (Charge_Vs_PathlengthTOB[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTOB[elepos])->getTH2S());
 
       if (Charge_Vs_PathlengthTIDP[elepos] == nullptr) {
         edm::LogError("SiStripGainFromCalibTree")
             << "Harvesting: could not retrieve " << cvpTIDP.c_str() << ", statistics will not be summed!" << std::endl;
       } else
         (Charge_Vs_PathlengthTIDP[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTIDP[elepos])->getTH2S());
 
       if (Charge_Vs_PathlengthTIDM[elepos] == nullptr) {
         edm::LogError("SiStripGainFromCalibTree")
             << "Harvesting: could not retrieve " << cvpTIDM.c_str() << ", statistics will not be summed!" << std::endl;
       } else
         (Charge_Vs_PathlengthTIDM[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTIDM[elepos])->getTH2S());
 
       if (Charge_Vs_PathlengthTECP1[elepos] == nullptr) {
         edm::LogError("SiStripGainFromCalibTree")
             << "Harvesting: could not retrieve " << cvpTECP1.c_str() << ", statistics will not be summed!" << std::endl;
       } else
         (Charge_Vs_PathlengthTECP1[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTECP1[elepos])->getTH2S());
 
       if (Charge_Vs_PathlengthTECP2[elepos] == nullptr) {
         edm::LogError("SiStripGainFromCalibTree")
             << "Harvesting: could not retrieve " << cvpTECP2.c_str() << ", statistics will not be summed!" << std::endl;
       } else
         (Charge_Vs_PathlengthTECP2[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTECP2[elepos])->getTH2S());
 
       if (Charge_Vs_PathlengthTECM1[elepos] == nullptr) {
         edm::LogError("SiStripGainFromCalibTree")
             << "Harvesting: could not retrieve " << cvpTECM1.c_str() << ", statistics will not be summed!" << std::endl;
       } else
         (Charge_Vs_PathlengthTECM1[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTECM1[elepos])->getTH2S());
 
       if (Charge_Vs_PathlengthTECM2[elepos] == nullptr) {
         edm::LogError("SiStripGainFromCalibTree")
             << "Harvesting: could not retrieve " << cvpTECM2.c_str() << ", statistics will not be summed!" << std::endl;
       } else
         (Charge_Vs_PathlengthTECM2[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTECM2[elepos])->getTH2S());
 
       // Gather Charge monitoring histograms
       std::vector<std::pair<std::string, std::string>> tags =
           APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "");
       for (unsigned int i = 0; i < tags.size(); i++) {
         std::string tag = DQM_dir + "/" + (tags[i]).first + stag;
         Charge_1[elepos].push_back(APVGain::APVmon(0, 0, 0, 0, dbe->get(tag)));
         if ((Charge_1[elepos][i]).getMonitor() == nullptr) {
           edm::LogError("SiStripGainFromCalibTree")
               << "Harvesting: could not retrieve " << tag.c_str() << ", statistics will not be summed!" << std::endl;
         } else
           (Charge_1[Harvest][i]).getMonitor()->getTH1D()->Add((Charge_1[elepos][i]).getMonitor()->getTH1D());
       }
 
       tags = APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "woG2");
       for (unsigned int i = 0; i < tags.size(); i++) {
         std::string tag = DQM_dir + "/" + (tags[i]).first + stag;
         Charge_2[elepos].push_back(APVGain::APVmon(0, 0, 0, 0, dbe->get(tag)));
         if ((Charge_2[elepos][i]).getMonitor() == nullptr) {
           edm::LogError("SiStripGainFromCalibTree")
               << "Harvesting: could not retrieve " << tag.c_str() << ", statistics will not be summed!" << std::endl;
         } else
           (Charge_2[Harvest][i]).getMonitor()->getTH1D()->Add((Charge_2[elepos][i]).getMonitor()->getTH1D());
       }
 
       tags = APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "woG1");
       for (unsigned int i = 0; i < tags.size(); i++) {
         std::string tag = DQM_dir + "/" + (tags[i]).first + stag;
         Charge_3[elepos].push_back(APVGain::APVmon(0, 0, 0, 0, dbe->get(tag)));
         if ((Charge_3[elepos][i]).getMonitor() == nullptr) {
           edm::LogError("SiStripGainFromCalibTree")
               << "Harvesting: could not retrieve " << tag.c_str() << ", statistics will not be summed!" << std::endl;
         } else
           (Charge_3[Harvest][i]).getMonitor()->getTH1D()->Add((Charge_3[elepos][i]).getMonitor()->getTH1D());
       }
 
       tags = APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "woG1G2");
       for (unsigned int i = 0; i < tags.size(); i++) {
         std::string tag = DQM_dir + "/" + (tags[i]).first + stag;
         Charge_4[elepos].push_back(APVGain::APVmon(0, 0, 0, 0, dbe->get(tag)));
         if ((Charge_4[elepos][i]).getMonitor() == nullptr) {
           edm::LogError("SiStripGainFromCalibTree")
               << "Harvesting: could not retrieve " << tag.c_str() << ", statistics will not be summed!" << std::endl;
         } else
           (Charge_4[Harvest][i]).getMonitor()->getTH1D()->Add((Charge_4[elepos][i]).getMonitor()->getTH1D());
       }
     }
   }
 }

◆ bookDQMHistos()

void SiStripGainFromCalibTree::bookDQMHistos	(	const char *	dqm_dir,
		const char *	tag
	)

private

Definition at line 479 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob().

                                                                                  {
   edm::LogInfo("SiStripGainFromCalibTree")
       << "Setting " << dqm_dir << "in DQM and booking histograms for tag " << tag << std::endl;
 
   if (strcmp(booked_dir_.c_str(), dqm_dir) != 0) {
     booked_dir_ = dqm_dir;
     dbe->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 cviA     = std::string("Charge_Vs_Index_Absolute")  + 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 = (m_harvestingMode && AlgoMode == "PCL") ? Harvest : statCollectionFromMode(tag);
 
   // 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(NStripAPVs + 1);
   for (int a = 0; a <= 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.;
 
   Charge_Vs_Index[elepos] = dbe->book2S(cvi.c_str(), cvi.c_str(), NStripAPVs, &binXarray[0], 687, binYarray.data());
   //Charge_Vs_Index_Absolute[elepos]  = dbe->book2S(cviA.c_str()    , cviA.c_str()    , 88625, 0   , 88624,1000,0,4000);
   Charge_Vs_PathlengthTIB[elepos] = dbe->book2S(cvpTIB.c_str(), cvpTIB.c_str(), 20, 0.3, 1.3, 250, 0, 2000);
   Charge_Vs_PathlengthTOB[elepos] = dbe->book2S(cvpTOB.c_str(), cvpTOB.c_str(), 20, 0.3, 1.3, 250, 0, 2000);
   Charge_Vs_PathlengthTIDP[elepos] = dbe->book2S(cvpTIDP.c_str(), cvpTIDP.c_str(), 20, 0.3, 1.3, 250, 0, 2000);
   Charge_Vs_PathlengthTIDM[elepos] = dbe->book2S(cvpTIDM.c_str(), cvpTIDM.c_str(), 20, 0.3, 1.3, 250, 0, 2000);
   Charge_Vs_PathlengthTECP1[elepos] = dbe->book2S(cvpTECP1.c_str(), cvpTECP1.c_str(), 20, 0.3, 1.3, 250, 0, 2000);
   Charge_Vs_PathlengthTECP2[elepos] = dbe->book2S(cvpTECP2.c_str(), cvpTECP2.c_str(), 20, 0.3, 1.3, 250, 0, 2000);
   Charge_Vs_PathlengthTECM1[elepos] = dbe->book2S(cvpTECM1.c_str(), cvpTECM1.c_str(), 20, 0.3, 1.3, 250, 0, 2000);
   Charge_Vs_PathlengthTECM2[elepos] = dbe->book2S(cvpTECM2.c_str(), cvpTECM2.c_str(), 20, 0.3, 1.3, 250, 0, 2000);
 
   //Book Charge monitoring histograms
   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;
     MonitorElement* monitor = dbe->book1DD(htag.c_str(), (hnames[i]).second.c_str(), 100, 0., 1000.);
     int thick = APVGain::thickness((hnames[i]).first);
     int id = APVGain::subdetectorId((hnames[i]).first);
     int side = APVGain::subdetectorSide((hnames[i]).first);
     int plane = APVGain::subdetectorPlane((hnames[i]).first);
     Charge_1[elepos].push_back(APVGain::APVmon(thick, id, side, plane, monitor));
   }
 
   hnames = APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "woG2");
   for (unsigned int i = 0; i < hnames.size(); i++) {
     std::string htag = (hnames[i]).first + stag;
     MonitorElement* monitor = dbe->book1DD(htag.c_str(), (hnames[i]).second.c_str(), 100, 0., 1000.);
     int thick = APVGain::thickness((hnames[i]).first);
     int id = APVGain::subdetectorId((hnames[i]).first);
     int side = APVGain::subdetectorSide((hnames[i]).first);
     int plane = APVGain::subdetectorPlane((hnames[i]).first);
     Charge_2[elepos].push_back(APVGain::APVmon(thick, id, side, plane, monitor));
   }
 
   hnames = APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "woG1");
   for (unsigned int i = 0; i < hnames.size(); i++) {
     std::string htag = (hnames[i]).first + stag;
     MonitorElement* monitor = dbe->book1DD(htag.c_str(), (hnames[i]).second.c_str(), 100, 0., 1000.);
     int thick = APVGain::thickness((hnames[i]).first);
     int id = APVGain::subdetectorId((hnames[i]).first);
     int side = APVGain::subdetectorSide((hnames[i]).first);
     int plane = APVGain::subdetectorPlane((hnames[i]).first);
     Charge_3[elepos].push_back(APVGain::APVmon(thick, id, side, plane, monitor));
   }
 
   hnames = APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "woG1G2");
   for (unsigned int i = 0; i < hnames.size(); i++) {
     std::string htag = (hnames[i]).first + stag;
     MonitorElement* monitor = dbe->book1DD(htag.c_str(), (hnames[i]).second.c_str(), 100, 0., 1000.);
     int thick = APVGain::thickness((hnames[i]).first);
     int id = APVGain::subdetectorId((hnames[i]).first);
     int side = APVGain::subdetectorSide((hnames[i]).first);
     int plane = APVGain::subdetectorPlane((hnames[i]).first);
     Charge_4[elepos].push_back(APVGain::APVmon(thick, id, side, plane, monitor));
   }
 
   //Book validation histograms
   if (m_harvestingMode) {
     int MPVbin = 300;
     float MPVmin = 0.;
     float MPVmax = 600.;
 
     MPV_Vs_EtaTIB = dbe->book2DD("MPV_vs_EtaTIB", "MPV vs Eta TIB", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
     MPV_Vs_EtaTID = dbe->book2DD("MPV_vs_EtaTID", "MPV vs Eta TID", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
     MPV_Vs_EtaTOB = dbe->book2DD("MPV_vs_EtaTOB", "MPV vs Eta TOB", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
     MPV_Vs_EtaTEC = dbe->book2DD("MPV_vs_EtaTEC", "MPV vs Eta TEC", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
     MPV_Vs_EtaTECthin = dbe->book2DD("MPV_vs_EtaTEC1", "MPV vs Eta TEC-thin", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
     MPV_Vs_EtaTECthick = dbe->book2DD("MPV_vs_EtaTEC2", "MPV vs Eta TEC-thick", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
 
     MPV_Vs_PhiTIB = dbe->book2DD("MPV_vs_PhiTIB", "MPV vs Phi TIB", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
     MPV_Vs_PhiTID = dbe->book2DD("MPV_vs_PhiTID", "MPV vs Phi TID", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
     MPV_Vs_PhiTOB = dbe->book2DD("MPV_vs_PhiTOB", "MPV vs Phi TOB", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
     MPV_Vs_PhiTEC = dbe->book2DD("MPV_vs_PhiTEC", "MPV vs Phi TEC", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
     MPV_Vs_PhiTECthin = dbe->book2DD("MPV_vs_PhiTEC1", "MPV vs Phi TEC-thin", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
     MPV_Vs_PhiTECthick = dbe->book2DD("MPV_vs_PhiTEC2", "MPV vs Phi TEC-thick", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
 
     NoMPVfit = dbe->book2DD("NoMPVfit", "Modules with bad Landau Fit", 350, -350, 350, 240, 0, 120);
     NoMPVmasked = dbe->book2DD("NoMPVmasked", "Masked Modules", 350, -350, 350, 240, 0, 120);
 
     Gains = dbe->book1DD("Gains", "Gains", 300, 0, 2);
     MPVs = dbe->book1DD("MPVs", "MPVs", MPVbin, MPVmin, MPVmax);
     MPVs320 = dbe->book1DD("MPV_320", "MPV 320 thickness", MPVbin, MPVmin, MPVmax);
     MPVs500 = dbe->book1DD("MPV_500", "MPV 500 thickness", MPVbin, MPVmin, MPVmax);
     MPVsTIB = dbe->book1DD("MPV_TIB", "MPV TIB", MPVbin, MPVmin, MPVmax);
     MPVsTID = dbe->book1DD("MPV_TID", "MPV TID", MPVbin, MPVmin, MPVmax);
     MPVsTIDP = dbe->book1DD("MPV_TIDP", "MPV TIDP", MPVbin, MPVmin, MPVmax);
     MPVsTIDM = dbe->book1DD("MPV_TIDM", "MPV TIDM", MPVbin, MPVmin, MPVmax);
     MPVsTOB = dbe->book1DD("MPV_TOB", "MPV TOB", MPVbin, MPVmin, MPVmax);
     MPVsTEC = dbe->book1DD("MPV_TEC", "MPV TEC", MPVbin, MPVmin, MPVmax);
     MPVsTECP = dbe->book1DD("MPV_TECP", "MPV TECP", MPVbin, MPVmin, MPVmax);
     MPVsTECM = dbe->book1DD("MPV_TECM", "MPV TECM", MPVbin, MPVmin, MPVmax);
     MPVsTECthin = dbe->book1DD("MPV_TEC1", "MPV TEC1", MPVbin, MPVmin, MPVmax);
     MPVsTECthick = dbe->book1DD("MPV_TEC2", "MPV TEC2", MPVbin, MPVmin, MPVmax);
     MPVsTECP1 = dbe->book1DD("MPV_TECP1", "MPV TECP1", MPVbin, MPVmin, MPVmax);
     MPVsTECP2 = dbe->book1DD("MPV_TECP2", "MPV TECP2", MPVbin, MPVmin, MPVmax);
     MPVsTECM1 = dbe->book1DD("MPV_TECM1", "MPV TECM1", MPVbin, MPVmin, MPVmax);
     MPVsTECM2 = dbe->book1DD("MPV_TECM2", "MPV TECM2", MPVbin, MPVmin, MPVmax);
 
     MPVError = dbe->book1DD("MPVError", "MPV Error", 150, 0, 150);
     MPVErrorVsMPV = dbe->book2DD("MPVErrorVsMPV", "MPV Error vs MPV", 300, 0, 600, 150, 0, 150);
     MPVErrorVsEta = dbe->book2DD("MPVErrorVsEta", "MPV Error vs Eta", 50, -3.0, 3.0, 150, 0, 150);
     MPVErrorVsPhi = dbe->book2DD("MPVErrorVsPhi", "MPV Error vs Phi", 50, -3.4, 3.4, 150, 0, 150);
     MPVErrorVsN = dbe->book2DD("MPVErrorVsN", "MPV Error vs N", 500, 0, 1000, 150, 0, 150);
 
     DiffWRTPrevGainTIB = dbe->book1DD("DiffWRTPrevGainTIB", "Diff w.r.t. PrevGain TIB", 250, 0, 2);
     DiffWRTPrevGainTID = dbe->book1DD("DiffWRTPrevGainTID", "Diff w.r.t. PrevGain TID", 250, 0, 2);
     DiffWRTPrevGainTOB = dbe->book1DD("DiffWRTPrevGainTOB", "Diff w.r.t. PrevGain TOB", 250, 0, 2);
     DiffWRTPrevGainTEC = dbe->book1DD("DiffWRTPrevGainTEC", "Diff w.r.t. PrevGain TEC", 250, 0, 2);
 
     GainVsPrevGainTIB = dbe->book2DD("GainVsPrevGainTIB", "Gain vs PrevGain TIB", 100, 0, 2, 100, 0, 2);
     GainVsPrevGainTID = dbe->book2DD("GainVsPrevGainTID", "Gain vs PrevGain TID", 100, 0, 2, 100, 0, 2);
     GainVsPrevGainTOB = dbe->book2DD("GainVsPrevGainTOB", "Gain vs PrevGain TOB", 100, 0, 2, 100, 0, 2);
     GainVsPrevGainTEC = dbe->book2DD("GainVsPrevGainTEC", "Gain vs PrevGain TEC", 100, 0, 2, 100, 0, 2);
 
     std::vector<std::pair<std::string, std::string>> hnames =
         APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "newG2");
     for (unsigned int i = 0; i < hnames.size(); i++) {
       MonitorElement* monitor = dbe->book1DD((hnames[i]).first.c_str(), (hnames[i]).second.c_str(), 100, 0., 1000.);
       int thick = APVGain::thickness((hnames[i]).first);
       int id = APVGain::subdetectorId((hnames[i]).first);
       int side = APVGain::subdetectorSide((hnames[i]).first);
       int plane = APVGain::subdetectorPlane((hnames[i]).first);
       newCharge.push_back(APVGain::APVmon(thick, id, side, plane, monitor));
     }
   }
 }

◆ connect()

template<typename T >

edm::Handle<T> SiStripGainFromCalibTree::connect	(	const T *&	ptr,
		edm::EDGetTokenT< T >	token,
		const edm::Event &	evt
	)

inlineprivate

Definition at line 121 of file SiStripGainFromCalibTree.cc.

References edm::Event::getByToken(), patZpeak::handle, and unpackBuffers-CaloStage2::token.

Referenced by algoAnalyze(), and o2o_db_cfgmap.DbManagerDAQ::update_hashmap().

                                                                                            {
     edm::Handle<T> handle;
     evt.getByToken(token, handle);
     ptr = handle.product();
     return handle;  //return handle to keep alive pointer (safety first)
   }

◆ getNewObject()

std::unique_ptr< SiStripApvGain > SiStripGainFromCalibTree::getNewObject ( )

overrideprivatevirtual

Implements ConditionDBWriter< SiStripApvGain >.

Definition at line 1760 of file SiStripGainFromCalibTree.cc.

References a, sistrip::APV, APVsCollOrdered, m_harvestingMode, getGTfromDQMFile::obj, produceTagFilter(), FastTimerService_cff::range, and ConditionDBWriter< SiStripApvGain >::setDoStore().

                                                                      {
   auto obj = std::make_unique<SiStripApvGain>();
   if (!m_harvestingMode)
     return obj;
 
   if (!produceTagFilter()) {
     edm::LogWarning("SiStripGainFromCalibTree")
         << "getNewObject -> will not produce a paylaod because produceTagFilter returned false " << endl;
     setDoStore(false);
     return obj;
   }
 
   std::vector<float>* theSiStripVector = nullptr;
   unsigned int PreviousDetId = 0;
   for (unsigned int a = 0; a < APVsCollOrdered.size(); a++) {
     stAPVGain* APV = APVsCollOrdered[a];
     if (APV == nullptr) {
       printf("Bug\n");
       continue;
     }
     if (APV->SubDet <= 2)
       continue;
     if (APV->DetId != PreviousDetId) {
       if (theSiStripVector != nullptr) {
         SiStripApvGain::Range range(theSiStripVector->begin(), theSiStripVector->end());
         if (!obj->put(PreviousDetId, range))
           printf("Bug to put detId = %i\n", PreviousDetId);
       }
       theSiStripVector = new std::vector<float>;
       PreviousDetId = APV->DetId;
     }
     if (theSiStripVector != nullptr) {
       theSiStripVector->push_back(APV->Gain);
     }
   }
   if (theSiStripVector != nullptr) {
     SiStripApvGain::Range range(theSiStripVector->begin(), theSiStripVector->end());
     if (!obj->put(PreviousDetId, range))
       printf("Bug to put detId = %i\n", PreviousDetId);
   }
 
   if (theSiStripVector != nullptr)
     delete theSiStripVector;
 
   return obj;
 }

◆ getPeakOfLandau()

void SiStripGainFromCalibTree::getPeakOfLandau	(	TH1 *	InputHisto,
		double *	FitResults,
		double	LowRange = `50`,
		double	HighRange = `5400`
	)

private

Definition at line 1085 of file SiStripGainFromCalibTree.cc.

References MinNrEntries.

Referenced by algoComputeMPVandGain().

                                                                                                                      {
   FitResults[0] = -0.5;  //MPV
   FitResults[1] = 0;     //MPV error
   FitResults[2] = -0.5;  //Width
   FitResults[3] = 0;     //Width error
   FitResults[4] = -0.5;  //Fit Chi2/NDF
   FitResults[5] = 0;     //Normalization
 
   if (InputHisto->GetEntries() < MinNrEntries)
     return;
 
   // perform fit with standard landau
   TF1* MyLandau = new TF1("MyLandau", "landau", LowRange, HighRange);
   MyLandau->SetParameter(1, 300);
   InputHisto->Fit(MyLandau, "0QR WW");
 
   // MPV is parameter 1 (0=constant, 1=MPV, 2=Sigma)
   FitResults[0] = MyLandau->GetParameter(1);                      //MPV
   FitResults[1] = MyLandau->GetParError(1);                       //MPV error
   FitResults[2] = MyLandau->GetParameter(2);                      //Width
   FitResults[3] = MyLandau->GetParError(2);                       //Width error
   FitResults[4] = MyLandau->GetChisquare() / MyLandau->GetNDF();  //Fit Chi2/NDF
   FitResults[5] = MyLandau->GetParameter(0);
 
   delete MyLandau;
 }

◆ isBFieldConsistentWithMode()

bool SiStripGainFromCalibTree::isBFieldConsistentWithMode ( const edm::EventSetup & iSetup ) const

private

Definition at line 797 of file SiStripGainFromCalibTree.cc.

References edm::EventSetup::getData(), m_calibrationMode, MagFieldCurrentTh, submitPVValidationJobs::runInfo, and runInfoToken_.

Referenced by algoBeginJob(), and algoBeginRun().

                                                                                            {
   const auto& runInfo = iSetup.getData(runInfoToken_);
 
   double average_current = runInfo.m_avg_current;
   bool isOn = (average_current > MagFieldCurrentTh);
   bool is0T = (m_calibrationMode.substr(m_calibrationMode.length() - 2) == "0T");
 
   return ((isOn && !is0T) || (!isOn && is0T));
 }

◆ IsGoodLandauFit()

bool SiStripGainFromCalibTree::IsGoodLandauFit ( double * FitResults )

private

Definition at line 1112 of file SiStripGainFromCalibTree.cc.

Referenced by algoComputeMPVandGain().

                                                                  {
   if (FitResults[0] <= 0)
     return false;
   //   if(FitResults[1] > MaxMPVError   )return false;
   //   if(FitResults[4] > MaxChi2OverNDF)return false;
   return true;
 }

◆ MakeCalibrationMap()

void SiStripGainFromCalibTree::MakeCalibrationMap ( )

private

Definition at line 1817 of file SiStripGainFromCalibTree.cc.

References sistrip::APV, APVsColl, m_calibrationPath, RandomServiceHelper::t1, and useCalibration.

Referenced by algoBeginJob().

                                                   {
   if (!useCalibration)
     return;
 
   TChain* t1 = new TChain("SiStripCalib/APVGain");
   t1->Add(m_calibrationPath.c_str());
 
   unsigned int tree_DetId;
   unsigned char tree_APVId;
   double tree_Gain;
 
   t1->SetBranchAddress("DetId", &tree_DetId);
   t1->SetBranchAddress("APVId", &tree_APVId);
   t1->SetBranchAddress("Gain", &tree_Gain);
 
   for (unsigned int ientry = 0; ientry < t1->GetEntries(); ientry++) {
     t1->GetEntry(ientry);
     stAPVGain* APV = APVsColl[(tree_DetId << 4) | (unsigned int)tree_APVId];
     APV->CalibGain = tree_Gain;
   }
 
   delete t1;
 }

◆ merge()

void SiStripGainFromCalibTree::merge	(	TH2 *	A,
		TH2 *	B
	)

private

Definition at line 347 of file SiStripGainFromCalibTree.cc.

References A, B, and x.

Referenced by algoEndRun().

                                                    {
   if (A->GetNbinsX() == B->GetNbinsX()) {
     A->Add(B);
   } else {
     for (int x = 0; x <= B->GetNbinsX() + 1; x++) {
       for (int y = 0; y <= B->GetNbinsY() + 1; y++) {
         A->SetBinContent(x, y, A->GetBinContent(x, y) + B->GetBinContent(x, y));
       }
     }
   }
 }

◆ processEvent()

void SiStripGainFromCalibTree::processEvent ( )

private

Definition at line 1120 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

                                             {
   edm::LogInfo("SiStripGainFromCalibTree") << "Processing run " << runnumber << " and event " << eventnumber << " for "
                                            << m_calibrationMode << " calibration." << std::endl;
 
   if (runnumber < SRun)
     SRun = runnumber;
   if (runnumber > ERun)
     ERun = runnumber;
 
   NEvent++;
   NTrack += (*trackp).size();
 
   int elepos = statCollectionFromMode(m_calibrationMode.c_str());
 
   unsigned int FirstAmplitude = 0;
   for (unsigned int i = 0; i < (*chargeoverpath).size(); i++) {
     FirstAmplitude += (*nstrips)[i];
     int TI = (*trackindex)[i];
 
     //printf("%i - %i - %i %i %i\n", (int)(*rawid)[i], (int)(*firststrip)[i]/128, (int)(*farfromedge)[i], (int)(*overlapping)[i], (int)(*saturation )[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;
 
     stAPVGain* APV =
         APVsColl[((*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;
 
     //printf("detId=%7i run=%7i event=%9i charge=%5i cs=%3i\n",(*rawid)[i],runnumber,eventnumber,(*charge)[i],(*nstrips)[i]);
 
     //double trans = atan2((*localdiry)[i],(*localdirx)[i])*(180/3.14159265);
     //double alpha = acos ((*localdirx)[i] / sqrt( pow((*localdirx)[i],2) +  pow((*localdirz)[i],2) ) ) * (180/3.14159265);
     //double beta  = acos ((*localdiry)[i] / sqrt( pow((*localdirx)[i],2) +  pow((*localdirz)[i],2) ) ) * (180/3.14159265);
 
     //printf("NStrip = %i : Charge = %i --> Path = %f  --> ChargeOverPath=%f\n",(*nstrips)[i],(*charge)[i],(*path)[i],(*chargeoverpath)[i]);
     //printf("Amplitudes: ");
     //for(unsigned int a=0;a<(*nstrips)[i];a++){printf("%i ",(*amplitude)[FirstAmplitude+a]);}
     //printf("\n");
 
     if (APV->SubDet > 2) {
       NClusterStrip++;
     } else {
       NClusterPixel++;
     }
 
     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
     }
 
     //printf("ChargeDifference = %i Vs %i with Gain = %f\n",(*charge)[i],Charge,APV->CalibGain);
 
     double ClusterChargeOverPath = ((double)Charge) / (*path)[i];
     if (APV->SubDet > 2) {
       if (Validation) {
         ClusterChargeOverPath /= (*gainused)[i];
       }
       if (OldGainRemoving) {
         ClusterChargeOverPath *= (*gainused)[i];
       }
     }
 
     // keep pixel cluster charge processing until here
     if (APV->SubDet <= 2)
       continue;
 
     (Charge_Vs_Index[elepos])->Fill(APV->Index, ClusterChargeOverPath);
 
     // Compute the charge for monitoring and fill the relative 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
     }
     std::vector<APVGain::APVmon>& v1 = Charge_1[elepos];
     std::vector<MonitorElement*> cmon1 = APVGain::FetchMonitor(v1, (*rawid)[i], tTopo_);
     for (unsigned int m = 0; m < cmon1.size(); m++)
       cmon1[m]->Fill(((double)mCharge1) / (*path)[i]);
 
     std::vector<APVGain::APVmon>& v2 = Charge_2[elepos];
     std::vector<MonitorElement*> cmon2 = APVGain::FetchMonitor(v2, (*rawid)[i], tTopo_);
     for (unsigned int m = 0; m < cmon2.size(); m++)
       cmon2[m]->Fill(((double)mCharge2) / (*path)[i]);
 
     std::vector<APVGain::APVmon>& v3 = Charge_3[elepos];
     std::vector<MonitorElement*> cmon3 = APVGain::FetchMonitor(v3, (*rawid)[i], tTopo_);
     for (unsigned int m = 0; m < cmon3.size(); m++)
       cmon3[m]->Fill(((double)mCharge3) / (*path)[i]);
 
     std::vector<APVGain::APVmon>& v4 = Charge_4[elepos];
     std::vector<MonitorElement*> cmon4 = APVGain::FetchMonitor(v4, (*rawid)[i], tTopo_);
     for (unsigned int m = 0; m < cmon4.size(); m++)
       cmon4[m]->Fill(((double)mCharge4) / (*path)[i]);
 
     // Fill Charge Vs pathLenght histograms
     if (APV->SubDet == StripSubdetector::TIB) {
       (Charge_Vs_PathlengthTIB[elepos])->Fill((*path)[i], Charge);  // TIB
 
     } else if (APV->SubDet == StripSubdetector::TOB) {
       (Charge_Vs_PathlengthTOB[elepos])->Fill((*path)[i], Charge);  // TOB
 
     } else if (APV->SubDet == StripSubdetector::TID) {
       if (APV->Eta < 0) {
         (Charge_Vs_PathlengthTIDM[elepos])->Fill((*path)[i], Charge);
       }  // TID minus
       else if (APV->Eta > 0) {
         (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) {
           (Charge_Vs_PathlengthTECM1[elepos])->Fill((*path)[i], Charge);
         }  // TEC minus, type 1
         else if (APV->Thickness > 0.04) {
           (Charge_Vs_PathlengthTECM2[elepos])->Fill((*path)[i], Charge);
         }  // TEC minus, type 2
       } else if (APV->Eta > 0) {
         if (APV->Thickness < 0.04) {
           (Charge_Vs_PathlengthTECP1[elepos])->Fill((*path)[i], Charge);
         }  // TEC plus, type 1
         else if (APV->Thickness > 0.04) {
           (Charge_Vs_PathlengthTECP2[elepos])->Fill((*path)[i], Charge);
         }  // TEC plus, type 2
       }
     }
 
   }  // END OF ON-CLUSTER LOOP
 }  //END OF processEvent()

◆ produceTagFilter()

bool SiStripGainFromCalibTree::produceTagFilter ( )

private

Definition at line 1733 of file SiStripGainFromCalibTree.cc.

References AlgoMode, BAD, Charge_Vs_Index, GOOD, Harvest, funct::integral(), m_calibrationMode, NStripAPVs, statCollectionFromMode(), tagCondition_GoodFrac, and tagCondition_NClusters.

Referenced by getNewObject().

                                                 {
   // The goal of this function is to check wether or not there is enough statistics
   // to produce a meaningful tag for the DB
   int elepos = (AlgoMode == "PCL") ? Harvest : statCollectionFromMode(m_calibrationMode.c_str());
   if (Charge_Vs_Index[elepos] == nullptr) {
     edm::LogError("SiStripGainFromCalibTree")
         << "produceTagFilter -> Return false: could not retrieve the " << m_calibrationMode.c_str() << " statistics.\n"
         << "Please check if input contains " << m_calibrationMode.c_str() << " data." << std::endl;
     return false;
   }
 
   float integral = (Charge_Vs_Index[elepos])->getTH2S()->Integral();
   if ((Charge_Vs_Index[elepos])->getTH2S()->Integral(0, NStripAPVs + 1, 0, 99999) < tagCondition_NClusters) {
     edm::LogWarning("SiStripGainFromCalibTree")
         << "calibrationMode  -> " << m_calibrationMode << "\n"
         << "produceTagFilter -> Return false: Statistics is too low : " << integral << endl;
     return false;
   }
   if ((1.0 * GOOD) / (GOOD + BAD) < tagCondition_GoodFrac) {
     edm::LogWarning("SiStripGainFromCalibTree")
         << "calibrationMode  -> " << m_calibrationMode << "\n"
         << "produceTagFilter ->  Return false: ratio of GOOD/TOTAL is too low: " << (1.0 * GOOD) / (GOOD + BAD) << endl;
     return false;
   }
   return true;
 }

◆ qualityMonitor()

void SiStripGainFromCalibTree::qualityMonitor ( )

private

Definition at line 1472 of file SiStripGainFromCalibTree.cc.

References a, AlgoMode, sistrip::APV, APVsCollOrdered, newFWLiteAna::bin, Charge_Vs_Index, DiffWRTPrevGainTEC, DiffWRTPrevGainTIB, DiffWRTPrevGainTID, DiffWRTPrevGainTOB, MillePedeFileConverter_cfg::e, APVGain::FetchMonitor(), dqm::impl::MonitorElement::Fill(), IntegrityClient_cfi::Gain, Gains, GainVsPrevGainTEC, GainVsPrevGainTIB, GainVsPrevGainTID, GainVsPrevGainTOB, h, Harvest, m_calibrationMode, MPV_Vs_EtaTEC, MPV_Vs_EtaTECthick, MPV_Vs_EtaTECthin, MPV_Vs_EtaTIB, MPV_Vs_EtaTID, MPV_Vs_EtaTOB, MPV_Vs_PhiTEC, MPV_Vs_PhiTECthick, MPV_Vs_PhiTECthin, MPV_Vs_PhiTIB, MPV_Vs_PhiTID, MPV_Vs_PhiTOB, MPVError, MPVErrorVsEta, MPVErrorVsMPV, MPVErrorVsN, MPVErrorVsPhi, MPVs, MPVs320, MPVs500, MPVsTEC, MPVsTECM, MPVsTECM1, MPVsTECM2, MPVsTECP, MPVsTECP1, MPVsTECP2, MPVsTECthick, MPVsTECthin, MPVsTIB, MPVsTID, MPVsTIDM, MPVsTIDP, MPVsTOB, newCharge, NoMPVfit, NoMPVmasked, VtxSmearedParameters_cfi::Phi, dttmaxenums::R, statCollectionFromMode(), ntupleEnum::SubDet, tTopo_, and z.

Referenced by algoEndJob().

                                               {
   int elepos = (AlgoMode == "PCL") ? Harvest : statCollectionFromMode(m_calibrationMode.c_str());
 
   for (unsigned int a = 0; a < APVsCollOrdered.size(); a++) {
     stAPVGain* APV = APVsCollOrdered[a];
     if (APV == nullptr)
       continue;
 
     unsigned int Index = APV->Index;
     unsigned int SubDet = APV->SubDet;
     unsigned int DetId = APV->DetId;
     float z = APV->z;
     float Eta = APV->Eta;
     float R = APV->R;
     float Phi = APV->Phi;
     float Thickness = APV->Thickness;
     double FitMPV = APV->FitMPV;
     double FitMPVErr = APV->FitMPVErr;
     double Gain = APV->Gain;
     double NEntries = APV->NEntries;
     double PreviousGain = APV->PreviousGain;
 
     if (SubDet < 3)
       continue;  // avoid to loop over Pixel det id
 
     if (Gain != 1.) {
       std::vector<MonitorElement*> charge_histos = APVGain::FetchMonitor(newCharge, DetId, tTopo_);
       TH2S* chvsidx = (Charge_Vs_Index[elepos])->getTH2S();
       int bin = chvsidx->GetXaxis()->FindBin(Index);
       TH1D* Proj = chvsidx->ProjectionY("proj", bin, bin);
       for (int binId = 0; binId < Proj->GetXaxis()->GetNbins(); binId++) {
         double new_charge = Proj->GetXaxis()->GetBinCenter(binId) / Gain;
         if (Proj->GetBinContent(binId) != 0.) {
           for (unsigned int h = 0; h < charge_histos.size(); h++) {
             TH1D* chisto = (charge_histos[h])->getTH1D();
             for (int e = 0; e < Proj->GetBinContent(binId); e++)
               chisto->Fill(new_charge);
           }
         }
       }
     }
 
     if (FitMPV <= 0.) {  // No fit of MPV
       if (APV->isMasked)
         NoMPVmasked->Fill(z, R);
       else
         NoMPVfit->Fill(z, R);
 
     } else {  // Fit of MPV
       if (FitMPV > 0.)
         Gains->Fill(Gain);
 
       MPVs->Fill(FitMPV);
       if (Thickness < 0.04)
         MPVs320->Fill(FitMPV);
       if (Thickness > 0.04)
         MPVs500->Fill(FitMPV);
 
       MPVError->Fill(FitMPVErr);
       MPVErrorVsMPV->Fill(FitMPV, FitMPVErr);
       MPVErrorVsEta->Fill(Eta, FitMPVErr);
       MPVErrorVsPhi->Fill(Phi, FitMPVErr);
       MPVErrorVsN->Fill(NEntries, FitMPVErr);
 
       if (SubDet == 3) {
         MPV_Vs_EtaTIB->Fill(Eta, FitMPV);
         MPV_Vs_PhiTIB->Fill(Phi, FitMPV);
         MPVsTIB->Fill(FitMPV);
 
       } else if (SubDet == 4) {
         MPV_Vs_EtaTID->Fill(Eta, FitMPV);
         MPV_Vs_PhiTID->Fill(Phi, FitMPV);
         MPVsTID->Fill(FitMPV);
         if (Eta < 0.)
           MPVsTIDM->Fill(FitMPV);
         if (Eta > 0.)
           MPVsTIDP->Fill(FitMPV);
 
       } else if (SubDet == 5) {
         MPV_Vs_EtaTOB->Fill(Eta, FitMPV);
         MPV_Vs_PhiTOB->Fill(Phi, FitMPV);
         MPVsTOB->Fill(FitMPV);
 
       } else if (SubDet == 6) {
         MPV_Vs_EtaTEC->Fill(Eta, FitMPV);
         MPV_Vs_PhiTEC->Fill(Phi, FitMPV);
         MPVsTEC->Fill(FitMPV);
         if (Eta < 0.)
           MPVsTECM->Fill(FitMPV);
         if (Eta > 0.)
           MPVsTECP->Fill(FitMPV);
         if (Thickness < 0.04) {
           MPV_Vs_EtaTECthin->Fill(Eta, FitMPV);
           MPV_Vs_PhiTECthin->Fill(Phi, FitMPV);
           MPVsTECthin->Fill(FitMPV);
           if (Eta > 0.)
             MPVsTECP1->Fill(FitMPV);
           if (Eta < 0.)
             MPVsTECM1->Fill(FitMPV);
         }
         if (Thickness > 0.04) {
           MPV_Vs_EtaTECthick->Fill(Eta, FitMPV);
           MPV_Vs_PhiTECthick->Fill(Phi, FitMPV);
           MPVsTECthick->Fill(FitMPV);
           if (Eta > 0.)
             MPVsTECP2->Fill(FitMPV);
           if (Eta < 0.)
             MPVsTECM2->Fill(FitMPV);
         }
       }
     }
 
     if (SubDet == 3 && PreviousGain != 0.)
       DiffWRTPrevGainTIB->Fill(Gain / PreviousGain);
     else if (SubDet == 4 && PreviousGain != 0.)
       DiffWRTPrevGainTID->Fill(Gain / PreviousGain);
     else if (SubDet == 5 && PreviousGain != 0.)
       DiffWRTPrevGainTOB->Fill(Gain / PreviousGain);
     else if (SubDet == 6 && PreviousGain != 0.)
       DiffWRTPrevGainTEC->Fill(Gain / PreviousGain);
 
     if (SubDet == 3)
       GainVsPrevGainTIB->Fill(PreviousGain, Gain);
     else if (SubDet == 4)
       GainVsPrevGainTID->Fill(PreviousGain, Gain);
     else if (SubDet == 5)
       GainVsPrevGainTOB->Fill(PreviousGain, Gain);
     else if (SubDet == 6)
       GainVsPrevGainTEC->Fill(PreviousGain, Gain);
   }
 }

◆ statCollectionFromMode()

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

inlineprivate

Definition at line 333 of file SiStripGainFromCalibTree.cc.

References relativeConstraints::empty, None, AlCaHLTBitMon_QueryRunRegistry::string, and makeGlobalPositionRcd_cfg::tag.

Referenced by algoBeginJob(), algoComputeMPVandGain(), algoEndJob(), algoEndRun(), bookDQMHistos(), processEvent(), produceTagFilter(), qualityMonitor(), and storeOnTree().

                                                                                  {
   std::vector<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;
 }

◆ storeOnTree()

void SiStripGainFromCalibTree::storeOnTree ( TFileService * tfs )

private

Definition at line 1604 of file SiStripGainFromCalibTree.cc.

References a, sistrip::APV, APVsCollOrdered, BAD, Charge_Vs_Index, Gains, GOOD, m_calibrationMode, TFileService::make(), MASKED, NClusterStrip, NEvent, NStripAPVs, NTrack, OutputGains, statCollectionFromMode(), mps_setup::stdout, and tfs.

Referenced by algoEndJob().

                                                             {
   unsigned int tree_Index;
   unsigned int tree_Bin;
   unsigned int tree_DetId;
   unsigned char tree_APVId;
   unsigned char tree_SubDet;
   float tree_x;
   float tree_y;
   float tree_z;
   float tree_Eta;
   float tree_R;
   float tree_Phi;
   float tree_Thickness;
   float tree_FitMPV;
   float tree_FitMPVErr;
   float tree_FitWidth;
   float tree_FitWidthErr;
   float tree_FitChi2NDF;
   float tree_FitNorm;
   double tree_Gain;
   double tree_PrevGain;
   double tree_PrevGainTick;
   double tree_NEntries;
   bool tree_isMasked;
 
   TTree* MyTree;
   MyTree = tfs->make<TTree>("APVGain", "APVGain");
   MyTree->Branch("Index", &tree_Index, "Index/i");
   MyTree->Branch("Bin", &tree_Bin, "Bin/i");
   MyTree->Branch("DetId", &tree_DetId, "DetId/i");
   MyTree->Branch("APVId", &tree_APVId, "APVId/b");
   MyTree->Branch("SubDet", &tree_SubDet, "SubDet/b");
   MyTree->Branch("x", &tree_x, "x/F");
   MyTree->Branch("y", &tree_y, "y/F");
   MyTree->Branch("z", &tree_z, "z/F");
   MyTree->Branch("Eta", &tree_Eta, "Eta/F");
   MyTree->Branch("R", &tree_R, "R/F");
   MyTree->Branch("Phi", &tree_Phi, "Phi/F");
   MyTree->Branch("Thickness", &tree_Thickness, "Thickness/F");
   MyTree->Branch("FitMPV", &tree_FitMPV, "FitMPV/F");
   MyTree->Branch("FitMPVErr", &tree_FitMPVErr, "FitMPVErr/F");
   MyTree->Branch("FitWidth", &tree_FitWidth, "FitWidth/F");
   MyTree->Branch("FitWidthErr", &tree_FitWidthErr, "FitWidthErr/F");
   MyTree->Branch("FitChi2NDF", &tree_FitChi2NDF, "FitChi2NDF/F");
   MyTree->Branch("FitNorm", &tree_FitNorm, "FitNorm/F");
   MyTree->Branch("Gain", &tree_Gain, "Gain/D");
   MyTree->Branch("PrevGain", &tree_PrevGain, "PrevGain/D");
   MyTree->Branch("PrevGainTick", &tree_PrevGainTick, "PrevGainTick/D");
   MyTree->Branch("NEntries", &tree_NEntries, "NEntries/D");
   MyTree->Branch("isMasked", &tree_isMasked, "isMasked/O");
 
   FILE* Gains = stdout;
   fprintf(Gains, "NEvents   = %i\n", NEvent);
   fprintf(Gains, "NTracks   = %i\n", NTrack);
   //fprintf(Gains,"NClustersPixel = %i\n",NClusterPixel);
   fprintf(Gains, "NClustersStrip = %i\n", NClusterStrip);
   //fprintf(Gains,"Number of Pixel Dets = %lu\n",static_cast<unsigned long>(NPixelDets));
   fprintf(Gains, "Number of Strip APVs = %lu\n", static_cast<unsigned long>(NStripAPVs));
   fprintf(Gains,
           "GoodFits = %i BadFits = %i ratio = %f%%   (MASKED=%i)\n",
           GOOD,
           BAD,
           (100.0 * GOOD) / (GOOD + BAD),
           MASKED);
 
   Gains = fopen(OutputGains.c_str(), "w");
   fprintf(Gains, "NEvents   = %i\n", NEvent);
   fprintf(Gains, "NTracks   = %i\n", NTrack);
   //fprintf(Gains,"NClustersPixel = %i\n",NClusterPixel);
   fprintf(Gains, "NClustersStrip = %i\n", NClusterStrip);
   fprintf(Gains, "Number of Strip APVs = %lu\n", static_cast<unsigned long>(NStripAPVs));
   //fprintf(Gains,"Number of Pixel Dets = %lu\n",static_cast<unsigned long>(NPixelDets));
   fprintf(Gains,
           "GoodFits = %i BadFits = %i ratio = %f%%   (MASKED=%i)\n",
           GOOD,
           BAD,
           (100.0 * GOOD) / (GOOD + BAD),
           MASKED);
 
   int elepos = statCollectionFromMode(m_calibrationMode.c_str());
 
   for (unsigned int a = 0; a < APVsCollOrdered.size(); a++) {
     stAPVGain* APV = APVsCollOrdered[a];
     if (APV == nullptr)
       continue;
     //     printf(      "%i | %i | PreviousGain = %7.5f NewGain = %7.5f (#clusters=%8.0f)\n", APV->DetId,APV->APVId,APV->PreviousGain,APV->Gain, APV->NEntries);
     fprintf(Gains,
             "%i | %i | PreviousGain = %7.5f(tick) x %7.5f(particle) NewGain (particle) = %7.5f (#clusters=%8.0f)\n",
             APV->DetId,
             APV->APVId,
             APV->PreviousGainTick,
             APV->PreviousGain,
             APV->Gain,
             APV->NEntries);
 
     tree_Index = APV->Index;
     tree_Bin = (Charge_Vs_Index[elepos])->getTH2S()->GetXaxis()->FindBin(APV->Index);
     tree_DetId = APV->DetId;
     tree_APVId = APV->APVId;
     tree_SubDet = APV->SubDet;
     tree_x = APV->x;
     tree_y = APV->y;
     tree_z = APV->z;
     tree_Eta = APV->Eta;
     tree_R = APV->R;
     tree_Phi = APV->Phi;
     tree_Thickness = APV->Thickness;
     tree_FitMPV = APV->FitMPV;
     tree_FitMPVErr = APV->FitMPVErr;
     tree_FitWidth = APV->FitWidth;
     tree_FitWidthErr = APV->FitWidthErr;
     tree_FitChi2NDF = APV->FitChi2;
     tree_FitNorm = APV->FitNorm;
     tree_Gain = APV->Gain;
     tree_PrevGain = APV->PreviousGain;
     tree_PrevGainTick = APV->PreviousGainTick;
     tree_NEntries = APV->NEntries;
     tree_isMasked = APV->isMasked;
 
     if (tree_DetId == 402673324) {
       printf("%i | %i : %f --> %f  (%f)\n", tree_DetId, tree_APVId, tree_PrevGain, tree_Gain, tree_NEntries);
     }
 
     MyTree->Fill();
   }
   if (Gains)
     fclose(Gains);
 }

◆ swapBFieldMode()

void SiStripGainFromCalibTree::swapBFieldMode ( void )

private

Definition at line 807 of file SiStripGainFromCalibTree.cc.

References m_calibrationMode.

Referenced by algoBeginJob(), and algoBeginRun().

                                               {
   if (m_calibrationMode.substr(m_calibrationMode.length() - 2) == "0T") {
     m_calibrationMode.erase(m_calibrationMode.length() - 2, 2);
   } else {
     m_calibrationMode.append("0T");
   }
 }

Member Data Documentation

◆ AlgoMode

std::string SiStripGainFromCalibTree::AlgoMode

private

Definition at line 163 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoBeginJob(), algoBeginRun(), algoComputeMPVandGain(), algoEndJob(), algoEndRun(), bookDQMHistos(), produceTagFilter(), qualityMonitor(), and SiStripGainFromCalibTree().

◆ AllowSaturation

bool SiStripGainFromCalibTree::AllowSaturation

private

Definition at line 144 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ amplitude

const std::vector<unsigned char>* SiStripGainFromCalibTree::amplitude = nullptr

private

Definition at line 307 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ amplitude_token_

edm::EDGetTokenT<std::vector<unsigned char> > SiStripGainFromCalibTree::amplitude_token_

private

Definition at line 308 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ APVsColl

std::unordered_map<unsigned int, stAPVGain*> SiStripGainFromCalibTree::APVsColl

private

Definition at line 323 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), algoComputeMPVandGain(), MakeCalibrationMap(), processEvent(), and ~SiStripGainFromCalibTree().

◆ APVsCollOrdered

std::vector<stAPVGain*> SiStripGainFromCalibTree::APVsCollOrdered

private

Definition at line 322 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), algoBeginRun(), getNewObject(), qualityMonitor(), storeOnTree(), and ~SiStripGainFromCalibTree().

◆ BAD

unsigned int SiStripGainFromCalibTree::BAD

private

Definition at line 253 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), algoComputeMPVandGain(), produceTagFilter(), and storeOnTree().

◆ booked_dir_

std::string SiStripGainFromCalibTree::booked_dir_

private

Definition at line 170 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos().

◆ CalibPrefix_

string SiStripGainFromCalibTree::CalibPrefix_

private

Definition at line 318 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

◆ CalibrationLevel

int SiStripGainFromCalibTree::CalibrationLevel

private

Definition at line 148 of file SiStripGainFromCalibTree.cc.

Referenced by algoComputeMPVandGain(), and SiStripGainFromCalibTree().

◆ CalibSuffix_

string SiStripGainFromCalibTree::CalibSuffix_

private

Definition at line 319 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

◆ charge

const std::vector<unsigned int>* SiStripGainFromCalibTree::charge = nullptr

private

Definition at line 301 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ Charge_1

std::array<std::vector<APVGain::APVmon>, 7> SiStripGainFromCalibTree::Charge_1

private

Charge per cm per layer / wheel

Definition at line 173 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndRun(), bookDQMHistos(), and processEvent().

◆ Charge_2

std::array<std::vector<APVGain::APVmon>, 7> SiStripGainFromCalibTree::Charge_2

private

Charge per cm per layer / wheel without G2

Definition at line 174 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndRun(), bookDQMHistos(), and processEvent().

◆ Charge_3

std::array<std::vector<APVGain::APVmon>, 7> SiStripGainFromCalibTree::Charge_3

private

Charge per cm per layer / wheel without G1

Definition at line 175 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndRun(), bookDQMHistos(), and processEvent().

◆ Charge_4

std::array<std::vector<APVGain::APVmon>, 7> SiStripGainFromCalibTree::Charge_4

private

Charge per cm per layer / wheel without G1 and G1

Definition at line 176 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndRun(), bookDQMHistos(), and processEvent().

◆ charge_token_

edm::EDGetTokenT<std::vector<unsigned int> > SiStripGainFromCalibTree::charge_token_

private

Definition at line 302 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ Charge_Vs_Index

std::vector<MonitorElement*> SiStripGainFromCalibTree::Charge_Vs_Index

private

Charge per cm for each detector id

Definition at line 172 of file SiStripGainFromCalibTree.cc.

Referenced by algoComputeMPVandGain(), algoEndJob(), algoEndRun(), bookDQMHistos(), processEvent(), produceTagFilter(), qualityMonitor(), SiStripGainFromCalibTree(), and storeOnTree().

◆ Charge_Vs_PathlengthTECM1

std::vector<MonitorElement*> SiStripGainFromCalibTree::Charge_Vs_PathlengthTECM1

private

Charge vs pathlength in TECP thin

Definition at line 184 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), algoEndRun(), bookDQMHistos(), processEvent(), and SiStripGainFromCalibTree().

◆ Charge_Vs_PathlengthTECM2

std::vector<MonitorElement*> SiStripGainFromCalibTree::Charge_Vs_PathlengthTECM2

private

Charge vs pathlength in TECP thick

Definition at line 185 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), algoEndRun(), bookDQMHistos(), processEvent(), and SiStripGainFromCalibTree().

◆ Charge_Vs_PathlengthTECP1

std::vector<MonitorElement*> SiStripGainFromCalibTree::Charge_Vs_PathlengthTECP1

private

Charge vs pathlength in TECP thin

Definition at line 182 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), algoEndRun(), bookDQMHistos(), processEvent(), and SiStripGainFromCalibTree().

◆ Charge_Vs_PathlengthTECP2

std::vector<MonitorElement*> SiStripGainFromCalibTree::Charge_Vs_PathlengthTECP2

private

Charge vs pathlength in TECP thick

Definition at line 183 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), algoEndRun(), bookDQMHistos(), processEvent(), and SiStripGainFromCalibTree().

◆ Charge_Vs_PathlengthTIB

std::vector<MonitorElement*> SiStripGainFromCalibTree::Charge_Vs_PathlengthTIB

private

Charge vs pathlength in TIB

Definition at line 178 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), algoEndRun(), bookDQMHistos(), processEvent(), and SiStripGainFromCalibTree().

◆ Charge_Vs_PathlengthTIDM

std::vector<MonitorElement*> SiStripGainFromCalibTree::Charge_Vs_PathlengthTIDM

private

Charge vs pathlength in TIDM

Definition at line 181 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), algoEndRun(), bookDQMHistos(), processEvent(), and SiStripGainFromCalibTree().

◆ Charge_Vs_PathlengthTIDP

std::vector<MonitorElement*> SiStripGainFromCalibTree::Charge_Vs_PathlengthTIDP

private

Charge vs pathlength in TIDP

Definition at line 180 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), algoEndRun(), bookDQMHistos(), processEvent(), and SiStripGainFromCalibTree().

◆ Charge_Vs_PathlengthTOB

std::vector<MonitorElement*> SiStripGainFromCalibTree::Charge_Vs_PathlengthTOB

private

Charge vs pathlength in TOB

Definition at line 179 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), algoEndRun(), bookDQMHistos(), processEvent(), and SiStripGainFromCalibTree().

◆ chargeoverpath

const std::vector<double>* SiStripGainFromCalibTree::chargeoverpath = nullptr

private

Definition at line 305 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ chargeoverpath_token_

edm::EDGetTokenT<std::vector<double> > SiStripGainFromCalibTree::chargeoverpath_token_

private

Definition at line 306 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ dbe

DQMStore* SiStripGainFromCalibTree::dbe

private

Definition at line 131 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndRun(), bookDQMHistos(), and SiStripGainFromCalibTree().

◆ DiffWRTPrevGainTEC

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTEC

private

ratio gain / PreviousGain for TEC disks

Definition at line 235 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ DiffWRTPrevGainTIB

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTIB

private

ratio Gain / PreviousGain for TIB layers

Definition at line 232 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ DiffWRTPrevGainTID

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTID

private

ratio Gain / PreviousGain for TID disks

Definition at line 233 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ DiffWRTPrevGainTOB

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTOB

private

ratio Gain / PreviousGain for TOB layers

Definition at line 234 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ doChargeMonitorPerPlane

bool SiStripGainFromCalibTree::doChargeMonitorPerPlane

private

Charge monitor per detector plane

Definition at line 154 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndRun(), bookDQMHistos(), and SiStripGainFromCalibTree().

◆ dqm_tag_

std::vector<string> SiStripGainFromCalibTree::dqm_tag_

private

Definition at line 169 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and SiStripGainFromCalibTree().

◆ ERun

unsigned int SiStripGainFromCalibTree::ERun

private

Definition at line 251 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and processEvent().

◆ eventnumber

unsigned int SiStripGainFromCalibTree::eventnumber = 0

private

Definition at line 259 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ EventPrefix_

string SiStripGainFromCalibTree::EventPrefix_

private

Definition at line 314 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

◆ EventSuffix_

string SiStripGainFromCalibTree::EventSuffix_

private

Definition at line 315 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

◆ farfromedge

const std::vector<bool>* SiStripGainFromCalibTree::farfromedge = nullptr

private

Definition at line 299 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ farfromedge_token_

edm::EDGetTokenT<std::vector<bool> > SiStripGainFromCalibTree::farfromedge_token_

private

Definition at line 300 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ FirstSetOfConstants

bool SiStripGainFromCalibTree::FirstSetOfConstants

private

Definition at line 145 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ firststrip

const std::vector<unsigned short>* SiStripGainFromCalibTree::firststrip = nullptr

private

Definition at line 291 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ firststrip_token_

edm::EDGetTokenT<std::vector<unsigned short> > SiStripGainFromCalibTree::firststrip_token_

private

Definition at line 292 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ Gains

MonitorElement* SiStripGainFromCalibTree::Gains

private

distribution of gain factors

Definition at line 207 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), qualityMonitor(), and storeOnTree().

◆ gainToken_

edm::ESGetToken<SiStripGain, SiStripGainRcd> SiStripGainFromCalibTree::gainToken_

private

Definition at line 328 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginRun(), and SiStripGainFromCalibTree().

◆ gainused

const std::vector<double>* SiStripGainFromCalibTree::gainused = nullptr

private

Definition at line 309 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ gainused_token_

edm::EDGetTokenT<std::vector<double> > SiStripGainFromCalibTree::gainused_token_

private

Definition at line 310 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ gainusedTick

const std::vector<double>* SiStripGainFromCalibTree::gainusedTick = nullptr

private

Definition at line 311 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ gainusedTick_token_

edm::EDGetTokenT<std::vector<double> > SiStripGainFromCalibTree::gainusedTick_token_

private

Definition at line 312 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ GainVsPrevGainTEC

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTEC

private

Gain vs PreviousGain for TEC

Definition at line 240 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ GainVsPrevGainTIB

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTIB

private

Gain vs PreviousGain for TIB

Definition at line 237 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ GainVsPrevGainTID

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTID

private

Gain vs PreviousGain for TID

Definition at line 238 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ GainVsPrevGainTOB

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTOB

private

Gain vs PreviousGain for TOB

Definition at line 239 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ GOOD

unsigned int SiStripGainFromCalibTree::GOOD

private

Definition at line 252 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), algoComputeMPVandGain(), produceTagFilter(), and storeOnTree().

◆ localdirx

const std::vector<double>* SiStripGainFromCalibTree::localdirx = nullptr

private

Definition at line 285 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ localdirx_token_

edm::EDGetTokenT<std::vector<double> > SiStripGainFromCalibTree::localdirx_token_

private

Definition at line 286 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ localdiry

const std::vector<double>* SiStripGainFromCalibTree::localdiry = nullptr

private

Definition at line 287 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ localdiry_token_

edm::EDGetTokenT<std::vector<double> > SiStripGainFromCalibTree::localdiry_token_

private

Definition at line 288 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ localdirz

const std::vector<double>* SiStripGainFromCalibTree::localdirz = nullptr

private

Definition at line 289 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ localdirz_token_

edm::EDGetTokenT<std::vector<double> > SiStripGainFromCalibTree::localdirz_token_

private

Definition at line 290 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ m_calibrationMode

std::string SiStripGainFromCalibTree::m_calibrationMode

private

Type of statistics for the calibration

Definition at line 155 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), algoBeginJob(), algoBeginRun(), algoComputeMPVandGain(), algoEndJob(), algoEndRun(), isBFieldConsistentWithMode(), processEvent(), produceTagFilter(), qualityMonitor(), SiStripGainFromCalibTree(), storeOnTree(), and swapBFieldMode().

◆ m_calibrationPath

std::string SiStripGainFromCalibTree::m_calibrationPath

private

Definition at line 156 of file SiStripGainFromCalibTree.cc.

Referenced by MakeCalibrationMap(), and SiStripGainFromCalibTree().

◆ m_DQMdir

std::string SiStripGainFromCalibTree::m_DQMdir

private

DQM folder hosting the charge statistics and the monitor plots

Definition at line 157 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), algoEndRun(), and SiStripGainFromCalibTree().

◆ m_harvestingMode

bool SiStripGainFromCalibTree::m_harvestingMode

private

Definition at line 152 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoBeginJob(), algoBeginRun(), algoEndJob(), algoEndRun(), bookDQMHistos(), getNewObject(), and SiStripGainFromCalibTree().

◆ m_splitDQMstat

bool SiStripGainFromCalibTree::m_splitDQMstat

private

Definition at line 153 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), algoEndRun(), and SiStripGainFromCalibTree().

◆ MagFieldCurrentTh

double SiStripGainFromCalibTree::MagFieldCurrentTh

private

Definition at line 132 of file SiStripGainFromCalibTree.cc.

Referenced by isBFieldConsistentWithMode(), and SiStripGainFromCalibTree().

◆ MASKED

unsigned int SiStripGainFromCalibTree::MASKED

private

Definition at line 254 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), algoComputeMPVandGain(), and storeOnTree().

◆ MaxChi2OverNDF

double SiStripGainFromCalibTree::MaxChi2OverNDF

private

Definition at line 135 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree().

◆ MaxMPVError

double SiStripGainFromCalibTree::MaxMPVError

private

Definition at line 134 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree().

◆ MaxNrStrips

unsigned int SiStripGainFromCalibTree::MaxNrStrips

private

Definition at line 140 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ MaxTrackChiOverNdf

double SiStripGainFromCalibTree::MaxTrackChiOverNdf

private

Definition at line 142 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ MaxTrackEta

double SiStripGainFromCalibTree::MaxTrackEta

private

Definition at line 139 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ MaxTrackingIteration

int SiStripGainFromCalibTree::MaxTrackingIteration

private

Definition at line 143 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ MaxTrackMomentum

double SiStripGainFromCalibTree::MaxTrackMomentum

private

Definition at line 137 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ MinNrEntries

double SiStripGainFromCalibTree::MinNrEntries

private

Definition at line 133 of file SiStripGainFromCalibTree.cc.

Referenced by getPeakOfLandau(), and SiStripGainFromCalibTree().

◆ MinTrackEta

double SiStripGainFromCalibTree::MinTrackEta

private

Definition at line 138 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ MinTrackHits

unsigned int SiStripGainFromCalibTree::MinTrackHits

private

Definition at line 141 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ MinTrackMomentum

double SiStripGainFromCalibTree::MinTrackMomentum

private

Definition at line 136 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ MPV_Vs_EtaTEC

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTEC

private

MPV vs Eta for TEC planes

Definition at line 193 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPV_Vs_EtaTECthick

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTECthick

private

MPV vs Eta for TEC tick planes

Definition at line 195 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPV_Vs_EtaTECthin

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTECthin

private

MPV vs Eta for TEC thin planes

Definition at line 194 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPV_Vs_EtaTIB

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTIB

private

MPV vs Eta for TIB planes

Definition at line 190 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPV_Vs_EtaTID

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTID

private

MPV vs Eta for TID planes

Definition at line 191 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPV_Vs_EtaTOB

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTOB

private

MPV vs Eta for TOB planes

Definition at line 192 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPV_Vs_PhiTEC

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTEC

private

MPV vs Phi for TEC planes

Definition at line 200 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPV_Vs_PhiTECthick

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTECthick

private

MPV vs Phi for TID tick planes

Definition at line 202 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPV_Vs_PhiTECthin

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTECthin

private

MPV vs Phi for TEC thin planes

Definition at line 201 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPV_Vs_PhiTIB

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTIB

private

MPV vs Phi for TIB planes

Definition at line 197 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPV_Vs_PhiTID

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTID

private

MPV vs Phi for TID planes

Definition at line 198 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPV_Vs_PhiTOB

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTOB

private

MPV vs Phi for TOB planes

Definition at line 199 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVError

MonitorElement* SiStripGainFromCalibTree::MPVError

private

error of Landau fit

Definition at line 226 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVErrorVsEta

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsEta

private

error of Landau fit vs Eta

Definition at line 228 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVErrorVsMPV

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsMPV

private

error of Landau fit vs MPV

Definition at line 227 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVErrorVsN

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsN

private

error of landau fit vs number of entries

Definition at line 230 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVErrorVsPhi

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsPhi

private

error of Landau fit vs Phi

Definition at line 229 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVs

MonitorElement* SiStripGainFromCalibTree::MPVs

private

distribution of MPVs

Definition at line 208 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVs320

MonitorElement* SiStripGainFromCalibTree::MPVs320

private

distribution of MPVs for thin sensors

Definition at line 209 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVs500

MonitorElement* SiStripGainFromCalibTree::MPVs500

private

distribution of MPVs for tick sensors

Definition at line 210 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTEC

MonitorElement* SiStripGainFromCalibTree::MPVsTEC

private

distribution of MPVs for TEC disks

Definition at line 216 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTECM

MonitorElement* SiStripGainFromCalibTree::MPVsTECM

private

distribution of MPVs for TECM disks

Definition at line 218 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTECM1

MonitorElement* SiStripGainFromCalibTree::MPVsTECM1

private

distribution of MPVs for TECM thin sensors

Definition at line 223 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTECM2

MonitorElement* SiStripGainFromCalibTree::MPVsTECM2

private

distribution of MPVs for TECM tick sensors

Definition at line 224 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTECP

MonitorElement* SiStripGainFromCalibTree::MPVsTECP

private

distribution of MPVs for TECP disks

Definition at line 217 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTECP1

MonitorElement* SiStripGainFromCalibTree::MPVsTECP1

private

distribution of MPVs for TECP thin sensors

Definition at line 221 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTECP2

MonitorElement* SiStripGainFromCalibTree::MPVsTECP2

private

distribution of MPVs for TECP tick sensors

Definition at line 222 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTECthick

MonitorElement* SiStripGainFromCalibTree::MPVsTECthick

private

distribution of MPVs for TEC tick sensors

Definition at line 220 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTECthin

MonitorElement* SiStripGainFromCalibTree::MPVsTECthin

private

distribution of MPVs for TEC thin sensors

Definition at line 219 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTIB

MonitorElement* SiStripGainFromCalibTree::MPVsTIB

private

distribution of MPVs for TIB planes

Definition at line 211 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTID

MonitorElement* SiStripGainFromCalibTree::MPVsTID

private

distribution of MPVs for TID disks

Definition at line 212 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTIDM

MonitorElement* SiStripGainFromCalibTree::MPVsTIDM

private

distribution of MPVs for TIDM disks

Definition at line 214 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTIDP

MonitorElement* SiStripGainFromCalibTree::MPVsTIDP

private

distribution of MPVs for TIDP disks

Definition at line 213 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ MPVsTOB

MonitorElement* SiStripGainFromCalibTree::MPVsTOB

private

distribution of MPVs for TOB planes

Definition at line 215 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ NClusterPixel

unsigned int SiStripGainFromCalibTree::NClusterPixel

private

Definition at line 247 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and processEvent().

◆ NClusterStrip

unsigned int SiStripGainFromCalibTree::NClusterStrip

private

Definition at line 246 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), algoEndJob(), processEvent(), and storeOnTree().

◆ NEvent

unsigned int SiStripGainFromCalibTree::NEvent

private

Definition at line 244 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), algoBeginJob(), processEvent(), and storeOnTree().

◆ newCharge

std::vector<APVGain::APVmon> SiStripGainFromCalibTree::newCharge

private

Definition at line 242 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ NoMPVfit

MonitorElement* SiStripGainFromCalibTree::NoMPVfit

private

R,Z map of missing APV calibration (no fit)

Definition at line 205 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ NoMPVmasked

MonitorElement* SiStripGainFromCalibTree::NoMPVmasked

private

R,Z map of missing APV calibration (masked modules)

Definition at line 204 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

◆ NPixelDets

int SiStripGainFromCalibTree::NPixelDets

private

Definition at line 249 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob().

◆ NStripAPVs

int SiStripGainFromCalibTree::NStripAPVs

private

Definition at line 248 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), algoEndJob(), bookDQMHistos(), produceTagFilter(), and storeOnTree().

◆ nstrips

const std::vector<unsigned short>* SiStripGainFromCalibTree::nstrips = nullptr

private

Definition at line 293 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ nstrips_token_

edm::EDGetTokenT<std::vector<unsigned short> > SiStripGainFromCalibTree::nstrips_token_

private

Definition at line 294 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ NTrack

unsigned int SiStripGainFromCalibTree::NTrack

private

Definition at line 245 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), processEvent(), and storeOnTree().

◆ OldGainRemoving

bool SiStripGainFromCalibTree::OldGainRemoving

private

Definition at line 147 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ OutputGains

std::string SiStripGainFromCalibTree::OutputGains

private

Definition at line 164 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree(), and storeOnTree().

◆ overlapping

const std::vector<bool>* SiStripGainFromCalibTree::overlapping = nullptr

private

Definition at line 297 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ overlapping_token_

edm::EDGetTokenT<std::vector<bool> > SiStripGainFromCalibTree::overlapping_token_

private

Definition at line 298 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ path

const std::vector<double>* SiStripGainFromCalibTree::path = nullptr

private

Definition at line 303 of file SiStripGainFromCalibTree.cc.

Referenced by python.rootplot.rootmath.Target::__repr__(), algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ path_token_

edm::EDGetTokenT<std::vector<double> > SiStripGainFromCalibTree::path_token_

private

Definition at line 304 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ qualityToken_

edm::ESGetToken<SiStripQuality, SiStripQualityRcd> SiStripGainFromCalibTree::qualityToken_

private

Definition at line 329 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginRun(), and SiStripGainFromCalibTree().

◆ rawid

const std::vector<unsigned int>* SiStripGainFromCalibTree::rawid = nullptr

private

Definition at line 283 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ rawid_token_

edm::EDGetTokenT<std::vector<unsigned int> > SiStripGainFromCalibTree::rawid_token_

private

Definition at line 284 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ runInfoToken_

edm::ESGetToken<RunInfo, RunInfoRcd> SiStripGainFromCalibTree::runInfoToken_

private

Definition at line 327 of file SiStripGainFromCalibTree.cc.

Referenced by isBFieldConsistentWithMode(), and SiStripGainFromCalibTree().

◆ runnumber

unsigned int SiStripGainFromCalibTree::runnumber = 0

private

Definition at line 260 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ saturation

const std::vector<bool>* SiStripGainFromCalibTree::saturation = nullptr

private

Definition at line 295 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ saturation_token_

edm::EDGetTokenT<std::vector<bool> > SiStripGainFromCalibTree::saturation_token_

private

Definition at line 296 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ saveSummary

bool SiStripGainFromCalibTree::saveSummary

private

Definition at line 150 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), and SiStripGainFromCalibTree().

◆ SRun

unsigned int SiStripGainFromCalibTree::SRun

private

Definition at line 250 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and processEvent().

◆ tagCondition_GoodFrac

double SiStripGainFromCalibTree::tagCondition_GoodFrac

private

Definition at line 161 of file SiStripGainFromCalibTree.cc.

Referenced by produceTagFilter(), and SiStripGainFromCalibTree().

◆ tagCondition_NClusters

double SiStripGainFromCalibTree::tagCondition_NClusters

private

Definition at line 160 of file SiStripGainFromCalibTree.cc.

Referenced by produceTagFilter(), and SiStripGainFromCalibTree().

◆ tfs

TFileService* SiStripGainFromCalibTree::tfs

private

Definition at line 130 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), and storeOnTree().

◆ tkGeomToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiStripGainFromCalibTree::tkGeomToken_

private

Definition at line 326 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and SiStripGainFromCalibTree().

◆ trackalgo

const std::vector<int>* SiStripGainFromCalibTree::trackalgo = nullptr

private

Definition at line 277 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ trackalgo_token_

edm::EDGetTokenT<std::vector<int> > SiStripGainFromCalibTree::trackalgo_token_

private

Definition at line 278 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ trackchi2ndof

const std::vector<double>* SiStripGainFromCalibTree::trackchi2ndof = nullptr

private

Definition at line 265 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ trackchi2ndof_token_

edm::EDGetTokenT<std::vector<double> > SiStripGainFromCalibTree::trackchi2ndof_token_

private

Definition at line 266 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ tracketa

const std::vector<double>* SiStripGainFromCalibTree::tracketa = nullptr

private

Definition at line 271 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ tracketa_token_

edm::EDGetTokenT<std::vector<double> > SiStripGainFromCalibTree::tracketa_token_

private

Definition at line 272 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ trackhitsvalid

const std::vector<unsigned int>* SiStripGainFromCalibTree::trackhitsvalid = nullptr

private

Definition at line 275 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ trackhitsvalid_token_

edm::EDGetTokenT<std::vector<unsigned int> > SiStripGainFromCalibTree::trackhitsvalid_token_

private

Definition at line 276 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ trackindex

const std::vector<int>* SiStripGainFromCalibTree::trackindex = nullptr

private

Definition at line 281 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ trackindex_token_

edm::EDGetTokenT<std::vector<int> > SiStripGainFromCalibTree::trackindex_token_

private

Definition at line 282 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ trackp

const std::vector<float>* SiStripGainFromCalibTree::trackp = nullptr

private

Definition at line 267 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), algoAnalyzeTheTree(), and processEvent().

◆ trackp_token_

edm::EDGetTokenT<std::vector<float> > SiStripGainFromCalibTree::trackp_token_

private

Definition at line 268 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ trackphi

const std::vector<double>* SiStripGainFromCalibTree::trackphi = nullptr

private

Definition at line 273 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ trackphi_token_

edm::EDGetTokenT<std::vector<double> > SiStripGainFromCalibTree::trackphi_token_

private

Definition at line 274 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ TrackPrefix_

string SiStripGainFromCalibTree::TrackPrefix_

private

Definition at line 316 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

◆ trackpt

const std::vector<float>* SiStripGainFromCalibTree::trackpt = nullptr

private

Definition at line 269 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ trackpt_token_

edm::EDGetTokenT<std::vector<float> > SiStripGainFromCalibTree::trackpt_token_

private

Definition at line 270 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ TrackSuffix_

string SiStripGainFromCalibTree::TrackSuffix_

private

Definition at line 317 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

◆ TrigTech

const std::vector<bool>* SiStripGainFromCalibTree::TrigTech = nullptr

private

Definition at line 261 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

◆ TrigTech_token_

edm::EDGetTokenT<std::vector<bool> > SiStripGainFromCalibTree::TrigTech_token_

private

Definition at line 262 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

◆ tTopo_

const TrackerTopology* SiStripGainFromCalibTree::tTopo_ = nullptr

private

Definition at line 330 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), processEvent(), and qualityMonitor().

◆ tTopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiStripGainFromCalibTree::tTopoToken_

private

Definition at line 325 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and SiStripGainFromCalibTree().

◆ useCalibration

bool SiStripGainFromCalibTree::useCalibration

private

Definition at line 151 of file SiStripGainFromCalibTree.cc.

Referenced by MakeCalibrationMap(), processEvent(), and SiStripGainFromCalibTree().

◆ Validation

bool SiStripGainFromCalibTree::Validation

private

Definition at line 146 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

◆ VChargeHisto

std::vector<std::string> SiStripGainFromCalibTree::VChargeHisto

private

Charge monitor plots to be output

Definition at line 158 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndRun(), bookDQMHistos(), and SiStripGainFromCalibTree().

◆ VInputFiles

vector<string> SiStripGainFromCalibTree::VInputFiles

private

Definition at line 165 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ DQMStore

◆ MonitorElement

Constructor & Destructor Documentation

◆ SiStripGainFromCalibTree()

◆ ~SiStripGainFromCalibTree()

Member Function Documentation

◆ algoAnalyze()

◆ algoAnalyzeTheTree()

◆ algoBeginJob()

◆ algoBeginRun()

◆ algoComputeMPVandGain()

◆ algoEndJob()

◆ algoEndRun()

◆ bookDQMHistos()

◆ connect()

◆ getNewObject()

◆ getPeakOfLandau()

◆ isBFieldConsistentWithMode()

◆ IsGoodLandauFit()

◆ MakeCalibrationMap()

◆ merge()

◆ processEvent()

◆ produceTagFilter()

◆ qualityMonitor()

◆ statCollectionFromMode()

◆ storeOnTree()

◆ swapBFieldMode()

Member Data Documentation

◆ AlgoMode

◆ AllowSaturation

◆ amplitude

◆ amplitude_token_

◆ APVsColl

◆ APVsCollOrdered

◆ BAD

◆ booked_dir_

◆ CalibPrefix_

◆ CalibrationLevel

◆ CalibSuffix_

◆ charge

◆ Charge_1

◆ Charge_2

◆ Charge_3

◆ Charge_4

◆ charge_token_

◆ Charge_Vs_Index

◆ Charge_Vs_PathlengthTECM1

◆ Charge_Vs_PathlengthTECM2

◆ Charge_Vs_PathlengthTECP1

◆ Charge_Vs_PathlengthTECP2

◆ Charge_Vs_PathlengthTIB

◆ Charge_Vs_PathlengthTIDM

◆ Charge_Vs_PathlengthTIDP

◆ Charge_Vs_PathlengthTOB

◆ chargeoverpath

◆ chargeoverpath_token_

◆ dbe

◆ DiffWRTPrevGainTEC

◆ DiffWRTPrevGainTIB

◆ DiffWRTPrevGainTID

◆ DiffWRTPrevGainTOB

◆ doChargeMonitorPerPlane

◆ dqm_tag_

◆ ERun

◆ eventnumber

◆ EventPrefix_

◆ EventSuffix_

◆ farfromedge

◆ farfromedge_token_

◆ FirstSetOfConstants

◆ firststrip

◆ firststrip_token_

◆ Gains