#include <memory>
#include <string>
#include <iostream>
#include "TH1F.h"
#include "TH2F.h"
#include "TFile.h"
#include "math.h"
#include "TMath.h"
#include "TF1.h"
#include "CalibCalorimetry/HcalStandardModules/interface/HFLightCal.h"
#include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
#include "CalibFormats/HcalObjects/interface/HcalDbService.h"
#include "CalibFormats/HcalObjects/interface/HcalDbRecord.h"
#include "CondFormats/HcalObjects/interface/HcalQIEShape.h"
#include "CondFormats/HcalObjects/interface/HcalQIECoder.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"

Functions
Double_t	FitFun (Double_t x, Double_t par)

void	HistSpecs (TH1F *hist, Double_t &mean, Double_t &rms, Double_t range=4)

Variables
Int_t	eventN =0

Int_t	itsmax [26][36][2]

Int_t	itspinmax [8][3]

Int_t	Nev

Int_t	runN =0

Function Documentation

Double_t FitFun	(	Double_t *	x,
		Double_t *	par
	)

Definition at line 185 of file HFLightCal.cc.

References reco::ParticleMasses::k0, Nev, funct::pow(), diffTwoXMLs::r1, diffTwoXMLs::r2, and mathSSE::sqrt().

Referenced by HFLightCal::endJob().

                                             { 
 // Spectra fit function: Pedestal Gaussian + asymmetric 1PE + 2PE +3PE peaks
 
   Double_t sum,xx,A0,C0,r0,sigma0,mean1,sigma1,A1,C1,r1,mean2,sigma2,A2,C2,r2,mean3,sigma3,A3,C3,r3;
 
   const Double_t k0=2.0, k1=1.6, k2=2.0;
 
   xx=x[0];
   sigma0 = par[2];
   A0 = 2*Nev/(2+2*par[0]+par[0]*par[0]+pow(par[0],3)/3+pow(par[0],4)/12+
           pow(par[0],5)/60+pow(par[0],6)/360);
   r0 = ((xx-par[1])/sigma0);
   C0 = 1/(sigma0* TMath::Exp(-k0*k0/2)/k0 +
       sigma0*sqrt(2*3.14159)*0.5*(1+TMath::Erf(k0/1.41421)));
   //sum = 1/(sqrt(2*3.14159)*par[2])*A0*TMath::Exp(-0.5*r0*r0);
   if(r0 < k0) sum = C0*A0*TMath::Exp(-0.5*r0*r0);
   else sum = C0*A0*TMath::Exp(0.5*k0*k0-k0*r0);
 
   mean1 = par[1]+par[3];
   sigma1 = par[4];
   //sigma1 = 1.547+0.125*par[3]+0.004042*par[3]*par[3];
   //sigma1 = (sigma1+(9.1347e-3+3.845e-2*par[3])*par[4]*2.0)*par[2];
   A1 = A0*par[0];
   C1 = 1/(sigma1* TMath::Exp(-k1*k1/2)/k1 +
       sigma1*sqrt(2*3.14159)*0.5*(1+TMath::Erf(k1/1.41421)));
   r1 = ((xx-mean1)/sigma1);
   if(r1 < k1) sum += C1*A1*TMath::Exp(-0.5*r1*r1);
   else sum += C1*A1*TMath::Exp(0.5*k1*k1-k1*r1);
 
   mean2 = 2*par[3]+par[1];
   sigma2 = sqrt(2*sigma1*sigma1 - pow(par[2],2));
   //A2 = A0*par[5]*par[0]*par[0]/2;
   A2 = A0*par[0]*par[0]/2;
   C2 = 1/(sigma2* TMath::Exp(-k2*k2/2)/k2 +
       sigma2*sqrt(2*3.14159)*0.5*(1+TMath::Erf(k2/1.41421)));
   r2 = ((xx-mean2)/sigma2);
   if(r2 < k2) sum += C2*A2*TMath::Exp(-0.5*r2*r2);
   else sum += C2*A2*TMath::Exp(0.5*k2*k2-k2*r2);
 
   mean3 = 3*par[3]+par[1];
   sigma3 = sqrt(3*sigma1*sigma1 - 2*pow(par[2],2));
   A3 = A0*par[0]*par[0]*par[0]/6;
   C3 = 1/(sigma3*sqrt(2*3.14159));
   r3 = ((xx-mean3)/sigma3);
   sum += C3*A3*TMath::Exp(-0.5*r3*r3);
 
   return sum;
 }

void HistSpecs	(	TH1F *	hist,
		Double_t &	mean,
		Double_t &	rms,
		Double_t	range = `4`
	)

Definition at line 169 of file HFLightCal.cc.

Referenced by HFLightCal::endJob().

                                                                             {
   Double_t xmin,xmax;
   mean=hist->GetMean();
   rms=hist->GetRMS();
   xmin=hist->GetXaxis()->GetXmin();
   xmax=hist->GetXaxis()->GetXmax();
   hist->SetAxisRange(mean-range*rms-100,mean+range*rms+100);
   mean=hist->GetMean();
   rms=hist->GetRMS();
   hist->SetAxisRange(mean-range*rms-100,mean+range*rms+100);
   mean=hist->GetMean();
   rms=hist->GetRMS();
   hist->SetAxisRange(xmin,xmax);
   return;
 }