extern void printGenParticle(const MCParticlePointer genParticle) {
printParticle(genParticle);
cout << "GenParticle Information" << endl;
cout << setw(30) << "pdgid" << setw(30) << "status" << setw(30) << "mother index" << setw(30) << "energy" << endl;
cout << setw(30) << genParticle->pdgId() << setw(30) << genParticle->status() << setw(30)
<< genParticle->motherIndex() << setw(30) << genParticle->energy() << endl;
}
extern void printMuon(const MuonPointer muon) {
printParticle(muon);
cout << "Muon Information" << endl;
cout << setw(30) << "is global" << setw(30) << "PF RelIso DR=0.4" << setw(30) << "----" << endl;
cout << setw(30) << muon->isGlobal() << setw(30) << muon->PFRelIso04() << setw(30) << "" << endl;
}
extern void printElectron(const ElectronPointer electron) {
printParticle(electron);
cout << "Electron Information" << endl;
cout << setw(30) << "isTight" << setw(30) << "isNonIso" << setw(30) << "isConv" << endl;
cout << setw(30) << electron->isTightElectron() << setw(30) << electron->isTightNonIsoElectron() << setw(30) << electron->isTightConversionElectron() << endl;
cout << setw(30) << "sigma_{ieta ieta}" << setw(30) << "|Delta phi_{in}|" << setw(30) << "|Delta eta_{in}|" << setw(30) << "HadOverEm" << endl;
cout << setw(30) << electron->sigmaIEtaIEta() << setw(30) << fabs(electron->dPhiIn()) << setw(30) << fabs(electron->dEtaIn()) << setw(30) << electron->HadOverEm() << endl;
cout << setw(30) << "Isolation" << setw(30) << "superClusterEta" << endl;
cout << setw(30) << electron->PFRelIso03DeltaBeta() << setw(30) << electron->superClusterEta() << endl;
}
static void display2() {
static GLfloat xRot = 0.0f, yRot = 0.0f, zRot = 0.0f;
glClear( GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT );
/* We don't want to modify the projection matrix. */
// glMatrixMode( GL_MODELVIEW );
glLoadIdentity( );
/* Move down the z-axis. */
glTranslatef( deltaX, deltaY, deltaZ );
if( rotateX ) {
//xRot += 45.0f;
if( ++xRot > 360.0f ) {
xRot = 0.0f;
}
//rotateX = !rotateX;
}
if( rotateY ) {
//yRot += 45.0f;
if( ++yRot > 360.0f ) {
yRot = 0.0f;
}
//rotateY = !rotateY;
}
if( rotateZ ) {
//zRot += 45.0f;
if( ++zRot > 360.0f ) {
zRot = 0.0f;
}
//rotateZ = !rotateZ;
}
/* Rotate. */
glRotatef( xRot, 1.0f, 0.0f, 0.0f );
glRotatef( yRot, 0.0f, 1.0f, 0.0f );
glRotatef( zRot, 0.0f, 0.0f, 1.0f ); /* Rotate. */
//printCube( cube0, index );
printParticle(particles, isovalue);
glutSwapBuffers();
glutPostRedisplay();
}
extern void printJet(const JetPointer jet) {
printParticle(jet);
cout << "Jet Information" << endl;
cout << setw(30) << "JECUnc" << endl;
cout << setw(30) << jet->JECUnc() << endl;
cout << setw(30) << "ptRaw" << setw(30) << "pxRaw" << setw(30) << "pyRaw" << setw(30) << "pzRaw" << endl;
cout << setw(30) << sqrt(jet->PxRaw() * jet->PxRaw() + jet->PyRaw() * jet->PyRaw()) << setw(30) << jet->PxRaw()
<< setw(30) << jet->PyRaw() << setw(30) << jet->PzRaw() << endl;
cout << setw(30) << "L1OffJEC" << setw(30) << "L2L3ResJEC" << setw(30) << "L2RelJEC" << setw(30) << "L3AbsJEC" << endl;
cout << setw(30) << jet->L1OffJEC() << setw(30) << jet->L2L3ResJEC() << setw(30) << jet->L2RelJEC() << setw(30) << jet->L3AbsJEC() << endl;;
cout << setw(30) << "emf" << setw(30) << "n90Hits" << setw(30) << "fHPD" << setw(30) << "B tag(SSV)" << endl;
cout << setw(30) << jet->emf() << setw(30) << jet->n90Hits() << setw(30) << jet->fHPD() << setw(30)
<< jet->isBJet(BtagAlgorithm::CombinedSecondaryVertexV2, BtagAlgorithm::MEDIUM) << endl << endl;
if (jet->getUsedAlgorithm() == JetAlgorithm::CA08PF || jet->getUsedAlgorithm() == JetAlgorithm::PF2PAT) {
printPFJetExtra(jet);
}
}
