1 |
|
// $Id$ |
2 |
|
|
3 |
+ |
#include "MitProd/TreeFiller/interface/FillerConversionElectrons.h" |
4 |
|
#include "FWCore/MessageLogger/interface/MessageLogger.h" |
4 |
– |
#include "FWCore/Framework/interface/ESHandle.h" |
5 |
|
#include "DataFormats/Common/interface/Handle.h" |
6 |
– |
#include "FWCore/ServiceRegistry/interface/Service.h" |
6 |
|
#include "DataFormats/TrackReco/interface/Track.h" |
8 |
– |
#include "MitAna/DataTree/interface/Electron.h" |
7 |
|
#include "DataFormats/TrackReco/interface/TrackFwd.h" |
10 |
– |
|
11 |
– |
|
12 |
– |
#include "MitAna/DataTree/interface/Track.h" |
8 |
|
#include "MitAna/DataTree/interface/Names.h" |
14 |
– |
#include "MitProd/TreeFiller/interface/FillerConversionElectrons.h" |
9 |
|
|
10 |
|
using namespace std; |
11 |
|
using namespace edm; |
12 |
|
using namespace mithep; |
13 |
|
|
14 |
< |
//------------------------------------------------------------------------------------------------- |
15 |
< |
FillerConversionElectrons::FillerConversionElectrons(const edm::ParameterSet& cfg, bool active, const mithep::TrackCol* conversionInOutTracks, const mithep::TrackCol* conversionOutInTracks, const mithep::TrackMap* conversionInOutTrackMap, const mithep::TrackMap* conversionOutInTrackMap) : |
16 |
< |
BaseFiller(cfg, "ConversionElectrons", active), |
17 |
< |
mitName_(Conf().getUntrackedParameter<string>("mitName","")), |
18 |
< |
conversionElectrons_(new mithep::ElectronArr), |
14 |
> |
//-------------------------------------------------------------------------------------------------- |
15 |
> |
FillerConversionElectrons::FillerConversionElectrons(const edm::ParameterSet &cfg, bool active, |
16 |
> |
const TrackCol *convInOutTracks, |
17 |
> |
const TrackCol *convOutInTracks, |
18 |
> |
const TrackMap *convInOutTrackMap, |
19 |
> |
const TrackMap *convOutInTrackMap) : |
20 |
> |
BaseFiller(cfg,"ConversionElectrons",active), |
21 |
> |
mitName_(Conf().getUntrackedParameter<string>("mitName","ConversionElectrons")), |
22 |
> |
conversionElectrons_(new mithep::ElectronArr(16)), |
23 |
|
convElectronMap_(new mithep::ConversionElectronMap), |
24 |
< |
conversionInOutTracks_(conversionInOutTracks), |
25 |
< |
conversionOutInTracks_(conversionOutInTracks), |
26 |
< |
conversionInOutTrackMap_(conversionInOutTrackMap), |
27 |
< |
conversionOutInTrackMap_(conversionOutInTrackMap) |
24 |
> |
conversionInOutTracks_(convInOutTracks), |
25 |
> |
conversionOutInTracks_(convOutInTracks), |
26 |
> |
conversionInOutTrackMap_(convInOutTrackMap), |
27 |
> |
conversionOutInTrackMap_(convOutInTrackMap) |
28 |
|
{ |
29 |
+ |
// Constructor. |
30 |
|
} |
31 |
|
|
32 |
< |
//------------------------------------------------------------------------------------------------- |
32 |
> |
//-------------------------------------------------------------------------------------------------- |
33 |
|
FillerConversionElectrons::~FillerConversionElectrons() |
34 |
|
{ |
35 |
+ |
// Destructor. |
36 |
+ |
|
37 |
+ |
delete conversionElectrons_; |
38 |
|
delete convElectronMap_; |
39 |
|
} |
40 |
|
|
41 |
< |
//------------------------------------------------------------------------------------------------- |
41 |
> |
//-------------------------------------------------------------------------------------------------- |
42 |
|
void FillerConversionElectrons::BookDataBlock(TreeWriter &tws) |
43 |
|
{ |
44 |
+ |
// Add converted electron branch to tree. |
45 |
+ |
|
46 |
|
tws.AddBranch(mitName_.c_str(),&conversionElectrons_); |
47 |
|
} |
48 |
|
|
49 |
< |
//------------------------------------------------------------------------------------------------- |
49 |
> |
//-------------------------------------------------------------------------------------------------- |
50 |
|
void FillerConversionElectrons::FillDataBlock(const edm::Event &event, |
51 |
< |
const edm::EventSetup &setup) |
51 |
> |
const edm::EventSetup &setup) |
52 |
|
{ |
53 |
< |
// Fill global muon information. Links to the various standalone components are created in the |
50 |
< |
// link resolution. |
53 |
> |
// Fill converted electron array. |
54 |
|
|
55 |
|
conversionElectrons_->Reset(); |
56 |
|
convElectronMap_->Reset(); |
57 |
|
FillFromTracks(conversionInOutTracks_, conversionInOutTrackMap_); |
58 |
|
FillFromTracks(conversionOutInTracks_, conversionOutInTrackMap_); |
59 |
|
conversionElectrons_->Trim(); |
57 |
– |
|
58 |
– |
|
59 |
– |
} |
60 |
– |
|
61 |
– |
void FillerConversionElectrons::FillFromTracks(const mithep::TrackCol* tracks, const TrackMap* trackMap) { |
62 |
– |
for (UInt_t i=0; i<tracks->GetEntries(); i++) { |
63 |
– |
mithep::Track* track = (Track*)tracks->At(i); |
64 |
– |
mithep::Electron* electron = conversionElectrons_->AddNew(); |
65 |
– |
electron->SetTrackerTrack(track); |
66 |
– |
convElectronMap_->Add(trackMap->GetEdmRef(track), electron); |
67 |
– |
} |
60 |
|
} |
61 |
|
|
62 |
+ |
//-------------------------------------------------------------------------------------------------- |
63 |
+ |
void FillerConversionElectrons::FillFromTracks(const TrackCol *tracks, |
64 |
+ |
const TrackMap *trackMap) |
65 |
+ |
{ |
66 |
+ |
// Obtain information from tracks. |
67 |
|
|
68 |
+ |
for (UInt_t i=0; i<tracks->GetEntries(); ++i) { |
69 |
+ |
mithep::Track *track = const_cast<Track*>(tracks->At(i)); |
70 |
+ |
mithep::Electron *electron = conversionElectrons_->AddNew(); |
71 |
+ |
electron->SetTrackerTrk(track); |
72 |
+ |
convElectronMap_->Add(trackMap->GetEdm(track), electron); |
73 |
+ |
} |
74 |
+ |
} |