ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/FastOpenGlDisplayer/src/Frog_Analyzer.cc
Revision: 1.23
Committed: Fri Jul 11 04:38:40 2008 UTC (16 years, 10 months ago) by querten
Content type: text/plain
Branch: MAIN
CVS Tags: Version_1_02
Changes since 1.22: +2 -2 lines
Log Message: 
Add Run and Event Number

File Contents

# User Rev Content
1 querten 1.1 // -*- C++ -*-
2     //
3     // Package: Frog_Analyzer
4     // Class: Frog_Analyzer
5     //
6     /**\class Frog_Analyzer Frog_Analyzer.cc Visualisation/Frog/src/Frog_Analyzer.cc
7    
8     Description: <one line class summary>
9    
10     Implementation:
11     <Notes on implementation>
12     */
13     //
14     // Original Author: Loic QUERTENMONT
15     // Created: Fri Oct 26 07:22:12 CEST 2007
16 querten 1.23 // $Id: Frog_Analyzer.cc,v 1.22 2008/07/09 08:49:50 querten Exp $
17 querten 1.1 //
18     //
19    
20    
21     // system include files
22     #include <memory>
23    
24     // user include files
25     #include "FWCore/Framework/interface/Frameworkfwd.h"
26     #include "FWCore/Framework/interface/EDAnalyzer.h"
27    
28     #include "FWCore/Framework/interface/Event.h"
29     #include "FWCore/Framework/interface/MakerMacros.h"
30    
31     #include "FWCore/ParameterSet/interface/ParameterSet.h"
32     #include "FWCore/ServiceRegistry/interface/Service.h"
33     #include "FWCore/Framework/interface/ESHandle.h"
34    
35     #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
36     #include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"
37     #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
38     #include "Geometry/TrackerNumberingBuilder/interface/GeometricDet.h"
39     #include "Geometry/CommonTopologies/interface/PixelTopology.h"
40     #include "Geometry/CommonTopologies/interface/StripTopology.h"
41     #include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetType.h"
42     #include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetType.h"
43     #include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
44     #include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
45     #include "DataFormats/GeometrySurface/interface/BoundSurface.h"
46     #include "DataFormats/DetId/interface/DetId.h"
47    
48     #include "Geometry/CaloEventSetup/interface/CaloTopologyRecord.h"
49     #include "Geometry/Records/interface/IdealGeometryRecord.h"
50     #include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
51     #include "Geometry/CaloGeometry/interface/CaloGeometry.h"
52     #include "Geometry/EcalBarrelAlgo/interface/EcalBarrelGeometry.h"
53     #include "Geometry/EcalEndcapAlgo/interface/EcalEndcapGeometry.h"
54     #include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
55     #include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
56    
57     #include "Geometry/CommonDetUnit/interface/TrackingGeometry.h"
58     #include "Geometry/Records/interface/MuonGeometryRecord.h"
59     #include "Geometry/DTGeometry/interface/DTGeometry.h"
60     #include "Geometry/CSCGeometry/interface/CSCGeometry.h"
61     #include "Geometry/RPCGeometry/interface/RPCGeometry.h"
62    
63    
64     #include "DataFormats/CaloRecHit/interface/CaloRecHit.h"
65     #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
66     #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
67    
68     #include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
69     #include "SimDataFormats/TrackingHit/interface/PSimHit.h"
70     #include "SimDataFormats/CrossingFrame/interface/CrossingFrame.h"
71     #include "SimDataFormats/CrossingFrame/interface/MixCollection.h"
72     #include "SimDataFormats/CaloHit/interface/PCaloHit.h"
73     #include "SimDataFormats/Track/interface/SimTrack.h"
74     #include "SimDataFormats/Vertex/interface/SimVertex.h"
75    
76     #include "DataFormats/Provenance/interface/BranchDescription.h"
77     #include "DataFormats/Provenance/interface/Provenance.h"
78     #include "DataFormats/Candidate/interface/Candidate.h"
79     #include "DataFormats/HepMCCandidate/interface/GenParticleCandidate.h"
80     #include "DataFormats/RecoCandidate/interface/RecoChargedCandidate.h"
81     #include "DataFormats/TrackReco/interface/Track.h"
82 roberfro 1.5 #include "TrackingTools/PatternTools/interface/Trajectory.h"
83 querten 1.6 #include "TrackingTools/PatternTools/interface/TrajTrackAssociation.h"
84    
85 querten 1.1
86 querten 1.2 #include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"
87     #include "DataFormats/CSCRecHit/interface/CSCSegmentCollection.h"
88    
89     #include "DataFormats/DTRecHit/interface/DTRecHitCollection.h"
90     #include "DataFormats/CSCRecHit/interface/CSCRecHit2DCollection.h"
91     #include "DataFormats/RPCRecHit/interface/RPCRecHitCollection.h"
92    
93 querten 1.1 #include "DataFormats/Math/interface/Point3D.h"
94     #include "DataFormats/GeometrySurface/interface/TrapezoidalPlaneBounds.h"
95     #include "DataFormats/GeometrySurface/interface/RectangularPlaneBounds.h"
96    
97 querten 1.11 #include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
98     #include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
99     #include "DataFormats/SiPixelDetId/interface/PXBDetId.h"
100     #include "DataFormats/SiPixelDetId/interface/PXFDetId.h"
101     #include "DataFormats/SiStripDetId/interface/TIBDetId.h"
102     #include "DataFormats/SiStripDetId/interface/TIDDetId.h"
103     #include "DataFormats/SiStripDetId/interface/TOBDetId.h"
104     #include "DataFormats/SiStripDetId/interface/TECDetId.h"
105    
106     #include "DataFormats/MuonDetId/interface/DTChamberId.h"
107     #include "DataFormats/MuonDetId/interface/CSCDetId.h"
108     #include "DataFormats/MuonDetId/interface/RPCDetId.h"
109    
110     #include "DataFormats/EcalDetId/interface/EBDetId.h"
111     #include "DataFormats/EcalDetId/interface/EEDetId.h"
112     #include "DataFormats/EcalDetId/interface/ESDetId.h"
113     #include "DataFormats/HcalDetId/interface/HcalDetId.h"
114    
115     #include "Visualisation/Frog/soft/Includes/FROG/FROG_Element_Base_With_DetId.h"
116     #include "Visualisation/Frog/soft/Includes/FROG/FROG_Element_Tools.h"
117     #include "Visualisation/Frog/soft/Includes/FROG/FROG_Element_Tools.cpp"
118 roberfro 1.13 #include "Visualisation/Frog/soft/Includes/FROG/FROG_Element_Event_Track.h"
119    
120     #include "Visualisation/Frog/soft/Includes/FROG/FROG_Events.h"
121     #include "Visualisation/Frog/soft/Includes/FROG/FROG_Events.cpp"
122    
123     #include "Visualisation/Frog/soft/Includes/FROG/FROG_Geometry.h"
124     #include "Visualisation/Frog/soft/Includes/FROG/FROG_Geometry.cpp"
125 querten 1.11
126 roberfro 1.13 #include "Visualisation/Frog/soft/Includes/FROG/FROG_DetId.h"
127 querten 1.18 #include "Visualisation/Frog/soft/Includes/FROG/FROG_ReadCards.cpp"
128 querten 1.1
129 querten 1.11
130 querten 1.1 using namespace edm;
131     using namespace std;
132    
133     //
134     // class decleration
135     //
136    
137     class Frog_Analyzer : public edm::EDAnalyzer {
138     public:
139     explicit Frog_Analyzer(const edm::ParameterSet&);
140     ~Frog_Analyzer();
141    
142    
143     private:
144     virtual void beginJob(const edm::EventSetup& iSetup);
145     virtual void analyze (const edm::Event&, const edm::EventSetup&);
146     virtual void endJob ();
147    
148 querten 1.15 FROG_Element_Base* events_prim;
149 querten 1.1
150 querten 1.11 std::string OutputFile;
151 querten 1.1
152     std::vector<InputTag> SimTrackProducers;
153     std::vector<InputTag> SimVertexProducers;
154     std::vector<InputTag> SimHitProducers;
155    
156 querten 1.9 std::vector<InputTag> TrackProducers;
157 roberfro 1.5 std::vector<InputTag> TrajectoryProducers;
158 querten 1.1 std::vector<InputTag> EcalRecHitProducers;
159     std::vector<InputTag> HcalHBHERecHitProducers;
160     std::vector<InputTag> HcalHORecHitProducers;
161     std::vector<InputTag> HcalHFRecHitProducers;
162 querten 1.2 std::vector<InputTag> DTSegmentProducers;
163     std::vector<InputTag> CSCSegmentProducers;
164    
165     std::vector<InputTag> RPCHitsProducers;
166 querten 1.1
167 querten 1.22 // ----------member data ---------------------------
168    
169     bool ProduceGeom;
170     int NEventsInVisFile;
171     int NEvents;
172 querten 1.1
173 querten 1.22 unsigned int SRun;
174     unsigned int SEvent;
175 querten 1.1 };
176    
177     //
178     // constructors and destructor
179     //
180     Frog_Analyzer::Frog_Analyzer(const edm::ParameterSet& iConfig)
181     {
182     OutputFile = iConfig.getParameter<std::string >("OutputFile");
183    
184     SimTrackProducers = iConfig.getParameter<std::vector<InputTag> >("SimTrackProducers");
185     SimVertexProducers = iConfig.getParameter<std::vector<InputTag> >("SimVertexProducers");
186     SimHitProducers = iConfig.getParameter<std::vector<InputTag> >("SimHitProducers");
187    
188 querten 1.10 TrackProducers = iConfig.getParameter<std::vector<InputTag> >("TrackProducers");
189 querten 1.9 TrajectoryProducers = iConfig.getParameter<std::vector<InputTag> >("TrajectoryProducers");
190 querten 1.1 EcalRecHitProducers = iConfig.getParameter<std::vector<InputTag> >("EcalRecHitProducers");
191     HcalHBHERecHitProducers = iConfig.getParameter<std::vector<InputTag> >("HcalHBHERecHitProducers");
192     HcalHORecHitProducers = iConfig.getParameter<std::vector<InputTag> >("HcalHORecHitProducers");
193     HcalHFRecHitProducers = iConfig.getParameter<std::vector<InputTag> >("HcalHFRecHitProducers");
194    
195 querten 1.2 DTSegmentProducers = iConfig.getParameter<std::vector<InputTag> >("DTSegmentProducers");
196     CSCSegmentProducers = iConfig.getParameter<std::vector<InputTag> >("CSCSegmentProducers");
197    
198     RPCHitsProducers = iConfig.getParameter<std::vector<InputTag> >("RPCHitsProducers");
199    
200 querten 1.22 ProduceGeom = iConfig.getParameter<bool >("ProduceGeom" );
201     NEventsInVisFile = iConfig.getParameter<int >("NEventsInVisFile" );
202    
203     NEvents = 0;
204 querten 1.1 }
205    
206    
207     Frog_Analyzer::~Frog_Analyzer()
208     {
209     }
210    
211     // ------------ method called once each job just before starting event loop ------------
212     void
213     Frog_Analyzer::beginJob(const edm::EventSetup& iSetup)
214     {
215 querten 1.22 if(ProduceGeom){
216     DetId Detid;
217     unsigned int SubDet;
218    
219     unsigned int Frog_DetId;
220     FROG_Element_Base* FEB = NULL;
221     FROG_Element_Base_With_DetId* FEB_Temp = NULL;
222    
223     FROG_Element_Base* prim = new FROG_Element_Base(C_PRIMARY);
224     FROG_Element_Base* geom = new FROG_Element_Base(C_GEOMETRY); prim->addDaughter(geom);
225     FROG_Element_Base_With_DetId* tracker = new FROG_Element_Base_With_DetId(C_FEB_DETID, DETID_TRK); geom->addDaughter(tracker);
226     FROG_Element_Base_With_DetId* muon = new FROG_Element_Base_With_DetId(C_FEB_DETID, DETID_MUON); geom->addDaughter(muon);
227     FROG_Element_Base_With_DetId* ecal = new FROG_Element_Base_With_DetId(C_FEB_DETID, DETID_ECAL); geom->addDaughter(ecal);
228     FROG_Element_Base_With_DetId* hcal = new FROG_Element_Base_With_DetId(C_FEB_DETID, DETID_HCAL); geom->addDaughter(hcal);
229    
230    
231     // ### TRACKER GEOMETRY ###
232    
233     edm::ESHandle<TrackerGeometry> tkGeom;
234     iSetup.get<TrackerDigiGeometryRecord>().get( tkGeom );
235     vector<GeomDet*> TrackerDets = tkGeom->dets();
236    
237     for(unsigned int i=0;i<TrackerDets.size();i++){
238     Detid = TrackerDets[i]->geographicalId();
239     SubDet = Detid.subdetId();
240     if(SubDet<1 || SubDet>6) continue;
241 querten 1.11
242 querten 1.1
243 querten 1.22 GeomDet* DetUnit = TrackerDets[i];
244     if(!DetUnit)continue;
245     const BoundPlane plane = DetUnit->surface();
246     const TrapezoidalPlaneBounds* trapezoidalBounds( dynamic_cast<const TrapezoidalPlaneBounds*>(&(plane.bounds())));
247     const RectangularPlaneBounds* rectangularBounds( dynamic_cast<const RectangularPlaneBounds*>(&(plane.bounds())));
248    
249     float width = 0;
250     float length = 0;
251     float thickness = 0;
252     float TrapezoidalParam = 0;
253    
254     if(trapezoidalBounds){
255     std::vector<float> const & parameters = (*trapezoidalBounds).parameters();
256     width = parameters[0]*2;
257     length = parameters[3]*2;
258     thickness = (*trapezoidalBounds).thickness();
259     TrapezoidalParam = parameters[1]/parameters[0];
260     }else if(rectangularBounds){
261     width = DetUnit->surface().bounds().width();
262     length = DetUnit->surface().bounds().length();
263     thickness = DetUnit->surface().bounds().thickness();
264     TrapezoidalParam = 1;
265     }
266 querten 1.1
267 querten 1.22 Surface::GlobalPoint WidthVector = plane.toGlobal( LocalPoint(width/2, 0, 0) );
268     Surface::GlobalPoint LengthVector = plane.toGlobal( LocalPoint(0, length/2, 0) );
269     Surface::GlobalPoint ThickVector = plane.toGlobal( LocalPoint(0, 0, thickness/2) );
270    
271     GlobalVector Pos = GlobalVector(DetUnit->position().basicVector());
272    
273     FROG_Element_Geom_TrackerMod* mod = new FROG_Element_Geom_TrackerMod(
274     Detid.rawId(), TrapezoidalParam,
275     Pos.x(), Pos.y(), Pos.z(),
276     WidthVector.x() -Pos.x(), WidthVector.y() -Pos.y(), WidthVector.z() -Pos.z(),
277     LengthVector.x()-Pos.x(), LengthVector.y()-Pos.y(), LengthVector.z()-Pos.z(),
278     ThickVector.x() -Pos.x(), ThickVector.y() -Pos.y(), ThickVector.z() -Pos.z());
279 querten 1.1
280    
281 querten 1.22 Frog_DetId = DETID_TRK;
282     FEB = tracker;
283 querten 1.11
284 querten 1.22 Frog_DetId += SubDet*100000;
285 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
286     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
287     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
288    
289 querten 1.22 if(SubDet==PixelSubdetector::PixelBarrel){
290     PXBDetId detid = PXBDetId(Detid);
291 querten 1.11
292 querten 1.22 Frog_DetId += detid.layer()*10000;
293     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
294     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
295     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
296 querten 1.11
297 querten 1.22 Frog_DetId += detid.ladder()*100;
298     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
299     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
300     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
301     }else if(SubDet==PixelSubdetector::PixelEndcap){
302     PXFDetId detid = PXFDetId(Detid);
303 querten 1.11
304 querten 1.22 Frog_DetId += detid.side()*10000;
305     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
306     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
307     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
308 querten 1.11
309 querten 1.22 Frog_DetId += detid.disk()*1000;
310     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
311     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
312     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
313 querten 1.11
314 querten 1.22 Frog_DetId += detid.blade()*10;
315     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
316     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
317     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
318 querten 1.11
319 querten 1.22 Frog_DetId += detid.panel()*1;
320     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
321     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
322     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
323     }else if(SubDet==StripSubdetector::TIB){
324     TIBDetId detid = TIBDetId(Detid);
325 querten 1.11
326 querten 1.22 Frog_DetId += detid.layer()*10000;
327     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
328     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
329     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
330 querten 1.11
331 querten 1.22 Frog_DetId += detid.stringNumber()*100;
332     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
333     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
334     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
335     }else if(SubDet==StripSubdetector::TID){
336     TIDDetId detid = TIDDetId(Detid);
337 querten 1.11
338 querten 1.22 Frog_DetId += detid.side()*10000;
339     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
340     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
341     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
342 querten 1.11
343 querten 1.22 Frog_DetId += detid.wheel()*1000;
344     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
345     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
346     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
347    
348     Frog_DetId += detid.ring()*10;
349     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
350     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
351     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
352     }else if(SubDet==StripSubdetector::TOB){
353     TOBDetId detid = TOBDetId(Detid);
354 querten 1.11
355 querten 1.22 Frog_DetId += detid.layer()*10000;
356     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
357     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
358     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
359 querten 1.11
360 querten 1.22 Frog_DetId += detid.rodNumber()*100;
361     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
362     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
363     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
364     }else if(SubDet==StripSubdetector::TEC){
365     TECDetId detid = TECDetId(Detid);
366 querten 1.11
367 querten 1.22 Frog_DetId += detid.side()*10000;
368     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
369     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
370     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
371 querten 1.11
372 querten 1.22 Frog_DetId += detid.wheel()*1000;
373     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
374     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
375     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
376 querten 1.1
377 querten 1.22 Frog_DetId += detid.ring()*10;
378 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
379     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
380     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
381 querten 1.22
382     Frog_DetId += detid.petalNumber();
383 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
384     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
385     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
386 querten 1.22 }else{
387     if(mod)delete mod;
388     continue;
389 querten 1.11 }
390 querten 1.22 FEB->addDaughter(mod);
391     }
392 querten 1.11
393    
394 querten 1.22 // ### MUON GEOMETRY ###
395 querten 1.11
396 querten 1.22 edm::ESHandle<DTGeometry> DtGeom;
397     iSetup.get<MuonGeometryRecord>().get( DtGeom );
398     const vector<GeomDet*> DtDets = DtGeom->dets();
399    
400     edm::ESHandle<CSCGeometry> CscGeom;
401     iSetup.get<MuonGeometryRecord>().get( CscGeom );
402     const vector<GeomDet*> CscDets = CscGeom->dets();
403    
404     edm::ESHandle<RPCGeometry> RpcGeom;
405     iSetup.get<MuonGeometryRecord>().get( RpcGeom );
406     const vector<GeomDet*> RpcDets = RpcGeom->dets();
407    
408     vector<GeomDet*> MuonDets;
409     for(unsigned int i=0;i<DtDets.size() ;i++){MuonDets.push_back(DtDets [i]);}
410     for(unsigned int i=0;i<CscDets.size();i++){MuonDets.push_back(CscDets[i]);}
411     for(unsigned int i=0;i<RpcDets.size();i++){MuonDets.push_back(RpcDets[i]);}
412 querten 1.1
413 querten 1.22 for(unsigned int i=0;i<MuonDets.size();i++)
414     {
415     Detid = DetId(MuonDets[i]->geographicalId());
416     SubDet = Detid.subdetId();
417 querten 1.1
418 querten 1.22 GeomDet* DetUnit = MuonDets[i];
419     if(!DetUnit)continue;
420     const BoundPlane plane = DetUnit->surface();
421     const TrapezoidalPlaneBounds* trapezoidalBounds( dynamic_cast<const TrapezoidalPlaneBounds*>(&(plane.bounds())));
422     const RectangularPlaneBounds* rectangularBounds( dynamic_cast<const RectangularPlaneBounds*>(&(plane.bounds())));
423    
424     float width = 0;
425     float length = 0;
426     float thickness = 0;
427     float TrapezoidalParam = 0;
428    
429     if(trapezoidalBounds)
430     {
431     std::vector<float> const & parameters = (*trapezoidalBounds).parameters();
432     width = parameters[0]*2;
433     length = parameters[3]*2;
434     thickness = (*trapezoidalBounds).thickness();
435     TrapezoidalParam = parameters[1]/parameters[0];
436     }else if(rectangularBounds){
437     width = DetUnit->surface().bounds().width();
438     length = DetUnit->surface().bounds().length();
439     thickness = DetUnit->surface().bounds().thickness();
440     TrapezoidalParam = 1;
441     }
442 querten 1.1
443 querten 1.22 Surface::GlobalPoint WidthVector = plane.toGlobal( LocalPoint(width/2, 0, 0) );
444     Surface::GlobalPoint LengthVector = plane.toGlobal( LocalPoint(0, length/2, 0) );
445     Surface::GlobalPoint ThickVector = plane.toGlobal( LocalPoint(0, 0, thickness/2) );
446    
447     GlobalVector Pos = GlobalVector(DetUnit->position().basicVector());
448    
449     FROG_Element_Geom_TrackerMod* mod = new FROG_Element_Geom_TrackerMod(
450     Detid.rawId(), TrapezoidalParam,
451     Pos.x(), Pos.y(), Pos.z(),
452     WidthVector.x() -Pos.x(), WidthVector.y() -Pos.y(), WidthVector.z() -Pos.z(),
453     LengthVector.x()-Pos.x(), LengthVector.y()-Pos.y(), LengthVector.z()-Pos.z(),
454     ThickVector.x() -Pos.x(), ThickVector.y() -Pos.y(), ThickVector.z() -Pos.z());
455 querten 1.7
456 querten 1.22 Frog_DetId = DETID_MUON;
457     FEB = muon;
458 querten 1.7
459 querten 1.11 Frog_DetId += SubDet*100000;
460     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
461     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
462     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
463    
464     if(SubDet==1){
465 querten 1.22 DTChamberId detid = DTChamberId(Detid);
466 querten 1.11
467 querten 1.22 Frog_DetId += (detid.wheel()+3)*10000;
468 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
469     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
470     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
471    
472 querten 1.22 Frog_DetId += detid.station()*100;
473 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
474     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
475     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
476    
477 querten 1.22 Frog_DetId += detid.sector()*1;
478 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
479     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
480     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
481     }else if(SubDet==2){
482 querten 1.22 CSCDetId detid = CSCDetId(Detid);
483    
484     Frog_DetId += detid.endcap()*10000;
485 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
486     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
487     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
488 querten 1.22
489     Frog_DetId += detid.station()*1000;
490 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
491     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
492     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
493 querten 1.22
494     Frog_DetId += detid.ring()*100;
495 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
496     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
497     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
498 querten 1.22
499     Frog_DetId += detid.chamber()*1;
500 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
501     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
502     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
503 querten 1.22 }else if(SubDet==3){
504     RPCDetId detid = RPCDetId(Detid);
505 querten 1.11
506 querten 1.22 Frog_DetId += (detid.region()+2)*10000;
507 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
508     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
509     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
510    
511 querten 1.22 // Barrel
512     if(detid.region()==0){
513     Frog_DetId += (detid.ring()+3)*1000;
514     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
515     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
516     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
517     // Endcap
518     }else{
519     Frog_DetId += detid.ring()*1000;
520     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
521     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
522     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
523     }
524    
525     Frog_DetId += detid.station()*100;
526 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
527     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
528     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
529    
530 querten 1.22 Frog_DetId += detid.sector()*1;
531 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
532     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
533     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
534     }
535    
536 querten 1.22 FEB->addDaughter(mod);
537     }
538    
539    
540    
541     // ### CALO GEOMETRY ###
542    
543     edm::ESHandle<CaloGeometry> CaloGeom;
544     iSetup.get<IdealGeometryRecord>().get( CaloGeom );
545     const vector<DetId> CaloDets = CaloGeom->getValidDetIds();
546 querten 1.11
547 querten 1.22 for(unsigned int i=0;i<CaloDets.size();i++)
548     {
549     Detid = CaloDets[i];
550     SubDet = Detid.subdetId();
551    
552     if(Detid.det()==DetId::Ecal){
553     if(SubDet<1 || SubDet>3) continue;
554 querten 1.11
555 querten 1.22 Frog_DetId = DETID_ECAL;
556     FEB = ecal;
557 querten 1.11
558 querten 1.22 Frog_DetId += SubDet*100000;
559 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
560     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
561     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
562    
563 querten 1.22 if(SubDet==1){
564     EBDetId detid = EBDetId(Detid);
565    
566     Frog_DetId += ((detid.zside()==1)+1)*10000;
567     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
568     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
569     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
570    
571     Frog_DetId += (detid.tower_ieta()*detid.zside())*100;
572     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
573     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
574     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
575    
576     Frog_DetId += detid.tower_iphi()*1;
577     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
578     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
579     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
580     }else if(SubDet==2){
581     EEDetId detid = EEDetId(Detid);
582    
583     Frog_DetId += ((detid.zside()==1)+1)*10000;
584     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
585     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
586     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
587    
588     Frog_DetId += detid.iquadrant()*1000;
589     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
590     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
591     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
592    
593     Frog_DetId += detid.isc();
594     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
595     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
596     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
597     }else if(SubDet==3){
598     ESDetId detid = ESDetId(Detid);
599    
600     Frog_DetId += ((detid.zside()==1)+1)*10000;
601     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
602     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
603     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
604    
605     Frog_DetId += detid.plane()*1000;
606     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
607     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
608     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
609    
610     Frog_DetId += detid.six()*50;
611     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
612     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
613     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
614    
615     Frog_DetId += detid.siy();
616     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
617     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
618     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
619     }
620     }else if(Detid.det()==DetId::Hcal){
621     if(SubDet<1 || SubDet>4) continue;
622    
623     Frog_DetId = DETID_HCAL;
624     FEB = hcal;
625    
626     Frog_DetId += SubDet*100000;
627 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
628     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
629     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
630    
631 querten 1.22 if(true){
632     HcalDetId detid = HcalDetId(Detid);
633    
634     Frog_DetId += ((detid.zside()==1)+1)*10000;
635     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
636     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
637     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
638    
639     Frog_DetId += detid.ietaAbs()*100;
640     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
641     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
642     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
643    
644     Frog_DetId += detid.iphi()*1;
645     FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
646     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
647     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
648     }
649    
650     }else if(Detid.det() == DetId::Calo){
651     continue;
652     /*
653     Frog_DetId = DETID_CALO + 1000000;
654     FEB = hcal;
655    
656     Frog_DetId += SubDet*100000;
657 querten 1.11 FEB_Temp = dynamic_cast<FROG_Element_Base_With_DetId*>(FEB);
658     FEB = FEB_Temp->getDaughterWithDetId(Frog_DetId);
659     if(FEB==NULL){FEB = new FROG_Element_Base_With_DetId(C_FEB_DETID, Frog_DetId); FEB_Temp->addDaughter(FEB);}
660 querten 1.22 */
661 querten 1.11 }
662    
663 querten 1.22 const CaloCellGeometry* CellGeom = CaloGeom->getGeometry(Detid);
664     GlobalPoint CellPos = CellGeom->getPosition();
665     const CaloCellGeometry::CornersVec CellCorners = CellGeom->getCorners();
666    
667     float cX = (CellCorners[0].x() + CellCorners[2].x())/2;
668     float cY = (CellCorners[0].y() + CellCorners[2].y())/2;
669     float cZ = (CellCorners[0].z() + CellCorners[2].z())/2;
670    
671     float wX = (CellCorners[1].x() - CellCorners[0].x())/2;
672     float wY = (CellCorners[1].y() - CellCorners[0].y())/2;
673     float wZ = (CellCorners[1].z() - CellCorners[0].z())/2;
674    
675     float hX = (CellCorners[3].x() - CellCorners[0].x())/2;
676     float hY = (CellCorners[3].y() - CellCorners[0].y())/2;
677     float hZ = (CellCorners[3].z() - CellCorners[0].z())/2;
678    
679     float F = sqrt( pow(CellCorners[4].x()+CellCorners[6].x(),2) + pow(CellCorners[4].y()+CellCorners[6].y(),2) + pow(CellCorners[4].z()+CellCorners[6].z(),2) );
680     F /= sqrt( pow(CellCorners[0].x()+CellCorners[2].x(),2) + pow(CellCorners[0].y()+CellCorners[2].y(),2) + pow(CellCorners[0].z()+CellCorners[2].z(),2) );
681    
682     FROG_Element_Geom_CaloMod* mod = new FROG_Element_Geom_CaloMod(
683     Detid.rawId(), F,
684     cX , cY , cZ ,
685     wX , wY , wZ ,
686     hX , hY , hZ );
687 querten 1.11
688 querten 1.22 FEB->addDaughter(mod);
689 querten 1.11 }
690    
691 querten 1.22 // ### Save .geom ###
692 querten 1.11
693 querten 1.22 FROG_Geometry* CMS_Geom = new FROG_Geometry(prim);
694     CMS_Geom->Save("CMS.geom");
695 querten 1.11 }
696    
697 querten 1.15 events_prim = new FROG_Element_Base(C_PRIMARY);
698 querten 1.1 }
699    
700     // ------------ method called once each job just after ending the event loop ------------
701     void
702 roberfro 1.13 Frog_Analyzer::endJob()
703 querten 1.22 {
704     if(NEvents>0){
705     char OutputFileName[1024];
706     sprintf(OutputFileName,"%s",OutputFile.substr(0,OutputFile.size()-4).c_str());
707     if(NEventsInVisFile>=0)sprintf(OutputFileName,"%s_R%i_E%i.vis",OutputFileName,SRun,SEvent);
708     sprintf(OutputFileName,"%s.vis",OutputFileName);
709    
710     FROG_Events* events = new FROG_Events(events_prim);
711     events->Save((char*) OutputFileName);
712     delete events;
713     delete events_prim;
714    
715     events_prim = new FROG_Element_Base(C_PRIMARY);
716     NEvents=0;
717     }
718 querten 1.1 }
719    
720    
721    
722     //
723     // member functions
724     //
725    
726     // ------------ method called to for each event ------------
727     void
728     Frog_Analyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
729     {
730 querten 1.22 if(NEvents == 0){
731     SRun = iEvent.id().run();
732     SEvent = iEvent.id().event();
733     }
734    
735    
736 querten 1.23 FROG_Element_Event* event = new FROG_Element_Event(iEvent.id().run(), iEvent.id().event() );
737 querten 1.15 events_prim->addDaughter(event);
738 querten 1.1
739 querten 1.2 // access the tracker
740     edm::ESHandle<TrackerGeometry> tkGeom;
741     iSetup.get<TrackerDigiGeometryRecord>().get( tkGeom );
742 querten 1.1
743 querten 1.2 edm::ESHandle<DTGeometry> DtGeom;
744     iSetup.get<MuonGeometryRecord>().get( DtGeom );
745    
746     edm::ESHandle<CSCGeometry> CscGeom;
747     iSetup.get<MuonGeometryRecord>().get( CscGeom );
748 querten 1.1
749 querten 1.2 edm::ESHandle<RPCGeometry> RpcGeom;
750     iSetup.get<MuonGeometryRecord>().get( RpcGeom );
751 querten 1.1
752 querten 1.2 // ### SIMTRACK ###
753 querten 1.1 for(unsigned int i=0;i<SimTrackProducers.size();i++){
754     edm::Handle<std::vector< SimTrack > > h_SimTracks;
755     iEvent.getByLabel(SimTrackProducers[i], h_SimTracks);
756     std::vector< SimTrack > SimTrackColl = *h_SimTracks.product();
757    
758 querten 1.15 FROG_Element_Base_With_DetId* frogSimTrkColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_SIMTRK + i);
759     event->addDaughter(frogSimTrkColl);
760    
761 querten 1.1 for ( unsigned int a = 0; a < SimTrackColl.size(); ++a ) {
762     SimTrack simTrack =SimTrackColl[a];
763    
764 roberfro 1.13 FROG_Element_Event_Sim_Track* frogSimTrk = new FROG_Element_Event_Sim_Track(simTrack.trackId(),
765     simTrack.momentum().x(), simTrack.momentum().y(), simTrack.momentum().z(),
766     simTrack.momentum().e(), simTrack.type(), simTrack.vertIndex(),
767     simTrack.charge() );
768 querten 1.15 frogSimTrkColl->addDaughter(frogSimTrk);
769 querten 1.1 }
770     }
771    
772 querten 1.2 // ### SIMVERTEX ###
773 querten 1.1 for(unsigned int i=0;i<SimVertexProducers.size();i++){
774     edm::Handle<std::vector< SimVertex > > h_Vertex;
775     iEvent.getByLabel(SimVertexProducers[i], h_Vertex);
776     std::vector< SimVertex > VertexColl = *h_Vertex.product();
777    
778 querten 1.15 FROG_Element_Base_With_DetId* frogSimVtxColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_SIMVTX + i);
779     event->addDaughter(frogSimVtxColl);
780    
781 querten 1.1 for (unsigned int b = 0; b < VertexColl.size(); ++b ) {
782     SimVertex Vertex = VertexColl[b];
783    
784 roberfro 1.13 FROG_Element_Event_Sim_Vertex* frogSimVtx = new FROG_Element_Event_Sim_Vertex(
785     Vertex.position().x(), Vertex.position().y(), Vertex.position().z(),
786     Vertex.parentIndex());
787 querten 1.1
788 querten 1.15 frogSimVtxColl->addDaughter(frogSimVtx);
789 querten 1.1 }
790     }
791    
792 roberfro 1.13
793 querten 1.2 // ### SIMHIT ###
794 querten 1.15
795 querten 1.1 for(unsigned int i=0;i<SimHitProducers.size();i++){
796     edm::Handle<std::vector< PSimHit > > h_Hits;
797     iEvent.getByLabel(SimHitProducers[i], h_Hits);
798     std::vector< PSimHit > Hits = *h_Hits.product();
799    
800 querten 1.15 FROG_Element_Base_With_DetId* frogSimHitColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_SIMHIT + i);
801     event->addDaughter(frogSimHitColl);
802    
803 querten 1.1 for(unsigned int h=0; h<Hits.size(); h++)
804     {
805 querten 1.2 DetId theDetUnitId(Hits[h].detUnitId());
806     const GeomDet* theDet = tkGeom->idToDet(theDetUnitId);
807 querten 1.1
808 querten 1.15
809     FROG_Element_Event_Sim_Hit* frogSimHit = new FROG_Element_Event_Sim_Hit(
810     theDet->surface().toGlobal(Hits[h].localPosition()).x(),theDet->surface().toGlobal(Hits[h].localPosition()).y(),theDet->surface().toGlobal(Hits[h].localPosition()).z(),
811     Hits[h].energyLoss(), Hits[h].processType() );
812     frogSimHitColl->addDaughter(frogSimHit);
813 querten 1.1 }
814     }
815 querten 1.15
816 querten 1.1
817    
818 querten 1.9 // ### TRAJECTORY ###
819 roberfro 1.5 for(unsigned int i=0;i<TrajectoryProducers.size();i++){
820 querten 1.6 Handle<TrajTrackAssociationCollection> trajTrackAssociationHandle;
821     iEvent.getByLabel(TrajectoryProducers[i], trajTrackAssociationHandle);
822     const TrajTrackAssociationCollection TrajToTrackMap = *trajTrackAssociationHandle.product();
823    
824 querten 1.15 FROG_Element_Base_With_DetId* frogTrkColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_TRK + i);
825     event->addDaughter(frogTrkColl);
826    
827 querten 1.6 for(TrajTrackAssociationCollection::const_iterator it = TrajToTrackMap.begin(); it!=TrajToTrackMap.end(); ++it) {
828     Trajectory recoTraj = *it->key;
829 querten 1.15 FROG_Element_Event_Track* frogTrk = new FROG_Element_Event_Track(i,it->val->p(),it->val->pt(),it->key->chiSquared());
830     frogTrkColl->addDaughter(frogTrk);
831 querten 1.1
832 roberfro 1.5 std::vector<TrajectoryMeasurement> measurements = recoTraj.measurements();
833     for(unsigned int h=0;h<measurements.size();h++){
834    
835     TrajectoryMeasurement::ConstRecHitPointer h_it = measurements[h].recHit();
836 querten 1.1 if(!h_it->isValid() )continue;
837     DetId detId = h_it->geographicalId();
838 roberfro 1.5 GlobalPoint globalPos = measurements[h].updatedState().globalPosition();
839 querten 1.1
840 querten 1.15 FROG_Element_Event_Hit* frogHit = new FROG_Element_Event_Hit(detId.rawId(), globalPos.x(),globalPos.y(), globalPos.z(), -1);
841     frogTrk->addDaughter(frogHit);
842 querten 1.1 }
843     }
844     }
845    
846 querten 1.9 // ### TRACK ###
847    
848     for(unsigned int i=0;i<TrackProducers.size();i++){
849     edm::Handle<std::vector< reco::Track > > h_Tracks;
850     iEvent.getByLabel(TrackProducers[i], h_Tracks);
851     std::vector< reco::Track > TrackColl = *h_Tracks.product();
852    
853 querten 1.15 FROG_Element_Base_With_DetId* frogTrkColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_TRK + TrajectoryProducers.size()+i);
854     event->addDaughter(frogTrkColl);
855    
856 querten 1.11 //printf("Size of the TrackCOllection : %i = %i\n",i,TrackColl.size());
857 querten 1.9 for ( unsigned int t = 0; t < TrackColl.size(); ++t ) {
858 querten 1.15 reco::Track recoTrack =TrackColl[t];
859 querten 1.9
860 querten 1.15 FROG_Element_Event_Track* frogTrk = new FROG_Element_Event_Track(TrajectoryProducers.size()+i,recoTrack.p(),recoTrack.pt(), recoTrack.chi2());
861     frogTrkColl->addDaughter(frogTrk);
862 querten 1.9
863 querten 1.15 for(unsigned int h=0;h<recoTrack.recHitsSize();h++){
864     TrackingRecHitRef h_it = recoTrack.recHit(h);
865     if(!h_it->isValid() )continue;
866     DetId detId = h_it->geographicalId();
867     const GeomDet * theDet = NULL;
868    
869     if(detId.det()==DetId::Tracker) theDet = tkGeom->idToDet(detId);
870     if(detId.det()==DetId::Muon && detId.subdetId() == 1) theDet = DtGeom->idToDet(detId);
871     if(detId.det()==DetId::Muon && detId.subdetId() == 2) theDet = CscGeom->idToDet(detId);
872     if(detId.det()==DetId::Muon && detId.subdetId() == 3) theDet = RpcGeom->idToDet(detId);
873    
874     LocalPoint localPos = h_it->localPosition();
875    
876     FROG_Element_Event_Hit* frogHit = new FROG_Element_Event_Hit(
877     detId.rawId(),
878     theDet->surface().toGlobal(localPos).x(),
879     theDet->surface().toGlobal(localPos).y(),
880     theDet->surface().toGlobal(localPos).z(),
881     -1);
882     frogTrk->addDaughter(frogHit);
883     }
884 querten 1.9 }
885     }
886    
887    
888 querten 1.1
889 querten 1.2 // ### ECALRecHits ###
890 querten 1.1 for(unsigned int i=0;i<EcalRecHitProducers.size();i++){
891     edm::Handle<EcalRecHitCollection > h_Ecal_RecHits;
892     iEvent.getByLabel(EcalRecHitProducers[i], h_Ecal_RecHits);
893     EcalRecHitCollection Ecal_RecHits = *h_Ecal_RecHits.product();
894    
895 querten 1.15 FROG_Element_Base_With_DetId* frogHitColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_ECALHIT + i);
896     event->addDaughter(frogHitColl);
897    
898 querten 1.1 for(unsigned int eh=0;eh<Ecal_RecHits.size();eh++){
899 roberfro 1.14 FROG_Element_Event_CaloHit* frogCaloHit = new FROG_Element_Event_CaloHit((Ecal_RecHits[eh].detid()).rawId(),
900     Ecal_RecHits[eh].energy(), Ecal_RecHits[eh].time());
901 querten 1.15 frogHitColl->addDaughter(frogCaloHit);
902 querten 1.1 }
903     }
904    
905    
906 querten 1.2 // ### HCALRecHits ###
907 querten 1.1 for(unsigned int i=0;i<HcalHBHERecHitProducers.size();i++){
908     edm::Handle<HBHERecHitCollection > h_HcalHBHE_RecHits;
909     iEvent.getByLabel(HcalHBHERecHitProducers[i], h_HcalHBHE_RecHits);
910     HBHERecHitCollection HcalHBHE_RecHits = *h_HcalHBHE_RecHits.product();
911    
912 querten 1.15 FROG_Element_Base_With_DetId* frogHitColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_HCALHIT + i);
913     event->addDaughter(frogHitColl);
914    
915 querten 1.1 for(unsigned int hh=0;hh<HcalHBHE_RecHits.size();hh++){
916 roberfro 1.14 FROG_Element_Event_CaloHit* frogCaloHit = new FROG_Element_Event_CaloHit((HcalHBHE_RecHits[hh].detid()).rawId(),
917     HcalHBHE_RecHits[hh].energy(), HcalHBHE_RecHits[hh].time());
918 querten 1.15 frogHitColl->addDaughter(frogCaloHit);
919 querten 1.1 }
920     }
921    
922     for(unsigned int i=0;i<HcalHORecHitProducers.size();i++){
923     edm::Handle<HORecHitCollection > h_HcalHO_RecHits;
924     iEvent.getByLabel(HcalHORecHitProducers[i], h_HcalHO_RecHits);
925     HORecHitCollection HcalHO_RecHits = *h_HcalHO_RecHits.product();
926    
927 querten 1.15 FROG_Element_Base_With_DetId* frogHitColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_HCALHIT + HcalHBHERecHitProducers.size() + i);
928     event->addDaughter(frogHitColl);
929    
930 querten 1.1 for(unsigned int hh=0;hh<HcalHO_RecHits.size();hh++){
931 querten 1.15 FROG_Element_Event_CaloHit* frogCaloHit = new FROG_Element_Event_CaloHit((HcalHO_RecHits[hh].detid()).rawId(),
932     HcalHO_RecHits[hh].energy(), HcalHO_RecHits[hh].time());
933     frogHitColl->addDaughter(frogCaloHit);
934 querten 1.1 }
935     }
936    
937     for(unsigned int i=0;i<HcalHFRecHitProducers.size();i++){
938     edm::Handle<HFRecHitCollection > h_HcalHF_RecHits;
939     iEvent.getByLabel(HcalHFRecHitProducers[i], h_HcalHF_RecHits);
940     HFRecHitCollection HcalHF_RecHits = *h_HcalHF_RecHits.product();
941    
942 querten 1.15 FROG_Element_Base_With_DetId* frogHitColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_HCALHIT + HcalHBHERecHitProducers.size() + HcalHORecHitProducers.size() + i);
943     event->addDaughter(frogHitColl);
944    
945 querten 1.1 for(unsigned int hh=0;hh<HcalHF_RecHits.size();hh++){
946 querten 1.15 FROG_Element_Event_CaloHit* frogCaloHit = new FROG_Element_Event_CaloHit((HcalHF_RecHits[hh].detid()).rawId(),
947     HcalHF_RecHits[hh].energy(), HcalHF_RecHits[hh].time());
948     frogHitColl->addDaughter(frogCaloHit);
949 querten 1.1 }
950     }
951    
952 querten 1.2 // ### Muon Segments ###
953     for(unsigned int i=0;i<CSCSegmentProducers.size();i++){
954     edm::Handle<CSCSegmentCollection > h_CSC_Segments;
955     iEvent.getByLabel(CSCSegmentProducers[i], h_CSC_Segments);
956     CSCSegmentCollection CSC_Segments = *h_CSC_Segments.product();
957    
958 querten 1.15 FROG_Element_Base_With_DetId* frogSegColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_MUONSEG + i);
959     event->addDaughter(frogSegColl);
960    
961 querten 1.2 for(unsigned int s=0;s<CSC_Segments.size();s++){
962     DetId theDetUnitId = CSC_Segments[s].geographicalId();
963     const GeomDet* theDet = CscGeom->idToDet(theDetUnitId);
964    
965 querten 1.15 FROG_Element_Event_Segment* frogSeg = new FROG_Element_Event_Segment(theDetUnitId.rawId(),
966     theDet->surface().toGlobal(CSC_Segments[s].localPosition()).x() , theDet->surface().toGlobal(CSC_Segments[s].localPosition()).y() , theDet->surface().toGlobal(CSC_Segments[s].localPosition()).z() ,
967     theDet->surface().toGlobal(CSC_Segments[s].localDirection()).x(), theDet->surface().toGlobal(CSC_Segments[s].localDirection()).y(), theDet->surface().toGlobal(CSC_Segments[s].localDirection()).z() );
968     frogSegColl->addDaughter(frogSeg);
969 querten 1.2 }
970     }
971    
972     for(unsigned int i=0;i<DTSegmentProducers.size();i++){
973     edm::Handle<DTRecSegment4DCollection > h_DT_Segments;
974     iEvent.getByLabel(DTSegmentProducers[i], h_DT_Segments);
975     DTRecSegment4DCollection DT_Segments = *h_DT_Segments.product();
976    
977 querten 1.15 FROG_Element_Base_With_DetId* frogSegColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_MUONSEG + CSCSegmentProducers.size() + i);
978     event->addDaughter(frogSegColl);
979    
980 querten 1.2 for(unsigned int s=0;s<DT_Segments.size();s++){
981     DetId theDetUnitId = DT_Segments[s].geographicalId();
982     const GeomDet* theDet = DtGeom->idToDet(theDetUnitId);
983    
984 querten 1.15 FROG_Element_Event_Segment* frogSeg = new FROG_Element_Event_Segment(theDetUnitId.rawId(),
985     theDet->surface().toGlobal(DT_Segments[s].localPosition()).x() , theDet->surface().toGlobal(DT_Segments[s].localPosition()).y() , theDet->surface().toGlobal(DT_Segments[s].localPosition()).z() ,
986     theDet->surface().toGlobal(DT_Segments[s].localDirection()).x(), theDet->surface().toGlobal(DT_Segments[s].localDirection()).y(), theDet->surface().toGlobal(DT_Segments[s].localDirection()).z() );
987     frogSegColl->addDaughter(frogSeg);
988 querten 1.2 }
989     }
990    
991     // ### Muon Hits ###
992     for(unsigned int i=0;i<RPCHitsProducers.size();i++){
993     edm::Handle<RPCRecHitCollection > h_RPC_Hits;
994     iEvent.getByLabel(RPCHitsProducers[i], h_RPC_Hits);
995     RPCRecHitCollection RPC_Hits = *h_RPC_Hits.product();
996    
997 querten 1.15 FROG_Element_Base_With_DetId* frogHitColl = new FROG_Element_Base_With_DetId(C_FEB_DETID, EVTID_MUONHIT + i);
998     event->addDaughter(frogHitColl);
999    
1000 querten 1.2 for(unsigned int h=0;h<RPC_Hits.size();h++){
1001     DetId theDetUnitId = RPC_Hits[h].geographicalId();
1002     const GeomDet* theDet = RpcGeom->idToDet(theDetUnitId);
1003    
1004 querten 1.15 FROG_Element_Event_Hit* frogHit = new FROG_Element_Event_Hit(theDetUnitId.rawId(),
1005     theDet->surface().toGlobal(RPC_Hits[h].localPosition()).x() , theDet->surface().toGlobal(RPC_Hits[h].localPosition()).y() , theDet->surface().toGlobal(RPC_Hits[h].localPosition()).z(), -1 );
1006     frogHitColl->addDaughter(frogHit);
1007 querten 1.2 }
1008     }
1009    
1010 querten 1.22 NEvents++;
1011     if(NEvents == NEventsInVisFile){
1012     char OutputFileName[1024];
1013     sprintf(OutputFileName,"%s",OutputFile.substr(0,OutputFile.size()-4).c_str());
1014     if(NEventsInVisFile>=0)sprintf(OutputFileName,"%s_R%i_E%i.vis",OutputFileName,SRun,SEvent);
1015     sprintf(OutputFileName,"%s.vis",OutputFileName);
1016    
1017     FROG_Events* events = new FROG_Events(events_prim);
1018     events->Save((char*) OutputFileName);
1019     delete events;
1020     delete events_prim;
1021 roberfro 1.13
1022 querten 1.22 events_prim = new FROG_Element_Base(C_PRIMARY);
1023     NEvents=0;
1024     }
1025 querten 1.15
1026 querten 1.1 }
1027    
1028    
1029     //define this as a plug-in
1030     DEFINE_FWK_MODULE(Frog_Analyzer);
1031    
1032