EgammaCoreTools. More...

#include <PiZeroAnalyzer.h>

Inheritance diagram for PiZeroAnalyzer:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	PiZeroAnalyzer (const edm::ParameterSet &)

virtual	~PiZeroAnalyzer ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
void	makePizero (const edm::EventSetup &es, const edm::Handle< EcalRecHitCollection > eb, const edm::Handle< EcalRecHitCollection > ee)

Private Attributes
edm::EDGetTokenT < edm::SortedCollection < EcalRecHit, edm::StrictWeakOrdering < EcalRecHit > > >	barrelEcalHits_token_

int	clusEtaSize_

int	clusPhiSize_

double	clusSeedThr_
	parameters needed for pizero finding More...

std::stringstream	currentFolder_

double	cutStep_

edm::EDGetTokenT < edm::SortedCollection < EcalRecHit, edm::StrictWeakOrdering < EcalRecHit > > >	endcapEcalHits_token_

std::string	fName_

MonitorElement *	hIsoPi0EB_

MonitorElement *	hMinvPi0EB_

MonitorElement *	hPt1Pi0EB_

MonitorElement *	hPt2Pi0EB_

MonitorElement *	hPtPi0EB_

double	minPhoEtCut_

int	nEvt_

int	numberOfSteps_

bool	ParameterLogWeighted_

double	ParameterT0_barl_

double	ParameterW0_

double	ParameterX0_

edm::ParameterSet	posCalcParameters_

unsigned int	prescaleFactor_

double	seleMinvMaxPi0_

double	seleMinvMinPi0_

double	selePi0BeltDeta_

double	selePi0BeltDR_

double	selePi0Iso_

double	selePtGammaOne_

double	selePtGammaTwo_

double	selePtPi0_

double	seleS4S9GammaOne_

double	seleS4S9GammaTwo_

double	seleXtalMinEnergy_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

EgammaCoreTools.

$Id: PiZeroAnalyzer authors: Nancy Marinelli, U. of Notre Dame, US Jamie Antonelli, U. of Notre Dame, US

Definition at line 77 of file PiZeroAnalyzer.h.

Constructor & Destructor Documentation

PiZeroAnalyzer::PiZeroAnalyzer ( const edm::ParameterSet & pset )

explicit

Definition at line 20 of file PiZeroAnalyzer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   fName_              = pset.getUntrackedParameter<std::string>("Name");
   prescaleFactor_     = pset.getUntrackedParameter<int>("prescaleFactor",1);
 
   barrelEcalHits_token_ = consumes<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > >(pset.getParameter<edm::InputTag>("barrelEcalHits"));
   endcapEcalHits_token_ = consumes<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > >(pset.getParameter<edm::InputTag>("endcapEcalHits"));
 
   nEvt_=0;
 
   // parameters for Pizero finding
   seleXtalMinEnergy_    = pset.getParameter<double> ("seleXtalMinEnergy");
   clusSeedThr_          = pset.getParameter<double> ("clusSeedThr");
   clusEtaSize_          = pset.getParameter<int> ("clusEtaSize");
   clusPhiSize_          = pset.getParameter<int> ("clusPhiSize");
   ParameterLogWeighted_ = pset.getParameter<bool> ("ParameterLogWeighted");
   ParameterX0_          = pset.getParameter<double> ("ParameterX0");
   ParameterT0_barl_     = pset.getParameter<double> ("ParameterT0_barl");
   ParameterW0_          = pset.getParameter<double> ("ParameterW0");
 
   selePtGammaOne_       = pset.getParameter<double> ("selePtGammaOne");
   selePtGammaTwo_       = pset.getParameter<double> ("selePtGammaTwo");
   seleS4S9GammaOne_     = pset.getParameter<double> ("seleS4S9GammaOne");
   seleS4S9GammaTwo_     = pset.getParameter<double> ("seleS4S9GammaTwo");
   selePtPi0_            = pset.getParameter<double> ("selePtPi0");
   selePi0Iso_           = pset.getParameter<double> ("selePi0Iso");
   selePi0BeltDR_        = pset.getParameter<double> ("selePi0BeltDR");
   selePi0BeltDeta_      = pset.getParameter<double> ("selePi0BeltDeta");
   seleMinvMaxPi0_       = pset.getParameter<double> ("seleMinvMaxPi0");
   seleMinvMinPi0_       = pset.getParameter<double> ("seleMinvMinPi0");
 
   posCalcParameters_    = pset.getParameter<edm::ParameterSet>("posCalcParameters");
 }

PiZeroAnalyzer::~PiZeroAnalyzer ( )

virtual

Definition at line 54 of file PiZeroAnalyzer.cc.

54 {

55 }

Member Function Documentation

void PiZeroAnalyzer::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	esup
	)

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 81 of file PiZeroAnalyzer.cc.

References edm::Event::getByToken(), and edm::EventBase::id().

 {
   using namespace edm;
 
   if (nEvt_% prescaleFactor_ ) return;
   nEvt_++;
   LogInfo("PiZeroAnalyzer") << "PiZeroAnalyzer Analyzing event number: " << e.id() << " Global Counter " << nEvt_ <<"\n";
 
   // Get EcalRecHits
   bool validEcalRecHits=true;
   Handle<EcalRecHitCollection> barrelHitHandle;
   EcalRecHitCollection barrelRecHits;
   e.getByToken(barrelEcalHits_token_, barrelHitHandle);
   if (!barrelHitHandle.isValid()) {
     edm::LogError("PhotonProducer") << "Error! Can't get the product: barrelEcalHits_token_";
     validEcalRecHits=false;
   }
 
   Handle<EcalRecHitCollection> endcapHitHandle;
   e.getByToken(endcapEcalHits_token_, endcapHitHandle);
   EcalRecHitCollection endcapRecHits;
   if (!endcapHitHandle.isValid()) {
     edm::LogError("PhotonProducer") << "Error! Can't get the product: endcapEcalHits_token_";
     validEcalRecHits=false;
   }
 
   if (validEcalRecHits) makePizero(esup,  barrelHitHandle, endcapHitHandle);
 }

void PiZeroAnalyzer::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 57 of file PiZeroAnalyzer.cc.

References DQMStore::IBooker::book1D(), MonitorElement::setAxisTitle(), and DQMStore::IBooker::setCurrentFolder().

 {
   currentFolder_.str("");
   currentFolder_ << "Egamma/PiZeroAnalyzer/";
   ibooker.setCurrentFolder(currentFolder_.str());
 
   hMinvPi0EB_ = ibooker.book1D("Pi0InvmassEB","Pi0 Invariant Mass in EB",100,0.,0.5);
   hMinvPi0EB_->setAxisTitle("Inv Mass [GeV] ",1);
 
   hPt1Pi0EB_ = ibooker.book1D("Pt1Pi0EB","Pt 1st most energetic Pi0 photon in EB",100,0.,20.);
   hPt1Pi0EB_->setAxisTitle("1st photon Pt [GeV] ",1);
 
   hPt2Pi0EB_ = ibooker.book1D("Pt2Pi0EB","Pt 2nd most energetic Pi0 photon in EB",100,0.,20.);
   hPt2Pi0EB_->setAxisTitle("2nd photon Pt [GeV] ",1);
 
   hPtPi0EB_ = ibooker.book1D("PtPi0EB","Pi0 Pt in EB",100,0.,20.);
   hPtPi0EB_->setAxisTitle("Pi0 Pt [GeV] ",1);
 
   hIsoPi0EB_ = ibooker.book1D("IsoPi0EB","Pi0 Iso in EB",50,0.,1.);
   hIsoPi0EB_->setAxisTitle("Pi0 Iso",1);
 }

void PiZeroAnalyzer::makePizero	(	const edm::EventSetup &	es,
		const edm::Handle< EcalRecHitCollection >	eb,
		const edm::Handle< EcalRecHitCollection >	ee
	)

private

Definition at line 110 of file PiZeroAnalyzer.cc.

