---

Hector Class Reference

#include <SimTransport/HectorProducer/interface/Hector.h>

List of all members.

Public Member Functions
void	add (const HepMC::GenEvent *ev, const edm::EventSetup &es)
HepMC::GenEvent *	addPartToHepMC (HepMC::GenEvent *event)
void	clear ()
void	clearApertureFlags ()
void	filterD1 ()
void	filterFP420 ()
void	filterZDC ()
int	getDirect (unsigned int part_n) const
	Hector (const edm::ParameterSet &ps, bool verbosity, bool FP420Transport, bool ZDCTransport)
void	print () const
virtual	~Hector ()
Private Attributes
string	beam1filename
string	beam2filename
double	etacut
double	lengthd1
double	lengthfp420
double	lengthzdc
H_BeamLine *	m_beamlineD11
H_BeamLine *	m_beamlineD12
H_BeamLine *	m_beamlineFP4201
H_BeamLine *	m_beamlineFP4202
H_BeamLine *	m_beamlineZDC1
H_BeamLine *	m_beamlineZDC2
std::map< unsigned int, H_BeamParticle * >	m_beamPart
std::map< unsigned int, int >	m_direct
std::map< unsigned int, double >	m_eAtTrPoint
std::map< unsigned int, double >	m_eta
bool	m_FP420Transport
std::map< unsigned int, bool >	m_isCharged
std::map< unsigned int, bool >	m_isStoppedd1
std::map< unsigned int, bool >	m_isStoppedfp420
std::map< unsigned int, bool >	m_isStoppedzdc
std::map< unsigned int, int >	m_pdg
std::map< unsigned int, double >	m_pz
H_RecRPObject *	m_rp420_b
H_RecRPObject *	m_rp420_f
float	m_rpp420_b
float	m_rpp420_f
float	m_rppd1
float	m_rppzdc
double	m_sig_e
double	m_sigmaSTX
double	m_sigmaSTY
bool	m_smearAng
bool	m_smearE
std::map< unsigned int, double >	m_TxAtTrPoint
std::map< unsigned int, double >	m_TyAtTrPoint
bool	m_verbosity
std::map< unsigned int, double >	m_xAtTrPoint
std::map< unsigned int, double >	m_yAtTrPoint
bool	m_ZDCTransport
edm::ESHandle< ParticleDataTable >	pdt

---

Detailed Description

Definition at line 34 of file Hector.h.

---

Constructor & Destructor Documentation

Hector::Hector	(	const edm::ParameterSet &	ps,
		bool	verbosity,
		bool	FP420Transport,
		bool	ZDCTransport
	)

Definition at line 17 of file Hector.cc.

References b1, b2, beam1filename, beam2filename, GenMuonPlsPt100GeV_cfg::cout, e, lat::endl(), etacut, edm::ParameterSet::getParameter(), lengthd1, lengthfp420, lengthzdc, m_beamlineD11, m_beamlineD12, m_beamlineFP4201, m_beamlineFP4202, m_beamlineZDC1, m_beamlineZDC2, m_FP420Transport, m_rpp420_b, m_rpp420_f, m_rppd1, m_rppzdc, m_sig_e, m_sigmaSTX, m_sigmaSTY, m_smearAng, m_smearE, m_verbosity, m_ZDCTransport, alivecheck_mergeAndRegister::msg, and cms::Exception::what().

                                                                                                   : 
  m_verbosity(verbosity), 
  m_FP420Transport(FP420Transport),
  m_ZDCTransport(ZDCTransport)
{
  
  edm::LogInfo ("Hector") << "===================================================================\n"  
                          << "=== Hector Constructor start                       =====\n";
  
  // Create LHC beam line
  edm::ParameterSet hector_par = param.getParameter<edm::ParameterSet>("Hector");
  
  // User definitons
  lengthfp420    = hector_par.getParameter<double>("BeamLineLengthFP420" );
  m_rpp420_f     = (float) hector_par.getParameter<double>("RP420f");
  m_rpp420_b     = (float) hector_par.getParameter<double>("RP420b");
  lengthzdc      = hector_par.getParameter<double>("BeamLineLengthZDC" );
  lengthd1       = hector_par.getParameter<double>("BeamLineLengthD1" );
  beam1filename  = hector_par.getParameter<string>("Beam1");
  beam2filename  = hector_par.getParameter<string>("Beam2");  
  m_rppzdc       = (float) lengthzdc ;
  m_rppd1        = (float) lengthd1 ;
  m_smearAng     = hector_par.getParameter<bool>("smearAng");
  m_sigmaSTX     = hector_par.getParameter<double>("sigmaSTX" );
  m_sigmaSTY     = hector_par.getParameter<double>("sigmaSTY" );
  m_smearE       = hector_par.getParameter<bool>("smearEnergy");
  m_sig_e        = hector_par.getParameter<double>("sigmaEnergy");
  etacut         = hector_par.getParameter<double>("EtaCutForHector" );
  
  edm::LogInfo ("Hector") << "Hector parameters: \n" 
                          << "   lengthfp420:    " << lengthfp420 << "\n"
                          << "   m_rpp420_f:    " << m_rpp420_f << "\n"
                          << "   m_rpp420_b:    " << m_rpp420_b << "\n"
                          << "   lengthzdc:    " << lengthzdc << "\n"
                          << "   lengthd1:    " << lengthd1 << "\n";
  
  if(m_verbosity) {
    cout << "===================================================================" << endl;  
    cout << "=== Hector:    m_FP420Transport=" << m_FP420Transport <<"m_ZDCTransport " << m_ZDCTransport << endl;
    cout << "=== lengthfp420: " << lengthfp420 << endl;
    cout << "=== m_rpp420_f: " << m_rpp420_f << endl;
    cout << "=== m_rpp420_b: " << m_rpp420_b << endl;
    cout << "=== lengthzdc: " << lengthzdc << endl;
    cout << "=== lengthd1: " << lengthd1 << endl;
    cout << "===================================================================" << endl;
  }  
  edm::FileInPath b1(beam1filename.c_str());
  edm::FileInPath b2(beam2filename.c_str());
  
  // construct beam line for FP420:                                                                                           .
  if(m_FP420Transport && lengthfp420>0. ) {
    m_beamlineFP4201 = new H_BeamLine(  1, lengthfp420 + 0.1 ); // (direction, length)
    m_beamlineFP4202 = new H_BeamLine( -1, lengthfp420 + 0.1 ); //
    try {
      m_beamlineFP4201->fill( b1.fullPath(),  1, "IP5" );
      m_beamlineFP4202->fill( b2.fullPath(), -1, "IP5" );
    } catch ( const edm::Exception& e ) {
      std::string msg = e.what();
      msg += " caught in Hector... \nERROR: Could not locate SimTransport/HectorData data files.";
      edm::LogError ("DataNotFound") << msg;
    }
    m_beamlineFP4201->offsetElements( 120, -0.097 );
    m_beamlineFP4202->offsetElements( 120, +0.097 );
    m_beamlineFP4201->calcMatrix();
    m_beamlineFP4202->calcMatrix();
  }  
  else{
    cout << "=== Hector: WARNING: lengthfp420=  " << lengthfp420 << endl;
  }
  
  
  if (m_ZDCTransport && lengthzdc>0. && lengthd1>0.) {
    // construct beam line for ZDC:                                                                                           .
    m_beamlineZDC1 = new H_BeamLine(  1, lengthzdc + 0.1 ); // (direction, length)
    m_beamlineZDC2 = new H_BeamLine( -1, lengthzdc + 0.1 ); //
    try {
      m_beamlineZDC1->fill( b1.fullPath(),  1, "IP5" );
      m_beamlineZDC2->fill( b2.fullPath(), -1, "IP5" );
    } catch ( const edm::Exception& e ) {
      std::string msg = e.what();
      msg += " caught in Hector... \nERROR: Could not locate SimTransport/HectorData data files.";
      edm::LogError ("DataNotFound") << msg;
    }
    m_beamlineZDC1->offsetElements( 120, -0.097 );
    m_beamlineZDC2->offsetElements( 120, +0.097 );
    m_beamlineZDC1->calcMatrix();
    m_beamlineZDC2->calcMatrix();
    
    
    // construct beam line for D1:                                                                                           .
    m_beamlineD11 = new H_BeamLine(  1, lengthd1 + 0.1 ); // (direction, length)
    m_beamlineD12 = new H_BeamLine( -1, lengthd1 + 0.1 ); //
    try {
      m_beamlineD11->fill( b1.fullPath(),  1, "IP5" );
      m_beamlineD12->fill( b2.fullPath(), -1, "IP5" );
    } catch ( const edm::Exception& e ) {
      std::string msg = e.what();
      msg += " caught in Hector... \nERROR: Could not locate SimTransport/HectorData data files.";
      edm::LogError ("DataNotFound") << msg;
    }
    m_beamlineD11->offsetElements( 120, -0.097 );
    m_beamlineD12->offsetElements( 120, +0.097 );
    m_beamlineD11->calcMatrix();
    m_beamlineD12->calcMatrix();
  }  
  else{
    cout << "=== Hector: WARNING: lengthzdc=  " << lengthzdc << "lengthd1=  " << lengthd1 << endl;
  }
  
  edm::LogInfo ("Hector") << "===================================================================\n";
  
}

Hector::~Hector

(

)

[virtual]

Definition at line 130 of file Hector.cc.

References it, m_beamlineD11, m_beamlineD12, m_beamlineFP4201, m_beamlineFP4202, m_beamlineZDC1, m_beamlineZDC2, and m_beamPart.

               {
  
  edm::LogInfo ("Hector") << "===================================================================\n"  
                          << "=== Start delete Hector                                       =====\n";
  for (std::map<unsigned int,H_BeamParticle*>::iterator it = m_beamPart.begin(); it != m_beamPart.end(); ++it ) {
    delete (*it).second;
  }
  //
  delete m_beamlineFP4201;
  delete m_beamlineFP4202;
  delete m_beamlineZDC1;
  delete m_beamlineZDC2;
  delete m_beamlineD11;
  delete m_beamlineD12;
  
  edm::LogInfo ("Hector") << "===================================================================\n";  
}

---

Member Function Documentation

void Hector::add	(	const HepMC::GenEvent *	ev,
		const edm::EventSetup &	es
	)

Adds the stable protons from the event ev to a beamline

Definition at line 177 of file Hector.cc.

References funct::abs(), GenMuonPlsPt100GeV_cfg::cout, lat::endl(), etacut, g, edm::EventSetup::getData(), parsecf::pyparsing::line(), m_beamPart, m_direct, m_eta, m_isCharged, m_pdg, m_pz, m_verbosity, pdt, and funct::sqrt().

Referenced by HectorProducer::produce().

                                                                         {
  
  H_BeamParticle * h_p;
  unsigned int line;
  
  //  unsigned int npart = ev->nStable();
  
  /*  for (HepMC::GenEvent::particle_const_iterator eventParticle =evt->particles_begin();
      eventParticle != evt->particles_end();
      eventParticle = find_if(++eventParticle, evt->particles_end(), ( (*eventParticle)->status() == 1 ) ) ) {*/
  //    for ( HepMC::GenVertex::particle_iterator  eventParticle = (*vitr)->particles_begin(HepMC::descendants);
  for (HepMC::GenEvent::particle_const_iterator eventParticle =evt->particles_begin();
       eventParticle != evt->particles_end();
       ++eventParticle ) {
    //    if ( (*eventParticle)->status() == 1 && !(*eventParticle)->end_vertex()) {
    if ( (*eventParticle)->status() == 1 ) {
      //      if ( abs( (*eventParticle)->pdg_id() ) == 2212 ) { // if it's proton
      if ( abs( (*eventParticle)->momentum().eta())>etacut){
        line = (*eventParticle)->barcode();
        if ( m_beamPart.find(line) == m_beamPart.end() ) {
          double charge=1.;
          m_isCharged[line] = false;// neutrals
          HepMC::GenParticle * g = (*eventParticle);    
          iSetup.getData( pdt );
          const ParticleData * part = pdt->particle( g->pdg_id() );
          if (part){
            charge = part->charge();
            if(m_verbosity) cout <<"charge=  " << charge << endl;
          }
          if(charge !=0) m_isCharged[line] = true;//charged
          double mass = (*eventParticle)->generatedMass();
          if(m_verbosity) cout << "=== Hector:add:    status=  "<< g->status() <<"mass =  " << mass <<"charge=  " << charge << endl;

          //      h_p = new H_BeamParticle();

          //      H_BeamParticle myparticle(mass,charge);
          h_p = new H_BeamParticle(mass,charge);

          double px,py,pz,pt;
          px = (*eventParticle)->momentum().px();         
          py = (*eventParticle)->momentum().py();         
          pz = (*eventParticle)->momentum().pz();         
          
          h_p->set4Momentum( px, py, pz, (*eventParticle)->momentum().e() );
          
          pt = sqrt( (px*px) + (py*py) );
          
          // Clears H_BeamParticle::positions and sets the initial one
          h_p->setPosition( (*eventParticle)->production_vertex()->position().x(), (*eventParticle)->production_vertex()->position().y(), std::atan2( px, pt ), std::atan2( py, pt ), (*eventParticle)->production_vertex()->position().z() );
          


          m_beamPart[line] = h_p;
          m_direct[line] = 0;
          m_direct[line] = ( pz > 0 ) ? 1 : -1;
          
          m_eta[line] = (*eventParticle)->momentum().eta();
          m_pdg[line] = (*eventParticle)->pdg_id();
          m_pz[line]  = (*eventParticle)->momentum().pz();
          
          
          if(m_verbosity) {
            cout << "=== Hector:add:            pz=  "<< pz << endl;
            cout << "=== Hector:add:            pt=  "<<pt  << endl;
            cout << "=== Hector:add:   m_isCharged[line]=  "<< m_isCharged[line]  << endl;
          } 
        }// if find line
      }// if 2212 or eta 8.2
    }// if status
  }// for loop
  
}

HepMC::GenEvent * Hector::addPartToHepMC

(

HepMC::GenEvent *

event

)

Return vector of the particle lines (HepMC::GenParticle::barcode()) in a beamline

Definition at line 495 of file Hector.cc.

References funct::cos(), GenMuonPlsPt100GeV_cfg::cout, lat::endl(), relval_parameters_module::energy, it, lengthd1, parsecf::pyparsing::line(), m_beamPart, m_direct, m_eAtTrPoint, m_FP420Transport, m_isStoppedd1, m_isStoppedfp420, m_isStoppedzdc, m_rpp420_b, m_rpp420_f, m_TxAtTrPoint, m_TyAtTrPoint, m_verbosity, m_xAtTrPoint, m_yAtTrPoint, m_ZDCTransport, pi, edm::second(), funct::sin(), funct::sqrt(), and theta.

Referenced by HectorProducer::produce().

                                                           {
  
  unsigned int line;
  //  H_BeamParticle * part;
  HepMC::GenParticle * gpart;
  long double tx,ty,theta,fi,energy,time = 0;
  std::map< unsigned int, H_BeamParticle* >::iterator it;
  
  
  for (it = m_beamPart.begin(); it != m_beamPart.end(); ++it ) {
    line = (*it).first;
    if(!m_FP420Transport) m_isStoppedfp420[line] = true;
    if(!m_ZDCTransport) {m_isStoppedzdc[line] = false;m_isStoppedd1[line] = true;}
    if(m_verbosity) {
      cout << "=== Hector:addPartToHepMC: isStoppedfp420=" << (*m_isStoppedfp420.find(line)).second <<  endl;
      cout << "=== Hector:addPartToHepMC: isStoppedzdc=" << (*m_isStoppedzdc.find(line)).second <<  endl;
      cout << "=== Hector:addPartToHepMC: isStoppedd1=" << (*m_isStoppedd1.find(line)).second <<  endl;
    }
    if (!((*m_isStoppedfp420.find(line)).second) || (!((*m_isStoppedd1.find(line)).second) && ((*m_isStoppedzdc.find(line)).second))){
      //  if ( !((*m_isStoppedfp420.find(line)).second) || !((*m_isStoppedd1.find(line)).second)  ) {
      gpart = evt->barcode_to_particle( line );
      if ( gpart ) {
        tx     = (*m_TxAtTrPoint.find(line)).second / 1000000.;
        ty     = (*m_TyAtTrPoint.find(line)).second / 1000000.;
        theta  = sqrt((tx*tx) + (ty*ty));
        double ddd = 0.;
        if ( !((*m_isStoppedfp420.find(line)).second)) {
          if( (*m_direct.find( line )).second >0 ) {
            ddd = m_rpp420_f;
          }
          else if((*m_direct.find( line )).second <0 ) {
            ddd = m_rpp420_b;
            theta= pi-theta;
          }
        }
        else {
          ddd = lengthd1;
          if((*m_direct.find( line )).second <0 ) theta= pi-theta;
        }
        
        fi     = std::atan2(tx,ty); // tx, ty never == 0?
        energy = (*m_eAtTrPoint.find(line)).second;
        
        //      HepMC::GenEvent::vertex_iterator v_it;
        //      time   = ( *evt->vertices_begin() )->position().t(); // does time important?
        //long double time_buf = 0;
        //for (v_it = evt->vertices_begin(); v_it != evt->vertices_end(); ++v_it)  { // since no operator--
        //        time_buf = ( *v_it )->position().t();
        //  time = (  time > time_buf ) ? time : time_buf;
        //      }

        //      time = ( ddd*1000. - gpart->production_vertex()->position().z()*mm )/c_light;

        time = ( ddd*1000. - gpart->production_vertex()->position().z()*mm ); // mm
        //                               decaylength/(p.Beta()*p.Gamma()*c_light);
        //      time = ( ddd*1000. - gpart->production_vertex()->position().z()*mm )/300.; // nsec

        if(ddd != 0.) {
          if(m_verbosity) {
            std::cout<<"=========Hector:: x= "<< (*(m_xAtTrPoint.find(line))).second*0.001<<std::endl;
            std::cout<<"=========Hector:: y= "<< (*(m_yAtTrPoint.find(line))).second*0.001<<std::endl;
            std::cout<<"=========Hector:: z= "<< ddd * (*(m_direct.find( line ))).second*1000.<<std::endl;
            std::cout<<"=========Hector:: t= "<< time  <<std::endl;
          }
          
          
          
          HepMC::GenVertex * vert = new HepMC::GenVertex( HepMC::FourVector( (*(m_xAtTrPoint.find(line))).second*0.001,
                                                                             (*(m_yAtTrPoint.find(line))).second*0.001,
                                                                             ddd * (*(m_direct.find( line ))).second*1000.,
                                                                             time + .001*time ) );
          
          
          
          gpart->set_status( 2 );
          vert->add_particle_in( gpart );
          vert->add_particle_out( new HepMC::GenParticle( HepMC::FourVector(energy*std::sin(theta)*std::sin(fi),
                                                                            energy*std::sin(theta)*std::cos(fi),
                                                                            energy*std::cos(theta),
                                                                            energy ),
                                                          gpart->pdg_id(), 1, gpart->flow() ) );
          evt->add_vertex( vert );
          
        if(m_verbosity) std::cout << "Hector::TRANSPORTED pz= "<< gpart->momentum().pz()  << " eta= "<< gpart->momentum().eta()  << " status= "<< gpart->status()  <<std::endl;
          
          
        }// ddd
      }// if gpart
    }// if !isStopped
    
    else {
      gpart = evt->barcode_to_particle( line );
      if ( gpart ) {
        HepMC::GenVertex * vert= new HepMC::GenVertex();
        gpart->set_status( 2 );
        vert->add_particle_in( gpart );
        vert->add_particle_out( gpart );
        evt->add_vertex( vert );
        if(m_verbosity) std::cout << "Hector::NON-transp. pz= "<< gpart->momentum().pz()  << " eta= "<< gpart->momentum().eta()  << " status= "<< gpart->status()  <<std::endl;
      }
    }

  }//for 
  //  cout << "=== Hector:addPartToHepMC: end " << endl;
  
  //  if(m_verbosity){
  //   std::cout<<"============================================================================"<<std::endl;
  //  std::cout<<"=========Hector::addPartToHepMC print:"<<std::endl;
  //  evt->print();
  //}
  
  
  
  return evt;
} 

void Hector::clear

(

void

)

Clears BeamParticle, prepares Hector for a next Aperture check or/and a next event

Definition at line 154 of file Hector.cc.

References it, m_beamPart, m_direct, m_eta, m_isCharged, m_pdg, and m_pz.

Referenced by HectorProducer::produce().

                  {
  for ( std::map<unsigned int,H_BeamParticle*>::iterator it = m_beamPart.begin(); it != m_beamPart.end(); ++it ) {
    delete (*it).second;
  };
  m_beamPart.clear();
  m_direct.clear();
  m_eta.clear();
  m_pdg.clear();
  m_pz.clear();
  m_isCharged.clear();  
}

void Hector::clearApertureFlags

(

)

Clears ApertureFlags, prepares Hector for a next event

Definition at line 148 of file Hector.cc.

References m_isStoppedd1, m_isStoppedfp420, and m_isStoppedzdc.

Referenced by HectorProducer::produce().

                               {
  m_isStoppedfp420.clear();
  m_isStoppedzdc.clear();
  m_isStoppedd1.clear();
}

void Hector::filterD1

(

)

propagate the particles through a beamline to ZDC

Definition at line 402 of file Hector.cc.

References GenMuonPlsPt100GeV_cfg::cout, direction, lat::endl(), it, lengthd1, parsecf::pyparsing::line(), m_beamlineD11, m_beamlineD12, m_beamPart, m_direct, m_eAtTrPoint, m_eta, m_isCharged, m_isStoppedd1, m_isStoppedzdc, m_pdg, m_pz, m_sig_e, m_sigmaSTX, m_sigmaSTY, m_smearAng, m_smearE, m_TxAtTrPoint, m_TyAtTrPoint, m_verbosity, m_xAtTrPoint, m_yAtTrPoint, and edm::second().

Referenced by HectorProducer::produce().

                     {
  unsigned int line;
  H_BeamParticle * part;
  std::map< unsigned int, H_BeamParticle* >::iterator it;
  
  bool is_stop_d1;
  int direction;
  
  float x1_d1;
  float y1_d1;
  
  if ( m_beamPart.size() && lengthd1>0.) {
    
    for (it = m_beamPart.begin(); it != m_beamPart.end(); ++it ) {
      line = (*it).first;
      part = (*it).second;
      if(m_verbosity) cout << "=== Hector:filterD1:  isStopped=" << (*m_isStoppedzdc.find(line)).second <<  endl;
      if ( ((*m_isStoppedzdc.find(line)).second) || !((*m_isCharged.find( line )).second) ) {
        
        if(m_verbosity) cout << "filterD1:pdg_id=" << (*m_pdg.find( line )).second << " momentum().eta=" << (*m_eta.find( line )).second<< " pz=" << (*m_pz.find( line )).second<< endl;
        
        direction = (*m_direct.find( line )).second;
        if(m_verbosity) cout << "filterD1:direction=" << direction << endl;
        if (m_smearAng) {
          if ( m_sigmaSTX>0. && m_sigmaSTY>0.) {
            part->smearAng(m_sigmaSTX,m_sigmaSTY);
          }
          else {
            part->smearAng(); 
          }
        }
        if (m_smearE) {
          if ( m_sig_e ) {
            part->smearE(m_sig_e);
          }
          else {
            part->smearE(); 
          }
        }
        if ( direction == 1 ) {
          part->computePath( m_beamlineD11, 1 );
          is_stop_d1 = part->stopped( m_beamlineD11 );
          m_isStoppedd1[line] = is_stop_d1;
          if(m_verbosity) cout << "=== Hector:filterD1:    pos         is_stop_d1=  "<< is_stop_d1 << endl;
        }
        else  if ( direction == -1 ){
          part->computePath( m_beamlineD12, -1 );
          is_stop_d1 = part->stopped( m_beamlineD12 );
          m_isStoppedd1[line] = is_stop_d1;
          if(m_verbosity) cout << "=== Hector:filterD1:    neg         is_stop_d1=  "<< is_stop_d1 << endl;
        }
        else {
          is_stop_d1 = true;
          m_isStoppedd1[line] = is_stop_d1;
          if(m_verbosity) cout << "=== Hector:filterD1:    0         is_stop_d1=  "<< endl;
        }
        //propagating
        if (!is_stop_d1 ) {
          part->propagate( lengthd1 );
          x1_d1 = part->getX();
          y1_d1 = part->getY();
          m_xAtTrPoint[line]  = x1_d1;
          m_yAtTrPoint[line]  = y1_d1;
          m_TxAtTrPoint[line] = part->getTX();
          m_TyAtTrPoint[line] = part->getTY();
          m_eAtTrPoint[line]  = part->getE();
        }
      }// if stopzdc
      else {
        m_isStoppedd1[line] = false;// not stopped in propagating to ZDC and therefore in  propagation to D1 too.
        if(m_verbosity) cout << "=== Hector:filterD1:  isStopped=" << (*m_isStoppedd1.find(line)).second <<  endl;
      }
      
    } // for (it = m_beamPart.begin(); it != m_beamPart.end(); it++ ) 
  } // if ( m_beamPart.size() )
  
}

void Hector::filterFP420

(

)

propagate the particles through a beamline to FP420

Definition at line 254 of file Hector.cc.

References GenMuonPlsPt100GeV_cfg::cout, direction, lat::endl(), it, lengthfp420, parsecf::pyparsing::line(), m_beamlineFP4201, m_beamlineFP4202, m_beamPart, m_direct, m_eAtTrPoint, m_eta, m_isCharged, m_isStoppedfp420, m_pdg, m_pz, m_rpp420_b, m_rpp420_f, m_sig_e, m_sigmaSTX, m_sigmaSTY, m_smearAng, m_smearE, m_TxAtTrPoint, m_TyAtTrPoint, m_verbosity, m_xAtTrPoint, m_yAtTrPoint, and edm::second().

Referenced by HectorProducer::produce().

                        {
  unsigned int line;
  H_BeamParticle * part;
  std::map< unsigned int, H_BeamParticle* >::iterator it;
  
  bool is_stop;
  int direction;
  
  float x1_420;
  float y1_420;
  
  if ( m_beamPart.size() && lengthfp420>0. ) {
    
    for (it = m_beamPart.begin(); it != m_beamPart.end(); ++it ) {
      line = (*it).first;
      part = (*it).second;
      
      if(m_verbosity) {
        cout << "filterFP420:pdg_id=" << (*m_pdg.find( line )).second << " momentum().eta=" << (*m_eta.find( line )).second<< " pz=" << (*m_pz.find( line )).second<< endl;
      }
      if ( (*m_isCharged.find( line )).second ) {
        direction = (*m_direct.find( line )).second;
        if (m_smearAng) {
          if ( m_sigmaSTX>0. && m_sigmaSTY>0.) {
            part->smearAng(m_sigmaSTX,m_sigmaSTY);
          }
          else {
            part->smearAng(); 
          }
        }
        if (m_smearE) {
          if ( m_sig_e ) {
            part->smearE(m_sig_e);
          }
          else {
            part->smearE(); 
          }
        }
        if ( direction == 1 ) {
          part->computePath( m_beamlineFP4201, 1 );
          is_stop = part->stopped( m_beamlineFP4201 );
          if(m_verbosity) cout << "=== Hector:filterFP420:   pos          is_stop=  "<< is_stop << endl;
        }
        else if ( direction == -1 ){
          part->computePath( m_beamlineFP4202, -1 );
          is_stop = part->stopped( m_beamlineFP4202 );
          if(m_verbosity) cout << "=== Hector:filterFP420:   neg          is_stop=  "<< is_stop << endl;
        }
        else {
          is_stop = true;
          if(m_verbosity) cout << "=== Hector:filterFP420:   0          is_stop=  "<< is_stop << endl;
        }
        
        //propagating
        m_isStoppedfp420[line] = is_stop;
        if(m_verbosity) cout << "=== Hector:filterFP420:  isStopped=" << (*m_isStoppedfp420.find(line)).second <<  endl;
        
        if (!is_stop) {
          if ( direction == 1 ) part->propagate( m_rpp420_f );
          if ( direction == -1 ) part->propagate( m_rpp420_b );
          x1_420 = part->getX();
          y1_420 = part->getY();
          if(m_verbosity) cout << "=== Hector:filterFP420: x=  "<< x1_420 <<" y= " << y1_420 << endl;
          
          m_xAtTrPoint[line]  = x1_420;
          m_yAtTrPoint[line]  = y1_420;
          m_TxAtTrPoint[line] = part->getTX();
          m_TyAtTrPoint[line] = part->getTY();
          m_eAtTrPoint[line]  = part->getE();
          
        }
      }// if isCharged
      else {
        m_isStoppedfp420[line] = true;// imply that neutral particles stopped to reach 420m
        if(m_verbosity) cout << "=== Hector:filterFP420:  isStopped=" << (*m_isStoppedfp420.find(line)).second <<  endl;
      }
      
    } // for (it = m_beamPart.begin(); it != m_beamPart.end(); it++ ) 
  } // if ( m_beamPart.size() )
  
}

void Hector::filterZDC

(

)

propagate the particles through a beamline to ZDC

Definition at line 336 of file Hector.cc.

References GenMuonPlsPt100GeV_cfg::cout, direction, lat::endl(), it, lengthzdc, parsecf::pyparsing::line(), m_beamlineZDC1, m_beamlineZDC2, m_beamPart, m_direct, m_eta, m_isCharged, m_isStoppedfp420, m_isStoppedzdc, m_pdg, m_pz, m_sig_e, m_sigmaSTX, m_sigmaSTY, m_smearAng, m_smearE, m_verbosity, and edm::second().

Referenced by HectorProducer::produce().

                      {
  unsigned int line;
  H_BeamParticle * part;
  std::map< unsigned int, H_BeamParticle* >::iterator it;
  
  bool is_stop_zdc;
  int direction;
  
  if ( m_beamPart.size() && lengthzdc>0. ) {
    
    for (it = m_beamPart.begin(); it != m_beamPart.end(); ++it ) {
      line = (*it).first;
      part = (*it).second;
      if(m_verbosity) cout << "=== Hector:filterZDC:  isStopped=" << (*m_isStoppedfp420.find(line)).second <<  endl;
      if ( ((*m_isStoppedfp420.find(line)).second) && ((*m_isCharged.find(line)).second) ) {
        
        if(m_verbosity) cout << "filterZDC:pdg_id=" << (*m_pdg.find( line )).second << " momentum().eta=" << (*m_eta.find( line )).second<< " pz=" << (*m_pz.find( line )).second<< endl;
        
        direction = (*m_direct.find( line )).second;
        if(m_verbosity) cout << "filterZDC:direction=" << direction << endl;
        if (m_smearAng) {
          if ( m_sigmaSTX>0. && m_sigmaSTY>0.) {
            part->smearAng(m_sigmaSTX,m_sigmaSTY);
          }
          else {
            part->smearAng(); 
          }
        }
        if (m_smearE) {
          if ( m_sig_e ) {
            part->smearE(m_sig_e);
          }
          else {
            part->smearE(); 
          }
        }
        if ( direction == 1 ) {
          part->computePath( m_beamlineZDC1, 1 );
          is_stop_zdc = part->stopped( m_beamlineZDC1 );
          m_isStoppedzdc[line] = is_stop_zdc;
          if(m_verbosity) cout << "=== Hector:filterZDC:     pos        is_stop_zdc=  "<< is_stop_zdc << endl;
        }
        else if ( direction == -1 ){
          part->computePath( m_beamlineZDC2, -1 );
          is_stop_zdc = part->stopped( m_beamlineZDC2 );
          m_isStoppedzdc[line] = is_stop_zdc;
          if(m_verbosity) cout << "=== Hector:filterZDC:    neg         is_stop_zdc=  "<< is_stop_zdc << endl;
        }
        else {
          m_isStoppedzdc[line] = true;
          if(m_verbosity) cout << "=== Hector:filterZDC:    0         is_stop_zdc=  "<< endl;
        }
      }// if stopfp420 charged particles
      else if ( ((*m_isCharged.find(line)).second) ){
        m_isStoppedzdc[line] = false;// not stopped in propagating to FP420 and therefore in  propagation to ZDC too.
        if(m_verbosity) cout << "=== Hector:filterZDC:  isStopped=" << (*m_isStoppedzdc.find(line)).second <<  endl;
      }
      else {
        m_isStoppedzdc[line] = true;// neutrals particles considered as stopped in propagating to ZDC
        if(m_verbosity) cout << "=== Hector:filterZDC:  isStopped=" << (*m_isStoppedzdc.find(line)).second <<  endl;
      }
      
    } // for (it = m_beamPart.begin(); it != m_beamPart.end(); it++ ) 
  } // if ( m_beamPart.size() )
  
}

int Hector::getDirect

(

unsigned int

part_n

)

const

Definition at line 480 of file Hector.cc.

References it, and m_direct.

                                                  {
  std::map<unsigned int, int>::const_iterator it = m_direct.find( part_n );
  if ( it != m_direct.end() ){
    return (*it).second;
  }
  return 0;
}

void Hector::print

(

void

)

const

Prints properties of all particles in a beamline

Definition at line 488 of file Hector.cc.

References it, and m_beamPart.

                         {
  for (std::map<unsigned int,H_BeamParticle*>::const_iterator it = m_beamPart.begin(); it != m_beamPart.end(); ++it ) {
    (*it).second->printProperties();
  };
}

---

Member Data Documentation

string Hector::beam1filename [private]

Definition at line 121 of file Hector.h.

Referenced by Hector().

string Hector::beam2filename [private]

Definition at line 122 of file Hector.h.

Referenced by Hector().

double Hector::etacut [private]

Definition at line 78 of file Hector.h.

Referenced by add(), and Hector().

double Hector::lengthd1 [private]

Definition at line 76 of file Hector.h.

Referenced by addPartToHepMC(), filterD1(), and Hector().

double Hector::lengthfp420 [private]

Definition at line 74 of file Hector.h.

Referenced by filterFP420(), and Hector().

double Hector::lengthzdc [private]

Definition at line 75 of file Hector.h.

Referenced by filterZDC(), and Hector().

H_BeamLine* Hector::m_beamlineD11 [private]

Definition at line 98 of file Hector.h.

Referenced by filterD1(), Hector(), and ~Hector().

H_BeamLine* Hector::m_beamlineD12 [private]

Definition at line 99 of file Hector.h.

Referenced by filterD1(), Hector(), and ~Hector().

H_BeamLine* Hector::m_beamlineFP4201 [private]

Definition at line 94 of file Hector.h.

Referenced by filterFP420(), Hector(), and ~Hector().

H_BeamLine* Hector::m_beamlineFP4202 [private]

Definition at line 95 of file Hector.h.

Referenced by filterFP420(), Hector(), and ~Hector().

H_BeamLine* Hector::m_beamlineZDC1 [private]

Definition at line 96 of file Hector.h.

Referenced by filterZDC(), Hector(), and ~Hector().

H_BeamLine* Hector::m_beamlineZDC2 [private]

Definition at line 97 of file Hector.h.

Referenced by filterZDC(), Hector(), and ~Hector().

std::map<unsigned int, H_BeamParticle*> Hector::m_beamPart [private]

Definition at line 105 of file Hector.h.

Referenced by add(), addPartToHepMC(), clear(), filterD1(), filterFP420(), filterZDC(), print(), and ~Hector().

std::map<unsigned int, int> Hector::m_direct [private]

Definition at line 106 of file Hector.h.

Referenced by add(), addPartToHepMC(), clear(), filterD1(), filterFP420(), filterZDC(), and getDirect().

std::map<unsigned int, double> Hector::m_eAtTrPoint [private]

Definition at line 114 of file Hector.h.

Referenced by addPartToHepMC(), filterD1(), and filterFP420().

std::map<unsigned int, double> Hector::m_eta [private]

Definition at line 116 of file Hector.h.

Referenced by add(), clear(), filterD1(), filterFP420(), and filterZDC().

bool Hector::m_FP420Transport [private]

Definition at line 125 of file Hector.h.

Referenced by addPartToHepMC(), and Hector().

std::map<unsigned int, bool> Hector::m_isCharged [private]

Definition at line 119 of file Hector.h.

Referenced by add(), clear(), filterD1(), filterFP420(), and filterZDC().

std::map<unsigned int, bool> Hector::m_isStoppedd1 [private]

Definition at line 109 of file Hector.h.

Referenced by addPartToHepMC(), clearApertureFlags(), and filterD1().

std::map<unsigned int, bool> Hector::m_isStoppedfp420 [private]

Definition at line 107 of file Hector.h.

Referenced by addPartToHepMC(), clearApertureFlags(), filterFP420(), and filterZDC().

std::map<unsigned int, bool> Hector::m_isStoppedzdc [private]

Definition at line 108 of file Hector.h.

Referenced by addPartToHepMC(), clearApertureFlags(), filterD1(), and filterZDC().

std::map<unsigned int, int> Hector::m_pdg [private]

Definition at line 117 of file Hector.h.

Referenced by add(), clear(), filterD1(), filterFP420(), and filterZDC().

std::map<unsigned int, double> Hector::m_pz [private]

Definition at line 118 of file Hector.h.

Referenced by add(), clear(), filterD1(), filterFP420(), and filterZDC().

H_RecRPObject* Hector::m_rp420_b [private]

Definition at line 103 of file Hector.h.

H_RecRPObject* Hector::m_rp420_f [private]

Definition at line 102 of file Hector.h.

float Hector::m_rpp420_b [private]

Definition at line 86 of file Hector.h.

Referenced by addPartToHepMC(), filterFP420(), and Hector().

float Hector::m_rpp420_f [private]

Definition at line 85 of file Hector.h.

Referenced by addPartToHepMC(), filterFP420(), and Hector().

float Hector::m_rppd1 [private]

Definition at line 88 of file Hector.h.

Referenced by Hector().

float Hector::m_rppzdc [private]

Definition at line 87 of file Hector.h.

Referenced by Hector().

double Hector::m_sig_e [private]

Definition at line 80 of file Hector.h.

Referenced by filterD1(), filterFP420(), filterZDC(), and Hector().

double Hector::m_sigmaSTX [private]

Definition at line 82 of file Hector.h.

Referenced by filterD1(), filterFP420(), filterZDC(), and Hector().

double Hector::m_sigmaSTY [private]

Definition at line 83 of file Hector.h.

Referenced by filterD1(), filterFP420(), filterZDC(), and Hector().

bool Hector::m_smearAng [private]

Definition at line 79 of file Hector.h.

Referenced by filterD1(), filterFP420(), filterZDC(), and Hector().

bool Hector::m_smearE [private]

Definition at line 81 of file Hector.h.

Referenced by filterD1(), filterFP420(), filterZDC(), and Hector().

std::map<unsigned int, double> Hector::m_TxAtTrPoint [private]

Definition at line 112 of file Hector.h.

Referenced by addPartToHepMC(), filterD1(), and filterFP420().

std::map<unsigned int, double> Hector::m_TyAtTrPoint [private]

Definition at line 113 of file Hector.h.

Referenced by addPartToHepMC(), filterD1(), and filterFP420().

bool Hector::m_verbosity [private]

Definition at line 124 of file Hector.h.

Referenced by add(), addPartToHepMC(), filterD1(), filterFP420(), filterZDC(), and Hector().

std::map<unsigned int, double> Hector::m_xAtTrPoint [private]

Definition at line 110 of file Hector.h.

Referenced by addPartToHepMC(), filterD1(), and filterFP420().

std::map<unsigned int, double> Hector::m_yAtTrPoint [private]

Definition at line 111 of file Hector.h.

Referenced by addPartToHepMC(), filterD1(), and filterFP420().

bool Hector::m_ZDCTransport [private]

Definition at line 126 of file Hector.h.

Referenced by addPartToHepMC(), and Hector().

edm::ESHandle< ParticleDataTable > Hector::pdt [private]

Definition at line 91 of file Hector.h.

Referenced by add().

---

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

• SimTransport/HectorProducer/interface/Hector.h
• SimTransport/HectorProducer/src/Hector.cc

---