...
How to get flavor tagging information inside Marlin processor:
Code Block |
---|
PIDHandler pidh( colJet ) ;// get collectionLCCollectionjet collection LCCollection* colJet = evt->getCollection("VertexJets"); // get PIDHandler associated with the jet collection PIDHandler pidh( colJet ); // get algorithm ID associated with LCFIPlus int algo = pidh.getAlgorithmID( "lcfiplus" ); // get index number for flavor tagging int ibtag = pidh.getParameterIndex(algo, "BTag"); int ictag = pidh.getParameterIndex(algo, "CTag"); // loop over jets to extract flavor tagging information for(int i=0; i < colJet->getNumberOfElements(); i++) { ReconstructedParticle *part = dynamic_cast<ReconstructedParticle*>( colJet->getElementAt( i ) ); const ParticleID &pid = pidh.getParticleID(part, algo); cout << "btag = " << pid.getParameters()[ibtag] << endl; cout << "ctag = " << pid.getParameters()[ictag] << endl; } |