#ifndef analysis_h #define analysis_h #include "VisMassCycleTree.h" #include "TString.h" //#include "TEfficiency.h" #include "TH2F.h" #include "TFile.h" #include "TStyle.h" #include "TROOT.h" class analysis : public VisMassCycleTree{ public: analysis(TTree *tree=0, TString inputtag=""):VisMassCycleTree(tree),InputTag(inputtag) { initialize();} virtual ~analysis(){ terminate();} void runAnalysis(TString filename, bool update); private: bool checkZtautauEventSel(); bool checkZtautauEventSelBeforeControl(); bool isInA(); bool isInB(); bool isInC(); bool isInD(); TString InputTag; void write(TString name, bool update); // fillers void FillHistosWithoutTrigger(); void FillHistosPassTrigger_tau12loose(); void FillHistosPassTrigger_tau12loose1p(); void FillHistosPassTrigger_tau12loose3p(); void FillHistosPassTrigger_tau16loose(); void FillHistosPassTrigger_tau16loose1p(); void FillHistosPassTrigger_tau16loose3p(); void FillHistosPassTrigger_tau16medium(); void FillHistosPassTrigger_tau16medium1p(); void FillHistosPassTrigger_tau16medium3p(); // the histograms ------------- // // all taus ------------------- // no selection ---------- TH2F *h2_tau_vmpt_A; TH2F *h2_tau_vmpt_B; TH2F *h2_tau_vmpt_C; TH2F *h2_tau_vmpt_D; // ----------------------- // EF_tau12_loose TH2F *h2_tau_vmpt_EF_tau12_loose_A; TH2F *h2_tau_vmpt_EF_tau12_loose_B; TH2F *h2_tau_vmpt_EF_tau12_loose_C; TH2F *h2_tau_vmpt_EF_tau12_loose_D; // EF_tau16_loose TH2F *h2_tau_vmpt_EF_tau16_loose_A; TH2F *h2_tau_vmpt_EF_tau16_loose_B; TH2F *h2_tau_vmpt_EF_tau16_loose_C; TH2F *h2_tau_vmpt_EF_tau16_loose_D; // EF_tau16_medium TH2F *h2_tau_vmpt_EF_tau16_medium_A; TH2F *h2_tau_vmpt_EF_tau16_medium_B; TH2F *h2_tau_vmpt_EF_tau16_medium_C; TH2F *h2_tau_vmpt_EF_tau16_medium_D; // for future reference // w-control region // // no selection ---------- TH2F *h2_tau1p_vmpt_A; TH2F *h2_tau1p_vmpt_B; TH2F *h2_tau1p_vmpt_C; TH2F *h2_tau1p_vmpt_D; // ----------------------- // EF_tau12_loose TH2F *h2_tau1p_vmpt_EF_tau12_loose_A; TH2F *h2_tau1p_vmpt_EF_tau12_loose_B; TH2F *h2_tau1p_vmpt_EF_tau12_loose_C; TH2F *h2_tau1p_vmpt_EF_tau12_loose_D; // EF_tau16_loose TH2F *h2_tau1p_vmpt_EF_tau16_loose_A; TH2F *h2_tau1p_vmpt_EF_tau16_loose_B; TH2F *h2_tau1p_vmpt_EF_tau16_loose_C; TH2F *h2_tau1p_vmpt_EF_tau16_loose_D; // EF_tau16_medium TH2F *h2_tau1p_vmpt_EF_tau16_medium_A; TH2F *h2_tau1p_vmpt_EF_tau16_medium_B; TH2F *h2_tau1p_vmpt_EF_tau16_medium_C; TH2F *h2_tau1p_vmpt_EF_tau16_medium_D; // for future reference // w-control region // no selection ---------- TH2F *h2_tau3p_vmpt_A; TH2F *h2_tau3p_vmpt_B; TH2F *h2_tau3p_vmpt_C; TH2F *h2_tau3p_vmpt_D; // ----------------------- // EF_tau12_loose TH2F *h2_tau3p_vmpt_EF_tau12_loose_A; TH2F *h2_tau3p_vmpt_EF_tau12_loose_B; TH2F *h2_tau3p_vmpt_EF_tau12_loose_C; TH2F *h2_tau3p_vmpt_EF_tau12_loose_D; // EF_tau16_loose TH2F *h2_tau3p_vmpt_EF_tau16_loose_A; TH2F *h2_tau3p_vmpt_EF_tau16_loose_B; TH2F *h2_tau3p_vmpt_EF_tau16_loose_C; TH2F *h2_tau3p_vmpt_EF_tau16_loose_D; // EF_tau16_medium TH2F *h2_tau3p_vmpt_EF_tau16_medium_A; TH2F *h2_tau3p_vmpt_EF_tau16_medium_B; TH2F *h2_tau3p_vmpt_EF_tau16_medium_C; TH2F *h2_tau3p_vmpt_EF_tau16_medium_D; // for future reference // w-control region void initialize() { int b_vm = 40; int b_pt = 40; double min_vm = 0.0; double max_vm = 200.0; double min_pt = 0.0; double max_pt = 200.0; h2_tau_vmpt_A = new TH2F("h2_tau_vmpt_"+InputTag+"_A","h2_tau_vmpt_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_B = new TH2F("h2_tau_vmpt_"+InputTag+"_B","h2_tau_vmpt_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_C = new TH2F("h2_tau_vmpt_"+InputTag+"_C","h2_tau_vmpt_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_D = new TH2F("h2_tau_vmpt_"+InputTag+"_D","h2_tau_vmpt_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); // EF_tau12_loose h2_tau_vmpt_EF_tau12_loose_A = new TH2F("h2_tau_vmpt_EF_tau12_loose_"+InputTag+"_A","h2_tau_vmpt_EF_tau12_loose_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_EF_tau12_loose_B = new TH2F("h2_tau_vmpt_EF_tau12_loose_"+InputTag+"_B","h2_tau_vmpt_EF_tau12_loose_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_EF_tau12_loose_C = new TH2F("h2_tau_vmpt_EF_tau12_loose_"+InputTag+"_C","h2_tau_vmpt_EF_tau12_loose_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_EF_tau12_loose_D = new TH2F("h2_tau_vmpt_EF_tau12_loose_"+InputTag+"_D","h2_tau_vmpt_EF_tau12_loose_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); // EF_tau16_loose h2_tau_vmpt_EF_tau16_loose_A = new TH2F("h2_tau_vmpt_EF_tau16_loose_"+InputTag+"_A","h2_tau_vmpt_EF_tau16_loose_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_EF_tau16_loose_B = new TH2F("h2_tau_vmpt_EF_tau16_loose_"+InputTag+"_B","h2_tau_vmpt_EF_tau16_loose_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_EF_tau16_loose_C = new TH2F("h2_tau_vmpt_EF_tau16_loose_"+InputTag+"_C","h2_tau_vmpt_EF_tau16_loose_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_EF_tau16_loose_D = new TH2F("h2_tau_vmpt_EF_tau16_loose_"+InputTag+"_D","h2_tau_vmpt_EF_tau16_loose_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); // EF_tau16_medium h2_tau_vmpt_EF_tau16_medium_A = new TH2F("h2_tau_vmpt_EF_tau16_medium_"+InputTag+"_A","h2_tau_vmpt_EF_tau16_medium_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_EF_tau16_medium_B = new TH2F("h2_tau_vmpt_EF_tau16_medium_"+InputTag+"_B","h2_tau_vmpt_EF_tau16_medium_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_EF_tau16_medium_C = new TH2F("h2_tau_vmpt_EF_tau16_medium_"+InputTag+"_C","h2_tau_vmpt_EF_tau16_medium_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau_vmpt_EF_tau16_medium_D = new TH2F("h2_tau_vmpt_EF_tau16_medium_"+InputTag+"_D","h2_tau_vmpt_EF_tau16_medium_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_A = new TH2F("h2_tau1p_vmpt_"+InputTag+"_A","h2_tau1p_vmpt_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_B = new TH2F("h2_tau1p_vmpt_"+InputTag+"_B","h2_tau1p_vmpt_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_C = new TH2F("h2_tau1p_vmpt_"+InputTag+"_C","h2_tau1p_vmpt_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_D = new TH2F("h2_tau1p_vmpt_"+InputTag+"_D","h2_tau1p_vmpt_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); // EF_tau12_loose h2_tau1p_vmpt_EF_tau12_loose_A = new TH2F("h2_tau1p_vmpt_EF_tau12_loose_"+InputTag+"_A","h2_tau1p_vmpt_EF_tau12_loose_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_EF_tau12_loose_B = new TH2F("h2_tau1p_vmpt_EF_tau12_loose_"+InputTag+"_B","h2_tau1p_vmpt_EF_tau12_loose_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_EF_tau12_loose_C = new TH2F("h2_tau1p_vmpt_EF_tau12_loose_"+InputTag+"_C","h2_tau1p_vmpt_EF_tau12_loose_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_EF_tau12_loose_D = new TH2F("h2_tau1p_vmpt_EF_tau12_loose_"+InputTag+"_D","h2_tau1p_vmpt_EF_tau12_loose_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); // EF_tau16_loose h2_tau1p_vmpt_EF_tau16_loose_A = new TH2F("h2_tau1p_vmpt_EF_tau16_loose_"+InputTag+"_A","h2_tau1p_vmpt_EF_tau16_loose_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_EF_tau16_loose_B = new TH2F("h2_tau1p_vmpt_EF_tau16_loose_"+InputTag+"_B","h2_tau1p_vmpt_EF_tau16_loose_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_EF_tau16_loose_C = new TH2F("h2_tau1p_vmpt_EF_tau16_loose_"+InputTag+"_C","h2_tau1p_vmpt_EF_tau16_loose_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_EF_tau16_loose_D = new TH2F("h2_tau1p_vmpt_EF_tau16_loose_"+InputTag+"_D","h2_tau1p_vmpt_EF_tau16_loose_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); // EF_tau16_medium h2_tau1p_vmpt_EF_tau16_medium_A = new TH2F("h2_tau1p_vmpt_EF_tau16_medium_"+InputTag+"_A","h2_tau1p_vmpt_EF_tau16_medium_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_EF_tau16_medium_B = new TH2F("h2_tau1p_vmpt_EF_tau16_medium_"+InputTag+"_B","h2_tau1p_vmpt_EF_tau16_medium_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_EF_tau16_medium_C = new TH2F("h2_tau1p_vmpt_EF_tau16_medium_"+InputTag+"_C","h2_tau1p_vmpt_EF_tau16_medium_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau1p_vmpt_EF_tau16_medium_D = new TH2F("h2_tau1p_vmpt_EF_tau16_medium_"+InputTag+"_D","h2_tau1p_vmpt_EF_tau16_medium_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_A = new TH2F("h2_tau3p_vmpt_"+InputTag+"_A","h2_tau3p_vmpt_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_B = new TH2F("h2_tau3p_vmpt_"+InputTag+"_B","h2_tau3p_vmpt_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_C = new TH2F("h2_tau3p_vmpt_"+InputTag+"_C","h2_tau3p_vmpt_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_D = new TH2F("h2_tau3p_vmpt_"+InputTag+"_D","h2_tau3p_vmpt_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); // EF_tau12_loose h2_tau3p_vmpt_EF_tau12_loose_A = new TH2F("h2_tau3p_vmpt_EF_tau12_loose_"+InputTag+"_A","h2_tau3p_vmpt_EF_tau12_loose_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_EF_tau12_loose_B = new TH2F("h2_tau3p_vmpt_EF_tau12_loose_"+InputTag+"_B","h2_tau3p_vmpt_EF_tau12_loose_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_EF_tau12_loose_C = new TH2F("h2_tau3p_vmpt_EF_tau12_loose_"+InputTag+"_C","h2_tau3p_vmpt_EF_tau12_loose_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_EF_tau12_loose_D = new TH2F("h2_tau3p_vmpt_EF_tau12_loose_"+InputTag+"_D","h2_tau3p_vmpt_EF_tau12_loose_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); // EF_tau16_loose h2_tau3p_vmpt_EF_tau16_loose_A = new TH2F("h2_tau3p_vmpt_EF_tau16_loose_"+InputTag+"_A","h2_tau3p_vmpt_EF_tau16_loose_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_EF_tau16_loose_B = new TH2F("h2_tau3p_vmpt_EF_tau16_loose_"+InputTag+"_B","h2_tau3p_vmpt_EF_tau16_loose_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_EF_tau16_loose_C = new TH2F("h2_tau3p_vmpt_EF_tau16_loose_"+InputTag+"_C","h2_tau3p_vmpt_EF_tau16_loose_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_EF_tau16_loose_D = new TH2F("h2_tau3p_vmpt_EF_tau16_loose_"+InputTag+"_D","h2_tau3p_vmpt_EF_tau16_loose_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); // EF_tau16_medium h2_tau3p_vmpt_EF_tau16_medium_A = new TH2F("h2_tau3p_vmpt_EF_tau16_medium_"+InputTag+"_A","h2_tau3p_vmpt_EF_tau16_medium_"+InputTag+"_A",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_EF_tau16_medium_B = new TH2F("h2_tau3p_vmpt_EF_tau16_medium_"+InputTag+"_B","h2_tau3p_vmpt_EF_tau16_medium_"+InputTag+"_B",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_EF_tau16_medium_C = new TH2F("h2_tau3p_vmpt_EF_tau16_medium_"+InputTag+"_C","h2_tau3p_vmpt_EF_tau16_medium_"+InputTag+"_C",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); h2_tau3p_vmpt_EF_tau16_medium_D = new TH2F("h2_tau3p_vmpt_EF_tau16_medium_"+InputTag+"_D","h2_tau3p_vmpt_EF_tau16_medium_"+InputTag+"_D",b_vm,min_vm,max_vm,b_pt,min_pt,max_pt); } void terminate() { // no selection ---------- h2_tau_vmpt_A->Delete(); h2_tau_vmpt_B->Delete(); h2_tau_vmpt_C->Delete(); h2_tau_vmpt_D->Delete(); // ----------------------- // EF_tau12_loose h2_tau_vmpt_EF_tau12_loose_A->Delete(); h2_tau_vmpt_EF_tau12_loose_B->Delete(); h2_tau_vmpt_EF_tau12_loose_C->Delete(); h2_tau_vmpt_EF_tau12_loose_D->Delete(); // EF_tau16_loose h2_tau_vmpt_EF_tau16_loose_A->Delete(); h2_tau_vmpt_EF_tau16_loose_B->Delete(); h2_tau_vmpt_EF_tau16_loose_C->Delete(); h2_tau_vmpt_EF_tau16_loose_D->Delete(); // EF_tau16_medium h2_tau_vmpt_EF_tau16_medium_A->Delete(); h2_tau_vmpt_EF_tau16_medium_B->Delete(); h2_tau_vmpt_EF_tau16_medium_C->Delete(); h2_tau_vmpt_EF_tau16_medium_D->Delete(); // no selection ---------- h2_tau1p_vmpt_A->Delete(); h2_tau1p_vmpt_B->Delete(); h2_tau1p_vmpt_C->Delete(); h2_tau1p_vmpt_D->Delete(); // ----------------------- // EF_tau12_loose h2_tau1p_vmpt_EF_tau12_loose_A->Delete(); h2_tau1p_vmpt_EF_tau12_loose_B->Delete(); h2_tau1p_vmpt_EF_tau12_loose_C->Delete(); h2_tau1p_vmpt_EF_tau12_loose_D->Delete(); // EF_tau16_loose h2_tau1p_vmpt_EF_tau16_loose_A->Delete(); h2_tau1p_vmpt_EF_tau16_loose_B->Delete(); h2_tau1p_vmpt_EF_tau16_loose_C->Delete(); h2_tau1p_vmpt_EF_tau16_loose_D->Delete(); // EF_tau16_medium h2_tau1p_vmpt_EF_tau16_medium_A->Delete(); h2_tau1p_vmpt_EF_tau16_medium_B->Delete(); h2_tau1p_vmpt_EF_tau16_medium_C->Delete(); h2_tau1p_vmpt_EF_tau16_medium_D->Delete(); // no selection ---------- h2_tau3p_vmpt_A->Delete(); h2_tau3p_vmpt_B->Delete(); h2_tau3p_vmpt_C->Delete(); h2_tau3p_vmpt_D->Delete(); // ----------------------- // EF_tau12_loose h2_tau3p_vmpt_EF_tau12_loose_A->Delete(); h2_tau3p_vmpt_EF_tau12_loose_B->Delete(); h2_tau3p_vmpt_EF_tau12_loose_C->Delete(); h2_tau3p_vmpt_EF_tau12_loose_D->Delete(); // EF_tau16_loose h2_tau3p_vmpt_EF_tau16_loose_A->Delete(); h2_tau3p_vmpt_EF_tau16_loose_B->Delete(); h2_tau3p_vmpt_EF_tau16_loose_C->Delete(); h2_tau3p_vmpt_EF_tau16_loose_D->Delete(); // EF_tau16_medium h2_tau3p_vmpt_EF_tau16_medium_A->Delete(); h2_tau3p_vmpt_EF_tau16_medium_B->Delete(); h2_tau3p_vmpt_EF_tau16_medium_C->Delete(); h2_tau3p_vmpt_EF_tau16_medium_D->Delete(); } }; #endif