AnalyzerWithCentrality Class Reference

Inheritance diagram for AnalyzerWithCentrality:

Public Member Functions
	AnalyzerWithCentrality (const edm::ParameterSet &)
	~AnalyzerWithCentrality ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	endJob ()

Private Attributes
CentralityProvider *	centrality_
edm::Service< TFileService >	fs
TH1D *	h1
TNtuple *	nt

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 27 of file AnalyzerWithCentrality.cc.

Constructor & Destructor Documentation

AnalyzerWithCentrality::AnalyzerWithCentrality ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 57 of file AnalyzerWithCentrality.cc.

References fs, h1, TFileDirectory::make(), and nt.

                                                                             : 
centrality_(0)
{
   //now do what ever initialization is needed
   h1 = fs->make<TH1D>("h1","histogram",100,0,100);
   nt = fs->make<TNtuple>("hi","hi","hf:hft:hftp:hftm:eb:ee:eep:eem:npix:et:zdc:zdcp:zdcm:bin:trig");
}

AnalyzerWithCentrality::~AnalyzerWithCentrality

(

)

Definition at line 65 of file AnalyzerWithCentrality.cc.

{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

void AnalyzerWithCentrality::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 80 of file AnalyzerWithCentrality.cc.

References newFWLiteAna::bin, CentralityProvider::bMean(), CentralityBins::bMeanOfBin(), CentralityProvider::bSigma(), CentralityBins::bSigmaOfBin(), centrality_, CentralityProvider::centralityValue(), gather_cfg::cout, reco::Centrality::EtEBSum(), reco::Centrality::EtEESum(), reco::Centrality::EtEESumMinus(), reco::Centrality::EtEESumPlus(), reco::Centrality::EtHFhitSum(), reco::Centrality::EtHFtowerSum(), reco::Centrality::EtHFtowerSumMinus(), reco::Centrality::EtHFtowerSumPlus(), reco::Centrality::EtMidRapiditySum(), CentralityProvider::getBin(), CentralityProvider::getNbins(), h1, reco::Centrality::multiplicityPixel(), pileupCalc::nbins, CentralityProvider::NcollMean(), CentralityBins::NcollMeanOfBin(), CentralityProvider::NcollSigma(), CentralityBins::NcollSigmaOfBin(), CentralityProvider::newEvent(), CentralityProvider::NhardMean(), CentralityProvider::NhardSigma(), CentralityProvider::NpartMean(), CentralityBins::NpartMeanOfBin(), CentralityProvider::NpartSigma(), CentralityBins::NpartSigmaOfBin(), nt, CentralityProvider::print(), CentralityProvider::raw(), reco::Centrality::zdcSum(), reco::Centrality::zdcSumMinus(), and reco::Centrality::zdcSumPlus().

{
   using namespace edm;
   if(!centrality_) centrality_ = new CentralityProvider(iSetup);
   centrality_->newEvent(iEvent,iSetup);

   double hf = centrality_->raw()->EtHFhitSum();
   double hft = centrality_->raw()->EtHFtowerSum();
   double hftp = centrality_->raw()->EtHFtowerSumPlus();
   double hftm = centrality_->raw()->EtHFtowerSumMinus();
   double eb = centrality_->raw()->EtEBSum();
   double ee = centrality_->raw()->EtEESum();
   double eep = centrality_->raw()->EtEESumPlus();
   double eem = centrality_->raw()->EtEESumMinus();
   double zdc = centrality_->raw()->zdcSum();
   double zdcm = centrality_->raw()->zdcSumMinus();
   double zdcp = centrality_->raw()->zdcSumPlus();
   double npix = centrality_->raw()->multiplicityPixel();
   double et = centrality_->raw()->EtMidRapiditySum();

   cout<<"Centrality variables in the event:"<<endl;
   cout<<"Total energy in HF hits : "<<hf<<endl;
   cout<<"Asymmetry of HF towers : "<<fabs(hftp-hftm)/(hftp+hftm)<<endl;
   cout<<"Total energy in EE basic clusters : "<<eep+eem<<endl;
   cout<<"Total energy in EB basic clusters : "<<eb<<endl;
   
   centrality_->print();
   
   cout<<"Centrality of the event : "<<centrality_->centralityValue()<<endl;

   int bin = centrality_->getBin();
   cout<<"a"<<endl;
   nt->Fill(hf,hft,hftp,hftm,eb,ee,eep,eem,npix,et,zdc,zdcp,zdcm,bin,1);
   cout<<"b"<<endl;

   h1->Fill(bin);
   cout<<"c"<<endl;

   int nbins = centrality_->getNbins(); 
   cout<<"d"<<endl;
   int binsize = 100./nbins;
   cout<<"e"<<endl;
   char* binName = Form("%d to % d",bin*binsize,(bin+1)*binsize);
   cout<<"The event falls into centrality bin : "<<binName<<" id : "<<bin<<endl;

   double npartMean = centrality_->NpartMean();
   double npartSigma = centrality_->NpartSigma();
   cout<<"Npart Mean : "<<npartMean<<"   Variance : "<<npartSigma<<endl;

   // or, alternatively,
   npartMean = centrality_->NpartMeanOfBin(bin);
   npartSigma = centrality_->NpartSigmaOfBin(bin);

   double ncollMean = centrality_->NcollMean();
   double ncollSigma = centrality_->NcollSigma();
   cout<<"Ncoll Mean : "<<ncollMean<<"   Variance : "<<ncollSigma<<endl;

   // or, alternatively,
   ncollMean = centrality_->NcollMeanOfBin(bin);
   ncollSigma = centrality_->NcollSigmaOfBin(bin);


   double nhardMean = centrality_->NhardMean();
   double nhardSigma = centrality_->NhardSigma();
   cout<<"Nhard Mean : "<<nhardMean<<"   Variance : "<<nhardSigma<<endl;

   double bMean = centrality_->bMean();
   double bSigma = centrality_->bSigma();
   cout<<"b Mean : "<<bMean<<"   Variance : "<<bSigma<<endl;

   bMean = centrality_->bMeanOfBin(bin);
   bSigma = centrality_->bSigmaOfBin(bin);


}

void AnalyzerWithCentrality::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 159 of file AnalyzerWithCentrality.cc.

{
}

void AnalyzerWithCentrality::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 165 of file AnalyzerWithCentrality.cc.

                               {
}

Member Data Documentation

CentralityProvider* AnalyzerWithCentrality::centrality_

private

Definition at line 40 of file AnalyzerWithCentrality.cc.

Referenced by analyze().

edm::Service<TFileService> AnalyzerWithCentrality::fs

private

Definition at line 41 of file AnalyzerWithCentrality.cc.

Referenced by AnalyzerWithCentrality().

TH1D* AnalyzerWithCentrality::h1

private

Definition at line 42 of file AnalyzerWithCentrality.cc.

Referenced by utils.Chi2::absval(), analyze(), AnalyzerWithCentrality(), utils.Chi2::check_filled_bins(), utils.KS::do_test(), utils.Chi2::do_test(), utils.BinToBin::do_test(), utils.BinToBin1percent::do_test(), utils.StatisticalTest::get_rank(), and utils.StatisticalTest::set_operands().

TNtuple* AnalyzerWithCentrality::nt

private

Definition at line 43 of file AnalyzerWithCentrality.cc.

Referenced by analyze(), and AnalyzerWithCentrality().