|
|
|
#include <HFLightCalRand.h>
Public Member Functions | |
| virtual void | analyze (const edm::Event &fEvent, const edm::EventSetup &fSetup) |
| virtual void | beginJob () |
| virtual void | endJob (void) |
| HFLightCalRand (const edm::ParameterSet &fConfiguration) | |
| virtual | ~HFLightCalRand () |
Private Attributes | |
| edm::InputTag | hcalCalibDigiCollectionTag_ |
| edm::InputTag | hfDigiCollectionTag_ |
| std::string | histfile |
| TH2F * | hnpemapM |
| TH2F * | hnpemapP |
| TH1F * | hnpevar |
| TH1F * | hped [26][36][2] |
| TH1F * | hpedmean |
| TH1F * | hpedpin [8][3] |
| TH1F * | hpedrms |
| TH2F * | hsignalmapM |
| TH2F * | hsignalmapP |
| TH1F * | hsignalmean |
| TH1F * | hsignalrms |
| TH2F * | hsignalRMSmapM |
| TH2F * | hsignalRMSmapP |
| TH1F * | hsp [26][36][2] |
| TH1F * | hspe [26][36][2] |
| TH1F * | hspepin [8][3] |
| TH1F * | hspes |
| TH1F * | hsppin [8][3] |
| TH1F * | htmax |
| TH1F * | htmean |
| TH1F * | hts [26][36][2] |
| TH1F * | htsm [26][36][2] |
| TH1F * | htsmpin [8][3] |
| TH1F * | htspin [8][3] |
| TFile * | mFile |
| FILE * | preFile |
| std::string | prefile |
| std::string | textfile |
| FILE * | tFile |
Definition at line 12 of file HFLightCalRand.h.
| HFLightCalRand::HFLightCalRand | ( | const edm::ParameterSet & | fConfiguration | ) |
Definition at line 43 of file HFLightCalRand.cc.
References edm::ParameterSet::getUntrackedParameter(), histfile, and textfile.
: hfDigiCollectionTag_(fConfiguration.getParameter<edm::InputTag>("hfDigiCollectionTag")), hcalCalibDigiCollectionTag_(fConfiguration.getParameter<edm::InputTag>("hcalCalibDigiCollectionTag")) { //std::string histfile = fConfiguration.getUntrackedParameter<string>("rootFile"); histfile = fConfiguration.getUntrackedParameter<string>("rootFile"); textfile = fConfiguration.getUntrackedParameter<string>("textFile"); }
| HFLightCalRand::~HFLightCalRand | ( | ) | [virtual] |
Definition at line 52 of file HFLightCalRand.cc.
{
//delete mFile;
}
| void HFLightCalRand::analyze | ( | const edm::Event & | fEvent, |
| const edm::EventSetup & | fSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 373 of file HFLightCalRand.cc.
References ecalMGPA::adc(), calib, gather_cfg::cout, HcalDetId::depth(), edm::EventID::event(), EventN, HcalObjRepresent::Fill(), edm::Event::getByLabel(), hcalCalibDigiCollectionTag_, HcalForward, hfDigiCollectionTag_, hped, hsp, hspe, htmax, htmean, hts, htsm, edm::EventBase::id(), HFDataFrame::id(), HcalDetId::ieta(), HcalDetId::iphi(), edm::EventID::run(), runNumb, convertSQLiteXML::runNumber, HFDataFrame::size(), and HcalDetId::subdet().
{
// event ID
edm::EventID eventId = fEvent.id();
int runNumber = eventId.run ();
int eventNumber = eventId.event ();
if (runNumb==0) runNumb=runNumber;
EventN++;
if (verbose) std::cout << "========================================="<<std::endl
<< "run/event: "<<runNumber<<'/'<<eventNumber<<std::endl;
Double_t buf[20];
Double_t maxADC,signal,ped=0,meant;
Int_t maxisample=0,i1=3,i2=6;
// HF PIN-diodes
edm::Handle<HcalCalibDigiCollection> calib;
fEvent.getByLabel(hcalCalibDigiCollectionTag_, calib);
if (verbose) std::cout<<"Analysis-> total CAL digis= "<<calib->size()<<std::endl;
/* COMMENTED OUT by J. Mans (7-28-2008) as major changes needed with new Calib DetId
re-commented out by R.Ofierzynski (11.Nov.2008) - to be able to provide a consistent code for CMSSW_3_0_0_pre3:
major changes are needed for the new Calib DetId which does not have the old methods any more
for (unsigned j = 0; j < calib->size (); ++j) {
const HcalCalibDataFrame digi = (*calib)[j];
HcalElectronicsId elecId = digi.elecId();
HcalCalibDetId calibId = digi.id();
if (verbose) std::cout<<calibId.sectorString().c_str()<<" "<<calibId.rbx()<<" "<<elecId.fiberChanId()<<std::endl;
int isector = calibId.rbx()-1;
int ipin = elecId.fiberChanId();
int iside = -1;
if (calibId.sectorString() == "HFP") iside = 0;
else if (calibId.sectorString() == "HFM") iside = 4;
if (iside != -1) {
maxADC=-99;
for (int isample = 0; isample < digi.size(); ++isample) {
int adc = digi[isample].adc();
int capid = digi[isample].capid ();
double linear_ADC = digi[isample].nominal_fC();
if (verbose) std::cout<<"PIN linear_ADC = "<<linear_ADC<<std::endl;
htspin[isector+iside][ipin]->Fill(isample,linear_ADC);
buf[isample]=linear_ADC;
if (maxADC<linear_ADC) {
maxADC=linear_ADC;
maxisample=isample;
}
}
i1=maxisample-1;
i2=maxisample+2;
if (i1<0) {i1=0;i2=3;}
else if (i2>9) {i1=6;i2=9;}
if (i1==0) ped=buf[8]+buf[9]+buf[6]+buf[7];
else if (i1==1) ped=buf[8]+buf[9]+buf[6]+buf[7];
else if (i1==2) ped=buf[0]+buf[1]+buf[6]+buf[7];
else if (i1==3) ped=buf[0]+buf[1]+buf[2]+buf[7];
else if (i1>=4) ped=buf[0]+buf[1]+buf[2]+buf[3];
signal=0;
for (ii=0;ii<4;ii++) signal+=TMath::Max(buf[ii+i1],ped/4.0);
hsppin[isector+iside][ipin]->Fill(signal);
hspepin[isector+iside][ipin]->Fill(signal);
hpedpin[isector+iside][ipin]->Fill(ped);
// Mean PIN signal time estimation
ped=ped/4;
meant=0;
for (ii=0;ii<4;ii++) meant+=(TMath::Max(buf[ii+i1],ped)-ped)*(ii+i1);
if (signal-ped*4>0) meant/=(signal-ped*4);
else meant=i1+1;
htsmpin[isector+iside][ipin]->Fill(meant);
}
}
*/
// HF
edm::Handle<HFDigiCollection> hf_digi;
fEvent.getByLabel(hfDigiCollectionTag_, hf_digi);
if (verbose) std::cout<<"Analysis-> total HF digis= "<<hf_digi->size()<<std::endl;
for (unsigned ihit = 0; ihit < hf_digi->size (); ++ihit) {
const HFDataFrame& frame = (*hf_digi)[ihit];
HcalDetId detId = frame.id();
int ieta = detId.ieta();
int iphi = detId.iphi();
int depth = detId.depth();
if (verbose) std::cout <<"HF digi # " <<ihit<<": eta/phi/depth: "
<<ieta<<'/'<<iphi<<'/'<< depth << std::endl;
if (ieta>0) ieta = ieta-29;
else ieta = 13-ieta-29;
for (int ii=0; ii<10; ii++) buf[ii]=0;
maxADC=-99;
for (int isample = 0; isample < frame.size(); ++isample) {
int adc = frame[isample].adc();
int capid = frame[isample].capid ();
double linear_ADC = frame[isample].nominal_fC();
double nominal_fC = detId.subdet () == HcalForward ? 2.6 * linear_ADC : linear_ADC;
if (verbose) std::cout << "Analysis-> HF sample # " << isample
<< ", capid=" << capid
<< ": ADC=" << adc
<< ", linearized ADC=" << linear_ADC
<< ", nominal fC=" << nominal_fC << std::endl;
hts[ieta][(iphi-1)/2][depth-1]->Fill(isample,linear_ADC);
buf[isample]=linear_ADC;
if (maxADC<linear_ADC) {
maxADC=linear_ADC;
maxisample=isample;
}
}
// Signal is four capIDs around maxTS, Pedestal is four capID off the signal
maxADC=-99;
for (int ibf=0; ibf<frame.size()-1; ibf++) {
Double_t sumbuf=0;
for (int jbf=0; jbf<2; jbf++) {
sumbuf += buf[ibf+jbf];
if (ibf+jbf<2) sumbuf -= (buf[ibf+jbf+4]+buf[ibf+jbf+8])/2.0;
else if (ibf+jbf<4) sumbuf -= buf[ibf+jbf+4];
else if (ibf+jbf<6) sumbuf -= (buf[ibf+jbf-4]+buf[ibf+jbf+4])/2.0;
else if (ibf+jbf<8) sumbuf -= buf[ibf+jbf-4];
else if (ibf+jbf<10) sumbuf -= (buf[ibf+jbf-4]+buf[ibf+jbf-8])/2.0;
}
if (sumbuf>maxADC) {
maxADC=sumbuf;
if (buf[ibf]>buf[ibf+1]) maxisample=ibf;
else maxisample=ibf+1;
//maxisample=ibf;
}
}
htmax->Fill(maxisample);
i1=maxisample-1;
i2=maxisample+2;
if (i1<0) {i1=0;i2=3;}
else if (i2>9) {i1=6;i2=9;}
else if (i2<9) {
if (buf[i1]<buf[i2+1]) {i1=i1+1;i2=i2+1;}
}
signal=buf[i1]+buf[i1+1]+buf[i1+2]+buf[i1+3];
hsp[ieta][(iphi-1)/2][depth-1]->Fill(signal);
hspe[ieta][(iphi-1)/2][depth-1]->Fill(signal);
if (i1==0) ped=(buf[4]+buf[8])/2.0+(buf[5]+buf[9])/2.0+buf[6]+buf[7];
else if (i1==1) ped=(buf[0]+buf[8])/2.0+(buf[5]+buf[9])/2.0+buf[6]+buf[7];
else if (i1==2) ped=(buf[0]+buf[8])/2.0+(buf[1]+buf[9])/2.0+buf[6]+buf[7];
else if (i1==3) ped=(buf[0]+buf[8])/2.0+(buf[1]+buf[9])/2.0+buf[2]+buf[7];
else if (i1==4) ped=(buf[0]+buf[8])/2.0+(buf[1]+buf[9])/2.0+buf[2]+buf[3];
else if (i1==5) ped=(buf[0]+buf[4])/2.0+(buf[1]+buf[9])/2.0+buf[2]+buf[3];
else if (i1==6) ped=(buf[0]+buf[4])/2.0+(buf[1]+buf[5])/2.0+buf[2]+buf[3];
/*
if (i1==0) ped=buf[8]+buf[9]+buf[6]+buf[7];
else if (i1==1) ped=buf[0]+buf[9]+buf[6]+buf[7];
else if (i1==2) ped=buf[0]+buf[1]+buf[6]+buf[7];
else if (i1==3) ped=buf[0]+buf[1]+buf[2]+buf[7];
else if (i1>=4) ped=buf[0]+buf[1]+buf[2]+buf[3];
*/
hped[ieta][(iphi-1)/2][depth-1]->Fill(ped);
// Mean signal time estimation
ped=ped/4;
meant=(buf[i1]-ped)*i1+(buf[i1+1]-ped)*(i1+1)+(buf[i1+2]-ped)*(i1+2)+(buf[i1+3]-ped)*(i1+3);
meant /= (buf[i1]-ped)+(buf[i1+1]-ped)+(buf[i1+2]-ped)+(buf[i1+3]-ped);
htmean->Fill(meant);
htsm[ieta][(iphi-1)/2][depth-1]->Fill(meant);
}
}
| void HFLightCalRand::beginJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 56 of file HFLightCalRand.cc.
References gather_cfg::cout, histfile, hnpemapM, hnpemapP, hnpevar, hped, hpedmean, hpedpin, hpedrms, hsignalmapM, hsignalmapP, hsignalmean, hsignalrms, hsignalRMSmapM, hsignalRMSmapP, hsp, hspe, hspepin, hspes, hsppin, htmax, htmean, hts, htsm, htsmpin, htspin, i, j, gen::k, mFile, NULL, textfile, and tFile.
{
char htit[64];
mFile = new TFile (histfile.c_str(),"RECREATE");
if ((tFile = fopen(textfile.c_str(),"w"))==NULL)printf("\nNo textfile open\n\n");
// General Histos
htmax = new TH1F("htmax","Max TS",10,-0.5,9.5);
htmean = new TH1F("htmean","Mean signal TS",100,0,10);
hsignalmean = new TH1F("hsignalmean","Mean ADC 4maxTS",1201,-25,30000);
hsignalrms = new TH1F("hsignalrms","RMS ADC 4maxTS",500,0,500);
hpedmean = new TH1F("hpedmean","Mean ADC 4lowTS",200,-10,90);
hpedrms = new TH1F("hpedrms","RMS ADC 4lowTS",200,0,100);
hspes = new TH1F("hspes","SPE if measured",200,0,40);
hnpevar = new TH1F("hnpevar","~N PE input",500,0,500);
hsignalmapP = new TH2F("hsignalmapP","Mean(Response) - Mean(Pedestal) HFP;#eta;#phi",26,28.5,41.5,36,0,72);
hsignalRMSmapP = new TH2F("hsignalRMSmapP","RMS Response HFP;#eta;#phi",26,28.5,41.5,36,0,72);
hnpemapP = new TH2F("hnpemapP","~N PE input HFP;#eta;#phi",26,28.5,41.5,36,0,72);
hnpemapP->SetOption("COLZ");hsignalmapP->SetOption("COLZ");hsignalRMSmapP->SetOption("COLZ");
hsignalmapM = new TH2F("hsignalmapM","Mean(Response) - Mean(Pedestal) HFM;#eta;#phi",26,-41.5,-28.5,36,0,72);
hsignalRMSmapM = new TH2F("hsignalRMSmapM","RMS Response HFM;#eta;#phi",26,-41.5,-28.5,36,0,72);
hnpemapM = new TH2F("hnpemapM","~N PE input HFM;#eta;#phi",26,-41.5,-28.5,36,0,72);
hnpemapM->SetOption("COLZ");hsignalmapM->SetOption("COLZ");hsignalRMSmapM->SetOption("COLZ");
// Channel-by-channel histos
for (int i=0;i<13;i++) for (int j=0;j<36;j++) for (int k=0;k<2;k++) {
if (i>10 && j%2==0) continue;
sprintf(htit,"ts_+%d_%d_%d",i+29,j*2+1,k+1);
hts[i][j][k] = new TH1F(htit,htit,10,-0.5,9.5); // TimeSlices (pulse shape)
sprintf(htit,"tsmean_+%d_%d_%d",i+29,j*2+1,k+1);
htsm[i][j][k] = new TH1F(htit,htit,100,0,10); // Mean signal time estimated from TS
sprintf(htit,"sp_+%d_%d_%d",i+29,j*2+1,k+1);
hsp[i][j][k] = new TH1F(htit,htit,1201,-25,30000); // Big-scale spectrum (linear ADC)
sprintf(htit,"spe_+%d_%d_%d",i+29,j*2+1,k+1);
hspe[i][j][k] = new TH1F(htit,htit,200,-9.5,190.5); // Small-scale spectrum (linear ADC)
sprintf(htit,"ped_+%d_%d_%d",i+29,j*2+1,k+1);
hped[i][j][k] = new TH1F(htit,htit,200,-9.5,190.5); // Pedestal spectrum
sprintf(htit,"ts_-%d_%d_%d",i+29,j*2+1,k+1);
hts[i+13][j][k] = new TH1F(htit,htit,10,-0.5,9.5);
sprintf(htit,"tsmean_-%d_%d_%d",i+29,j*2+1,k+1);
htsm[i+13][j][k] = new TH1F(htit,htit,100,0,10);
sprintf(htit,"sp_-%d_%d_%d",i+29,j*2+1,k+1);
hsp[i+13][j][k] = new TH1F(htit,htit,1201,-25,30000);
sprintf(htit,"spe_-%d_%d_%d",i+29,j*2+1,k+1);
hspe[i+13][j][k] = new TH1F(htit,htit,200,-9.5,190.5);
sprintf(htit,"ped_-%d_%d_%d",i+29,j*2+1,k+1);
hped[i+13][j][k] = new TH1F(htit,htit,200,-9.5,190.5);
}
// PIN-diodes histos
for (int i=0;i<4;i++) for (int j=0;j<3;j++) {
sprintf(htit,"ts_PIN%d_+Q%d",j+1,i+1);
htspin[i][j] = new TH1F(htit,htit,10,-0.5,9.5);
sprintf(htit,"sp_PIN%d_+Q%d",j+1,i+1);
hsppin[i][j] = new TH1F(htit,htit,1601,-25,40000);
sprintf(htit,"spe_PIN%d_+Q%d",j+1,i+1);
hspepin[i][j] = new TH1F(htit,htit,200,-9.5,190.5);
sprintf(htit,"ped_PIN%d_+Q%d",j+1,i+1);
hpedpin[i][j] = new TH1F(htit,htit,200,-9.5,190.5);
sprintf(htit,"tsmean_PIN%d_+Q%d",j+1,i+1);
htsmpin[i][j] = new TH1F(htit,htit,100,0,10);
sprintf(htit,"ts_PIN%d_-Q%d",j+1,i+1);
htspin[i+4][j] = new TH1F(htit,htit,10,-0.5,9.5);
sprintf(htit,"sp_PIN%d_-Q%d",j+1,i+1);
hsppin[i+4][j] = new TH1F(htit,htit,1601,-25,40000);
sprintf(htit,"spe_PIN%d_-Q%d",j+1,i+1);
hspepin[i+4][j] = new TH1F(htit,htit,200,-9.5,190.5);
sprintf(htit,"ped_PIN%d_-Q%d",j+1,i+1);
hpedpin[i+4][j] = new TH1F(htit,htit,200,-9.5,190.5);
sprintf(htit,"tsmean_PIN%d_-Q%d",j+1,i+1);
htsmpin[i+4][j] = new TH1F(htit,htit,100,0,10);
}
std::cout<<std::endl<<"beginJob: histfile="<<histfile.c_str()<<" textfile="<<textfile.c_str()<<std::endl;
return;
}
| void HFLightCalRand::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 193 of file HFLightCalRand.cc.
References gather_cfg::cout, alignCSCRings::e, EventN, Fit3Peak(), HistSpec(), hnpemapM, hnpemapP, hnpevar, hped, hpedmean, hpedpin, hpedrms, hsignalmapM, hsignalmapP, hsignalmean, hsignalrms, hsignalRMSmapM, hsignalRMSmapP, hsp, hspe, hspepin, hspes, hsppin, hts, htspin, i, j, gen::k, timingPdfMaker::mean, mFile, NEvents, plotscripts::rms(), runNumb, and tFile.
{
Double_t mean,rms,meanped,rmsped,npevar;
Double_t par[5],parm[5],dspe=0,dnpe;
Int_t tsmax,intspe;
fprintf(tFile,"#RunN %d Events processed %d",runNumb,EventN);
std::cout<<"endJob: histos processing..."<<std::endl;
std::cout<<"RunN= "<<runNumb<<" Events processed= "<<EventN<<std::endl;
for (int i=0;i<26;i++) for (int j=0;j<36;j++) for (int k=0;k<2;k++) {
if (i>10 && i<13 && j%2==0) continue;
if (i>23 && j%2==0) continue;
meanped=rmsped=mean=rms=0;
if (hsp[i][j][k]->Integral()>0) {
HistSpec(hped[i][j][k],meanped,rmsped);
HistSpec(hsp[i][j][k],mean,rms);
if (hspe[i][j][k]->Integral()>hsp[i][j][k]->Integral()*0.9 || mean<100) {
HistSpec(hspe[i][j][k],mean,rms);
}
hsignalmean->Fill(mean);
hsignalrms->Fill(rms);
hpedmean->Fill(meanped);
hpedrms->Fill(rmsped);
}
}
meanped=hpedmean->GetMean();
rmsped=hpedrms->GetMean();
mean=hsignalmean->GetMean();
rms=hsignalrms->GetMean();
fprintf(tFile," MeanInput=<%.2f> [linADCcount] RMS=<%.2f>\n",mean,rms);
fprintf(tFile,"#eta/phi/depth sum4maxTS RMS ~N_PE sum4lowTS RMS maxTS SPE +/- Err Comment\n");
TF1* fPed = new TF1("fPed","gaus",0,120);
fPed->SetNpx(200);
TF1 *fTot = new TF1("fTot",Fit3Peak ,0,200,5);
fTot->SetNpx(800);
for (int i=0;i<26;i++) for (int j=0;j<36;j++) for (int k=0;k<2;k++) {
if (i>10 && i<13 && j%2==0) continue;
if (i>23 && j%2==0) continue;
HistSpec(hped[i][j][k],meanped,rmsped);
HistSpec(hsp[i][j][k],mean,rms);
par[3]=0;
if (hspe[i][j][k]->Integral()>hsp[i][j][k]->Integral()*0.9 || mean<100) {
HistSpec(hspe[i][j][k],mean,rms);
if (hspe[i][j][k]->Integral(1,(int) (meanped+3*rmsped+12))/NEvents>0.1) {
//if (hspe[i][j][k]->Integral()>100 && mean-meanped<100) {
if (mean+rms*3-meanped-rmsped*3>1 && rmsped>0) { // SPE fit if low intensity>0
par[1] = meanped;
par[2] = rmsped;
par[0] = hped[i][j][k]->GetMaximum();
fPed->SetParameters(par);
hped[i][j][k]->Fit(fPed,"BQ0");
fPed->GetParameters(&par[0]);
hped[i][j][k]->Fit(fPed,"B0Q","",par[1]-par[2]*3,par[1]+par[2]*3);
fPed->GetParameters(par);
hped[i][j][k]->Fit(fPed,"BLIQ","",par[1]-par[2]*3,par[1]+par[2]*3);
fPed->GetParameters(&par[0]);
fPed->GetParameters(&parm[0]);
NEvents = (int) hspe[i][j][k]->Integral();
par[0]=0.1;
par[3]=10;
par[4]=6;
par[2]=par[2]/0.93;
fTot->SetParameters(par);
fTot->SetParLimits(0,0,2);
//fTot->FixParameter(1,par[1]);
//fTot->FixParameter(2,par[2]);
fTot->SetParLimits(1,par[1]-4,par[1]+4);
fTot->SetParLimits(2,par[2]*0.9,par[2]);
fTot->SetParLimits(3,1.2,50);
//fTot->SetParLimits(4,par[2]*1.1,100);
par[4]=0;
fTot->SetParLimits(4,-1.64,1.64);
hspe[i][j][k]->Fit(fTot,"BLEQ","");
fTot->GetParameters(par);
dspe=fTot->GetParError(3);
dnpe=fTot->GetParError(0);
if (par[3]<1.21 || dnpe>par[0]) par[3]=-1;
else if (par[0]>1.96 || par[3]>49) par[3]=0;
else {
hspes->Fill(par[3]);
}
}
}
}
// NPE
npevar=0;
if (par[3]>0) npevar=par[0]; // NPE from SPE fit
else { // NPE from high intensity signal
if (hspe[i][j][k]->Integral()>hsp[i][j][k]->Integral()*0.98) {
HistSpec(hspe[i][j][k],mean,rms,3);
}
else {
HistSpec(hsp[i][j][k],mean,rms,3);
}
if (rmsped>0) {
if ((rms*rms-rmsped*rmsped/0.8836)>1 && mean>meanped+2) {
npevar=(mean-meanped-2)*(mean-meanped-2)/(rms*rms-rmsped*rmsped/0.8836);
}
else if (mean<100) {
intspe=int(hspe[i][j][k]->Integral());
hspe[i][j][k]->SetAxisRange(meanped+2+rmsped*4,300);
npevar=hspe[i][j][k]->Integral()/intspe;
if (npevar>0.01) npevar=-1;
else npevar=0;
hspe[i][j][k]->SetAxisRange(-20,300);
}
}
}
if (npevar>5.0e-5) hnpevar->Fill(npevar);
if (i<13) {
hsignalmapP->Fill(i+28.6+k/2.0,j*2+1,mean-meanped-1.8);
hsignalRMSmapP->Fill(i+28.6+k/2.0,j*2+1,rms);
if (npevar>0) hnpemapP->Fill(i+28.6+k/2.0,j*2+1,npevar);
fprintf(tFile,"%3d%4d%5d %11.2f%8.2f",i+29,j*2+1,k+1,mean,rms);
}
else {
fprintf(tFile,"%3d%4d%5d %11.2f%8.2f",13-i-29,j*2+1,k+1,mean,rms);
hsignalmapM->Fill(13-i-28.6-k/2.0,j*2+1,mean-meanped-1.8);
hsignalRMSmapM->Fill(13-i-28.6-k/2.0,j*2+1,rms);
if (npevar>0) hnpemapM->Fill(13-i-28.6-k/2.0,j*2+1,npevar);
}
if (npevar>0) fprintf(tFile," %9.4f",npevar);
else fprintf(tFile," 0 ");
fprintf(tFile," %8.2f%8.2f",meanped,rmsped);
tsmax=hts[i][j][k]->GetMaximumBin()-1;
fprintf(tFile," %4d",tsmax);
if (par[3]>0 && par[3]<99) fprintf(tFile,"%8.2f%7.2f",par[3]+1,dspe);
else if (npevar>0) fprintf(tFile,"%8.2f 0 ",(mean-meanped-1)/npevar);
else fprintf(tFile," 0 0 ");
// Diagnostics
fprintf(tFile," ");
if (hsp[i][j][k]->GetEntries()<=0) fprintf(tFile,"NoSignal\n");
else if (hsp[i][j][k]->GetEntries()<=10) fprintf(tFile,"Nev<10\n");
else {
if (hsp[i][j][k]->Integral()<=10 || mean>12000) fprintf(tFile,"SignalOffRange\n");
else {
if (hsp[i][j][k]->Integral()<100) fprintf(tFile,"Nev<100/");
if (npevar>0 && par[3]>0 && (npevar*NEvents<10 || npevar<0.001))
fprintf(tFile,"LowSignal/");
else if (npevar==0 && fabs(mean-meanped)<3) fprintf(tFile,"LowSignal/");
if (par[3]<0) fprintf(tFile,"BadFit/");
else if (par[3]==0) fprintf(tFile,"NoSPEFit/");
else if (par[3]>0 && npevar>1) fprintf(tFile,"NPE>1/");
if (npevar<0) fprintf(tFile,"Problem/");
if (mean<2) fprintf(tFile,"LowMean/");
if (rms<0.5) fprintf(tFile,"LowRMS/");
if (meanped<-1) fprintf(tFile,"LowPed/");
else if (meanped>25) fprintf(tFile,"HighPed/");
if (rmsped<0.5 && rmsped>0) fprintf(tFile,"NarrowPed/");
else if (rmsped>10) fprintf(tFile,"WidePed/");
if (hped[i][j][k]->GetBinContent(201)>10) fprintf(tFile,"PedOffRange");
fprintf(tFile,"-\n");
}
}
}
for (int i=0;i<8;i++) for (int j=0;j<3;j++) {
HistSpec(hpedpin[i][j],meanped,rmsped);
HistSpec(hsppin[i][j],mean,rms);
if (hspepin[i][j]->Integral()>hsppin[i][j]->Integral()*0.9 || mean<100) {
HistSpec(hspepin[i][j],mean,rms);
}
if (i<4) fprintf(tFile," PIN%d +Q%d %12.2f %6.2f",j+1,i+1,mean,rms);
else fprintf(tFile," PIN%d -Q%d %12.2f %6.2f",j+1,i-3,mean,rms);
fprintf(tFile," %15.2f %6.2f",meanped,rmsped);
tsmax=htspin[i][j]->GetMaximumBin()-1;
fprintf(tFile," %4d\n",tsmax);
}
mFile->Write();
mFile->Close();
fclose(tFile);
std::cout<<std::endl<<" --endJob-- done"<<std::endl;
return;
}
Definition at line 47 of file HFLightCalRand.h.
Referenced by analyze().
Definition at line 46 of file HFLightCalRand.h.
Referenced by analyze().
std::string HFLightCalRand::histfile [private] |
Definition at line 27 of file HFLightCalRand.h.
Referenced by beginJob(), and HFLightCalRand().
TH2F * HFLightCalRand::hnpemapM [private] |
Definition at line 38 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH2F* HFLightCalRand::hnpemapP [private] |
Definition at line 38 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH1F * HFLightCalRand::hnpevar [private] |
Definition at line 39 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH1F* HFLightCalRand::hped[26][36][2] [private] |
Definition at line 37 of file HFLightCalRand.h.
Referenced by analyze(), beginJob(), and endJob().
TH1F * HFLightCalRand::hpedmean [private] |
Definition at line 39 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH1F* HFLightCalRand::hpedpin[8][3] [private] |
Definition at line 43 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH1F * HFLightCalRand::hpedrms [private] |
Definition at line 39 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH2F * HFLightCalRand::hsignalmapM [private] |
Definition at line 38 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH2F * HFLightCalRand::hsignalmapP [private] |
Definition at line 38 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH1F* HFLightCalRand::hsignalmean [private] |
Definition at line 39 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH1F * HFLightCalRand::hsignalrms [private] |
Definition at line 39 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH2F * HFLightCalRand::hsignalRMSmapM [private] |
Definition at line 38 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH2F * HFLightCalRand::hsignalRMSmapP [private] |
Definition at line 38 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH1F* HFLightCalRand::hsp[26][36][2] [private] |
Definition at line 35 of file HFLightCalRand.h.
Referenced by analyze(), beginJob(), and endJob().
TH1F* HFLightCalRand::hspe[26][36][2] [private] |
Definition at line 36 of file HFLightCalRand.h.
Referenced by analyze(), beginJob(), and endJob().
TH1F* HFLightCalRand::hspepin[8][3] [private] |
Definition at line 42 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH1F * HFLightCalRand::hspes [private] |
Definition at line 39 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH1F* HFLightCalRand::hsppin[8][3] [private] |
Definition at line 41 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TH1F * HFLightCalRand::htmax [private] |
Definition at line 39 of file HFLightCalRand.h.
Referenced by analyze(), and beginJob().
TH1F * HFLightCalRand::htmean [private] |
Definition at line 39 of file HFLightCalRand.h.
Referenced by analyze(), and beginJob().
TH1F* HFLightCalRand::hts[26][36][2] [private] |
Definition at line 33 of file HFLightCalRand.h.
Referenced by analyze(), beginJob(), and endJob().
TH1F* HFLightCalRand::htsm[26][36][2] [private] |
Definition at line 34 of file HFLightCalRand.h.
Referenced by analyze(), and beginJob().
TH1F* HFLightCalRand::htsmpin[8][3] [private] |
Definition at line 44 of file HFLightCalRand.h.
Referenced by beginJob().
TH1F* HFLightCalRand::htspin[8][3] [private] |
Definition at line 40 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
TFile* HFLightCalRand::mFile [private] |
Definition at line 30 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
FILE* HFLightCalRand::preFile [private] |
Definition at line 32 of file HFLightCalRand.h.
std::string HFLightCalRand::prefile [private] |
Definition at line 29 of file HFLightCalRand.h.
std::string HFLightCalRand::textfile [private] |
Definition at line 28 of file HFLightCalRand.h.
Referenced by beginJob(), and HFLightCalRand().
FILE* HFLightCalRand::tFile [private] |
Definition at line 31 of file HFLightCalRand.h.
Referenced by beginJob(), and endJob().
1.7.3