#include <HFPreLightCal.h>

Inheritance diagram for HFPreLightCal:

Public Member Functions
virtual void	analyze (const edm::Event &fEvent, const edm::EventSetup &fSetup)
virtual void	beginJob ()
virtual void	endJob (void)
	HFPreLightCal (const edm::ParameterSet &fConfiguration)
virtual	~HFPreLightCal ()
Private Attributes
edm::InputTag	hcalCalibDigiCollectionTag_
edm::InputTag	hfDigiCollectionTag_
std::string	histfile
TH1F *	hts [26][36][2]
TH1F *	htsmax
TH1F *	htspin [8][3]
TH1F *	htspinmax
TFile *	mFile
std::string	textfile
FILE *	tFile

Detailed Description

Definition at line 12 of file HFPreLightCal.h.

Constructor & Destructor Documentation

HFPreLightCal::HFPreLightCal ( const edm::ParameterSet & fConfiguration )

Definition at line 38 of file HFPreLightCal.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>("rootPreFile");
  textfile = fConfiguration.getUntrackedParameter<string>("textPreFile");
}

HFPreLightCal::~HFPreLightCal ( ) [virtual]

Definition at line 47 of file HFPreLightCal.cc.

                               {
  //delete mFile;
}

Member Function Documentation

void HFPreLightCal::analyze	(	const edm::Event &	fEvent,
		const edm::EventSetup &	fSetup
	)		`[virtual]`

Implements edm::EDAnalyzer.

Definition at line 127 of file HFPreLightCal.cc.

References ecalMGPA::adc(), calib, gather_cfg::cout, HcalDetId::depth(), edm::EventID::event(), event_NN, HcalObjRepresent::Fill(), edm::Event::getByLabel(), hcalCalibDigiCollectionTag_, HcalForward, hfDigiCollectionTag_, hts, edm::EventBase::id(), HFDataFrame::id(), HcalDetId::ieta(), HcalDetId::iphi(), edm::EventID::run(), run_NN, convertSQLiteXML::runNumber, HFDataFrame::size(), and HcalDetId::subdet().

                                                                               {

  // event ID
  edm::EventID eventId = fEvent.id();
  int runNumber = eventId.run ();
  int eventNumber = eventId.event ();
  if (run_NN==0) run_NN=runNumber;
  event_NN++;
  if (verbose) std::cout << "========================================="<<std::endl
                         << "run/event: "<<runNumber<<'/'<<eventNumber<<std::endl;

  // 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 

  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) {
      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);
      }
    }
  }
  */  
  // 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 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);
    }
  }
  return;
}

void HFPreLightCal::beginJob ( void ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 51 of file HFPreLightCal.cc.

References gather_cfg::cout, cmsRelvalreport::exit, histfile, hts, htsmax, htspin, htspinmax, i, j, gen::k, mFile, NULL, textfile, and tFile.

                             {

  char htit[64];
  std::cout<<std::endl<<"HFPreLightCal beginJob: --> ";

  mFile = new TFile (histfile.c_str(),"RECREATE");
  if ((tFile = fopen(textfile.c_str(),"w"))==NULL) {
    printf("\nNo textfile open\n\n");
    std::cout<<"Problem with output Pre-textFILE => exit"<<std::endl;
    exit(1);
  }
  // Histos
  htsmax = new TH1F("htsmax","Max TS",100,0,10);
  htspinmax = new TH1F("htspinmax","Max TS PIN",100,0,10);
  // 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,"tspre_+%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,"tspre_-%d_%d_%d",i+29,j*2+1,k+1);
    hts[i+13][j][k] = new TH1F(htit,htit,10,-0.5,9.5); // TimeSlices (pulse shape)
  } 
  // PIN-diodes histos
  for (int i=0;i<4;i++) for (int j=0;j<3;j++) {
    sprintf(htit,"tspre_PIN%d_+Q%d",j+1,i+1);
    htspin[i][j] = new TH1F(htit,htit,10,-0.5,9.5);
    sprintf(htit,"tspre_PIN%d_-Q%d",j+1,i+1);
    htspin[i+4][j] = new TH1F(htit,htit,10,-0.5,9.5);
  }
  std::cout<<"histfile="<<histfile.c_str()<<"  textfile="<<textfile.c_str()<<std::endl;
  return;
}

void HFPreLightCal::endJob ( void ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 84 of file HFPreLightCal.cc.

References cont, gather_cfg::cout, event_NN, hts, htsmax, htspin, htspinmax, i, cuy::ii, j, gen::k, mFile, and tFile.

{
  Double_t sum,cont;
  Int_t tsmax;

  std::cout<<std::endl<<"HFPreLightCal endJob --> ";

  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;
    sum=tsmax=0;
    for (int ii=1; ii<=10; ii++) {
      cont = hts[i][j][k]->GetBinContent(ii);
      if (ii<3) cont=cont-(hts[i][j][k]->GetBinContent(ii+4)+hts[i][j][k]->GetBinContent(ii+8))/2;
      else if (ii<5) cont=cont-hts[i][j][k]->GetBinContent(ii+4);
      else if (ii<7) cont=cont-(hts[i][j][k]->GetBinContent(ii-4)+hts[i][j][k]->GetBinContent(ii+4))/2;
      else if (ii<9) cont=cont-hts[i][j][k]->GetBinContent(ii-4);
      else cont=cont-(hts[i][j][k]->GetBinContent(ii-4)+hts[i][j][k]->GetBinContent(ii-8))/2;
      if (cont>sum) {
        sum = cont;
        tsmax=ii;
      }
    }
    htsmax->Fill(tsmax);
    if (i<13) fprintf(tFile," %d  %d  %d  %d\n",i+29,j*2+1,k+1,tsmax);
    else      fprintf(tFile," %d  %d  %d  %d\n",13-i-29,j*2+1,k+1,tsmax);
  }

  for (int i=0;i<8;i++) for (int j=0;j<3;j++) {
    sum=tsmax=0;
    tsmax = htspin[i][j]->GetMaximumBin();
    htspinmax->Fill(tsmax);
    if (i<4) fprintf(tFile,"%d  %d  %d\n",i+1,j+1,tsmax);
    else     fprintf(tFile,"%d  %d  %d\n",-i+3,j+1,tsmax);
  } 

  mFile->Write();
  mFile->Close();
  fclose(tFile);
  std::cout<<" Nevents = "<<event_NN<<std::endl;
  return;
}