ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/FastOpenGlDisplayer/src/OpenGLDisplayer.cc
(Generate patch)

Comparing UserCode/FastOpenGlDisplayer/src/OpenGLDisplayer.cc (file contents):
Revision 1.9 by querten, Fri Mar 21 15:11:29 2008 UTC vs.
Revision 1.15 by querten, Sat Apr 5 07:40:33 2008 UTC

# Line 54 | Line 54
54   #include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
55   #include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
56  
57 <
57 > #include "DataFormats/CaloRecHit/interface/CaloRecHit.h"
58 > #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
59 > #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
60  
61   #include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
62   #include "SimDataFormats/TrackingHit/interface/PSimHit.h"
# Line 98 | Line 100 | class OpenGLDisplayer : public edm::EDAn
100  
101  
102        MySimEvents* MyEvents;
103 <      Geometry*    Geom;
103 >      Geometry*    Geom_Tracker;
104 >      Geometry*    Geom_ECAL;
105 >      Geometry*    Geom_HCAL;
106  
107 +      std::string              OutputFile;
108 +
109 +      std::vector<InputTag> SimTrackProducers;
110 +      std::vector<InputTag> SimVertexProducers;
111 +      std::vector<InputTag> SimHitProducers;
112 +
113 +      std::vector<InputTag> TrackProducers;
114 +      std::vector<InputTag> EcalRecHitProducers;
115 +      std::vector<InputTag> HcalHBHERecHitProducers;
116 +      std::vector<InputTag> HcalHORecHitProducers;
117 +      std::vector<InputTag> HcalHFRecHitProducers;
118  
104 –      std::vector<std::string> SimHitSubdetectors;
119  
106 –      std::string              OutputFile;
120  
121  
122        // ----------member data ---------------------------
# Line 114 | Line 127 | class OpenGLDisplayer : public edm::EDAn
127   //
128   OpenGLDisplayer::OpenGLDisplayer(const edm::ParameterSet& iConfig)
129   {
130 <    SimHitSubdetectors  = iConfig.getParameter<std::vector<std::string> >("SimHitSubdetectors");
130 >    OutputFile              = iConfig.getParameter<std::string >("OutputFile");
131 >
132 >    SimTrackProducers       = iConfig.getParameter<std::vector<InputTag> >("SimTrackProducers");
133 >    SimVertexProducers      = iConfig.getParameter<std::vector<InputTag> >("SimVertexProducers");
134 >    SimHitProducers         = iConfig.getParameter<std::vector<InputTag> >("SimHitProducers");
135 >
136 >    TrackProducers          = iConfig.getParameter<std::vector<InputTag> >("TrackProducers");
137 >    EcalRecHitProducers     = iConfig.getParameter<std::vector<InputTag> >("EcalRecHitProducers");
138 >    HcalHBHERecHitProducers = iConfig.getParameter<std::vector<InputTag> >("HcalHBHERecHitProducers");
139 >    HcalHORecHitProducers   = iConfig.getParameter<std::vector<InputTag> >("HcalHORecHitProducers");
140 >    HcalHFRecHitProducers   = iConfig.getParameter<std::vector<InputTag> >("HcalHFRecHitProducers");
141  
119 –    OutputFile          = iConfig.getParameter<std::string >("OutputFile");
142   }
143  
144  
# Line 131 | Line 153 | OpenGLDisplayer::beginJob(const edm::Eve
153     DetId  Detid;
154     int    SubDet;
155  
156 <   MyEvents = new MySimEvents();
157 <   Geom     = new Geometry();
156 >   MyEvents     = new MySimEvents();
157 >   Geom_Tracker = new Geometry();
158 >   Geom_ECAL    = new Geometry();
159 >   Geom_HCAL    = new Geometry();
160 >
161  
162     edm::ESHandle<TrackerGeometry> tkGeom;
163     iSetup.get<TrackerDigiGeometryRecord>().get( tkGeom );
# Line 173 | Line 198 | OpenGLDisplayer::beginJob(const edm::Eve
198  
199        GlobalVector Pos =  GlobalVector(DetUnit->position().basicVector());
200  
201 <      Geom->Add_TrackerDet(Detid.rawId(), TrapezoidalParam,
201 >      Geom_Tracker->Add_TrackerDet(Detid.rawId(), TrapezoidalParam,
202             Pos.x(),                  Pos.y(),                  Pos.z(),
203             WidthVector.x() -Pos.x(), WidthVector.y() -Pos.y(), WidthVector.z() -Pos.z(),
204             LengthVector.x()-Pos.x(), LengthVector.y()-Pos.y(), LengthVector.z()-Pos.z(),
# Line 193 | Line 218 | OpenGLDisplayer::beginJob(const edm::Eve
218     {
219        Detid  = CaloDets[i];
220        SubDet = Detid.subdetId();
196 –    
197 –       if(! (Detid.det()==DetId::Ecal||Detid.det()==DetId::Hcal)) continue;
221  
222 <      //EcalBarrel=1, EcalEndcap=2, EcalPreshower=3, EcalTriggerTower=4, EcalLaserPnDiode=5
223 < //      if(SubDet!=EcalBarrel)continue;
222 >      Geometry* Geom_temp = NULL;
223 >      if(Detid.det()==DetId::Ecal){
224 >        Geom_temp = Geom_ECAL;
225 >
226 >        if(SubDet<1 || SubDet>3){
227 >                printf("Don't Save EcalTriggerTower or EcalLaserPnDiode\n");
228 >                continue;
229 >        }
230 >      }else if(Detid.det()==DetId::Hcal){
231 >        Geom_temp = Geom_HCAL;
232 >
233 >        if(SubDet<1 || SubDet>4){
234 >                printf("Don't Save HcalEmpty, HcalTriggerTower, HcalOther\n");
235 >                continue;
236 >        }
237 >      }else{
238 >        continue;
239 >      }
240  
241        const CaloCellGeometry* CellGeom = CaloGeom->getGeometry(Detid);
242        GlobalPoint CellPos =  CellGeom->getPosition();
204 –
243        const CaloCellGeometry::CornersVec CellCorners =  CellGeom->getCorners();
244  
245 <      float Coord[24];
246 <      for(unsigned int i=0;i<8;i++){
247 <         Coord[i*3+0] = CellCorners[i].x();
248 <         Coord[i*3+1] = CellCorners[i].y();
249 <         Coord[i*3+2] = CellCorners[i].z();
250 <      }
251 <
252 <      Geom->Add_CaloDet(Detid.rawId(),
253 <                        CellPos.x(),   CellPos.y(),   CellPos.z(),
254 <                        Coord[0 ]  ,   Coord[1 ]  ,   Coord[2 ]  ,
217 <                        Coord[3 ]  ,   Coord[4 ]  ,   Coord[5 ]  ,
218 <                        Coord[6 ]  ,   Coord[7 ]  ,   Coord[8 ]  ,
219 <                        Coord[9 ]  ,   Coord[10]  ,   Coord[11]  ,
220 <                        Coord[12]  ,   Coord[13]  ,   Coord[14]  ,
221 <                        Coord[15]  ,   Coord[16]  ,   Coord[17]  ,
222 <                        Coord[18]  ,   Coord[19]  ,   Coord[20]  ,
223 <                        Coord[21]  ,   Coord[22]  ,   Coord[23]  );
224 <                        
225 <      printf("Det = %i Subdet = %i --> %6.2f %6.2f %6.2f - %i\n",  Detid.det(), SubDet, CellPos.x(), CellPos.y(), CellPos.z(), CellCorners.size());
245 >      Geom_temp->Add_CaloDet(Detid.rawId(),
246 >                        CellPos.x()         ,   CellPos.y()         ,   CellPos.z(),
247 >                        CellCorners[0].x()  ,   CellCorners[0].y()  ,   CellCorners[0].z()  ,
248 >                        CellCorners[1].x()  ,   CellCorners[1].y()  ,   CellCorners[1].z()  ,
249 >                        CellCorners[2].x()  ,   CellCorners[2].y()  ,   CellCorners[2].z()  ,
250 >                        CellCorners[3].x()  ,   CellCorners[3].y()  ,   CellCorners[3].z()  ,
251 >                        CellCorners[4].x()  ,   CellCorners[4].y()  ,   CellCorners[4].z()  ,
252 >                        CellCorners[5].x()  ,   CellCorners[5].y()  ,   CellCorners[5].z()  ,
253 >                        CellCorners[6].x()  ,   CellCorners[6].y()  ,   CellCorners[6].z()  ,
254 >                        CellCorners[7].x()  ,   CellCorners[7].y()  ,   CellCorners[7].z()  );
255     }
256  
257 +   Geom_Tracker->Save("Tracker.geom");    
258 +   Geom_ECAL->Save("ECAL.geom");
259 +   Geom_HCAL->Save("HCAL.geom");
260  
229 –   Geom->Save("Tracker.geom");    
261  
262  
263   }
# Line 236 | Line 267 | void
267   OpenGLDisplayer::endJob() {
268    
269     MyEvents->Save((char*) OutputFile.c_str());
270 <   MyEvents->Load((char*) OutputFile.c_str());
270 > //   MyEvents->Load((char*) OutputFile.c_str());
271   }
272  
273  
# Line 249 | Line 280 | OpenGLDisplayer::endJob() {
280   void
281   OpenGLDisplayer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
282   {
283 <
253 <        edm::Handle<std::vector< SimTrack > > h_SimTracks;
254 <        iEvent.getByLabel("g4SimHits", h_SimTracks);
255 <        std::vector< SimTrack > SimTrackColl = *h_SimTracks.product();
256 <
257 <        edm::Handle<std::vector< SimVertex > > h_Vertex;
258 <        iEvent.getByLabel("g4SimHits", h_Vertex);
259 <        std::vector< SimVertex > VertexColl = *h_Vertex.product();
283 >        MySimEvent* MyEvent = new MySimEvent;
284  
285  
286  
# Line 269 | Line 293 | OpenGLDisplayer::analyze(const edm::Even
293          TrackingGeometry::DetContainer theDetUnits = theTracker.dets();
294  
295  
296 +        // Save Data
297  
298 +        for(unsigned int i=0;i<SimTrackProducers.size();i++){
299 +           edm::Handle<std::vector< SimTrack > > h_SimTracks;
300 +           iEvent.getByLabel(SimTrackProducers[i], h_SimTracks);
301 +           std::vector< SimTrack > SimTrackColl = *h_SimTracks.product();
302  
303 <        // Save the data
275 <
276 <        MySimEvent* MyEvent = new MySimEvent;
303 >           MySimEvent* MyEvent = new MySimEvent;
304    
305 <        for ( unsigned int a = 0; a < SimTrackColl.size(); ++a ) {
305 >           for ( unsigned int a = 0; a < SimTrackColl.size(); ++a ) {
306                  MySimTrack  MysimTrack;
307                  SimTrack    simTrack    =SimTrackColl[a];
308  
# Line 289 | Line 316 | OpenGLDisplayer::analyze(const edm::Even
316                  MysimTrack.charge       =simTrack.charge();
317  
318                  MyEvent->MySimTrackCollection.push_back(MysimTrack);
319 <        }
319 >           }
320 >        }
321 >
322  
323 +        for(unsigned int i=0;i<SimVertexProducers.size();i++){
324 +           edm::Handle<std::vector< SimVertex > > h_Vertex;
325 +           iEvent.getByLabel(SimVertexProducers[i], h_Vertex);
326 +           std::vector< SimVertex > VertexColl = *h_Vertex.product();
327  
328 <        for (unsigned int b = 0; b < VertexColl.size(); ++b ) {
328 >           for (unsigned int b = 0; b < VertexColl.size(); ++b ) {
329                  MySimVertex     MyVertex;
330                  SimVertex       Vertex  =  VertexColl[b];
331  
# Line 302 | Line 335 | OpenGLDisplayer::analyze(const edm::Even
335                  MyVertex.z              =Vertex.position().z();
336  
337                  MyEvent->MySimVertexCollection.push_back(MyVertex);
338 <        }
338 >           }
339 >        }
340  
341 <        for(unsigned int i=0; i<SimHitSubdetectors.size(); i++)
308 <        {
341 >        for(unsigned int i=0;i<SimHitProducers.size();i++){
342             edm::Handle<std::vector< PSimHit > > h_Hits;
343 <           iEvent.getByLabel("g4SimHits",SimHitSubdetectors[i].c_str(), h_Hits);
343 >           iEvent.getByLabel(SimHitProducers[i], h_Hits);
344             std::vector< PSimHit > Hits = *h_Hits.product();
345  
346             for(unsigned int h=0; h<Hits.size(); h++)
# Line 323 | Line 356 | OpenGLDisplayer::analyze(const edm::Even
356                  Hit.dEdX        = Hits[h].energyLoss();
357  
358                  MyEvent->MyPSimHitCollection.push_back(Hit);
326 –
359             }
360          }
361  
330 –        edm::Handle<std::vector< reco::Track > > h_Tracks;
331 –        iEvent.getByLabel("ctfWithMaterialTracks", h_Tracks);
332 –        std::vector< reco::Track > TrackColl = *h_Tracks.product();
362  
363 <        for ( unsigned int t = 0; t < TrackColl.size(); ++t ) {
363 >        for(unsigned int i=0;i<TrackProducers.size();i++){
364 >           edm::Handle<std::vector< reco::Track > > h_Tracks;
365 >           iEvent.getByLabel(TrackProducers[i], h_Tracks);
366 >           std::vector< reco::Track > TrackColl = *h_Tracks.product();
367 >
368 >           for ( unsigned int t = 0; t < TrackColl.size(); ++t ) {
369                  MyRecoTrack  MyrecoTrack;
370                  reco::Track  recoTrack  =TrackColl[t];
371  
338 – //                MyrecoTrack.Hits = new MyRecoHit[MyrecoTrack.N];
339 –
372                  for(unsigned int h=0;h<recoTrack.recHitsSize();h++){
373                        TrackingRecHitRef h_it     = recoTrack.recHit(h);
374                        if(!h_it->isValid() )continue;
# Line 351 | Line 383 | OpenGLDisplayer::analyze(const edm::Even
383                        hit.DetId  = detId.rawId();
384                        hit.Charge = -1;
385                        MyrecoTrack.Hits.push_back(hit);
354 –
355 –
356 – //                    printf("%8.2f %8.2f %8.2f\n",theDet->surface().toGlobal(localPos).x(),theDet->surface().toGlobal(localPos).y(),theDet->surface().toGlobal(localPos).z());
386                  }
387                  MyEvent->MyRecoTrackCollection.push_back(MyrecoTrack);
388 +           }
389          }
390  
361 –        MyEvents->Events.push_back(MyEvent);  
391  
392 +        for(unsigned int i=0;i<EcalRecHitProducers.size();i++){
393 +           edm::Handle<EcalRecHitCollection > h_Ecal_RecHits;
394 +           iEvent.getByLabel(EcalRecHitProducers[i], h_Ecal_RecHits);
395 +           EcalRecHitCollection Ecal_RecHits = *h_Ecal_RecHits.product();
396 +
397 +           for(unsigned int eh=0;eh<Ecal_RecHits.size();eh++){
398 +                MyCaloHit temp_EcalHit;
399 +                temp_EcalHit.E     = Ecal_RecHits[eh].energy();
400 +                temp_EcalHit.t     = Ecal_RecHits[eh].time();
401 +                temp_EcalHit.DetId = (Ecal_RecHits[eh].detid()).rawId();
402 +        
403 +                MyEvent->MyEcalCaloHitCollection.push_back(temp_EcalHit);
404 +           }
405 +        }
406 +
407 +
408 +        for(unsigned int i=0;i<HcalHBHERecHitProducers.size();i++){
409 +           edm::Handle<HBHERecHitCollection > h_HcalHBHE_RecHits;
410 +           iEvent.getByLabel(HcalHBHERecHitProducers[i], h_HcalHBHE_RecHits);
411 +           HBHERecHitCollection HcalHBHE_RecHits = *h_HcalHBHE_RecHits.product();
412 +
413 +           for(unsigned int hh=0;hh<HcalHBHE_RecHits.size();hh++){
414 +                MyCaloHit temp_HcalHBHEHit;
415 +                temp_HcalHBHEHit.E      = HcalHBHE_RecHits[hh].energy();
416 +                temp_HcalHBHEHit.t      = HcalHBHE_RecHits[hh].time();
417 +                temp_HcalHBHEHit.DetId = (HcalHBHE_RecHits[hh].detid()).rawId();
418 +
419 +                MyEvent->MyHcalCaloHitCollection.push_back(temp_HcalHBHEHit);
420 +           }
421 +        }
422 +
423 +
424 +        for(unsigned int i=0;i<HcalHORecHitProducers.size();i++){
425 +           edm::Handle<HORecHitCollection > h_HcalHO_RecHits;
426 +           iEvent.getByLabel(HcalHORecHitProducers[i], h_HcalHO_RecHits);
427 +           HORecHitCollection HcalHO_RecHits = *h_HcalHO_RecHits.product();
428 +
429 +           for(unsigned int hh=0;hh<HcalHO_RecHits.size();hh++){
430 +                MyCaloHit temp_HcalHOHit;
431 +                temp_HcalHOHit.E      = HcalHO_RecHits[hh].energy();
432 +                temp_HcalHOHit.t      = HcalHO_RecHits[hh].time();
433 +                temp_HcalHOHit.DetId = (HcalHO_RecHits[hh].detid()).rawId();
434 +
435 +                MyEvent->MyHcalCaloHitCollection.push_back(temp_HcalHOHit);
436 +           }
437 +        }
438 +
439 +
440 +        for(unsigned int i=0;i<HcalHFRecHitProducers.size();i++){
441 +           edm::Handle<HFRecHitCollection > h_HcalHF_RecHits;
442 +           iEvent.getByLabel(HcalHFRecHitProducers[i], h_HcalHF_RecHits);
443 +           HFRecHitCollection HcalHF_RecHits = *h_HcalHF_RecHits.product();
444 +
445 +           for(unsigned int hh=0;hh<HcalHF_RecHits.size();hh++){
446 +                MyCaloHit temp_HcalHFHit;
447 +                temp_HcalHFHit.E      = HcalHF_RecHits[hh].energy();
448 +                temp_HcalHFHit.t      = HcalHF_RecHits[hh].time();
449 +                temp_HcalHFHit.DetId = (HcalHF_RecHits[hh].detid()).rawId();
450 +
451 +                MyEvent->MyHcalCaloHitCollection.push_back(temp_HcalHFHit);
452 +           }
453 +        }
454 +
455 +        MyEvents->Events.push_back(MyEvent);  
456   }
457  
458  

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines