subprepdatjob-data-e250-350Gflav-4jets.sh subproclst-data8-delay-e250-350Gflav-4jets
flavortag-all-batch-revtx-350-4jets.xml revertex-all-batch.xml
runitallbatch2 suballbatch2 allSME250p80m30.lst ffhzzE250p80m30.lst
[noric02] ~/sidhome/lcfi $ diff /u/ey/homer/sidhome/lcfi/build/v01-16-03/MarlinFastJet/v00-01/src/{FastJetProcessor.cpp,FastJetProcessor.cpp.original} 50d49 < // string jetparts = _lcJetOutName + string("Parts"); 52,53d50 < // registerOutputCollection(LCIO::RECONSTRUCTEDPARTICLE, "jetOut", "The identified jets", jetparts, "JetOut"); < registerOutputCollection(LCIO::RECONSTRUCTEDPARTICLE, "jetpartsOut", "The identified jet components", _lcJetPartsName, "JetParts"); 390,395d386 < // cout << "step1 - creating vector of particles" << endl; < // create output collection with particles used in it < IMPL::LCCollectionVec* lccJetsOutParts = new IMPL::LCCollectionVec(LCIO::RECONSTRUCTEDPARTICLE); < < int nelem = 0; < 404d394 < evt->addCollection(lccJetsOutParts, _lcJetPartsName); // store empty collection 411c401 < nelem = convertFromRecParticle(particleIn); --- > convertFromRecParticle(particleIn); 413,418d402 < if (nelem < 2) < { < _statsNrSkippedEmptyEvents++; < evt->addCollection(lccJetsOutParts, _lcJetPartsName); // store empty collection < throw DataNotAvailableException("Collection is there, but its empty after quality cuts!"); < } 441d424 < 443,444c426 < // if (_reconstructedPars->getNumberOfElements() < (int)_nrJets) < if (nelem < (int)_nrJets) --- > if (_reconstructedPars->getNumberOfElements() < (int)_nrJets) 491,546d472 < // cout << "step2 - beginning loop over jets" << endl; < for (it=jets.begin(); it != jets.end(); it++) < { < // create a reconstructed particle for this jet, and add all the containing particles to it < // ReconstructedParticle* recparts = getRecPar( (*it), cs.constituents(*it) ); < for (unsigned int n = 0; n < cs.constituents(*it).size(); ++n) < { < // cout << "step3 - adding particle" << endl; < ReconstructedParticle* p = dynamic_cast< ReconstructedParticle* > (_reconstructedPars->getElementAt(cs.constituents(*it)[n].user_index())) ; < < // create a true ReconstructedParticle < ReconstructedParticleImpl* recp = new ReconstructedParticleImpl(); < < // save the particle's parameters < recp->setEnergy( p->getEnergy() ); < recp->setMass( p->getMass() ); < recp->setCharge( p->getCharge() ); < recp->setMomentum( p->getMomentum() ); < recp->setType( p->getType() ); < recp->setCovMatrix( p->getCovMatrix() ); < recp->setMass( p->getMass() ); < recp->setCharge( p->getCharge() ); < recp->setReferencePoint( p->getReferencePoint() ); < for ( unsigned int k = 0 ; k < p->getParticleIDs().size() ; k++ ) < { < ParticleID *pid = p->getParticleIDs()[k] ; < IMPL::ParticleIDImpl* implPID = new IMPL::ParticleIDImpl() ; < implPID->setType( pid->getType() ) ; < implPID->setPDG( pid->getPDG() ) ; < implPID->setLikelihood( pid->getLikelihood() ) ; < implPID->setAlgorithmType( pid->getAlgorithmType() ) ; < for ( unsigned int l = 0 ; l < pid->getParameters().size() ; l++ ) < { implPID->addParameter( pid->getParameters()[l] ) ; } < recp->addParticleID( implPID ) ; < } < recp->setParticleIDUsed( p->getParticleIDUsed() ); < recp->setGoodnessOfPID( p->getGoodnessOfPID() ); < for ( unsigned int k = 0 ; k < p->getParticles().size() ; k++ ) < { recp->addParticle( p->getParticles()[k] ); } < for ( unsigned int k = 0 ; k < p->getClusters().size() ; k++ ) < { recp->addCluster( p->getClusters()[k] ); } < for ( unsigned int k = 0 ; k < p->getTracks().size() ; k++ ) < { recp->addTrack( p->getTracks()[k] ); } < recp->setStartVertex( p->getStartVertex() ); < < lccJetsOutParts->addElement( recp ); < < < } < } < < // cout << "step4 -- adding collection" << endl; < // string jetparts = _lcJetOutName + string("Parts"); < evt->addCollection(lccJetsOutParts, _lcJetPartsName); < // cout << "step5 --- past new code" << endl; < 560d485 < // cout << "step6" << endl; 698c623 < int FastJetProcessor::convertFromRecParticle(LCCollection* recCol) --- > void FastJetProcessor::convertFromRecParticle(LCCollection* recCol) 700d624 < int nelem=0; 705,720c629,634 < // add some quality cuts to resolve a background overlay problem - Homer < double px = par->getMomentum()[0]; < double py = par->getMomentum()[1]; < double pz = par->getMomentum()[2]; < < if ((fabs(pz)/sqrt(px*px+py*py+pz*pz)) < 0.99 && par->getEnergy() < 1000.) { < < nelem++; < < _pjList.push_back( < fastjet::PseudoJet( par->getMomentum()[0], < par->getMomentum()[1], < par->getMomentum()[2], < par->getEnergy() ) ); < _pjList.back().set_user_index(i); // save the id of this recParticle < } --- > _pjList.push_back( > fastjet::PseudoJet( par->getMomentum()[0], > par->getMomentum()[1], > par->getMomentum()[2], > par->getEnergy() ) ); > _pjList.back().set_user_index(i); // save the id of this recParticle 722d635
IsolatedHighPElectronIdentifier.java
variables_id_c1.png variables_id_c2.png variables_id_c3.png
factory->AddVariable( "evisjets", "Evis JETS", "GeV", 'F' );
factory->AddVariable( "mvisjets", "Mvis JETS", "GeV", 'F' );
factory->AddVariable( "ptvisjets", "PTvis JETS", "GeV/c", 'F' );
factory->AddVariable( "blike", "B Likeness", "", 'F' );
factory->AddVariable( "clike", "C likeness", "", 'F' );
factory->AddVariable( "nelec", "Number of HE Electrons", 'F' );
factory->AddVariable( "hmass", "Higgs Mass","GeV/c2", 'F' );
factory->AddVariable( "ejet1", "Z energy","GeV", 'F' );
factory->AddVariable( "ejet2", "Zstr energy", "GeV", 'F' );
factory->AddVariable( "cjet1", "Z cos theta","", 'F' );
factory->AddVariable( "cjet2", "Zstr cos theta", "", 'F' );
factory->AddVariable( "mjet1", "Z mass", "GeV/c2", 'F' );
factory->AddVariable( "mjet2", "Zstr mass", "GeV/c2", 'F' );
factory->AddVariable( "dphi", "azimuthal separation of Z jets", "", 'F' );
factory->AddVariable( "jetthrust", "thrust", "1", 'F' );
factory->AddVariable( "nchtrks", "Number Charged Tracks", "", 'F' );
factory->AddVariable( "tjpart1", "Number of electrons", "", 'F' );
factory->AddVariable( "tjpart2", "Number of muons", "", 'F' );
if ( EvisJETS<140.0){ if (nTrks>10 && EvisJETS>60.0 && PTvisJETS>5.0 && hmass>50. && hmass<140. && fabs(fabs(dph)-3.14159)>0.15 && ej1>40. && jetthrust<0.99) presel=1; // #4 LOW EVIS } else { if (nTrks>20 && EvisJETS<220.0 && PTvisJETS>5.0 && hmass>100. && hmass<150. && fabs(fabs(dph)-3.14159)>0.05 && ej1<130. && jetthrust<0.99) presel=1; // #10 }
// common0410
[neal@localhost weights]$ grep -A 8 -i optimal- ../hzzv14-presel-v4-new-common0410results.txt --- Classifier ( #signal, #backgr.) Optimal-cut S/sqrt(S+B) NSig NBkg EffSig EffBkg -- --------------------------------------------------------------------------------------------- --- Cuts: ( 664.0744,1312202.5) -0.0050 0 0 0 0 0 --- Likelihood: ( 664.0744,1312202.5) 1.0000 0.728991 644.0848 779980.1 0.9699 0.5944 --- Fisher: ( 664.0744,1312202.5) 0.0120 1.70335 258.2099 22721.13 0.3888 0.01732 --- BDTG: ( 664.0744,1312202.5) -0.9823 1.17754 538.5353 208621.3 0.811 0.159 --- BDT: ( 664.0744,1312202.5) -0.0462 3.42709 100.6751 762.2891 0.1516 0.0005809 -- ---------------------------------------------------------------------------------------------
isigall= 1345.984602 isigpresel= 1460.101057 igood= 114.116455 ibad= 762.500000 intot= 77692005.984714 ipre= 876.616455 cuts__ (Cut Name ): all others 2f 4f 6f aa 1f 3f 5f q e+mu tau nu Signif. cut #0 (all ): 7.769e+07 625000 14507062 1019938 0 54073784 0 7464875 0 948 91 45 262 0.15 Delta(sig*BR) = 6.548594 +/- 0.046879 cut #1 (20.<PTvisJETS ): 1.283e+07 37500 5525512 746188 0 2968250 0 3546875 0 806 69 44 259 0.33 Delta(sig*BR) = 3.042092 +/- 0.024172 cut #2 (EvisJETS<220. ): 4.637e+06 25000 3287550 514400 0 130196 0 678875 0 608 42 41 219 0.42 Delta(sig*BR) = 2.365748 +/- 0.020248 cut #3 (50.<hmass<140. ): 3.306e+06 0 2542562 378675 0 60426 0 323875 0 536 38 35 206 0.45 Delta(sig*BR) = 2.233245 +/- 0.019503 cut #4 (||dph|-3.14159|>0.15 ): 2.772e+06 0 2083912 364150 0 57097 0 266000 0 514 36 33 202 0.47 Delta(sig*BR) = 2.120439 +/- 0.018965 cut #5 (nTrks>10 ): 1.426e+06 0 1136375 169088 0 30283 0 89875 0 466 10 19 177 0.56 Delta(sig*BR) = 1.777889 +/- 0.017994 cut #6 (jetthrust<0.99 ): 1.325e+06 0 1044700 165588 0 30033 0 84125 0 463 10 19 173 0.58 Delta(sig*BR) = 1.730704 +/- 0.017728 cut #7 (ej1>40. ): 1.313e+06 0 1035112 164625 0 29215 0 83250 0 462 10 19 173 0.58 Delta(sig*BR) = 1.725415 +/- 0.017704 cut #8 (MVA ): 8.766e+02 0 138 625 0 0 0 0 0 51 0 1 62 3.85 Delta(sig*BR) = 0.259452 +/- 0.015456