IParticleHandle_Jet::Imp Class Reference

Collaboration diagram for IParticleHandle_Jet::Imp:

Public Member Functions
double	pt () const
double	phi () const
double	eta () const
double	energy () const
double	transverseEnergy () const
double	coneR () const
double	energyForLengthAndCuts () const
void	createShapeFromJetParameters (const IParticleCollHandle_Jet *collHandleJet, const double &coneR, const double &eta, const double &phi, const double &energy, const SbVec3f &origin)
void	updateConeHeightParameters (SoCone *, SoTranslation *, const double &energy) const
void	updateConeHeightParameters () const
void	rerandomiseMaterial ()

Public Attributes
IParticleHandle_Jet *	theclass = nullptr
const IParticleCollHandle_Jet *	theCollHandle = nullptr
const xAOD::Jet *	m_jet = nullptr
SoSeparator *	sep = nullptr
SoCone *	cone = nullptr
SoMaterial *	m_randomMat = nullptr
SoSwitch *	m_bTagged = nullptr
SoSwitch *	m_bTaggingCollSwitch = nullptr
double	m_bTagWeightMV1 = 0.0
double	m_bTagWeightMV2c20 = 0.0
double	m_bTagWeightMV2c10 = 0.0
double	m_JetFitterCombNN_pb = 0.0
double	m_JetFitterCombNN_pc = 0.0
double	m_JetFitterCombNN_pu = 0.0
bool	considerTransverseEnergies = true
double	coneRPar = -1
double	scale = (10.0 * SYSTEM_OF_UNITS::m) / (100.0*SYSTEM_OF_UNITS::GeV)
double	maxR = 0.0 * SYSTEM_OF_UNITS::m

Detailed Description

Definition at line 51 of file IParticleHandle_Jet.cxx.

Member Function Documentation

coneR()

double IParticleHandle_Jet::Imp::coneR ( ) const

inline

Definition at line 97 of file IParticleHandle_Jet.cxx.

{ return coneRPar > 0 ? coneRPar : 0.4; }

createShapeFromJetParameters()

void IParticleHandle_Jet::Imp::createShapeFromJetParameters	(	const IParticleCollHandle_Jet *	collHandleJet,
		const double &	coneR,
		const double &	eta,
		const double &	phi,
		const double &	energy,
		const SbVec3f &	origin
	)

Definition at line 245 of file IParticleHandle_Jet.cxx.

