SiStripGainFromCalibTree Class Reference

Inheritance diagram for SiStripGainFromCalibTree:

Classes
class	isEqual

Public Member Functions
	SiStripGainFromCalibTree (const edm::ParameterSet &)
	~SiStripGainFromCalibTree ()
Public Member Functions inherited from ConditionDBWriter< SiStripApvGain >
	ConditionDBWriter (const edm::ParameterSet &iConfig)
virtual	~ConditionDBWriter ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Types
enum	statistic_type {   None =-1, StdBunch, StdBunch0T, FaABunch,   FaABunch0T, IsoBunch, IsoBunch0T, Harvest }

Private Member Functions
virtual void	algoAnalyze (const edm::Event &, const edm::EventSetup &) override
void	algoAnalyzeTheTree ()
virtual void	algoBeginJob (const edm::EventSetup &iSetup) override
virtual void	algoBeginRun (const edm::Run &run, const edm::EventSetup &iSetup) override
void	algoComputeMPVandGain ()
virtual void	algoEndJob () override
virtual 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)
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 ()
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_
__gnu_cxx::hash_map< unsigned int, stAPVGain *, __gnu_cxx::hash< unsigned int >, isEqual >	APVsColl
std::vector< stAPVGain * >	APVsCollOrdered
unsigned int	BAD
string	CalibPrefix_
int	CalibrationLevel
string	CalibSuffix_
const std::vector< unsigned int > *	charge =0
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
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_
const std::vector< double > *	gainused =0
edm::EDGetTokenT< std::vector < double > >	gainused_token_
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_
string	m_calibrationMode
string	m_calibrationPath
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
unsigned int	NClusterPixel
unsigned int	NClusterStrip
unsigned int	NEvent
int	NPixelDets
int	NStripAPVs
const std::vector< unsigned short > *	nstrips =0
edm::EDGetTokenT< std::vector < unsigned short > >	nstrips_token_
unsigned int	NTrack
bool	OldGainRemoving
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_
bool	useCalibration
bool	Validation
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 &)
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 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 119 of file SiStripGainFromCalibTree.cc.

Member Enumeration Documentation

enum SiStripGainFromCalibTree::statistic_type

private

Enumerator
None
StdBunch
StdBunch0T
FaABunch
FaABunch0T
IsoBunch
IsoBunch0T
Harvest

Definition at line 196 of file SiStripGainFromCalibTree.cc.

{None=-1, StdBunch, StdBunch0T, FaABunch, FaABunch0T, IsoBunch, IsoBunch0T, Harvest};

Constructor & Destructor Documentation

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

explicit

Definition at line 300 of file SiStripGainFromCalibTree.cc.

References AlgoMode, AllowSaturation, amplitude_token_, CalibPrefix_, CalibrationLevel, CalibSuffix_, charge_token_, Charge_Vs_Index, Charge_Vs_PathlengthTECM1, Charge_Vs_PathlengthTECM2, Charge_Vs_PathlengthTECP1, Charge_Vs_PathlengthTECP2, Charge_Vs_PathlengthTIB, Charge_Vs_PathlengthTIDM, Charge_Vs_PathlengthTIDP, Charge_Vs_PathlengthTOB, chargeoverpath_token_, dbe, dqm_tag_, EventPrefix_, EventSuffix_, farfromedge_token_, FirstSetOfConstants, firststrip_token_, gainused_token_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_25ns10e33_v2_cff::InputTag, diffTwoXMLs::label, localdirx_token_, localdiry_token_, localdirz_token_, m_calibrationMode, m_calibrationPath, m_DQMdir, m_harvestingMode, m_splitDQMstat, MagFieldCurrentTh, MaxChi2OverNDF, MaxMPVError, MaxNrStrips, MaxTrackChiOverNdf, MaxTrackEta, MaxTrackingIteration, MaxTrackMomentum, MinNrEntries, MinTrackEta, MinTrackHits, MinTrackMomentum, nstrips_token_, OldGainRemoving, cppFunctionSkipper::operator, OutputGains, overlapping_token_, path_token_, rawid_token_, saturation_token_, saveSummary, AlCaHLTBitMon_QueryRunRegistry::string, tagCondition_GoodFrac, tagCondition_NClusters, trackalgo_token_, trackchi2ndof_token_, tracketa_token_, trackhitsvalid_token_, trackindex_token_, trackp_token_, trackphi_token_, TrackPrefix_, trackpt_token_, TrackSuffix_, TrigTech_token_, useCalibration, Validation, and VInputFiles.

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


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


        // 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(), 0);
        //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(), 0);
        Charge_Vs_PathlengthTOB.insert( Charge_Vs_PathlengthTOB.begin(), dqm_tag_.size(), 0);
        Charge_Vs_PathlengthTIDP.insert( Charge_Vs_PathlengthTIDP.begin(), dqm_tag_.size(), 0);
        Charge_Vs_PathlengthTIDM.insert( Charge_Vs_PathlengthTIDM.begin(), dqm_tag_.size(), 0);
        Charge_Vs_PathlengthTECP1.insert( Charge_Vs_PathlengthTECP1.begin(), dqm_tag_.size(), 0);
        Charge_Vs_PathlengthTECP2.insert( Charge_Vs_PathlengthTECP2.begin(), dqm_tag_.size(), 0);
        Charge_Vs_PathlengthTECM1.insert( Charge_Vs_PathlengthTECM1.begin(), dqm_tag_.size(), 0);
        Charge_Vs_PathlengthTECM2.insert( Charge_Vs_PathlengthTECM2.begin(), dqm_tag_.size(), 0);

        


        // 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_)); 

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

SiStripGainFromCalibTree::~SiStripGainFromCalibTree

(

)

Definition at line 1248 of file SiStripGainFromCalibTree.cc.

{ 
}

Member Function Documentation

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

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 1274 of file SiStripGainFromCalibTree.cc.

References AlgoMode, amplitude, amplitude_token_, charge, charge_token_, chargeoverpath, chargeoverpath_token_, connect(), edm::EventID::event(), eventnumber, farfromedge, farfromedge_token_, firststrip, firststrip_token_, gainused, gainused_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(trackalgo     , trackalgo_token_     , iEvent);
     
    processEvent();
}

void SiStripGainFromCalibTree::algoAnalyzeTheTree ( )

private

Definition at line 933 of file SiStripGainFromCalibTree.cc.

References amplitude, CalibPrefix_, CalibSuffix_, charge, chargeoverpath, eventnumber, EventPrefix_, EventSuffix_, farfromedge, firststrip, gainused, i, localdirx, localdiry, localdirz, m_calibrationMode, NEvent, nstrips, NULL, overlapping, path, processEvent(), rawid, runnumber, saturation, trackalgo, trackchi2ndof, tracketa, trackhitsvalid, trackindex, trackp, trackphi, TrackPrefix_, trackpt, TrackSuffix_, MainPageGenerator::tree, 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   , NULL);
        tree->SetBranchAddress((EventPrefix_ + "run"            + EventSuffix_).c_str(), &runnumber     , NULL);
        tree->SetBranchAddress((EventPrefix_ + "TrigTech"       + EventSuffix_).c_str(), &TrigTech      , NULL);

        tree->SetBranchAddress((TrackPrefix_ + "chi2ndof"       + TrackSuffix_).c_str(), &trackchi2ndof , NULL);
        tree->SetBranchAddress((TrackPrefix_ + "momentum"       + TrackSuffix_).c_str(), &trackp        , NULL);
        tree->SetBranchAddress((TrackPrefix_ + "pt"             + TrackSuffix_).c_str(), &trackpt       , NULL);
        tree->SetBranchAddress((TrackPrefix_ + "eta"            + TrackSuffix_).c_str(), &tracketa      , NULL);
        tree->SetBranchAddress((TrackPrefix_ + "phi"            + TrackSuffix_).c_str(), &trackphi      , NULL);
        tree->SetBranchAddress((TrackPrefix_ + "hitsvalid"      + TrackSuffix_).c_str(), &trackhitsvalid, NULL);
                tree->SetBranchAddress((TrackPrefix_ + "algo"           + TrackSuffix_).c_str(), &trackalgo     , NULL);

        tree->SetBranchAddress((CalibPrefix_ + "trackindex"     + CalibSuffix_).c_str(), &trackindex    , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "rawid"          + CalibSuffix_).c_str(), &rawid         , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "localdirx"      + CalibSuffix_).c_str(), &localdirx     , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "localdiry"      + CalibSuffix_).c_str(), &localdiry     , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "localdirz"      + CalibSuffix_).c_str(), &localdirz     , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "firststrip"     + CalibSuffix_).c_str(), &firststrip    , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "nstrips"        + CalibSuffix_).c_str(), &nstrips       , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "saturation"     + CalibSuffix_).c_str(), &saturation    , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "overlapping"    + CalibSuffix_).c_str(), &overlapping   , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "farfromedge"    + CalibSuffix_).c_str(), &farfromedge   , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "charge"         + CalibSuffix_).c_str(), &charge        , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "path"           + CalibSuffix_).c_str(), &path          , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "chargeoverpath" + CalibSuffix_).c_str(), &chargeoverpath, NULL);
        tree->SetBranchAddress((CalibPrefix_ + "amplitude"      + CalibSuffix_).c_str(), &amplitude     , NULL);
        tree->SetBranchAddress((CalibPrefix_ + "gainused"       + CalibSuffix_).c_str(), &gainused      , NULL);


        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 451 of file SiStripGainFromCalibTree.cc.

References AlgoMode, sistrip::APV, stAPVGain::APVId, APVsColl, APVsCollOrdered, BAD, stAPVGain::Bin, bookDQMHistos(), stAPVGain::DetId, 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, i, stAPVGain::Index, isBFieldConsistentWithMode(), stAPVGain::isMasked, j, m_calibrationMode, m_DQMdir, m_harvestingMode, m_splitDQMstat, MakeCalibrationMap(), MASKED, NClusterPixel, NClusterStrip, PixelTopology::ncolumns(), stAPVGain::NEntries, NEvent, NPixelDets, PixelTopology::nrows(), NStripAPVs, StripTopology::nstrips(), NTrack, or, stAPVGain::Phi, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, stAPVGain::PreviousGain, stAPVGain::PreviousGainTick, stAPVGain::R, DetId::rawId(), SRun, statCollectionFromMode(), AlCaHLTBitMon_QueryRunRegistry::string, stAPVGain::SubDet, DetId::subdetId(), swapBFieldMode(), StripSubdetector::TEC, stAPVGain::Thickness, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, 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<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 || 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->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 606 of file SiStripGainFromCalibTree.cc.

References a, AlgoMode, sistrip::APV, stAPVGain::APVId, APVsCollOrdered, stAPVGain::DetId, cmsRelvalreport::exit, edm::EventSetup::get(), isBFieldConsistentWithMode(), stAPVGain::isMasked, edm::ESHandleBase::isValid(), m_calibrationMode, m_harvestingMode, stAPVGain::PreviousGain, stAPVGain::PreviousGainTick, 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];

        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),1);
        if(APV->PreviousGainTick!=1 and newPreviousGainTick!=APV->PreviousGainTick)edm::LogWarning("SiStripGainFromCalibTree")<< "WARNING: TickMarkGain in the global tag changed\n";
        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 985 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, i, stAPVGain::Index, IsGoodLandauFit(), stAPVGain::isMasked, m_calibrationMode, MASKED, stAPVGain::NEntries, NULL, stAPVGain::PreviousGain, statCollectionFromMode(), and stAPVGain::SubDet.

Referenced by algoEndJob().

                                                     {
    unsigned int I=0;
    TH1F* Proj = NULL;
    double FitResults[6];
    double MPVmean = 300;

        int elepos = (AlgoMode == "PCL")? Harvest : statCollectionFromMode(m_calibrationMode.c_str());

        if ( Charge_Vs_Index[elepos]==0 ) {
            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(__gnu_cxx::hash_map<unsigned int, stAPVGain*,  __gnu_cxx::hash<unsigned int>, isEqual >::iterator 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
                __gnu_cxx::hash_map<unsigned int, stAPVGain*,  __gnu_cxx::hash<unsigned int>, isEqual >::iterator tmpit;
                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 755 of file SiStripGainFromCalibTree.cc.

References algoAnalyzeTheTree(), algoComputeMPVandGain(), AlgoMode, Charge_Vs_Index, Charge_Vs_PathlengthTECM1, Charge_Vs_PathlengthTECM2, Charge_Vs_PathlengthTECP1, Charge_Vs_PathlengthTECP2, Charge_Vs_PathlengthTIB, Charge_Vs_PathlengthTIDM, Charge_Vs_PathlengthTIDP, Charge_Vs_PathlengthTOB, Harvest, m_calibrationMode, m_harvestingMode, TFileService::make(), NClusterPixel, NClusterStrip, NPixelDets, NStripAPVs, cppFunctionSkipper::operator, or, saveSummary, statCollectionFromMode(), storeOnTree(), and tfs.

                                     {
    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();
   
    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]!=0 )           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]!=0 )   tfs->make<TH2S> ( *(Charge_Vs_PathlengthTIB[elepos])->getTH2S() );
            if( Charge_Vs_PathlengthTOB[elepos]!=0 )   tfs->make<TH2S> ( *(Charge_Vs_PathlengthTOB[elepos])->getTH2S() );
            if( Charge_Vs_PathlengthTIDP[elepos]!=0 )  tfs->make<TH2S> ( *(Charge_Vs_PathlengthTIDP[elepos])->getTH2S() );
            if( Charge_Vs_PathlengthTIDM[elepos]!=0 )  tfs->make<TH2S> ( *(Charge_Vs_PathlengthTIDM[elepos])->getTH2S() );
            if( Charge_Vs_PathlengthTECP1[elepos]!=0 ) tfs->make<TH2S> ( *(Charge_Vs_PathlengthTECP1[elepos])->getTH2S() );
            if( Charge_Vs_PathlengthTECP2[elepos]!=0 ) tfs->make<TH2S> ( *(Charge_Vs_PathlengthTECP2[elepos])->getTH2S() );
            if( Charge_Vs_PathlengthTECM1[elepos]!=0 ) tfs->make<TH2S> ( *(Charge_Vs_PathlengthTECM1[elepos])->getTH2S() );
            if( Charge_Vs_PathlengthTECM2[elepos]!=0 ) 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 648 of file SiStripGainFromCalibTree.cc.

References AlgoMode, Charge_Vs_Index, Charge_Vs_PathlengthTECM1, Charge_Vs_PathlengthTECM2, Charge_Vs_PathlengthTECP1, Charge_Vs_PathlengthTECP2, Charge_Vs_PathlengthTIB, Charge_Vs_PathlengthTIDM, Charge_Vs_PathlengthTIDP, Charge_Vs_PathlengthTOB, dbe, Harvest, m_calibrationMode, m_DQMdir, m_harvestingMode, m_splitDQMstat, merge(), statCollectionFromMode(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                       {
    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.size()!=0 && 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.c_str());
                //Charge_Vs_Index_Absolute[elepos]  = dbe->get(cviA.c_str());
                Charge_Vs_PathlengthTIB[elepos]   = dbe->get(cvpTIB.c_str());
                Charge_Vs_PathlengthTOB[elepos]   = dbe->get(cvpTOB.c_str());
                Charge_Vs_PathlengthTIDP[elepos]  = dbe->get(cvpTIDP.c_str());
                Charge_Vs_PathlengthTIDM[elepos]  = dbe->get(cvpTIDM.c_str());
                Charge_Vs_PathlengthTECP1[elepos] = dbe->get(cvpTECP1.c_str());
                Charge_Vs_PathlengthTECP2[elepos] = dbe->get(cvpTECP2.c_str());
                Charge_Vs_PathlengthTECM1[elepos] = dbe->get(cvpTECM1.c_str());
                Charge_Vs_PathlengthTECM2[elepos] = dbe->get(cvpTECM2.c_str());

                if (Charge_Vs_Index[elepos]==0) {
                    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]==0) {
                    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]==0) {
                    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]==0) {
                    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]==0) {
                    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]==0) {
                    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]==0) {
                    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]==0) {
                    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]==0) {
                    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());
            }
    }
}

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

private

Definition at line 411 of file SiStripGainFromCalibTree.cc.

References AlgoMode, Charge_Vs_Index, Charge_Vs_PathlengthTECM1, Charge_Vs_PathlengthTECM2, Charge_Vs_PathlengthTECP1, Charge_Vs_PathlengthTECP2, Charge_Vs_PathlengthTIB, Charge_Vs_PathlengthTIDM, Charge_Vs_PathlengthTIDP, Charge_Vs_PathlengthTOB, dbe, Harvest, m_harvestingMode, statCollectionFromMode(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by algoBeginJob().

{
    static std::string booked_dir = "";

    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.size()!=0 && 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);

    Charge_Vs_Index[elepos]           = dbe->book2S(cvi.c_str()     , cvi.c_str()     , 88625, 0   , 88624,2000,0,4000);
    //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);
}

template<typename T >

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

inlineprivate

Definition at line 152 of file SiStripGainFromCalibTree.cc.

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

Referenced by Vispa.Gui.ZoomableScrollArea.ZoomableScrollArea::__init__(), Vispa.Views.PropertyView.BooleanProperty::__init__(), Vispa.Gui.FindDialog.FindDialog::_addScript(), Vispa.Gui.FindDialog.FindDialog::_addStringProperty(), Vispa.Main.Application.Application::_connectSignals(), Vispa.Plugins.ConfigEditor.CodeTableView.CodeTableView::_createItem(), Vispa.Gui.BoxContentDialog.BoxContentDialog::addButton(), Vispa.Gui.ToolBoxContainer.ToolBoxContainer::addWidget(), algoAnalyze(), Vispa.Views.PropertyView.PropertyView::append(), Vispa.Views.PropertyView.PropertyView::appendAddRow(), Vispa.Main.Application.Application::createAction(), Vispa.Views.PropertyView.TextEditWithButtonProperty::createButton(), Vispa.Views.LineDecayView.LineDecayView::createLineDecayContainer(), Vispa.Views.PropertyView.TextEditWithButtonProperty::createLineEdit(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.PropertyView.TextEditWithButtonProperty::createTextEdit(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filterDialog(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Gui.PortWidget.PortWidget::mouseMoveEvent(), Vispa.Views.PropertyView.BooleanProperty::setChecked(), Vispa.Main.SplitterTab.SplitterTab::setController(), Vispa.Plugins.Browser.BrowserTab.BrowserTab::setController(), Vispa.Views.PropertyView.BooleanProperty::setReadOnly(), Vispa.Views.PropertyView.DropDownProperty::setReadOnly(), Vispa.Views.PropertyView.TextEditWithButtonProperty::setReadOnly(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setTab(), Vispa.Views.PropertyView.IntegerProperty::setValue(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView(), Vispa.Plugins.EdmBrowser.EdmBrowserTabController.EdmBrowserTabController::updateViewMenu(), and Vispa.Plugins.ConfigEditor.ConfigEditorTabController.ConfigEditorTabController::updateViewMenu().

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

SiStripApvGain * SiStripGainFromCalibTree::getNewObject ( )

overrideprivatevirtual

Implements ConditionDBWriter< SiStripApvGain >.

Definition at line 1211 of file SiStripGainFromCalibTree.cc.

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

{
    SiStripApvGain* obj = new 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 = NULL;
    unsigned int PreviousDetId = 0; 
    for(unsigned int a=0;a<APVsCollOrdered.size();a++){
        stAPVGain* APV = APVsCollOrdered[a];
        if(APV==NULL){ printf("Bug\n"); continue; }
        if(APV->SubDet<=2)continue;
        if(APV->DetId != PreviousDetId){
            if(theSiStripVector!=NULL){
                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!=NULL){
        SiStripApvGain::Range range(theSiStripVector->begin(),theSiStripVector->end());
        if ( !obj->put(PreviousDetId,range) )  printf("Bug to put detId = %i\n",PreviousDetId);
    }

    return obj;
}

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

private

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

References edm::EventSetup::get(), m_calibrationMode, MagFieldCurrentTh, and edm::ESHandle< class >::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 822 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 1252 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;
    }
}

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

private

Definition at line 287 of file SiStripGainFromCalibTree.cc.

References x().

Referenced by algoEndRun().

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

void SiStripGainFromCalibTree::processEvent ( )

private

Definition at line 829 of file SiStripGainFromCalibTree.cc.

References AllowSaturation, sistrip::APV, APVsColl, stAPVGain::CalibGain, Charge_Vs_Index, Charge_Vs_PathlengthTECM1, Charge_Vs_PathlengthTECM2, Charge_Vs_PathlengthTECP1, Charge_Vs_PathlengthTECP2, Charge_Vs_PathlengthTIB, Charge_Vs_PathlengthTIDM, Charge_Vs_PathlengthTIDP, Charge_Vs_PathlengthTOB, ERun, stAPVGain::Eta, eventnumber, edm::false, farfromedge, HcalObjRepresent::Fill(), FirstSetOfConstants, firststrip, i, stAPVGain::Index, m_calibrationMode, MaxNrStrips, MaxTrackChiOverNdf, MaxTrackEta, MaxTrackingIteration, MaxTrackMomentum, MinTrackEta, MinTrackHits, MinTrackMomentum, NClusterPixel, NClusterStrip, NEvent, nstrips, NTrack, OldGainRemoving, overlapping, path, stAPVGain::PreviousGain, runnumber, alignCSCRings::s, saturation, SRun, statCollectionFromMode(), stAPVGain::SubDet, StripSubdetector::TEC, stAPVGain::Thickness, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, trackalgo, trackchi2ndof, tracketa, trackhitsvalid, trackp, funct::true, useCalibration, and Validation.

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];}
        }
        //(Charge_Vs_Index_Absolute[elepos])->Fill(APV->Index,Charge);   
        (Charge_Vs_Index[elepos])         ->Fill(APV->Index,ClusterChargeOverPath);

        if(APV->SubDet==StripSubdetector::TIB){ (Charge_Vs_PathlengthTIB[elepos])  ->Fill((*path)[i],Charge); 
        }else if(APV->SubDet==StripSubdetector::TOB){ (Charge_Vs_PathlengthTOB[elepos])  ->Fill((*path)[i],Charge);
        }else if(APV->SubDet==StripSubdetector::TID){
            if(APV->Eta<0){           (Charge_Vs_PathlengthTIDM[elepos]) ->Fill((*path)[i],Charge);
            }else if(APV->Eta>0){                      (Charge_Vs_PathlengthTIDP[elepos]) ->Fill((*path)[i],Charge);
            }
        }else if(APV->SubDet==StripSubdetector::TEC){
            if(APV->Eta<0){
                if(APV->Thickness<0.04){          (Charge_Vs_PathlengthTECM1[elepos])->Fill((*path)[i],Charge);
                }else if(APV->Thickness>0.04){          (Charge_Vs_PathlengthTECM2[elepos])->Fill((*path)[i],Charge);
                }
            }else if(APV->Eta>0){
                if(APV->Thickness<0.04){          (Charge_Vs_PathlengthTECP1[elepos])->Fill((*path)[i],Charge);
                }else if(APV->Thickness>0.04){          (Charge_Vs_PathlengthTECP2[elepos])->Fill((*path)[i],Charge);
                }
            }
        }

    }// END OF ON-CLUSTER LOOP
}//END OF processEvent()

