#include <MyPackages/HcalCalibrator/src/HcalCalibrator.cc>

Inheritance diagram for HcalCalibrator:

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

	~HcalCalibrator () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

void	beginJob () override

void	endJob () override

Private Attributes
bool	allowMissingInputs_

bool	mApplyPhiSymCorFlag

int	mCalibAbsIEtaMax

int	mCalibAbsIEtaMin

std::string	mCalibMethod

std::string	mCalibType

bool	mCombinePhi

int	mHbClusterSize

int	mHeClusterSize

std::string	mHistoFileName

std::string	mInputFileList

double	mMaxConeDist

double	mMaxEOverP

double	mMaxEtThirdJet

double	mMaxProbeJetEmFrac

double	mMaxTagJetAbsEta

double	mMaxTagJetEmFrac

double	mMaxTargetE

double	mMaxTrkEmE

double	mMinCellE

double	mMinDPhiDiJets

double	mMinEOverP

double	mMinProbeJetAbsEta

double	mMinTagJetEt

double	mMinTargetE

std::string	mOutputCorCoefFileName

std::string	mOutputFile

std::string	mPhiSymCorFileName

bool	mSumDepths

bool	mSumSmallDepths

const CaloGeometry *	mTheCaloGeometry

const HcalTopology *	mTheHcalTopology

bool	mUseConeClustering

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

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

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

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

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

void	consumesMany (const TypeToGet &id)

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

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

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

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

Detailed Description

Description: <one line="" class="" summary>="">

Implementation:

This is an interface to run the the hcal calibration code for isolated tracsk and dijets. It takes tha parameters for the calibration settings from a python file and passes them to the actual calibration code in "endJob()".

Definition at line 22 of file HcalCalibrator.h.

Constructor & Destructor Documentation

HcalCalibrator::HcalCalibrator ( const edm::ParameterSet & conf )

explicit

Definition at line 64 of file HcalCalibrator.cc.

     :
 
       mInputFileList(conf.getUntrackedParameter<string>("inputFileList")),
       //  mOutputFile(conf.getUntrackedParameter<string>("outputFile")),
 
       mCalibType(conf.getUntrackedParameter<string>("calibType")),
       mCalibMethod(conf.getUntrackedParameter<string>("calibMethod")),
       mMinTargetE(conf.getUntrackedParameter<double>("minTargetE")),
       mMaxTargetE(conf.getUntrackedParameter<double>("maxTargetE")),
       mMinCellE(conf.getUntrackedParameter<double>("minCellE")),
       mMinEOverP(conf.getUntrackedParameter<double>("minEOverP")),
       mMaxEOverP(conf.getUntrackedParameter<double>("maxEOverP")),
       mMaxTrkEmE(conf.getUntrackedParameter<double>("maxTrkEmE")),
       mMaxEtThirdJet(conf.getUntrackedParameter<double>("maxEtThirdJet")),
       mMinDPhiDiJets(conf.getUntrackedParameter<double>("minDPhiDiJets")),
       mSumDepths(conf.getUntrackedParameter<bool>("sumDepths")),
       mSumSmallDepths(conf.getUntrackedParameter<bool>("sumSmallDepths")),
       mCombinePhi(conf.getUntrackedParameter<bool>("combinePhi")),
       mHbClusterSize(conf.getUntrackedParameter<int>("hbClusterSize")),
       mHeClusterSize(conf.getUntrackedParameter<int>("heClusterSize")),
 
       mUseConeClustering(conf.getUntrackedParameter<bool>("useConeClustering")),
       mMaxConeDist(conf.getUntrackedParameter<double>("maxConeDist")),
 
       mCalibAbsIEtaMax(conf.getUntrackedParameter<int>("calibAbsIEtaMax")),
       mCalibAbsIEtaMin(conf.getUntrackedParameter<int>("calibAbsIEtaMin")),
       mMaxProbeJetEmFrac(conf.getUntrackedParameter<double>("maxProbeJetEmFrac")),
       mMaxTagJetEmFrac(conf.getUntrackedParameter<double>("maxTagJetEmFrac")),
       mMaxTagJetAbsEta(conf.getUntrackedParameter<double>("maxTagJetAbsEta")),
       mMinTagJetEt(conf.getUntrackedParameter<double>("minTagJetEt")),
       mMinProbeJetAbsEta(conf.getUntrackedParameter<double>("minProbeJetAbsEta")),
       mPhiSymCorFileName(conf.getUntrackedParameter<string>("phiSymCorFileName")),
       mApplyPhiSymCorFlag(conf.getUntrackedParameter<bool>("applyPhiSymCorFlag")),
       mOutputCorCoefFileName(conf.getUntrackedParameter<string>("outputCorCoefFileName")),
       mHistoFileName(conf.getUntrackedParameter<string>("histoFileName"))
 
 {}

HcalCalibrator::~HcalCalibrator ( )

override

Definition at line 105 of file HcalCalibrator.cc.

105 {}

Member Function Documentation

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

overrideprivate

Definition at line 109 of file HcalCalibrator.cc.

References edm::EventSetup::get(), mTheCaloGeometry, mTheHcalTopology, edm::ESHandle< T >::product(), and PVValHelper::pT.

                                                                                 {
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   mTheCaloGeometry = pG.product();
   edm::ESHandle<HcalTopology> pT;
   iSetup.get<HcalRecNumberingRecord>().get(pT);
   mTheHcalTopology = pT.product();
 }

void HcalCalibrator::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 120 of file HcalCalibrator.cc.

                               {
   //  ESHandle<CaloGeometry> theGeometry;
   //  ESHandle<CaloSubdetectorGeometry> theEndcapGeometry_handle, theBarrelGeometry_handle;
   //  evtSetup.get<CaloGeometryRecord>().get( theGeometry );
 
   // edm::ESHandle<CaloGeometry> pG;
   // evtSetup.get<CaloGeometryRecord>().get(pG);
   // mTheCaloGeometry = pG.product();
 }

void HcalCalibrator::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 132 of file HcalCalibrator.cc.

References gather_cfg::cout, DEFINE_FWK_MODULE, MillePedeFileConverter_cfg::fileName, PatBasicFWLiteAnalyzer::files, merge::inputFiles, mApplyPhiSymCorFlag, mCalibAbsIEtaMax, mCalibAbsIEtaMin, mCalibMethod, mCalibType, mCombinePhi, mHbClusterSize, mHeClusterSize, mHistoFileName, mInputFileList, mMaxConeDist, mMaxEOverP, mMaxEtThirdJet, mMaxProbeJetEmFrac, mMaxTagJetAbsEta, mMaxTagJetEmFrac, mMaxTargetE, mMaxTrkEmE, mMinCellE, mMinDPhiDiJets, mMinEOverP, mMinProbeJetAbsEta, mMinTagJetEt, mMinTargetE, mOutputCorCoefFileName, mPhiSymCorFileName, mSumDepths, mSumSmallDepths, mTheCaloGeometry, mTheHcalTopology, mUseConeClustering, hcalCalib::SetApplyPhiSymCorFlag(), hcalCalib::SetCalibAbsIEtaMax(), hcalCalib::SetCalibAbsIEtaMin(), hcalCalib::SetCalibMethod(), hcalCalib::SetCalibType(), hcalCalib::SetCaloGeometry(), hcalCalib::SetCombinePhiFlag(), hcalCalib::SetConeMaxDist(), hcalCalib::SetHbClusterSize(), hcalCalib::SetHeClusterSize(), hcalCalib::SetHistoFileName(), hcalCalib::SetMaxEOverP(), hcalCalib::SetMaxEtThirdJet(), hcalCalib::SetMaxProbeJetEmFrac(), hcalCalib::SetMaxTagJetAbsEta(), hcalCalib::SetMaxTagJetEmFrac(), hcalCalib::SetMaxTargetE(), hcalCalib::SetMaxTrkEmE(), hcalCalib::SetMinCellE(), hcalCalib::SetMinDPhiDiJets(), hcalCalib::SetMinEOverP(), hcalCalib::SetMinProbeJetAbsEta(), hcalCalib::SetMinTagJetEt(), hcalCalib::SetMinTargetE(), hcalCalib::SetOutputCorCoefFileName(), hcalCalib::SetPhiSymCorFileName(), hcalCalib::SetSumDepthsFlag(), hcalCalib::SetSumSmallDepthsFlag(), and hcalCalib::SetUseConeClustering().

Referenced by o2olib.O2ORunMgr::executeJob().

                             {
   if (mCalibType != "DI_JET" && mCalibType != "ISO_TRACK") {
     std::cout << "\n\nUnknown calibration type " << mCalibType << std::endl;
     std::cout << "Please select ISO_TRACK or DI_JET in the python file." << std::endl;
     return;
   }
 
   if (mCalibMethod != "L3" && mCalibMethod != "MATRIX_INV_OF_ETA_AVE" && mCalibMethod != "L3_AND_MTRX_INV") {
     std::cout << "\n\nUnknown calibration method " << mCalibMethod << std::endl;
     std::cout << "Supported methods for IsoTrack calibration are: L3, MATRIX_INV_OF_ETA_AVE, L3_AND_MTRX_INV"
               << std::endl;
     std::cout << "For DiJets the supported method is L3" << std::endl;
     return;
   }
 
   if (mCalibType == "DI_JET" && mCalibMethod != "L3") {
     std::cout << "\n\nDiJet calibration can use only the L3 method. Please change the python file." << std::endl;
     return;
   }
 
   if (mCalibAbsIEtaMin < 1 || mCalibAbsIEtaMax > 41 || mCalibAbsIEtaMin > mCalibAbsIEtaMax) {
     std::cout
         << "\n\nInvalid ABS(iEta) calibration range. Check calibAbsIEtaMin and calibAbsIEtaMax in the python file."
         << std::endl;
     return;
   }
 
   hcalCalib* calibrator = new hcalCalib();
 
   // set the parameters controlling the calibratoration
 
   calibrator->SetCalibType(mCalibType);
   calibrator->SetCalibMethod(mCalibMethod);
   calibrator->SetMinTargetE(mMinTargetE);
   calibrator->SetMaxTargetE(mMaxTargetE);
   calibrator->SetMaxEtThirdJet(mMaxEtThirdJet);
   calibrator->SetMinDPhiDiJets(mMinDPhiDiJets);
   calibrator->SetSumDepthsFlag(mSumDepths);
   calibrator->SetSumSmallDepthsFlag(mSumSmallDepths);
   calibrator->SetCombinePhiFlag(mCombinePhi);
   calibrator->SetMinCellE(mMinCellE);
   calibrator->SetMinEOverP(mMinEOverP);
   calibrator->SetMaxEOverP(mMaxEOverP);
   calibrator->SetMaxTrkEmE(mMaxTrkEmE);
   calibrator->SetHbClusterSize(mHbClusterSize);
   calibrator->SetHeClusterSize(mHeClusterSize);
 
   calibrator->SetUseConeClustering(mUseConeClustering);
   calibrator->SetConeMaxDist(mMaxConeDist);
 
   calibrator->SetCalibAbsIEtaMax(mCalibAbsIEtaMax);
   calibrator->SetCalibAbsIEtaMin(mCalibAbsIEtaMin);
   calibrator->SetMaxProbeJetEmFrac(mMaxProbeJetEmFrac);
   calibrator->SetMaxTagJetEmFrac(mMaxTagJetEmFrac);
   calibrator->SetMaxTagJetAbsEta(mMaxTagJetAbsEta);
 
   calibrator->SetMinTagJetEt(mMinTagJetEt);
 
   calibrator->SetMinProbeJetAbsEta(mMinProbeJetAbsEta);
   calibrator->SetApplyPhiSymCorFlag(mApplyPhiSymCorFlag);
   calibrator->SetPhiSymCorFileName(mPhiSymCorFileName);
   calibrator->SetOutputCorCoefFileName(mOutputCorCoefFileName);
 
   calibrator->SetHistoFileName(mHistoFileName);
 
   calibrator->SetCaloGeometry(mTheCaloGeometry, mTheHcalTopology);
 
   std::ifstream inputFileList;  // contains list of input root files
 
   TString files = mInputFileList;
   inputFileList.open(files.Data());
 
   std::vector<TString> inputFiles;
   while (!inputFileList.eof()) {
     TString fileName;
     inputFileList >> fileName;
     if (!fileName.BeginsWith("#") && !fileName.Contains(" ") && fileName != "")
       inputFiles.push_back(fileName);
   }
   inputFileList.close();
 
   std::cout << "\nInput files for processing:" << std::endl;
   for (std::vector<TString>::iterator it = inputFiles.begin(); it != inputFiles.end(); ++it) {
     std::cout << "file: " << it->Data() << std::endl;
   }
   std::cout << std::endl;
 
   TChain* fChain = new TChain("hcalCalibTree");
 
   for (std::vector<TString>::iterator f_it = inputFiles.begin(); f_it != inputFiles.end(); ++f_it) {
     fChain->Add(f_it->Data());
   }
 
   fChain->Process(calibrator);
 
   if (fChain)
     delete fChain;
   delete calibrator;
 
   return;
 }

Member Data Documentation

bool HcalCalibrator::allowMissingInputs_

private

Definition at line 75 of file HcalCalibrator.h.

bool HcalCalibrator::mApplyPhiSymCorFlag

private

Definition at line 67 of file HcalCalibrator.h.

Referenced by endJob().

int HcalCalibrator::mCalibAbsIEtaMax

private

Definition at line 57 of file HcalCalibrator.h.

Referenced by endJob().

int HcalCalibrator::mCalibAbsIEtaMin

private

Definition at line 58 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mCalibMethod

private

Definition at line 38 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mCalibType

private

Definition at line 37 of file HcalCalibrator.h.

Referenced by endJob().

bool HcalCalibrator::mCombinePhi

private

Definition at line 50 of file HcalCalibrator.h.

Referenced by endJob().

int HcalCalibrator::mHbClusterSize

private

Definition at line 51 of file HcalCalibrator.h.

Referenced by endJob().

int HcalCalibrator::mHeClusterSize

private

Definition at line 52 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mHistoFileName

private

Definition at line 70 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mInputFileList

private

Definition at line 35 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxConeDist

private

Definition at line 55 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxEOverP

private

Definition at line 43 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxEtThirdJet

private

Definition at line 46 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxProbeJetEmFrac

private

Definition at line 60 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxTagJetAbsEta

private

Definition at line 62 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxTagJetEmFrac

private

Definition at line 61 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxTargetE

private

Definition at line 40 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxTrkEmE

private

Definition at line 44 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinCellE

private

Definition at line 41 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinDPhiDiJets

private

Definition at line 47 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinEOverP

private

Definition at line 42 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinProbeJetAbsEta

private

Definition at line 64 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinTagJetEt

private

Definition at line 63 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinTargetE

private

Definition at line 39 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mOutputCorCoefFileName

private

Definition at line 69 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mOutputFile

private

Definition at line 34 of file HcalCalibrator.h.

std::string HcalCalibrator::mPhiSymCorFileName

private

Definition at line 66 of file HcalCalibrator.h.

Referenced by endJob().

bool HcalCalibrator::mSumDepths

private

Definition at line 48 of file HcalCalibrator.h.

Referenced by endJob().

bool HcalCalibrator::mSumSmallDepths

private

Definition at line 49 of file HcalCalibrator.h.

Referenced by endJob().

const CaloGeometry* HcalCalibrator::mTheCaloGeometry

private

Definition at line 72 of file HcalCalibrator.h.

Referenced by analyze(), and endJob().

const HcalTopology* HcalCalibrator::mTheHcalTopology

private

Definition at line 73 of file HcalCalibrator.h.

Referenced by analyze(), and endJob().

bool HcalCalibrator::mUseConeClustering

private

Definition at line 54 of file HcalCalibrator.h.

Referenced by endJob().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation