|
|
|
#include <HcalPedestalsAnalysis.h>
Public Member Functions | |
| void | analyze (const edm::Event &event, const edm::EventSetup &eventSetup) |
| HcalPedestalsAnalysis (const edm::ParameterSet &ps) | |
| virtual | ~HcalPedestalsAnalysis () |
Private Attributes | |
| std::vector< NewPedBunch > | Bunches |
| bool | dumpXML |
| bool | firsttime |
| int | firstTS |
| TH1F * | HBMeans |
| TH1F * | HBWidths |
| TH1F * | HEMeans |
| TH1F * | HEWidths |
| TH1F * | HFMeans |
| TH1F * | HFWidths |
| bool | hiSaveFlag |
| TH1F * | HOMeans |
| TH1F * | HOWidths |
| int | lastTS |
| std::string | pedsADCfilename |
| std::string | pedsfCfilename |
| HcalPedestals * | rawPedsItem |
| HcalPedestals * | rawPedsItemfc |
| HcalPedestalWidths * | rawWidthsItem |
| HcalPedestalWidths * | rawWidthsItemfc |
| std::string | ROOTfilename |
| int | runnum |
| TFile * | theFile |
| bool | verboseflag |
| std::string | widthsADCfilename |
| std::string | widthsfCfilename |
| std::string | XMLfilename |
| std::string | XMLtag |
Definition at line 68 of file HcalPedestalsAnalysis.h.
| HcalPedestalsAnalysis::HcalPedestalsAnalysis | ( | const edm::ParameterSet & | ps | ) |
Definition at line 10 of file HcalPedestalsAnalysis.cc.
References dumpXML, firsttime, firstTS, edm::ParameterSet::getUntrackedParameter(), hiSaveFlag, lastTS, rawPedsItem, rawPedsItemfc, rawWidthsItem, rawWidthsItemfc, and verboseflag.
{
hiSaveFlag = ps.getUntrackedParameter<bool>("hiSaveFlag", false);
dumpXML = ps.getUntrackedParameter<bool>("dumpXML", true);
verboseflag = ps.getUntrackedParameter<bool>("verbose", false);
firstTS = ps.getUntrackedParameter<int>("firstTS", 0);
lastTS = ps.getUntrackedParameter<int>("lastTS", 9);
firsttime = true;
rawPedsItem = new HcalPedestals(true);
rawWidthsItem = new HcalPedestalWidths(true);
rawPedsItemfc = new HcalPedestals(false);
rawWidthsItemfc = new HcalPedestalWidths(false);
}
| HcalPedestalsAnalysis::~HcalPedestalsAnalysis | ( | ) | [virtual] |
Definition at line 26 of file HcalPedestalsAnalysis.cc.
References HcalCondObjectContainer< Item >::addValues(), Bunches, gather_cfg::cout, CastorDbASCIIIO::dumpObject(), dumpXML, HBMeans, HBWidths, HEMeans, HEWidths, HFMeans, HFWidths, hiSaveFlag, HOMeans, HOWidths, i, pedsADCfilename, pedsfCfilename, rawPedsItem, rawPedsItemfc, rawWidthsItem, rawWidthsItemfc, runnum, HcalPedestalWidth::setSigma(), theFile, verboseflag, widthsADCfilename, widthsfCfilename, XMLfilename, and XMLtag.
{
//Calculate pedestal constants
std::cout << "Calculating Pedestal constants...\n";
std::vector<NewPedBunch>::iterator bunch_it;
for(bunch_it=Bunches.begin(); bunch_it != Bunches.end(); bunch_it++)
{
if(bunch_it->usedflag){
if(verboseflag) std::cout << "Analyzing channel " << bunch_it->detid << std::endl;
//pedestal constant is the mean
if(bunch_it->num[0][0]!=0) bunch_it->cap[0] /= bunch_it->num[0][0];
if(bunch_it->num[1][1]!=0) bunch_it->cap[1] /= bunch_it->num[1][1];
if(bunch_it->num[2][2]!=0) bunch_it->cap[2] /= bunch_it->num[2][2];
if(bunch_it->num[3][3]!=0) bunch_it->cap[3] /= bunch_it->num[3][3];
if(bunch_it->num[0][0]!=0) bunch_it->capfc[0] /= bunch_it->num[0][0];
if(bunch_it->num[1][1]!=0) bunch_it->capfc[1] /= bunch_it->num[1][1];
if(bunch_it->num[2][2]!=0) bunch_it->capfc[2] /= bunch_it->num[2][2];
if(bunch_it->num[3][3]!=0) bunch_it->capfc[3] /= bunch_it->num[3][3];
//widths are the covariance matrix--assumed symmetric
bunch_it->sig[0][0] = (bunch_it->prod[0][0]/bunch_it->num[0][0])-(bunch_it->cap[0])*(bunch_it->cap[0]);
bunch_it->sig[1][1] = (bunch_it->prod[1][1]/bunch_it->num[1][1])-(bunch_it->cap[1])*(bunch_it->cap[1]);
bunch_it->sig[2][2] = (bunch_it->prod[2][2]/bunch_it->num[2][2])-(bunch_it->cap[2])*(bunch_it->cap[2]);
bunch_it->sig[3][3] = (bunch_it->prod[3][3]/bunch_it->num[3][3])-(bunch_it->cap[3])*(bunch_it->cap[3]);
bunch_it->sig[0][1] = (bunch_it->prod[0][1])/(bunch_it->num[0][1])-(bunch_it->cap[0]*bunch_it->cap[1]);
bunch_it->sig[0][2] = (bunch_it->prod[0][2])/(bunch_it->num[0][2])-(bunch_it->cap[0]*bunch_it->cap[2]);
bunch_it->sig[0][3] = (bunch_it->prod[3][0])/(bunch_it->num[3][0])-(bunch_it->cap[0]*bunch_it->cap[3]); // sig03 MISNAMED in object!
bunch_it->sig[1][2] = (bunch_it->prod[1][2])/(bunch_it->num[1][2])-(bunch_it->cap[1]*bunch_it->cap[2]);
bunch_it->sig[1][3] = (bunch_it->prod[1][3])/(bunch_it->num[1][3])-(bunch_it->cap[1]*bunch_it->cap[3]);
bunch_it->sig[2][3] = (bunch_it->prod[2][3])/(bunch_it->num[2][3])-(bunch_it->cap[2]*bunch_it->cap[3]);
bunch_it->sigfc[0][0] = (bunch_it->prodfc[0][0]/bunch_it->num[0][0])-(bunch_it->capfc[0])*(bunch_it->capfc[0]);
bunch_it->sigfc[1][1] = (bunch_it->prodfc[1][1]/bunch_it->num[1][1])-(bunch_it->capfc[1])*(bunch_it->capfc[1]);
bunch_it->sigfc[2][2] = (bunch_it->prodfc[2][2]/bunch_it->num[2][2])-(bunch_it->capfc[2])*(bunch_it->capfc[2]);
bunch_it->sigfc[3][3] = (bunch_it->prodfc[3][3]/bunch_it->num[3][3])-(bunch_it->capfc[3])*(bunch_it->capfc[3]);
bunch_it->sigfc[0][1] = (bunch_it->prodfc[0][1]/(bunch_it->num[0][1]))-(bunch_it->capfc[0]*bunch_it->capfc[1]);
bunch_it->sigfc[0][2] = (bunch_it->prodfc[0][2]/(bunch_it->num[0][2]))-(bunch_it->capfc[0]*bunch_it->capfc[2]);
bunch_it->sigfc[0][3] = (bunch_it->prodfc[3][0]/(bunch_it->num[3][0]))-(bunch_it->capfc[0]*bunch_it->capfc[3]); //sig03 MISNAMED in object!
bunch_it->sigfc[1][2] = (bunch_it->prodfc[1][2]/(bunch_it->num[1][2]))-(bunch_it->capfc[1]*bunch_it->capfc[2]);
bunch_it->sigfc[1][3] = (bunch_it->prodfc[1][3]/(bunch_it->num[1][3]))-(bunch_it->capfc[1]*bunch_it->capfc[3]);
bunch_it->sigfc[2][3] = (bunch_it->prodfc[2][3]/(bunch_it->num[2][3]))-(bunch_it->capfc[2]*bunch_it->capfc[3]);
if(bunch_it->detid.subdet() == 1){
for(int i = 0; i != 4; i++){
HBMeans->Fill(bunch_it->cap[i]);
HBWidths->Fill(bunch_it->sig[i][i]);
}
}
if(bunch_it->detid.subdet() == 2){
for(int i = 0; i != 4; i++){
HEMeans->Fill(bunch_it->cap[i]);
HEWidths->Fill(bunch_it->sig[i][i]);
}
}
if(bunch_it->detid.subdet() == 3){
for(int i = 0; i != 4; i++){
HOMeans->Fill(bunch_it->cap[i]);
HOWidths->Fill(bunch_it->sig[i][i]);
}
}
if(bunch_it->detid.subdet() == 4){
for(int i = 0; i != 4; i++){
HFMeans->Fill(bunch_it->cap[i]);
HFWidths->Fill(bunch_it->sig[i][i]);
}
}
const HcalPedestal item(bunch_it->detid, bunch_it->cap[0], bunch_it->cap[1], bunch_it->cap[2], bunch_it->cap[3]);
rawPedsItem->addValues(item);
HcalPedestalWidth widthsp(bunch_it->detid);
widthsp.setSigma(0,0,bunch_it->sig[0][0]);
widthsp.setSigma(0,1,bunch_it->sig[0][1]);
widthsp.setSigma(0,2,bunch_it->sig[0][2]);
widthsp.setSigma(0,3,bunch_it->sig[0][3]);
widthsp.setSigma(1,1,bunch_it->sig[1][1]);
widthsp.setSigma(1,2,bunch_it->sig[1][2]);
widthsp.setSigma(1,3,bunch_it->sig[1][3]);
widthsp.setSigma(2,2,bunch_it->sig[2][2]);
widthsp.setSigma(2,3,bunch_it->sig[2][3]);
widthsp.setSigma(3,3,bunch_it->sig[3][3]);
rawWidthsItem->addValues(widthsp);
const HcalPedestal itemfc(bunch_it->detid, bunch_it->capfc[0], bunch_it->capfc[1],
bunch_it->capfc[2], bunch_it->capfc[3]);
rawPedsItemfc->addValues(itemfc);
HcalPedestalWidth widthspfc(bunch_it->detid);
widthspfc.setSigma(0,0,bunch_it->sigfc[0][0]);
widthspfc.setSigma(0,1,bunch_it->sigfc[0][1]);
widthspfc.setSigma(0,2,bunch_it->sigfc[0][2]);
widthspfc.setSigma(0,3,bunch_it->sigfc[0][3]);
widthspfc.setSigma(1,1,bunch_it->sigfc[1][1]);
widthspfc.setSigma(1,2,bunch_it->sigfc[1][2]);
widthspfc.setSigma(1,3,bunch_it->sigfc[1][3]);
widthspfc.setSigma(2,2,bunch_it->sigfc[2][2]);
widthspfc.setSigma(2,3,bunch_it->sigfc[2][3]);
widthspfc.setSigma(3,3,bunch_it->sigfc[3][3]);
rawWidthsItemfc->addValues(widthspfc);
}
}
// dump the resulting list of pedestals into a file
std::ofstream outStream1(pedsADCfilename.c_str());
HcalDbASCIIIO::dumpObject (outStream1, (*rawPedsItem) );
std::ofstream outStream2(widthsADCfilename.c_str());
HcalDbASCIIIO::dumpObject (outStream2, (*rawWidthsItem) );
std::ofstream outStream3(pedsfCfilename.c_str());
HcalDbASCIIIO::dumpObject (outStream3, (*rawPedsItemfc) );
std::ofstream outStream4(widthsfCfilename.c_str());
HcalDbASCIIIO::dumpObject (outStream4, (*rawWidthsItemfc) );
if(dumpXML){
std::ofstream outStream5(XMLfilename.c_str());
HcalDbXml::dumpObject (outStream5, runnum, 0, 2147483647, XMLtag, 1, (*rawPedsItem), (*rawWidthsItem));
}
if(hiSaveFlag){
theFile->Write();
}else{
theFile->cd();
theFile->cd("HB");
HBMeans->Write();
HBWidths->Write();
theFile->cd();
theFile->cd("HF");
HFMeans->Write();
HFWidths->Write();
theFile->cd();
theFile->cd("HE");
HEMeans->Write();
HEWidths->Write();
theFile->cd();
theFile->cd("HO");
HOMeans->Write();
HOWidths->Write();
}
std::cout << "Writing ROOT file... ";
theFile->Close();
std::cout << "ROOT file closed.\n";
delete rawPedsItem;
delete rawWidthsItem;
delete rawPedsItemfc;
delete rawWidthsItemfc;
}
| void HcalPedestalsAnalysis::analyze | ( | const edm::Event & | event, |
| const edm::EventSetup & | eventSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 175 of file HcalPedestalsAnalysis.cc.
References a, HcalQIESample::adc(), HcalElectronicsMap::allPrecisionId(), Bunches, NewPedBunch::cap, NewPedBunch::capfc, HcalQIESample::capid(), HcalQIECoder::charge(), NewPedBunch::detid, firsttime, firstTS, edm::EventSetup::get(), edm::Event::getByType(), HBMeans, HBWidths, HEMeans, HEWidths, HFMeans, HFWidths, HOMeans, HOWidths, i, HODataFrame::id(), HBHEDataFrame::id(), edm::EventBase::id(), HFDataFrame::id(), j, lastTS, NewPedBunch::num, pedsADCfilename, pedsfCfilename, NewPedBunch::prod, NewPedBunch::prodfc, edm::ESHandle< T >::product(), DetId::rawId(), ROOTfilename, edm::EventID::run(), runnum, HFDataFrame::sample(), HBHEDataFrame::sample(), HODataFrame::sample(), NewPedBunch::sig, NewPedBunch::sigfc, HODataFrame::size(), HBHEDataFrame::size(), HFDataFrame::size(), theFile, NewPedBunch::usedflag, widthsADCfilename, widthsfCfilename, XMLfilename, and XMLtag.
{
using namespace edm;
using namespace std;
edm::Handle<HBHEDigiCollection> hbhe; e.getByType(hbhe);
edm::Handle<HODigiCollection> ho; e.getByType(ho);
edm::Handle<HFDigiCollection> hf; e.getByType(hf);
edm::ESHandle<HcalDbService> conditions;
iSetup.get<HcalDbRecord>().get(conditions);
const HcalQIEShape* shape = conditions->getHcalShape();
if(firsttime)
{
runnum = e.id().run();
std::string runnum_string;
std::stringstream tempstringout;
tempstringout << runnum;
runnum_string = tempstringout.str();
ROOTfilename = runnum_string + "-peds_ADC.root";
pedsADCfilename = runnum_string + "-peds_ADC.txt";
pedsfCfilename = runnum_string + "-peds_fC.txt";
widthsADCfilename = runnum_string + "-widths_ADC.txt";
widthsfCfilename = runnum_string + "-widths_fC.txt";
XMLfilename = runnum_string + "-peds_ADC_complete.xml";
XMLtag = "Hcal_pedestals_" + runnum_string;
theFile = new TFile(ROOTfilename.c_str(), "RECREATE");
theFile->cd();
// Create sub-directories
theFile->mkdir("HB");
theFile->mkdir("HE");
theFile->mkdir("HF");
theFile->mkdir("HO");
theFile->cd();
HBMeans = new TH1F("All Ped Means HB","All Ped Means HB", 100, 0, 9);
HBWidths = new TH1F("All Ped Widths HB","All Ped Widths HB", 100, 0, 3);
HEMeans = new TH1F("All Ped Means HE","All Ped Means HE", 100, 0, 9);
HEWidths = new TH1F("All Ped Widths HE","All Ped Widths HE", 100, 0, 3);
HFMeans = new TH1F("All Ped Means HF","All Ped Means HF", 100, 0, 9);
HFWidths = new TH1F("All Ped Widths HF","All Ped Widths HF", 100, 0, 3);
HOMeans = new TH1F("All Ped Means HO","All Ped Means HO", 100, 0, 9);
HOWidths = new TH1F("All Ped Widths HO","All Ped Widths HO", 100, 0, 3);
edm::ESHandle<HcalElectronicsMap> refEMap;
iSetup.get<HcalElectronicsMapRcd>().get(refEMap);
const HcalElectronicsMap* myRefEMap = refEMap.product();
std::vector<HcalGenericDetId> listEMap = myRefEMap->allPrecisionId();
for (std::vector<HcalGenericDetId>::const_iterator it = listEMap.begin(); it != listEMap.end(); it++)
{
HcalGenericDetId mygenid(it->rawId());
if(mygenid.isHcalDetId())
{
NewPedBunch a;
HcalDetId chanid(mygenid.rawId());
a.detid = chanid;
a.usedflag = false;
for(int i = 0; i != 4; i++)
{
a.cap[i] = 0;
a.capfc[i] = 0;
for(int j = 0; j != 4; j++)
{
a.sig[i][j] = 0;
a.sigfc[i][j] = 0;
a.prod[i][j] = 0;
a.prodfc[i][j] = 0;
a.num[i][j] = 0;
}
}
Bunches.push_back(a);
}
}
firsttime = false;
}
std::vector<NewPedBunch>::iterator bunch_it;
for(HBHEDigiCollection::const_iterator j = hbhe->begin(); j != hbhe->end(); j++)
{
const HBHEDataFrame digi = (const HBHEDataFrame)(*j);
for(bunch_it = Bunches.begin(); bunch_it != Bunches.end(); bunch_it++)
if(bunch_it->detid.rawId() == digi.id().rawId()) break;
bunch_it->usedflag = true;
for(int ts = firstTS; ts != lastTS+1; ts++)
{
if(digi.sample(ts).adc() > 15) continue;
const HcalQIECoder* coder = conditions->getHcalCoder(digi.id().rawId());
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts).capid()] += 1;
bunch_it->cap[digi.sample(ts).capid()] += digi.sample(ts).adc();
double charge1 = coder->charge(*shape, digi.sample(ts).adc(), digi.sample(ts).capid());
bunch_it->capfc[digi.sample(ts).capid()] += charge1;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts).capid()] += (digi.sample(ts).adc() * digi.sample(ts).adc());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts).capid()] += charge1 * charge1;
if((ts+1 < digi.size()) && (ts+1 < lastTS)){
if(digi.sample(ts+1).adc() > 15) continue;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += digi.sample(ts).adc()*digi.sample(ts+1).adc();
double charge2 = coder->charge(*shape, digi.sample(ts+1).adc(), digi.sample(ts+1).capid());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += charge1*charge2;
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += 1;
}
if((ts+2 < digi.size()) && (ts+2 < lastTS)){
if(digi.sample(ts+2).adc() > 15) continue;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += digi.sample(ts).adc()*digi.sample(ts+2).adc();
double charge2 = coder->charge(*shape, digi.sample(ts+2).adc(), digi.sample(ts+2).capid());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += charge1*charge2;
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += 1;
}
if((ts+3 < digi.size()) && (ts+3 < lastTS)){
if(digi.sample(ts+3).adc() > 15) continue;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += digi.sample(ts).adc()*digi.sample(ts+3).adc();
double charge2 = coder->charge(*shape, digi.sample(ts+3).adc(), digi.sample(ts+3).capid());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += charge1*charge2;
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += 1;
}
}
}
for(HODigiCollection::const_iterator j = ho->begin(); j != ho->end(); j++)
{
const HODataFrame digi = (const HODataFrame)(*j);
for(bunch_it = Bunches.begin(); bunch_it != Bunches.end(); bunch_it++)
if(bunch_it->detid.rawId() == digi.id().rawId()) break;
bunch_it->usedflag = true;
for(int ts = firstTS; ts <= lastTS; ts++)
{
if(digi.sample(ts).adc() > 15) continue;
const HcalQIECoder* coder = conditions->getHcalCoder(digi.id().rawId());
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts).capid()] += 1;
bunch_it->cap[digi.sample(ts).capid()] += digi.sample(ts).adc();
double charge1 = coder->charge(*shape, digi.sample(ts).adc(), digi.sample(ts).capid());
bunch_it->capfc[digi.sample(ts).capid()] += charge1;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts).capid()] += (digi.sample(ts).adc() * digi.sample(ts).adc());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts).capid()] += charge1 * charge1;
if((ts+1 < digi.size()) && (ts+1 < lastTS)){
if(digi.sample(ts+1).adc() > 15) continue;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += digi.sample(ts).adc()*digi.sample(ts+1).adc();
double charge2 = coder->charge(*shape, digi.sample(ts+1).adc(), digi.sample(ts+1).capid());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += charge1*charge2;
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += 1;
}
if((ts+2 < digi.size()) && (ts+2 < lastTS)){
if(digi.sample(ts+2).adc() > 15) continue;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += digi.sample(ts).adc()*digi.sample(ts+2).adc();
double charge2 = coder->charge(*shape, digi.sample(ts+2).adc(), digi.sample(ts+2).capid());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += charge1*charge2;
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += 1;
}
if((ts+3 < digi.size()) && (ts+3 < lastTS)){
if(digi.sample(ts+3).adc() > 15) continue;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += digi.sample(ts).adc()*digi.sample(ts+3).adc();
double charge2 = coder->charge(*shape, digi.sample(ts+3).adc(), digi.sample(ts+3).capid());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += charge1*charge2;
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += 1;
}
}
}
for(HFDigiCollection::const_iterator j = hf->begin(); j != hf->end(); j++)
{
const HFDataFrame digi = (const HFDataFrame)(*j);
for(bunch_it = Bunches.begin(); bunch_it != Bunches.end(); bunch_it++)
if(bunch_it->detid.rawId() == digi.id().rawId()) break;
bunch_it->usedflag = true;
for(int ts = firstTS; ts <= lastTS; ts++)
{
if(digi.sample(ts).adc() > 15) continue;
const HcalQIECoder* coder = conditions->getHcalCoder(digi.id().rawId());
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts).capid()] += 1;
bunch_it->cap[digi.sample(ts).capid()] += digi.sample(ts).adc();
double charge1 = coder->charge(*shape, digi.sample(ts).adc(), digi.sample(ts).capid());
bunch_it->capfc[digi.sample(ts).capid()] += charge1;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts).capid()] += (digi.sample(ts).adc() * digi.sample(ts).adc());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts).capid()] += charge1 * charge1;
if((ts+1 < digi.size()) && (ts+1 < lastTS)){
if(digi.sample(ts+1).adc() > 15) continue;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += digi.sample(ts).adc()*digi.sample(ts+1).adc();
double charge2 = coder->charge(*shape, digi.sample(ts+1).adc(), digi.sample(ts+1).capid());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += charge1*charge2;
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += 1;
}
if((ts+2 < digi.size()) && (ts+2 < lastTS)){
if(digi.sample(ts+2).adc() > 15) continue;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += digi.sample(ts).adc()*digi.sample(ts+2).adc();
double charge2 = coder->charge(*shape, digi.sample(ts+2).adc(), digi.sample(ts+2).capid());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += charge1*charge2;
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += 1;
}
if((ts+3 < digi.size()) && (ts+3 < lastTS)){
if(digi.sample(ts+3).adc() > 15) continue;
bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += digi.sample(ts).adc()*digi.sample(ts+3).adc();
double charge2 = coder->charge(*shape, digi.sample(ts+3).adc(), digi.sample(ts+3).capid());
bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += charge1*charge2;
bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += 1;
}
}
}
//this is the last brace
}
std::vector<NewPedBunch> HcalPedestalsAnalysis::Bunches [private] |
Definition at line 80 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
bool HcalPedestalsAnalysis::dumpXML [private] |
Definition at line 83 of file HcalPedestalsAnalysis.h.
Referenced by HcalPedestalsAnalysis(), and ~HcalPedestalsAnalysis().
bool HcalPedestalsAnalysis::firsttime [private] |
Definition at line 106 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and HcalPedestalsAnalysis().
int HcalPedestalsAnalysis::firstTS [private] |
Definition at line 86 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and HcalPedestalsAnalysis().
TH1F* HcalPedestalsAnalysis::HBMeans [private] |
Definition at line 96 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
TH1F* HcalPedestalsAnalysis::HBWidths [private] |
Definition at line 97 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
TH1F* HcalPedestalsAnalysis::HEMeans [private] |
Definition at line 98 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
TH1F* HcalPedestalsAnalysis::HEWidths [private] |
Definition at line 99 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
TH1F* HcalPedestalsAnalysis::HFMeans [private] |
Definition at line 100 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
TH1F* HcalPedestalsAnalysis::HFWidths [private] |
Definition at line 101 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
bool HcalPedestalsAnalysis::hiSaveFlag [private] |
Definition at line 82 of file HcalPedestalsAnalysis.h.
Referenced by HcalPedestalsAnalysis(), and ~HcalPedestalsAnalysis().
TH1F* HcalPedestalsAnalysis::HOMeans [private] |
Definition at line 102 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
TH1F* HcalPedestalsAnalysis::HOWidths [private] |
Definition at line 103 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
int HcalPedestalsAnalysis::lastTS [private] |
Definition at line 87 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and HcalPedestalsAnalysis().
std::string HcalPedestalsAnalysis::pedsADCfilename [private] |
Definition at line 89 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
std::string HcalPedestalsAnalysis::pedsfCfilename [private] |
Definition at line 90 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
HcalPedestals* HcalPedestalsAnalysis::rawPedsItem [private] |
Definition at line 107 of file HcalPedestalsAnalysis.h.
Referenced by HcalPedestalsAnalysis(), and ~HcalPedestalsAnalysis().
Definition at line 109 of file HcalPedestalsAnalysis.h.
Referenced by HcalPedestalsAnalysis(), and ~HcalPedestalsAnalysis().
Definition at line 108 of file HcalPedestalsAnalysis.h.
Referenced by HcalPedestalsAnalysis(), and ~HcalPedestalsAnalysis().
Definition at line 110 of file HcalPedestalsAnalysis.h.
Referenced by HcalPedestalsAnalysis(), and ~HcalPedestalsAnalysis().
std::string HcalPedestalsAnalysis::ROOTfilename [private] |
Definition at line 88 of file HcalPedestalsAnalysis.h.
Referenced by analyze().
int HcalPedestalsAnalysis::runnum [private] |
Definition at line 85 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
TFile* HcalPedestalsAnalysis::theFile [private] |
Definition at line 105 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
bool HcalPedestalsAnalysis::verboseflag [private] |
Definition at line 84 of file HcalPedestalsAnalysis.h.
Referenced by HcalPedestalsAnalysis(), and ~HcalPedestalsAnalysis().
std::string HcalPedestalsAnalysis::widthsADCfilename [private] |
Definition at line 91 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
std::string HcalPedestalsAnalysis::widthsfCfilename [private] |
Definition at line 92 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
std::string HcalPedestalsAnalysis::XMLfilename [private] |
Definition at line 93 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
std::string HcalPedestalsAnalysis::XMLtag [private] |
Definition at line 94 of file HcalPedestalsAnalysis.h.
Referenced by analyze(), and ~HcalPedestalsAnalysis().
1.7.3