Inheritance diagram for SiStripGainFromCalibTree:

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
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from 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 93 of file SiStripGainFromCalibTree.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 340 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 1647 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 1682 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 1221 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 631 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, 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 791 of file SiStripGainFromCalibTree.cc.

References a, AlgoMode, sistrip::APV, stAPVGain::APVId, APVsCollOrdered, stAPVGain::DetId, cmsRelvalreport::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 1274 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 991 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 839 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 457 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob().

 {
     edm::LogInfo("SiStripGainFromCalibTree") << "Setting " << dqm_dir << "in DQM and booking histograms for tag "
                                              << tag << std::endl;
 
     if ( strcmp(booked_dir_.c_str(),dqm_dir)!=0 ) {
         booked_dir_ = dqm_dir;
         dbe->setCurrentFolder(dqm_dir);
     }
 
     std::string stag(tag);
     if(!stag.empty() && stag[0]!='_') stag.insert(0,1,'_');
 
     std::string cvi      = std::string("Charge_Vs_Index") + stag;
     //std::string cviA     = std::string("Charge_Vs_Index_Absolute")  + stag;
     std::string cvpTIB   = std::string("Charge_Vs_PathlengthTIB")   + stag;
     std::string cvpTOB   = std::string("Charge_Vs_PathlengthTOB")   + stag;
     std::string cvpTIDP  = std::string("Charge_Vs_PathlengthTIDP")  + stag;
     std::string cvpTIDM  = std::string("Charge_Vs_PathlengthTIDM")  + stag;
     std::string cvpTECP1 = std::string("Charge_Vs_PathlengthTECP1") + stag;
     std::string cvpTECP2 = std::string("Charge_Vs_PathlengthTECP2") + stag;
     std::string cvpTECM1 = std::string("Charge_Vs_PathlengthTECM1") + stag;
     std::string cvpTECM2 = std::string("Charge_Vs_PathlengthTECM2") + stag;
 
     int elepos = (m_harvestingMode && AlgoMode=="PCL")? Harvest : statCollectionFromMode(tag);
 
     // The cluster charge is stored by exploiting a non uniform binning in order 
     // reduce the histogram memory size. The bin width is relaxed with a falling
     // exponential function and the bin boundaries are stored in the binYarray. 
     // The binXarray is used to provide as many bins as the APVs.
     //
     // More details about this implementations are here:
     // https://indico.cern.ch/event/649344/contributions/2672267/attachments/1498323/2332518/OptimizeChHisto.pdf
 
     std::vector<float> binXarray;
     binXarray.reserve( NStripAPVs+1 );
     for(int a=0;a<=NStripAPVs;a++){
        binXarray.push_back( (float)a );
     }
  
     std::array<float,688> binYarray;
     double p0 = 5.445;
     double p1 = 0.002113;
     double p2 = 69.01576;
     double y = 0.;
     for(int b=0;b<687;b++) {
        binYarray[b] = y;
        if(y<=902.) y = y + 2.;
        else y = ( p0 - log(exp(p0-p1*y) - p2*p1)) / p1;
     }
     binYarray[687] = 4000.;
 
     Charge_Vs_Index[elepos]           = dbe->book2S(cvi.c_str()     , cvi.c_str()     , NStripAPVs, &binXarray[0], 687, binYarray.data());
     //Charge_Vs_Index_Absolute[elepos]  = dbe->book2S(cviA.c_str()    , cviA.c_str()    , 88625, 0   , 88624,1000,0,4000);
     Charge_Vs_PathlengthTIB[elepos]   = dbe->book2S(cvpTIB.c_str()  , cvpTIB.c_str()  , 20   , 0.3 , 1.3  , 250,0,2000);
     Charge_Vs_PathlengthTOB[elepos]   = dbe->book2S(cvpTOB.c_str()  , cvpTOB.c_str()  , 20   , 0.3 , 1.3  , 250,0,2000);
     Charge_Vs_PathlengthTIDP[elepos]  = dbe->book2S(cvpTIDP.c_str() , cvpTIDP.c_str() , 20   , 0.3 , 1.3  , 250,0,2000);
     Charge_Vs_PathlengthTIDM[elepos]  = dbe->book2S(cvpTIDM.c_str() , cvpTIDM.c_str() , 20   , 0.3 , 1.3  , 250,0,2000);
     Charge_Vs_PathlengthTECP1[elepos] = dbe->book2S(cvpTECP1.c_str(), cvpTECP1.c_str(), 20   , 0.3 , 1.3  , 250,0,2000);
     Charge_Vs_PathlengthTECP2[elepos] = dbe->book2S(cvpTECP2.c_str(), cvpTECP2.c_str(), 20   , 0.3 , 1.3  , 250,0,2000);
     Charge_Vs_PathlengthTECM1[elepos] = dbe->book2S(cvpTECM1.c_str(), cvpTECM1.c_str(), 20   , 0.3 , 1.3  , 250,0,2000);
     Charge_Vs_PathlengthTECM2[elepos] = dbe->book2S(cvpTECM2.c_str(), cvpTECM2.c_str(), 20   , 0.3 , 1.3  , 250,0,2000);
 
     //Book Charge monitoring histograms
     std::vector<std::pair<std::string,std::string>> hnames = APVGain::monHnames(VChargeHisto,doChargeMonitorPerPlane,"");
     for (unsigned int i=0;i<hnames.size();i++){
         std::string htag = (hnames[i]).first + stag;
         MonitorElement* monitor = dbe->book1DD( htag.c_str(), (hnames[i]).second.c_str(), 100   , 0. , 1000. );
         int 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 127 of file SiStripGainFromCalibTree.cc.

References edm::Event::getByToken(), cmsBatch::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 1608 of file SiStripGainFromCalibTree.cc.

References a, sistrip::APV, APVsCollOrdered, stAPVGain::DetId, stAPVGain::Gain, m_harvestingMode, haddnano::obj, produceTagFilter(), 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 1038 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 772 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 1065 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 1658 of file SiStripGainFromCalibTree.cc.

References sistrip::APV, APVsColl, stAPVGain::CalibGain, m_calibrationPath, 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 327 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 1072 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 1578 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 1352 of file SiStripGainFromCalibTree.cc.

References a, AlgoMode, sistrip::APV, APVsCollOrdered, stringResolutionProvider_cfi::bin, Charge_Vs_Index, stAPVGain::DetId, DiffWRTPrevGainTEC, DiffWRTPrevGainTIB, DiffWRTPrevGainTID, DiffWRTPrevGainTOB, MillePedeFileConverter_cfg::e, stAPVGain::Eta, reco::tau::disc::Eta(), APVGain::FetchMonitor(), MonitorElement::Fill(), stAPVGain::FitMPV, stAPVGain::FitMPVErr, stAPVGain::Gain, g4SimHits_cfi::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, colinearityKinematic::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 314 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 1462 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 783 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 169 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::AllowSaturation

private

Definition at line 150 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

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

private

Definition at line 294 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 294 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 309 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 308 of file SiStripGainFromCalibTree.cc.

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

unsigned int SiStripGainFromCalibTree::BAD

private

Definition at line 261 of file SiStripGainFromCalibTree.cc.

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

std::string SiStripGainFromCalibTree::booked_dir_

private

Definition at line 176 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos().

string SiStripGainFromCalibTree::CalibPrefix_

private

Definition at line 303 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

int SiStripGainFromCalibTree::CalibrationLevel

private

Definition at line 154 of file SiStripGainFromCalibTree.cc.

Referenced by algoComputeMPVandGain(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::CalibSuffix_

private

Definition at line 304 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 291 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 179 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 180 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 181 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 182 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 291 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 178 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 190 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 191 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 188 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 189 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 184 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 187 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 186 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 185 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 293 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 293 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

DQMStore* SiStripGainFromCalibTree::dbe

private

Definition at line 137 of file SiStripGainFromCalibTree.cc.

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

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTEC

private

ratio gain / PreviousGain for TEC disks

Definition at line 242 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTIB

private

ratio Gain / PreviousGain for TIB layers

Definition at line 239 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTID

private

ratio Gain / PreviousGain for TID disks

Definition at line 240 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::DiffWRTPrevGainTOB

private

ratio Gain / PreviousGain for TOB layers

Definition at line 241 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

bool SiStripGainFromCalibTree::doChargeMonitorPerPlane

private

Charge monitor per detector plane

Definition at line 160 of file SiStripGainFromCalibTree.cc.

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

std::vector<string> SiStripGainFromCalibTree::dqm_tag_

private

Definition at line 175 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::ERun

private

Definition at line 259 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and processEvent().

unsigned int SiStripGainFromCalibTree::eventnumber =0

private

Definition at line 267 of file SiStripGainFromCalibTree.cc.

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

string SiStripGainFromCalibTree::EventPrefix_

private

Definition at line 299 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::EventSuffix_

private

Definition at line 300 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 290 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 290 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

bool SiStripGainFromCalibTree::FirstSetOfConstants

private

Definition at line 151 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

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

private

Definition at line 286 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 286 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

MonitorElement* SiStripGainFromCalibTree::Gains

private

distribution of gain factors

Definition at line 214 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 295 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 295 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 296 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 296 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTEC

private

Gain vs PreviousGain for TEC

Definition at line 247 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTIB

private

Gain vs PreviousGain for TIB

Definition at line 244 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTID

private

Gain vs PreviousGain for TID

Definition at line 245 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::GainVsPrevGainTOB

private

Gain vs PreviousGain for TOB

Definition at line 246 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

unsigned int SiStripGainFromCalibTree::GOOD

private

Definition at line 260 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 283 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 283 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 284 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 284 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 285 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 285 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

std::string SiStripGainFromCalibTree::m_calibrationMode

private

Type of statistics for the calibration

Definition at line 161 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 162 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 163 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::m_harvestingMode

private

Definition at line 158 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::m_splitDQMstat

private

Definition at line 159 of file SiStripGainFromCalibTree.cc.

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

double SiStripGainFromCalibTree::MagFieldCurrentTh

private

Definition at line 138 of file SiStripGainFromCalibTree.cc.

Referenced by isBFieldConsistentWithMode(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::MASKED

private

Definition at line 262 of file SiStripGainFromCalibTree.cc.

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

double SiStripGainFromCalibTree::MaxChi2OverNDF

private

Definition at line 141 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxMPVError

private

Definition at line 140 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::MaxNrStrips

private

Definition at line 146 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxTrackChiOverNdf

private

Definition at line 148 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxTrackEta

private

Definition at line 145 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

int SiStripGainFromCalibTree::MaxTrackingIteration

private

Definition at line 149 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxTrackMomentum

private

Definition at line 143 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MinNrEntries

private

Definition at line 139 of file SiStripGainFromCalibTree.cc.

Referenced by getPeakOfLandau(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MinTrackEta

private

Definition at line 144 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::MinTrackHits

private

Definition at line 147 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MinTrackMomentum

private

Definition at line 142 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTEC

private

MPV vs Eta for TEC planes

Definition at line 200 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 202 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 201 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTIB

private

MPV vs Eta for TIB planes

Definition at line 197 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTID

private

MPV vs Eta for TID planes

Definition at line 198 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_EtaTOB

private

MPV vs Eta for TOB planes

Definition at line 199 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTEC

private

MPV vs Phi for TEC planes

Definition at line 207 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 209 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 208 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTIB

private

MPV vs Phi for TIB planes

Definition at line 204 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTID

private

MPV vs Phi for TID planes

Definition at line 205 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPV_Vs_PhiTOB

private

MPV vs Phi for TOB planes

Definition at line 206 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVError

private

error of Landau fit

Definition at line 233 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsEta

private

error of Landau fit vs Eta

Definition at line 235 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsMPV

private

error of Landau fit vs MPV

Definition at line 234 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsN

private

error of landau fit vs number of entries

Definition at line 237 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVErrorVsPhi

private

error of Landau fit vs Phi

Definition at line 236 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVs

private

distribution of MPVs

Definition at line 215 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVs320

private

distribution of MPVs for thin sensors

Definition at line 216 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVs500

private

distribution of MPVs for tick sensors

Definition at line 217 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTEC

private

distribution of MPVs for TEC disks

Definition at line 223 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECM

private

distribution of MPVs for TECM disks

Definition at line 225 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECM1

private

distribution of MPVs for TECM thin sensors

Definition at line 230 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECM2

private

distribution of MPVs for TECM tick sensors

Definition at line 231 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECP

private

distribution of MPVs for TECP disks

Definition at line 224 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECP1

private

distribution of MPVs for TECP thin sensors

Definition at line 228 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECP2

private

distribution of MPVs for TECP tick sensors

Definition at line 229 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECthick

private

distribution of MPVs for TEC tick sensors

Definition at line 227 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTECthin

private

distribution of MPVs for TEC thin sensors

Definition at line 226 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTIB

private

distribution of MPVs for TIB planes

Definition at line 218 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTID

private

distribution of MPVs for TID disks

Definition at line 219 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTIDM

private

distribution of MPVs for TIDM disks

Definition at line 221 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTIDP

private

distribution of MPVs for TIDP disks

Definition at line 220 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

MonitorElement* SiStripGainFromCalibTree::MPVsTOB

private

distribution of MPVs for TOB planes

Definition at line 222 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

unsigned int SiStripGainFromCalibTree::NClusterPixel

private

Definition at line 255 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and processEvent().

unsigned int SiStripGainFromCalibTree::NClusterStrip

private

Definition at line 254 of file SiStripGainFromCalibTree.cc.

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

unsigned int SiStripGainFromCalibTree::NEvent

private

Definition at line 252 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 249 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 212 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 211 of file SiStripGainFromCalibTree.cc.

Referenced by bookDQMHistos(), and qualityMonitor().

int SiStripGainFromCalibTree::NPixelDets

private

Definition at line 257 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob().

int SiStripGainFromCalibTree::NStripAPVs

private

Definition at line 256 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 287 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 287 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::NTrack

private

Definition at line 253 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::OldGainRemoving

private

Definition at line 153 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

std::string SiStripGainFromCalibTree::OutputGains

private

Definition at line 170 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree(), and storeOnTree().

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

private

Definition at line 289 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 289 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 292 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 292 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 282 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 282 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::runnumber =0

private

Definition at line 268 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 288 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 288 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

bool SiStripGainFromCalibTree::saveSummary

private

Definition at line 156 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::SRun

private

Definition at line 258 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and processEvent().

double SiStripGainFromCalibTree::tagCondition_GoodFrac

private

Definition at line 167 of file SiStripGainFromCalibTree.cc.

Referenced by produceTagFilter(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::tagCondition_NClusters

private

Definition at line 166 of file SiStripGainFromCalibTree.cc.

Referenced by produceTagFilter(), and SiStripGainFromCalibTree().

TFileService* SiStripGainFromCalibTree::tfs

private

Definition at line 136 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob().

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

private

Definition at line 278 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 278 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 272 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 272 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 275 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 275 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 277 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 281 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 281 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 273 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 273 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 276 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 301 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 274 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 274 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::TrackSuffix_

private

Definition at line 302 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 269 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 269 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

const TrackerTopology* SiStripGainFromCalibTree::tTopo_

private

Definition at line 310 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::useCalibration

private

Definition at line 157 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::Validation

private

Definition at line 152 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 164 of file SiStripGainFromCalibTree.cc.

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

vector<string> SiStripGainFromCalibTree::VInputFiles

private

Definition at line 171 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation