TString name_ttEThad = "ttEThad_" + str_ieta + "_" + str_iphi;
TString name_ttETtot = "ttETtot_" + str_ieta + "_" + str_iphi;
_ttETem[ieta][iphi] = new TH1F(name_ttETem, name_ttETem, 48, -2, 10);
_ttEThad[ieta][iphi] = new TH1F(name_ttEThad, name_ttEThad, 48, -2, 10);
_ttETtot[ieta][iphi] = new TH1F(name_ttETtot, name_ttETtot, 48, -2, 10);
}
}
TTDir2->cd("..");
TDirectory *TTDir3 = new TDirectory("TTDir3", "Energy deposited in each of 1280 Trigger Towers");
TTDir3->cd();
for(int i = 0; i<1280 ; i++)
{
char buffer[33];
sprintf(buffer,"%d", i);
TString str = buffer;
TString name_totTTet = "totTTet_" + str;
_totTTet[i] = new TH1F(name_totTTet, name_totTTet, 48, -2, 10);
}
TTDir3->cd("..");
hist_L1JetLead_TTlead = new TH2F("hist_L1JetLead_TTlead", "Leading TT of Leading L1 jet vs leading L1 jet", 100, 0, 500, 100, 0, 150);//
hist_2ndL1JetLead_TTlead = new TH2F("hist_2ndL1JetLead_TTlead", "Leading TT of 2nd Leading L1 jet vs 2nd leading L1 jet", 100, 0, 500, 100, 0, 150);
hist_JetLead_TTmax = new TH2F("hist_JetLead_TTmax", "Leading TT of Leading jet vs leading jet", 100, 0, 500, 100, 0, 150);// hist_2ndJetLead_TTmax = new TH2F("hist_2ndJetLead_TTmax", "Leading TT of Leading jet vs 2nd leading jet", 100, 0, 500, 100, 0, 150);
evtno=0;//
}// TTplots::begin()
////////////////////////////////////////processEvent()////////////////////////////////////////////////////
bool TTplots::processEvent(cafe::Event& event)
{
bool passed = true;
evtno++;
cout << "-->Now beinning of process event " << evtno < //============================== L1 Trigger Towers=====================================
cafe::Collection L1CalEMTowers=event.getL1CalEMTowers();
cafe::Collection L1CalHADTowers=event.getL1CalHADTowers();
cafe::Collection L1CalTotalTowers=event.getL1CalTotalTowers();
//------------------------------map-----------------------------------
TString name_ttETtot = "ttETtot_" + str_ieta + "_" + str_iphi;
_ttETem[ieta][iphi] = new TH1F(name_ttETem, name_ttETem, 48, -2, 10);
_ttEThad[ieta][iphi] = new TH1F(name_ttEThad, name_ttEThad, 48, -2, 10);
_ttETtot[ieta][iphi] = new TH1F(name_ttETtot, name_ttETtot, 48, -2, 10);
}
}
TTDir2->cd("..");
TDirectory *TTDir3 = new TDirectory("TTDir3", "Energy deposited in each of 1280 Trigger Towers");
TTDir3->cd();
for(int i = 0; i<1280 ; i++)
{
char buffer[33];
sprintf(buffer,"%d", i);
TString str = buffer;
TString name_totTTet = "totTTet_" + str;
_totTTet[i] = new TH1F(name_totTTet, name_totTTet, 48, -2, 10);
}
TTDir3->cd("..");
hist_L1JetLead_TTlead = new TH2F("hist_L1JetLead_TTlead", "Leading TT of Leading L1 jet vs leading L1 jet", 100, 0, 500, 100, 0, 150);//
hist_2ndL1JetLead_TTlead = new TH2F("hist_2ndL1JetLead_TTlead", "Leading TT of 2nd Leading L1 jet vs 2nd leading L1 jet", 100, 0, 500, 100, 0, 150);
hist_JetLead_TTmax = new TH2F("hist_JetLead_TTmax", "Leading TT of Leading jet vs leading jet", 100, 0, 500, 100, 0, 150);// hist_2ndJetLead_TTmax = new TH2F("hist_2ndJetLead_TTmax", "Leading TT of Leading jet vs 2nd leading jet", 100, 0, 500, 100, 0, 150);
evtno=0;//
}// TTplots::begin()
////////////////////////////////////////processEvent()////////////////////////////////////////////////////
bool TTplots::processEvent(cafe::Event& event)
{
bool passed = true;
evtno++;
cout << "-->Now beinning of process event " << evtno <
cafe::Collection
cafe::Collection
cafe::Collection
//------------------------------map-----------------------------------
Comments
Post a Comment
https://gengwg.blogspot.com/