Changeset 228 in svn
- Timestamp:
- Feb 2, 2009, 12:44:43 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Delphes.cpp
r212 r228 15 15 #include "TApplication.h" 16 16 #include "TStopwatch.h" 17 18 #include "Utilities/ExRootAnalysis/interface/ExRootTreeReader.h" 19 #include "Utilities/ExRootAnalysis/interface/ExRootTreeWriter.h" 20 #include "Utilities/ExRootAnalysis/interface/ExRootTreeBranch.h" 21 22 #include "interface/DataConverter.h" 23 #include "interface/HEPEVTConverter.h" 24 #include "interface/LHEFConverter.h" 25 #include "interface/STDHEPConverter.h" 26 27 #include "interface/SmearUtil.h" 28 #include "interface/BFieldProp.h" 29 #include "interface/TriggerUtil.h" 30 #include "interface/VeryForward.h" 31 #include "interface/JetUtils.h" 32 #include "interface/FrogUtil.h" 17 #include "TFile.h" 18 19 #include "ExRootTreeReader.h" 20 #include "ExRootTreeWriter.h" 21 #include "ExRootTreeBranch.h" 22 23 #include "DataConverter.h" 24 #include "HEPEVTConverter.h" 25 #include "LHEFConverter.h" 26 #include "STDHEPConverter.h" 27 28 #include "SmearUtil.h" 29 #include "BFieldProp.h" 30 #include "TriggerUtil.h" 31 #include "VeryForward.h" 32 #include "JetsUtil.h" 33 #include "FrogUtil.h" 33 34 34 35 #include <vector> … … 54 55 { 55 56 int appargc = 2; 56 char *appName = "Delphes"; 57 char *appargv[] = {appName, "-b"}; 57 char *appName= new char[20]; 58 char *appOpt= new char[20]; 59 sprintf(appName,"Delphes"); 60 sprintf(appOpt,"-b"); 61 char *appargv[] = {appName,appOpt}; 58 62 TApplication app(appName, &appargc, appargv); 59 63 delete [] appName; 64 delete [] appOpt; 65 66 60 67 if(argc != 4 && argc != 3 && argc != 5) { 61 68 cout << " Usage: " << argv[0] << " input_file output_file [detector_card] [trigger_card] " << endl; … … 66 73 exit(1); 67 74 } 75 76 77 68 78 69 79 cout << endl << endl; … … 163 173 << right << setw(7) <<"**"<<""<<endl; 164 174 cout <<"** **"<< endl; 165 175 166 176 //read the trigger input file 167 177 string TrigDatacard("data/trigger.dat"); … … 182 192 183 193 //Propagation of tracks in the B field 184 TrackPropagation *TRACP = new TrackPropagation(DetDatacard); 194 TrackPropagation *TRACP = new TrackPropagation(DET); 195 //TrackPropagation *TRACP = new TrackPropagation(DetDatacard); 185 196 186 197 //Jet information 187 JetsUtil *JETRUN = new JetsUtil(DetDatacard); 198 JetsUtil *JETRUN = new JetsUtil(DET); 199 //JetsUtil *JETRUN = new JetsUtil(DetDatacard); 188 200 189 201 //VFD information 190 VeryForward * VFD = new VeryForward(DetDatacard); 202 VeryForward * VFD = new VeryForward(DET); 203 //VeryForward * VFD = new VeryForward(DetDatacard); 191 204 192 205 // data converters 193 DataConverter *converter= 0;206 DataConverter *converter=NULL; 194 207 cout <<"** **"<<endl; 195 208 cout <<"** ####### Start convertion to TRoot format ######## **"<< endl; … … 285 298 treeReader->ReadEntry(entry); 286 299 treeWriter->Clear(); 287 if((entry % 100) == 0 && entry > 0 ) {300 if((entry % 100) == 0 && entry > 0 ) { 288 301 cout << left << setw(52) <<"** Processing element # "<<"" 289 302 << left << setw(15) << entry <<"" … … 516 529 treeWriterT->Write(); 517 530 delete treeWriterT; 531 delete treeReaderT; 518 532 } // trigger 519 533 triggerwatch.Stop(); … … 522 536 // 3.2 ************** FROG display 523 537 frogwatch.Start(); 524 if(DET->FLAG_frog == 1) 525 { 538 if(DET->FLAG_frog == 1) { 526 539 cout <<"** **"<<endl; 527 540 cout <<"** ################## Start FROG ################# **"<< endl; 528 541 529 FrogDisplay *FROG = new FrogDisplay(); 530 FROG->BuidEvents(outputfilename,DET->NEvents_Frog); 531 FROG->BuildGeom(DetDatacard); 532 } 533 frogwatch.Stop(); 542 FrogDisplay *FROG = new FrogDisplay(DET); 543 FROG->BuildEvents(outputfilename); 544 FROG->BuildGeom(); 545 delete FROG; 546 } 547 frogwatch.Stop(); 534 548 535 549 … … 591 605 delete JETRUN; 592 606 delete VFD; 593 if(converter)delete converter;607 delete converter; 594 608 595 609 // todo("TODO");
Note:
See TracChangeset
for help on using the changeset viewer.