void proof(Int_t runNum=-1,Int_t NumCores=10,TString OutPutName="histograms/temp.root") { //////////////////////////////////////////////////////////// // open the root files and put the together into a TChain // //////////////////////////////////////////////////////////// if (runNum != -99){ ch = new TChain("caltree"); std::cout << "Creating TChain... "; if (runNum ==-1){ // if no runNum give make a chain manually ch->Add("./rootfiles/run-0573-??.root"); } else { stringstream ss; ss<<"./rootfiles/run-"<<setfill('0')<<setw(4)<<runNum<<"-??.root"; ch->Add(ss.str().c_str()); if (OutPutName == "histograms/temp.root"){ ss.str(""); ss<<"histograms/HistogramsFromRun"<<runNum<<".root"; OutPutName=ss.str(); } } } std::cout << "Done." << std::endl; /////////////////////////////////////////////////////////////// // look for settings objects in each of the files // // in the chain. Set their names as settings0 1 2... // // for saving in the final histogram file // /////////////////////////////////////////////////////////////// vector <R00TLeSettings *> ListOfSettings; TObjArray *fileElements=ch->GetListOfFiles(); TIter next(fileElements); TChainElement *chEl=0; int FileSuccessCount=0; int count=0; while (( chEl=(TChainElement*)next() )) { TFile f(chEl->GetTitle()); if (f.IsOpen()){ FileSuccessCount++; }else { cout<<"Warning could not open "<<chEl->GetTitle()<<endl; } if ( f.GetListOfKeys()->Contains("TheSettings")) { stringstream ss; ss<<"Settings"<<count; TString newName = ss.str().c_str(); R00TLeSettings *temp = (R00TLeSettings*)f.Get("TheSettings"); temp->SetName(newName); ListOfSettings.push_back(temp); } else { cout<<"No Seetings object found in "<<chEl->GetTitle()<<endl; exit(1); } count++; } if ( FileSuccessCount == 0 ){ cout<<"No files were successfully opened"<<endl; return; } std::cout << "\nList of files:" << std::endl; ch->GetListOfFiles()->Print(); cout<<endl; TString user =gSystem->Getenv("R00TLe_User"); TString install =gSystem->Getenv("R00TLeInstall"); TString src = install+"/users/"+user+"/src/Analoop.cc+g"; TString user =gSystem->Getenv("R00TLe_User"); TString install =gSystem->Getenv("R00TLeInstall"); TString src = install+"/users/"+user+"/src/Analoop.cc+g"; TString soFile=install+"/users/"+user+"/src/Analoop_cc.so"; ///////////////////////////////////////// // compile and load Analoop in to ROOT // ///////////////////////////////////////// gROOT->ProcessLine(".L " +src); //////////////////////////////////////// // Set up the PROOF Sandbox directory // //////////////////////////////////////// TString sandbox = install + "/users/"+user+"/PROOF/"; gEnv->SetValue("ProofLite.SandBox",sandbox); ///////////////// // start PROOF // ///////////////// std::stringstream ss; ss<<NumCores; TString workerString = "workers=" + TString(ss.str().c_str()); TProof * p =TProof::Open(workerString); p->SetParameter("PROOF_UseTreeCache",1); p->SetParameter("PROOF_CacheSize",1024.*1024.); ////////////////////////////////////////////////////// // load the analoop library on each worker thread // ////////////////////////////////////////////////////// TString so_cmd = "gSystem->Load(\"" +soFile+ "\")"; p->Exec(so_cmd); ///////////////////////////////////////////////////////////// // make instance of analoop that will be passed into PROOF // ///////////////////////////////////////////////////////////// TSelector* myAnaloop = TSelector::GetSelector(src); ////////////////////////////////////////////////////////////// // Put the settings object from the first file in the input // // list of the selector // ////////////////////////////////////////////////////////////// TList * list = new TList(); list->Add(ListOfSettings[0]); myAnaloop->SetInputList(list); //////////////////////////////////////////////////////////////// // Magic lines that loop over the tchain and analyzes it with // // the analoop program using PROOF on different worker // // nodes // //////////////////////////////////////////////////////////////// ch->SetProof(kTRUE); ch->SetCacheSize(1024.*1024.); ch->AddBranchToCache("*",kTRUE); ch->Process(myAnaloop,OutPutName);//Deep Magicks //Preform Clean up for PROOF p->CleanupSession(gProof->GetSessionTag()); p->ClearCache(); p->Close(); delete ch; }
void InitAndLoadLibs(Int_t runMode=kMyRunModeLocal, Int_t workers=0,Bool_t debug=0) { // Loads libs and par files + custom task and classes // Custom stuff to be loaded // listToLoad->Add(new TObjString("$ALICE_ROOT/ANALYSIS/AliCentralitySelectionTask.cxx+")); listToLoad->Add(new TObjString("/Users/mfloris/Work/AliSoft/AliRoot-trunk/PWGPP/background/AliHistoListWrapper.cxx+")); // FIXME listToLoad->Add(new TObjString("AliAnalysisMultPbTrackHistoManager.cxx+")); listToLoad->Add(new TObjString("AliAnalysisMultPbCentralitySelector.cxx+")); listToLoad->Add(new TObjString("AliAnalysisTaskMultPbTracks.cxx+")); // listToLoad->Add(new TObjString("$ALICE_ROOT/PWG0/multPbPb/AliAnalysisMultPbTrackHistoManager.cxx+")); // listToLoad->Add(new TObjString("$ALICE_ROOT/PWG0/multPbPb/AliAnalysisMultPbCentralitySelector.cxx+")); // listToLoad->Add(new TObjString("$ALICE_ROOT/PWG0/multPbPb/AliAnalysisTaskMultPbTracks.cxx+")); if (runMode == kMyRunModeCAF) { cout << "Init in CAF mode" << endl; gEnv->SetValue("XSec.GSI.DelegProxy", "2"); // TProof::Mgr("alice-caf.cern.ch")->SetROOTVersion("VO_ALICE@ROOT::v5-30-03-1"); TProof::Mgr("alice-caf.cern.ch")->SetROOTVersion("current"); TProof * p = TProof::Open("alice-caf.cern.ch", workers>0 ? Form("workers=%d",workers) : "1x"); // TProof * p = TProof::Open("skaf.saske.sk", workers>0 ? Form("workers=%d",workers) : ""); p->Exec("TObject *o = gEnv->GetTable()->FindObject(\"Proof.UseMergers\"); gEnv->GetTable()->Remove(o);", kTRUE); // TProof::Mgr("alice-caf.cern.ch")->SetROOTVersion("VO_ALICE@ROOT::v5-28-00f"); gProof->EnablePackage("VO_ALICE@AliRoot::v5-02-12-AN"); gSystem->Load("libCore"); gSystem->Load("libTree"); gSystem->Load("libGeom"); gSystem->Load("libVMC"); gSystem->Load("libPhysics"); gSystem->Load("libMinuit"); gSystem->Load("libSTEERBase"); gSystem->Load("libESD"); gSystem->Load("libAOD"); gSystem->Load("libANALYSIS"); gSystem->Load("libOADB"); gSystem->Load("libANALYSISalice"); // Enable the needed package // gProof->UploadPackage("$ALICE_ROOT/obj/STEERBase"); // gProof->EnablePackage("$ALICE_ROOT/obj/STEERBase"); // gProof->UploadPackage("$ALICE_ROOT/obj/ESD"); // gProof->EnablePackage("$ALICE_ROOT/obj/ESD"); // gProof->UploadPackage("$ALICE_ROOT/obj/AOD"); // gProof->EnablePackage("$ALICE_ROOT/obj/AOD"); // gProof->UploadPackage("$ALICE_ROOT/obj/ANALYSIS"); // gProof->EnablePackage("$ALICE_ROOT/obj/ANALYSIS"); // gProof->UploadPackage("$ALICE_ROOT/obj/OADB"); // gProof->EnablePackage("$ALICE_ROOT/obj/OADB"); // gProof->UploadPackage("$ALICE_ROOT/obj/ANALYSISalice"); // gProof->EnablePackage("$ALICE_ROOT/obj/ANALYSISalice"); // gProof->UploadPackage("$ALICE_ROOT/obj/PWG0base"); // gProof->EnablePackage("$ALICE_ROOT/obj/PWG0base"); // gROOT->ProcessLine(gSystem->ExpandPathName(".include $ALICE_ROOT/include")); // gROOT->ProcessLine(gSystem->ExpandPathName(".include $ALICE_ROOT/PWG0/multPb")); // gROOT->ProcessLine(gSystem->ExpandPathName(".include $ALICE_ROOT/PWG1/background")); } else if (runMode == kMyRunModeProofLite) { cout << "Init in CAF mode" << endl; gEnv->SetValue("XSec.GSI.DelegProxy", "2"); TProof * p = TProof::Open(""); // TProof * p = TProof::Open("skaf.saske.sk", workers>0 ? Form("workers=%d",workers) : ""); // p->Exec("TObject *o = gEnv->GetTable()->FindObject(\"Proof.UseMergers\"); gEnv->GetTable()->Remove(o);", kTRUE); // TProof::Mgr("alice-caf.cern.ch")->SetROOTVersion("VO_ALICE@ROOT::v5-28-00f"); // TProof::Mgr("alice-caf.cern.ch")->SetROOTVersion("5.28/00f"); gProof->UploadPackage("$ALICE_ROOT/ANALYSIS/macros/AliRootProofLite.par"); gProof->EnablePackage("AliRootProofLite"); } else { cout << "Init in Local or Grid mode" << endl; gSystem->Load("libCore"); gSystem->Load("libTree"); gSystem->Load("libGeom"); gSystem->Load("libVMC"); gSystem->Load("libPhysics"); gSystem->Load("libMinuit"); gSystem->Load("libSTEERBase"); gSystem->Load("libESD"); gSystem->Load("libAOD"); gSystem->Load("libANALYSIS"); gSystem->Load("libOADB"); gSystem->Load("libANALYSISalice"); // Use AliRoot includes to compile our task gROOT->ProcessLine(".include $ALICE_ROOT/include"); // gSystem->Load("libVMC"); // gSystem->Load("libTree"); // gSystem->Load("libSTEERBase"); // gSystem->Load("libESD"); // gSystem->Load("libAOD"); // gSystem->Load("libANALYSIS"); // gSystem->Load("libANALYSISalice"); // gSystem->Load("libPWG0base"); gROOT->ProcessLine(gSystem->ExpandPathName(".include $ALICE_ROOT/PWG0/multPb")); gROOT->ProcessLine(gSystem->ExpandPathName(".include $ALICE_ROOT/PWGPP/background")); // gROOT->ProcessLine(gSystem->ExpandPathName(".include $ALICE_ROOT/PWGPP/background/")); } // Load helper classes TIterator * iter = listToLoad->MakeIterator(); TObjString * name = 0; while (name = (TObjString *)iter->Next()) { gSystem->ExpandPathName(name->String()); cout << name->String().Data(); if (runMode == kMyRunModeCAF || runMode == kMyRunModeProofLite) { gProof->Load(name->String()+(debug?"+g":"")); } else { gROOT->LoadMacro(name->String()+(debug?"+g":"")); } } }
void runAnalysis(TString sample="all", TString proc="all", TString sys="all", TString updown="both", bool getLogs=false) { TString desdir = "/uscms/home/algomez/work/CMSSW_4_2_4/src/Yumiceva/TreeAnalyzer/test/backup/"; //TString desdir = "/eos/uscms/store/user/algomez/FourTop/resultsTreeAnalyzer/results13/"; TProof *p = TProof::Open("lite://", desdir ,desdir); //p->AddDynamicPath(""); p->Exec("gSystem->Load(\"/uscms/home/algomez/work/CMSSW_4_2_4/lib/slc5_amd64_gcc434/libYumicevaTop7TeV.so\")"); p->Exec("gSystem->Load(\"/uscms/home/algomez/work/CMSSW_4_2_4/lib/slc5_amd64_gcc434/libCondFormatsJetMETObjects.so\")"); p->AddIncludePath("/uscms/home/algomez/work/CMSSW_4_2_4/src/"); p->Archive(" ",desdir); //p->AddInput(new TNamed("PROOF_OUTPUTFILE_LOCATION", "LOCAL")); if (sample=="tttt"|| sample=="ttttSM"||sample=="all") { TDSet *mc_tttt_SM = new TDSet("top","*","/PATNtupleMaker"); mc_tttt_SM->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_SM-Tuple-PATskim.root"); if (proc =="all" ) { mc_tttt_SM->Process("Analyzer.C+","sample=tttt_SM"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_tttt_SM->Process("Analyzer.C+", sys+"UP sample=tttt_SM"); if (updown == "both" || updown == "down") mc_tttt_SM->Process("Analyzer.C+", sys+"DOWN sample=tttt_SM"); } } if (sample=="tttt"||sample=="tttt400"||sample=="all") { TDSet *mc_tttt_Gh400 = new TDSet("top","*","/PATNtupleMaker"); mc_tttt_Gh400->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh400_50k_1-Tuple-PATskim.root"); mc_tttt_Gh400->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh400_50k_2-Tuple-PATskim.root"); if (proc =="all" ) { mc_tttt_Gh400->Process("Analyzer.C+","sample=tttt_Gh400"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_tttt_Gh400->Process("Analyzer.C+", sys+"UP sample=tttt_Gh400"); if (updown == "both" || updown == "down") mc_tttt_Gh400->Process("Analyzer.C+", sys+"DOWN sample=tttt_Gh400"); } } if (sample=="tttt"||sample=="tttt500"||sample=="all") { TDSet *mc_tttt_Gh500 = new TDSet("top","*","/PATNtupleMaker"); mc_tttt_Gh500->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh500_50k_1-Tuple-PATskim.root"); mc_tttt_Gh500->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh500_50k_2-Tuple-PATskim.root"); if (proc =="all" ) { mc_tttt_Gh500->Process("Analyzer.C+","sample=tttt_Gh500"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_tttt_Gh500->Process("Analyzer.C+", sys+"UP sample=tttt_Gh500"); if (updown == "both" || updown == "down") mc_tttt_Gh500->Process("Analyzer.C+", sys+"DOWN sample=tttt_Gh500"); } } if (sample=="tttt"||sample=="tttt600"||sample=="all") { TDSet *mc_tttt_Gh600 = new TDSet("top","*","/PATNtupleMaker"); //mc_tttt_Gh600->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh600_50k_1-Tuple-PATskim.root"); mc_tttt_Gh600->Add("/uscms/home/algomez/work/CMSSW_4_2_4/src/Yumiceva/Top7TeV/test/tttt_Gh600_test_FASTSIM_FULLSIM/tttt_Gh600-Tuple-PATskim.root"); //mc_tttt_Gh600->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh600_50k_2-Tuple-PATskim.root"); if (proc =="all" ) { mc_tttt_Gh600->Process("Analyzer.C+","sample=tttt_Gh600"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_tttt_Gh600->Process("Analyzer.C+", sys+"UP sample=tttt_Gh600"); if (updown == "both" || updown == "down") mc_tttt_Gh600->Process("Analyzer.C+", sys+"DOWN sample=tttt_Gh600"); } } if (sample=="tttt"||sample=="tttt700"||sample=="all") { TDSet *mc_tttt_Gh700 = new TDSet("top","*","/PATNtupleMaker"); mc_tttt_Gh700->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh700_50k_1-Tuple-PATskim.root"); mc_tttt_Gh700->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh700_50k_2-Tuple-PATskim.root"); if (proc =="all" ) { mc_tttt_Gh700->Process("Analyzer.C+","sample=tttt_Gh700"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_tttt_Gh700->Process("Analyzer.C+", sys+"UP sample=tttt_Gh700"); if (updown == "both" || updown == "down") mc_tttt_Gh700->Process("Analyzer.C+", sys+"DOWN sample=tttt_Gh700"); } } if (sample=="tttt"||sample=="tttt800"||sample=="all") { TDSet *mc_tttt_Gh800 = new TDSet("top","*","/PATNtupleMaker"); mc_tttt_Gh800->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh800_50k_1-Tuple-PATskim.root"); mc_tttt_Gh800->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh800_50k_2-Tuple-PATskim.root"); if (proc =="all" ) { mc_tttt_Gh800->Process("Analyzer.C+","sample=tttt_Gh800"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_tttt_Gh800->Process("Analyzer.C+", sys+"UP sample=tttt_Gh800"); if (updown == "both" || updown == "down") mc_tttt_Gh800->Process("Analyzer.C+", sys+"DOWN sample=tttt_Gh800"); } } if (sample=="tttt"||sample=="tttt900"||sample=="all") { TDSet *mc_tttt_Gh900 = new TDSet("top","*","/PATNtupleMaker"); mc_tttt_Gh900->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh900_50k_1-Tuple-PATskim.root"); mc_tttt_Gh900->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh900_50k_2-Tuple-PATskim.root"); if (proc =="all" ) { mc_tttt_Gh900->Process("Analyzer.C+","sample=tttt_Gh900"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_tttt_Gh900->Process("Analyzer.C+", sys+"UP sample=tttt_Gh900"); if (updown == "both" || updown == "down") mc_tttt_Gh900->Process("Analyzer.C+", sys+"DOWN sample=tttt_Gh900"); } } if (sample=="tttt"||sample=="tttt1000"||sample=="all") { TDSet *mc_tttt_Gh1000 = new TDSet("top","*","/PATNtupleMaker"); mc_tttt_Gh1000->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh1000_50k_1-Tuple-PATskim.root"); mc_tttt_Gh1000->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh1000_50k_2-Tuple-PATskim.root"); if (proc =="all" ) { mc_tttt_Gh1000->Process("Analyzer.C+","sample=tttt_Gh1000"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_tttt_Gh1000->Process("Analyzer.C+", sys+"UP sample=tttt_Gh1000"); if (updown == "both" || updown == "down") mc_tttt_Gh1000->Process("Analyzer.C+", sys+"DOWN sample=tttt_Gh1000"); } } if (sample=="tttt"||sample=="ttttUED6"||sample=="all") { TDSet *mc_tttt_UED6 = new TDSet("top","*","/PATNtupleMaker"); mc_tttt_UED6->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_UED6-Tuple-PATskim.root"); if (proc =="all" ) { mc_tttt_UED6->Process("Analyzer.C+","sample=ttttUED6"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_tttt_UED6->Process("Analyzer.C+", sys+"UP sample=ttttUED6"); if (updown == "both" || updown == "down") mc_tttt_UED6->Process("Analyzer.C+", sys+"DOWN sample=ttttUED6"); } } /*if (sample=="tttt"||sample=="tttt1200"||sample=="all") { TDSet *mc_tttt_Gh1200 = new TDSet("top","*","/PATNtupleMaker"); mc_tttt_Gh1200->Add("/eos/uscms/store/user/algomez/FourTop/Ntuples/tttt_Gh1200-Tuple-PATskim.root"); if (proc =="all" ) { mc_tttt_Gh1200->Process("Analyzer.C+","sample=tttt_Gh1200"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_tttt_Gh1200->Process("Analyzer.C+", sys+"UP sample=tttt_Gh1200"); if (updown == "both" || updown == "down") mc_tttt_Gh1200->Process("Analyzer.C+", sys+"DOWN sample=tttt_Gh1200"); } }*/ if (sample=="MC"||sample=="ttbar"||sample=="all") { TDSet *mc_ttbar = new TDSet("top","*","/PATNtupleMaker"); mc_ttbar->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/TTbar_Fall11_Mu.root"); ///uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/TTbar_Mu.root"); if (proc =="all" ) { mc_ttbar->Process("Analyzer.C+","sample=ttbar"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_ttbar->Process("Analyzer.C+", sys+"UP sample=ttbar"); if (updown == "both" || updown == "down") mc_ttbar->Process("Analyzer.C+", sys+"DOWN sample=ttbar"); } } if (sample=="TTbar"||sample=="all") { TDSet *mc_TTbar = new TDSet("top","*","/PATNtupleMaker"); mc_TTbar->Add("/eos/uscms/store/user/algomez/ttbar/ttbar_pythia8-Tuple-PATskim.root"); mc_TTbar->Add("/eos/uscms/store/user/algomez/ttbar/ttbar_pythia8_100k-Tuple-PATskim.root"); if (proc =="all" ) { mc_TTbar->Process("Analyzer.C+","sample=TTbar"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_TTbar->Process("Analyzer.C+", sys+"UP sample=TTbar"); if (updown == "both" || updown == "down") mc_TTbar->Process("Analyzer.C+", sys+"DOWN sample=TTbar"); } } if (sample=="TTbar_POWHEG"||sample=="all") { TDSet *mc_ttbar = new TDSet("top","*","/PATNtupleMaker"); mc_ttbar->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/TT_Powheg_Mu.root"); if (proc =="all" ) mc_ttbar->Process("Analyzer.C+","sample=ttbar_powheg"); if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_ttbar->Process("Analyzer.C+", sys+"UP sample=ttbar_powheg"); if (updown == "both" || updown == "down") mc_ttbar->Process("Analyzer.C+", sys+"DOWN sample=ttbar_powheg"); } } if (sample=="TTbarM"||sample=="MC") { TDSet *mc_ttbarm = new TDSet("top","*","/PATNtupleMaker"); mc_ttbarm->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/TTbar_matchup_Mu.root"); if (proc =="all" ) {mc_ttbarm->Process("Analyzer.C+","sample=ttbarMatchup");} } if (sample=="TTbarM"||sample=="MC") { TDSet *mc_ttbarm = new TDSet("top","*","/PATNtupleMaker"); mc_ttbarm->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/TTbar_matchdown_Mu.root"); if (proc =="all" ) {mc_ttbarm->Process("Analyzer.C+","sample=ttbarMatchdown");} } if (sample=="TTbarQ"||sample=="MC") { TDSet *mc_ttbarm = new TDSet("top","*","/PATNtupleMaker"); mc_ttbarm->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/TTbar_scaleup_Mu.root"); if (proc =="all" ) {mc_ttbarm->Process("Analyzer.C+","sample=ttbarScaleup");} } if (sample=="TTbarQ"||sample=="MC") { TDSet *mc_ttbarm = new TDSet("top","*","/PATNtupleMaker"); mc_ttbarm->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/TTbar_scaledown_Mu.root"); if (proc =="all" ) {mc_ttbarm->Process("Analyzer.C+","sample=ttbarScaledown");} } if (sample=="WJets"||sample=="all") { TDSet *mc_WJets = new TDSet("top","*","/PATNtupleMaker"); mc_WJets->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/WJets_Mu.root"); if (proc =="all" ) { mc_WJets->Process("Analyzer.C+","sample=WJets"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_WJets->Process("Analyzer.C+", sys+"UP sample=WJets"); if (updown == "both" || updown == "down") mc_WJets->Process("Analyzer.C+", sys+"DOWN sample=WJets"); } } if (sample=="MC"||sample=="Wjets"||sample=="Wbb"||sample=="all") { TDSet *mc_Wbb = new TDSet("top","*","/PATNtupleMaker"); mc_Wbb->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/WJets_Mu.root"); if (proc =="all" ) { mc_Wbb->Process("Analyzer.C+","sample=Wbb"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_Wbb->Process("Analyzer.C+", sys+"UP sample=Wbb"); if (updown == "both" || updown == "down") mc_Wbb->Process("Analyzer.C+", sys+"DOWN sample=Wbb"); } } if (sample=="MC"||sample=="Wjets"||sample=="Wcc"||sample=="all") { TDSet *mc_Wcc = new TDSet("top","*","/PATNtupleMaker"); mc_Wcc->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/WJets_Mu.root"); if (proc =="all" ) { mc_Wcc->Process("Analyzer.C+","sample=Wcc"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_Wcc->Process("Analyzer.C+", sys+"UP sample=Wcc"); if (updown == "both" || updown == "down") mc_Wcc->Process("Analyzer.C+", sys+"DOWN sample=Wcc"); } } if (sample=="MC"||sample=="Wjets"||sample=="Wqq"||sample=="all") { TDSet *mc_Wqq = new TDSet("top","*","/PATNtupleMaker"); mc_Wqq->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/WJets_Mu.root"); if (proc =="all" ) { mc_Wqq->Process("Analyzer.C+","sample=Wqq"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_Wqq->Process("Analyzer.C+", sys+"UP sample=Wqq"); if (updown == "both" || updown == "down") mc_Wqq->Process("Analyzer.C+", sys+"DOWN sample=Wqq"); } } if (sample=="QCD"||sample=="all") { TDSet *mc_QCD = new TDSet("top","*","/PATNtupleMaker"); mc_QCD->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/QCD_Mu.root"); if (proc =="all" ) { mc_QCD->Process("Analyzer.C+","sample=QCD"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_QCD->Process("Analyzer.C+", sys+"UP sample=QCD"); if (updown == "both" || updown == "down") mc_QCD->Process("Analyzer.C+", sys+"DOWN sample=QCD"); } } if (sample=="MC"||sample=="ST"||sample=="STsch"||sample=="all") { TDSet *mc_STsch = new TDSet("top","*","/PATNtupleMaker"); mc_STsch->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STsch_Mu.root"); if (proc =="all" ) { mc_STsch->Process("Analyzer.C+","sample=STsch"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_STsch->Process("Analyzer.C+", sys+"UP sample=STsch"); if (updown == "both" || updown == "down") mc_STsch->Process("Analyzer.C+", sys+"DOWN sample=STsch"); } } if (sample=="MC"||sample=="ST"||sample=="STtch"||sample=="all") { TDSet *mc_STtch = new TDSet("top","*","/PATNtupleMaker"); mc_STtch->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STtch_Mu.root"); if (proc =="all" ) { mc_STtch->Process("Analyzer.C+","sample=STtch"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_STtch->Process("Analyzer.C+", sys+"UP sample=STtch"); if (updown == "both" || updown == "down") mc_STtch->Process("Analyzer.C+", sys+"DOWN sample=STtch"); } } if (sample=="MC"||sample=="ST"||sample=="STtWch"||sample=="all") { TDSet *mc_STtWch = new TDSet("top","*","/PATNtupleMaker"); mc_STtWch->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STtWch_Mu.root"); if (proc =="all" ) { mc_STtWch->Process("Analyzer.C+","sample=STtWch"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_STtWch->Process("Analyzer.C+", sys+"UP sample=STtWch"); if (updown == "both" || updown == "down") mc_STtWch->Process("Analyzer.C+", sys+"DOWN sample=STtWch"); } } if (sample=="MC"||sample=="ST"||sample=="STsch_bar"||sample=="all") { TDSet *mc_STsch_bar = new TDSet("top","*","/PATNtupleMaker"); mc_STsch_bar->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STsch_bar_Mu.root"); if (proc =="all" ) { mc_STsch_bar->Process("Analyzer.C+","sample=STsch_bar"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_STsch_bar->Process("Analyzer.C+", sys+"UP sample=STsch_bar"); if (updown == "both" || updown == "down") mc_STsch_bar->Process("Analyzer.C+", sys+"DOWN sample=STsch_bar"); } } if (sample=="MC"||sample=="ST"||sample=="STtch_bar"||sample=="all") { TDSet *mc_STtch_bar = new TDSet("top","*","/PATNtupleMaker"); mc_STtch_bar->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STtch_bar_Mu.root"); if (proc =="all" ) { mc_STtch_bar->Process("Analyzer.C+","sample=STtch_bar"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_STtch_bar->Process("Analyzer.C+", sys+"UP sample=STtch_bar"); if (updown == "both" || updown == "down") mc_STtch_bar->Process("Analyzer.C+", sys+"DOWN sample=STtch_bar"); } } if (sample=="MC"||sample=="ST"||sample=="STtWch_bar"||sample=="all") { TDSet *mc_STtWch_bar = new TDSet("top","*","/PATNtupleMaker"); mc_STtWch_bar->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STtWch_bar_Mu.root"); if (proc =="all" ) { mc_STtWch_bar->Process("Analyzer.C+","sample=STtWch_bar"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_STtWch_bar->Process("Analyzer.C+", sys+"UP sample=STtWch_bar"); if (updown == "both" || updown == "down") mc_STtWch_bar->Process("Analyzer.C+", sys+"DOWN sample=STtWch_bar"); } } if (sample=="MC"||sample=="diboson"||sample=="WW"||sample=="all") { TDSet *mc_WW = new TDSet("top","*","/PATNtupleMaker"); mc_WW->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/WW_Mu.root"); if (proc =="all" ) { mc_WW->Process("Analyzer.C+","sample=WW"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_WW->Process("Analyzer.C+", sys+"UP sample=WW"); if (updown == "both" || updown == "down") mc_WW->Process("Analyzer.C+", sys+"DOWN sample=WW"); } } if (sample=="MC"||sample=="diboson"||sample=="WZ"||sample=="all") { TDSet *mc_WZ = new TDSet("top","*","/PATNtupleMaker"); mc_WZ->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/WZ_Mu.root"); if (proc =="all" ) { mc_WZ->Process("Analyzer.C+","sample=WZ"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_WZ->Process("Analyzer.C+", sys+"UP sample=WZ"); if (updown == "both" || updown == "down") mc_WZ->Process("Analyzer.C+", sys+"DOWN sample=WZ"); } } if (sample=="MC"||sample=="ZJets"||sample=="all") { TDSet *mc_ZJets = new TDSet("top","*","/PATNtupleMaker"); mc_ZJets->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/ZJets_Mu.root"); if (proc =="all" ) { mc_ZJets->Process("Analyzer.C+","sample=ZJets"); } if ( proc =="sys" ) { if (updown == "both" || updown == "up") mc_ZJets->Process("Analyzer.C+", sys+"UP sample=ZJets"); if (updown == "both" || updown == "down") mc_ZJets->Process("Analyzer.C+", sys+"DOWN sample=ZJets"); } } if (sample=="data"||sample=="all") { TDSet *data = new TDSet("top","*","/PATNtupleMaker"); data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_May10ReReco_Mu.root"); data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v4_Run_165088_166967_Mu.root"); data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v4_Run_167039_167913_Mu.root"); //data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Aug05ReReco_Mu.root"); //data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v6_Run_172620_173198_Mu.root"); //data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v6_Run_173236_173692_Mu.root"); //data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Run2011B-PromptReco-v1_175860_178380_Mu.root"); //data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Run2011B-PromptReco-v1_178420_180252_Mu.root"); data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Aug05ReReco_Mu.root"); data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v6_Run_172620_173198_Mu.root"); data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v6_Run_173236_173692_Mu.root"); data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Run2011B-PromptReco-v1_175860_177452_Mu.root"); data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Run2011B-PromptReco-v1_177718_178380_Mu.root"); data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Run2011B-PromptReco-v1_178420_180252_Mu.root "); data->Process("Analyzer.C+","sample=data"); // get log files if (getLogs) { logList = p->GetManager()->GetSessionLogs()->GetListOfLogs(); for (int i=1; i< logList->GetSize(); ++i) { logElem = ( TProofLogElem* ) logList->At( i ); macro = logElem->GetMacro(); macro->SaveSource("data_muons_"+TString(Form("%i",i))+".stdout"); } } } if (sample=="dataMissing"||sample=="all") { TDSet *data = new TDSet("top","*","/PATNtupleMaker"); data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/missing.root "); data->Process("Analyzer.C+","sample=dataMissing"); // get log files if (getLogs) { logList = p->GetManager()->GetSessionLogs()->GetListOfLogs(); for (int i=1; i< logList->GetSize(); ++i) { logElem = ( TProofLogElem* ) logList->At( i ); macro = logElem->GetMacro(); macro->SaveSource("data_missing_"+TString(Form("%i",i))+".stdout"); } } } if (sample=="dataQCD2"||sample=="all") { TDSet *dataQCD2 = new TDSet("top","*","/PATNtupleMaker"); dataQCD2->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011/42x/Run2011_May10ReReco_Mu.root"); dataQCD2->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011/42x/Run2011_PromptReco1_Mu.root"); dataQCD2->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011/42x/Run2011_PromptReco2_Mu.root"); dataQCD2->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011/42x/SingleMu_Aug05ReReco_Mu.root"); dataQCD2->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011/42x/SingleMu_PromptReco-v6_Mu.root"); dataQCD2->Process("Analyzer.C+","QCD2 sample=data"); } }