void RecoInfoTree::Loop()
{
BookHistograms();
if (fChain == 0) return;
Long64_t nentries = fChain->GetEntriesFast();
Long64_t nbytes = 0, nb = 0;
for (Long64_t jentry=0; jentry<nentries;jentry++) {
Long64_t ientry = LoadTree(jentry);
if (ientry < 0) break;
nb = fChain->GetEntry(jentry); nbytes += nb;
if(jentry%1000==0) std::cout<<" processing event "<<jentry<<std::endl;
for(int i=0; i<(int)Muon_n ; i++){
//if(Muon_n>0){
TLorentzVector* p4 = (TLorentzVector*) Muon_4Momentum->At(i);
if(p4->Pt()<10) continue ;
double correctedIso = ( (*Muon_particleIso)[i] +
TMath::Max(double( (*Muon_neutralIso)[i] + (*Muon_photonIso)[i]-
TMath::Max((0.5*(*Muon_dbeta)[i]),0.0)),0.0 ))/p4->Pt();
Muon_Iso_->Fill(correctedIso);
Muon_Significance3D_->Fill(TMath::Abs((*Muon_Significance3D)[i]));
Muon_NewVz_->Fill((*Muon_NewVz)[i]);
Muon_normChi2_->Fill((*Muon_normChi2)[i]);
Muon_validHits_->Fill((*Muon_validHits)[i]);
Muon_tkHits_->Fill((*Muon_tkHits)[i]);
Muon_pixHits_->Fill((*Muon_pixHits)[i]);
Muon_numberOfMatches_->Fill((*Muon_numberOfMatches)[i]);
}
for(int i=0; i<(int)Electron_n;i++ ){
TLorentzVector* p4 = (TLorentzVector*) Electron_4Momentum->At(i);
if(p4->Pt()<10) continue ;
double correctedIso = ( (*Electron_particleIso)[i] +
TMath::Max(double( (*Electron_neutralIso)[i] + (*Electron_photonIso)[i]-
TMath::Max((0.5*(*Electron_dbeta)[i]),0.0)),0.0 ))/p4->Pt();
Electron_Iso_->Fill(correctedIso);
Electron_MVAID_new_->Fill((*Electron_MVAID_new)[i]);
Electron_HoE_->Fill((*Electron_HoE)[i]);
Electron_SigmaIetaIeta_->Fill((*Electron_SigmaIetaIeta)[i]);
Electron_dEtaIn_->Fill((*Electron_dEtaIn)[i]);
Electron_dPhiIn_->Fill((*Electron_dPhiIn)[i]);
Electron_missing_hits_->Fill((*Electron_missing_hits)[i]);
Electron_Significance3D_->Fill(TMath::Abs((*Electron_Significance3D)[i]));
}
}
WriteHistograms();
}
int main (int argc,char **argv)
{
CheckOpts(argc,argv);
GetFQHN();
if (strlen(ERASE) > 0)
{
EraseAverages();
exit(0);
}
ReadAverages();
SummarizeAverages();
if (strlen(FLNAME) == 0)
{
if (TIMESTAMPS)
{
if ((NOW = time((time_t *)NULL)) == -1)
{
Verbose("Couldn't read system clock\n");
}
sprintf(FLNAME,"cfenvgraphs-%s-%s",CanonifyName(VFQNAME),ctime(&NOW));
}
else
{
sprintf(FLNAME,"cfenvgraphs-snapshot-%s",CanonifyName(VFQNAME));
}
}
Verbose("Creating sub-directory %s\n",FLNAME);
if (mkdir(FLNAME,0755) == -1)
{
Verbose("Writing to existing directory\n");
}
if (chdir(FLNAME))
{
perror("chdir");
exit(0);
}
Verbose("Writing data to sub-directory %s: \n x,y1,y2,y3...\n ",FLNAME);
if (NOWOPT)
{
MagnifyNow();
}
else
{
WriteGraphFiles();
WriteHistograms();
DiskArrivals();
PeerIntermittency();
}
return 0;
}
