CTPPSHector Class Reference

#include <CTPPSHector.h>

Public Types
typedef CLHEP::HepLorentzVector	LorentzVector

Public Member Functions
void	add (const HepMC::GenEvent *ev, const edm::EventSetup &es, CLHEP::HepRandomEngine *engine)
HepMC::GenEvent *	addPartToHepMC (HepMC::GenEvent *event)
void	ApplyBeamCorrection (LorentzVector &, CLHEP::HepRandomEngine *engine)
void	clear ()
void	clearApertureFlags ()
	CTPPSHector (const edm::ParameterSet &ps, bool verbosity, bool CTPPSTransport)
void	filterCTPPS (TRandom3 *)
double	get_BeamEnergy ()
double	get_BeamMomentum ()
std::vector< LHCTransportLink > &	getCorrespondenceMap ()
int	getDirect (unsigned int part_n) const
void	LorentzBoost (LorentzVector &p_out, const string &frame)
void	print () const
void	set_BeamEnergy (double e)
virtual	~CTPPSHector ()

Private Attributes
string	beam1filename
string	beam2filename
double	etacut
double	fBeamEnergy
double	fBeamMomentum
bool	fCrossAngleCorr
double	fCrossingAngle
double	fMomentumMin
double	fVtxMeanX
double	fVtxMeanY
double	fVtxMeanZ
double	lengthctpps
float	m_b_ctpps_b
std::unique_ptr< H_BeamLine >	m_beamlineCTPPS1
std::unique_ptr< H_BeamLine >	m_beamlineCTPPS2
std::map< unsigned int, H_BeamParticle * >	m_beamPart
bool	m_CTPPSTransport
std::map< unsigned int, int >	m_direct
std::map< unsigned int, double >	m_eAtTrPoint
std::map< unsigned int, double >	m_eta
float	m_f_ctpps_f
std::map< unsigned int, bool >	m_isCharged
std::map< unsigned int, bool >	m_isStoppedctpps
std::map< unsigned int, int >	m_pdg
std::map< unsigned int, double >	m_pz
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
int	NEvent
edm::ESHandle< ParticleDataTable >	pdt
std::vector< LHCTransportLink >	theCorrespondenceMap

Detailed Description

Definition at line 39 of file CTPPSHector.h.

Member Typedef Documentation

typedef CLHEP::HepLorentzVector CTPPSHector::LorentzVector

Definition at line 46 of file CTPPSHector.h.

Constructor & Destructor Documentation

CTPPSHector::CTPPSHector	(	const edm::ParameterSet &	ps,
		bool	verbosity,
		bool	CTPPSTransport
	)

Definition at line 20 of file CTPPSHector.cc.

References beam1filename, beam2filename, etacut, fBeamEnergy, fCrossAngleCorr, fCrossingAngle, objects.autophobj::float, fMomentumMin, fVtxMeanX, fVtxMeanY, fVtxMeanZ, edm::ParameterSet::getParameter(), lengthctpps, LogDebug, m_b_ctpps_b, m_beamlineCTPPS1, m_beamlineCTPPS2, m_CTPPSTransport, m_f_ctpps_f, m_sig_e, m_sigmaSTX, m_sigmaSTY, m_smearAng, m_smearE, m_verbosity, and theCorrespondenceMap.

                                                                                          : 
    m_smearAng(false),m_sig_e(0.),m_smearE(false),m_sigmaSTX(0.),m_sigmaSTY(0.),
    fCrossAngleCorr(false),fCrossingAngle(0.),fBeamMomentum(0),fBeamEnergy(0),
    fVtxMeanX(0.),fVtxMeanY(0.),fVtxMeanZ(0.),fMomentumMin(0.),
    m_verbosity(verbosity), 
    m_CTPPSTransport(CTPPSTransport),NEvent(0)

{
    // Create LHC beam line
    edm::ParameterSet hector_par = param.getParameter<edm::ParameterSet>("CTPPSHector");

    // User definitons
    lengthctpps     = hector_par.getParameter<double>("BeamLineLengthCTPPS" );
    m_f_ctpps_f     = (float) hector_par.getParameter<double>("CTPPSf");
    m_b_ctpps_b     = (float) hector_par.getParameter<double>("CTPPSb");

    beam1filename   = hector_par.getParameter<string>("Beam1");
    beam2filename   = hector_par.getParameter<string>("Beam2");  
    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" );
    //CTPPS
    fCrossAngleCorr = hector_par.getParameter<bool>("CrossAngleCorr");
    fCrossingAngle  = hector_par.getParameter<double>("CrossingAngle");
    fBeamEnergy     = hector_par.getParameter<double>("BeamEnergy"); // beam energy in GeV
    fVtxMeanX       = hector_par.getParameter<double>("VtxMeanX");
    fVtxMeanY       = hector_par.getParameter<double>("VtxMeanY");
    fVtxMeanZ       = hector_par.getParameter<double>("VtxMeanZ");
    fMomentumMin    = hector_par.getParameter<double>("MomentumMin"); 

    theCorrespondenceMap.clear();

    if(m_verbosity) {
        edm::LogInfo("CTPPSHectorSetup") << "===================================================================\n"  
            << " * * * * * * * * * * * * * * * * * * * * * * * * * * * *           \n"  
            << " *                                                         *       \n"  
            << " *                   --<--<--  A fast simulator --<--<--     *     \n"  
            << " *                 | --<--<--     of particle   --<--<--     *     \n"  
            << " *  ----HECTOR----<                                          *     \n"  
            << " *                 | -->-->-- transport through-->-->--      *     \n"   
            << " *                   -->-->-- generic beamlines -->-->--     *     \n"  
            << " *                                                           *     \n"   
            << " * JINST 2:P09005 (2007)                                     *     \n"  
            << " *      X Rouby, J de Favereau, K Piotrzkowski (CP3)         *     \n"  
            << " *       http://www.fynu.ucl.ac.be/hector.html               *     \n"  
            << " *                                                           *     \n"  
            << " * Center for Cosmology, Particle Physics and Phenomenology  *     \n"  
            << " *              Universite catholique de Louvain             *     \n"  
            << " *                 Louvain-la-Neuve, Belgium                 *     \n"  
            << " *                                                         *       \n"  
            << " * * * * * * * * * * * * * * * * * * * * * * * * * * * *           \n"   
            << " CTPPSHector configuration: \n" 
            << " m_CTPPSTransport   = " << m_CTPPSTransport << "\n"
            << " lengthctpps      = " << lengthctpps << "\n"
            << " m_f_ctpps_f      =  " << m_f_ctpps_f << "\n"
            << " m_b_ctpps_b      =  " << m_b_ctpps_b << "\n"
            << "===================================================================\n";
    }  
    edm::FileInPath b1(beam1filename.c_str());
    edm::FileInPath b2(beam2filename.c_str());

    // construct beam line for CTPPS (forward 1 backward 2):                                                                                           
    if(m_CTPPSTransport && lengthctpps>0. ) {
        m_beamlineCTPPS1 = std::unique_ptr<H_BeamLine>(new H_BeamLine( -1, lengthctpps + 0.1 )); // (direction, length)
        m_beamlineCTPPS2 = std::unique_ptr<H_BeamLine>(new H_BeamLine( 1, lengthctpps + 0.1 )); //
        m_beamlineCTPPS1->fill( b2.fullPath(), 1, "IP5" );
        m_beamlineCTPPS2->fill( b1.fullPath(), 1, "IP5" );
        m_beamlineCTPPS1->offsetElements( 120, 0.097 );
        m_beamlineCTPPS2->offsetElements( 120, 0.097 );
        m_beamlineCTPPS1->calcMatrix();
        m_beamlineCTPPS2->calcMatrix();
    } else {
        if ( m_verbosity ) LogDebug("CTPPSHectorSetup") << "CTPPSHector: WARNING: lengthctpps=  " << lengthctpps;
    } 
}

CTPPSHector::~CTPPSHector ( )

virtual

Definition at line 99 of file CTPPSHector.cc.

References m_beamPart.

                         {

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

}

Member Function Documentation

void CTPPSHector::add	(	const HepMC::GenEvent *	ev,
		const edm::EventSetup &	es,
		CLHEP::HepRandomEngine *	engine
	)

Adds the stable protons from the event ev to a beamline

Definition at line 123 of file CTPPSHector.cc.

References funct::abs(), ApplyBeamCorrection(), ALCARECOTkAlJpsiMuMu_cff::charge, cm_to_m, cm_to_um, MillePedeFileConverter_cfg::e, etacut, fCrossAngleCorr, fMomentumMin, fVtxMeanX, fVtxMeanY, fVtxMeanZ, g, GenParticle::GenParticle, edm::EventSetup::getData(), geometryCSVtoXML::line, LogDebug, LorentzBoost(), m_beamPart, m_direct, m_eta, m_isCharged, m_pdg, m_pz, m_verbosity, ResonanceBuilder::mass, NULL, pdt, funct::pow(), and mathSSE::sqrt().

Referenced by CTPPSHectorProducer::produce(), and counter.Counter::register().

                                                                                                             {

    H_BeamParticle* h_p  = NULL;
    unsigned int line;

    for (HepMC::GenEvent::particle_const_iterator eventParticle =evt->particles_begin();
            eventParticle != evt->particles_end();
            ++eventParticle ) {
        if ( (*eventParticle)->status() == 1 && (*eventParticle)->pdg_id()==2212 ){
            if ( abs( (*eventParticle)->momentum().eta())>etacut && abs( (*eventParticle)->momentum().pz())>fMomentumMin){
                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(charge !=0) m_isCharged[line] = true;//charged
                    double mass = (*eventParticle)->generatedMass();

                    h_p = new H_BeamParticle(mass,charge);

                    double px,py,pz,e;
                    double TXforPosition=0.0, TYforPosition=0.0;//urad

                    px = (*eventParticle)->momentum().px();   
                    py = (*eventParticle)->momentum().py();   
                    pz = (*eventParticle)->momentum().pz();   

                    e = sqrt(pow(mass,2)+pow(px,2)+pow(py,2)+pow(pz,2));

                    // Apply Beam and Crossing Angle Corrections
                    LorentzVector p_out(px,py,pz,e);
                    ApplyBeamCorrection(p_out, engine);
                    if (fCrossAngleCorr) LorentzBoost(const_cast<LorentzVector&>(p_out),"LAB");

                    // from mm to cm        
                    double XforPosition = (*eventParticle)->production_vertex()->position().x()/cm;//cm
                    double YforPosition = (*eventParticle)->production_vertex()->position().y()/cm;//cm
                    double ZforPosition = (*eventParticle)->production_vertex()->position().z()/cm;//cm

                    if(m_verbosity) LogDebug("CTPPSHectorEventProcessing") << " fVtxMeanX: " << fVtxMeanX << " fVtxMeanY: " << fVtxMeanY << " fVtxMeanZ: "  << fVtxMeanZ ;
                    // It is important to set the Position before the 4Momentum otherwise HECTOR resets variables
                    h_p->setPosition(-(XforPosition-fVtxMeanX)*cm_to_um,(YforPosition-fVtxMeanY)*cm_to_um,TXforPosition,TYforPosition,-(ZforPosition)*cm_to_m);
                    
                    h_p->set4Momentum( -p_out.px(), p_out.py(), -p_out.pz(), p_out.e() );

                    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) { 
                        LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector:add: barcode = " << line 
                            << " status = " << g->status() 
                            << " PDG Id = " << g->pdg_id() 
                            << " mass = " << mass 
                            << " pz = " << pz 
                            << " charge = " << charge 
                            << " m_isCharged[line] = " << m_isCharged[line];
                    } 
                }// if find line
            }// if eta > 8.2
        }// if status
    }// for loop

}

HepMC::GenEvent * CTPPSHector::addPartToHepMC

(

HepMC::GenEvent *

event

)

Definition at line 345 of file CTPPSHector.cc.

References funct::cos(), GenParticle::GenParticle, geometryCSVtoXML::line, LogDebug, m_b_ctpps_b, m_beamPart, m_CTPPSTransport, m_direct, m_eAtTrPoint, m_f_ctpps_f, m_isStoppedctpps, m_TxAtTrPoint, m_TyAtTrPoint, m_verbosity, m_xAtTrPoint, m_yAtTrPoint, NEvent, pi, ProtonMassSQ, edm::second(), funct::sin(), mathSSE::sqrt(), theCorrespondenceMap, theta(), and ntuplemaker::time.

Referenced by get_BeamMomentum(), and CTPPSHectorProducer::produce().

                                                                {
    NEvent++;
    theCorrespondenceMap.clear();

    unsigned int line;

    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_CTPPSTransport) m_isStoppedctpps[line] = true;
        if(m_verbosity) {
            LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector:addPartToHepMC: barcode = " << line << "\n"
                << "CTPPSHector:addPartToHepMC: isStoppedctpps=" << (*m_isStoppedctpps.find(line)).second;
        }
        if (!((*m_isStoppedctpps.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.;
                long double fi_  = 0.; 
                if ( !((*m_isStoppedctpps.find(line)).second)) {
                    if( (*m_direct.find( line )).second >0 ) {
                        ddd = m_f_ctpps_f;
                        fi_    = std::atan2(tx,ty); // tx, ty never == 0?
                    }
                    else if((*m_direct.find( line )).second <0 ) {
                        ddd = m_b_ctpps_b;
                        theta= CLHEP::pi-theta;
                        fi_    = std::atan2(tx,ty); // tx, ty never == 0?
                    }
                } 
                fi = fi_; 
                energy = (*m_eAtTrPoint.find(line)).second;
                time = ( ddd*meter - gpart->production_vertex()->position().z()*mm ); // mm

                if(ddd != 0.) {
                    if(m_verbosity) {
                        LogDebug("CTPPSHectorEventProcessing") <<"CTPPSHector:: x= "<< (*(m_xAtTrPoint.find(line))).second*0.001<< "\n"
                            <<"CTPPSHector:: y= "<< (*(m_yAtTrPoint.find(line))).second*0.001<< "\n"
                            <<"CTPPSHector:: z= "<< ddd * (*(m_direct.find( line ))).second*1000. << "\n"
                            <<"CTPPSHector:: t= "<< time;
                    }

                    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 );
                    double Pmom = sqrt(energy*energy - ProtonMassSQ);
                    vert->add_particle_out( new HepMC::GenParticle( HepMC::FourVector(Pmom*std::sin(theta)*std::sin(fi),
                                    Pmom*std::sin(theta)*std::cos(fi),
                                    Pmom*std::cos(theta),
                                    energy ),gpart->pdg_id(), 1, gpart->flow() ) );
                    evt->add_vertex( vert );

                    int ingoing = (*vert->particles_in_const_begin())->barcode();
                    int outgoing = (*vert->particles_out_const_begin())->barcode();
                    LHCTransportLink theLink(ingoing,outgoing);
                    if (m_verbosity) LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector:addPartToHepMC: LHCTransportLink " << theLink;
                    theCorrespondenceMap.push_back(theLink);

                    if(m_verbosity) LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector::TRANSPORTED pz= " << gpart->momentum().pz()  
                        << " eta= "<< gpart->momentum().eta() << " status= "<< gpart->status();
                }// ddd
            }// if gpart
        }// if !isStopped

        else {
            gpart = evt->barcode_to_particle( line );
            if ( gpart ) {
                gpart->set_status( 2 );
                if(m_verbosity) LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector::NON-transp. pz= " << gpart->momentum().pz()  
                    << " eta= "<< gpart->momentum().eta()  << " status= "<< gpart->status();
            }
        }
    }//for 

    return evt;
} 

void CTPPSHector::ApplyBeamCorrection	(	LorentzVector &	p_out,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 281 of file CTPPSHector.cc.

References funct::cos(), MillePedeFileConverter_cfg::e, m_sig_e, m_sigmaSTX, m_sigmaSTY, m_smearAng, m_smearE, AlCaHLTBitMon_ParallelJobs::p, pi, ProtonMassSQ, funct::sin(), mathSSE::sqrt(), and theta().

Referenced by add(), and get_BeamMomentum().

{

    double microrad = 1.e-6;
    double theta = p_out.theta(); if (p_out.pz()<0) theta=CLHEP::pi-theta;
    double dtheta_x = (double)(m_smearAng)?CLHEP::RandGauss::shoot(engine,0.,m_sigmaSTX):0;
    double dtheta_y = (double)(m_smearAng)?CLHEP::RandGauss::shoot(engine,0.,m_sigmaSTY):0;
    double denergy  = (double)(m_smearE)?CLHEP::RandGauss::shoot(engine,0.,m_sig_e):0.;

    double p = sqrt((p_out.px())*(p_out.px())+(p_out.py())*(p_out.py())+(p_out.pz())*(p_out.pz()));
    double px = p*sin(theta+dtheta_x*microrad)*cos(p_out.phi());
    double py = p*sin(theta+dtheta_y*microrad)*sin(p_out.phi());
    double pz = p*(cos(theta)+denergy);

    if (p_out.pz()<0) pz*=-1;

    double e  = sqrt(px*px+py*py+pz*pz+ProtonMassSQ);
    p_out.setPx(px);
    p_out.setPy(py);
    p_out.setPz(pz);
    p_out.setE(e);

}

void CTPPSHector::clear

(

void

)

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

Definition at line 111 of file CTPPSHector.cc.

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

Referenced by CTPPSHectorProducer::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 CTPPSHector::clearApertureFlags

(

)

Clears ApertureFlags, prepares CTPPSHector for a next event

Definition at line 107 of file CTPPSHector.cc.

References m_isStoppedctpps.

Referenced by CTPPSHectorProducer::produce().

                                    {
    m_isStoppedctpps.clear();
}

void CTPPSHector::filterCTPPS

(

TRandom3 *

rootEngine

)

propagate the particles through a beamline to CTPPS

Definition at line 197 of file CTPPSHector.cc.

References lengthctpps, geometryCSVtoXML::line, LogDebug, m_b_ctpps_b, m_beamlineCTPPS1, m_beamlineCTPPS2, m_beamPart, m_direct, m_eAtTrPoint, m_f_ctpps_f, m_isCharged, m_isStoppedctpps, m_TxAtTrPoint, m_TyAtTrPoint, m_verbosity, m_xAtTrPoint, m_yAtTrPoint, NULL, and edm::second().

Referenced by CTPPSHectorProducer::produce().

                                                 {

    unsigned int line;
    H_BeamParticle * part = NULL;
 
    std::map< unsigned int, H_BeamParticle* >::iterator it;

    bool is_stop;
    int direction;

    float x1_ctpps;
    float y1_ctpps;

    if ( m_beamPart.size() && lengthctpps>0. ) {

        for (it = m_beamPart.begin(); it != m_beamPart.end(); ++it ) {
            line = (*it).first;
            part = (*it).second;

            if(m_verbosity) LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector:filterCTPPS: barcode = " << line;
            if ( (*m_isCharged.find( line )).second ) {
                direction = (*m_direct.find( line )).second;
                if ( direction == 1 && m_beamlineCTPPS1 != 0 ) {
                    
            part->computePath(&*m_beamlineCTPPS1);

                    is_stop = part->stopped(&* m_beamlineCTPPS1 );
                    if(m_verbosity) LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector:filterCTPPS: barcode = " << line << " positive is_stop=  "<< is_stop;
                }
                else if ( direction == -1 && m_beamlineCTPPS2 != 0 ){

                    part->computePath(&*m_beamlineCTPPS2 );

                    is_stop = part->stopped(&*m_beamlineCTPPS2 );
                    if(m_verbosity) LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector:filterCTPPS: barcode = " << line << " negative is_stop=  "<< is_stop;
                }
                else {
                    is_stop = true;
                    if(m_verbosity) LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector:filterCTPPS: barcode = " << line << " 0      is_stop=  "<< is_stop;
                }

                //propagating
                m_isStoppedctpps[line] = is_stop;
                if(m_verbosity) LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector:filterCTPPS: barcode = " << line << " isStopped=" << (*m_isStoppedctpps.find(line)).second;

                if (!is_stop) {
                    if ( direction == 1 ) part->propagate( m_f_ctpps_f ); 
                    if ( direction == -1 ) part->propagate( m_b_ctpps_b );  
                    x1_ctpps = -part->getX()/millimeter;
                    y1_ctpps = part->getY()/millimeter;
                    if(m_verbosity) LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector:filterCTPPS: barcode = " << line << " x=  "<< x1_ctpps <<" y= " << y1_ctpps;

                    m_xAtTrPoint[line]  = x1_ctpps;
                    m_yAtTrPoint[line]  = y1_ctpps;
                    m_TxAtTrPoint[line] = -part->getTX();
                    m_TyAtTrPoint[line] = part->getTY();
                    m_eAtTrPoint[line]  = part->getE();

                }
            }// if isCharged
            else {
                m_isStoppedctpps[line] = true;// imply that neutral particles stopped to reach 420m
                if(m_verbosity) LogDebug("CTPPSHectorEventProcessing") << "CTPPSHector:filterCTPPS: barcode = " << line << " isStopped=" << (*m_isStoppedctpps.find(line)).second;
            }

        } // for (it = m_beamPart.begin(); it != m_beamPart.end(); it++ ) 
    } // if ( m_beamPart.size() )

}//

double CTPPSHector::get_BeamEnergy ( )

inline

Definition at line 63 of file CTPPSHector.h.

References fBeamEnergy.

{return fBeamEnergy;};

double CTPPSHector::get_BeamMomentum ( )

inline

Definition at line 65 of file CTPPSHector.h.

References addPartToHepMC(), ApplyBeamCorrection(), fBeamMomentum, getDirect(), and print().

{return fBeamMomentum;};

std::vector<LHCTransportLink>& CTPPSHector::getCorrespondenceMap ( )

inline

Definition at line 76 of file CTPPSHector.h.

References theCorrespondenceMap.

Referenced by CTPPSHectorProducer::produce().

{ return theCorrespondenceMap; }

int CTPPSHector::getDirect

(

unsigned int

part_n

)

const

Definition at line 267 of file CTPPSHector.cc.

References m_direct.

Referenced by get_BeamMomentum().

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

void CTPPSHector::LorentzBoost	(	LorentzVector &	p_out,
		const string &	frame
	)

Definition at line 305 of file CTPPSHector.cc.

References funct::cos(), fCrossingAngle, p4, funct::sin(), and funct::tan().

Referenced by add().

{
    // Use a matrix
    double microrad = 1.e-6;
    TMatrixD tmpboost(4,4);
    double alpha_ = 0.;
    double phi_  = fCrossingAngle*microrad;
    if (p_out.pz()<0) phi_*=-1;
    tmpboost(0,0) = 1./cos(phi_);
    tmpboost(0,1) = - cos(alpha_)*sin(phi_);
    tmpboost(0,2) = - tan(phi_)*sin(phi_);
    tmpboost(0,3) = - sin(alpha_)*sin(phi_);
    tmpboost(1,0) = - cos(alpha_)*tan(phi_);
    tmpboost(1,1) = 1.;
    tmpboost(1,2) = cos(alpha_)*tan(phi_);
    tmpboost(1,3) = 0.;
    tmpboost(2,0) = 0.;
    tmpboost(2,1) = - cos(alpha_)*sin(phi_);
    tmpboost(2,2) = cos(phi_);
    tmpboost(2,3) = - sin(alpha_)*sin(phi_);
    tmpboost(3,0) = - sin(alpha_)*tan(phi_);
    tmpboost(3,1) = 0.;
    tmpboost(3,2) = sin(alpha_)*tan(phi_);
    tmpboost(3,3) = 1.;

    if(frame=="LAB") tmpboost.Invert();

    TMatrixD p4(4,1);
    p4(0,0) = p_out.e();
    p4(1,0) = p_out.px();
    p4(2,0) = p_out.py();
    p4(3,0) = p_out.pz();
    TMatrixD p4lab(4,1);
    p4lab = tmpboost * p4;
    p_out.setPx(p4lab(1,0));
    p_out.setPy(p4lab(2,0));
    p_out.setPz(p4lab(3,0));
    p_out.setE(p4lab(0,0));
}

void CTPPSHector::print

(

void

)

const

Prints properties of all particles in a beamline

Definition at line 275 of file CTPPSHector.cc.

References m_beamPart.

Referenced by get_BeamMomentum().

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

void CTPPSHector::set_BeamEnergy ( double  e )

inline

Definition at line 61 of file CTPPSHector.h.

References MillePedeFileConverter_cfg::e, fBeamEnergy, fBeamMomentum, ProtonMassSQ, and mathSSE::sqrt().

{fBeamEnergy=e;fBeamMomentum = sqrt(fBeamEnergy*fBeamEnergy - ProtonMassSQ);};

Member Data Documentation

string CTPPSHector::beam1filename

private

Definition at line 121 of file CTPPSHector.h.

Referenced by CTPPSHector().

string CTPPSHector::beam2filename

private

Definition at line 122 of file CTPPSHector.h.

Referenced by CTPPSHector().

double CTPPSHector::etacut

private

Definition at line 82 of file CTPPSHector.h.

Referenced by add(), and CTPPSHector().

double CTPPSHector::fBeamEnergy

private

Definition at line 95 of file CTPPSHector.h.

Referenced by CTPPSHector(), get_BeamEnergy(), and set_BeamEnergy().

double CTPPSHector::fBeamMomentum

private

Definition at line 94 of file CTPPSHector.h.

Referenced by get_BeamMomentum(), and set_BeamEnergy().

bool CTPPSHector::fCrossAngleCorr

private

Definition at line 92 of file CTPPSHector.h.

Referenced by add(), and CTPPSHector().

double CTPPSHector::fCrossingAngle

private

Definition at line 93 of file CTPPSHector.h.

Referenced by CTPPSHector(), and LorentzBoost().

double CTPPSHector::fMomentumMin

private

Definition at line 99 of file CTPPSHector.h.

Referenced by add(), and CTPPSHector().

double CTPPSHector::fVtxMeanX

private

Definition at line 96 of file CTPPSHector.h.

Referenced by add(), and CTPPSHector().

double CTPPSHector::fVtxMeanY

private

Definition at line 97 of file CTPPSHector.h.

Referenced by add(), and CTPPSHector().

double CTPPSHector::fVtxMeanZ

private

Definition at line 98 of file CTPPSHector.h.

Referenced by add(), and CTPPSHector().

double CTPPSHector::lengthctpps

private

Definition at line 80 of file CTPPSHector.h.

Referenced by CTPPSHector(), and filterCTPPS().

float CTPPSHector::m_b_ctpps_b

private

Definition at line 89 of file CTPPSHector.h.

Referenced by addPartToHepMC(), CTPPSHector(), and filterCTPPS().

std::unique_ptr<H_BeamLine> CTPPSHector::m_beamlineCTPPS1

private

Definition at line 104 of file CTPPSHector.h.

Referenced by CTPPSHector(), and filterCTPPS().

std::unique_ptr<H_BeamLine> CTPPSHector::m_beamlineCTPPS2

private

Definition at line 105 of file CTPPSHector.h.

Referenced by CTPPSHector(), and filterCTPPS().

std::map<unsigned int, H_BeamParticle*> CTPPSHector::m_beamPart

private

Definition at line 107 of file CTPPSHector.h.

Referenced by add(), addPartToHepMC(), clear(), filterCTPPS(), print(), and ~CTPPSHector().

bool CTPPSHector::m_CTPPSTransport

private

Definition at line 125 of file CTPPSHector.h.

Referenced by addPartToHepMC(), and CTPPSHector().

std::map<unsigned int, int> CTPPSHector::m_direct

private

Definition at line 108 of file CTPPSHector.h.

Referenced by add(), addPartToHepMC(), clear(), filterCTPPS(), and getDirect().

std::map<unsigned int, double> CTPPSHector::m_eAtTrPoint

private

Definition at line 114 of file CTPPSHector.h.

Referenced by addPartToHepMC(), and filterCTPPS().

std::map<unsigned int, double> CTPPSHector::m_eta

private

Definition at line 116 of file CTPPSHector.h.

Referenced by add(), and clear().

float CTPPSHector::m_f_ctpps_f

private

Definition at line 88 of file CTPPSHector.h.

Referenced by addPartToHepMC(), CTPPSHector(), and filterCTPPS().

std::map<unsigned int, bool> CTPPSHector::m_isCharged

private

Definition at line 119 of file CTPPSHector.h.

Referenced by add(), clear(), and filterCTPPS().

std::map<unsigned int, bool> CTPPSHector::m_isStoppedctpps

private

Definition at line 109 of file CTPPSHector.h.

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

std::map<unsigned int, int> CTPPSHector::m_pdg

private

Definition at line 117 of file CTPPSHector.h.

Referenced by add(), and clear().

std::map<unsigned int, double> CTPPSHector::m_pz

private

Definition at line 118 of file CTPPSHector.h.

Referenced by add(), and clear().

double CTPPSHector::m_sig_e

private

Definition at line 84 of file CTPPSHector.h.

Referenced by ApplyBeamCorrection(), and CTPPSHector().

double CTPPSHector::m_sigmaSTX

private

Definition at line 86 of file CTPPSHector.h.

Referenced by ApplyBeamCorrection(), and CTPPSHector().

double CTPPSHector::m_sigmaSTY

private

Definition at line 87 of file CTPPSHector.h.

Referenced by ApplyBeamCorrection(), and CTPPSHector().

bool CTPPSHector::m_smearAng

private

Definition at line 83 of file CTPPSHector.h.

Referenced by ApplyBeamCorrection(), and CTPPSHector().

bool CTPPSHector::m_smearE

private

Definition at line 85 of file CTPPSHector.h.

Referenced by ApplyBeamCorrection(), and CTPPSHector().

std::map<unsigned int, double> CTPPSHector::m_TxAtTrPoint

private

Definition at line 112 of file CTPPSHector.h.

Referenced by addPartToHepMC(), and filterCTPPS().

std::map<unsigned int, double> CTPPSHector::m_TyAtTrPoint

private

Definition at line 113 of file CTPPSHector.h.

Referenced by addPartToHepMC(), and filterCTPPS().

bool CTPPSHector::m_verbosity

private

Definition at line 124 of file CTPPSHector.h.

Referenced by add(), addPartToHepMC(), CTPPSHector(), and filterCTPPS().

std::map<unsigned int, double> CTPPSHector::m_xAtTrPoint

private

Definition at line 110 of file CTPPSHector.h.

Referenced by addPartToHepMC(), and filterCTPPS().

std::map<unsigned int, double> CTPPSHector::m_yAtTrPoint

private

Definition at line 111 of file CTPPSHector.h.

Referenced by addPartToHepMC(), and filterCTPPS().

int CTPPSHector::NEvent

private

Definition at line 128 of file CTPPSHector.h.

Referenced by addPartToHepMC().

edm::ESHandle< ParticleDataTable > CTPPSHector::pdt

private

Definition at line 101 of file CTPPSHector.h.

Referenced by add().

std::vector<LHCTransportLink> CTPPSHector::theCorrespondenceMap

private

Definition at line 127 of file CTPPSHector.h.

Referenced by addPartToHepMC(), CTPPSHector(), and getCorrespondenceMap().