bool SiStripGainFromCalibTree::produceTagFilter ( )

private

Definition at line 1181 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]==0 ) {
        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; 
}

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

inlineprivate

Definition at line 274 of file SiStripGainFromCalibTree.cc.

References FSQ::None, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by algoBeginJob(), algoComputeMPVandGain(), algoEndJob(), algoEndRun(), bookDQMHistos(), processEvent(), produceTagFilter(), 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)=="") return 0;  // return StdBunch calibration mode for backward compatibility

    return None;
}

void SiStripGainFromCalibTree::storeOnTree ( TFileService *  tfs )

private

Definition at line 1065 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, GOOD, stAPVGain::Index, stAPVGain::isMasked, m_calibrationMode, TFileService::make(), MASKED, NClusterPixel, NClusterStrip, stAPVGain::NEntries, NEvent, NPixelDets, NStripAPVs, NTrack, NULL, OutputGains, stAPVGain::Phi, stAPVGain::PreviousGain, stAPVGain::PreviousGainTick, stAPVGain::R, statCollectionFromMode(), 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==NULL)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 598 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 192 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::AllowSaturation

private

Definition at line 175 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

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

private

Definition at line 253 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 253 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

__gnu_cxx::hash_map<unsigned int, stAPVGain*, __gnu_cxx::hash<unsigned int>, isEqual > SiStripGainFromCalibTree::APVsColl

private

Definition at line 270 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 269 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), algoBeginRun(), getNewObject(), and storeOnTree().

unsigned int SiStripGainFromCalibTree::BAD

private

Definition at line 220 of file SiStripGainFromCalibTree.cc.

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

string SiStripGainFromCalibTree::CalibPrefix_

private

Definition at line 260 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

int SiStripGainFromCalibTree::CalibrationLevel

private

Definition at line 179 of file SiStripGainFromCalibTree.cc.

Referenced by algoComputeMPVandGain(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::CalibSuffix_

private

Definition at line 261 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 250 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 250 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 200 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 208 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 209 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 206 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 207 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 202 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 205 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 204 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 203 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 252 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 252 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

DQMStore* SiStripGainFromCalibTree::dbe

private

Definition at line 162 of file SiStripGainFromCalibTree.cc.

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

std::vector<string> SiStripGainFromCalibTree::dqm_tag_

private

Definition at line 198 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::ERun

private

Definition at line 218 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and processEvent().

unsigned int SiStripGainFromCalibTree::eventnumber =0

private

Definition at line 226 of file SiStripGainFromCalibTree.cc.

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

string SiStripGainFromCalibTree::EventPrefix_

private

Definition at line 256 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::EventSuffix_

private

Definition at line 257 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 249 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 249 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

bool SiStripGainFromCalibTree::FirstSetOfConstants

private

Definition at line 176 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

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

private

Definition at line 245 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 245 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 254 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 254 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::GOOD

private

Definition at line 219 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 242 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 242 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 243 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 243 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 244 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 244 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::m_calibrationMode

private

Definition at line 185 of file SiStripGainFromCalibTree.cc.

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

string SiStripGainFromCalibTree::m_calibrationPath

private

Definition at line 186 of file SiStripGainFromCalibTree.cc.

Referenced by MakeCalibrationMap(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::m_DQMdir

private

Definition at line 187 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::m_harvestingMode

private

Definition at line 183 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::m_splitDQMstat

private

Definition at line 184 of file SiStripGainFromCalibTree.cc.

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

double SiStripGainFromCalibTree::MagFieldCurrentTh

private

Definition at line 163 of file SiStripGainFromCalibTree.cc.

Referenced by isBFieldConsistentWithMode(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::MASKED

private

Definition at line 221 of file SiStripGainFromCalibTree.cc.

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

double SiStripGainFromCalibTree::MaxChi2OverNDF

private

Definition at line 166 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxMPVError

private

Definition at line 165 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::MaxNrStrips

private

Definition at line 171 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxTrackChiOverNdf

private

Definition at line 173 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxTrackEta

private

Definition at line 170 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

int SiStripGainFromCalibTree::MaxTrackingIteration

private

Definition at line 174 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MaxTrackMomentum

private

Definition at line 168 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MinNrEntries

private

Definition at line 164 of file SiStripGainFromCalibTree.cc.

Referenced by getPeakOfLandau(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MinTrackEta

private

Definition at line 169 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::MinTrackHits

private

Definition at line 172 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::MinTrackMomentum

private

Definition at line 167 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::NClusterPixel

private

Definition at line 214 of file SiStripGainFromCalibTree.cc.

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

unsigned int SiStripGainFromCalibTree::NClusterStrip

private

Definition at line 213 of file SiStripGainFromCalibTree.cc.

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

unsigned int SiStripGainFromCalibTree::NEvent

private

Definition at line 211 of file SiStripGainFromCalibTree.cc.

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

int SiStripGainFromCalibTree::NPixelDets

private

Definition at line 216 of file SiStripGainFromCalibTree.cc.

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

int SiStripGainFromCalibTree::NStripAPVs

private

Definition at line 215 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 246 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 246 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::NTrack

private

Definition at line 212 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::OldGainRemoving

private

Definition at line 178 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

std::string SiStripGainFromCalibTree::OutputGains

private

Definition at line 193 of file SiStripGainFromCalibTree.cc.

Referenced by SiStripGainFromCalibTree(), and storeOnTree().

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

private

Definition at line 248 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 248 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

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

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 241 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 241 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::runnumber =0

private

Definition at line 227 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 247 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 247 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

bool SiStripGainFromCalibTree::saveSummary

private

Definition at line 181 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob(), and SiStripGainFromCalibTree().

unsigned int SiStripGainFromCalibTree::SRun

private

Definition at line 217 of file SiStripGainFromCalibTree.cc.

Referenced by algoBeginJob(), and processEvent().

double SiStripGainFromCalibTree::tagCondition_GoodFrac

private

Definition at line 190 of file SiStripGainFromCalibTree.cc.

Referenced by produceTagFilter(), and SiStripGainFromCalibTree().

double SiStripGainFromCalibTree::tagCondition_NClusters

private

Definition at line 189 of file SiStripGainFromCalibTree.cc.

Referenced by produceTagFilter(), and SiStripGainFromCalibTree().

TFileService* SiStripGainFromCalibTree::tfs

private

Definition at line 161 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob().

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

private

Definition at line 237 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 237 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 231 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 231 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 234 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 234 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 236 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 236 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 240 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 240 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 232 of file SiStripGainFromCalibTree.cc.

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

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

private

Definition at line 232 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

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

private

Definition at line 235 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 235 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::TrackPrefix_

private

Definition at line 258 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 233 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 233 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

string SiStripGainFromCalibTree::TrackSuffix_

private

Definition at line 259 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().

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

private

Definition at line 228 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and algoAnalyzeTheTree().

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

private

Definition at line 228 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyze(), and SiStripGainFromCalibTree().

bool SiStripGainFromCalibTree::useCalibration

private

Definition at line 182 of file SiStripGainFromCalibTree.cc.

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

bool SiStripGainFromCalibTree::Validation

private

Definition at line 177 of file SiStripGainFromCalibTree.cc.

Referenced by processEvent(), and SiStripGainFromCalibTree().

vector<string> SiStripGainFromCalibTree::VInputFiles

private

Definition at line 194 of file SiStripGainFromCalibTree.cc.

Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().