Pi0EtaToGGRecAlg Class Reference

#include <Pi0EtaToGGRecAlg.h>

Inheritance diagram for Pi0EtaToGGRecAlg:

Public Member Functions
	Pi0EtaToGGRecAlg (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()

Detailed Description

Definition at line 16 of file Pi0EtaToGGRecAlg.h.

Constructor & Destructor Documentation

Pi0EtaToGGRecAlg()

Pi0EtaToGGRecAlg::Pi0EtaToGGRecAlg	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Individual photons

Definition at line 15 of file misc/Pi0EtaToGGRecAlg.cxx.

                                                                                  :
   Algorithm(name, pSvcLocator)
{
   //Declare the properties
   declareProperty("Pi0MinMass",       m_minmass_cut  = 0.07);
   declareProperty("Pi0MaxMass",       m_maxmass_cut  = 0.18);
   declareProperty("Pi0ChisqCut",      m_chisq_cut    = 20.0);
 
   /// Individual photons
   declareProperty("PhotonUseBarrelEmc",    m_useBarrel        = true);
   declareProperty("PhotonUseEndcapEmc",    m_useEndcap        = true);
   declareProperty("PhotonMinEnergy",       m_minEnergy        = 0.02);
   declareProperty("PhotonDeltaAngleCut",   m_angled_cut       = 20.0);
   declareProperty("PhotonDeltaThetaCut",   m_thed_cut         = 20.0);
   declareProperty("PhotonDeltaPhiCut",     m_phid_cut         = 20.0);
 
}

Member Function Documentation

execute()

StatusCode Pi0EtaToGGRecAlg::execute

(

)

Perform pi0 1C kinfit

Fill number of pi0s in EvtRecEvent

Definition at line 44 of file misc/Pi0EtaToGGRecAlg.cxx.

                                     {
 
   MsgStream log(msgSvc(), name());
   log << MSG::INFO << "in execute()" <<endreq;
 
   // get event header, eventlist and tracklist from TDS
   SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
   SmartDataPtr<EvtRecEvent> recEvt(eventSvc(), EventModel::EvtRec::EvtRecEvent);
   SmartDataPtr<EvtRecTrackCol> recTrkCol(eventSvc(),  EventModel::EvtRec::EvtRecTrackCol);
 
   bool save2TDS = false;
   SmartDataPtr<EvtRecPi0Col> recPi0Col(eventSvc(), EventModel::EvtRec::EvtRecPi0Col);
   if ( !recPi0Col ) {
      recPi0Col = new EvtRecPi0Col;
      save2TDS = true;
   }
 
   KinematicFit * kmfit = KinematicFit::instance();
 
   //
   //  Set up the pi0's list
   //
   for(int i1 = recEvt->totalCharged(); i1 < (recEvt->totalTracks()-1); i1++) {
      EvtRecTrack* g1Trk = *(recTrkCol->begin() + i1);
      if ( !validPhoton(g1Trk,recEvt,recTrkCol) ) continue;
      RecEmcShower* g1Shower = g1Trk->emcShower();
      HepLorentzVector g1P4 = getP4(g1Shower);
 
      for(int i2 = i1+1; i2 < recEvt->totalTracks(); i2++) {
     EvtRecTrack* g2Trk = *(recTrkCol->begin() + i2);
     if ( !validPhoton(g2Trk,recEvt,recTrkCol) ) continue;
     RecEmcShower* g2Shower = g2Trk->emcShower();
     HepLorentzVector g2P4 = getP4(g2Shower);
 
     HepLorentzVector p2g = g1P4 + g2P4;
     if ( p2g.m() < m_minmass_cut ) continue;
     if ( p2g.m() > m_maxmass_cut ) continue;
 
     /// Perform pi0 1C kinfit 
     kmfit->init();
     kmfit->setIterNumber(1);
     kmfit->AddTrack(0, 0.0, g1Shower);
     kmfit->AddTrack(1, 0.0, g2Shower);
     kmfit->AddResonance(0, xmpi0, 0, 1);
     if ( !kmfit->Fit(0) ) continue;
     kmfit->BuildVirtualParticle(0);
 
     double pi0_chisq = kmfit->chisq(0);
     if ( pi0_chisq > m_chisq_cut ) continue;
 
     double e1 = (kmfit->pfit(0)).e();
     double e2 = (kmfit->pfit(1)).e();
     HepLorentzVector ppi0 = kmfit->pfit(0) + kmfit->pfit(1);
 
     EvtRecPi0* recPi0 = new EvtRecPi0();
     recPi0->setRawMass(p2g.restMass());
     recPi0->setChisq(pi0_chisq);
     recPi0->setFcos((e1-e2)/ppi0.rho());
 
     if ( g1P4.e() >= g2P4.e() ) {
        recPi0->setHiPfit(kmfit->pfit(0));
        recPi0->setLoPfit(kmfit->pfit(1));
        recPi0->setHiEnGamma(g1Trk);
        recPi0->setLoEnGamma(g2Trk);
     }
     else {
        recPi0->setHiPfit(kmfit->pfit(1));
        recPi0->setLoPfit(kmfit->pfit(0));
        recPi0->setHiEnGamma(g2Trk);
        recPi0->setLoEnGamma(g1Trk);
     }
 
     std::cout << "p(g1) = " <<  g1P4 << std::endl; 
     std::cout << "p(g2) = " <<  g2P4 << std::endl; 
     std::cout << "chi2(pi0) = " << pi0_chisq << std::endl; 
     std::cout << "m(gg) before fill = " << p2g.restMass() << std::endl << std::endl;
 
     recPi0Col->push_back(recPi0);
      }
   }
 
   /// Fill number of pi0s in EvtRecEvent
   recEvt->setNumberOfPi0(recPi0Col->size());
 
   if ( save2TDS ) {
      StatusCode sc = eventSvc()->registerObject(EventModel::EvtRec::EvtRecPi0Col, recPi0Col);
      if ( sc.isFailure() ) {
     log << MSG::ERROR << "could not register EvtRecPi0Col in TDS" <<endreq;
     return StatusCode::FAILURE;
      }
   }
 
   return StatusCode::SUCCESS;
}

finalize()

StatusCode Pi0EtaToGGRecAlg::finalize

(

)

Definition at line 141 of file misc/Pi0EtaToGGRecAlg.cxx.

                                      {
 
   MsgStream log(msgSvc(), name());
   log << MSG::INFO << "in finalize()" << endreq;
 
   return StatusCode::SUCCESS;
}

initialize()

StatusCode Pi0EtaToGGRecAlg::initialize

(

)

Definition at line 34 of file misc/Pi0EtaToGGRecAlg.cxx.

                                        {
 
   MsgStream log(msgSvc(), name());
   log << MSG::INFO << "in initialize()" <<endreq;
 
   return StatusCode::SUCCESS;
}

---

The documentation for this class was generated from the following files:

• Pi0EtaToGGRecAlg.h
• misc/Pi0EtaToGGRecAlg.cxx
• Pi0EtaToGGRecAlg.cxx