EgammaCoreTools. More...

#include <PiZeroAnalyzer.h>

Inheritance diagram for PiZeroAnalyzer:

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

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

	PiZeroAnalyzer (const edm::ParameterSet &)

	~PiZeroAnalyzer () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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_

const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_

const edm::ESGetToken< CaloTopology, CaloTopologyRecord >	caloTopologyToken_

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 DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

EgammaCoreTools.

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

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

Definition at line 75 of file PiZeroAnalyzer.h.

Constructor & Destructor Documentation

◆ PiZeroAnalyzer()

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

explicit

Definition at line 20 of file PiZeroAnalyzer.cc.

References deDxTools::esConsumes().

     : caloGeometryToken_{esConsumes()}, caloTopologyToken_{esConsumes()} {
   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::~PiZeroAnalyzer ( )

override

Definition at line 56 of file PiZeroAnalyzer.cc.

56 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 81 of file PiZeroAnalyzer.cc.

References barrelEcalHits_token_, HLT_2022v15_cff::barrelRecHits, MillePedeFileConverter_cfg::e, endcapEcalHits_token_, HLT_2022v15_cff::endcapRecHits, edm::HandleBase::isValid(), makePizero(), nEvt_, and prescaleFactor_.

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

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 58 of file PiZeroAnalyzer.cc.

References dqm::implementation::IBooker::book1D(), currentFolder_, hIsoPi0EB_, hMinvPi0EB_, hPt1Pi0EB_, hPt2Pi0EB_, hPtPi0EB_, dqm::impl::MonitorElement::setAxisTitle(), and dqm::implementation::NavigatorBase::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);
 }

◆ makePizero()

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

private

Definition at line 112 of file PiZeroAnalyzer.cc.

References edm::SortedCollection< T, SORT >::begin(), PositionCalc::Calculate_Location(), caloGeometryToken_, caloTopologyToken_, clusEtaSize_, clusPhiSize_, clusSeedThr_, funct::cos(), gather_cfg::cout, DetId::Ecal, EcalBarrel, EcalPreshower, edm::SortedCollection< T, SORT >::end(), HCALHighEnergyHPDFilter_cfi::energy, JetChargeProducer_cfi::exp, dqm::impl::MonitorElement::Fill(), edm::EventSetup::getHandle(), CaloSubdetectorTopology::getWindow(), hIsoPi0EB_, hMinvPi0EB_, hPt1Pi0EB_, hPt2Pi0EB_, hPtPi0EB_, mps_fire::i, l1ctLayer2EG_cff::id, EBDetId::ieta(), LEDCalibrationChannels::ieta, EBDetId::iphi(), LEDCalibrationChannels::iphi, dqmiolumiharvest::j, dqmdumpme::k, posCalcParameters_, edm::Handle< T >::product(), DetachedQuadStep_cff::seeds, seleMinvMaxPi0_, seleMinvMinPi0_, selePi0BeltDeta_, selePi0BeltDR_, selePi0Iso_, selePtGammaOne_, selePtGammaTwo_, selePtPi0_, seleS4S9GammaOne_, seleS4S9GammaTwo_, seleXtalMinEnergy_, funct::sin(), jetUpdater_cfi::sort, mathSSE::sqrt(), x, and y.

Referenced by analyze().

                                                                             {
   const EcalRecHitCollection* hitCollection_p = rhEB.product();
 
   //edm::ESHandle<CaloGeometry> geoHandle;
   //es.get<CaloGeometryRecord>().get(geoHandle);
   auto geoHandle = es.getHandle(caloGeometryToken_);
 
   //edm::ESHandle<CaloTopology> theCaloTopology;
   //es.get<CaloTopologyRecord>().get(theCaloTopology);
   auto theCaloTopology = es.getHandle(caloTopologyToken_);
 
   const CaloSubdetectorTopology* topology_p;
   const CaloSubdetectorGeometry* geometry_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
   const CaloSubdetectorGeometry* 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(), [](auto& x, auto& y) { return (x.energy() > y.energy()); });
   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

◆ barrelEcalHits_token_

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

private

Definition at line 94 of file PiZeroAnalyzer.h.

Referenced by analyze().

◆ caloGeometryToken_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> PiZeroAnalyzer::caloGeometryToken_

private

Definition at line 97 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ caloTopologyToken_

const edm::ESGetToken<CaloTopology, CaloTopologyRecord> PiZeroAnalyzer::caloTopologyToken_

private

Definition at line 98 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ clusEtaSize_

int PiZeroAnalyzer::clusEtaSize_

private

Definition at line 107 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ clusPhiSize_

int PiZeroAnalyzer::clusPhiSize_

private

Definition at line 108 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ clusSeedThr_

double PiZeroAnalyzer::clusSeedThr_

private

parameters needed for pizero finding

Definition at line 106 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ currentFolder_

std::stringstream PiZeroAnalyzer::currentFolder_

private

Definition at line 128 of file PiZeroAnalyzer.h.

Referenced by bookHistograms().

◆ cutStep_

double PiZeroAnalyzer::cutStep_

private

Definition at line 102 of file PiZeroAnalyzer.h.

◆ endcapEcalHits_token_

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

private

Definition at line 95 of file PiZeroAnalyzer.h.

Referenced by analyze().

◆ fName_

std::string PiZeroAnalyzer::fName_

private

Definition at line 87 of file PiZeroAnalyzer.h.

◆ hIsoPi0EB_

MonitorElement* PiZeroAnalyzer::hIsoPi0EB_

private

Definition at line 133 of file PiZeroAnalyzer.h.

Referenced by bookHistograms(), and makePizero().

◆ hMinvPi0EB_

MonitorElement* PiZeroAnalyzer::hMinvPi0EB_

private

Definition at line 130 of file PiZeroAnalyzer.h.

Referenced by bookHistograms(), and makePizero().

◆ hPt1Pi0EB_

MonitorElement* PiZeroAnalyzer::hPt1Pi0EB_

private

Definition at line 131 of file PiZeroAnalyzer.h.

Referenced by bookHistograms(), and makePizero().

◆ hPt2Pi0EB_

MonitorElement* PiZeroAnalyzer::hPt2Pi0EB_

private

Definition at line 132 of file PiZeroAnalyzer.h.

Referenced by bookHistograms(), and makePizero().

◆ hPtPi0EB_

MonitorElement* PiZeroAnalyzer::hPtPi0EB_

private

Definition at line 134 of file PiZeroAnalyzer.h.

Referenced by bookHistograms(), and makePizero().

◆ minPhoEtCut_

double PiZeroAnalyzer::minPhoEtCut_

private

Definition at line 100 of file PiZeroAnalyzer.h.

◆ nEvt_

int PiZeroAnalyzer::nEvt_

private

Definition at line 90 of file PiZeroAnalyzer.h.

Referenced by analyze().

◆ numberOfSteps_

int PiZeroAnalyzer::numberOfSteps_

private

Definition at line 103 of file PiZeroAnalyzer.h.

◆ ParameterLogWeighted_

bool PiZeroAnalyzer::ParameterLogWeighted_

private

Definition at line 112 of file PiZeroAnalyzer.h.

◆ ParameterT0_barl_

double PiZeroAnalyzer::ParameterT0_barl_

private

Definition at line 114 of file PiZeroAnalyzer.h.

◆ ParameterW0_

double PiZeroAnalyzer::ParameterW0_

private

Definition at line 115 of file PiZeroAnalyzer.h.

◆ ParameterX0_

double PiZeroAnalyzer::ParameterX0_

private

Definition at line 113 of file PiZeroAnalyzer.h.

◆ posCalcParameters_

edm::ParameterSet PiZeroAnalyzer::posCalcParameters_

private

Definition at line 92 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ prescaleFactor_

unsigned int PiZeroAnalyzer::prescaleFactor_

private

Definition at line 88 of file PiZeroAnalyzer.h.

Referenced by analyze().

◆ seleMinvMaxPi0_

double PiZeroAnalyzer::seleMinvMaxPi0_

private

Definition at line 125 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ seleMinvMinPi0_

double PiZeroAnalyzer::seleMinvMinPi0_

private

Definition at line 126 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ selePi0BeltDeta_

double PiZeroAnalyzer::selePi0BeltDeta_

private

Definition at line 123 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ selePi0BeltDR_

double PiZeroAnalyzer::selePi0BeltDR_

private

Definition at line 122 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ selePi0Iso_

double PiZeroAnalyzer::selePi0Iso_

private

Definition at line 124 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ selePtGammaOne_

double PiZeroAnalyzer::selePtGammaOne_

private

Definition at line 117 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ selePtGammaTwo_

double PiZeroAnalyzer::selePtGammaTwo_

private

Definition at line 118 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ selePtPi0_

double PiZeroAnalyzer::selePtPi0_

private

Definition at line 119 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ seleS4S9GammaOne_

double PiZeroAnalyzer::seleS4S9GammaOne_

private

Definition at line 120 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ seleS4S9GammaTwo_

double PiZeroAnalyzer::seleS4S9GammaTwo_

private

Definition at line 121 of file PiZeroAnalyzer.h.

Referenced by makePizero().

◆ seleXtalMinEnergy_

double PiZeroAnalyzer::seleXtalMinEnergy_

private

Definition at line 110 of file PiZeroAnalyzer.h.

Referenced by makePizero().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PiZeroAnalyzer()

◆ ~PiZeroAnalyzer()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ makePizero()

Member Data Documentation

◆ barrelEcalHits_token_

◆ caloGeometryToken_

◆ caloTopologyToken_

◆ clusEtaSize_

◆ clusPhiSize_

◆ clusSeedThr_

◆ currentFolder_

◆ cutStep_

◆ endcapEcalHits_token_

◆ fName_

◆ hIsoPi0EB_

◆ hMinvPi0EB_

◆ hPt1Pi0EB_

◆ hPt2Pi0EB_

◆ hPtPi0EB_

◆ minPhoEtCut_

◆ nEvt_

◆ numberOfSteps_

◆ ParameterLogWeighted_

◆ ParameterT0_barl_

◆ ParameterW0_

◆ ParameterX0_

◆ posCalcParameters_

◆ prescaleFactor_

◆ seleMinvMaxPi0_

◆ seleMinvMinPi0_

◆ selePi0BeltDeta_

◆ selePi0BeltDR_

◆ selePi0Iso_

◆ selePtGammaOne_

◆ selePtGammaTwo_

◆ selePtPi0_

◆ seleS4S9GammaOne_

◆ seleS4S9GammaTwo_

◆ seleXtalMinEnergy_