#include <PFRootEventManagerColin.h>

Inheritance diagram for PFRootEventManagerColin:

Public Types
enum	Mode { Neutral =0, HIGH_E_TAUS =1 }

Public Types inherited from PFRootEventManager
typedef std::map< int, int >	EventToEntry

typedef std::map< int, EventToEntry >	LumisMap

typedef std::map< int, LumisMap >	RunsMap

enum	Verbosity { SHUTUP = 0, VERBOSE }

enum	View_t { XY = 0, RZ = 1, EPE = 2, EPH = 3, NViews = 4 }
	viewport definition More...

Public Member Functions
	PFRootEventManagerColin (const char *file)

bool	processEntry (int entry)
	process one entry (pass the TTree entry) More...

bool	processHIGH_E_TAUS ()

bool	processNeutral ()

void	readSpecificOptions (const char *file)

void	write ()

	~PFRootEventManagerColin ()

Public Member Functions inherited from PFRootEventManager
const reco::PFBlockCollection &	blocks () const

int	chargeValue (const int &pdgId) const

const reco::PFSimParticle &	closestParticle (reco::PFTrajectoryPoint::LayerType layer, double eta, double phi, double &peta, double &pphi, double &pe) const
	find the closest PFSimParticle to a point (eta,phi) in a given detector More...

void	clustering ()
	read data from testbeam tree More...

void	connect (const char *infilename="")
	open the root file and connect to the tree More...

bool	countChargedAndPhotons () const

bool	eventAccepted () const
	returns true if the event is accepted(have a look at the function implementation) More...

int	eventNumber ()

int	eventToEntry (int run, int lumi, int event) const

std::string	expand (const std::string &oldString) const

void	fillClusterMask (std::vector< bool > &mask, const reco::PFClusterCollection &clusters) const
	cluster mask set to true for rechits inside TCutG More...

void	fillOutEventWithBlocks (const reco::PFBlockCollection &blocks)
	fills outEvent with blocks More...

void	fillOutEventWithCaloTowers (const CaloTowerCollection &cts)
	fills outEvent with calo towers More...

void	fillOutEventWithClusters (const reco::PFClusterCollection &clusters)
	fills OutEvent with clusters More...

void	fillOutEventWithPFCandidates (const reco::PFCandidateCollection &pfCandidates)
	fills OutEvent with candidates More...

void	fillOutEventWithSimParticles (const reco::PFSimParticleCollection &ptcs)
	fills OutEvent with sim particles More...

void	fillRecHitMask (std::vector< bool > &mask, const reco::PFRecHitCollection &rechits) const
	rechit mask set to true for rechits inside TCutG More...

void	fillTrackMask (std::vector< bool > &mask, const reco::PFRecTrackCollection &tracks) const
	track mask set to true for rechits inside TCutG More...

void	fillTrackMask (std::vector< bool > &mask, const reco::GsfPFRecTrackCollection &tracks) const

std::string	getGenParticleName (int partId, std::string &latexStringName) const
	get name of genParticle More...

bool	highPtJet (double ptMin) const
	returns true if there is at least one jet with pT>pTmin More...

bool	highPtPFCandidate (double ptMin, reco::PFCandidate::ParticleType type=reco::PFCandidate::X) const
	returns true if there is a PFCandidate of a given type over a given pT More...

void	initializeEventInformation ()

bool	isHadronicTau () const
	study the sim event to check if the tau decay is hadronic More...

void	mcTruthMatching (std::ostream &out, const reco::PFCandidateCollection &candidates, std::vector< std::list< simMatch > > &candSimMatchTrack, std::vector< std::list< simMatch > > &candSimMatchEcal) const

void	particleFlow ()
	performs particle flow More...

void	pfCandCompare (int)
	compare particle flow More...

	PFRootEventManager ()
	default constructor More...

	PFRootEventManager (const char *file)

void	PreprocessRecHits (reco::PFRecHitCollection &rechits, bool findNeighbours)
	preprocess a rechit vector from a given rechit branch More...

void	PreprocessRecTracks (reco::PFRecTrackCollection &rectracks)
	preprocess a rectrack vector from a given rectrack branch More...

void	PreprocessRecTracks (reco::GsfPFRecTrackCollection &rectracks)

void	print (std::ostream &out=std::cout, int maxNLines=-1) const
	print information More...

void	printCluster (const reco::PFCluster &cluster, std::ostream &out=std::cout) const

void	printClusters (const reco::PFClusterCollection &clusters, std::ostream &out=std::cout) const
	print clusters More...

void	printGenParticles (std::ostream &out=std::cout, int maxNLines=-1) const
	print the HepMC truth More...

void	printMCCalib (std::ofstream &out) const
	print calibration information More...

void	printRecHit (const reco::PFRecHit &rh, unsigned index, const char *seed=" ", std::ostream &out=std::cout) const

void	printRecHits (const reco::PFRecHitCollection &rechits, const PFClusterAlgo &clusterAlgo, std::ostream &out=std::cout) const
	print rechits More...

virtual bool	processEvent (int run, int lumi, int event)
	process one event (pass the CMS event number) More...

void	readCMSSWJets ()

bool	readFromSimulation (int entry)
	read data from simulation tree More...

void	readOptions (const char *file, bool refresh=true, bool reconnect=false)

void	reconstructCaloJets ()
	reconstruct calo jets More...

void	reconstructFWLiteJets (const reco::CandidatePtrVector &Candidates, std::vector< ProtoJet > &output)
	used by the reconstruct*Jets functions More...

void	reconstructGenJets ()
	reconstruct gen jets More...

void	reconstructPFJets ()
	reconstruct pf jets More...

void	reset ()
	reset before next event More...

void	setRecHitNeigbours (reco::PFRecHit &rh, const std::map< unsigned, unsigned > &detId2index)

edm::InputTag	stringToTag (const std::vector< std::string > &tagname)
	returns an InputTag from a vector of strings More...

double	tauBenchmark (const reco::PFCandidateCollection &candidates)
	COLIN need to get rid of this mess. More...

bool	trackInsideGCut (const reco::PFTrack &track) const
	is PFTrack inside cut G ? yes if at least one trajectory point is inside. More...

TTree *	tree ()
	get tree More...

virtual	~PFRootEventManager ()
	destructor More...

Public Attributes
int	mode_

NeutralEvent *	neutralEvent_

TTree *	outTreeMy_

TauEvent *	tauEvent_

Public Attributes inherited from PFRootEventManager
std::ofstream *	calibFile_

boost::shared_ptr < PFEnergyCalibration >	calibration_

std::vector< ProtoJet >	caloJets_
	calo Jets More...

std::vector< reco::CaloJet >	caloJetsCMSSW_

edm::Handle< std::vector < reco::CaloJet > >	caloJetsHandle_
	CMSSW calo Jets. More...

edm::InputTag	caloJetsTag_

reco::CaloMETCollection	caloMets_
	Calo MET. More...

reco::CaloMETCollection	caloMetsCMSSW_

edm::Handle < reco::CaloMETCollection >	caloMetsHandle_
	CMSSW Calo MET. More...

edm::InputTag	caloMetsTag_

CaloTowerCollection	caloTowers_

edm::Handle< CaloTowerCollection >	caloTowersHandle_
	input collection of calotowers More...

reco::CandidatePtrVector	caloTowersPtrs_

edm::InputTag	caloTowersTag_

PFClusterAlgo	clusterAlgoECAL_

PFClusterAlgo	clusterAlgoHCAL_
	clustering algorithm for HCAL More...

PFClusterAlgo	clusterAlgoHFEM_
	clustering algorithm for HF, electro-magnetic layer More...

PFClusterAlgo	clusterAlgoHFHAD_
	clustering algorithm for HF, hadronic layer More...

PFClusterAlgo	clusterAlgoPS_
	clustering algorithm for PS More...

boost::shared_ptr < pftools::PFClusterCalibration >	clusterCalibration_

std::auto_ptr < reco::PFClusterCollection >	clustersECAL_

edm::Handle < reco::PFClusterCollection >	clustersECALHandle_
	clusters ECAL More...

edm::InputTag	clustersECALTag_

std::auto_ptr < reco::PFClusterCollection >	clustersHCAL_

edm::Handle < reco::PFClusterCollection >	clustersHCALHandle_
	clusters HCAL More...

edm::InputTag	clustersHCALTag_

std::auto_ptr < reco::PFClusterCollection >	clustersHFEM_

edm::Handle < reco::PFClusterCollection >	clustersHFEMHandle_
	clusters HCAL More...

edm::InputTag	clustersHFEMTag_

std::auto_ptr < reco::PFClusterCollection >	clustersHFHAD_

edm::Handle < reco::PFClusterCollection >	clustersHFHADHandle_
	clusters HCAL More...

edm::InputTag	clustersHFHADTag_

std::auto_ptr < reco::PFClusterCollection >	clustersPS_

edm::Handle < reco::PFClusterCollection >	clustersPSHandle_
	clusters PS More...

edm::InputTag	clustersPSTag_

reco::GsfPFRecTrackCollection	convBremGsfrecTracks_

edm::Handle < reco::GsfPFRecTrackCollection >	convBremGsfrecTracksHandle_
	reconstructed secondary GSF tracks More...

edm::InputTag	convBremGsfrecTracksTag_

reco::PFConversionCollection	conversion_

edm::Handle < reco::PFConversionCollection >	conversionHandle_
	conversions More...

edm::InputTag	conversionTag_

std::vector< reco::CaloJet >	corrcaloJetsCMSSW_

edm::Handle< std::vector < reco::CaloJet > >	corrcaloJetsHandle_
	CMSSW corrected calo Jets. More...

edm::InputTag	corrcaloJetsTag_

bool	debug_
	debug printouts for this PFRootEventManager on/off More...

double	DeltaMETcut

double	DeltaPhicut

reco::PFRecTrackCollection	displacedRecTracks_

edm::Handle < reco::PFRecTrackCollection >	displacedRecTracksHandle_

edm::InputTag	displacedRecTracksTag_

bool	doClustering_

bool	doCompare_
	comparison with pf CMSSW More...

bool	doJets_
	jets on/off More...

bool	doMet_
	MET on/off. More...

bool	doParticleFlow_
	particle flow on/off More...

bool	doPFCandidateBenchmark_

bool	doPFJetBenchmark_
	PFJet benchmark on/off. More...

bool	doPFMETBenchmark_
	PFMET benchmark on/off. More...

bool	doTauBenchmark_
	tau benchmark on/off More...

fwlite::ChainEvent *	ev_
	NEW: input event. More...

edm::EventAuxiliary *	eventAuxiliary_
	event auxiliary information More...

TBranch *	eventAuxiliaryBranch_

bool	fastsim_
	Fastsim or fullsim. More...

TFile *	file_
	input file More...

bool	filterHadronicTaus_

unsigned	filterNParticles_

std::vector< int >	filterTaus_

bool	findRecHitNeighbours_
	find rechit neighbours ? More...

reco::GenJetCollection	genJets_
	gen Jets More...

reco::GenJetCollection	genJetsCMSSW_

edm::Handle < reco::GenJetCollection >	genJetsHandle_
	CMSSW gen Jets. More...

edm::InputTag	genJetsTag_

reco::GenParticleCollection	genParticlesCMSSW_

reco::GenParticleRefVector	genParticlesforJets_

edm::Handle < reco::GenParticleRefVector >	genParticlesforJetsHandle_
	input collection of gen particles More...

reco::CandidatePtrVector	genParticlesforJetsPtrs_

edm::InputTag	genParticlesforJetsTag_

edm::Handle < reco::GenParticleCollection >	genParticlesforMETHandle_
	CMSSW GenParticles. More...

edm::InputTag	genParticlesforMETTag_

reco::GsfPFRecTrackCollection	gsfrecTracks_

edm::Handle < reco::GsfPFRecTrackCollection >	gsfrecTracksHandle_
	reconstructed GSF tracks More...

edm::InputTag	gsfrecTracksTag_

TH1F *	h_deltaETvisible_MCEHT_
	output histo dET ( EHT - MC) More...

TH1F *	h_deltaETvisible_MCPF_
	output histo dET ( PF - MC) More...

int	iEvent_
	current event More...

std::vector< std::string >	inFileNames_
	input file names More...

bool	JECinCaloMet_
	propagate the Jet Energy Corrections to the caloMET on/off More...

PFJetAlgorithm	jetAlgo_

int	jetAlgoType_
	jet algo type More...

FWLiteJetProducer	jetMaker_
	wrapper to official jet algorithms More...

bool	jetsDebug_
	debug printouts for jet algo on/off More...

RunsMap	mapEventToEntry_

edm::HepMCProduct	MCTruth_

edm::Handle< edm::HepMCProduct >	MCTruthHandle_
	MC truth. More...

edm::InputTag	MCTruthTag_

double	MET1cut
	PFMET Benchmark. More...

std::auto_ptr< METManager >	metManager_

reco::MuonCollection	muons_

edm::Handle< reco::MuonCollection >	muonsHandle_
	muons More...

edm::InputTag	muonsTag_

IO *	options_
	options file parser More...

EventColin *	outEvent_

TFile *	outFile_
	output file More...

std::string	outFileName_
	output filename More...

TTree *	outTree_
	output tree More...

PFAlgo	pfAlgo_
	particle flow algorithm More...

PFBlockAlgo	pfBlockAlgo_
	algorithm for building the particle flow blocks More...

std::auto_ptr < reco::PFBlockCollection >	pfBlocks_
	reconstructed pfblocks More...

reco::PFCandidateCollection	pfCandCMSSW_

edm::Handle < reco::PFCandidateCollection >	pfCandidateHandle_
	CMSSW PF candidates. More...

PFCandidateManager	pfCandidateManager_

std::auto_ptr < reco::PFCandidateCollection >	pfCandidates_
	reconstructed pfCandidates More...

reco::CandidatePtrVector	pfCandidatesPtrs_

edm::InputTag	pfCandidateTag_

PFJetBenchmark	PFJetBenchmark_
	PFJet Benchmark. More...

reco::PFJetCollection	pfJets_
	PF Jets. More...

reco::PFJetCollection	pfJetsCMSSW_

edm::Handle < reco::PFJetCollection >	pfJetsHandle_
	CMSSW PF Jets. More...

edm::InputTag	pfJetsTag_

reco::PFMETCollection	pfMets_
	PF MET. More...

reco::PFMETCollection	pfMetsCMSSW_

edm::Handle < reco::PFMETCollection >	pfMetsHandle_
	CMSSW PF MET. More...

edm::InputTag	pfMetsTag_

reco::PFDisplacedTrackerVertexCollection	pfNuclearTrackerVertex_

edm::Handle < reco::PFDisplacedTrackerVertexCollection >	pfNuclearTrackerVertexHandle_
	PFDisplacedVertex. More...

edm::InputTag	pfNuclearTrackerVertexTag_

reco::VertexCollection	primaryVertices_

edm::Handle < reco::VertexCollection >	primaryVerticesHandle_
	reconstructed primary vertices More...

edm::InputTag	primaryVerticesTag_

bool	printClusters_
	print clusters yes/no More...

double	printClustersEMin_

bool	printGenParticles_
	print MC truth yes/no More...

double	printGenParticlesPtMin_

bool	printMCTruthMatching_

bool	printPFBlocks_
	print PFBlocks yes/no More...

bool	printPFCandidates_
	print PFCandidates yes/no More...

double	printPFCandidatesPtMin_

bool	printPFJets_
	print PFJets yes/no More...

double	printPFJetsPtMin_

bool	printRecHits_
	print rechits yes/no More...

double	printRecHitsEMin_

bool	printSimParticles_
	print true particles yes/no More...

double	printSimParticlesPtMin_

reco::PFRecHitCollection	rechitsCLEANED_

std::vector< edm::Handle < reco::PFRecHitCollection > >	rechitsCLEANEDHandles_

std::vector< edm::InputTag >	rechitsCLEANEDTags_

std::vector < reco::PFRecHitCollection >	rechitsCLEANEDV_
	rechits HF CLEANED More...

reco::PFRecHitCollection	rechitsECAL_

edm::Handle < reco::PFRecHitCollection >	rechitsECALHandle_
	rechits ECAL More...

edm::InputTag	rechitsECALTag_

reco::PFRecHitCollection	rechitsHCAL_

edm::Handle < reco::PFRecHitCollection >	rechitsHCALHandle_
	rechits HCAL More...

edm::InputTag	rechitsHCALTag_

reco::PFRecHitCollection	rechitsHFEM_

edm::Handle < reco::PFRecHitCollection >	rechitsHFEMHandle_
	rechits HF EM More...

edm::InputTag	rechitsHFEMTag_

reco::PFRecHitCollection	rechitsHFHAD_

edm::Handle < reco::PFRecHitCollection >	rechitsHFHADHandle_
	rechits HF HAD More...

edm::InputTag	rechitsHFHADTag_

reco::PFRecHitCollection	rechitsPS_

edm::Handle < reco::PFRecHitCollection >	rechitsPSHandle_
	rechits PS More...

edm::InputTag	rechitsPSTag_

reco::PFRecTrackCollection	recTracks_

edm::Handle < reco::PFRecTrackCollection >	recTracksHandle_
	reconstructed tracks More...

edm::InputTag	recTracksTag_

reco::TrackCollection	stdTracks_

edm::Handle < reco::TrackCollection >	stdTracksHandle_
	standard reconstructed tracks More...

edm::InputTag	stdTracksTag_

bool	tauBenchmarkDebug_
	tau benchmark debug More...

reco::METCollection	tcMets_
	TCMET. More...

reco::METCollection	tcMetsCMSSW_

edm::Handle< reco::METCollection >	tcMetsHandle_
	CMSSW TCMET. More...

edm::InputTag	tcMetsTag_

boost::shared_ptr < PFEnergyCalibrationHF >	thepfEnergyCalibrationHF_

TTree *	tree_
	input tree More...

reco::PFSimParticleCollection	trueParticles_

edm::Handle < reco::PFSimParticleCollection >	trueParticlesHandle_
	true particles More...

edm::InputTag	trueParticlesTag_

bool	useAtHLT
	Use HLT tracking. More...

bool	useConvBremGsfTracks_
	Use Secondary Gsf Tracks. More...

bool	useConvBremPFRecTracks_
	Use Conv Brem KF Tracks. More...

bool	usePFConversions_
	Use of conversions in PFAlgo. More...

bool	usePFNuclearInteractions_
	Use of PFDisplacedVertex in PFAlgo. More...

bool	usePFV0s_
	Use of V0 in PFAlgo. More...

reco::PFV0Collection	v0_

edm::Handle< reco::PFV0Collection >	v0Handle_
	V0. More...

edm::InputTag	v0Tag_

int	verbosity_
	verbosity More...

Detailed Description

Definition at line 63 of file PFRootEventManagerColin.h.

Member Enumeration Documentation

enum PFRootEventManagerColin::Mode

Enumerator
Neutral
HIGH_E_TAUS

Definition at line 80 of file PFRootEventManagerColin.h.

             {
     Neutral=0,
     HIGH_E_TAUS=1
   };

Constructor & Destructor Documentation

PFRootEventManagerColin::PFRootEventManagerColin ( const char * file )

Definition at line 11 of file PFRootEventManagerColin.cc.

References neutralEvent_, outTreeMy_, readSpecificOptions(), and tauEvent_.

   : PFRootEventManager(file) {
   
   tauEvent_ = 0;
   neutralEvent_ = 0;
   outTreeMy_ = 0;   
   
   //   readOptions(file, false, false);
   
   // book histos here
   //   neutralEvent_ = new NeutralEvent();  
 
 
   //   tauEvent_ = new TauEvent();  
   //   outTree_ = new TTree("Tau","");
   //   outTree_->Branch("event","TauEvent", &tauEvent_,32000,2);
 
   readSpecificOptions(file);
 
 }

PFRootEventManagerColin::~PFRootEventManagerColin ( )

Definition at line 32 of file PFRootEventManagerColin.cc.

                                                   {
 
   //   delete event_;
   //   delete outTree_;
 }

Member Function Documentation

bool PFRootEventManagerColin::processEntry ( int entry )

virtual

process one entry (pass the TTree entry)

Reimplemented from PFRootEventManager.

Definition at line 77 of file PFRootEventManagerColin.cc.

References ExpressReco_HICollisions_FallBack::cerr, HIGH_E_TAUS, mode_, Neutral, outTreeMy_, PFRootEventManager::processEntry(), processHIGH_E_TAUS(), processNeutral(), NeutralEvent::rCode, TauEvent::reset(), and tauEvent_.

                                                     {
 
   tauEvent_->reset();
 
   if( ! PFRootEventManager::processEntry(entry) ) {
     // cerr<<"event was not accepted"<<endl;
     // print();
     tauEvent_->rCode = 10;
     return false; // event not accepted
   }
 
   bool rvalue = false;
   switch(mode_) {
   case Neutral:
     // cout<<"colin: process Neutral"<<endl;
     rvalue = processNeutral();
     break;
   case HIGH_E_TAUS:
     // cout<<"colin: process highETaus"<<endl;
     rvalue = processHIGH_E_TAUS();
     break;
   default:
     cerr<<"colin: undefined mode"<<endl;
     assert(0);
   }
   outTreeMy_->Fill();
 
 
   outTreeMy_->Fill();
 
 
   outTreeMy_->Fill();
 
 
   return rvalue;
 }

bool PFRootEventManagerColin::processHIGH_E_TAUS ( )

Definition at line 192 of file PFRootEventManagerColin.cc.

References abs, reco::PFBlock::associatedElements(), Association::block, DeDxDiscriminatorTools::charge(), reco::PFTrack::charge(), ExpressReco_HICollisions_FallBack::chi2, TauEvent::chi2ECAL, reco::PFTrajectoryPoint::ClosestApproach, PFRootEventManager::clustersECAL_, PFRootEventManager::clustersHCAL_, reco::PFSimParticle::daughterIds(), reco::PFBlockElement::ECAL, NeutralEvent::eECAL, NeutralEvent::eHCAL, reco::PFBlock::elements(), asciidump::elements, NeutralEvent::eNeutral, NeutralEvent::etaECAL, NeutralEvent::etaHCAL, NeutralEvent::etaNeutral, TauEvent::etaTrack, reco::PFBlockElement::HCAL, i, getHLTprescales::index, edm::Ref< C, T, F >::isNull(), reco::PFBlock::linkData(), NeutralEvent::nECAL, NeutralEvent::nHCAL, reco::PFSimParticle::pdgCode(), PFRootEventManager::pfBlocks_, TauEvent::pHadron, TauEvent::pTrack, TauEvent::ptTrack, NeutralEvent::rCode, PFRootEventManager::recTracks_, edm::OwnVector< T, P >::size(), tauEvent_, reco::PFBlockElement::TRACK, and PFRootEventManager::trueParticles_.

Referenced by processEntry().

                                                  {
 
 
   // true info 
   // 1 charged hadron, 2 photons
   // save charged part mom, save sum E photons
   
 
   int iHadron  = -1;
   int iPi0 = -1;
   unsigned nStableChargedHadrons=0;
   unsigned nPi0=0;
   for(unsigned i=0; i<trueParticles_.size(); i++) {
     
     const reco::PFSimParticle& part = trueParticles_[i];
     
     int pdgCode = part.pdgCode();
     double charge = part.charge();
 
     if( std::abs(pdgCode) > 100 &&
         charge !=0 && 
         part.daughterIds().empty() ) {
       nStableChargedHadrons++;
       iHadron = i;
     }
     else if( std::abs(pdgCode)==111) {
       nPi0++;
       iPi0 = i; 
     } 
 
     // cout<<i<<" "<<part<<endl;
   }
 
 
   // one has to select 1 charged and 2 photons 
   // to use this function.
 
   // even after filtering events with one stable charged particle,
   // this particle can be a lepton (eg leptonic pion decay)
   if( nStableChargedHadrons==0 ) {
     tauEvent_->rCode = 4; 
     return false;
   }
   assert( nStableChargedHadrons==1 );
 
 
   
   double pHadron = trueParticles_[iHadron].extrapolatedPoint(reco::PFTrajectoryPoint::ClosestApproach ).momentum().P(); 
   tauEvent_->pHadron = pHadron;
 
   //  tauEvent_->eEcalHadron = trueParticles_[iHadron].ecalEnergy();
   
   if(nPi0 == 1) {
     math::XYZTLorentzVector pi0mom =  trueParticles_[iPi0].extrapolatedPoint(reco::PFTrajectoryPoint::ClosestApproach ).momentum();
     tauEvent_->eNeutral = pi0mom.E();
     tauEvent_->etaNeutral = pi0mom.Eta();
   }
   else {
     tauEvent_->eNeutral = 0;
   }
 
   //   if( tauEvent_->eNeutral > 0.1* tauEvent_->pHadron ) {
   //     print();
   //   }
 
 
   // check that there is 
   // only one track
   // 0 or 1 ecal cluster
   // 0 or 1 hcal cluster
 
   if( recTracks_.size() != 1 ) {
     //     cout<<"more than 1 track"<<endl;
     tauEvent_->rCode = 1;
     return false;
   }
   if( clustersECAL_->size() > 1 ) {
     //     cout<<"more than 1 ecal cluster"<<endl;
     tauEvent_->rCode = 2;
     // return false;
   }
   if( clustersHCAL_->size() > 1 ) {
     //     cout<<"more than 1 hcal cluster"<<endl;
     tauEvent_->rCode = 3;
     return false;
   }
   // save track mom + neutral info.
 
   tauEvent_->pTrack = recTracks_[0].extrapolatedPoint(reco::PFTrajectoryPoint::ClosestApproach ).momentum().P();
   tauEvent_->ptTrack = recTracks_[0].extrapolatedPoint(reco::PFTrajectoryPoint::ClosestApproach ).momentum().Pt();
   tauEvent_->etaTrack = recTracks_[0].extrapolatedPoint(reco::PFTrajectoryPoint::ClosestApproach ).momentum().Eta();
 
   // access blocks
 
   // take the track block
 
   // look for the closest associated ecal and hcal clusters
 
   // fill the tree
 
     
 
 
   for(unsigned i=0; i<pfBlocks_->size(); i++) {
     const reco::PFBlock& block = (*pfBlocks_)[i];
     
     const edm::OwnVector< reco::PFBlockElement >& 
       elements = block.elements();
     
     // look for the track
     int iTrack = -1;
     unsigned nTracks = 0;
     for(unsigned ie=0; ie<elements.size(); ie++) {
       if(elements[ie].type() == reco::PFBlockElement::TRACK  ) {
         iTrack = ie;
         nTracks++;
       }
     }
     
     if(nTracks!=1) continue; // no track, or too many tracks in the block
     
     std::multimap<double, unsigned> sortedElems;
     block.associatedElements( iTrack, 
                               block.linkData(),
                               sortedElems );
     
     tauEvent_->nECAL=0;
     tauEvent_->nHCAL=0;
     
     typedef std::multimap<double, unsigned>::iterator IE;
     for(IE ie = sortedElems.begin(); ie != sortedElems.end(); ++ie ) {
       
       
       double chi2 = ie->first;
       unsigned index = ie->second;
       
       reco::PFBlockElement::Type type = elements[index].type();
       
       reco::PFClusterRef clusterRef = elements[index].clusterRef();
 
 
       
 
       if( type == reco::PFBlockElement::ECAL ) {
         if(!tauEvent_->nECAL ) { // closest ecal
           assert( !clusterRef.isNull() );
           tauEvent_->eECAL = clusterRef->energy();
           tauEvent_->etaECAL = clusterRef->position().Eta();
           tauEvent_->chi2ECAL = chi2;
           tauEvent_->nECAL++;
         }
       }
 
 
       else if( type == reco::PFBlockElement::HCAL ) {
         if(!tauEvent_->nHCAL ) { // closest hcal
           assert( !clusterRef.isNull() );
           tauEvent_->eHCAL = clusterRef->energy();
           tauEvent_->etaHCAL = clusterRef->position().Eta();
           tauEvent_->nHCAL++;
         }
       } 
     } // eles associated to the track
   } // blocks
 
 
 
 
   return false;    
 }

bool PFRootEventManagerColin::processNeutral ( )

Definition at line 117 of file PFRootEventManagerColin.cc.

References ExpressReco_HICollisions_FallBack::cerr, NeutralEvent::eECAL, NeutralEvent::eHCAL, NeutralEvent::eNeutral, eta(), NeutralEvent::etaNeutral, configurableAnalysis::GenParticle, edm::HepMCProduct::GetEvent(), i, PFRootEventManager::MCTruth_, NeutralEvent::nECAL, neutralEvent_, NeutralEvent::nHCAL, outTreeMy_, phi, and NeutralEvent::reset().

Referenced by processEntry().

                                              {
   //   else {
   //     cerr<<"event accepted"<<endl;
   //   }
   
   //   if( ! ( (*clustersECAL_).size() <= 1 && 
   //      (*clustersHCAL_).size() <= 1 ) ) {
   //     cerr<<"wrong number of ECAL or HCAL clusters :"
   //    <<(*clustersECAL_).size()<<","<<(*clustersHCAL_).size()<<endl;
   //     return false; 
   //   }
   // 1 HCAL cluster
   
   neutralEvent_->reset();
   
   // particle
 
   const HepMC::GenEvent* myGenEvent = MCTruth_.GetEvent();
   if(!myGenEvent) {
     assert(0);
   }
 
   if( myGenEvent->particles_size() != 1 ) {
     cerr<<"wrong number of particles:"
         <<myGenEvent->particles_size()<<endl;
     return 0;
   }
 
   // take first particle  
   const HepMC::GenParticle* particle = *(myGenEvent->particles_begin() );
   
   // and check that it's a K0L 
   if( particle->pdg_id() != 130 ) {
     cerr<<"not a K0L : "<<particle->pdg_id()<<endl;
     return false;
   }
 
   neutralEvent_->eNeutral = particle->momentum().e();
 
   double eta =  particle->momentum().eta();
   double phi =  particle->momentum().phi();
   neutralEvent_->etaNeutral = eta;
   
 
   neutralEvent_->nECAL = (*clustersECAL_).size();
 
   // look for the closest ECAL cluster from the particle.
 
   double minDist2 = 9999999;
   // int iClosest = -1;
   for( unsigned i=0; i<(*clustersECAL_).size(); ++i) {
     double deta = (*clustersECAL_)[i].position().Eta() - eta;
     double dphi = (*clustersECAL_)[i].position().Phi() - phi;
     double dist2 = deta*deta + dphi*dphi;
     if(dist2 < minDist2) {
       minDist2 = dist2;
       neutralEvent_->eECAL = (*clustersECAL_)[i].energy();
     }
   }
 
 
   neutralEvent_->nHCAL = (*clustersHCAL_).size();
   if( (*clustersHCAL_).size() == 1 ) {
     neutralEvent_->eHCAL = (*clustersHCAL_)[0].energy();
   }
 
   
   outTreeMy_->Fill();
   
   if(  neutralEvent_->nECAL<1 && neutralEvent_->eNeutral<1 )
     return true;
   else return false;
 }

void PFRootEventManagerColin::readSpecificOptions ( const char * file )

virtual

Reimplemented from PFRootEventManager.

Definition at line 39 of file PFRootEventManagerColin.cc.

References ExpressReco_HICollisions_FallBack::cerr, gather_cfg::cout, cmsRelvalreport::exit, IO::GetOpt(), HIGH_E_TAUS, mode_, Neutral, neutralEvent_, PFRootEventManager::options_, PFRootEventManager::outFile_, outTreeMy_, and tauEvent_.

Referenced by PFRootEventManagerColin().

                                                                   {
 
 
 
   cout<<"calling PFRootEventManagerColin::readSpecificOptions"<<endl; 
   //   PFRootEventManager::readOptions(file, refresh, reconnect);
 
   
   options_->GetOpt("colin", "mode", mode_);
 
   if(outTreeMy_) delete outTreeMy_;
   
   outFile_->cd();
   switch(mode_) {
   case Neutral:
     cout<<"colin: Neutral mode"<<endl;
     neutralEvent_ = new NeutralEvent();  
     outTreeMy_ = new TTree("Neutral","");
     outTreeMy_->Branch("event","NeutralEvent", &neutralEvent_,32000,2);
     gDirectory->ls();
     break;
   case HIGH_E_TAUS:
     cout<<"colin: highETaus mode"<<endl;
     tauEvent_ = new TauEvent();  
     outTreeMy_ = new TTree("Tau","");
     outTreeMy_->Branch("event","TauEvent", &tauEvent_,32000,2);
     gDirectory->ls();
     break;
   default:
     cerr<<"colin: undefined mode"<<endl;
     exit(1);
   }
 }

void PFRootEventManagerColin::write ( void )

virtual

Reimplemented from PFRootEventManager.

Definition at line 366 of file PFRootEventManagerColin.cc.

References PFRootEventManager::outFile_, outTreeMy_, and PFRootEventManager::write().

                                     {
   // write histos here
   outFile_->cd();
   outTreeMy_->Write();
 
   PFRootEventManager::write();
 }

Member Data Documentation

int PFRootEventManagerColin::mode_

Definition at line 88 of file PFRootEventManagerColin.h.

Referenced by processEntry(), and readSpecificOptions().

NeutralEvent* PFRootEventManagerColin::neutralEvent_

Definition at line 86 of file PFRootEventManagerColin.h.

Referenced by PFRootEventManagerColin(), processNeutral(), and readSpecificOptions().

TTree* PFRootEventManagerColin::outTreeMy_

Definition at line 85 of file PFRootEventManagerColin.h.

Referenced by PFRootEventManagerColin(), processEntry(), processNeutral(), readSpecificOptions(), and write().

TauEvent* PFRootEventManagerColin::tauEvent_

Definition at line 87 of file PFRootEventManagerColin.h.

Referenced by PFRootEventManagerColin(), processEntry(), processHIGH_E_TAUS(), and readSpecificOptions().

Public Types

Public Member Functions

Public Attributes

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation