HFPreLightCal Class Reference

#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 ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Attributes
std::string	histfile
TH1F *	hts [26][36][2]
TH1F *	htsmax
TH1F *	htspin [8][3]
TH1F *	htspinmax
TFile *	mFile
std::string	textfile
FILE *	tFile

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 11 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().

                                                                   {
  //std::string histfile = fConfiguration.getUntrackedParameter<string>("rootFile");
  histfile = fConfiguration.getUntrackedParameter<string>("rootPreFile");
  textfile = fConfiguration.getUntrackedParameter<string>("textPreFile");
}

HFPreLightCal::~HFPreLightCal ( )

virtual

Definition at line 44 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 124 of file HFPreLightCal.cc.

References ecalMGPA::adc(), calib, gather_cfg::cout, HcalDetId::depth(), edm::EventID::event(), event_NN, HcalObjRepresent::Fill(), edm::Event::getByType(), HcalForward, HFDataFrame::id(), edm::EventBase::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.getByType(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.getByType(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 48 of file HFPreLightCal.cc.

References gather_cfg::cout, cmsRelvalreport::exit, i, j, gen::k, and NULL.

                             {

  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 81 of file HFPreLightCal.cc.

References cont, gather_cfg::cout, event_NN, i, j, and gen::k.

{
  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;
}

Member Data Documentation

std::string HFPreLightCal::histfile

private

Definition at line 26 of file HFPreLightCal.h.

TH1F* HFPreLightCal::hts[26][36][2]

private

Definition at line 30 of file HFPreLightCal.h.

TH1F* HFPreLightCal::htsmax

private

Definition at line 31 of file HFPreLightCal.h.

TH1F* HFPreLightCal::htspin[8][3]

private

Definition at line 32 of file HFPreLightCal.h.

TH1F * HFPreLightCal::htspinmax

private

Definition at line 31 of file HFPreLightCal.h.

TFile* HFPreLightCal::mFile

private

Definition at line 28 of file HFPreLightCal.h.

std::string HFPreLightCal::textfile

private

Definition at line 27 of file HFPreLightCal.h.

FILE* HFPreLightCal::tFile

private

Definition at line 29 of file HFPreLightCal.h.