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

Inheritance diagram for HcalCalibrator:

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

	~HcalCalibrator ()

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 (const std::string &iProcessName, std::vector< const char * > &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 Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob ()

virtual void	endJob ()

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

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

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 25 of file HcalCalibrator.h.

Constructor & Destructor Documentation

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

explicit

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

Definition at line 115 of file HcalCalibrator.cc.

115 {

116 }

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 121 of file HcalCalibrator.cc.

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

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

privatevirtual

Reimplemented from edm::EDAnalyzer.

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

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 149 of file HcalCalibrator.cc.

References gather_cfg::cout, MillePedeFileConverter_cfg::fileName, MainPageGenerator::files, diJetCalib::inputFileList, 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().

                        {
  
 
   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 83 of file HcalCalibrator.h.

bool HcalCalibrator::mApplyPhiSymCorFlag

private

Definition at line 73 of file HcalCalibrator.h.

Referenced by endJob().

int HcalCalibrator::mCalibAbsIEtaMax

private

Definition at line 62 of file HcalCalibrator.h.

Referenced by endJob().

int HcalCalibrator::mCalibAbsIEtaMin

private

Definition at line 63 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mCalibMethod

private

Definition at line 43 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mCalibType

private

Definition at line 42 of file HcalCalibrator.h.

Referenced by endJob().

bool HcalCalibrator::mCombinePhi

private

Definition at line 55 of file HcalCalibrator.h.

Referenced by endJob().

int HcalCalibrator::mHbClusterSize

private

Definition at line 56 of file HcalCalibrator.h.

Referenced by endJob().

int HcalCalibrator::mHeClusterSize

private

Definition at line 57 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mHistoFileName

private

Definition at line 76 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mInputFileList

private

Definition at line 40 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxConeDist

private

Definition at line 60 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxEOverP

private

Definition at line 48 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxEtThirdJet

private

Definition at line 51 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxProbeJetEmFrac

private

Definition at line 65 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxTagJetAbsEta

private

Definition at line 67 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxTagJetEmFrac

private

Definition at line 66 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxTargetE

private

Definition at line 45 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMaxTrkEmE

private

Definition at line 49 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinCellE

private

Definition at line 46 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinDPhiDiJets

private

Definition at line 52 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinEOverP

private

Definition at line 47 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinProbeJetAbsEta

private

Definition at line 69 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinTagJetEt

private

Definition at line 68 of file HcalCalibrator.h.

Referenced by endJob().

double HcalCalibrator::mMinTargetE

private

Definition at line 44 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mOutputCorCoefFileName

private

Definition at line 75 of file HcalCalibrator.h.

Referenced by endJob().

std::string HcalCalibrator::mOutputFile

private

Definition at line 39 of file HcalCalibrator.h.

std::string HcalCalibrator::mPhiSymCorFileName

private

Definition at line 72 of file HcalCalibrator.h.

Referenced by endJob().

bool HcalCalibrator::mSumDepths

private

Definition at line 53 of file HcalCalibrator.h.

Referenced by endJob().

bool HcalCalibrator::mSumSmallDepths

private

Definition at line 54 of file HcalCalibrator.h.

Referenced by endJob().

const CaloGeometry* HcalCalibrator::mTheCaloGeometry

private

Definition at line 79 of file HcalCalibrator.h.

Referenced by analyze(), and endJob().

const HcalTopology* HcalCalibrator::mTheHcalTopology

private

Definition at line 80 of file HcalCalibrator.h.

Referenced by analyze(), and endJob().

bool HcalCalibrator::mUseConeClustering

private

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