Changeset 313 in svn for trunk/src/LHCOConverter.cc

Timestamp:

Mar 10, 2009, 12:57:20 PM (15 years ago)

Author:

severine ovyn

Message:

OK etmis and plus joli

File:

• trunk/src/LHCOConverter.cc (modified) (6 diffs)

trunk/src/LHCOConverter.cc

@@ -63 +63 @@
         const unsigned int r_find = outputfilename_.rfind(".root");
         if(r_find > outputfilename_.size()) {
-           cout << "ERROR: the extension of the input file is not '.root'. Exiting...\n"; valid_=false; 
+           cerr << "ERROR: the extension of the input file is not '.root'. Exiting...\n"; valid_=false; 
         }
         else {
@@ -69 +69 @@
                 ofstream outfile( outputfilename_.c_str());
                 outfile.close();
-                cout << "INFO: " << outputfilename_ << " has been created\n";
+                cout << left  << setw(20) <<"**        INFO: "<<""
+                              << left  << setw(29) << outputfilename_ <<""
+                              << right << setw(20) <<"created          **"<<endl;
+
         }
     }
@@ -127 +130 @@
 
   Long64_t Nevents = analysisTree->GetEntries();
-  cout << "** TTree 'Analysis' contains " << Nevents << " events" << endl;
   Nevents = min(Nevents,triggerTree->GetEntries());
   if (Nevents != analysisTree->GetEntries()) 
@@ -155 +157 @@
 
     // 6 MET
-    BranchReader(branchETmis,line,lhcoETmisID);
+    BranchReaderETmis(branchETmis,line);
 
   } // event loop
@@ -173 +175 @@
     for (unsigned int i=0; i< branch_size; i++) {
         particle = (TRootParticle*) branch->At(i);
-        double temp =  sqrt(particle->PT*particle->PT + particle->Pz * particle->Pz);
         TLorentzVector pmu;
         pmu.SetPtEtaPhiE(particle->PT, particle->Eta, particle->Phi, particle->E);
@@ -217 +218 @@
     // return outfile;
 }
+
+//ostringstream LHCOConverter::BranchReader(const TClonesArray* branch, unsigned int& line, unsigned short int lhcoID) const {
+void LHCOConverter::BranchReaderETmis(const TClonesArray* branch, unsigned int& line) const {
+    //ostringstream outfile;
+    ofstream outfile( outputfilename_.c_str(),ios_base::app);
+    TRootETmis * particle = 0;
+    particle = (TRootETmis*) branch->At(0);
+    outfile << fixed << setprecision(3)
+            << setw(3) << line++    // line counter
+            << setw(4) << 6 << " "   // object ID in lhco format
+            << setw(7) << 0.000 <<" "
+            << setw(7) << particle->Phi <<" "
+            << setw(7) << particle->ET  <<" "
+            << setw(7) << 0. <<" "       // invariant mass
+            << setw(7) << 0  <<" "       // number of tracks
+            << setw(7) << 0.  <<" "
+            << setw(7) << 0. <<" "      // E_had/E_em
+            << setw(7) << 0.0 << setw(7) << 0.0  // dummy/dummy
+            << endl;
+
+    outfile.close();
+    // return outfile;
+}
+
+