References PositionCalc::Calculate_Location(), funct::cos(), gather_cfg::cout, DetId::Ecal, EcalBarrel, EcalPreshower, relval_parameters_module::energy, create_public_lumi_plots::exp, edm::EventSetup::get(), CaloSubdetectorTopology::getWindow(), i, EBDetId::ieta(), EBDetId::iphi(), j, relval_2017::k, edm::Handle< T >::product(), funct::sin(), and mathSSE::sqrt().

                                                                                                                                                  {
   const EcalRecHitCollection *hitCollection_p = rhEB.product();
 
   edm::ESHandle<CaloGeometry> geoHandle;
   es.get<CaloGeometryRecord>().get(geoHandle);
 
   edm::ESHandle<CaloTopology> theCaloTopology;
   es.get<CaloTopologyRecord>().get(theCaloTopology);
 
   const CaloSubdetectorGeometry *geometry_p;
   const CaloSubdetectorTopology *topology_p;
   const CaloSubdetectorGeometry *geometryES_p;
   geometry_p = geoHandle->getSubdetectorGeometry(DetId::Ecal,EcalBarrel);
   geometryES_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
 
   // Parameters for the position calculation:
   PositionCalc posCalculator_ = PositionCalc(posCalcParameters_);
   //
   std::map<DetId, EcalRecHit> recHitsEB_map;
   //
   std::vector<EcalRecHit> seeds;
 
   seeds.clear();
   //
   vector<EBDetId> usedXtals;
   usedXtals.clear();
   //
   EcalRecHitCollection::const_iterator itb;
   //
   static const int MAXCLUS = 2000;
   int nClus=0;
   vector<float> eClus;
   vector<float> etClus;
   vector<float> etaClus;
   vector<float> phiClus;
   vector<EBDetId> max_hit;
   vector< vector<EcalRecHit> > RecHitsCluster;
   vector<float> s4s9Clus;
 
   // find cluster seeds in EB
   for(itb=rhEB->begin(); itb!=rhEB->end(); ++itb){
     EBDetId id(itb->id());
     double energy = itb->energy();
     if (energy > seleXtalMinEnergy_) {
       std::pair<DetId, EcalRecHit> map_entry(itb->id(), *itb);
       recHitsEB_map.insert(map_entry);
     }
     if (energy > clusSeedThr_) seeds.push_back(*itb);
   } // Eb rechits
 
   sort(seeds.begin(), seeds.end(), ecalRecHitLess());
   for (std::vector<EcalRecHit>::iterator itseed=seeds.begin(); itseed!=seeds.end(); itseed++) {
     EBDetId seed_id = itseed->id();
     std::vector<EBDetId>::const_iterator usedIds;
 
     bool seedAlreadyUsed=false;
     for(usedIds=usedXtals.begin(); usedIds!=usedXtals.end(); usedIds++){
       if(*usedIds==seed_id){
     seedAlreadyUsed=true;
     //cout<< " Seed with energy "<<itseed->energy()<<" was used !"<<endl;
     break;
       }
     }
     if(seedAlreadyUsed)continue;
     topology_p = theCaloTopology->getSubdetectorTopology(DetId::Ecal,EcalBarrel);
     std::vector<DetId> clus_v = topology_p->getWindow(seed_id,clusEtaSize_,clusPhiSize_);
     //std::vector<DetId> clus_used;
     std::vector<std::pair<DetId, float> > clus_used;
 
     vector<EcalRecHit> RecHitsInWindow;
 
     double simple_energy = 0;
 
     for (std::vector<DetId>::iterator det=clus_v.begin(); det!=clus_v.end(); det++) {
       // EBDetId EBdet = *det;
       //      cout<<" det "<< EBdet<<" ieta "<<EBdet.ieta()<<" iphi "<<EBdet.iphi()<<endl;
       bool  HitAlreadyUsed=false;
       for(usedIds=usedXtals.begin(); usedIds!=usedXtals.end(); usedIds++){
     if(*usedIds==*det){
       HitAlreadyUsed=true;
       break;
     }
       }
       if(HitAlreadyUsed)continue;
       if (recHitsEB_map.find(*det) != recHitsEB_map.end()){
     //      cout<<" Used det "<< EBdet<<endl;
     std::map<DetId, EcalRecHit>::iterator aHit;
     aHit = recHitsEB_map.find(*det);
     usedXtals.push_back(*det);
     RecHitsInWindow.push_back(aHit->second);
     clus_used.push_back( std::pair<DetId, float>(*det, 1.) );
     simple_energy = simple_energy + aHit->second.energy();
       }
     }
 
     math::XYZPoint clus_pos = posCalculator_.Calculate_Location(clus_used,hitCollection_p,geometry_p,geometryES_p);
     float theta_s = 2. * atan(exp(-clus_pos.eta()));
     float p0x_s = simple_energy * sin(theta_s) * cos(clus_pos.phi());
     float p0y_s = simple_energy * sin(theta_s) * sin(clus_pos.phi());
     //      float p0z_s = simple_energy * cos(theta_s);
     float et_s = sqrt( p0x_s*p0x_s + p0y_s*p0y_s);
 
     //cout << "       Simple Clustering: E,Et,px,py,pz: "<<simple_energy<<" "<<et_s<<" "<<p0x_s<<" "<<p0y_s<<" "<<endl;
 
     eClus.push_back(simple_energy);
     etClus.push_back(et_s);
     etaClus.push_back(clus_pos.eta());
     phiClus.push_back(clus_pos.phi());
     max_hit.push_back(seed_id);
     RecHitsCluster.push_back(RecHitsInWindow);
     //Compute S4/S9 variable
     //We are not sure to have 9 RecHits so need to check eta and phi:
     float s4s9_[4];
     for(int i=0;i<4;i++)s4s9_[i]= itseed->energy();
     for(unsigned int j=0; j<RecHitsInWindow.size();j++){
       //cout << " Simple cluster rh, ieta, iphi : "<<((EBDetId)RecHitsInWindow[j].id()).ieta()<<" "<<((EBDetId)RecHitsInWindow[j].id()).iphi()<<endl;
       if((((EBDetId)RecHitsInWindow[j].id()).ieta() == seed_id.ieta()-1 && seed_id.ieta()!=1 ) || ( seed_id.ieta()==1 && (((EBDetId)RecHitsInWindow[j].id()).ieta() == seed_id.ieta()-2))){
     if(((EBDetId)RecHitsInWindow[j].id()).iphi() == seed_id.iphi()-1 ||((EBDetId)RecHitsInWindow[j].id()).iphi()-360 == seed_id.iphi()-1 ){
       s4s9_[0]+=RecHitsInWindow[j].energy();
     }else{
       if(((EBDetId)RecHitsInWindow[j].id()).iphi() == seed_id.iphi()){
         s4s9_[0]+=RecHitsInWindow[j].energy();
         s4s9_[1]+=RecHitsInWindow[j].energy();
       }else{
         if(((EBDetId)RecHitsInWindow[j].id()).iphi() == seed_id.iphi()+1 ||((EBDetId)RecHitsInWindow[j].id()).iphi()-360 == seed_id.iphi()+1 ){
           s4s9_[1]+=RecHitsInWindow[j].energy();
         }
       }
     }
       }else{
     if(((EBDetId)RecHitsInWindow[j].id()).ieta() == seed_id.ieta()){
       if(((EBDetId)RecHitsInWindow[j].id()).iphi() == seed_id.iphi()-1 ||((EBDetId)RecHitsInWindow[j].id()).iphi()-360 == seed_id.iphi()-1 ){
         s4s9_[0]+=RecHitsInWindow[j].energy();
         s4s9_[3]+=RecHitsInWindow[j].energy();
       }else{
         if(((EBDetId)RecHitsInWindow[j].id()).iphi() == seed_id.iphi()+1 ||((EBDetId)RecHitsInWindow[j].id()).iphi()-360 == seed_id.iphi()+1 ){
           s4s9_[1]+=RecHitsInWindow[j].energy();
           s4s9_[2]+=RecHitsInWindow[j].energy();
         }
       }
     }else{
       if((((EBDetId)RecHitsInWindow[j].id()).ieta() == seed_id.ieta()+1 && seed_id.ieta()!=-1 ) || ( seed_id.ieta()==-1 && (((EBDetId)RecHitsInWindow[j].id()).ieta() == seed_id.ieta()+2))){
         if(((EBDetId)RecHitsInWindow[j].id()).iphi() == seed_id.iphi()-1 ||((EBDetId)RecHitsInWindow[j].id()).iphi()-360 == seed_id.iphi()-1 ){
           s4s9_[3]+=RecHitsInWindow[j].energy();
         }else{
           if(((EBDetId)RecHitsInWindow[j].id()).iphi() == seed_id.iphi()){
         s4s9_[2]+=RecHitsInWindow[j].energy();
         s4s9_[3]+=RecHitsInWindow[j].energy();
           }else{
         if(((EBDetId)RecHitsInWindow[j].id()).iphi() == seed_id.iphi()+1 ||((EBDetId)RecHitsInWindow[j].id()).iphi()-360 == seed_id.iphi()+1 ){
           s4s9_[2]+=RecHitsInWindow[j].energy();
         }
           }
         }
       }else{
         cout<<" (EBDetId)RecHitsInWindow[j].id()).ieta() "<<((EBDetId)RecHitsInWindow[j].id()).ieta()<<" seed_id.ieta() "<<seed_id.ieta()<<endl;
         cout<<" Problem with S4 calculation "<<endl;return;
       }
     }
       }
     }
     s4s9Clus.push_back(*max_element( s4s9_,s4s9_+4)/simple_energy);
     //    cout<<" s4s9Clus[0] "<<s4s9_[0]/simple_energy<<" s4s9Clus[1] "<<s4s9_[1]/simple_energy<<" s4s9Clus[2] "<<s4s9_[2]/simple_energy<<" s4s9Clus[3] "<<s4s9_[3]/simple_energy<<endl;
     //    cout<<" Max "<<*max_element( s4s9_,s4s9_+4)/simple_energy<<endl;
     nClus++;
     if (nClus == MAXCLUS) return;
   }  //  End loop over seed clusters
 
   // cout<< " Pi0 clusters: "<<nClus<<endl;
 
   // Selection, based on Simple clustering
   //pi0 candidates
   static const int MAXPI0S = 200;
   int npi0_s=0;
 
   vector<EBDetId> scXtals;
   scXtals.clear();
 
   if (nClus <= 1) return;
   for(Int_t i=0 ; i<nClus ; i++){
     for(Int_t j=i+1 ; j<nClus ; j++){
       //      cout<<" i "<<i<<"  etClus[i] "<<etClus[i]<<" j "<<j<<"  etClus[j] "<<etClus[j]<<endl;
       if( etClus[i]>selePtGammaOne_ && etClus[j]>selePtGammaTwo_ && s4s9Clus[i]>seleS4S9GammaOne_ && s4s9Clus[j]>seleS4S9GammaTwo_){
     float theta_0 = 2. * atan(exp(-etaClus[i]));
     float theta_1 = 2. * atan(exp(-etaClus[j]));
 
     float p0x = eClus[i] * sin(theta_0) * cos(phiClus[i]);
     float p1x = eClus[j] * sin(theta_1) * cos(phiClus[j]);
     float p0y = eClus[i] * sin(theta_0) * sin(phiClus[i]);
     float p1y = eClus[j] * sin(theta_1) * sin(phiClus[j]);
     float p0z = eClus[i] * cos(theta_0);
     float p1z = eClus[j] * cos(theta_1);
 
     float pt_pi0 = sqrt( (p0x+p1x)*(p0x+p1x) + (p0y+p1y)*(p0y+p1y));
     //      cout<<" pt_pi0 "<<pt_pi0<<endl;
     if (pt_pi0 < selePtPi0_)continue;
     float m_inv = sqrt ( (eClus[i] + eClus[j])*(eClus[i] + eClus[j]) - (p0x+p1x)*(p0x+p1x) - (p0y+p1y)*(p0y+p1y) - (p0z+p1z)*(p0z+p1z) );
     if ( (m_inv<seleMinvMaxPi0_) && (m_inv>seleMinvMinPi0_) ){
 
       //New Loop on cluster to measure isolation:
       vector<int> IsoClus;
       IsoClus.clear();
       float Iso = 0;
       TVector3 pi0vect = TVector3((p0x+p1x), (p0y+p1y), (p0z+p1z));
       for(Int_t k=0 ; k<nClus ; k++){
         if(k==i || k==j)continue;
         TVector3 Clusvect = TVector3(eClus[k] * sin(2. * atan(exp(-etaClus[k]))) * cos(phiClus[k]), eClus[k] * sin(2. * atan(exp(-etaClus[k]))) * sin(phiClus[k]) , eClus[k] * cos(2. * atan(exp(-etaClus[k]))));
         float dretaclpi0 = fabs(etaClus[k] - pi0vect.Eta());
         float drclpi0 = Clusvect.DeltaR(pi0vect);
 
         if((drclpi0<selePi0BeltDR_) && (dretaclpi0<selePi0BeltDeta_) ){
 
           Iso = Iso + etClus[k];
           IsoClus.push_back(k);
         }
       }
 
 
       if(Iso/pt_pi0<selePi0Iso_){
 
         hMinvPi0EB_->Fill(m_inv);
         hPt1Pi0EB_->Fill(etClus[i]);
         hPt2Pi0EB_->Fill(etClus[j]);
         hPtPi0EB_->Fill(pt_pi0);
         hIsoPi0EB_->Fill(Iso/pt_pi0);
 
 
         npi0_s++;
       }
 
       if(npi0_s == MAXPI0S) return;
     }
       }
     }
   }
 }

Member Data Documentation

edm::EDGetTokenT<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > > PiZeroAnalyzer::barrelEcalHits_token_

private

Definition at line 96 of file PiZeroAnalyzer.h.

int PiZeroAnalyzer::clusEtaSize_

private

Definition at line 105 of file PiZeroAnalyzer.h.

int PiZeroAnalyzer::clusPhiSize_

private

Definition at line 106 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::clusSeedThr_

private

parameters needed for pizero finding

Definition at line 104 of file PiZeroAnalyzer.h.

std::stringstream PiZeroAnalyzer::currentFolder_

private

Definition at line 126 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::cutStep_

private

Definition at line 100 of file PiZeroAnalyzer.h.

edm::EDGetTokenT<edm::SortedCollection<EcalRecHit,edm::StrictWeakOrdering<EcalRecHit> > > PiZeroAnalyzer::endcapEcalHits_token_

private

Definition at line 97 of file PiZeroAnalyzer.h.

std::string PiZeroAnalyzer::fName_

private

Definition at line 89 of file PiZeroAnalyzer.h.

MonitorElement* PiZeroAnalyzer::hIsoPi0EB_

private

Definition at line 131 of file PiZeroAnalyzer.h.

MonitorElement* PiZeroAnalyzer::hMinvPi0EB_

private

Definition at line 128 of file PiZeroAnalyzer.h.

MonitorElement* PiZeroAnalyzer::hPt1Pi0EB_

private

Definition at line 129 of file PiZeroAnalyzer.h.

MonitorElement* PiZeroAnalyzer::hPt2Pi0EB_

private

Definition at line 130 of file PiZeroAnalyzer.h.

MonitorElement* PiZeroAnalyzer::hPtPi0EB_

private

Definition at line 132 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::minPhoEtCut_

private

Definition at line 98 of file PiZeroAnalyzer.h.

int PiZeroAnalyzer::nEvt_

private

Definition at line 92 of file PiZeroAnalyzer.h.

int PiZeroAnalyzer::numberOfSteps_

private

Definition at line 101 of file PiZeroAnalyzer.h.

bool PiZeroAnalyzer::ParameterLogWeighted_

private

Definition at line 110 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::ParameterT0_barl_

private

Definition at line 112 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::ParameterW0_

private

Definition at line 113 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::ParameterX0_

private

Definition at line 111 of file PiZeroAnalyzer.h.

edm::ParameterSet PiZeroAnalyzer::posCalcParameters_

private

Definition at line 94 of file PiZeroAnalyzer.h.

unsigned int PiZeroAnalyzer::prescaleFactor_

private

Definition at line 90 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::seleMinvMaxPi0_

private

Definition at line 123 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::seleMinvMinPi0_

private

Definition at line 124 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::selePi0BeltDeta_

private

Definition at line 121 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::selePi0BeltDR_

private

Definition at line 120 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::selePi0Iso_

private

Definition at line 122 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::selePtGammaOne_

private

Definition at line 115 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::selePtGammaTwo_

private

Definition at line 116 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::selePtPi0_

private

Definition at line 117 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::seleS4S9GammaOne_

private

Definition at line 118 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::seleS4S9GammaTwo_

private

Definition at line 119 of file PiZeroAnalyzer.h.

double PiZeroAnalyzer::seleXtalMinEnergy_

private

Definition at line 108 of file PiZeroAnalyzer.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation