|
|
|
#include <DigiBXCorrHistogramMaker.h>
Public Member Functions | |
| void | beginRun (const unsigned int nrun) |
| void | book (const char *dirname, const std::map< int, std::string > &labels) |
| DigiBXCorrHistogramMaker (const edm::ParameterSet &iConfig) | |
| DigiBXCorrHistogramMaker (const int ncorbins=1000) | |
| void | fill (const T &he, const std::map< int, int > &ndigi, const edm::Handle< APVCyclePhaseCollection > &phase) |
| void | fillcorr (const T &he1, const T &he2, const std::map< int, int > &ndigi) |
| ~DigiBXCorrHistogramMaker () | |
Private Attributes | |
| std::map< unsigned int, int > | m_binmax |
| const bool | m_dbx3Histo |
| const bool | m_dbx3Histo3D |
| std::map< int, TProfile * > | m_digicorr |
| std::string | m_hitname |
| std::map< int, std::string > | m_labels |
| const int | m_nbins |
| int | m_ncorbins |
| std::map< int, TProfile * > | m_ndigivsbx |
| std::map< int, TH2F * > | m_ndigivsbx2D |
| std::map< int, TH2F * > | m_ndigivsbx2Dzoom |
| std::map< int, TH2F * > | m_ndigivsbx2Dzoom2 |
| std::map< int, TH2F * > | m_ndigivscycle |
| std::map< int, TProfile2D * > | m_ndigivscycle2dbx |
| std::map< int, TProfile2D * > | m_ndigivscycledbx |
| std::map< int, TProfile2D ** > | m_ndigivscycletime |
| std::map< int, TH2F * > | m_ndigivscyclezoom |
| std::map< int, TH2F * > | m_ndigivscyclezoom2 |
| std::map< int, TProfile * > | m_ndigivsdbx |
| std::map< int, TH2F * > | m_ndigivsdbx2D |
| std::map< int, TH2F * > | m_ndigivsdbx2Dzoom |
| std::map< int, TH2F * > | m_ndigivsdbx2Dzoom2 |
| std::map< int, TProfile2D * > | m_ndigivsdbx3zoom |
| std::map< int, TProfile * > | m_ndigivsdbxincycle |
| std::map< int, TH2F * > | m_ndigivsdbxincycle2D |
| std::map< int, TProfile2D * > | m_ndigivsdbxincycle3 |
| std::map< int, TH3F * > | m_ndigivsdbxincycle33D |
| std::map< int, TProfile * > | m_ndigivsdbxzoom |
| std::map< int, TProfile * > | m_ndigivsdbxzoom2 |
| std::map< int, TProfile * > | m_nmeandigivscycle |
| std::map< int, std::string > | m_phasepart |
| RunHistogramManager | m_rhm |
| const bool | m_runHisto |
| std::vector< int > | m_scalefact |
Definition at line 22 of file DigiBXCorrHistogramMaker.h.
| DigiBXCorrHistogramMaker< T >::DigiBXCorrHistogramMaker | ( | const int | ncorbins = 1000 | ) |
Definition at line 90 of file DigiBXCorrHistogramMaker.h.
: m_ncorbins(ncorbins), m_hitname("digi"), m_dbx3Histo(false), m_dbx3Histo3D(false), m_runHisto(true), m_labels(), m_binmax(), m_phasepart(), m_scalefact(), m_nbins(200), m_rhm() { }
| DigiBXCorrHistogramMaker< T >::DigiBXCorrHistogramMaker | ( | const edm::ParameterSet & | iConfig | ) |
Definition at line 96 of file DigiBXCorrHistogramMaker.h.
References edm::ParameterSet::getUntrackedParameter(), DigiBXCorrHistogramMaker< T >::m_binmax, and DigiBXCorrHistogramMaker< T >::m_phasepart.
: m_ncorbins(iConfig.getUntrackedParameter<int>("corrNbins")), m_hitname(iConfig.getUntrackedParameter<std::string>("hitName","digi")), m_dbx3Histo(iConfig.getUntrackedParameter<bool>("dbx3Histo",false)), m_dbx3Histo3D(iConfig.getUntrackedParameter<bool>("dbx3Histo3D",false)), m_runHisto(iConfig.getUntrackedParameter<bool>("runHisto",true)), m_labels(), m_scalefact(iConfig.getUntrackedParameter<std::vector<int> >("scaleFactors",std::vector<int>(1,5))), m_nbins(iConfig.getUntrackedParameter<int>("numberOfBins",200)), m_rhm() { std::vector<edm::ParameterSet> wantedsubds(iConfig.getUntrackedParameter<std::vector<edm::ParameterSet> >("wantedSubDets",std::vector<edm::ParameterSet>())); for(std::vector<edm::ParameterSet>::iterator ps=wantedsubds.begin();ps!=wantedsubds.end();++ps) { m_binmax[ps->getParameter<unsigned int>("detSelection")] = ps->getParameter<int>("binMax"); m_phasepart[ps->getParameter<unsigned int>("detSelection")] = ps->getUntrackedParameter<std::string>("phasePartition","None"); } }
| DigiBXCorrHistogramMaker< T >::~DigiBXCorrHistogramMaker | ( | ) | [inline] |
Definition at line 28 of file DigiBXCorrHistogramMaker.h.
{ };
| void DigiBXCorrHistogramMaker< T >::beginRun | ( | const unsigned int | nrun | ) |
Definition at line 318 of file DigiBXCorrHistogramMaker.h.
References i.
Referenced by MultiplicityTimeCorrelations::beginRun().
{
m_rhm.beginRun(nrun);
for(std::map<int,std::string>::const_iterator lab=m_labels.begin();lab!=m_labels.end();lab++) {
const int i = lab->first;
if(m_runHisto) {
if(m_ndigivscycletime[i]) {
(*m_ndigivscycletime[i])->GetXaxis()->SetTitle("Event 1 BX mod(70)"); (*m_ndigivscycletime[i])->GetYaxis()->SetTitle("time [Orb#]");
(*m_ndigivscycletime[i])->SetBit(TH1::kCanRebin);
}
}
}
}
| void DigiBXCorrHistogramMaker< T >::book | ( | const char * | dirname, |
| const std::map< int, std::string > & | labels | ||
| ) |
Definition at line 119 of file DigiBXCorrHistogramMaker.h.
References i, reco_calib_source_client_cfg::labels, TFileDirectory::make(), mergeVDriftHistosByStation::name, SiStripTKNumbers::nstrips(), indexGen::title, and interactiveExample::ui.
Referenced by MultiplicityTimeCorrelations::MultiplicityTimeCorrelations().
{
m_labels = labels;
edm::Service<TFileService> tfserv;
TFileDirectory subev = tfserv->mkdir(dirname);
SiStripTKNumbers trnumb;
edm::LogInfo("NumberOfBins") << "Number of Bins: " << m_nbins;
edm::LogInfo("ScaleFactors") << "x-axis range scale factors: ";
for(std::vector<int>::const_iterator sf=m_scalefact.begin();sf!=m_scalefact.end();++sf) {
edm::LogVerbatim("ScaleFactors") << *sf ;
}
edm::LogInfo("BinMaxValue") << "Setting bin max values";
for(std::map<int,std::string>::const_iterator lab=m_labels.begin();lab!=m_labels.end();lab++) {
const int i = lab->first; const std::string slab = lab->second; const unsigned int ui = i;
if(m_binmax.find(ui)==m_binmax.end()) {
edm::LogVerbatim("NotConfiguredBinMax") << "Bin max for " << lab->second
<< " not configured: " << trnumb.nstrips(i) << " used";
m_binmax[ui] = trnumb.nstrips(i);
}
edm::LogVerbatim("BinMaxValue") << "Bin max for " << lab->second << " is " << m_binmax[ui];
}
edm::LogInfo("PhasePartitions") << "Partitions for APV Cycle Phase";
for(std::map<int,std::string>::const_iterator lab=m_labels.begin();lab!=m_labels.end();lab++) {
const int i = lab->first; const std::string slab = lab->second; const unsigned int ui = i;
edm::LogVerbatim("PhasePartitions") << "Partition for " << lab->second << " is " << ((m_phasepart.find(ui)!=m_phasepart.end()) ? m_phasepart[ui] : "not found") ;
}
for(std::map<int,std::string>::const_iterator lab=m_labels.begin();lab!=m_labels.end();lab++) {
const int i = lab->first; const std::string slab = lab->second; const unsigned int ui = i;
char name[200];
char title [500];
// vs DBX
if(m_scalefact.size()>=1) {
sprintf(title,"%s %s multiplicity vs BX separation",slab.c_str(),m_hitname.c_str());
sprintf(name,"n%sdigivsdbx2D",slab.c_str());
m_ndigivsdbx2D[i] = subev.make<TH2F>(name,title,1000,-0.5,500000-0.5,m_nbins,0,m_binmax[ui]/(m_scalefact[0]*m_nbins)*m_nbins);
sprintf(name,"n%sdigivsdbx2Dzoom2",slab.c_str());
m_ndigivsdbx2Dzoom2[i] = subev.make<TH2F>(name,title,1000,-0.5,50000-0.5,m_nbins,0,m_binmax[ui]/(m_scalefact[0]*m_nbins)*m_nbins);
sprintf(name,"n%sdigivsdbx2Dzoom",slab.c_str());
m_ndigivsdbx2Dzoom[i] = subev.make<TH2F>(name,title,1000,-0.5,999.5,m_nbins,0,m_binmax[ui]/(m_scalefact[0]*m_nbins)*m_nbins);
m_ndigivsdbx2D[i]->GetXaxis()->SetTitle("#DeltaBX"); m_ndigivsdbx2D[i]->GetYaxis()->SetTitle("Number of Hits");
m_ndigivsdbx2Dzoom2[i]->GetXaxis()->SetTitle("#DeltaBX"); m_ndigivsdbx2Dzoom2[i]->GetYaxis()->SetTitle("Number of Hits");
m_ndigivsdbx2Dzoom[i]->GetXaxis()->SetTitle("#DeltaBX"); m_ndigivsdbx2Dzoom[i]->GetYaxis()->SetTitle("Number of Hits");
}
sprintf(title,"%s %s multiplicity vs BX separation",slab.c_str(),m_hitname.c_str());
sprintf(name,"n%sdigivsdbx",slab.c_str());
m_ndigivsdbx[i] = subev.make<TProfile>(name,title,1000,-0.5,500000.-0.5);
sprintf(name,"n%sdigivsdbxzoom2",slab.c_str());
m_ndigivsdbxzoom2[i] = subev.make<TProfile>(name,title,1000,-0.5,50000.-0.5);
sprintf(name,"n%sdigivsdbxzoom",slab.c_str());
m_ndigivsdbxzoom[i] = subev.make<TProfile>(name,title,1000,-0.5,999.5);
m_ndigivsdbx[i]->GetXaxis()->SetTitle("#DeltaBX"); m_ndigivsdbx[i]->GetYaxis()->SetTitle("Number of Hits");
m_ndigivsdbxzoom2[i]->GetXaxis()->SetTitle("#DeltaBX"); m_ndigivsdbxzoom2[i]->GetYaxis()->SetTitle("Number of Hits");
m_ndigivsdbxzoom[i]->GetXaxis()->SetTitle("#DeltaBX"); m_ndigivsdbxzoom[i]->GetYaxis()->SetTitle("Number of Hits");
sprintf(name,"n%sdigivsdbx3zoom",slab.c_str());
sprintf(title,"%s %s multiplicity vs Triplets BX separation",slab.c_str(),m_hitname.c_str());
m_ndigivsdbx3zoom[i] = subev.make<TProfile2D>(name,title,100,-0.5,999.5,100,-0.5,999.5);
m_ndigivsdbx3zoom[i]->GetXaxis()->SetTitle("#DeltaBX(n,n-1)"); m_ndigivsdbx3zoom[i]->GetYaxis()->SetTitle("#DeltaBX(n,n-2)");
sprintf(name,"%sdigicorr",slab.c_str());
sprintf(title,"%s %s DBX correlation",slab.c_str(),m_hitname.c_str());
m_digicorr[i] = subev.make<TProfile>(name,title,m_ncorbins,-0.5,m_ncorbins-0.5);
m_digicorr[i]->GetXaxis()->SetTitle("#DeltaBX"); m_digicorr[i]->GetYaxis()->SetTitle("Number of Hits");
// vs DBX w.r.t. cycle
if(m_scalefact.size()>=1) {
if(m_phasepart.find(ui)!=m_phasepart.end() && m_phasepart[ui]!="None") {
sprintf(name,"n%sdigivsdbxincycle",slab.c_str());
sprintf(title,"%s %s multiplicity vs BX separation w.r.t. cycle",slab.c_str(),m_hitname.c_str());
m_ndigivsdbxincycle[i] = subev.make<TProfile>(name,title,1000,-0.5,999.5);
m_ndigivsdbxincycle[i]->GetXaxis()->SetTitle("#DeltaBX w.r.t. cycle"); m_ndigivsdbxincycle[i]->GetYaxis()->SetTitle("Number of Hits");
sprintf(name,"n%sdigivsdbxincycle2D",slab.c_str());
sprintf(title,"%s %s multiplicity vs BX separation w.r.t. cycle",slab.c_str(),m_hitname.c_str());
m_ndigivsdbxincycle2D[i] = subev.make<TH2F>(name,title,1000,-0.5,999.5,m_nbins,0.,m_binmax[ui]/(m_scalefact[0]*m_nbins)*m_nbins);
m_ndigivsdbxincycle2D[i]->GetXaxis()->SetTitle("#DeltaBX w.r.t. cycle"); m_ndigivsdbxincycle2D[i]->GetYaxis()->SetTitle("Number of Hits");
if(m_dbx3Histo) {
sprintf(name,"n%sdigivsdbxincycle3",slab.c_str());
sprintf(title,"%s %s multiplicity vs Triplets BX separation w.r.t. cycle",slab.c_str(),m_hitname.c_str());
m_ndigivsdbxincycle3[i] = subev.make<TProfile2D>(name,title,2000,-0.5,1999.5,30,-0.5,2099.5);
m_ndigivsdbxincycle3[i]->GetXaxis()->SetTitle("#DeltaBX(n,n-1)");
m_ndigivsdbxincycle3[i]->GetYaxis()->SetTitle("#DeltaBX(n,n-2)-#DeltaBX(n,n-1)");
if(m_dbx3Histo3D) {
sprintf(name,"n%sdigivsdbxincycle33D",slab.c_str());
sprintf(title,"%s %s multiplicity vs Triplets BX separation w.r.t. cycle",slab.c_str(),m_hitname.c_str());
m_ndigivsdbxincycle33D[i] = subev.make<TH3F>(name,title,2000,-0.5,1999.5,30,-0.5,2099.5,50,0.,m_binmax[ui]/(m_scalefact[0]*50)*50);
m_ndigivsdbxincycle33D[i]->GetXaxis()->SetTitle("#DeltaBX(n,n-1)");
m_ndigivsdbxincycle33D[i]->GetYaxis()->SetTitle("#DeltaBX(n,n-2)-#DeltaBX(n,n-1)");
}
}
}
}
// vs absolute BX mod 70
if(m_phasepart.find(ui)!=m_phasepart.end() && m_phasepart[ui]!="None") {
sprintf(title,"%s Mean %s multiplicity vs BX mod(70)",slab.c_str(),m_hitname.c_str());
sprintf(name,"n%smeandigivscycle",slab.c_str());
m_nmeandigivscycle[i] =subev.make<TProfile>(name,title,70,-0.5,69.5);
m_nmeandigivscycle[i]->GetXaxis()->SetTitle("absolute BX mod(70)");
m_nmeandigivscycle[i]->GetYaxis()->SetTitle("Mean number of Hits");
sprintf(title,"%s %s multiplicity vs BX mod(70)",slab.c_str(),m_hitname.c_str());
if(m_scalefact.size()>=1) {
sprintf(name,"n%sdigivscycle",slab.c_str());
m_ndigivscycle[i] =subev.make<TH2F>(name,title,70,-0.5,69.5,m_nbins,0,m_binmax[ui]/(m_scalefact[0]*m_nbins)*m_nbins);
m_ndigivscycle[i]->GetXaxis()->SetTitle("absolute BX mod(70)"); m_ndigivscycle[i]->GetYaxis()->SetTitle("Number of Hits");
}
if(m_scalefact.size()>=2) {
sprintf(name,"n%sdigivscyclezoom",slab.c_str());
m_ndigivscyclezoom[i] =subev.make<TH2F>(name,title,70,-0.5,69.5,m_nbins,0,m_binmax[ui]/(m_scalefact[1]*m_nbins)*m_nbins);
m_ndigivscyclezoom[i]->GetXaxis()->SetTitle("absolute BX mod(70)");
m_ndigivscyclezoom[i]->GetYaxis()->SetTitle("Number of Hits");
}
if(m_scalefact.size()>=3) {
sprintf(name,"n%sdigivscyclezoom2",slab.c_str());
m_ndigivscyclezoom2[i] =subev.make<TH2F>(name,title,70,-0.5,69.5,m_nbins,0,m_binmax[ui]/(m_scalefact[2]*m_nbins)*m_nbins);
m_ndigivscyclezoom2[i]->GetXaxis()->SetTitle("absolute BX mod(70)");
m_ndigivscyclezoom2[i]->GetYaxis()->SetTitle("Number of Hits");
}
sprintf(name,"n%sdigivscycledbx",slab.c_str());
sprintf(title,"%s %s multiplicity vs BX mod(70) and DBX",slab.c_str(),m_hitname.c_str());
m_ndigivscycledbx[i] = subev.make<TProfile2D>(name,title,70,-0.5,69.5,1000,-0.5,999.5);
m_ndigivscycledbx[i]->GetXaxis()->SetTitle("Event 1 BX mod(70)"); m_ndigivscycledbx[i]->GetYaxis()->SetTitle("#DeltaBX event 2-1");
sprintf(name,"n%sdigivscycle2dbx",slab.c_str());
sprintf(title,"%s %s multiplicity vs BX mod(70) and DBX",slab.c_str(),m_hitname.c_str());
m_ndigivscycle2dbx[i] = subev.make<TProfile2D>(name,title,70,-0.5,69.5,1000,-0.5,999.5);
m_ndigivscycle2dbx[i]->GetXaxis()->SetTitle("Event 2 BX mod(70)"); m_ndigivscycle2dbx[i]->GetYaxis()->SetTitle("#DeltaBX event 2-1");
}
// Multiplicity in cycle vs time is booked also if the phase is not corrected
if(m_runHisto) {
sprintf(name,"n%sdigivscycletime",slab.c_str());
sprintf(title,"%s %s multiplicity vs BX mod(70) and Orbit",slab.c_str(),m_hitname.c_str());
m_ndigivscycletime[i] = m_rhm.makeTProfile2D(name,title,70,-0.5,69.5,90,0.,90.*11223);
// m_ndigivscycletime[i]->GetXaxis()->SetTitle("Event 1 BX mod(70)"); m_ndigivscycletime[i]->GetYaxis()->SetTitle("time [Orb#]");
// m_ndigivscycletime[i]->SetBit(TH1::kCanRebin);
}
// vs BX number
sprintf(title,"%s %s mean multiplicity vs BX",slab.c_str(),m_hitname.c_str());
sprintf(name,"n%sdigivsbx",slab.c_str());
m_ndigivsbx[i] =subev.make<TProfile>(name,title,3564,-0.5,3563.5);
m_ndigivsbx[i]->GetXaxis()->SetTitle("BX#"); m_ndigivsbx[i]->GetYaxis()->SetTitle("Mean Number of Hits");
sprintf(title,"%s %s multiplicity vs BX",slab.c_str(),m_hitname.c_str());
if(m_scalefact.size()>=1) {
sprintf(name,"n%sdigivsbx2D",slab.c_str());
m_ndigivsbx2D[i] =subev.make<TH2F>(name,title,3564,-0.5,3563.5,m_nbins,0,m_binmax[ui]/(m_scalefact[0]*m_nbins)*m_nbins);
m_ndigivsbx2D[i]->GetXaxis()->SetTitle("BX#"); m_ndigivsbx2D[i]->GetYaxis()->SetTitle("Number of Hits");
}
if(m_scalefact.size()>=2) {
sprintf(name,"n%sdigivsbx2Dzoom",slab.c_str());
m_ndigivsbx2Dzoom[i] =subev.make<TH2F>(name,title,3564,-0.5,3563.5,m_nbins,0,m_binmax[ui]/(m_scalefact[1]*m_nbins)*m_nbins);
m_ndigivsbx2Dzoom[i]->GetXaxis()->SetTitle("BX#"); m_ndigivsbx2Dzoom[i]->GetYaxis()->SetTitle("Number of Hits");
}
if(m_scalefact.size()>=3) {
sprintf(name,"n%sdigivsbx2Dzoom2",slab.c_str());
m_ndigivsbx2Dzoom2[i] =subev.make<TH2F>(name,title,3564,-0.5,3563.5,m_nbins,0,m_binmax[ui]/(m_scalefact[2]*m_nbins)*m_nbins);
m_ndigivsbx2Dzoom2[i]->GetXaxis()->SetTitle("BX#"); m_ndigivsbx2Dzoom2[i]->GetYaxis()->SetTitle("Number of Hits");
}
}
}
| void DigiBXCorrHistogramMaker< T >::fill | ( | const T & | he, |
| const std::map< int, int > & | ndigi, | ||
| const edm::Handle< APVCyclePhaseCollection > & | phase | ||
| ) |
Definition at line 337 of file DigiBXCorrHistogramMaker.h.
References edm::HandleBase::failedToGet(), HcalObjRepresent::Fill(), i, APVCyclePhaseCollection::invalid, edm::HandleBase::isValid(), APVCyclePhaseCollection::multiphase, APVCyclePhaseCollection::nopartition, and interactiveExample::ui.
Referenced by MultiplicityTimeCorrelations::analyze().
{
for(std::map<int,int>::const_iterator digi=ndigi.begin();digi!=ndigi.end();digi++) {
if(m_labels.find(digi->first) != m_labels.end()) {
const int i=digi->first; const unsigned int ui = i;
int thephase = APVCyclePhaseCollection::invalid;
if(m_phasepart.find(ui)!=m_phasepart.end() && m_phasepart[ui]!="None") {
if(!phase.failedToGet() && phase.isValid()) {
thephase = phase->getPhase(m_phasepart[ui]);
}
}
long long tbx = he.absoluteBX();
if(thephase!=APVCyclePhaseCollection::nopartition &&
thephase!=APVCyclePhaseCollection::multiphase &&
thephase!=APVCyclePhaseCollection::invalid) tbx -= thephase;
if(m_nmeandigivscycle.find(i)!=m_nmeandigivscycle.end()) m_nmeandigivscycle[i]->Fill(tbx%70,digi->second);
if(m_ndigivscycle.find(i)!=m_ndigivscycle.end()) m_ndigivscycle[i]->Fill(tbx%70,digi->second);
if(m_ndigivscyclezoom.find(i)!=m_ndigivscyclezoom.end()) m_ndigivscyclezoom[i]->Fill(tbx%70,digi->second);
if(m_ndigivscyclezoom2.find(i)!=m_ndigivscyclezoom2.end()) m_ndigivscyclezoom2[i]->Fill(tbx%70,digi->second);
if(m_runHisto) {
if(m_ndigivscycletime.find(i)!=m_ndigivscycletime.end()) {
if(m_ndigivscycletime[i]!=0 && (*m_ndigivscycletime[i])!=0 ) (*m_ndigivscycletime[i])->Fill(tbx%70,(int)he._orbit,digi->second);
}
}
m_ndigivsbx[i]->Fill(he.bx(),digi->second);
if(m_ndigivsbx2D.find(i)!=m_ndigivsbx2D.end()) m_ndigivsbx2D[i]->Fill(he.bx(),digi->second);
if(m_ndigivsbx2Dzoom.find(i)!=m_ndigivsbx2Dzoom.end()) m_ndigivsbx2Dzoom[i]->Fill(he.bx(),digi->second);
if(m_ndigivsbx2Dzoom2.find(i)!=m_ndigivsbx2Dzoom2.end()) m_ndigivsbx2Dzoom2[i]->Fill(he.bx(),digi->second);
if(he.depth()>0) {
long long dbx = he.deltaBX();
m_ndigivsdbx[i]->Fill(dbx,digi->second);
m_ndigivsdbxzoom[i]->Fill(dbx,digi->second);
m_ndigivsdbxzoom2[i]->Fill(dbx,digi->second);
if(m_ndigivsdbx2D.find(i)!=m_ndigivsdbx2D.end()) m_ndigivsdbx2D[i]->Fill(dbx,digi->second);
if(m_ndigivsdbx2Dzoom.find(i)!=m_ndigivsdbx2Dzoom.end()) m_ndigivsdbx2Dzoom[i]->Fill(dbx,digi->second);
if(m_ndigivsdbx2Dzoom2.find(i)!=m_ndigivsdbx2Dzoom2.end()) m_ndigivsdbx2Dzoom2[i]->Fill(dbx,digi->second);
if(thephase!=APVCyclePhaseCollection::nopartition &&
thephase!=APVCyclePhaseCollection::multiphase &&
thephase!=APVCyclePhaseCollection::invalid) {
long long dbxincycle = he.deltaBXinCycle(thephase);
if(m_ndigivsdbxincycle2D.find(i)!=m_ndigivsdbxincycle2D.end()) m_ndigivsdbxincycle2D[i]->Fill(dbxincycle,digi->second);
if(m_ndigivsdbxincycle.find(i)!=m_ndigivsdbxincycle.end()) m_ndigivsdbxincycle[i]->Fill(dbxincycle,digi->second);
}
long long prevtbx = he.absoluteBX(1);
if(thephase!=APVCyclePhaseCollection::nopartition &&
thephase!=APVCyclePhaseCollection::multiphase &&
thephase!=APVCyclePhaseCollection::invalid) prevtbx -= thephase;
if(m_ndigivscycledbx.find(i)!=m_ndigivscycledbx.end()) m_ndigivscycledbx[i]->Fill(prevtbx%70,dbx,digi->second);
if(m_ndigivscycle2dbx.find(i)!=m_ndigivscycle2dbx.end()) m_ndigivscycle2dbx[i]->Fill(tbx%70,dbx,digi->second);
if(he.depth()>1) {
long long dbx2 = he.deltaBX(2);
m_ndigivsdbx3zoom[i]->Fill(dbx2,dbx,digi->second);
if(thephase!=APVCyclePhaseCollection::nopartition &&
thephase!=APVCyclePhaseCollection::multiphase &&
thephase!=APVCyclePhaseCollection::invalid) {
long long dbxincycle = he.deltaBXinCycle(thephase);
long long dbxincycle2 = he.deltaBXinCycle(2,thephase);
if(m_dbx3Histo) {
if(m_ndigivsdbxincycle3.find(i)!=m_ndigivsdbxincycle3.end()) m_ndigivsdbxincycle3[i]->Fill(dbxincycle,dbxincycle2-dbxincycle,digi->second);
if(m_dbx3Histo3D) {
if(m_ndigivsdbxincycle33D.find(i)!=m_ndigivsdbxincycle33D.end()) m_ndigivsdbxincycle33D[i]->Fill(dbxincycle,dbxincycle2-dbxincycle,digi->second);
}
}
}
}
}
}
else {
edm::LogWarning("MissingKey") << " Key " << digi->first << " is missing ";
}
}
}
| void DigiBXCorrHistogramMaker< T >::fillcorr | ( | const T & | he1, |
| const T & | he2, | ||
| const std::map< int, int > & | ndigi | ||
| ) |
Definition at line 432 of file DigiBXCorrHistogramMaker.h.
References i.
{
for(std::map<int,int>::const_iterator digi=ndigi.begin();digi!=ndigi.end();digi++) {
if(m_labels.find(digi->first) != m_labels.end()) {
const int i=digi->first;
long long dbx = he2.deltaBX(he1);
m_digicorr[i]->Fill(dbx,digi->second);
}
else {
edm::LogWarning("MissingKey") << " Key " << digi->first << " is missing ";
}
}
}
std::map<unsigned int, int> DigiBXCorrHistogramMaker< T >::m_binmax [private] |
Definition at line 44 of file DigiBXCorrHistogramMaker.h.
Referenced by DigiBXCorrHistogramMaker< T >::DigiBXCorrHistogramMaker().
const bool DigiBXCorrHistogramMaker< T >::m_dbx3Histo [private] |
Definition at line 39 of file DigiBXCorrHistogramMaker.h.
const bool DigiBXCorrHistogramMaker< T >::m_dbx3Histo3D [private] |
Definition at line 40 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile*> DigiBXCorrHistogramMaker< T >::m_digicorr [private] |
Definition at line 84 of file DigiBXCorrHistogramMaker.h.
std::string DigiBXCorrHistogramMaker< T >::m_hitname [private] |
Definition at line 38 of file DigiBXCorrHistogramMaker.h.
std::map<int,std::string> DigiBXCorrHistogramMaker< T >::m_labels [private] |
Definition at line 43 of file DigiBXCorrHistogramMaker.h.
const int DigiBXCorrHistogramMaker< T >::m_nbins [private] |
Definition at line 47 of file DigiBXCorrHistogramMaker.h.
int DigiBXCorrHistogramMaker< T >::m_ncorbins [private] |
Definition at line 37 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile*> DigiBXCorrHistogramMaker< T >::m_ndigivsbx [private] |
Definition at line 64 of file DigiBXCorrHistogramMaker.h.
std::map<int,TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsbx2D [private] |
Definition at line 65 of file DigiBXCorrHistogramMaker.h.
std::map<int,TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsbx2Dzoom [private] |
Definition at line 66 of file DigiBXCorrHistogramMaker.h.
std::map<int,TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsbx2Dzoom2 [private] |
Definition at line 67 of file DigiBXCorrHistogramMaker.h.
std::map<int,TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivscycle [private] |
Definition at line 60 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile2D*> DigiBXCorrHistogramMaker< T >::m_ndigivscycle2dbx [private] |
Definition at line 72 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile2D*> DigiBXCorrHistogramMaker< T >::m_ndigivscycledbx [private] |
Definition at line 70 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile2D**> DigiBXCorrHistogramMaker< T >::m_ndigivscycletime [private] |
Definition at line 74 of file DigiBXCorrHistogramMaker.h.
std::map<int,TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivscyclezoom [private] |
Definition at line 61 of file DigiBXCorrHistogramMaker.h.
std::map<int,TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivscyclezoom2 [private] |
Definition at line 62 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbx [private] |
Definition at line 51 of file DigiBXCorrHistogramMaker.h.
std::map<int,TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbx2D [private] |
Definition at line 76 of file DigiBXCorrHistogramMaker.h.
std::map<int,TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbx2Dzoom [private] |
Definition at line 78 of file DigiBXCorrHistogramMaker.h.
std::map<int,TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbx2Dzoom2 [private] |
Definition at line 77 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile2D*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbx3zoom [private] |
Definition at line 80 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxincycle [private] |
Definition at line 55 of file DigiBXCorrHistogramMaker.h.
std::map<int,TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxincycle2D [private] |
Definition at line 56 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile2D*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxincycle3 [private] |
Definition at line 81 of file DigiBXCorrHistogramMaker.h.
std::map<int,TH3F*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxincycle33D [private] |
Definition at line 82 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxzoom [private] |
Definition at line 53 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxzoom2 [private] |
Definition at line 52 of file DigiBXCorrHistogramMaker.h.
std::map<int,TProfile*> DigiBXCorrHistogramMaker< T >::m_nmeandigivscycle [private] |
Definition at line 58 of file DigiBXCorrHistogramMaker.h.
std::map<int,std::string> DigiBXCorrHistogramMaker< T >::m_phasepart [private] |
Definition at line 45 of file DigiBXCorrHistogramMaker.h.
Referenced by DigiBXCorrHistogramMaker< T >::DigiBXCorrHistogramMaker().
RunHistogramManager DigiBXCorrHistogramMaker< T >::m_rhm [private] |
Definition at line 49 of file DigiBXCorrHistogramMaker.h.
const bool DigiBXCorrHistogramMaker< T >::m_runHisto [private] |
Definition at line 41 of file DigiBXCorrHistogramMaker.h.
std::vector<int> DigiBXCorrHistogramMaker< T >::m_scalefact [private] |
Definition at line 46 of file DigiBXCorrHistogramMaker.h.
1.7.3