user includes More...

#include <SiStripGainsPCLWorker.h>

Inheritance diagram for SiStripGainsPCLWorker:

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

Public Member Functions inherited from DQMEDAnalyzer
void	beginRun (edm::Run const &, edm::EventSetup const &) final

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer (void)

void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

	~EDAnalyzerBase () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

virtual	~EDConsumerBase () noexcept(false)

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

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

virtual void	beginJob ()

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

virtual void	checkBookAPVColls (const edm::EventSetup &setup)

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

virtual void	endJob ()

void	processEvent (const TrackerTopology *topo)

int	statCollectionFromMode (const char *tag) const

Private Attributes
bool	AllowSaturation

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

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

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

std::vector< std::shared_ptr< stAPVGain > >	APVsCollOrdered

const TrackerGeometry *	bareTkGeomPtr_

std::string	CalibPrefix_

std::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_

bool	doChargeMonitorPerPlane

std::vector< std::string >	dqm_tag_

unsigned int	ERun

unsigned int	eventnumber =0

std::string	EventPrefix_

std::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_

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_

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_DQMdir

unsigned int	MaxNrStrips

double	MaxTrackChiOverNdf

double	MaxTrackEta

int	MaxTrackingIteration

double	MaxTrackMomentum

double	MinTrackEta

unsigned int	MinTrackHits

double	MinTrackMomentum

unsigned int	NClusterPixel

unsigned int	NClusterStrip

unsigned int	NEvent

int	NPixelDets

int	NStripAPVs

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

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

unsigned int	NTrack

bool	OldGainRemoving

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_

unsigned int	SRun

edm::ESHandle< TrackerGeometry >	tkGeom_

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_

std::string	TrackPrefix_

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

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

std::string	TrackSuffix_

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

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

bool	useCalibration

bool	Validation

std::vector< std::string >	VChargeHisto

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

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

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

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

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

void	consumesMany (const TypeToGet &id)

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

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

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

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

Detailed Description

user includes

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

Definition at line 75 of file SiStripGainsPCLWorker.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 9 of file SiStripGainsPCLWorker.cc.

                                                                            :
   NEvent(0),
   NTrack(0),
   NClusterStrip(0),
   NClusterPixel(0),
   NStripAPVs(0),
   NPixelDets(0),
   SRun(1<<31),
   ERun(0),
   bareTkGeomPtr_(nullptr)
 {
   
   MinTrackMomentum        = iConfig.getUntrackedParameter<double>  ("minTrackMomentum"        ,  3.0);
   MaxTrackMomentum        = iConfig.getUntrackedParameter<double>  ("maxTrackMomentum"        ,  99999.0);
   MinTrackEta             = iConfig.getUntrackedParameter<double>  ("minTrackEta"             , -5.0);
   MaxTrackEta             = iConfig.getUntrackedParameter<double>  ("maxTrackEta"             ,  5.0);
   MaxNrStrips             = iConfig.getUntrackedParameter<unsigned>("maxNrStrips"             ,  2);
   MinTrackHits            = iConfig.getUntrackedParameter<unsigned>("MinTrackHits"            ,  8);
   MaxTrackChiOverNdf      = iConfig.getUntrackedParameter<double>  ("MaxTrackChiOverNdf"      ,  3);
   MaxTrackingIteration    = iConfig.getUntrackedParameter<int>     ("MaxTrackingIteration"    ,  7);
   AllowSaturation         = iConfig.getUntrackedParameter<bool>    ("AllowSaturation"         ,  false);
   FirstSetOfConstants     = iConfig.getUntrackedParameter<bool>    ("FirstSetOfConstants"     ,  true);
   Validation              = iConfig.getUntrackedParameter<bool>    ("Validation"              ,  false);
   OldGainRemoving         = iConfig.getUntrackedParameter<bool>    ("OldGainRemoving"         ,  false);
   useCalibration          = iConfig.getUntrackedParameter<bool>    ("UseCalibration"          ,  false);
   doChargeMonitorPerPlane = iConfig.getUntrackedParameter<bool>    ("doChargeMonitorPerPlane" ,  false);
   m_DQMdir                = iConfig.getUntrackedParameter<std::string>  ("DQMdir"             , "AlCaReco/SiStripGains");
   m_calibrationMode       = iConfig.getUntrackedParameter<std::string>  ("calibrationMode"    , "StdBunch");
   VChargeHisto            = iConfig.getUntrackedParameter<std::vector<std::string> >  ("ChargeHisto");
 
   //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_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");
 
   std::string label = swhallowgain_pset.getUntrackedParameter<std::string>("label");
   CalibPrefix_ = swhallowgain_pset.getUntrackedParameter<std::string>("prefix");
   CalibSuffix_ = swhallowgain_pset.getUntrackedParameter<std::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<std::string>("label");
   EventPrefix_ = evtinfo_pset.getUntrackedParameter<std::string>("prefix");
   EventSuffix_ = evtinfo_pset.getUntrackedParameter<std::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<std::string>("label");
   TrackPrefix_ = track_pset.getUntrackedParameter<std::string>("prefix");
   TrackSuffix_ = track_pset.getUntrackedParameter<std::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_));
   
 }

Member Function Documentation

void SiStripGainsPCLWorker::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 147 of file SiStripGainsPCLWorker.cc.

References amplitude, amplitude_token_, charge, charge_token_, chargeoverpath, chargeoverpath_token_, connect(), edm::EventID::event(), eventnumber, farfromedge, farfromedge_token_, firststrip, firststrip_token_, gainused, gainused_token_, gainusedTick, gainusedTick_token_, edm::EventSetup::get(), edm::EventBase::id(), localdirx, localdirx_token_, localdiry, localdiry_token_, localdirz, localdirz_token_, nstrips, nstrips_token_, overlapping, overlapping_token_, path, path_token_, processEvent(), edm::ESHandle< T >::product(), 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_.

 {
   using namespace edm;
 
   //  this->checkBookAPVColls(iSetup); // check whether APV colls are booked and do so if not yet done
 
   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);
  
   edm::ESHandle<TrackerTopology> TopoHandle;
   iSetup.get<TrackerTopologyRcd>().get( TopoHandle );
   const TrackerTopology* topo = TopoHandle.product();
 
   processEvent(topo);
 
 }

void SiStripGainsPCLWorker::beginJob ( void )

privatevirtual

Definition at line 356 of file SiStripGainsPCLWorker.cc.

357 {

358 }

void SiStripGainsPCLWorker::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	run,
		edm::EventSetup const &	es
	)

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 493 of file SiStripGainsPCLWorker.cc.

                                                                                                             {
 
   ibooker.cd();
   std::string dqm_dir = m_DQMdir;
   const char* tag = dqm_tag_[statCollectionFromMode(m_calibrationMode.c_str())].c_str();
 
   edm::LogInfo("SiStripGainsPCLWorker") << "Setting " << dqm_dir << "in DQM and booking histograms for tag "
                     << tag << std::endl;
   
   ibooker.setCurrentFolder(dqm_dir);
  
   std::string stag(tag);
   if(!stag.empty() && stag[0]!='_') stag.insert(0,1,'_');
   
   std::string cvi      = std::string("Charge_Vs_Index") + stag;
   std::string cvpTIB   = std::string("Charge_Vs_PathlengthTIB")   + stag;
   std::string cvpTOB   = std::string("Charge_Vs_PathlengthTOB")   + stag;
   std::string cvpTIDP  = std::string("Charge_Vs_PathlengthTIDP")  + stag;
   std::string cvpTIDM  = std::string("Charge_Vs_PathlengthTIDM")  + stag;
   std::string cvpTECP1 = std::string("Charge_Vs_PathlengthTECP1") + stag;
   std::string cvpTECP2 = std::string("Charge_Vs_PathlengthTECP2") + stag;
   std::string cvpTECM1 = std::string("Charge_Vs_PathlengthTECM1") + stag;
   std::string cvpTECM2 = std::string("Charge_Vs_PathlengthTECM2") + stag;
   
   int elepos = statCollectionFromMode(tag);
 
   // 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]           = ibooker.book2S(cvi.c_str()     , cvi.c_str()     , NStripAPVs, &binXarray[0], 687, binYarray.data());
   Charge_Vs_PathlengthTIB[elepos]   = ibooker.book2S(cvpTIB.c_str()  , cvpTIB.c_str()  , 20   , 0.3 , 1.3  , 250,0,2000);
   Charge_Vs_PathlengthTOB[elepos]   = ibooker.book2S(cvpTOB.c_str()  , cvpTOB.c_str()  , 20   , 0.3 , 1.3  , 250,0,2000);
   Charge_Vs_PathlengthTIDP[elepos]  = ibooker.book2S(cvpTIDP.c_str() , cvpTIDP.c_str() , 20   , 0.3 , 1.3  , 250,0,2000);
   Charge_Vs_PathlengthTIDM[elepos]  = ibooker.book2S(cvpTIDM.c_str() , cvpTIDM.c_str() , 20   , 0.3 , 1.3  , 250,0,2000);
   Charge_Vs_PathlengthTECP1[elepos] = ibooker.book2S(cvpTECP1.c_str(), cvpTECP1.c_str(), 20   , 0.3 , 1.3  , 250,0,2000);
   Charge_Vs_PathlengthTECP2[elepos] = ibooker.book2S(cvpTECP2.c_str(), cvpTECP2.c_str(), 20   , 0.3 , 1.3  , 250,0,2000);
   Charge_Vs_PathlengthTECM1[elepos] = ibooker.book2S(cvpTECM1.c_str(), cvpTECM1.c_str(), 20   , 0.3 , 1.3  , 250,0,2000);
   Charge_Vs_PathlengthTECM2[elepos] = ibooker.book2S(cvpTECM2.c_str(), cvpTECM2.c_str(), 20   , 0.3 , 1.3  , 250,0,2000);
 
   std::vector<std::pair<std::string,std::string>> hnames = APVGain::monHnames(VChargeHisto,doChargeMonitorPerPlane,"");
   for (unsigned int i=0;i<hnames.size();i++){
     std::string htag = (hnames[i]).first + stag;
     MonitorElement* monitor = ibooker.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 = ibooker.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 = ibooker.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 = ibooker.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) );
   }
 }

void SiStripGainsPCLWorker::checkBookAPVColls ( const edm::EventSetup & setup )

privatevirtual

Definition at line 363 of file SiStripGainsPCLWorker.cc.

References sistrip::APV, APVsColl, APVsCollOrdered, bareTkGeomPtr_, TrackerGeometry::dets(), edm::EventSetup::get(), mps_fire::i, PixelTopology::ncolumns(), NPixelDets, PixelTopology::nrows(), NStripAPVs, StripTopology::nstrips(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::rawId(), ntupleEnum::SubDet, DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, tkGeom_, and StripSubdetector::TOB.

Referenced by dqmBeginRun().

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

void SiStripGainsPCLWorker::dqmBeginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 107 of file SiStripGainsPCLWorker.cc.

References a, sistrip::APV, APVsCollOrdered, checkBookAPVColls(), cmsRelvalreport::exit, edm::EventSetup::get(), SiStripGain::getApvGain(), SiStripGain::getNumberOfTags(), SiStripGain::getRange(), SiStripQuality::IsApvBad(), edm::ESHandleBase::isValid(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, and mps_setup::stdout.

                                                                                 {
   
   using namespace edm;
 
   this->checkBookAPVColls(iSetup); // check whether APV colls are booked and do so if not yet done
   
   edm::ESHandle<SiStripGain> gainHandle;
   iSetup.get<SiStripGainRcd>().get(gainHandle);
   if(!gainHandle.isValid()){edm::LogError("SiStripGainPCLWorker")<< "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++){
 
     std::shared_ptr<stAPVGain> APV = APVsCollOrdered[a];
 
     if(APV->SubDet==PixelSubdetector::PixelBarrel || APV->SubDet==PixelSubdetector::PixelEndcap) continue;
     
     APV->isMasked      = SiStripQuality_->IsApvBad(APV->DetId,APV->APVId);
           
     if(gainHandle->getNumberOfTags()!=2){edm::LogError("SiStripGainPCLWorker")<< "NUMBER OF GAIN TAG IS EXPECTED TO BE 2\n";fflush(stdout);exit(0);};          
     float newPreviousGain = gainHandle->getApvGain(APV->APVId,gainHandle->getRange(APV->DetId, 1),1);
     if(APV->PreviousGain!=1 and newPreviousGain!=APV->PreviousGain)edm::LogWarning("SiStripGainPCLWorker")<< "WARNING: ParticleGain in the global tag changed\n";
     APV->PreviousGain = newPreviousGain;
     
     float newPreviousGainTick = gainHandle->getApvGain(APV->APVId,gainHandle->getRange(APV->DetId, 0),0);
     if(APV->PreviousGainTick!=1 and newPreviousGainTick!=APV->PreviousGainTick){
       edm::LogWarning("SiStripGainPCLWorker")<< "WARNING: TickMarkGain in the global tag changed\n"<< std::endl
                          <<" APV->SubDet: "<< APV->SubDet << " APV->APVId:" << APV->APVId << std::endl
                          <<" APV->PreviousGainTick: "<<APV->PreviousGainTick<<" newPreviousGainTick: "<<newPreviousGainTick<<std::endl;
     }
     APV->PreviousGainTick = newPreviousGainTick;      
   }
 
 }

void SiStripGainsPCLWorker::endJob ( void )

privatevirtual

Definition at line 479 of file SiStripGainsPCLWorker.cc.

Referenced by o2olib.O2ORunMgr::executeJob().

480 {

481 }

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

static

Definition at line 485 of file SiStripGainsPCLWorker.cc.

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

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

void SiStripGainsPCLWorker::processEvent ( const TrackerTopology * topo )

private

Definition at line 190 of file SiStripGainsPCLWorker.cc.

Referenced by analyze().

                                                                     {
 
   edm::LogInfo("SiStripGainsPCLWorker") << "Processing run " << runnumber 
                     << " and event " << eventnumber 
                     << std::endl;
 
   if(runnumber<SRun)SRun=runnumber;
   if(runnumber>ERun)ERun=runnumber;
   
   NEvent++;
   NTrack+=(*trackp).size();
 
   edm::LogInfo("SiStripGainsPCLWorker")
     <<"for mode"<< m_calibrationMode 
     <<" Nevent:"<<NEvent 
     <<" NTrack:"<<NTrack
     <<" NClusterStrip:"<<NClusterStrip
     <<" NClusterPixel:"<<NClusterPixel
     <<" NStripAPVs:"<<NStripAPVs
     <<" NPixelDets:"<<NPixelDets
     <<std::endl;
 
   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];
     
     if((*tracketa      )[TI]  < MinTrackEta          )continue;
     if((*tracketa      )[TI]  > MaxTrackEta          )continue;
     if((*trackp        )[TI]  < MinTrackMomentum     )continue;
     if((*trackp        )[TI]  > MaxTrackMomentum     )continue;
     if((*trackhitsvalid)[TI]  < MinTrackHits         )continue;
     if((*trackchi2ndof )[TI]  > MaxTrackChiOverNdf   )continue;
     if((*trackalgo     )[TI]  > MaxTrackingIteration )continue;
     
     std::shared_ptr<stAPVGain> APV = 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;
     
     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               
     }
     
     double ClusterChargeOverPath   =  ( (double) Charge )/(*path)[i] ;
     if(APV->SubDet>2){
       if(Validation)     {ClusterChargeOverPath/=(*gainused)[i];}
       if(OldGainRemoving){ClusterChargeOverPath*=(*gainused)[i];}
     }
 
     // keep processing of pixel cluster charge until here
     if(APV->SubDet<=2) continue;
       
     // real histogram for calibration
     (Charge_Vs_Index[elepos])->Fill(APV->Index,ClusterChargeOverPath);
     
     LogDebug("SiStripGainsPCLWorker") <<" for mode "<< m_calibrationMode << "\n"
                       <<" i "<< i
                       <<" NClusterStrip "<< NClusterStrip
                       <<" useCalibration "<< useCalibration
                       <<" FirstSetOfConstants "<< FirstSetOfConstants
                       <<" APV->DetId "<< APV->DetId
                       <<" APV->Index "<< APV->Index 
                       <<" Charge "<< Charge
                       <<" ClusterChargeOverPath "<< ClusterChargeOverPath
                       <<std::endl;
     
     // Fill monitoring histograms
     int mCharge1 = 0;
     int mCharge2 = 0;
     int mCharge3 = 0;
     int mCharge4 = 0;
     if(APV->SubDet>2) {
       for(unsigned int s=0;s<(*nstrips)[i];s++){
         int StripCharge =  (*amplitude)[FirstAmplitude-(*nstrips)[i]+s];
         if(StripCharge>1024)      StripCharge = 255;
         else if(StripCharge>254)  StripCharge = 254;
         mCharge1 += StripCharge;
         mCharge2 += StripCharge;
         mCharge3 += StripCharge;
         mCharge4 += StripCharge;
       }
       // Revome gains for monitoring
       mCharge2 *= (*gainused)[i];                         // remove G2
       mCharge3 *= (*gainusedTick)[i];                     // remove G1
       mCharge4 *= ( (*gainused)[i] * (*gainusedTick)[i]); // remove G1 and G2
     } 
     std::vector<APVGain::APVmon>& v1 = Charge_1[elepos];
     std::vector<MonitorElement*> cmon1 = APVGain::FetchMonitor(v1, (*rawid)[i], topo);
     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], topo);
     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], topo);
     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], topo);
     for(unsigned int m=0; m<cmon4.size(); m++) cmon4[m]->Fill(( (double) mCharge4 )/(*path)[i]);
 
 
     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
 
   LogDebug("SiStripGainsPCLWorker")<<" for mode"<< m_calibrationMode 
                    <<" entries in histogram:"<< (Charge_Vs_Index[elepos])->getTH2S()->GetEntries()
                    <<std::endl;
 
 }//END OF processEvent()

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

inlineprivate

Definition at line 190 of file SiStripGainsPCLWorker.h.

References dqm_tag_, None, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by bookHistograms(), and processEvent().

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

Member Data Documentation

bool SiStripGainsPCLWorker::AllowSaturation

private

Definition at line 129 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

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

private

Definition at line 173 of file SiStripGainsPCLWorker.h.

Referenced by analyze().

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

private

Definition at line 173 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

std::unordered_map<unsigned int, std::shared_ptr<stAPVGain> > SiStripGainsPCLWorker::APVsColl

private

Definition at line 185 of file SiStripGainsPCLWorker.h.

Referenced by checkBookAPVColls(), and processEvent().

std::vector<std::shared_ptr<stAPVGain> > SiStripGainsPCLWorker::APVsCollOrdered

private

Definition at line 184 of file SiStripGainsPCLWorker.h.

Referenced by checkBookAPVColls(), and dqmBeginRun().

const TrackerGeometry* SiStripGainsPCLWorker::bareTkGeomPtr_

private

Definition at line 141 of file SiStripGainsPCLWorker.h.

Referenced by checkBookAPVColls().

std::string SiStripGainsPCLWorker::CalibPrefix_

private

Definition at line 181 of file SiStripGainsPCLWorker.h.

Referenced by SiStripGainsPCLWorker().

std::string SiStripGainsPCLWorker::CalibSuffix_

private

Definition at line 182 of file SiStripGainsPCLWorker.h.

Referenced by SiStripGainsPCLWorker().

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

private

Definition at line 170 of file SiStripGainsPCLWorker.h.

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

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

private

Charge per cm per layer / wheel

Definition at line 96 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), and processEvent().

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

private

Charge per cm per layer / wheel without G2

Definition at line 97 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), and processEvent().

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

private

Charge per cm per layer / wheel without G1

Definition at line 98 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), and processEvent().

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

private

Charge per cm per layer / wheel without G1 and G1

Definition at line 99 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), and processEvent().

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

private

Definition at line 170 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Charge per cm for each detector id

Definition at line 95 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), processEvent(), and SiStripGainsPCLWorker().

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

private

Charge vs pathlength in TECP thin

Definition at line 107 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), processEvent(), and SiStripGainsPCLWorker().

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

private

Charge vs pathlength in TECP thick

Definition at line 108 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), processEvent(), and SiStripGainsPCLWorker().

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

private

Charge vs pathlength in TECP thin

Definition at line 105 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), processEvent(), and SiStripGainsPCLWorker().

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

private

Charge vs pathlength in TECP thick

Definition at line 106 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), processEvent(), and SiStripGainsPCLWorker().

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

private

Charge vs pathlength in TIB

Definition at line 101 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), processEvent(), and SiStripGainsPCLWorker().

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

private

Charge vs pathlength in TIDM

Definition at line 104 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), processEvent(), and SiStripGainsPCLWorker().

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

private

Charge vs pathlength in TIDP

Definition at line 103 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), processEvent(), and SiStripGainsPCLWorker().

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

private

Charge vs pathlength in TOB

Definition at line 102 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), processEvent(), and SiStripGainsPCLWorker().

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

private

Definition at line 172 of file SiStripGainsPCLWorker.h.

Referenced by analyze().

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

private

Definition at line 172 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

bool SiStripGainsPCLWorker::doChargeMonitorPerPlane

private

Charge monitor per detector plane

Definition at line 134 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), and SiStripGainsPCLWorker().

std::vector<std::string> SiStripGainsPCLWorker::dqm_tag_

private

Definition at line 91 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), SiStripGainsPCLWorker(), and statCollectionFromMode().

unsigned int SiStripGainsPCLWorker::ERun

private

Definition at line 119 of file SiStripGainsPCLWorker.h.

Referenced by processEvent().

unsigned int SiStripGainsPCLWorker::eventnumber =0

private

Definition at line 146 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

std::string SiStripGainsPCLWorker::EventPrefix_

private

Definition at line 177 of file SiStripGainsPCLWorker.h.

Referenced by SiStripGainsPCLWorker().

std::string SiStripGainsPCLWorker::EventSuffix_

private

Definition at line 178 of file SiStripGainsPCLWorker.h.

Referenced by SiStripGainsPCLWorker().

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

private

Definition at line 169 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 169 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

bool SiStripGainsPCLWorker::FirstSetOfConstants

private

Definition at line 130 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

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

private

Definition at line 165 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 165 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 174 of file SiStripGainsPCLWorker.h.

Referenced by analyze().

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

private

Definition at line 174 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 175 of file SiStripGainsPCLWorker.h.

Referenced by analyze().

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

private

Definition at line 175 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 162 of file SiStripGainsPCLWorker.h.

Referenced by analyze().

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

private

Definition at line 162 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 163 of file SiStripGainsPCLWorker.h.

Referenced by analyze().

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

private

Definition at line 163 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 164 of file SiStripGainsPCLWorker.h.

Referenced by analyze().

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

private

Definition at line 164 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

std::string SiStripGainsPCLWorker::m_calibrationMode

private

Type of statistics for the calibration

Definition at line 137 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), processEvent(), and SiStripGainsPCLWorker().

std::string SiStripGainsPCLWorker::m_DQMdir

private

DQM folder hosting the charge statistics and the monitor plots

Definition at line 136 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), and SiStripGainsPCLWorker().

unsigned int SiStripGainsPCLWorker::MaxNrStrips

private

Definition at line 125 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

double SiStripGainsPCLWorker::MaxTrackChiOverNdf

private

Definition at line 127 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

double SiStripGainsPCLWorker::MaxTrackEta

private

Definition at line 124 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

int SiStripGainsPCLWorker::MaxTrackingIteration

private

Definition at line 128 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

double SiStripGainsPCLWorker::MaxTrackMomentum

private

Definition at line 122 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

double SiStripGainsPCLWorker::MinTrackEta

private

Definition at line 123 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

unsigned int SiStripGainsPCLWorker::MinTrackHits

private

Definition at line 126 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

double SiStripGainsPCLWorker::MinTrackMomentum

private

Definition at line 121 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

unsigned int SiStripGainsPCLWorker::NClusterPixel

private

Definition at line 115 of file SiStripGainsPCLWorker.h.

Referenced by processEvent().

unsigned int SiStripGainsPCLWorker::NClusterStrip

private

Definition at line 114 of file SiStripGainsPCLWorker.h.

Referenced by processEvent().

unsigned int SiStripGainsPCLWorker::NEvent

private

Definition at line 112 of file SiStripGainsPCLWorker.h.

Referenced by processEvent().

int SiStripGainsPCLWorker::NPixelDets

private

Definition at line 117 of file SiStripGainsPCLWorker.h.

Referenced by checkBookAPVColls(), and processEvent().

int SiStripGainsPCLWorker::NStripAPVs

private

Definition at line 116 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), checkBookAPVColls(), and processEvent().

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

private

Definition at line 166 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 166 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

unsigned int SiStripGainsPCLWorker::NTrack

private

Definition at line 113 of file SiStripGainsPCLWorker.h.

Referenced by processEvent().

bool SiStripGainsPCLWorker::OldGainRemoving

private

Definition at line 132 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

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

private

Definition at line 168 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 168 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 171 of file SiStripGainsPCLWorker.h.

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

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

private

Definition at line 171 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 161 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 161 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

unsigned int SiStripGainsPCLWorker::runnumber =0

private

Definition at line 147 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 167 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 167 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

unsigned int SiStripGainsPCLWorker::SRun

private

Definition at line 118 of file SiStripGainsPCLWorker.h.

Referenced by processEvent().

edm::ESHandle<TrackerGeometry> SiStripGainsPCLWorker::tkGeom_

private

Definition at line 140 of file SiStripGainsPCLWorker.h.

Referenced by checkBookAPVColls().

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

private

Definition at line 157 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 157 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 151 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 151 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 154 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 154 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 156 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 156 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 160 of file SiStripGainsPCLWorker.h.

Referenced by analyze().

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

private

Definition at line 160 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 152 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and processEvent().

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

private

Definition at line 152 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

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

private

Definition at line 155 of file SiStripGainsPCLWorker.h.

Referenced by analyze().

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

private

Definition at line 155 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

std::string SiStripGainsPCLWorker::TrackPrefix_

private

Definition at line 179 of file SiStripGainsPCLWorker.h.

Referenced by SiStripGainsPCLWorker().

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

private

Definition at line 153 of file SiStripGainsPCLWorker.h.

Referenced by analyze().

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

private

Definition at line 153 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

std::string SiStripGainsPCLWorker::TrackSuffix_

private

Definition at line 180 of file SiStripGainsPCLWorker.h.

Referenced by SiStripGainsPCLWorker().

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

private

Definition at line 148 of file SiStripGainsPCLWorker.h.

Referenced by analyze().

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

private

Definition at line 148 of file SiStripGainsPCLWorker.h.

Referenced by analyze(), and SiStripGainsPCLWorker().

bool SiStripGainsPCLWorker::useCalibration

private

Definition at line 133 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

bool SiStripGainsPCLWorker::Validation

private

Definition at line 131 of file SiStripGainsPCLWorker.h.

Referenced by processEvent(), and SiStripGainsPCLWorker().

std::vector<std::string> SiStripGainsPCLWorker::VChargeHisto

private

Charge monitor plots to be output

Definition at line 138 of file SiStripGainsPCLWorker.h.

Referenced by bookHistograms(), and SiStripGainsPCLWorker().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation