CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
CTPPSPixGainCalibsESAnalyzer Class Reference
Inheritance diagram for CTPPSPixGainCalibsESAnalyzer:
edm::one::EDAnalyzer<> edm::one::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

void analyze (const edm::Event &e, const edm::EventSetup &c) override
 
 CTPPSPixGainCalibsESAnalyzer (edm::ParameterSet const &p)
 
 CTPPSPixGainCalibsESAnalyzer (int i)
 
void setReadablePlaneNames ()
 
 ~CTPPSPixGainCalibsESAnalyzer () override
 
- Public Member Functions inherited from edm::one::EDAnalyzer<>
 EDAnalyzer ()=default
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () 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
 
EDConsumerBaseoperator= (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 Attributes

std::map< uint32_t, std::string > detId_readable
 
std::string m_outfilename
 

Additional Inherited Members

- Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 15 of file CTPPSPixGainCalibsESAnalyzer.cc.

Constructor & Destructor Documentation

CTPPSPixGainCalibsESAnalyzer::CTPPSPixGainCalibsESAnalyzer ( edm::ParameterSet const &  p)
inlineexplicit

Definition at line 17 of file CTPPSPixGainCalibsESAnalyzer.cc.

References setReadablePlaneNames().

17  :
18  m_outfilename(p.getUntrackedParameter<std::string>("outputrootfile","output.root")) {
19  // std::cout<<"CTPPSPixGainCalibsESAnalyzer"<<std::endl;
21  }
CTPPSPixGainCalibsESAnalyzer::CTPPSPixGainCalibsESAnalyzer ( int  i)
inlineexplicit

Definition at line 22 of file CTPPSPixGainCalibsESAnalyzer.cc.

References setReadablePlaneNames().

22  {
23  //std::cout<<"CTPPSPixGainCalibsESAnalyzer "<<i<<std::endl;
25  }
CTPPSPixGainCalibsESAnalyzer::~CTPPSPixGainCalibsESAnalyzer ( )
inlineoverride

Definition at line 26 of file CTPPSPixGainCalibsESAnalyzer.cc.

References analyze(), EnergyCorrector::c, MillePedeFileConverter_cfg::e, and setReadablePlaneNames().

26  {
27  //std::cout<<"~CTPPSPixGainCalibsESAnalyzer "<<std::endl;
28  }

Member Function Documentation

void CTPPSPixGainCalibsESAnalyzer::analyze ( const edm::Event e,
const edm::EventSetup c 
)
override

Definition at line 65 of file CTPPSPixGainCalibsESAnalyzer.cc.

References DEFINE_FWK_MODULE, detId_readable, edm::EventID::event(), edm::eventsetup::heterocontainer::HCTypeTag::findType(), edm::EventSetup::get(), CTPPSPixelGainCalibrations::getCalibMap(), CTPPSPixelGainCalibration::getGain(), CTPPSPixelGainCalibrations::getGainCalibration(), CTPPSPixelGainCalibration::getIEnd(), CTPPSPixelGainCalibration::getNCols(), CTPPSPixelGainCalibration::getNRows(), CTPPSPixelGainCalibration::getPed(), edm::EventBase::id(), CTPPSPixelGainCalibration::isDead(), CTPPSPixelGainCalibration::isNoisy(), m_outfilename, hgcalPlots::ncols, edm::ESHandle< T >::product(), edm::EventID::run(), CTPPSPixelGainCalibrations::size(), AlCaHLTBitMon_QueryRunRegistry::string, and edm::eventsetup::EventSetupRecordKey::type().

Referenced by ~CTPPSPixGainCalibsESAnalyzer().

65  {
66  edm::LogPrint("CTPPSPixGainCalibsReader")<<"###CTPPSPixGainCalibsESAnalyzer::analyze";
67  edm::LogPrint("CTPPSPixGainCalibsReader") <<" I AM IN RUN NUMBER "<<e.id().run() ;
68  edm::LogPrint("CTPPSPixGainCalibsReader") <<" ---EVENT NUMBER "<<e.id().event() ;
70  if( recordKey.type() == edm::eventsetup::EventSetupRecordKey::TypeTag()) {
71  //record not found
72  edm::LogPrint("CTPPSPixGainCalibsReader") <<"Record \"CTPPSPixelGainCalibrationsRcd"<<"\" does not exist ";
73  }
74  //this part gets the handle of the event source and the record (i.e. the Database)
76  edm::LogPrint("CTPPSPixGainCalibsReader")<<"got eshandle";
77  context.get<CTPPSPixelGainCalibrationsRcd>().get(calhandle);
78  edm::LogPrint("CTPPSPixGainCalibsReader")<<"got context";
79  const CTPPSPixelGainCalibrations* pPixelGainCalibrations=calhandle.product();
80  edm::LogPrint("CTPPSPixGainCalibsReader")<<"got CTPPSPixelGainCalibrations* ";
81  edm::LogPrint("CTPPSPixGainCalibsReader")<< "print pointer address : " ;
82  edm::LogPrint("CTPPSPixGainCalibsReader") << pPixelGainCalibrations ;
83 
84  TFile myfile(m_outfilename.c_str(),"RECREATE");
85  myfile.cd();
86 
87  // the pPixelGainCalibrations object contains the map of detIds to pixel gains and pedestals for current run
88  // we get the map just to loop over the contents, but from here on it should be as the code (reconstruction etc) needs.
89  // Probably best to check that the key (detid) is in the list before calling the data
90 
91  edm::LogPrint("CTPPSPixGainCalibsReader") << "Size " << pPixelGainCalibrations->size();
92  const CTPPSPixelGainCalibrations::CalibMap & mymap = pPixelGainCalibrations->getCalibMap(); //just to get the keys?
93 
94  for (CTPPSPixelGainCalibrations::CalibMap::const_iterator it = mymap.begin(); it !=mymap.end() ; ++it) {
95  uint32_t detId = it->first;
96 
97  edm::LogPrint("CTPPSPixGainCalibsReader")<<"Address of detId = "<< (&detId) << " and of it = "<< (&it) << " and of it->first = "<< (&(it->first)) ;
98 
99  edm::LogPrint("CTPPSPixGainCalibsReader") << "Content of pPixelGainCalibrations for key: detId= " << detId ;
100  CTPPSPixelGainCalibration mycalibs0= pPixelGainCalibrations->getGainCalibration(detId);
101  const CTPPSPixelGainCalibration & mycalibs=it->second;
102 
103  edm::LogPrint("CTPPSPixGainCalibsReader")<<"Address of mycalibs0 = "<< (&mycalibs0) << " and of mycalibs = "<< (&mycalibs) << " and of it->second "<<(&(it->second));
104 
105 
106  std::string namep("pedsFromDB_" + detId_readable[detId]);
107  std::string nameg("gainsFromDB_" + detId_readable[detId]);
108  std::string tlp("Pedestals for "+detId_readable[detId]+"; column; row");
109  std::string tlg("Gains for "+detId_readable[detId]+"; column; row");
110  TH2D mypeds(namep.c_str(),tlp.c_str(),156,0.,156.,160,0.,160.);
111  TH2D mygains(nameg.c_str(),tlg.c_str(),156,0.,156.,160,0.,160.);
112 
113  int ncols = mycalibs.getNCols();
114  int npix = mycalibs.getIEnd();
115  int nrows = mycalibs.getNRows(); //should be == 160
116  edm::LogPrint("CTPPSPixGainCalibsReader")<<"Here ncols = "<<ncols<<" nrows ="<<nrows <<" npix="<<npix ;
117  for (int jrow=0; jrow < nrows ; ++jrow)
118  for(int icol=0; icol < ncols ; ++icol){
119  if(mycalibs.isDead(icol+jrow*ncols))
120  {edm::LogPrint("CTPPSPixGainCalibsReader")<<"Dead Pixel icol ="<< icol << " jrow ="<< jrow ; continue;}
121  if(mycalibs.isNoisy(icol+jrow*ncols))
122  {edm::LogPrint("CTPPSPixGainCalibsReader")<<"Noisy Pixel icol ="<< icol << " jrow ="<< jrow ;continue;}
123  mygains.Fill(icol,jrow,mycalibs.getGain(icol,jrow));
124  mypeds.Fill(icol,jrow,mycalibs.getPed(icol,jrow));
125  }
126 
127  mypeds.Write();
128  mygains.Write();
129  }
130  myfile.Write();
131  myfile.Close();
132 }
RunNumber_t run() const
Definition: EventID.h:39
EventNumber_t event() const
Definition: EventID.h:41
const CalibMap & getCalibMap() const
bool isDead(const uint32_t ipix) const
const CTPPSPixelGainCalibration & getGainCalibration(const uint32_t &detid) const
std::map< uint32_t, CTPPSPixelGainCalibration > CalibMap
static HCTypeTag findType(char const *iTypeName)
find a type based on the types name, if not found will return default HCTypeTag
std::map< uint32_t, std::string > detId_readable
heterocontainer::HCTypeTag TypeTag
edm::EventID id() const
Definition: EventBase.h:59
float getPed(const int &col, const int &row) const
bool isNoisy(const uint32_t ipix) const
T const * product() const
Definition: ESHandle.h:86
float getGain(const int &col, const int &row) const
void CTPPSPixGainCalibsESAnalyzer::setReadablePlaneNames ( )

Definition at line 37 of file CTPPSPixGainCalibsESAnalyzer.cc.

References detId_readable.

Referenced by CTPPSPixGainCalibsESAnalyzer(), and ~CTPPSPixGainCalibsESAnalyzer().

37  {
38  detId_readable[2014838784]="Arm_0_Sec_45_St_0_Pot_3_Plane_0";
39  detId_readable[2014904320]="Arm_0_Sec_45_St_0_Pot_3_Plane_1";
40  detId_readable[2014969856]="Arm_0_Sec_45_St_0_Pot_3_Plane_2";
41  detId_readable[2015035392]="Arm_0_Sec_45_St_0_Pot_3_Plane_3";
42  detId_readable[2015100928]="Arm_0_Sec_45_St_0_Pot_3_Plane_4";
43  detId_readable[2015166464]="Arm_0_Sec_45_St_0_Pot_3_Plane_5";
44  detId_readable[2023227392]="Arm_0_Sec_45_St_2_Pot_3_Plane_0";
45  detId_readable[2023292928]="Arm_0_Sec_45_St_2_Pot_3_Plane_1";
46  detId_readable[2023358464]="Arm_0_Sec_45_St_2_Pot_3_Plane_2";
47  detId_readable[2023424000]="Arm_0_Sec_45_St_2_Pot_3_Plane_3";
48  detId_readable[2023489536]="Arm_0_Sec_45_St_2_Pot_3_Plane_4";
49  detId_readable[2023555072]="Arm_0_Sec_45_St_2_Pot_3_Plane_5";
50  detId_readable[2031616000]="Arm_1_Sec_56_St_0_Pot_3_Plane_0";
51  detId_readable[2031681536]="Arm_1_Sec_56_St_0_Pot_3_Plane_1";
52  detId_readable[2031747072]="Arm_1_Sec_56_St_0_Pot_3_Plane_2";
53  detId_readable[2031812608]="Arm_1_Sec_56_St_0_Pot_3_Plane_3";
54  detId_readable[2031878144]="Arm_1_Sec_56_St_0_Pot_3_Plane_4";
55  detId_readable[2031943680]="Arm_1_Sec_56_St_0_Pot_3_Plane_5";
56  detId_readable[2040004608]="Arm_1_Sec_56_St_2_Pot_3_Plane_0";
57  detId_readable[2040070144]="Arm_1_Sec_56_St_2_Pot_3_Plane_1";
58  detId_readable[2040135680]="Arm_1_Sec_56_St_2_Pot_3_Plane_2";
59  detId_readable[2040201216]="Arm_1_Sec_56_St_2_Pot_3_Plane_3";
60  detId_readable[2040266752]="Arm_1_Sec_56_St_2_Pot_3_Plane_4";
61  detId_readable[2040332288]="Arm_1_Sec_56_St_2_Pot_3_Plane_5";
62 }
std::map< uint32_t, std::string > detId_readable

Member Data Documentation

std::map<uint32_t,std::string> CTPPSPixGainCalibsESAnalyzer::detId_readable
private

Definition at line 32 of file CTPPSPixGainCalibsESAnalyzer.cc.

Referenced by analyze(), and setReadablePlaneNames().

std::string CTPPSPixGainCalibsESAnalyzer::m_outfilename
private

Definition at line 33 of file CTPPSPixGainCalibsESAnalyzer.cc.

Referenced by analyze().