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 41 of file misc/Pi0EtaToGGRecAlg.cxx.

                                     {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "in execute()" << endmsg;
 
  // 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" << endmsg;
      return StatusCode::FAILURE;
    }
  }
 
  return StatusCode::SUCCESS;
}

finalize()

StatusCode Pi0EtaToGGRecAlg::finalize

(

)

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

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

initialize()

StatusCode Pi0EtaToGGRecAlg::initialize

(

)

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

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

---

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

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