{
    VP1Msg::messageVerbose("IParticleHandle_Jet::Imp::createShapeFromJetParameters()");
 
 
 
    sep = new SoSeparator();
    sep->ref();
 
    cone = new SoCone();
    cone->ref();
 
    // coneR IS the opening half-angle of the jet, in delta phi (easy) and in
    //  delta eta (trickier)
    // to try to get the jet extent right in eta, find theta corresponding
    //  to eta+R, eta-R, and take half of the difference:
    double thetaMax = 2.*atan(exp(-(eta+inputconeR)));
    double thetaMin = 2.*atan(exp(-(eta-inputconeR)));
    double deltaTheta = fabs(thetaMax-thetaMin);
    double etaScale = deltaTheta/(2.*inputconeR);
 
    // Translate by half cone height and flip 180 deg so point is at IP:
    SoTranslation *translate = new SoTranslation();
 
    SoRotationXYZ *flip = new SoRotationXYZ();
    flip->axis=SoRotationXYZ::Z;
    flip->angle=M_PI;
 
    SoRotationXYZ *ytoz = new SoRotationXYZ();
    ytoz->axis=SoRotationXYZ::X;
    ytoz->angle=M_PI/2.;
 
    // Cones should now be along Z-axis,point at IP.
    SoRotationXYZ *rotationPhi = new SoRotationXYZ();
    rotationPhi->axis=SoRotationXYZ::Z;
    rotationPhi->angle = phi+M_PI/2.; // starts from -y-axis in x-y plane
    // phi is measured from x-axis, so +M_PI/2
 
    SoRotationXYZ *rotationEta = new SoRotationXYZ();
    double signEta = fabs(eta)/eta;
    double theta = 2.*atan(signEta*exp(-fabs(eta)));
    if (theta<0.)theta+=M_PI;
    rotationEta->axis=SoRotationXYZ::X;
    rotationEta->angle=theta;
 
    //message("Eta: "+QString::number(eta)+" theta: "+QString::number(theta)+" phi: "+QString::number(phi));
 
    updateConeHeightParameters(cone, translate, energy);
 
    // play with Scale:
    SoScale* myScale = new SoScale();
    //message("Eta scaling factor is "+QString::number(etaScale));
    // maybe we need to squish along both axes...
    myScale->scaleFactor.setValue(etaScale,1.,etaScale);
 
    // translate to origin vertex rather than (0,0,0)
    SoTranslation *transvertex = new SoTranslation();
    transvertex->translation = origin;
 
    sep->addChild(transvertex);
    sep->addChild(rotationPhi);
    sep->addChild(rotationEta);//theta rotation done around x-axis
    sep->addChild(ytoz); // now it's along the z-axis
    sep->addChild(flip); // flip so tip at origin
    sep->addChild(translate); // back it up so base at origin // DO NOT MOVE THIS: ITS POSITION IS USED BY "updateConeHeightParameters(SoSeparator* sep,const double& energy)"
    sep->addChild(myScale); // squeeze jet according to eta
    //This is the point in the child sequence where we MAY add a random colour when appropriate.
    //Thus: The translation is the SIXTH child and the cone is the LAST child.
 
    //  sep->addChild(collHandleJet->collSettingsButton().defaultParameterMaterial());
 
    sep->addChild(collHandleJet->material());
    // ^^ FIXME - should rearrange so we don't need to reset material
 
    sep->addChild(cone);  // starts along y-axis // DO NOT MOVE THIS: ITS POSITION IS USED BY "updateConeHeightParameters(SoSeparator* sep,const double& energy)"
 
 
    // rendering b-tagged jets
    SoNode * lastChild = sep->getChild( sep->getNumChildren()-1 );
    if (lastChild->getTypeId() == SoCone::getClassTypeId())
    {
        VP1Msg::messageVerbose("adding b-tagging switch");
 
        m_bTagged = new SoSwitch();
        m_bTagged->whichChild=SO_SWITCH_ALL;
 
        m_bTagged->addChild(m_bTaggingCollSwitch);
 
        // add the b-tagging switch before the SoCone node
        sep->insertChild(m_bTagged, sep->getNumChildren()-1);
    }
 
 
    //return sep;
}

energy()

double IParticleHandle_Jet::Imp::energy ( ) const

inline

Definition at line 95 of file IParticleHandle_Jet.cxx.

{ return m_jet->e(); }

energyForLengthAndCuts()

double IParticleHandle_Jet::Imp::energyForLengthAndCuts ( ) const

inline

Definition at line 98 of file IParticleHandle_Jet.cxx.

{ return considerTransverseEnergies ? transverseEnergy() : energy(); }

eta()

double IParticleHandle_Jet::Imp::eta ( ) const

inline

Definition at line 94 of file IParticleHandle_Jet.cxx.

{ /*VP1Msg::messageVerbose("eta: " + QString::number(m_jet->eta()) );*/ return m_jet->eta(); }

phi()

double IParticleHandle_Jet::Imp::phi ( ) const

inline

Definition at line 93 of file IParticleHandle_Jet.cxx.

{ /*VP1Msg::messageVerbose("phi: " + QString::number(m_jet->phi()) );*/ return m_jet->phi(); }

pt()

double IParticleHandle_Jet::Imp::pt ( ) const

inline

Definition at line 92 of file IParticleHandle_Jet.cxx.

{ return m_jet->pt(); }

rerandomiseMaterial()

void IParticleHandle_Jet::Imp::rerandomiseMaterial

(

)

Definition at line 592 of file IParticleHandle_Jet.cxx.

{
    VP1Msg::messageVerbose("IParticleHandle_Jet::Imp::rerandomiseMaterial()");
 
    // TODO: Move to HSV system, instead of RGB, and make random colors more contrasting
 
 
    //Fixme: share this code with other systems!!
    if ( !m_randomMat ) { //We will anyway rerandomize it when we need it
        VP1Msg::messageVerbose("'m_randomMat not set. Returning.");
        return;
    }
 
    double r2 = 0.3*0.3;
    unsigned i(0);
    double r,g,b;
    bool ok;
    while (true) {
        r = (rand() / static_cast<double>(RAND_MAX));
        g = (rand() / static_cast<double>(RAND_MAX));
        b = (rand() / static_cast<double>(RAND_MAX));
        ok = true;
        //For now we make sure that we avoid black and red. This should be updated from bgd and highlight col automatically! (fixme).
        // -> and we should probably also make sure that tracks close in (eta,phi) are well separated in colour-space.
        if ( (r-1.0)*(r-1.0)+g*g+b*b < r2*0.5 )//avoid red (distance)
            ok = false;
        else if ( r*r/(r*r+g*g+b*b) > 0.8 )//avoid red (angle)
            ok = false;
        else if ( r*r+g*g+b*b < r2*2.0 )//avoid black
            ok = false;
        if (ok)
            break;
        ++i;
        if (i>50 ) {
            r2 *= 0.99;//To avoid problem in case we add too many forbidden spheres.
            if (i>1000) {
                //Just a safety
                VP1Msg::messageVerbose("IParticleHandle_Jet::Imp::rerandomiseMaterial() - Warning: Random colour could"
                        " not be selected such as to satisfy all separation criterias");
                break;
            }
        }
    }
 
    VP1MaterialButton::setMaterialParameters(m_randomMat,r, g, b, 0.1, 0.5 );
}

transverseEnergy()

double IParticleHandle_Jet::Imp::transverseEnergy ( ) const

inline

Definition at line 96 of file IParticleHandle_Jet.cxx.

{ return sin(2*atan(exp(-fabs(eta()))))*energy(); }

updateConeHeightParameters() [1/2]

void IParticleHandle_Jet::Imp::updateConeHeightParameters

(

)

const

Definition at line 376 of file IParticleHandle_Jet.cxx.

                                                              {
 
 
    if (!sep) {
        VP1Msg::messageVerbose("sep not defined. Returning.");
        return;
    }
 
    //NB: The translation is the SIXTH child and the cone is the LAST child.
    if (sep->getNumChildren()<6) {
        VP1Msg::messageVerbose("getNumChildren() < 6!!! Returning...");
        return;
    }
 
    const double energyJet = energyForLengthAndCuts();
 
    SoNode * sixthChild = sep->getChild(5);
    if (sixthChild->getTypeId()!=SoTranslation::getClassTypeId()) {
        VP1Msg::messageVerbose("6th child is not a translation!!! Returning...");
        return;
    }
    SoNode * lastChild = sep->getChild(sep->getNumChildren()-1);
    if (lastChild->getTypeId()!=SoCone::getClassTypeId()) {
        VP1Msg::messageVerbose("lastChild is not a cone!!! Returning...");
        return;
    }
    updateConeHeightParameters(static_cast<SoCone*>(lastChild),static_cast<SoTranslation*>(sixthChild), energyJet);
}

updateConeHeightParameters() [2/2]

void IParticleHandle_Jet::Imp::updateConeHeightParameters	(	SoCone *	cone,
		SoTranslation *	trans,
		const double &	energy
	)		const

Definition at line 343 of file IParticleHandle_Jet.cxx.

{
    VP1Msg::messageVerbose("IParticleHandle_Jet::Imp::updateConeHeightParameters()");
 
    double h( std::fabs(scale * energy) );
 
    // FIXME: check this! when maxR is set, what should that do?? Here it's only used as "h", but maybe it should set the max length of the cone's bottom radius?? Check with Ed!
    if (maxR > 0.0001) { // maxR > 0. is not univocally defined, because maxR is a floating point number
        VP1Msg::messageVerbose("maxR > 0.0001 ==> setting h = min(h, maxR)");
        h = std::min(h, maxR);
    }
 
    cone->height = h;
    cone->bottomRadius = tan(coneR()) * h;
    trans->translation = SbVec3f(0,-0.5*h, 0);
 
 
 
    //debug
    SbString strHeight, strRadius;
    (cone->height).get(strHeight);
    (cone->bottomRadius).get(strRadius);
//   std::cout << "input - energy: " << energy << " - scale: " << scale << " - maxR: " << maxR << " - h: " << h << " --- updated cone - height: " << strHeight.getString() << " - bottom radius: " << strRadius.getString() << std::endl;
 
    // you can also use the 'writeField()' method, direct to std output
    //SoDebug::writeField(&(cone->height));
 
    //  QString text = "updated cone - height: " + QString::number(cone->height) + " - bottom radius: " + QString::number(cone->bottomRadius) " - translation: " + trans->translation;
    //  VP1Msg::messageVerbose(text);
}

