126 |
|
#include "Visualisation/Frog/soft/Includes/FROG/FROG_DetId.h" |
127 |
|
#include "Visualisation/Frog/soft/Includes/FROG/FROG_ReadCards.cpp" |
128 |
|
|
129 |
+ |
#include "DataFormats/VertexReco/interface/NuclearInteraction.h" |
130 |
+ |
|
131 |
+ |
#include "Visualisation/Frog/soft/Includes/FROG/FROG_Element_Event_NuclInt.h" |
132 |
+ |
|
133 |
|
|
134 |
|
using namespace edm; |
135 |
|
using namespace std; |
156 |
|
std::vector<InputTag> SimTrackProducers; |
157 |
|
std::vector<InputTag> SimVertexProducers; |
158 |
|
std::vector<InputTag> SimHitProducers; |
159 |
+ |
std::vector<InputTag> NIProducers; |
160 |
|
|
161 |
|
std::vector<InputTag> TrackProducers; |
162 |
|
std::vector<InputTag> TrajectoryProducers; |
197 |
|
HcalHORecHitProducers = iConfig.getParameter<std::vector<InputTag> >("HcalHORecHitProducers"); |
198 |
|
HcalHFRecHitProducers = iConfig.getParameter<std::vector<InputTag> >("HcalHFRecHitProducers"); |
199 |
|
|
200 |
+ |
NIProducers = iConfig.getParameter<std::vector<InputTag> >("NIProducers"); |
201 |
+ |
|
202 |
|
DTSegmentProducers = iConfig.getParameter<std::vector<InputTag> >("DTSegmentProducers"); |
203 |
|
CSCSegmentProducers = iConfig.getParameter<std::vector<InputTag> >("CSCSegmentProducers"); |
204 |
|
|
740 |
|
} |
741 |
|
|
742 |
|
|
743 |
< |
FROG_Element_Event* event = new FROG_Element_Event(); |
743 |
> |
FROG_Element_Event* event = new FROG_Element_Event(iEvent.id().run(), iEvent.id().event() ); |
744 |
|
events_prim->addDaughter(event); |
745 |
|
|
746 |
|
// access the tracker |
820 |
|
} |
821 |
|
} |
822 |
|
|
823 |
+ |
// NUCLEAR INTERACTIONS |
824 |
+ |
for(unsigned int i=0;i<NIProducers.size();i++) { |
825 |
+ |
edm::Handle<std::vector<reco::NuclearInteraction> > h_NIs; |
826 |
+ |
iEvent.getByLabel(NIProducers[i], h_NIs); |
827 |
+ |
std::vector<reco::NuclearInteraction> NIs = *h_NIs.product(); |
828 |
+ |
FROG_Element_Base_With_DetId* frogNIColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_NUCLINT + i); |
829 |
+ |
event->addDaughter(frogNIColl); |
830 |
+ |
|
831 |
+ |
for(unsigned int h=0; h<NIs.size();h++){ |
832 |
+ |
reco::NuclearInteraction ni = NIs[i]; |
833 |
+ |
float vtx_x = ni.vertex().x(); |
834 |
+ |
float vtx_y = ni.vertex().y(); |
835 |
+ |
float vtx_z = ni.vertex().z(); |
836 |
+ |
float like = ni.likelihood(); |
837 |
+ |
FROG_Element_Event_NuclInt* frogNI = new FROG_Element_Event_NuclInt(vtx_x, vtx_y, vtx_z, like ); |
838 |
+ |
frogNIColl->addDaughter(frogNI); |
839 |
+ |
} |
840 |
+ |
} |
841 |
|
|
842 |
|
|
843 |
|
// ### TRAJECTORY ### |