search

veelken / SVfit

SVfit algorithm for tau pair mass reconstruction (CMSSW plugin version)
0 stars 0 forks source link

running the package #1

Open amithep opened 7 years ago

amithep commented 7 years ago

Hi, How to compile and run this package?

Amit

veelken commented 7 years ago

Hi Amit,

in the context of which analysis do you want to use SVfit ?

Cheers,

Christian

amithep commented 7 years ago

In Delphes analysis, I want to use this package. Please let me know if I can do that.

From: Christian Veelken notifications@github.com Sent: 07 November 2017 13:55:07 To: veelken/SVfit Cc: Amit Adhikary; Author Subject: Re: [veelken/SVfit] running the package (#1)

Hi Amit,

in the context of which analysis do you want to use SVfit ?

Cheers,

Christian

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/veelken/SVfit/issues/1#issuecomment-342409150, or mute the threadhttps://github.com/notifications/unsubscribe-auth/Af5ZDTaWNSMsR_eIaBVqLDvPvFjY0a3iks5s0BPigaJpZM4QUX_o.

amithep commented 7 years ago

Actually I want to reconstruct Z mass from, Z -> tau tau decay in Delphes.

veelken commented 7 years ago

Hi Amit,

the code can be used with Delphes. It is written to be used with the "scram" tool for compiling the code. In case you want to compile the code outside CMSSW, you would need to write Makefiles. Please not that the most recent version of SVfit is available via git clone https://github.com/svfit/ClassicSVfit TauAnalysis/ClassicSVfit git clone https://github.com/svfit/SVfitTF TauAnalysis/SVfitTF There exists a test program in https://github.com/SVfit/ClassicSVfit/blob/master/bin/testClassicSVfit.cc Note that you need to replace the line https://github.com/SVfit/ClassicSVfit/blob/master/bin/testClassicSVfit.cc#L49 by double kappa = 4.; // use 3 for emu, 4 for etau and mutau, 5 for tautau channel svFitAlgo.addLogM_fixed(true, kappa); using a value of kappa that depends on whether you are analyzing events with 2 leptons (where "lepton" means electron or muons), 1 lepton and 1 hadronic tau, or 2 hadronic taus.

Cheers,

Christian

veelken commented 7 years ago

Hi Amit,

I confirm that SVfit can reconstruct the Z mass in Z->tautau decays.

Cheeers,

Christian

amithep commented 7 years ago

Thank you very much. I will try with this and get back to you in case I face any problem.

amithep commented 6 years ago

Hi,

In the code "testClassicSVfit.cc" // define MET covariance TMatrixD covMET(2, 2);

What is a covMET matrix?

As you said, compiling SVfit requires writing makefile. Since I am a beginner I am trying. In the meantime I tried to run the test code "testClassicSVfit.cc", in a separate folder bringing all dependent header files (& run it with root -l testClassicSVfit.cc). But I got this error message >

Processing testClassicSVfit.cc... Error: class,struct,union or type ic_svFit not defined svFitAuxFunctions.h:110: Error: class,struct,union or type ic_svFit not defined /home/amit/root/cint/cint/lib/prec_stl/vector:495: unknown classic_svFit::SVfitQuantity, cannot cast Warning: Error occurred during reading source files Warning: Error occurred during dictionary source generation !!!Removing /home/amit/Desktop/SVfit/SVFit/AutoDict_vector_classic_svFitGraphPointcxx_ACLiC_dict.cxx /home/amit/Desktop/SVfit/SVFit/AutoDict_vector_classic_svFitGraphPointcxx_ACLiC_dict.h !!! Error: /home/amit/root/bin/rootcint: error loading headers... Error in : Dictionary generation failed! Warning: Unknown type 'std::vector' in function argument handled as int svFitAuxFunctions.h:110: Interpreter error recovered

I do not find any "ic_svFit" in any .h file, it is "classic_svFit".

veelken commented 6 years ago

Hi Amit,

the MET covariance matrix quantifies the experimental resolution on the missing Et. I suggest you use: covMET[0][0] = res res; covMET[0][1] = 0.; covMET[1][0] = 0.; covMET[1][1] = res res; i.e. assume that the matrix is diagonal. The symbol "res" is around 10-15 GeV for the ATLAS and CMS experiments (with Run 2 pileup conditions).

I am afraid I cannot help much with the compilation errors. Rather than trying to compile the code with interactive ROOT, I would suggest you write a Makefile to compile the code.

I hope that helps,

Christian

amithep commented 6 years ago

Thank you very much. I understand now. Yea, I am trying to make a makefile.

amithep commented 6 years ago

Hi, I am very sorry for bothering you again. Previously you said that svfit can be used outside CMSSW. But when I was trying to make files in TauAnalysis(ClassicSVfit & SVfitTF) dependent on "HadTauTFfromTGraph.cc", it requires "FWCore/ParameterSet/interface/FileInPath.h" which is inside CMSSW, which again requires other dependencies inside CMSSW. I am really confused now how to connect all the headers and sources in TauAnalysis without CMSSW.

veelken commented 6 years ago

Hi Amit, I suggest you compile https://github.com/veelken/ClassicSVfit/blob/master/interface/ClassicSVfit.h with #ifdef USE_SVFITTF not defined, so that the code gets compiled without the dependency on https://github.com/veelken/SVfitTF . Would that work for you ?

amithep commented 6 years ago

So only "ClassicSVfit" is enough for this calculation. Ok. Thanks.

amithep commented 6 years ago

Hi, Many thanks to you. I wrote a Makefile and linked the source and header files. So, for running the test code should I use 'root -l testClassicSVfit.cc' ?

veelken commented 6 years ago

Hi Amit, my recommendation is that you compile the code testClassicSVfit.cc as an executable in the Makefile (such that you can simply execute 'testClassicSVfit' in the terminal, i.e. without interactive ROOT).

amithep commented 6 years ago

ok.

amithep commented 6 years ago

I compiled the test code with the object files (.o) I created earlier in makefile. But it says 'cannot execute binary file'. Should I make an executable from all object files and execute code? (like ./exe testClassicSVfit.cc)

amithep commented 6 years ago

hi, I have a question. What is the mistag rate of a jet faking tau-tagged jet in CMS?

regards, amit

Get Outlook for Androidhttps://aka.ms/ghei36

From: Christian Veelken notifications@github.com Sent: Tuesday, November 7, 2017 1:55:07 PM To: veelken/SVfit Cc: Amit Adhikary; Author Subject: Re: [veelken/SVfit] running the package (#1)

Hi Amit,

in the context of which analysis do you want to use SVfit ?

Cheers,

Christian

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/veelken/SVfit/issues/1#issuecomment-342409150, or mute the threadhttps://github.com/notifications/unsubscribe-auth/Af5ZDTaWNSMsR_eIaBVqLDvPvFjY0a3iks5s0BPigaJpZM4QUX_o.