Member Data Documentation

cone

SoCone* IParticleHandle_Jet::Imp::cone = nullptr

Definition at line 63 of file IParticleHandle_Jet.cxx.

coneRPar

double IParticleHandle_Jet::Imp::coneRPar = -1

Definition at line 82 of file IParticleHandle_Jet.cxx.

considerTransverseEnergies

bool IParticleHandle_Jet::Imp::considerTransverseEnergies = true

Definition at line 81 of file IParticleHandle_Jet.cxx.

m_bTagged

SoSwitch* IParticleHandle_Jet::Imp::m_bTagged = nullptr

Definition at line 65 of file IParticleHandle_Jet.cxx.

m_bTaggingCollSwitch

SoSwitch* IParticleHandle_Jet::Imp::m_bTaggingCollSwitch = nullptr

Definition at line 66 of file IParticleHandle_Jet.cxx.

m_bTagWeightMV1

double IParticleHandle_Jet::Imp::m_bTagWeightMV1 = 0.0

Definition at line 70 of file IParticleHandle_Jet.cxx.

m_bTagWeightMV2c10

double IParticleHandle_Jet::Imp::m_bTagWeightMV2c10 = 0.0

Definition at line 72 of file IParticleHandle_Jet.cxx.

m_bTagWeightMV2c20

double IParticleHandle_Jet::Imp::m_bTagWeightMV2c20 = 0.0

Definition at line 71 of file IParticleHandle_Jet.cxx.

m_jet

const xAOD::Jet* IParticleHandle_Jet::Imp::m_jet = nullptr

Definition at line 60 of file IParticleHandle_Jet.cxx.

m_JetFitterCombNN_pb

double IParticleHandle_Jet::Imp::m_JetFitterCombNN_pb = 0.0

Definition at line 73 of file IParticleHandle_Jet.cxx.

m_JetFitterCombNN_pc

double IParticleHandle_Jet::Imp::m_JetFitterCombNN_pc = 0.0

Definition at line 74 of file IParticleHandle_Jet.cxx.

m_JetFitterCombNN_pu

double IParticleHandle_Jet::Imp::m_JetFitterCombNN_pu = 0.0

Definition at line 75 of file IParticleHandle_Jet.cxx.

m_randomMat

SoMaterial* IParticleHandle_Jet::Imp::m_randomMat = nullptr

Definition at line 64 of file IParticleHandle_Jet.cxx.

maxR

double IParticleHandle_Jet::Imp::maxR = 0.0 * SYSTEM_OF_UNITS::m

Definition at line 84 of file IParticleHandle_Jet.cxx.

scale

double IParticleHandle_Jet::Imp::scale = (10.0 * SYSTEM_OF_UNITS::m) / (100.0*SYSTEM_OF_UNITS::GeV)

Definition at line 83 of file IParticleHandle_Jet.cxx.

sep

SoSeparator* IParticleHandle_Jet::Imp::sep = nullptr

Definition at line 62 of file IParticleHandle_Jet.cxx.

theclass

IParticleHandle_Jet* IParticleHandle_Jet::Imp::theclass = nullptr

Definition at line 57 of file IParticleHandle_Jet.cxx.

theCollHandle

const IParticleCollHandle_Jet* IParticleHandle_Jet::Imp::theCollHandle = nullptr

Definition at line 58 of file IParticleHandle_Jet.cxx.

---

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

• IParticleHandle_Jet.cxx