#include <Eepipi.h>

Inheritance diagram for Eepipi:

Public Member Functions
	Eepipi (const string &name, ISvcLocator *pSvcLocator)

StatusCode	initialize ()

StatusCode	execute ()

StatusCode	finalize ()

StatusCode	getMaxEvent ()

	Eepipi (const string &name, ISvcLocator *pSvcLocator)

StatusCode	initialize ()

StatusCode	execute ()

StatusCode	finalize ()

StatusCode	getMaxEvent ()

Detailed Description

Definition at line 25 of file Generator/Eepipi/Eepipi-00-00-06/Eepipi/Eepipi.h.

Constructor & Destructor Documentation

◆ Eepipi() [1/2]

Eepipi::Eepipi	(	const string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 90 of file Eepipi/Eepipi-00-00-06/src/Eepipi.cxx.

                                                          :Algorithm( name, pSvcLocator )
{
  declareProperty("Ecms", m_Ecms = 3.65); // Ecm = sqrt(s) [GeV]
//  declareProperty("cosee", m_cosee = 0.99); // set costheta for e+ e- in final state 
  declareProperty("WriteMC", m_mctruth = 0); //write out the MC truth of final state to "pdata1.dat"
}

◆ Eepipi() [2/2]

Eepipi::Eepipi	(	const string &	name,
		ISvcLocator *	pSvcLocator
	)

Member Function Documentation

◆ execute() [1/2]

StatusCode Eepipi::execute ( )

Definition at line 129 of file Eepipi/Eepipi-00-00-06/src/Eepipi.cxx.

{ 
  MsgStream log(messageService(), name());
 
 
  int npart = 0;
  int pid1,pid2,pid3,pid4,pst1,pst2;
  pid1 = 11;
  pid2 =-11;
  pid3 = 211;
  pid4 =-211;
  pst1 = 1;
  pst2 = 1;
 
 
 
  // Fill event information
  GenEvent* evt = new GenEvent(1,1);
 
  GenVertex* prod_vtx = new GenVertex();
//  prod_vtx->add_particle_out( p );
  evt->add_vertex( prod_vtx );
              
  // incoming beam e- HepLorentzVector(px,py,pz,energy) required!
  GenParticle* p = new GenParticle( CLHEP::HepLorentzVector( MOMSET.p1[1][evtgen],  MOMSET.p1[2][evtgen], 
                                      MOMSET.p1[3][evtgen],  MOMSET.p1[0][evtgen]),
                              11, 3); 
  p->suggest_barcode( ++npart );
  prod_vtx->add_particle_in(p);
 
//  std::cout<<"incoming beam e+"<<endl;
  // incoming beam e+, e+ moving along the z-direction
  p = new GenParticle( CLHEP::HepLorentzVector( MOMSET.q1[1][evtgen],  MOMSET.q1[2][evtgen],
                                         MOMSET.q1[3][evtgen],  MOMSET.q1[4][evtgen]),
               -11, 3); 
 
  p->suggest_barcode( ++npart );
  prod_vtx->add_particle_in(p);
  
  // scattered lepton +
  p = new GenParticle( CLHEP::HepLorentzVector( MOMSET.p2[1][evtgen],  MOMSET.p2[2][evtgen], 
                                 MOMSET.p2[3][evtgen],  MOMSET.p2[0][evtgen]),
                     pid1,pst1); 
  p->suggest_barcode( ++npart );
  prod_vtx->add_particle_out(p);
 
  // scattered lepton -
  p = new GenParticle( CLHEP::HepLorentzVector( MOMSET.q2[1][evtgen],  MOMSET.q2[2][evtgen],
                                         MOMSET.q2[3][evtgen],  MOMSET.q2[0][evtgen]),
               pid2, pst2); 
  p->suggest_barcode( ++npart );
  prod_vtx->add_particle_out(p);
 
  // outgoing pi+
  p = new GenParticle( CLHEP::HepLorentzVector( MOMSET.p3[1][evtgen],  MOMSET.p3[2][evtgen], 
                                 MOMSET.p3[3][evtgen],  MOMSET.p3[0][evtgen]),
              pid3,pst1); 
  p->suggest_barcode( ++npart );
  prod_vtx->add_particle_out(p);
 
  //outgoing pi-
  p = new GenParticle( CLHEP::HepLorentzVector( MOMSET.q3[1][evtgen],  MOMSET.q3[2][evtgen],
                                         MOMSET.q3[3][evtgen],  MOMSET.q3[0][evtgen]),
               pid4, pst2); 
  p->suggest_barcode( ++npart );
  prod_vtx->add_particle_out(p);
 
 
  
  evtgen++;
 
  if( log.level() < MSG::INFO )
  {
    evt->print();  
  }
 
  // Check if the McCollection already exists
  SmartDataPtr<McGenEventCol> anMcCol(eventSvc(), "/Event/Gen");
  if (anMcCol!=0) 
  {
    // Add event to existing collection
    MsgStream log(messageService(), name());
    log << MSG::INFO << "Add McGenEvent to existing collection" << endreq;
    McGenEvent* mcEvent = new McGenEvent(evt);
    anMcCol->push_back(mcEvent);
  }
  else 
  {
    // Create Collection and add  to the transient store
    McGenEventCol *mcColl = new McGenEventCol;
    McGenEvent* mcEvent = new McGenEvent(evt);
    mcColl->push_back(mcEvent);
    StatusCode sc = eventSvc()->registerObject("/Event/Gen",mcColl);
    if (sc != StatusCode::SUCCESS) 
    {
      log << MSG::ERROR << "Could not register McGenEvent" << endreq;
      delete mcColl;
      delete evt;
      delete mcEvent;
      return StatusCode::FAILURE;
    }
    else 
    {
      //   log << MSG::INFO << "McGenEventCol created and " << npart <<" particles stored in McGenEvent" << endreq;
    }
  }
     
  return StatusCode::SUCCESS; 
}

◆ execute() [2/2]

StatusCode Eepipi::execute ( )

◆ finalize() [1/2]

StatusCode Eepipi::finalize ( )

Definition at line 239 of file Eepipi/Eepipi-00-00-06/src/Eepipi.cxx.

{
  MsgStream log(messageService(), name());
   char delcmd[300];
   strcpy(delcmd,"cat ");
   strcat(delcmd,"fresult.dat");
   system(delcmd);
 
  std::cout<< "Eepipi finalized" << endl;
  return StatusCode::SUCCESS;
}

◆ finalize() [2/2]

StatusCode Eepipi::finalize ( )

◆ getMaxEvent() [1/2]

StatusCode Eepipi::getMaxEvent ( )

Definition at line 251 of file Eepipi/Eepipi-00-00-06/src/Eepipi.cxx.

                               {
  IProperty* appPropMgr=0;
  StatusCode status = 
    serviceLocator()->getService("ApplicationMgr", IProperty::interfaceID(),
                 reinterpret_cast<IInterface*&>( appPropMgr ));
  if( status.isFailure() ) return status;
    
  IntegerProperty evtMax("EvtMax",0);
  status = appPropMgr->getProperty( &evtMax );
  if (status.isFailure()) return status;
  
  m_evtMax = evtMax.value();
  return status;
}

Referenced by initialize().

◆ getMaxEvent() [2/2]