Inheritance diagram for SiStripGainFromCalibTree:

Public Types
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

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 = 0

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 = 0

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 = 0

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 = 0

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

bool	FirstSetOfConstants

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

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

MonitorElement *	Gains

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

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

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

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

MonitorElement *	GainVsPrevGainTEC

MonitorElement *	GainVsPrevGainTIB

MonitorElement *	GainVsPrevGainTID

MonitorElement *	GainVsPrevGainTOB

unsigned int	GOOD

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

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

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

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

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

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 = 0

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

unsigned int	NTrack

bool	OldGainRemoving

std::string	OutputGains

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

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

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

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

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

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

unsigned int	runnumber = 0

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

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

bool	saveSummary

unsigned int	SRun

double	tagCondition_GoodFrac

double	tagCondition_NClusters

TFileService *	tfs

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

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

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

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

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

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

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

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

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

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

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

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

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

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

string	TrackPrefix_

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

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

string	TrackSuffix_

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

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

const TrackerTopology *	tTopo_

bool	useCalibration

bool	Validation

std::vector< std::string >	VChargeHisto

vector< string >	VInputFiles

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from 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)

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<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

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

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

Detailed Description

Definition at line 90 of file SiStripGainFromCalibTree.cc.

Member Typedef Documentation

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

Definition at line 93 of file SiStripGainFromCalibTree.cc.

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

Definition at line 92 of file SiStripGainFromCalibTree.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 355 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->();
   dbe->setVerbose(10);
 
   //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_));
 
   tTopo_ = nullptr;
 }

SiStripGainFromCalibTree::~SiStripGainFromCalibTree ( )

override

Definition at line 1801 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

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

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 1835 of file SiStripGainFromCalibTree.cc.

References AlgoMode, amplitude, amplitude_token_, charge, charge_token_, chargeoverpath, chargeoverpath_token_, connect(), DEFINE_FWK_MODULE, edm::EventID::event(), eventnumber, farfromedge, farfromedge_token_, firststrip, firststrip_token_, gainused, gainused_token_, gainusedTick, gainusedTick_token_, edm::EventBase::id(), localdirx, localdirx_token_, localdiry, localdiry_token_, localdirz, localdirz_token_, m_harvestingMode, nstrips, nstrips_token_, overlapping, overlapping_token_, path, path_token_, processEvent(), rawid, rawid_token_, edm::EventID::run(), 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();
 }

void SiStripGainFromCalibTree::algoAnalyzeTheTree ( )

private

Definition at line 1307 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
   }
 }

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

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 649 of file SiStripGainFromCalibTree.cc.

References AlgoMode, sistrip::APV, stAPVGain::APVId, APVsColl, APVsCollOrdered, BAD, stAPVGain::Bin, bookDQMHistos(), stAPVGain::DetId, TrackerGeometry::dets(), dqm_tag_, HTMLExport::elem(), ERun, stAPVGain::Eta, stAPVGain::FitChi2, stAPVGain::FitMPV, stAPVGain::FitMPVErr, stAPVGain::FitNorm, stAPVGain::FitWidth, stAPVGain::FitWidthErr, stAPVGain::Gain, edm::EventSetup::get(), GOOD, mps_fire::i, stAPVGain::Index, isBFieldConsistentWithMode(), stAPVGain::isMasked, dqmiolumiharvest::j, m_calibrationMode, m_DQMdir, m_harvestingMode, m_splitDQMstat, MakeCalibrationMap(), MASKED, NClusterPixel, NClusterStrip, PixelTopology::ncolumns(), stAPVGain::NEntries, NEvent, NPixelDets, PixelTopology::nrows(), NStripAPVs, StripTopology::nstrips(), NTrack, or, stAPVGain::Phi, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, stAPVGain::PreviousGain, stAPVGain::PreviousGainTick, stAPVGain::R, DetId::rawId(), SRun, statCollectionFromMode(), AlCaHLTBitMon_QueryRunRegistry::string, stAPVGain::SubDet, ntupleEnum::SubDet, DetId::subdetId(), swapBFieldMode(), StripSubdetector::TEC, stAPVGain::Thickness, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, tTopo_, stAPVGain::x, stAPVGain::y, and stAPVGain::z.

                                                                        {
   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());
   }
 
   edm::ESHandle<TrackerTopology> TopoHandle;
   iSetup.get<TrackerTopologyRcd>().get(TopoHandle);
   tTopo_ = TopoHandle.product();
 
   edm::ESHandle<TrackerGeometry> tkGeom;
   iSetup.get<TrackerDigiGeometryRecord>().get(tkGeom);
   auto const& Det = tkGeom->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;
 }

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

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 808 of file SiStripGainFromCalibTree.cc.

References a, AlgoMode, sistrip::APV, stAPVGain::APVId, APVsCollOrdered, stAPVGain::DetId, beamvalidation::exit(), edm::EventSetup::get(), SiStripGain::getApvGain(), SiStripGain::getNumberOfTags(), SiStripGain::getRange(), SiStripQuality::IsApvBad(), isBFieldConsistentWithMode(), stAPVGain::isMasked, edm::ESHandleBase::isValid(), m_calibrationMode, m_harvestingMode, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, stAPVGain::PreviousGain, stAPVGain::PreviousGainTick, mps_setup::stdout, stAPVGain::SubDet, 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;
     }
   }
 
   edm::ESHandle<SiStripGain> gainHandle;
   iSetup.get<SiStripGainRcd>().get(gainHandle);
   if (!gainHandle.isValid()) {
     edm::LogError("SiStripGainFromCalibTree") << "gainHandle is not valid\n";
     exit(0);
   }
 
   edm::ESHandle<SiStripQuality> SiStripQuality_;
   iSetup.get<SiStripQualityRcd>().get(SiStripQuality_);
 
   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);
     //      }
   }
 }

void SiStripGainFromCalibTree::algoComputeMPVandGain ( )

private

Definition at line 1363 of file SiStripGainFromCalibTree.cc.

References AlgoMode, sistrip::APV, stAPVGain::APVId, APVsColl, BAD, stAPVGain::Bin, CalibrationLevel, Charge_Vs_Index, stAPVGain::DetId, stAPVGain::FitChi2, stAPVGain::FitMPV, stAPVGain::FitMPVErr, stAPVGain::FitNorm, stAPVGain::FitWidth, stAPVGain::FitWidthErr, stAPVGain::Gain, getPeakOfLandau(), GOOD, Harvest, Exhume::I, mps_fire::i, stAPVGain::Index, IsGoodLandauFit(), stAPVGain::isMasked, m_calibrationMode, MASKED, stAPVGain::NEntries, stAPVGain::PreviousGain, statCollectionFromMode(), mps_setup::stdout, and stAPVGain::SubDet.

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");
 }

void SiStripGainFromCalibTree::algoEndJob ( )

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 1025 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);
   }
 }

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

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 865 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, 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, 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, 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, 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());
       }
     }
   }
 }

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

private

Definition at line 472 of file SiStripGainFromCalibTree.cc.

References a, AlgoMode, b, dqm::dqmstoreimpl::DQMStore::book1DD(), dqm::dqmstoreimpl::DQMStore::book2DD(), dqm::dqmstoreimpl::DQMStore::book2S(), booked_dir_, Charge_1, Charge_2, Charge_3, Charge_4, 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, dbe, DiffWRTPrevGainTEC, DiffWRTPrevGainTIB, DiffWRTPrevGainTID, DiffWRTPrevGainTOB, doChargeMonitorPerPlane, JetChargeProducer_cfi::exp, dqmdumpme::first, Gains, GainVsPrevGainTEC, GainVsPrevGainTIB, GainVsPrevGainTID, GainVsPrevGainTOB, Harvest, mps_fire::i, dqm-mbProfile::log, m_harvestingMode, APVGain::monHnames(), beam_dqm_sourceclient-live_cfg::monitor, 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, NStripAPVs, p1, p2, edm::second(), dqm::dqmstoreimpl::DQMStore::setCurrentFolder(), statCollectionFromMode(), AlCaHLTBitMon_QueryRunRegistry::string, APVGain::subdetectorId(), APVGain::subdetectorPlane(), APVGain::subdetectorSide(), VChargeHisto, and y.

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 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(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 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(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 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(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 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(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 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(id, side, plane, monitor));
     }
   }
 }

template<typename T >

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

inlineprivate

Definition at line 123 of file SiStripGainFromCalibTree.cc.

References edm::Event::getByToken(), patZpeak::handle, and edm::Handle< T >::product().

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

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

overrideprivatevirtual

Implements ConditionDBWriter< SiStripApvGain >.

Definition at line 1756 of file SiStripGainFromCalibTree.cc.

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

                                                                      {
   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;
     }
     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;
 }

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

private

Definition at line 1081 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;
 }

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

private

Definition at line 789 of file SiStripGainFromCalibTree.cc.

References edm::EventSetup::get(), RunInfo::m_avg_current, m_calibrationMode, MagFieldCurrentTh, and edm::ESHandle< T >::product().

Referenced by algoBeginJob(), and algoBeginRun().

                                                                                            {
   edm::ESHandle<RunInfo> runInfo;
   iSetup.get<RunInfoRcd>().get(runInfo);
 
   double average_current = runInfo.product()->m_avg_current;
   bool isOn = (average_current > MagFieldCurrentTh);
   bool is0T = (m_calibrationMode.substr(m_calibrationMode.length() - 2) == "0T");
 
   return ((isOn && !is0T) || (!isOn && is0T));
 }

bool SiStripGainFromCalibTree::IsGoodLandauFit ( double * FitResults )

private

Definition at line 1108 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;
 }

void SiStripGainFromCalibTree::MakeCalibrationMap ( )

private

Definition at line 1811 of file SiStripGainFromCalibTree.cc.

References sistrip::APV, APVsColl, stAPVGain::CalibGain, 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;
 }

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

private

Definition at line 343 of file SiStripGainFromCalibTree.cc.

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));
       }
     }
   }
 }

void SiStripGainFromCalibTree::processEvent ( )

private

Definition at line 1116 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()

bool SiStripGainFromCalibTree::produceTagFilter ( )

private

Definition at line 1729 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;
 }

void SiStripGainFromCalibTree::qualityMonitor ( )

private

Definition at line 1468 of file SiStripGainFromCalibTree.cc.

References a, AlgoMode, sistrip::APV, APVsCollOrdered, newFWLiteAna::bin, Charge_Vs_Index, stAPVGain::DetId, DiffWRTPrevGainTEC, DiffWRTPrevGainTIB, DiffWRTPrevGainTID, DiffWRTPrevGainTOB, MillePedeFileConverter_cfg::e, stAPVGain::Eta, APVGain::FetchMonitor(), dqm::impl::MonitorElement::Fill(), stAPVGain::FitMPV, stAPVGain::FitMPVErr, IntegrityClient_cfi::Gain, stAPVGain::Gain, Gains, GainVsPrevGainTEC, GainVsPrevGainTIB, GainVsPrevGainTID, GainVsPrevGainTOB, h, Harvest, stAPVGain::Index, stAPVGain::isMasked, 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, stAPVGain::NEntries, newCharge, NoMPVfit, NoMPVmasked, stAPVGain::Phi, VtxSmearedParameters_cfi::Phi, stAPVGain::PreviousGain, stAPVGain::R, dttmaxenums::R, statCollectionFromMode(), stAPVGain::SubDet, ntupleEnum::SubDet, stAPVGain::Thickness, tTopo_, z, and stAPVGain::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);
   }
 }

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

inlineprivate

Definition at line 329 of file SiStripGainFromCalibTree.cc.

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

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;
 }

void SiStripGainFromCalibTree::storeOnTree ( TFileService * tfs )

private

Definition at line 1600 of file SiStripGainFromCalibTree.cc.

References a, sistrip::APV, stAPVGain::APVId, APVsCollOrdered, BAD, Charge_Vs_Index, stAPVGain::DetId, stAPVGain::Eta, stAPVGain::FitChi2, stAPVGain::FitMPV, stAPVGain::FitMPVErr, stAPVGain::FitNorm, stAPVGain::FitWidth, stAPVGain::FitWidthErr, stAPVGain::Gain, Gains, GOOD, stAPVGain::Index, stAPVGain::isMasked, m_calibrationMode, TFileService::make(), MASKED, NClusterStrip, stAPVGain::NEntries, NEvent, NStripAPVs, NTrack, OutputGains, stAPVGain::Phi, stAPVGain::PreviousGain, stAPVGain::PreviousGainTick, stAPVGain::R, statCollectionFromMode(), mps_setup::stdout, stAPVGain::SubDet, stAPVGain::Thickness, stAPVGain::x, stAPVGain::y, and stAPVGain::z.

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);
 }

void SiStripGainFromCalibTree::swapBFieldMode ( void )

private

Definition at line 800 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

std::string SiStripGainFromCalibTree::AlgoMode

private

Definition at line 165 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::AllowSaturation

private

Definition at line 146 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

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

private

Definition at line 309 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 310 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 325 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 324 of file SiStripGainFromCalibTree.cc.

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

unsigned int SiStripGainFromCalibTree::BAD

private

Definition at line 255 of file SiStripGainFromCalibTree.cc.

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

std::string SiStripGainFromCalibTree::booked_dir_

private

Definition at line 172 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos().

string SiStripGainFromCalibTree::CalibPrefix_

private

Definition at line 320 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

int SiStripGainFromCalibTree::CalibrationLevel

private

Definition at line 150 of file SiStripGainFromCalibTree.cc.

Referenced by algoComputeMPVandGain(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::CalibSuffix_

private

Definition at line 321 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 303 of file SiStripGainFromCalibTree.cc.

Referenced by Vispa.Plugins.EdmBrowser.ParticleDataList.ParticleData::__repr__(), algoAnalyze(), and algoAnalyzeTheTree().

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

private

Charge per cm per layer / wheel

Definition at line 175 of file SiStripGainFromCalibTree.cc.

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

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

private

Charge per cm per layer / wheel without G2

Definition at line 176 of file SiStripGainFromCalibTree.cc.

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

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

private

Charge per cm per layer / wheel without G1

Definition at line 177 of file SiStripGainFromCalibTree.cc.

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

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

private

Charge per cm per layer / wheel without G1 and G1

Definition at line 178 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 304 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Charge per cm for each detector id

Definition at line 174 of file SiStripGainFromCalibTree.cc.

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

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

private

Charge vs pathlength in TECP thin

Definition at line 186 of file SiStripGainFromCalibTree.cc.

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

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

private

Charge vs pathlength in TECP thick

Definition at line 187 of file SiStripGainFromCalibTree.cc.

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

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

private

Charge vs pathlength in TECP thin

Definition at line 184 of file SiStripGainFromCalibTree.cc.

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

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

private

Charge vs pathlength in TECP thick

Definition at line 185 of file SiStripGainFromCalibTree.cc.

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

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

private

Charge vs pathlength in TIB

Definition at line 180 of file SiStripGainFromCalibTree.cc.

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

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

private

Charge vs pathlength in TIDM

Definition at line 183 of file SiStripGainFromCalibTree.cc.

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

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

private

Charge vs pathlength in TIDP

Definition at line 182 of file SiStripGainFromCalibTree.cc.

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

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

private

Charge vs pathlength in TOB

Definition at line 181 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 307 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 308 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

DQMStore* SiStripGainFromCalibTree::dbe

private

Definition at line 133 of file SiStripGainFromCalibTree.cc.

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

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTEC

private

ratio gain / PreviousGain for TEC disks

Definition at line 237 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTIB

private

ratio Gain / PreviousGain for TIB layers

Definition at line 234 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTID

private

ratio Gain / PreviousGain for TID disks

Definition at line 235 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTOB

private

ratio Gain / PreviousGain for TOB layers

Definition at line 236 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

bool SiStripGainFromCalibTree::doChargeMonitorPerPlane

private

Charge monitor per detector plane

Definition at line 156 of file SiStripGainFromCalibTree.cc.

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

std::vector<string> SiStripGainFromCalibTree::dqm_tag_

private

Definition at line 171 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::ERun

private

Definition at line 253 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and processEvent().

unsigned int SiStripGainFromCalibTree::eventnumber = 0

private

Definition at line 261 of file SiStripGainFromCalibTree.cc.

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

string SiStripGainFromCalibTree::EventPrefix_

private

Definition at line 316 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::EventSuffix_

private

Definition at line 317 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 301 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 302 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

bool SiStripGainFromCalibTree::FirstSetOfConstants

private

Definition at line 147 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

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

private

Definition at line 293 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 294 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

MonitorElement* SiStripGainFromCalibTree::Gains

private

distribution of gain factors

Definition at line 209 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 311 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 312 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 313 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 314 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTEC

private

Gain vs PreviousGain for TEC

Definition at line 242 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTIB

private

Gain vs PreviousGain for TIB

Definition at line 239 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTID

private

Gain vs PreviousGain for TID

Definition at line 240 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTOB

private

Gain vs PreviousGain for TOB

Definition at line 241 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

unsigned int SiStripGainFromCalibTree::GOOD

private

Definition at line 254 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 287 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 288 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 289 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 290 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 291 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 292 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

std::string SiStripGainFromCalibTree::m_calibrationMode

private

Type of statistics for the calibration

Definition at line 157 of file SiStripGainFromCalibTree.cc.

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

std::string SiStripGainFromCalibTree::m_calibrationPath

private

Definition at line 158 of file SiStripGainFromCalibTree.cc.

Referenced by MakeCalibrationMap(), and SiStripGainFromCalibTree().

std::string SiStripGainFromCalibTree::m_DQMdir

private

DQM folder hosting the charge statistics and the monitor plots

Definition at line 159 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::m_harvestingMode

private

Definition at line 154 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::m_splitDQMstat

private

Definition at line 155 of file SiStripGainFromCalibTree.cc.

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

double SiStripGainFromCalibTree::MagFieldCurrentTh

private

Definition at line 134 of file SiStripGainFromCalibTree.cc.

Referenced by isBFieldConsistentWithMode(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::MASKED

private

Definition at line 256 of file SiStripGainFromCalibTree.cc.

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

double SiStripGainFromCalibTree::MaxChi2OverNDF

private

Definition at line 137 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxMPVError

private

Definition at line 136 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::MaxNrStrips

private

Definition at line 142 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxTrackChiOverNdf

private

Definition at line 144 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxTrackEta

private

Definition at line 141 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

int SiStripGainFromCalibTree::MaxTrackingIteration

private

Definition at line 145 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxTrackMomentum

private

Definition at line 139 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MinNrEntries

private

Definition at line 135 of file SiStripGainFromCalibTree.cc.

Referenced by getPeakOfLandau(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MinTrackEta

private

Definition at line 140 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::MinTrackHits

private

Definition at line 143 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MinTrackMomentum

private

Definition at line 138 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTEC

private

MPV vs Eta for TEC planes

Definition at line 195 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTECthick

private

MPV vs Eta for TEC tick planes

Definition at line 197 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTECthin

private

MPV vs Eta for TEC thin planes

Definition at line 196 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTIB

private

MPV vs Eta for TIB planes

Definition at line 192 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTID

private

MPV vs Eta for TID planes

Definition at line 193 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTOB

private

MPV vs Eta for TOB planes

Definition at line 194 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTEC

private

MPV vs Phi for TEC planes

Definition at line 202 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTECthick

private

MPV vs Phi for TID tick planes

Definition at line 204 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTECthin

private

MPV vs Phi for TEC thin planes

Definition at line 203 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTIB

private

MPV vs Phi for TIB planes

Definition at line 199 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTID

private

MPV vs Phi for TID planes

Definition at line 200 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTOB

private

MPV vs Phi for TOB planes

Definition at line 201 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVError

private

error of Landau fit

Definition at line 228 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsEta

private

error of Landau fit vs Eta

Definition at line 230 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsMPV

private

error of Landau fit vs MPV

Definition at line 229 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsN

private

error of landau fit vs number of entries

Definition at line 232 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsPhi

private

error of Landau fit vs Phi

Definition at line 231 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVs

private

distribution of MPVs

Definition at line 210 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVs320

private

distribution of MPVs for thin sensors

Definition at line 211 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVs500

private

distribution of MPVs for tick sensors

Definition at line 212 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTEC

private

distribution of MPVs for TEC disks

Definition at line 218 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECM

private

distribution of MPVs for TECM disks

Definition at line 220 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECM1

private

distribution of MPVs for TECM thin sensors

Definition at line 225 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECM2

private

distribution of MPVs for TECM tick sensors

Definition at line 226 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECP

private

distribution of MPVs for TECP disks

Definition at line 219 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECP1

private

distribution of MPVs for TECP thin sensors

Definition at line 223 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECP2

private

distribution of MPVs for TECP tick sensors

Definition at line 224 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECthick

private

distribution of MPVs for TEC tick sensors

Definition at line 222 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECthin

private

distribution of MPVs for TEC thin sensors

Definition at line 221 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTIB

private

distribution of MPVs for TIB planes

Definition at line 213 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTID

private

distribution of MPVs for TID disks

Definition at line 214 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTIDM

private

distribution of MPVs for TIDM disks

Definition at line 216 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTIDP

private

distribution of MPVs for TIDP disks

Definition at line 215 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTOB

private

distribution of MPVs for TOB planes

Definition at line 217 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

unsigned int SiStripGainFromCalibTree::NClusterPixel

private

Definition at line 249 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and processEvent().

unsigned int SiStripGainFromCalibTree::NClusterStrip

private

Definition at line 248 of file SiStripGainFromCalibTree.cc.

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

unsigned int SiStripGainFromCalibTree::NEvent

private

Definition at line 246 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 244 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::NoMPVfit

private

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

Definition at line 207 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::NoMPVmasked

private

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

Definition at line 206 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

int SiStripGainFromCalibTree::NPixelDets

private

Definition at line 251 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob().

int SiStripGainFromCalibTree::NStripAPVs

private

Definition at line 250 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 295 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 296 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::NTrack

private

Definition at line 247 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::OldGainRemoving

private

Definition at line 149 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

std::string SiStripGainFromCalibTree::OutputGains

private

Definition at line 166 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree(), and storeOnTree().

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

private

Definition at line 299 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 300 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 305 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 306 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 285 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 286 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::runnumber = 0

private

Definition at line 262 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 297 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 298 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

bool SiStripGainFromCalibTree::saveSummary

private

Definition at line 152 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::SRun

private

Definition at line 252 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and processEvent().

double SiStripGainFromCalibTree::tagCondition_GoodFrac

private

Definition at line 163 of file SiStripGainFromCalibTree.cc.

Referenced by produceTagFilter(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::tagCondition_NClusters

private

Definition at line 162 of file SiStripGainFromCalibTree.cc.

Referenced by produceTagFilter(), and SiStripGainFromCalibTree().

TFileService* SiStripGainFromCalibTree::tfs

private

Definition at line 132 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob().

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

private

Definition at line 279 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 280 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 267 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 268 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 273 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 274 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 277 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 278 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 283 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 284 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 269 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 270 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 275 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 276 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::TrackPrefix_

private

Definition at line 318 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 271 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 272 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::TrackSuffix_

private

Definition at line 319 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 263 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 264 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

const TrackerTopology* SiStripGainFromCalibTree::tTopo_

private

Definition at line 326 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::useCalibration

private

Definition at line 153 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::Validation

private

Definition at line 148 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

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

private

Charge monitor plots to be output

Definition at line 160 of file SiStripGainFromCalibTree.cc.

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

vector<string> SiStripGainFromCalibTree::VInputFiles

private

Definition at line 167 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

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation